УДК 531/534:[57+61]
КОМПАРТМЕНТАЛЬНОЕ МОДЕЛИРОВАНИЕ ТРАНСПОРТНЫХ ПРОЦЕССОВ В КОРНЕ РАСТЕНИЯ, УЧИТЫВАЮЩЕЕ ПРИСУТСТВИЕ ПОГРАНИЧНЫХ СЛОЕВ
С.А. Логвенков1, 2, А.А. Штейн2
1 Кафедра высшей математики Национального исследовательского университета «Высшая школа экономики», Россия, 101000, Москва, ул. Мясницкая, 20, e-mail: [email protected]
2 Институт механики Московского государственного университета имени М.В. Ломоносова, Россия, 117899, Москва, Мичуринский проспект, 1, e-mail: [email protected]
Аннотация. На основе информации о структуре решения и асимптотических оценок в задаче о стационарном потоке через корень получена система приближенных алгебраических соотношений, аналогичная имеющимся компартментальным моделям. Предлагаемый метод имеет по сравнению с последними существенное преимущество, позволяя учитывать особенности анатомического строения корня и неоднородность распределения величин по его сечению. Это дает возможность сформулировать простые конечные зависимости, совпадающие с достаточной точностью с численным решением, получаемым в рамках континуальной модели. Применение предлагаемого подхода для решения конкретных задач оказывается более простым по сравнению как с численным решением, основывающимся на континуальной модели, так и с решением, получаемым асимптотическими методами.
Ключевые слова: поглощение воды растениями, корень растения, многофазные среды, клеточные системы, транспортные процессы, клеточные мембраны, компартментальные модели.
Введение
Изучение механизмов, участвующих в переносе веществ в корнях растений и обеспечивающих поддержание корневого давления для подъема воды
на необходимую высоту и снабжение растения питательными веществами, поступившими из почвы или окружающего раствора, оказывается невозможным без использования математических моделей.
По большей части для этой цели используются компартментальные модели, в которых различные цилиндрические слои корня рассматриваются как компартменты, отделенные друг от друга и от внешней среды мембранами, вообще говоря, с различными свойствами. При этом массообмен между компартментами описывается в рамках традиционных для неравновесной термодинамики линейных связей потоков с гидростатическими и осмотическими давлениями. Компартментальные модели, предложенные в работах [5-9], отличаясь друг от друга числом рассматриваемых компартментов и свойствами разделяющих их мембран, позволили находить хорошее совпадение с получаемыми в эксперименте зависимостями, в частности с зависимостью разности гидростатических давлений между ксилемой и окружающим раствором от потока через корень. Однако подобное рассмотрение остается в достаточной мере
© Логвенков С.А., Штейн А.А., 2013
Логвенков Сергей Алексеевич, к.ф.-м.н., доцент кафедры высшей математики, Москва
Штейн Александр Александрович, к.ф.-м.н., ведущий научный сотрудник Института механики, Москва
эмпирическим и оставляет в стороне основные особенности функционирования корня как распределенной механической системы. Вследствие этого модели такого рода неизбежно остаются привязанными к узкому кругу экспериментов и частных гипотез и мало пригодны для включения в полную гидромеханическую модель транспорта в растениях.
Недостатки компартментального моделирования привели к необходимости рассматривать более полные континуальные модели, учитывающие реальное анатомическое строение корня и включающие в себя основные физические механизмы, влияющие на перенос веществ [1, 2, 10]. Однако решение задач на основе континуальных моделей существенно сложнее, чем на основе компартментального моделирования. Последовательное использование асимптотических методов [4] оказалось громоздким и реализуемо лишь в очень простых случаях, охватывающих малую часть подлежащих решению задач.
В работе [1] авторами была получена компартментальная модель, основывающаяся не на схематическом представлении корня совокупностью мембран, а на осреднении континуальных уравнений раздельно в двух областях, граница между которыми задается областью поясков Каспари. Однако в каждой из этих областей осреднение по пространственной координате выполнялось стандартным методом, дающим хорошие результаты для гладких функций. Между тем, как показало численное решение [1, 10], в каждой из этих областей имеются пограничные слои, в которых осредняемые функции резко меняются, что может внести существенную погрешность в результаты, полученные методом осреднения, использованным в исследовании [1].
В предлагаемой работе на основе асимптотических оценок выведены уточненные осредненные соотношения, учитывающие информацию о структуре решения в задаче о стационарном потоке через корень, полученную на основе континуальной модели [1, 10]. Получены приближенные формулы, описывающие зависимость потока жидкости через корень от разности приложенных давлений и концентрации низкомолекулярного компонента от потока жидкости, а также распределение давлений, скоростей и концентраций по сечению корня.
Задача о стационарном потоке веществ через корень растения.
Континуальная модель
Методология континуального моделирования процессов массопереноса в корнях впервые предложена авторами в работах [1, 10]. В этих работах была разработана базовая континуальная модель, учитывающая реальное анатомическое строение корня, включающая в себя основные физические механизмы, влияющие на перенос веществ и позволяющая оценивать участие различных известных и предполагаемых механизмов в процессе поглощения жидкости корнем. Ткань корня моделируется жидким континуумом, заполняющим пористый каркас и состоящим из двух фаз, соответствующих внутриклеточной (симпласт) и внеклеточной (апопласт) жидкостям. Система уравнений и граничных условий подробно описана, например, в работе [1] и приводится здесь для случая стационарного осесимметричного течения без каких-либо обоснований и объяснений.
В предположении осевой симметрии корня решение ищется в области Г ^ т < T1, где г0 и г1 - координаты раздела корня с сосудами ксилемы с окружающей средой. Стационарные распределенные уравнения имеют следующий вид:
1 йухт т йт
т йт
(2)
(3)
(4)
(5)
р 1 й (1 -$ с2 ^т Р 2 ,
т йт
(6)
Здесь индексы 1 и 2 соответствуют внеклеточной и внутриклеточной фазам соответственно; р1 и р2 - объемные плотности фаз; v1 и у2 - скорости жидкости в апопласте и симпласте; У - межфазный поток жидкости; Ьр - объемный коэффициент гидравлической проводимости ткани; р1 и р2 - гидростатические давления в фазах;
Я = р0 ■Я0/д0 (Я - универсальная газовая постоянная, д0 - молярная масса
растворенного низкомолекулярного компонента, р0 - плотность воды); Т - абсолютная температура; с1 и с2 - массовые концентрации обобщенного низкомолекулярного компонента в каждой фазе; в1 и в2 - проводимости апопласта и симпласта; ^ - распределенный коэффициент отражения симпласта; к - межфазный активный поток вещества; X - проницаемость клеточных мембран по отношению к растворенному веществу; Д и ^2 - коэффициенты диффузии растворенных в апопласте и симпласте веществ. Коэффициенты системы (1)-(6) в дальнейшем предполагаются постоянными за исключением X и к, которые считаются кусочнопостоянными с разрывом в месте расположения поясков Каспари, соответствующем т = тс (т0 < тс < т1): при т > тс X = Х+, к = к+, при т < тс X = Х~, к = к .
Предполагается, что концентрация внешнего раствора, а также давления во внешнем растворе и сосудах ксилемы являются известными фиксированными величинами, в то время как концентрация в сосудах ксилемы подлежит определению наряду со скоростью поступления жидкости в нее.
Система граничных условий в пренебрежении потоком воды и растворенных веществ через внешнюю и внутреннюю границы симпласта по сравнению с потоками между симпластом и апопластом следующая (верхние индексы + и - соответствуют значениям функций по разные стороны от границы, определяющей пояски Каспари
т = тс):
(7)
т = тс: V-
йт йт
Приближенные формулы
Решение поставленной выше задачи при значениях коэффициентов, приведенных ниже, показало, что величины с2, р2, и ^2 испытывают
значительные изменения в пограничных слоях, примыкающих к внешней и внутренней границам поперечного сечения корня, а также слева и справа от поясков Каспари. Толщины пограничных слоев различаются для различных функций. Так, толщины пограничных слоев для концентрации в симпласте в зависимости от их локализации
имеют величины порядка Б2/V* и ^р2Б2/X* (V* и X* - характерные значения скорости и проницаемости клеточных мембран по отношению к растворенному веществу), а для давления в симпласте и скоростей - порядка ■\1р2$2/ £р [3].
При выбранных значениях коэффициентов функции с1 и р1 не испытывают резких изменений типа пограничного слоя. В области сечения, исключающей пограничные слои, величины , ^2, с1 и с2 меняются слабо. К этим заключениям можно прийти и
на основании более подробного асимптотического анализа. Перечисленные факты отражены на рис. 1. Заштрихованные области включают пограничные слои для всех величин, а буквенные обозначения соответствуют значениям функций на границе с пограничным слоем или значениям во всей области для слабо меняющихся величин.
Проинтегрируем уравнения (3) и (4) по интервалам, получаемым исключением из отрезков [т0; тс ] и [тс; т1] заштрихованных на рис. 1 областей, содержащих
пограничные слои. Учитывая, что величины и =-т^, и2 = -тv2 и с2 в каждой из областей практически постоянны, а также пренебрегая толщиной переходных слоев, получим приближенные равенства
В этих областях уравнение (1) принимает простой вид р2 - р1 + ЯТ(с1 - с2) = 0. Тогда значения величин на границах переходных областей оказываются связанными между собой соотношениями
(8)
Р2- Рх + ЯТ (сх- с2) = 0 Р* - Рі* + ЯТ(сх - С2 ) = 0
р** - р** + ЯТ(с° - С2+ ) = 0,
Р2+- Р° + ЯТ(С0 - С2+ ) = 0-
(9)
/г, * Р2 Рг С 2 /т2 ★★ 4- Р2 Рг + с2
Сх * р 1 р 1 П>, с0 ** + Р\ Р\ ГУ . 4*
=
^0
г = г
г = г
Г = г
Рис. 1. Расположение пограничных слоев для концентраций, давлений и скоростей (заштрихованные области); буквенные обозначения дают значения соответствующих функций на границе с пограничным слоем или значения во всей области для слабо меняющихся величин
Значения числовых параметров: Ьр = 4,2-10 6 кг-м 3-с 1-Па 1; жех1 = 0,15 МПа; ? = 0,89; г = 0,5 мм; г = 0,1 мм; гс = 0,35 мм; д/Ро = 0,07; р2/Ро = 0,2; к+ = 6,2 -10-5 кг-м~3-с-1; к- = -6,2 -10-9 кг-м-3-с-1; Х+ = 1,4 -10~2
X- = 2,1 -10-3 кг-м-3-с-1; Д = 3 -10-10 м2-с-1; Д = 5 -10-15 м2-с-1; р1 = в2 = 4 -10-18 м2-Па-1-с-1; ЯТ = 3,3 -107 Па.
Аналогичные рассуждения, проведенные для уравнений (5) и (6), дают следующие зависимости:
-3 -1
кг-м -с ;
10-18 м2-Па-1-с-1;
С2 °х X-,
к-
-
+ к+
с2 - с0 = —
2 0 Х+
(10)
Получим соотношение, связывающее значения давлений и концентраций в симпласте по разные стороны от переходных слоев, примыкающих слева и справа к пояскам Каспари при г = гс. Учитывая, что в этих слоях давление и концентрация
веществ в симпласте испытывают значительные изменения, в уравнении (4) можно пренебречь величиной скорости по сравнения со значениями производных давления и концентрации. Тогда после интегрирования уравнения (4) по переходной области получим
р2* - р2 = ^ят (с2+- с-).
(11)
Сложение уравнений (1) и (2) приводит к условию постоянства полного радиального потока жидкости. Отсюда, с учетом граничных условий (7), следует, что
РА+ +Р2и2+ = РЛ +Р2и- = РМг0)г0 = Р 2 У2(гс Ус .
(12)
Сложим уравнения (5) и (6). Интегрирование полученного уравнения на отрезке [г0 ; гс ] с учетом граничных условий (7) и уравнения (12) дает
(х - (1 - Ос2(гс )) () + Р2и2- ) = Р2П2гс ) .
(13)
добавить два уравнения, связывающие с2(гс) и —с другими неизвестными.
Для получения полной системы алгебраических уравнений необходимо
йс2( гс) йт
Рассмотрим асимптотику решения уравнения (6) в окрестности поясков Каспари. Учитывая сильное изменение концентрации в симпласте в этой области, можно пренебречь конвективным потоком растворенных веществ по сравнению с диффузионным. Чтобы проанализировать поведение решения в окрестности пограничных слоев слева и справа от т = тс, выполним преобразование растяжения соответствующих областей. Введем в каждой из них мелкомасштабную внутреннюю
е+ г - Гс + 1
переменную £ =—, где 8 = —
Р 2 ^2
- малые параметры. Знаки + и - указывают
Х~
на то, что слева и справа от поясков Каспари используются различные новые переменные, а пограничные слои имеют разную толщину. Переходя к новой переменной, получим в пограничных слоях уравнение вида (знаки + и - для простоты опущены)
г12
ч ё£2 ' Гс + 8£ ё£
2 +_ 2
с 2 + сі + — = 0 .
В первом приближении по малому параметру распределение концентрации в окрестности точки т = тс слева и справа имеет соответственно следующий вид:
к- ^ -^ с2 = сх + —+ Л~-ет + В • е т,
X
к+ ^ Х с2 = с0 + —+ Л+-еТ1 + В+ • е т .
2 0 Х+
Выполним сращивание решений в пограничных слоях с внешним решением. Будем использовать следующие условия сращивания [3]: Нш с2 = с- и Нш с2 = с+ .
к к+
Из этих условий получим, что Л+ = В ~ = 0, сх +—- = с-, с0 +-------= с+. Для определения
X X
оставшихся постоянных интегрирования Л~ и В+ будем использовать условия непрерывности концентрации с2 и диффузионного потока при т = тс:
сх + — + Л~= с0 + — + Л+, А- ~^ = -В\ — Х X- 0 Х+ А/ Р2А \Р2А
Отсюда следует
Теперь можно получить
Со Сх + Х+ X-
1 + ^ \Х+
X-
(14)
р2А (
йт
1+
С0 - сх +-----------------------
0 х Х+ X-
Уравнения (8)-(14) образуют полную систему уравнений для нахождения
Лс2{ гс)
+ + + ****** ± неизВестных щ, щ, Р2, Р2, Рг, Рі, Рі , с2, Сх , с2(тс) и
Расчеты и результаты
йт
В данном разделе в результате решения полученной системы уравнений будут найдены приближенные формулы, описывающие зависимость скорости течения от разности гидростатических давлений между окружающей средой и сосудами ксилемы, а также концентрации веществ в сосудах ксилемы от скорости потока на срезе корня.
Складывая четвертое уравнение (8) с третьим уравнением (9) и вычитая из этой суммы второе уравнение (8) и четвертое уравнение (9), получим равенство = -^.
Р 2 в1
Выполняя аналогичные действия с третьим уравнением (8), первым уравнением (9),
, . //"\\ и и
первым уравнением (8) и вторым уравнением (9), получим равенство = -^.
Р 2 в1
Учитывая полученные равенства и уравнение (12), находим, что и- = и+ и и- = и+. В дальнейшем при этих величинах будем опускать верхние индексы.
Введем обозначение и = — и1 +—и2, где р = р1 +р2. Тогда имеем выражения
и =■
Р1Р
Р2Р
^ - -и и и2 =------—------и . Складывая третье уравнение (8), четвертое
Р1Р1 + Р2Р2 Р1Р1 +Р2Р2
уравнение (8), первое уравнение (9) и уравнение (11), а затем вычитая из этой суммы четвертое уравнение (9), получаем
Ро - Рх =в“ в2
( (т л (т л >
1п и+ + 1п с и2 + (?
V V тс) V тоV
+ (^- (С2 - С2- ) + КТ (С0 - Сх )-
С учетом последних результатов это равенство можно переписать следующим образом:
- Рх = 1ПвТ|/То )р и + (? -1) ЯТ(с+ - с-)+ЯТ(Со - Сх).
Р1Р1+Р2Р2
(15)
Подставляя (14) в уравнение (13), найдем концентрацию низкомолекулярного компонента в сосудах ксилемы
с„ = •
и(1 -0 с0 + т- + т-\т- +Г
л/р2^-В2 (
с0 +
кг} X
+ г
у1р2Х ^2
Величину с- - с2 определим из уравнений (10):
(16)
с2 - с2 = с0 - сх +--------------------------------
2 2 0 х Х+
(17)
Алгебраические соотношения (15)—(17) дают интересующие нас зависимости концентрации в ксилеме сх и разности давлений между окружающей средой
и ксилемой Аp = р0 - рх от скорости и и, следовательно, от скорости вытекания
жидкости через срез корня 3 = 2uL / г12, где Ь - длина корня.
Численное решение и приближенные решения, полученные на основе двух компартментальных моделей (представленной соотношениями (15)-(17) и модели, сформулированной в работе [1]), сравнивались при одних и тех же значениях параметров. Результаты сравнения приведены на рис. 2 и 3. Видно, что полученные в настоящей работе соотношения с большей точностью совпадают с численным решением.
В частности, наклон асимптоты, к которой приближается график зависимости потока от разности гидростатических давлений на рис. 2, полученный с помощью компартментальной модели [1], заметно отличается от наклона асимптоты функции, полученной из точного решения, тогда как в модели, представленной в настоящей работе, такого отклонения не наблюдается. Указанное расхождение может привести к значительным ошибкам в определении гидравлической проводимости корней. Источник ошибки - в неучете присутствия пограничных слоев. Сравнение с экспериментом было выполнено для численного решения ранее [10] и здесь не проводится.
V -10-8, м/с
Ар, атм
Рис. 2. Зависимость скорости потока на срезе корня от разности гидростатических давлений между окружающей средой и ксилемой: 1 - расчет на основе континуальной модели; 2 - расчет на основе предлагаемой модели; 3 - расчет на основе компартментальной модели [1]
п, атм
v -10 8, м/с
Рис. 3. Зависимость осмотического давления в ксилеме от скорости потока на срезе корня: 1 - расчет на основе континуальной модели; 2 - расчет на основе предлагаемой модели; 3 - расчет на основе компартментальной модели [1]
Предлагаемый метод получения конечных соотношений, структурно сходных с традиционными компартментальными моделями, имеет по сравнению с последними существенное преимущество: он позволяет учитывать особенности анатомического строения корня и неоднородность распределения величин по сечению. Это дает возможность получать простые приближенные зависимости, совпадающие с достаточной точностью с численным решением. Использование таких зависимостей оказывается значительно проще использования распределенных моделей или применения асимптотических методов. В рамках предлагаемого подхода модель может быть в случае необходимости модифицирована с учетом дополнительных возможных механизмов.
Благодарности
Работа поддержана РФФИ (проект № 11-01-00774).
Список литературы
1. Логвенков С.А., Штейн А.А. Компартментальная модель поглощения воды корнями растений с учетом процессов на клеточном уровне // Российский журнал биомеханики. - 2008. - Т. 12, № 4. -C. 18-32.
2. Логвенков С.А., Штейн А.А. Влияние проницаемости поясков Каспари для воды и растворенных веществ на величину корневого давления. Математическое моделирование // Российский журнал биомеханики. - (Статья в настоящем выпуске).
3. Найфэ А. Введение в методы возмущений. - М.: Мир, 1984. - 535 с.
4. Юдина Е.Н. Асимптотическое исследование транспортных процессов в корне растения // Российский журнал биомеханики. - 2012. - Т. 16, № 2. - C. 59-67.
5. Dalton F.N., Raats P.A.C., Gardner W.R. Simultaneous uptake of water and solutes by plant roots // Agronomy Journal. - 1975. - Vol. 67. - P. 334-339.
6. Fiscus E.L. Determination of hydraulic and osmotic properties of soybean root system // Plant Physiology. -1977. - Vol. 59. - P. 1013-1020.
7. Ginsburg H. Model for iso-osmotic water flow in plant roots // Journal of Theoretical Biology. - 1971. -Vol. 32. - P. 147-158.
8. Miller D.M. Studies of root function in Zea mays. IV. Effects of applied pressure on the hydraulic conductivity and volume flow through excised root // Plant Physiology. - 1985. - Vol. 77. - P. 168-174.
9. Murphy R. Some compartmental models of the root: Steady-state behavior // Journal of Theoretical Biology. - 2000. - Vol. 207. - P. 557-576.
10. Stein A.A., Logvenkov S.A., Chalyuk A.T. Mathematical modelling of the plant root as a water-pumping cellular system // Mathematical modelling and computing in biology and medicine / ed. V. Capasso. -Bologna: Soc. Ed. Esculapio, 2003. - P. 206-212.
COMPARTMENT MODELLING OF TRANSPORT PROSESSES IN THE PLANT ROOT WITH ACCOUNT FOR THE PRESENCE OF BOUNDARY LAYERS
S.A. Logvenkov, A.A. Stein (Moscow, Russia)
Basing on the information about the structure of the solution and asymptotic estimates in the problem of steady flow across the root, a system of algebraic relations similar to the commonly used compartment models is obtained. As compared with these, the method proposed has an important advantage making it possible to take into account the characteristic features of the anatomical structure of the root and the non-uniformity of the parameter distribution over its cross-section. This enables us to formulate simple finite relationships fitting with sufficient accuracy with the numerical solution obtained within the framework of the continuum model. The application of the approach proposed to solving specific problems is simpler than both the numerical solution based on the continuum model and the solution obtained by asymptotic methods.
Key words: water absorption by plants, root plants, multiphase media, cellular systems, transport processes, cell membranes, compartmental models.
Получено 8 октября 2013