Компьютерное моделирование биотехнологических процессов и систем. Дворецкий Д.С - 43 стр.

UptoLike

ВУЗ: 

ТГТУ | Тамбов

Составители: 

Рубрика: 

0
*
2
2
*
=+
∂
∂
−
SS
R
r
S
D
; (4.43)
0
*
2
2
*
=−
∂
∂
−
CC
R
r
C
D
. (4.44)
Осложнение в решении этой системы уравнений состоит в том, что макро- и микробаланс опреде-
лены в различных системах координат. Эта проблема является общей при решении многих задач в био-
технологии, и часто решаются, используя аналитические приближения решения микробаланса. Однако
средства MatLab позволяют произвести формулировку задачи в единой системе координат. Порядок
решения данной задачи иллюстрирует программный код (листинг 1).
Листинг 1
% Сначала устанавливаем параметры решаемой задачи и геометрию области решения.
clear xfem fem1 fem2 appl1 appl2 appl3 elem
fem1.geom=rect2;
fem1.sdim={'x' 'r'};
fem1.mesh=meshinit(fem1,'hmax',{0.2,[],[3 0.01]});
fem2.geom=solid1([0 1]);
fem2.sdim='x';
fem2.mesh=meshinit(fem2,'hmax',0.02);
appl1.mode=flchedi2d;
appl1.dim='c1';
appl1.elemdefault='Lag2';
appl2.mode=flpdewc2d;
appl2.dim='lambda';
appl2.elemdefault='Lag2';
appl3.mode=flchecd1d;
appl3.dim='C';
appl3.elemdefault='Lag2';
% Задаем константы и переменные:
xfem.const={'gamma', 1.5,...
'D1', 1e-07,...
'D', 1e-06,...
'C0', 3,...
'k', 100,...
'epsil', 0.6,...
'Ap', 1000,...
'rp', 0.001,...
'u', 0.4};
elem{1}.elem=elcplextr;
elem{1}.src.g=2;
elem{1}.src.equ.meshp={{'x'}};
elem{1}.src.equ.var={'C_x',{'C'}};
elem{1}.src.equ.ind={1};
elem{1}.dst.g=1;
elem{1}.dst.bnd.ep={{'x'}};
elem{1}.dst.bnd.ind={3};
elem{2}.elem=elcplextr;
elem{2}.src.g=1;
elem{2}.src.bnd.meshp={{'x'}};
elem{2}.src.bnd.var={'Ndotn_x',{'D1*c1r/rp'}};
elem{2}.src.bnd.ind={3};
elem{2}.dst.g=2;

Страницы