�� sr !com.femlab.server.ModelFileHeader�D���% L tagst Lcom/femlab/util/FlStringList;L typesq ~ L vrsnt Lcom/femlab/util/FlVersion;xpw sr com.femlab.util.FlVersion��%�/B = I buildI majorL datet Ljava/lang/String;L extq ~ L nameq ~ L rcsq ~ L reactionExtq ~ L reactionNameq ~ xpw t q ~ t #COMSOL Reaction Engineering Lab 3.5t at COMSOL 3.5q ~ w [t $Name: $t $Date: 2008/12/03 17:02:19 $xur [Ljava.lang.String;��V�� {G xp t modelinfot xfemt guit g11t g12t g13t g7t g3t fem1uq ~ q ~ t femstructt guistructt drawq ~ q ~ q ~ q ~ t geomxsr com.femlab.api.client.ModelInfo� ^���% L descrq ~ L docURLq ~ [ imaget [Bxpw pt pxuq ~ t�`
COMSOL 3.5
a
0
603
$Name: $
$Date: 2008/12/03 17:02:19 $
1
geom
on
on
off
on
on
off
on
0.1
t
Darcy
cartesian
{"p"}
{"x","y","z"}
esdl
ES
{"shlag(2,'p')"}
4
2
2
off
_esdl
Lag2
H
static
off
ref
off
{"lm1"}
ideal
{1,1,1,1,1,1,1,1}
0
0
0
0
0
0
0
{1,2,1,3,4,4,4,5,5}
{"0"}
1
userdef
1
0.25
1
1
1
conductivity
0
1000
0.001
0
1
1
1
1
{1,2}
1e-5
9.81
0
SoluteTransSat
cartesian
{"c"}
{"x","y","z"}
esst
ES
{"shlag(2,'c')"}
4
2
2
off
_esst
Lag2
noncons
solid
time
off
{}
{}
{}
{}
ref
off
{"lm2"}
ideal
{1,1,1,1,1,1,1,1}
{"0"}
{"1"}
{"0"}
{1,2,1,3,2,4,2,1,1}
{"0"}
1
{"userdef"}
dthdt0
0
1
1
0.36
1
u_esdl
v_esdl
0
{"1"}
{"0"}
{"0"}
{"0"}
{"1"}
{"0"}
{"0"}
{"0"}
{"0"}
{"0"}
{"0"}
computed
computed
{"0"}
{"0"}
{"0"}
{"0"}
0
0
0
0
{1,2}
1e-5
FlPDEC
cartesian
{"S","S_t"}
{"x","y","z"}
c
{"shlag(2,'S')"}
4
2
2
off
_c
Lag2
off
ref
off
{"lm3","lm4"}
{"0"}
{"0"}
{"0"}
{"0"}
ideal
{1,1,1,1,1,1,1,1}
{"0"}
{"0"}
{"0"}
{"0"}
ideal
{"0"}
{"1"}
{"0"}
{"0"}
neu
{1,1,1,1,1,1,1,1,1}
1
{"0"}
{"0"}
{"0"}
{"0"}
ideal
{"0"}
{"0"}
{"0"}
{"rho_s/n"}
{1,2}
{"x","y"}
{"ref"}
{"shlag(2,'p')","shlag(2,'c')","shlag(2,'S')"}
4
2
2
on
1
general
SI
{1,1}
{"p","c","S"}
CSs_esdl/(rhof_esdl*g_esdl)
0
Ks_esdl*CKs_esdl
K_esdl
0
0
K_esdl
kaps_esdl*CKs_esdl
kap_esdl
0
0
kap_esdl
p/(rhof_esdl*g_esdl)+D_esdl
sqrt(px^2+py^2)
d(H_esdl,x)
d(H_esdl,y)
sqrt(Hx_esdl^2+Hy_esdl^2)
(-Kxx_esdl*(px+rhof_esdl*g_esdl*d(D_esdl,x))-Kxy_esdl*(py+rhof_esdl*g_esdl*d(D_esdl,y)))/(rhof_esdl*g_esdl)
(-Kyx_esdl*(px+rhof_esdl*g_esdl*d(D_esdl,x))-Kyy_esdl*(py+rhof_esdl*g_esdl*d(D_esdl,y)))/(rhof_esdl*g_esdl)
sqrt(u_esdl^2+v_esdl^2)
Qs_esst*cQ_c_esst+So_c_esst
1
alpha1_c_esst*u_esst^2/(U_esst+eps)+alpha2_c_esst*v_esst^2/(U_esst+eps)+theta_esst*tauL_c_esst*DmL_c_esst
(alpha1_c_esst-alpha2_c_esst)*u_esst*v_esst/(U_esst+eps)
(alpha1_c_esst-alpha2_c_esst)*u_esst*v_esst/(U_esst+eps)
alpha1_c_esst*v_esst^2/(U_esst+eps)+alpha2_c_esst*u_esst^2/(U_esst+eps)+theta_esst*tauL_c_esst*DmL_c_esst
U_esst*h*(u_esst^2+v_esst^2+eps)/(thDLxx_c_esst*u_esst^2+thDLxy_c_esst*u_esst*v_esst+thDLyx_c_esst*v_esst*u_esst+thDLyy_c_esst*v_esst^2)
cx
cy
sqrt(cx^2+cy^2)
sqrt(u_esst^2+v_esst^2)
u_esst/theta_s_esst
v_esst/theta_s_esst
sqrt(ua_esst^2+va_esst^2)
U_esst*c
-thDLxx_c_esst*cx-thDLxy_c_esst*cy
-thDLxx_c_esst*cx-thDLxy_c_esst*cy+u_esst*c
u_esst*c
-thDLyx_c_esst*cx-thDLyy_c_esst*cy
-thDLyx_c_esst*cx-thDLyy_c_esst*cy+v_esst*c
v_esst*c
sqrt(dflux_c_x_esst^2+dflux_c_y_esst^2)
sqrt((dflux_c_x_esst+u_esst*c)^2+(dflux_c_y_esst+v_esst*c)^2)
0
d(kP_c_esst*c,c)
kP_c_esst*c
1+rhob_c_esst*kP_c_esst/theta_esst
ua_esst/RF_c_esst
va_esst/RF_c_esst
sqrt(Sx^2+Sy^2)
sqrt(ga3x^2+ga3y^2)
0.25
1
1
1
1000
0.001
1
1
1
1
1
1
0.36
theta_s_esst
u_esdl
v_esdl
0
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
{0,0}
{1,2}
{"p","c","S"}
u_esdl*nx_esdl+v_esdl*ny_esdl
u_esdl*nx_esdl+v_esdl*ny_esdl
nx_esst*dflux_c_x_esst+ny_esst*dflux_c_y_esst
nx_esst*adflux_c_x_esst+ny_esst*adflux_c_y_esst
nx_esst*tflux_c_x_esst+ny_esst*tflux_c_y_esst
0
0
0
0
0
0
0
nx
ny
nx
ny
t!�
0
1
0
nx
ny
{0,0,0,0,0,0,0,0,0}
{1,2,1,3,4,5,4,6,6}
{"p","c","S"}
1000
{0,0,0,0,0,0,0,0}
{1,1,1,1,1,1,1,1}
1e-005
9.810000000000001
0
1e-005
{"B1","B2","B3","B4"}
{"R1"}
4.5e-8[m]
5e-6[m]
2600[kg/m^3]
205[m]
3.75[m]
1e-7[m/s]
11.25[m]
1e-5[m/s]
0.36
0.002
1000[kg/m^3]
0.001002[kg/m/s]
15[m]
13.975[m]
1.6e-5[kg/kg]
3.52e-6[kg/kg]
9.81[m/s^2]
10e-3[kg/m^3]
1600[kg/m^3]
1.38e-23[J/K]
293[K]
0.0002
sqrt(1.84e-5*K1*(1-n)^2/n^3)
sqrt(1.84e-5*K2*(1-n)^2/n^3)
dp/lp
4.03*K*T*log((1+(1-(R^2))^0.5)/R)^0.67
3*3.1416*visco*dc*u_esdl*(1.5*dp^2*lp)^0.33*(1/R)^0.67
3*3.1416*visco*ds*u_esdl*(1.5*dp^2*lp)^0.33*(1/R)^0.67
((1-(R^2))^0.5)^(-0.67)
(0.146*(rho_p-rho_f)*g*((dp^2)*lp)^0.67)/(visco*u_esdl*(1-(R^2))^0.5)
log((1+(1-(R^2))^0.5)/R)
A1*B1*C1
R^0.67*D1*E1
0.5*(dp/dc)^2*(3-(dp/(dp+dc)))
nd1+ng1+nI1
A1*B2*C1
0.5*(dp/ds)^2*(3-(dp/(dp+ds)))
nd2+ng1+nI2
(-1.5*((1-n)/dc)*alpha1*n01*u_esdl*Patt1)
(-1.5*((1-n)/ds)*alpha2*n02*u_esdl*Patt2)
(1-(S/Smax1))*(S<=Smax1)
(1-(S/Smax2))*(S<=Smax2)
{0}
{0}
esdl
weak
SI
uq ~ t��
off
time
stationary
auto
1.0E-6
25
tol
100
1
off
on
off
1.0
1.0E-4
10.0
off
on
0.0010
25
lumped
automatic
auto
on
1
4
tol
10
1
off
const
1.0
minimal
1.0
1.0E-2
10.0
0.0000000010
0.00000001
86400*range(0,10,3650)
0
1
0.1
11
step
none
tlist
free
bdf_ida
0.01
off
15
off
on
on
on
on
on
on
off
on
on
0.0010
5
1
0.1
0.75
linear
off
off
maybe
bweuler
0
off
6
0
300
0.0
0
lambda
0
step
none
off
1
off
0.0
0.0
0.0
on
auto
auto
auto
eliminate
auto
1000
on
auto
on
off
off
on
off
currgeom
l2
10000000
longest
weak
on
0
1
2
1
fltpft
2
1.7
0.5
0.5
on
off
off
sensitivity
adjoint
on
cholesky
analytic
500
10
1.0e-6
1.0e-6
1.0e-6
1.0e-6
0
off
c,S
c,S,p
on
init_expr_currsol_radio
Automatic
0
u_init_radio
Automatic
0
off
off
off
off
umfpack
linsolver
0.0
0.1
0.7
nd
on
off
nd
on
1.0E-8
on
off
tol
2
1.0E-6
400.0
10000
left
right
1
50
1
1.0
3
6
v
5000
512.0
off
2.0
on
droptol
gmres
0.02
100
0.8
on
on
shape
regular
on
2
2
off
Geom1
on
on
off
off
off
Geom1
0
off
off
off
2076121
COMSOL 3.5.0.603
1340916288241
1326317526059
off
html
off
a4
off
full
off
sr com.femlab.geom.Geom2�V���c/O xr com.femlab.geom.Geom�6��{ L geomAssoct Lcom/femlab/geom/GeomAssoc;L versionq ~ xpw q ~ w ur [B���T� xp G e o m 2 �����|�=-C ��6? � .@���� � 9@ .@���� � �? � B e z i e r C u r v e � .@ �? 9@ .@ �? A s s o c A t t r i b V e c t o r I n t $ xw xsq ~ &w q ~ w uq ~ * G e o m 2 �����|�=-C ��6? 9@ .@���� � @Z@ .@���� � �? � B e z i e r C u r v e 9@ .@ �? @Z@ .@ �? A s s o c A t t r i b V e c t o r I n t $ xw xsq ~ &w q ~ w uq ~ * G e o m 2 �����|�=-C ��6? @Z@ .@���� � �i@ .@���� � �? � B e z i e r C u r v e @Z@ .@ �? �i@ .@ �? A s s o c A t t r i b V e c t o r I n t $ xw xsq ~ &w q ~ w uq ~ * G e o m 2 �����|�=-C ��6? � �&@���� � �i@ �&@���� � �? � B e z i e r C u r v e � �&@ �? �i@ �&@ �? A s s o c A t t r i b V e c t o r I n t $ xw xsr com.femlab.geom.Rect�ȉ,E��6 D lxD lyL lxExprt !Lcom/femlab/geom/Geom$Expression;L lyExprq ~ 3xr com.femlab.geom.Prim2�Z��~�� D rotL baseq ~ L constrq ~ [ post [D[ posExprt "[Lcom/femlab/geom/Geom$Expression;L rotExprq ~ 3xq ~ &w q ~ w �uq ~ * � G e o m 2 �����|�=-C ��6? ���� � �i@ ���� � �i@ .@���� � .@���� � �? � �? � �? � �? � B e z i e r C u r v e �? �i@ �? B e z i e r C u r v e �i@ �? �i@ .@ �? B e z i e r C u r v e �i@ .@ �? .@ �? B e z i e r C u r v e .@ �? �? A s s o c A t t r i b V e c t o r I n t - xw xw t cornerur [D>���cZ xp ur "[Lcom.femlab.geom.Geom$Expression;�\2�Y�U xp sr com.femlab.geom.Geom$Expression OpI_h� D numScaleL exprStrq ~ xpw pw?� xsq ~ >w pw?� xw sq ~ >w pw?��F�R�9xxw @i� sq ~ >w pw?� xw@. sq ~ >w pw?� xxsq ~ &w q ~ w �uq ~ * � G e o m 2 �����|�=-C ��6? ���� � �&@���� � ������-@���� � 9@������-@���� � @Z@������-@���� � �i@ ���� � �i@ �&@���� � �i@������-@���� � �? ���� � �? � �? ���� � �? � �? � �? � �? � �? � �? � B e z i e r C u r v e �? �����&@ �? B e z i e r C u r v e �? �i@ �? B e z i e r C u r v e �����&@ �? ������-@ �? B e z i e r C u r v e �&@ �? �i@ �&@ �? B e z i e r C u r v e ������-@ �? 9@������-@ �? B e z i e r C u r v e 9@������-@ �? @Z@������-@ �? B e z i e r C u r v e @Z@������-@ �? �i@������-@ �? B e z i e r C u r v e �i@ �? �i@�����&@ �? B e z i e r C u r v e �i@�����&@ �? �i@������-@ �? A s s o c A t t r i b V e c t o r I n t 9 xw x