Научная статья на тему 'Исследование стационарных решений задачи динамики фитопланктона с учетом трансформации соединений фосфора, азота и кремния'

Исследование стационарных решений задачи динамики фитопланктона с учетом трансформации соединений фосфора, азота и кремния Текст научной статьи по специальности «Математика»

CC BY
110
19
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
ФИТОПЛАНКТОН / ФОСФОР / АЗОТ / КРЕМНИЙ / БИОГЕН / ХИМИКО-БИОЛОГИЧЕСКИЙ ИСТОЧНИК / УРАВНЕНИЕ КОНВЕКЦИИ-ДИФФУЗИИ-РЕАКЦИИ / ЗАДАЧА КОШИ ДЛЯ СИСТЕМЫ ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ / СТАЦИОНАРНЫЙ РЕЖИМ / PHYTOPLANKTON / PHOSPHORUS / NITROGEN / SILICON / BIOGEN / CHEMICAL-BIOLOGICAL SOURCE / CONVECTION-DIFFUSIONRESPONSE EQUATION / CAUCHY PROBLEM FOR SYSTEM OF ORDINARY DIFFERENTIAL EQUATIONS / STATIONARY MODE

Аннотация научной статьи по математике, автор научной работы — Белова Ю. В., Атаян А. М., Чистяков А. Е., Стражко А. В.

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

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

Похожие темы научных работ по математике , автор научной работы — Белова Ю. В., Атаян А. М., Чистяков А. Е., Стражко А. В.

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

Study on stationary solutions to the problem of phytoplankton dynamics considering transformation of phosphorus, nitrogen and silicon compounds

Introduction. The solution to the problem of transformation of phosphorus, nitrogen and silicon forms is studied. This problem arises under modeling phytoplankton dynamics in shallow-water bodies including the Azov Sea. The phytoplankton dynamics model is formulated as a boundary value problem for the system of diffusion-convection-response equations and takes into account the absorption and release of nutrients by phytoplankton, as well as the transition of nutrients from one compound to another. To calculate the initial conditions and parameters of the equations under which the steady-state regime occurs, the software is developed, which is based on the model describing changes in phytoplankton concentrations without considering current effects. This model is represented by a system of inhomogeneous differential equations. Based on the developed software, the initial conditions and parameters of the phytoplankton dynamics model in the Azov Sea are calculated experimentally.Materials and Methods. A 3D model of phytoplankton dynamics is considered taking into account the transformation of phosphorus, nitrogen and silicon compounds based on the system of nutrient transport equations. The case of a spatially uniform distribution of substances is considered to specify the parameters of the model at which the stationary modes occur. Because of simplification, a system of ordinary differential equations solved through the Runge-Kutta method is obtained.Research Results. The software is developed to specify the initial conditions and parameters of the phytoplankton dynamics model considering the transformation of phosphorus, nitrogen and silicon compounds. Several numerical experiments are performed under the assumption that the development of phytoplankton is limited by a single biogenic substance. As a result of the computational experiment, it can be seen that with the obtained values of the initial concentrations and parameters of the equations, stationary modes occur for the system of ordinary differential equations describing the case of the spatially uniform distribution of substances.Discussion and Conclusions. The mathematical model of the transformation of phosphorus, nitrogen and silicon forms in the problem of phytoplankton dynamics is studied. Stationary modes for the system of ordinary differential equations are obtained, for which the values of the system parameters and initial conditions are determined. The results obtained can be used in further simulation of the phytoplankton dynamics considering the transformation of phosphorus, nitrogen and silicon compounds with account for convection-diffusion, salinity, and temperature.

Текст научной работы на тему «Исследование стационарных решений задачи динамики фитопланктона с учетом трансформации соединений фосфора, азота и кремния»

МЕХАНИКА MECHANICS

УДК 519.6 https://doi.org/10.23947/1992-5980-2019-19-1-4-12

Исследование стационарных решений задачи динамики фитопланктона с учетом трансформации соединений фосфора, азота и кремния *

Ю. В. Белова1, А. М. Атаян2, А. Е. Чистяков3, А. В. Стражко4**

1,2,3,4 Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация

Study on stationary solutions to the problem of phytoplankton dynamics considering transformation of phosphorus, nitrogen and silicon compounds***

Yu.V. Belova1, A.M. Atayan2, A.E. Chistyakov3, A.V. Strazhko4**

1,2,3,4 Don State Technical University, Rostov-on-Don, Russian Federation

Ö О T3

M

"c

и (U

Ü £ Л

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

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

Introduction. The solution to the problem of transformation of phosphorus, nitrogen and silicon forms is studied. This problem arises under modeling phytoplankton dynamics in shallow-water bodies including the Azov Sea. The phytoplankton dynamics model is formulated as a boundary value problem for the system of diffusion-convection-response equations and takes into account the absorption and release of nutrients by phytoplankton, as well as the transition of nutrients from one compound to another. To calculate the initial conditions and parameters of the equations under which the steady-state regime occurs, the software is developed, which is based on the model describing changes in phytoplankton concentrations without considering current effects. This model is represented by a system of inhomogeneous differential equations. Based on the developed software, the initial conditions and parameters of the phytoplankton dynamics model in the Azov Sea are calculated experimentally.

Materials and Methods. A 3D model of phytoplankton dynamics is considered taking into account the transformation of phosphorus, nitrogen and silicon compounds based on the system of nutrient transport equations. The case of a spatially uniform distribution of substances is considered to specify the parameters of the model at which the stationary modes occur. Because of simplification, a system of ordinary differential equations solved through the Runge-Kutta method is obtained. Research Results. The software is developed to specify the initial conditions and parameters of the phytoplankton dynamics model considering the transformation of phosphorus, nitrogen and silicon compounds.

Several numerical experiments are performed under the assumption that the development of phytoplankton is limited by a single biogenic substance. As a result of the computational experiment, it can be seen that with the obtained values of the

fee

©

Работа выполнена при поддержке РНФ (проект № 17-11-01286). '* E-mail: [email protected], [email protected], [email protected], [email protected] '** The research is supported by the RSF (project no. 17-11-01286).

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

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

initial concentrations and parameters of the equations, stationary modes occur for the system of ordinary differential equations describing the case of the spatially uniform distribution of substances.

Discussion and Conclusions. The mathematical model of the transformation of phosphorus, nitrogen and silicon forms in the problem of phytoplankton dynamics is studied. Stationary modes for the system of ordinary differential equations are obtained, for which the values of the system parameters and initial conditions are determined. The results obtained can be used in further simulation of the phytoplankton dynamics considering the transformation of phosphorus, nitrogen and silicon compounds with account for convection-diffusion, salinity, and temperature.

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

Образец для цитирования. Белова, Ю. В. Исследование стационарных решений задачи динамики фитопланктона с учетом трансформации соединений фосфора, азота и кремния / Ю. В. Белова [и др.] // Вестник Донского гос. техн. ун-та. — 2019. — Т.19, №1. — С. 4-12. https://doi.org/10.23947/1992-5980-2Q19-19-1-4-12

Keywords: phytoplankton, phosphorus, nitrogen, silicon, biogen, chemical-biological source, convection-diffusion-response equation, Cauchy problem for system of ordinary differential equations, stationary mode.

For citation: Yu.V. Belova, et al. Study of stationary solutions to the problem of phytoplankton dynamics considering transformation of phosphorus, nitrogen and silicon compounds. Vestnik of DSTU, 2019, vol. 19, no. 1, pp. 4-12. https://doi.org/10.23947/1992-5980-2019-19-1-4-12

Введение. В связи с развитием крупных городов на побережье мелководных водоемов и речных систем, впадающих в данные водоемы, участились случаи эвтрофикации. Рост водорослей в водоемах вызван увеличением стока азотных и фосфорных соединений с прилегающих районов суши. Каждый водный объект является уникальным и требует подробного исследования. Для изучения водоемов прибегают к натурным исследованиям [1], а также математическому моделированию. Не преуменьшая роли натурных экспериментов можно сказать, что математическое моделирование является менее затратным, а также позволяет прогнозировать поведение экосистемы.

Для изучения Азовского моря разработана трехмерная модель гидродинамики [2, 3], включающая уравнения движения по трем пространственным направлениям. В работе [4] данная модель выполнена для случая динамического перестроения геометрии расчетной области за счет приливно -отливных явлений. Исследование точности данной модели приведено в работе [5]. В работах [6-8] приведена реконструкция экологической катастрофы, произошедшей в 2001 году, вызванной чрезмерной концентрацией водорослей в восточной части Азовского моря. В работе [9] предложены методы борьбы с заморными явлениями, возникающими в Азовском море. Работы [10-12] посвящены изучению динамики фито- и зоопланктона.

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

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

Материалы и методы. Модель основана на системе уравнений транспорта биогенных веществ [15, 16], вид которой для каждого модельного блока Fi имеет вид

дд

- + и — + v— + -W— = div ( k grad q )+ R

et dx ey cz ' q

«

S X a

<u

где q — концентрация i-ой компоненты, [мг/л]; ieM, M={F1, F2, F3, PO4, POP, DOP, NO3, NO2, NH4, Si}; {и,v,w} — компоненты вектора скорости водного потока, [м/с]; к - коэффициент турбулентного обмена, [м2/с]; Rq — функция-источник биогенных веществ, [мг/л-с].

В уравнении (1) индекс i указывает на вид субстанции (таблица 1).

Таблица 1

Биогенные вещества в модели динамики фитопланктона

Номер Обозначение Название

1 F1 зеленая водоросль Chlorella vulgaris

2 F2 синезеленая водоросль Aphanizomenon flos-aquae

3 F3 диатомовая водоросль Sceletonema costatum

4 PO4 фосфаты

5 POP взвешенный органический фосфор

6 DOP растворенный органический фосфор

7 NO3 нитраты

8 no2 нитриты

9 nh4 аммоний

10 Si растворенный неорганический кремний (кремниевые кислоты)

Химико-биологические реакции описаны следующими уравнениями

rf = cf(1~k—r)if -k—d4e -k—e , — = 13 ,

3

RPOP = SSPKFiDqFi KPDqPOP KPNqPOP , i=1

3

rdop = sspkfeqf + kpdqpop ~ kdnqdop , i=1

RPO = SSPCF (KF-~ 1) qFE — KPNqPOP + KDNqDOP , i=1

rnh4 = 23 SNCF (kfr~ -\ 4f + s3 SN (kfd + kfe ) ,

fn (чш-з > 4no2 > 4nh4 ) m

R = ± sC (K 1) f^ WW ^ ) ^ q + Kq Kq

Rno2 = sSNCF (KFR -^^Л---If + K42Vnh, -K234no2,

i=1 fn (4no, '4n02 '4nha ) qno2 + qno,

R C K 1 fx (qno3 'qm2' ) qno, K

Rno3 = s sncf \kfr 1I „ i \' qf +K23qno2 ,

i=1 fn (Чш, >Чш2 '4nha ) qno2 + qno,

RSi = SSCF3 (KF3R - 1) 4f, + SSKF3 d4f, •

Здесь KFR — удельная скорость дыхания фитопланктона; KFD — удельная скорость отмирания фитопланктона; KFE — удельная скорость экскреции фитопланктона; KPD — удельная скорость автолиза РОР; KPN — коэффициент фосфатофикации РОР; KDN — коэффициент фосфатофикации DOP; K42 — удельная скорость 3 окисления аммония до нитритов в процессе нитрификации; K23 — удельная скорость окисления нитритов до

и

й нитратов в процессе нитрификации, sP, sN, sS — нормировочные коэффициенты между содержанием N, P, Si в ТЗ

органическом веществе. 'й Скорость роста фитопланктона определяется выражениями:

cf:i = knfu m—n {fp (?PO4 ) • fn (4n0, 'яшг >4nh, )} , CF3 =knf3

min {fp ( 4po, ^ > fn (qno3 >4no-2 >4nh4 ) > fsi (4 s— )} ,

где KNF — максимальная удельная скорость роста фитопланктона.

(U

Ü С Л

Рис. 1. Модельная схема биогеохимической трансформации форм фосфора, азота и кремния

Функции, описывающие содержания биогенов

- для фосфора fr (qPOt ) = -

qPO4 ^ KPO4

где KPO — константа полунасыщения фосфатами;

- для кремния fs (qSi)=-

qSi

qSi + Ksi

где KSi — константа полунасыщения кремнием;

- для азота fN ( qm, qNOi ,qm ) = f( ( - qm2 , q^ ) + _fT ( q^ ) ,

(q^O, + qNO()exp(-K„JiqN^, )

qN

f(1) (q q q ) = y4NOi Чт^^УК f (2) (q s .N„4

Jn \4NO>4NO-,'4NHJ , ч > Jn \4NHJ

Кмо, + (Чно, + Чш-) 4 Кт4 + Чт,

где К — константа полунасыщения нитратами, К — константа полунасыщения аммонием, КрЛ —

коэффициент ингибирования аммония.

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

Ч^ (х, у, 2,0) = (х, у, 2), (х, у, г) е О ., / = 0, ¡еМ. (2)

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

- на боковой границе:

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

Ч = 0, на а, если ип < 0, ¡еМ; (3)

dq.

—L = 0, на ст , если ип > 0, i eM;

дп "

- на Eo — поверхность водоема:

dqL dz

= 0, ieM;

(4)

(5)

ей

и

S X а X <и

- на дне £„ :

7

k^ = Elq ,/e{Fb F2, F3}, dz

kdq- = e2q , /e{PO4, POP, DOP, NO3, NO2, NH4, Si}, (6)

dz '' '

где eu, e2i — скорости осаждения водорослей и питательных веществ на дно.

Стационарный режим. Рассмотрим случай пространственно-равномерного распределения субстанций (фитопланктона, форм фосфора, азота и кремния), тогда каждое из уравнений (1) упрощается и в результате получаем следующую систему обыкновенных дифференциальных уравнений (ОДУ):

dqF —

—7 = CF (1 - KFR )qF - KFDqF - KFEqF , i = 1,3, (7)

dt i i i

s K n — K П — K П

' р^ррчр. лpdhpop ^pn^pop •>

lpop dt ¡=i

dqDop _ s к n A- K n — K n

■ ¿J p^fmhf, ^pdhpop ^dnhd

dop dt ¡=i

qPO" ■ - A SPCF (KFR - 1) nF + KPNqP0P + KDNqD

_0

dt i-i

dq 3 f '2) (NH ) 3

ÜT - Ü (KFR - 1) ^ ( N03, N024 NH4) ^ + Ü—1 (KFD + KFE ) * - K42 qNH4

dnN02 - 3 (K -1) ff (N03, N02, NH4) nN02 n + Kn - Kn

dt 2SnC(Kfr -1) /nno2,NH4) ■ nN02 + 2 ^ + ^ -^

dnn0, ± r {„ NO3, NO2, NH4 ) nN03 „ -3 — A sNCF (^^^ -1)-7-f--3-nP + K23 ^^

dt 1—1 N ^ FR 'In (NO3, NO2, NH4) qNOl + nN03 ■ 2

= (KF3R - 1) qP3 + ssKFDqF, •

Решим систему обыкновенных дифференциальных уравнений методом Рунге -Кутты [15-17]. Проведем несколько численных экспериментов, предполагая, что развитие фитопланктона зависит от единственного лимитирующего вещества.

Результаты исследования. Для системы ОДУ (7) рассчитаем начальные условия и параметры уравнений, при которых наступают стационарные режимы. Возьмем начальные значения концентраций: qF (0) = 2,5 мг/л, qF (0) = 2,6 мг/л, qF (0) = 0,91 мг/л, qPOP (0) = 0,07 мг/л, qDOP (0) = 0,07 мг/л, qO (0) = 0,005

мг/л, qm (0) = 0,11 мг/л, qO (0) = 0,0178 мг/л, qO (0) = 0,304 мг/л, qsi (0) = 0,4 мг/л; коэффициенты: KF = 2,8 сут'1, KFR = 0,15 сут'1, KFD = 0,05 сут'1, KFE = 0,15 сут'1, KPD = 0,015 сут'1, KPN = 0,02 сут'1, KDN = 0,1 сут'1, K42 = 0,9 сут'1, K23 = 2,5 сут'1, Kpsi = 1,46 сут'1, sP = 0,01, sN = 0,016, sSi = 0,023, KO g = 0,024, Kn0i = 3,0, Knh4 = 2,0, Ks = 3,0.

.¡3 Полученные стационарные режимы системы ОДУ (7) в предположении, что развитие фитопланк-

СЛ

тона лимитируется единственным биогенным веществом —фосфором, азотом или кремнием — изобра-

М жены на рис. 2-4 соответственно. Рис. 2 описывает влияние фосфора на развитие различных видов фи-•д

и топланктона, рис. 3 описывает влияние азота на развитие различных видов фитопланктона, рис. 4 опи-<и

сывает влияние азота на развитие диатомовых водорослей.

а) зеленая водоросль (ChV), б) синезеленая водоросль (AF-A), в) диатомовая водоросль (SC), г) взвешенный органический фосфор (POP), д) растворенный органический фосфор (DOP), е) фосфаты (PO4)

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

концентраций и параметров уравнений наступают стационарные режимы для системы ОДУ (7), описывающей

случай пространственно-равномерного распределения субстанций. Полученные значения в дальнейшем будут

использованы для моделирования пространственно-неоднородного распределения субстанций, солености и

температуры с учетом движения водной среды [18]. ijS

К

Ch V, мг/л 10

AF-Aj мг/л

0 1 SC, мг/л A

NO2, мг/л

,

0,023 , .

,

xicr

,

9

а]

/

/

/

4 110 с 0 1 NH4, мг/л

в)

/

/

а)

1/

-

J t, 10®c 0 1

NOjj мг/л

,

,

,

0 1

,

и б)

4 t, Иге

г)

,

*\г

4 110 с

4 t, 10 с ' о 1

4 t, Иге

Ö О ТЗ

M

■¡з

и (U

Ü С Л

Рис. 3. Стационарный режим системы ОДУ в предположении, что развитие фитопланктона лимитируется азотом: а) зеленая водоросль (СНУ), б) синезеленая водоросль (ЛР-А), в) диатомовая водоросль (БС), г) аммоний (МИ4), д) нитриты (М02), е)

нитраты (М03)

.SC, мг/л

Юг

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

а)

/

/

/ 7

Si, мг/л

,

,

,

,

\ б)

4 110 с

Рис. 4. Стационарный режим системы ОДУ в предположении, что развитие фитопланктона (диатомовых водорослей) лимитируется кремнием: а) диатомовая водоросль (БС), б) кремний (Б)

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

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

Библиографический список

1. Якушев, Е. В. Комплексные океанологические исследования Азовского моря в 28-м рейсе научно-исследовательского судна «Акванавт» / Е. В. Якушев, А. И. Сухинов // Океанология. — 2003. — Т. 43. — № 1.

— С.44-53.

2. Сухинов, А. И. Численное моделирование экологического состояния Азовского моря с применением схем повышенного порядка точности на многопроцессорной вычислительной системе / А. И. Сухинов [и др.] // Компьютерные исследования и моделирование. — 2016. — Т. 8. — № 1. — С. 151-168.

3. Sukhinov, A.I., Sukhinov A.A. 3D model of diffusion-advection-aggregation suspensions in a water basins and its parallel realization. Parallel Computational Fluid Dynamics 2004: Multidisciplinary Applications — 2005. — pp. 223-230. DOI: 10.1016/B978-044452024-1/50029-4.

4. Sukhinov, A.I., Chistyakov, A.E., Shishenya, A.V., Timofeeva, E.F. Mathematical model for calculating coastal wave processes. Mathematical Models and Computer Simulations. — 2013. — Т. 5. — № 2. — pp. 122-129. DOI: 10.1134/S2070048213020087.

5. Сухинов, А. И. Предсказательное моделирование прибрежных гидрофизических процессов на многопроцессорной системе с использованием явных схем / А. И. Сухинов [и др.] // Математическое моделирование. — 2018. — Т. 30. — № 3. — С. 83-100.

6. Сухинов, А. И. Математическое моделирование условий формирования заморов в мелководных водоемах на многопроцессорной вычислительной системе / А. И. Сухинов [и др.] // Вычислительные методы и программирование. — 2013. — Т. 14. — № 1. — С. 103-112.

7. Sukhinov, A.I., Sukhinov A.A. Reconstruction of 2001 ecological disaster in the Azov sea on the basis of precise hydrophysics models. Parallel Computational Fluid Dynamics 2004: Multidisciplinary Applications — 2005. — pp. 231-238. DOI: 10.1016/B978-044452024-1/50030-0.

8. Debolskaya, E.I., Yakushev, E.V., Sukhinov, A.I. Formation of fish kills and anaerobic conditions in the sea of Azov. Water Resources. — 2005. — Т. 32. — № 2. — pp. 151-162. DOI: 10.1007/s11268-005-0020-5.

9. Никитина, А. В. Оптимальное управление устойчивым развитием при биологической реабилитации Азовского моря / А. В. Никитина [и др.] // Математическое моделирование. — 2016. — Т. 28. — № 7. — С. 96106.

10. Сухинов, А. И. Математическое моделирование процессов эвтрофикации в мелководных водоемах на многопроцессорной вычислительной системе / А. И. Сухинов, А.В. Никитина, А.Е. Чистяков // Вестник Южно-Уральского государственного университета. Серия: Вычислительная математика и информатика. — 2016. — Т. 5. — № 3. — С. 36 -53.

11. Никитина, А. В. Дифференциально-игровая модель предотвращения заморов в мелководных водоемах / А. В. Никитина, М. В. Пучкин, И. С. Семенов // Управление большими системами. — 2015. — Вып. 55.

— C. 343-361.

12. Сухинов, А. И. Математическая модель трансформации форм фосфора, азота и кремния в движущейся турбулентной водной среде в задачах динамики планктонных популяций / А. И. Сухинов, Ю. В. Белова // Инженерный вестник Дона. — 2015. — Т. 37. — № 3. — C. 50.

13. Дегтярева, Е. Е. Программная реализация трехмерной математической модели транспорта взвеси в мелководных акваториях [Электронный ресурс] / Е. Е. Дегтярева, Е. А. Проценко, А. Е. Чистяков // Инженерный вестник Дона. — 2012. — Т. 23. — № 4 — 2. -30 с. — Режим доступа : ivdon.ru/ru/magazine/archive/n4p2y2012/1283 (дата обращения : 12.12.2018).

14. Самарский, А. А. Теория разностных схем / А. А. Самарский — Москва : Наука, 1989. — 616 с.

15. Сухинов, А. И. Достаточные условия сходимости положительных решений линеаризованной двумерной задачи транспорта наносов / А. И. Сухинов, В. В. Сидорякина, А. А. Сухинов // Вестник Донского гос. техн. ун-та. — 2017. — Т. 17. — № 1 (88). — С. 5-17.

16. Самарский, А. А. Методы решения сеточных уравнений / А. А. Самарский, Е. С. Николаев. — Москва : Наука, 1978. — 532 с. ^

17. Марчук, Г. И. Математическое моделирование в проблеме окружающей среды / Г. И. Марчук. — Москва : Наука, 1982. — 319 с.

18. Белова, Ю. В. О четырехслойной итерационной схеме / Ю. В. Белова, А. Е. Чистяков,

У. А. Проценко //Вестник Донского гос. техн. ун-та. 2016. — Т. 16. — № 4 (87). — С. 146-149. 11

и s

(U

Поступила в редакцию 20.11.2018 Сдана в редакцию 21.11.2018 Запланирована в номер 11.01.2019

Received 20.11 .2018 Submitted 21.12.2018 Scheduled in the issue 11.01.2019

Об авторах:

Белова Юлия Валериевна,

младший научный сотрудник НИИ «Математическое моделирование и прогнозирование сложных систем» Донского государственного технического университета, (РФ, 344000, г. Ростов-на-Дону, пл. Гагарина, 1), ОЯСГО: https://orcid.org/0000-0002-2639-7451 [email protected]

Атаян Ася Михайловна,

ассистент кафедры «Программное обеспечение вычислительной техники и автоматизированных систем» Донского государственного технического университета, (РФ, 344000, г. Ростов-на-Дону, пл. Гагарина, 1),

ОКСГО: https://orcid.org/0000-0003-4629-1002 [email protected]

Authors:

Belova, Yuliya V.,

Junior research scholar, Research Institute for Mathematical Modeling and Forecasting of Complex Systems, Don State Technical University (1, Gagarin sq., Rostov-on-Don, 344000, RF),

ORCID: https://orcid.org/0000-0002-2639-7451 [email protected]

Atayan, Asya M.,

teaching assistant of the Computer and Automated Systems Software Department, Don State Technical University (1, Gagarin sq., Rostov-on-Don, 344000, RF), ORCID: https://orcid.org/0000-0003-4629-1002 [email protected]

Чистяков Александр Евгеньевич,

профессор кафедры «Программное обеспечение вычислительной техники и автоматизированных систем» Донского государственного технического университета, (РФ, 344000, г. Ростов-на-Дону, пл. Гагарина, 1), доктор физико-математических наук, профессор,

ORCID: https://orcid.org/0000-0002-8323-6005 cheese [email protected]

Стражко Александр Валентинович,

студент кафедры «Программное обеспечение вычислительной техники и автоматизированных систем» Донского государственного технического университета, (РФ, 344000, г. Ростов-на-Дону, пл. Гагарина, 1), ORCID: https://orcid.org/0000-0002-2449-8531 [email protected]

Chistyakov, Aleksandr E.,

professor of the Computer and Automated Systems Software Department, Don State Technical University (1, Gagarin sq., Rostov-on-Don, 344000, RF), Dr.Sci. (Phys.-Math.), professor, ORCID: https://orcid.org/0000-0002-8323-6005 cheese [email protected]

Strazhko, Aleksandr V.,

student of the Computer and Automated Systems Software Department, Don State Technical University (1, Gagarin sq., Rostov-on-Don, 344000, RF), ORCID: https://orcid.org/0000-0002-2449-8531 [email protected]

с о

T3

"c

и

(U >

Л £ Л

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