Научная статья на тему 'Моделирование процессов тепло- (массо-) переноса в микроканалах при чип-реализации электрофореза'

Моделирование процессов тепло- (массо-) переноса в микроканалах при чип-реализации электрофореза Текст научной статьи по специальности «Математика»

CC BY
62
8
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Научное приборостроение
ВАК
RSCI
Область наук

Аннотация научной статьи по математике, автор научной работы — Буляница А. Л.

В работе анализируется роль этапа математического моделирования процессов переноса тепла и вещества для решения общей задачи разработки микрофлюидной аналитической системы и проведения на ней химического и биологического анализа. Среди множества задач моделирования выделена следующая: исследование эффективности способов управления анализом многокомпонентных смесей химических и биологических веществ при чип-реализации электрофореза. Также рассмотрены сопутствующие вопросы, связанные с оптимизацией различных стадий анализа, включая обработку информативных сигналов.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Моделирование процессов тепло- (массо-) переноса в микроканалах при чип-реализации электрофореза»

ISSN 0868-5886 НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2005, том 15, № 3, c. 32-39

— Материалы научного семинара

МИКРОЧИПОВЫЕ ТЕХНОЛОГИИ В АНАЛИТИЧЕСКОЙ ХИМИИ

УДК 517.958+ 532.5+ 536.252 © А. Л. Буляница

МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ТЕПЛО- (МАССО-) ПЕРЕНОСА В МИКРОКАНАЛАХ ПРИ ЧИП-РЕАЛИЗАЦИИ ЭЛЕКТРОФОРЕЗА

В работе анализируется роль этапа математического моделирования процессов переноса тепла и вещества для решения общей задачи разработки микрофлюидной аналитической системы и проведения на ней химического и биологического анализа. Среди множества задач моделирования выделена следующая: исследование эффективности способов управления анализом многокомпонентных смесей химических и биологических веществ при чип-реализации электрофореза. Также рассмотрены сопутствующие вопросы, связанные с оптимизацией различных стадий анализа, включая обработку информативных сигналов.

ВВЕДЕНИЕ

Традиционно при решении задач проектирования прибора или разработки новой методики анализа математическое моделирование базовых процессов или явлений применяется как вспомогательная стадия работы. Тем самым результаты, достигнутые на стадии математического моделирования, крайне редко непосредственно проявляются в конечном результате. В частности, среди публикаций, касающихся разработки микрофлюидной аналитической системы, проводимой в Лаборатории информационно-измерительных био-и хемосенсорных микросистем Института аналитического приборостроения РАН [1-14], лишь меньшая их часть [3, 6, 12-14] достаточно подробно рассматривает постановку задач математического моделирования и используемые методы решения.

В контексте проводимых лабораторией работ основной задачей моделирования является исследование различных способов управления процессами анализа (разделения) смесей в каналах микрочипа при реализации электрофореза. Таким образом, главное внимание уделялось проверке эффективности традиционных и новых методов управления микропотоками вещества, выбору управляющих параметров/характеристик как конструктивных элементов системы, так и объектов анализа, а также поиску возможных оптимальных режимов.

Как правило, критерием оптимизации был максимум разрешения аналитических пиков компонент, и сама процедура решения также была хорошо известна в силу очевидной математической формализации задачи — поиск условного экстре -мума функции нескольких переменных при ограничениях в форме неравенств.

Следует выделить ограничения различной при-

роды (происхождения):

— ограничения, связанные с практическим проведением анализа (ограничение на максимальное время анализа, с чем могут быть связаны ограничение сверху длины микроканала и ограничение снизу величины управляющего продольного электрического поля);

— ограничение, связанное с тепловым действием электрического поля (ограничение сверху на величину тепловой мощности);

— ограничение, связанное с чувствительностью детектора (ограничение снизу на объем пробы или/и количество продукта реакции, если таковая осуществляется в микроканале);

— конструктивно-технологические ограничения по выбору размеров микроканала (ширина, соотношение глубины и ширины, угол наклона стенок и т. д.), связанные с уровнем развития технологии;

— финансово-экономические ограничения, не позволяющие принимать дорогостоящие решения (усложнение формы канала — переменность профиля, множественность изгибов и смешиваний, использование крайне сложных алгоритмов управления электрическими и тепловыми полями, что приведет к резкому удорожанию программно-математического обеспечения (ПМО), применение сверхдорогих реагентов при проведении химической реакции и т. п.).

Таким образом, практическая задача моделирования — поиск условно оптимальных конструктивных решений и режимов анализа для микрофлюидной аналитической системы при чип-реализации электрофореза.

БАЗОВЫЕ ПОЛОЖЕНИЯ

Основные базовые положения, в соответствии

с которыми осуществляется построение математических моделей процессов конвективно-диффузионного массопереноса в микроканалах чипа, достаточно традиционны. Среди признанных положений: 1) гипотеза сплошной среды, приемлемая для всех объектов, кроме самых крупных (например, клеток); 2) общая схема процесса — длинный тонкий канал прямоугольного или тра-пециидального сечения и короткая "пробка" анализируемой смеси; 3) система уравнений неразрывности и Навье—Стокса для описания распределения концентраций компонент; 4) система уравнений Пуассона—Больцмана для моделирования электрического потенциала, температурного поля и конвективного скоростного профиля.

Следует подчеркнуть, что распределение конвективной скорости по сечению канала в общем случае может отличаться от параболического, поскольку величины, входящие в соответствующее уравнение системы Пуассона—Больцмана (параметры канала, концентрация и вязкость буфера, температура и т. д.) могут быть различны.

Кроме того, одним из базовых положений, определяющим выбор метода решения, является отказ от получения точных численных решений (нестационарное пространственное распределение концентраций компонент) и использование интегральных характеристик (аналитических пиков). Такой подход допустим, поскольку детектор осуществляет суммирование откликов (по существу, концентраций) по всей глубине микроканала и иногда и по его ширине. Таким образом, анализируемыми характеристиками массопереноса будут параметры (центр тяжести, дисперсия, коэффициент асимметрии и т. д.) аналитических пиков компонент.

СПОСОБЫ УПРАВЛЕНИЯ АНАЛИЗОМ

Более подробно пять основных способов управления микропотоками вещества рассмотрены в работе [13], поэтому следует ограничиться их перечислением: 1) выбор топологии (геометрии) микроканала; 2) подбор управляющих напряжений для задания продольного (аксиального) электрического поля; 3) введение микрореактора с последующей регистрацией продукта; 4) управление состоянием поверхности микроканала и 5) создание управляемого температурного поля. Далее будут проиллюстрированы основные результаты моделирования указанных способов управления микроанализом при чип-реализации электрофореза.

двойной крест с варьированием ширины канала и длины пробки анализируемого вещества. Поскольку типичными размерами микроканала и пробы являются: длина Ь = 1-3 см, ширина 2к от 60 до 100 мкм, глубина Ь от 10 до 20 мкм (но не более полуширины), а объему пробы V от 50 до 300 пл соответствует длина пробки А от 80 до 500 мкм, очевидно, что базовая схема движения короткой пробки в длинном тонком канале подтверждается.

В первом приближении микроканал интерпретировался совокупностью плоских щелей. Параметры аналитического пика для базовой 2Б-мо-дели плоской щели вычислялись на основе метода моментов [6, 12]. Было показано [15, 16], что дифференциальное уравнение нестационарного пространственного распределения концентраций компоненты смеси (уравнение Навье—Стокса) может быть преобразовано для моментов к форме

дУп - в* д Уп

дt

дг 2

= пи*(1 -)уп-1 + п(п -1)В*у п = 0,1,2,3,...

п-2 з

(1)

где и , В — это нормированные на полуширину микроканала максимальная скорость конвективного движения и коэффициент диффузии соответственно; г — относительная ширина канала, отсчитываемая от оси [6].

Граничные условия для моментов любого порядка являются условиями второго рода, связанными с непроницаемостью стенки канала и осевой симметрией:

дУп (0, 0 = 0 дУп (1, 0 = 0

дг

дг

Начальные условия определяются исходным распределением вещества, которое мы полагаем соответствующим пробке относительной длины 2А с равномерным распределением вещества по сечению. Решение предполагает разложение конвективного скоростного профиля в ряд Фурье четным образом, т. е. вычисление коэффициентов В:

1 - = В0 + ^ В3 оо8(п ). Отметим, что

з =1

гипотеза о параболическом конвективном профиле была отвергнута в пользу более общей гипотезы — профиля с параметром клиновидности т. Тогда

Выбор топологии (геометрии) микроканала

В рамках указанного способа управления рассматривается оптимизация конструкции только для простейших топологий, а именно простой и

у1( г, 0 = и * В01 + ^

и В3

3 =1

В (П)'

-008(3) ,

и его среднее значение (среднее по сечению мик-

роканала положение центра тяжести пика) есть и*В0Х, где В0 = т/(т+1).

Средняя дисперсия пика есть

И=

2 + Ре2

В

1 2

=1(1

О X +

(

а!

3

3ре2'

2

В12

(2)

Жп )4

V У * *

где число Пекле определяется как Ре = и /О . Выражение (2) аналогично так называемой Тейлоровской дисперсии. Очевидно, что если максимальная скорость и коэффициент диффузии являются характеристиками анализируемого вещества, то управление положением центра тяжести и дисперсией пика осуществляется 3 параметрами: полушириной И (влияет на и , В*), объемом пробы (длина пробки А) и параметром клиновидности профиля т (последнее регулируется с помощью других методов управления).

На основе базовой модели плоской щели были решены некоторые сопровождающие задачи: канал трапециидального сечения был отмоделирован как совокупность плоских слоев разной ширины, а также проведено моделирование застойных зон (пристеночного обедненного веществом слоя) в соответствии с моделью неравномерной загрузки пробки [6]. Последние модели не подтверждаются экспериментами, поскольку выравнивание количества вещества по сечениям канала должно происходить, согласно расчетам, в течение нескольких секунд. Кроме того, эксперименты не подтверждают неравномерность загрузки. Таким образом, моделирование застойных зон должно осуществляться при воображаемой полупроницаемой стенке. Переход к трехмерной задаче как совокупности плоских щелей приводит к заниженным оценкам дисперсии, поскольку модель плоской щели не учитывает массообмен между слоями.

Расчет управляющих потенциалов

(продольного электрического поля)

На основе системы уравнений Кирхгофа и Ома ранее были выполнены [3] расчеты схемы управляющих потенциалов (топология: крест и двойной крест) с учетом нестационарности баланса потенциалов из-за изменения удельного сопротивления микроканала при движении смеси. При этом каналы рассматривались как делители напряжений. Исследовались два режима: ввод пробы и собственно анализ (разделение компонент). Для указанных двух режимов требовалось решить задачу оптимизации по критерию: минимум "ухода" вещества в боковые каналы при минимуме необходимых переключений управляющих потенциалов.

Основной принцип подбора потенциалов — минимизация "засасывания" буфера в пробу и "рассасывания" пробы в боковые каналы — обеспечивается подачей блокирующих потенциалов, примерно равных потенциалам в узлах (точках пересечения каналов). В расчетах было учтено переменное сечение микроканала вблизи от управляющих электродов. Была предложена удачная геометрическая аппроксимация переменного сечения канала, на основе которой были оценены параметры делителя напряжений. Подробный расчет геометрии каналов и баланса потенциалов также представлен в [3]. Так, при напряжении управляющих потенциалов для анализа и = 0, и2 = = 1500 В экспериментально найденные величины запирающих потенциалов будут и3 = = 5 В, и4 = 25 В; их расчетные значения [3] оказались соответственно и3 = 7.5 В и и4 = 29.7 В. Полученные оценки имели расхождение с экспериментально найденным оптимальным режимом управления не более 5 В.

Введение микрореактора

Основаниями предполагать эффективным введение микрореактора для последующей регистрации продукта можно считать: а) большую чувствительность детектирования к продукту, а не исходному реагенту; б) иные по сравнению с реагентом конвективные скорости и коэффициенты диффузии, что может способствовать изменению параметров пика.

Базовое уравнение Навье—Стокса для реагента и продукта реакции преобразовано с учетом химической реакции 1-го порядка. Очевидно, что указанное уравнение дополнительно содержит константу скорости химической реакции к, которая может быть управляющим параметром процесса анализа. Данную модификацию изначально предложил В.П. Андреев [17], использовав следующие допущения: 1) капилляр — круглый, субмиллиметрового диаметра; 2) профиль — конвективный, по скорости — параболический; 3) равенство конвективных скоростей продукта и реагента. Эти допущения неприемлемы в нашем случае, хотя сам предложенный подход позволил модифицировать уравнение (1) к уравнению моментов реагента и продукта.

Эти дифференциальные уравнения для реагента и продукта принимают вид соответственно

А - В дА д 1

д 2

= ПЩ (1 - Ыт1 )Ап-\ + п(п - 1) Ап-2 - А

(3)

+

У д

- В-

д2 У2

дг2

= пи

:(1 - кГ

1 )Уп2-1 + п(п - 1)В Уп2-2 + ^У,1

Здесь верхние индексы 1 соответствуют характеристикам реагента, 2 — продукта реакции первого порядка.

Граничными условиями являются условия второго рода; начальные условия соответствуют случаю равномерно заполненной "пробки" реагента относительной длины А и отсутствию продукта реакции. Видно, что первое уравнение в системе

1 1 —кг

(з) заменой момента уп ^ уп •е позволяет получить уравнение для "амплитуд" моментов реагента с решением, соответствующим отсутствию реакционного члена, т. е. уравнению (1).

Здесь подставлены нормированные скорости

(и) коэффициенты диффузии (Д) и параметры клиновидности (тг). Момент нулевого порядка (количество вещества) для реагента имеет вид

у1 = ехр(-к^, а для продукта реакции —

у02 = 1 - ехр(-к^. Так как необходимое для детектирования количество продукта ограничено снизу, то произведение к также должно быть ограничено снизу.

Средние по сечению положения центров тяжести пиков реагента и продукта определяются соответственно как

С1 = иВ1 и С2 =

= t(и 2 Вд - и1 В^ • ехр(-^)) /(1 - ехр(-^)) -

- (и 2 Вд - и1 В^)/ к.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Коэффициенты разложений конвективных профилей реагента и продукта реакции в ряды

Фурье по косинусам В3 (/ = 1,2) могут быть

различны. В случае совпадения конвективных профилей (то есть и1 = и2 и т1 = т2) будет равенство центров тяжести пиков реагента и продукта.

Выражение для средней дисперсии пика продукта реакции вычисляется по аналогии с теорией вероятностей:

22 < а >=< у22 >

/У02 -(< (у2)2 > /(У02)2).

иувилля, то момент первого порядка представлялся в форме

У2 (г, 0 = Щ0 ^) + ^Щ ^) оо8(к/г),

3

а среднее значение квадрата этого момента будет Wo2(t)+1 ).

23

Построение последовательности функций Щ(0 проводится по схеме:

а) ) = С 2(1 - ехр(-кО);

б) построение функций Щ производится с введением вспомогательных коэффициентов:

если к Ф В2п232, то

Жз (t) = у ехр (-В2п 2321) +

и 2 В

3

В П 2 32

+ а ехр (-к )+ в ехр (-(к + В1п 232)t), при этом

и1 В13

а = -

^ " и 2 В 2

ВП 3_

В2П232 - к

кщ В1. в = - 1 3

1

В1п232 (В2 - В1)п232 - к

и

7 =

и 2 В2 В2П232

-а-в;

если к = В2п2 32, то

и2В2

Щ (t) = 7ехр(-В2П3¿t) + - 2 2 .

В2П 3

и1 В1 2

+ (к----и 2 В 2 )t ехр(-к*) +

Вп 32

+ в ехр(-(к + В1п2 3 2)0

То есть требуется вычисление среднего момента второго порядка и среднего квадрата момента первого порядка. Поскольку уравнения решались методом разложения по собственным функциям соответствующей однородной задачи Штурма— Лиувилля, то момент первого порядка представлял-

при в = -

ки1 В)

1

Вп232 (В2 - В1П232 - к

и

У =

и2В2 _2 3

В2п231

-в .

Средний момент второго порядка для продукта реакции выражается громоздкой формулой:

< А >= 2и 2 Вд

2

2 X

и 2 В0 — - ЩВ0

1 1 - (1 + кХ)ехр(-кХ) и1 В1 - и2В0

2

к2

к

X -

1 - ехр(-кХ) к

„ 1 -ехр(-кХ) А2 „ . , чч . Ыч22-(2 + 2кХ + к2X2)ехр(-кХ) + 2В2 X - 2В2-—-- + — (1 - ехр(-кХ)) + (щ В1)2-1-2--——-- +

к 3 к2

1 - (1 + kX)ехр(-к*) Щ2 ^ (В1)2 1 - (1 + kX)exp(-kX) + 2 В,--1--> —-—----

1 к В^ П12 к

ки^у (В})2 2 у

П2 „г 4 ,-4

В1 1 П 1

1 - ехр(-к*) 1 - ехр(-(к + В1п212 )X)

к

22

к + В1п 1

У и 2 В1

и2В

2 - 1 X + у 1 - ехр(-В2П2/X) + ^ 1 - ехр(-кр + в 1- ехр(-(к + В1П1 У)

В2П212

В2П212

к + В1п 12

(4)

+

+

+

В случае совпадения конвективных профилей и коэффициентов диффузии продукта и реагента, дисперсии их пиков будут совпадать. Этот вывод также совпадает с выводом работы [17]. При этом управляющими параметрами являются как абсолютные значения, так и отношения одноименных величин: конвективных скоростей (и1/и2); параметров клиновидности (т^т2); коэффициентов диффузии продукта и реагента (В2/В1); а также константа скорости реакции к, длина (объем) пробки А и время анализа X.

Управление состоянием поверхности микроканала

Обоснование данного способа управления иллюстрируется следующей логической цепочкой:

2. Малая концентрация водоподобного буфера

1)

2)

3)

С0 (10 -10 М). Как следствие — относительная диэлектрическая проницаемость £= 80 (воде соответствует значение 81).

3. Величина дзета-потенциала <^0 лежит в пределах 10-100 мВ. В этой ситуации коэффициент в = (^0)/(ЯТ) < 1 и возможна линеаризация ~ (^ и Я обозначают соответственно постоянную Фарадея и универсальную газовую постоянную).

Тогда уравнение для скоростного профиля из системы Пуассона—Больцмана примет вид:

и (2) =

2Е^0££0 оЬ(Я)

П

1-

сЬ(Яг) еЬ(Я)

(5)

1) величина дзета-потенциала (потенциал на стенке микроканала) как граничное условие в уравнении из системы Пуассона—Больцмана влияет на профиль поперечного электростатического поля;

2) происходит трансформация профиля конвективной скорости (как средней скорости, так и параметра клиновидности);

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

3) меняются параметры (центр тяжести и дисперсия пиков) и, как следствие, разрешение соседних пиков.

В работе [6] представлено решение системы уравнений Пуассона—Больцмана при следующих допущениях.

1. Температура в каждом сечении канала одинакова Т = 300 К.

где в явном виде /(££0ЯТ), Е — на-

пряженность продольного управляющего электрического поля, П — динамический коэффициент вязкости, Е0 — электрическая постоянная (£0 =

107

=-- = 8.842 • 10-12 ф/

м). Решение уравнений

4пс

и обоснование соответствующих граничных условий также приведены в [6].

Аппроксимация (5) степенным профилем и(2) = итах (1 - 2т) проводится в соответствии с несколькими возможными критериями оптимизации, связанными с близостью двух профилей скорости. Предложенные критерии, позволяющие связать параметр Я, определяющий толщину двойного электрического слоя, и параметр клиновид-

ности конвективного профиля т обсуждены в работе [13].

Важнейшей особенностью аппроксимации (5) является то, что показатель клиновидности профиля m не зависит от величины дзета-потенциала, коэффициента динамической вязкости и от напряженности электрического поля. От этих параметров зависит абсолютная величина скорости, а не ее распределение по сечению канала. Однако при неприменимости линейной аппроксимации гиперболического синуса, то есть когда в— 1, что осуществляется при превышении дзета-потенциалом некоторого порога, не только средняя скорость, но и параметр клиновидности т будут зависеть и от величины дзета-потенциала, и от коэффициента динамической вязкости, и от напряженности электрического поля.

Создание управляемого температурного поля

Концепция данного нового способа управления заявлена в работе [13]. Ее обоснование также связано с уравнениями Пуассона—Больцмана, поскольку даже в линейную аппроксимацию (5) входит в явном виде температура. Эта концепция управления может быть реализована, исходя из трех алгоритмов различного уровня сложности и требующих различного ПМО. Во всех случаях основным критерием оптимального выбора управляющего температурного поля является наилучшее разрешение аналитических пиков компонент.

Были предложены и охарактеризованы следующие три основных алгоритма [13].

Алгоритм 1. Варьирование постоянной во всем канале температуры Т или применение микрогабаритных нагревателей для формирования продольных градиентов теплового поля.

Алгоритм 2. Задание начальной температуры и граничных условий на стенках канала с изменением Т естественным путем (действие управляющего электрического поля).

Алгоритм 3. Задание продольного градиента температуры для оптимального управления процессом разделения компонент с перспективой внедрения следящей системы управления температурой, "отслеживающей" прохождение быстрых компонент с более узкими пиками и меняющей температуру (скоростной профиль, коэффициент диффузии и, следовательно, закон формирования пика более медленных компонент) с целью сужения пика последних и достижения лучшего разрешения.

ОБРАБОТКА ИНФОРМАТИВНЫХ СИГНАЛОВ

Поскольку решение задачи поиска эффективного алгоритма обработки информативных сигналов влияет на качество анализа, то опосредовано ре-

зультаты данного направления работ представлены в [1-5, 7, 9, 10]. Однако подход к поиску алгоритма обработки информативных сигналов и его применение более детально описаны в работах [8,

11, 14].

Основания для применения специальных дополнительных алгоритмов обработки аналитических сигналов (пиков) в приборах электрофореза на микрочипах следующие: 1) переход к пиколит-ровым диапазонам анализируемых объемов (т.е. уменьшение анализируемых объемов на 3-5 порядков по сравнению с классическим капиллярным электрофорезом) приводит к уменьшению амплитуды информативного сигнала (уменьшению отношения сигнал/шум); 2) сокращение времени анализа (от 3 до 10 раз) приводит к сужению пиков. Таким образом, традиционная модель аналитического пика — кривая Гаусса оказывается неадекватна. Также следует заметить, что именно на модель гауссианы ориентированы традиционные критерии разрешения пиков.

Очевидно, что изменение отношения сигнал/шум существенно осложнит на практике алгоритм разрешения альтернативы "пик/шум": огрубление алгоритма (увеличение пикового порога) приведет к возрастанию ошибок первого рода — пропуск истинных пиков компонент, неизменность традиционных порогов приведет к резкому возрастанию ошибок второго рода — принятию шумов в качестве пиков (ложных). Кроме того, большие шумы существенно осложнят задачи выделения базовой линии, определение границ и центра тяжести пика. Сокращение времени анализа приведет к необходимости замены адекватных моделей пика на треугольную (кусочно-линейную) и параболическую [14]. Также заметим, что для узких островершинных (треугольных) пиков

8400

8000

7600

7200

6800

_|_._I_._I_._1_

1000 1050 1100 1150 1200 Такт

1250

Фрагмент электрофореграммы разделения смеси олигонуклеотидов Аш1х (10, 20, 30, 40, 50, 75 и 100 оснований). Один такт составляет 0.2 с

серьезной проблемой будет расчет площади пика.

Типичный фрагмент электрофореграммы на микрочипе представлен на рисунке. Очевидно наличие на рисунке нескольких ложных пиков (в области 1130 и 1200 тактов). Также следует отметить, что гауссова модель может быть применена лишь к последнему (правому пику). Пики 2-4 (отсчет слева), наверное, следует аппроксимировать треугольником, остальные — параболой.

Предложенный в работе [14] и подтвержденный результатами [8, 11] подход к обработке аналитических сигналов базировался на использовании цифрового фильтра (тангенсный фильтр Бат-терворта второго порядка). Алгоритм включал две стадии: а) префильтрация цифровым фильтром нижних частот (ФНЧ) с частотами среза 0.03-0.12 от частоты дискретизации (съема данных); б) обратный пересчет амплитуды пика, ширины, площади отфильтрованных сигналов к исходным по соответствующим таблицам.

ЗАКЛЮЧЕНИЕ

Благодаря использованию метода математического моделирования были получены следующие основные результаты.

1. Анализ, по крайней мере на уровне концепций, всех пяти предложенных способов управления процессом анализа смесей на микрочипе показал их обоснованность и эффективность.

2. Предложенная базовая модель плоской щели позволяет использовать метод моментов для оценки параметров аналитических пиков компонент и тем самым выбрать управление анализом, варьируя ширину микроканала, объем пробы и точку детектирования (или время анализа). В свою очередь, модификация модели на трехмерный случай позволит решать аналогичную задачу применительно к каналам прямоугольного и трапеции-дального сечения. В последних в качестве дополнительных управляющих параметров могут выступать глубина канала и угол наклона стенок.

3. На основе метода моментов в рамках базовой модели плоской щели получено явное выражение (3) для центра тяжести и средней дисперсии пиков продукта и реагента. Эти выражения содержат в качестве управляющих параметров нормированные значения максимальных конвективных скоростей, коэффициентов диффузии и константу скорости химической реакции первого порядка к. Выражение для дисперсии продукта имеет весьма громоздкий вид (4). В случае равенства скоростных профилей и коэффициентов диффузии продукта и реагента характеристики их аналитических пиков будут одинаковыми.

4. Модификация поверхности микроканала (величина дзета-потенциала) также позволяет управ-

лять анализом. Однако методика управления характеризуется следующими особенностями. Если величина дзета-потенциала не достигла некоторого минимального порога (порога применимости линейных аппроксимаций уравнений системы Пу-ассона—Больцмана), то: а) ее изменение не влияет на степень клиновидности профиля т, однако б) это изменение влияет на масштаб конвективной скорости, который прямо пропорционален величине дзета-потенциала (5).

5. Применимость температурных полей для управления анализом на микрочипе исследовалась ранее в работе [13]. В указанной работе оценивалось влияние вариации температуры только на толщину двойного электрического слоя (параметр Я) или на параметр клиновидности скоростного профиля. Выявлена низкая чувствительность т = т(Т). Однако способ управления может учитывать действие других тепловых эффектов — изменение динамической вязкости, коэффициента диффузии и т.д.

6. Фильтрация информативного сигнала, представленного на рисунке, показала наибольшую эффективность применения ФНЧ с относительными частотами среза 0.06-0.09: выбор большей частоты среза фильтра ухудшает качество подавления шумов, использование меньших частот приводит к очень сильным деформациям формы пиков — резкому уменьшению амплитуд для треугольных пиков и существенному расширению и асиммет-ризации параболических пиков.

Благодарности. Работа выполнена при финансовой поддержке в рамках гранта РФФИ № 03-0139003 ОБЕК_а "Теоретические и экспериментальные исследования явлений переноса и взаимодействия биологических объектов в микрофлюидных устройствах".

СПИСОК ЛИТЕРАТУРЫ

1. Беленький Б.Г. и др. Микрофлюидная аналитическая система с детектором лазер-индуцированной флуоресценции // Аллергология и иммунология. 2000. Т. 1, № 3. С. 101-102.

2. Евстрапов А.А. и др. Аналитические системы на основе микрочиповых технологий и методов капиллярного электрофореза // Тезисы докладов 1-й Всероссийской конференции "Аналитические приборы", 18-21 июня 2002, Санкт-Петербург, Россия. С. 40-41.

3. Евстрапов А.А. и др. Микроаналитические системы: портативные и лабораторные экспресс-анализаторы и биосенсоры на основе чип-реализации одноканального и мультиканального приборов капиллярного электрофореза и иммуноэлектрофореза. Отчет НИР-И, 2002,

СПб., ИАнП РАН, ГР № 01.20.02 04223. 65 с.

4. Евстрапов А.А. и др. Микрофлюидные аналитические системы на основе электрофоретиче-ских методов анализа // Новости науки и техники, Серия "Медицина, аллергия, астма и клиническая иммунология". 2001. № 1. С. 190-193.

5. Евстрапов А.А. и др. Экспресс-анализ олиго-нуклеотидов на планарном микрофлюидном чипе // Журнал аналитической химии. 2004. Т. 59, № 6. С. 587-594.

6. Буляница А.Л., Евстрапов А.А., Рудницкая Г.Е. Метод моментов при расчете параметров каналов в микроразмерных системах // Научное приборостроение. 2003. Т. 13, № 4. С. 28-40.

7. Евстрапов А.А. и др. Микрофлюидные аналитические системы на основе методов капиллярного электрофореза и микрочиповых технологий // Аллергия, астма и клиническая иммунология. 2003. Т. 7, № 9. С. 205-211.

8. Evstrapov A.A. et al. Estimation of data reduction effects on informative parameters of analytical signals, their consideration and compensation // Abstracts program 3rd Int. Symposium on Separations in BioSciencies SBS2003 "100 Years of Chromatography", 13-18 мая, 2003, Москва, Россия. С. 184.

9. Evstrapov A.A. et al. Devices and analytical systems on the basis of hydraulik microchips // Тезисы докладов Межд. конф. "Instrumentation in ecology and human safety", 4-6 ноября 2002, Санкт-Петербург, Россия. С. 45.

10. Евстрапов А.А. и др. Применение микрофлюидных аналитических систем в приборных комплексах биологического и экологического экспресс-анализа // Сборник материалов III Международного симпозиума "Аэрокосмические приборные технологии" (АПТ'04), 0204 июня 2004, Санкт-Петербург, Россия. С.316-317.

11. Буляница А.Л. и др. Оптимизация алгоритмов обработки электрофореграмм при анализе

биологических веществ на микрочипе // Каталог рефератов и статей Межд. форума "Аналитика и аналитики", 2-6 июня, 2003, Воронеж, Россия. Т. II. С. 592.

12. Bulianitsa A.L., Evstrapov A.A., Rudnitskaya G.E., Kurochkin V.E. Method of moments applied to assess analyte velocity profiles in electrophoretic separation of DNA fragments // Book of Abstract 2nd Black Sea Conference on Analytical Chemistry, Workshop on 1-st Marmara Analytical Chemistry, 14-17 Sept., 2003, Istambul, Turkey. Abstr. N P081, P. 135.

13. Буляница А.Л. Управление микропотоками вещества в канале микрофлюидного чипа с помощью регулируемых тепловых полей // Научное приборостроение. 2005. Т. 15, № 1. С. 56-61.

14. Евстрапов А.А. и др. Особенности применения алгоритмов цифровой фильтрации элек-трофореграмм при анализе веществ на микрочипе // Научное приборостроение. 2003. Т. 13, № 2. С. 57-63.

15. Туницкий Н.Н., Каминский В.А., Тимашев С.Ф. Методы физико-химической кинетики. М.: Химия, 1972. 198 с.

16. Андреев В.П., Брезгун А.А., Павленко И.В. Математическое моделирование процессов мас-сопереноса в противоточной распределительной хроматографии // Журнал физической химии. 1988. Т. LXII, N 9. C. 2448-2454.

17. Андреев В.П., Хидекель М.И. К математической модели проточно-инжекционного анализа. Часть 1: Химическая реакция первого порядка в прямой трубке. Препринт № 55. СПб.: ИАнП РАН, 1993. 21 с.

Институт аналитического приборостроения РАН,

Санкт-Петербург

Материал поступил в редакцию 17.03.2005.

MODELING PROCESSES OF HEAT (AND MASS) TRANSFER IN MICROCHANNELS FOR ELECTROPHORESIS ON A CHIP

A. L. Bulyanitsa

Institute for Analytical Instrumentation RAS, Saint-Petersburg

The paper analyzes the role of the stage of mathematical modeling of heat and mass transfer in the development of a microfluidic analytical system and performing chemical and biological analyses in such a system. Among numerous modeling problems, the main concern was with effective methods of analysis control for multicomponent mixtures of chemical and biological substances in on-chip electrophoresis. Also covered are related issues such as optimization of various analysis stages including signal processing.

i Надоели баннеры? Вы всегда можете отключить рекламу.