�� 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 fem0t fem1t g1t fem1t fem1.0q ~ q ~ t fem6t fem6.0q ~ q ~ t fem5t fem5.0q ~ q ~ uq ~ q ~ t femstructt guistructq ~ q ~ t drawt geomt mesht solutiont xmeshq ~ q ~ q ~ !q ~ "q ~ q ~ q ~ !q ~ "xsr 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
on
on
0.1
t
NavierStokes
cartesian
{"u","v","p","logk","logd","logw","phi","psi","nxw","nyw"}
{"x","y","z"}
chns
CHEM
{"shlag(2,'u')","shlag(2,'v')","shlag(1,'p')","shlag(2,'logk')","shlag(2,'logd')","shlag(2,'phi')","shlag(2,'psi')"}
2
off
_chns
Lagp2p1
time
Off
Off
k-epsilon
Off
Off
Off
pfnoncons
ref
off
{"lm1","lm2","lm3","lm4","lm5","lm6","lm7","lm8"}
ideal
0
0
{1,1,1,1}
uv
wall
ntotstress
totstress
cont
U0in
0
0
0
0
1
0
0
0
0
0
0.005
0.005
10
LTIT
open
1
0.01
0.05
dw
h/2
100
0
0
0
1
0
1
velocity
0
0
pi/2
h
{1,2,2,3}
1
1
1
0
0
1
0
0
0
0
power
0
0
0
1
1
0
0
0
0.5
0
pgc
0.25
0
sc
0.35
0
1
0.1
0
navierstokes
usesELM
0
Lagp2p1
0
Lagp2p1
0
0
0
0.25
1
gls
0.5
1
1
Lagp2p1
isoT
1000
1.2
1e-5
2e-5
7e-2
1
50
hmax_chns*0.5*1.2^2
0
1
{1}
10
1.44
1.92
1.0
1.3
0.09
5.5
0.42
{"x","y"}
{"ref"}
{"shlag(2,'u')","shlag(2,'v')","shlag(1,'p')","shlag(2,'logk')","shlag(2,'logd')","shlag(2,'phi')","shlag(2,'psi')"}
2
on
1
SI
{1}
{"u","v","p","logk","logd","phi","psi"}
sqrt(u^2+v^2)
vx-uy
ux+vy
etaT_chns/sigmad_chns+eta_chns
Cd1_chns*etaT_chns*exp(-logk)*pT_chns-Cd2_chns*rho_chns*exp(logd-logk)
exp(logd)
rho_chns*Cmu_chns*exp(2*logk-logd)
2*ux^2+2*vy^2+(uy+vx)^2
exp(logk)
etaT_chns/sigmak_chns+eta_chns
etaT_chns*exp(-logk)*pT_chns-rho_chns*exp(logd-logk)
-Dlogk_chns*(logkxx+logkyy)+rho_chns*(logkt+u*logkx+v*logky)-Dlogk_chns*(logkx^2+logky^2)-Slogk_chns
-Dlogd_chns*(logdxx+logdyy)+rho_chns*(logdt+u*logdx+v*logdy)-Dlogd_chns*(logdx^2+logdy^2)-Slogd_chns
3*sigma_chns*epsilon_chns/sqrt(8)
psi*lam_chns/epsilon_chns^2
min(max(0.5*(1+phi),0),1)
chi_chns*epsilon_chns^2
0.5*lam_chns*(phix^2+phiy^2)+0.25*lam_chns*(-1+phi^2)^2/epsilon_chns^2+dfdphi_chns*phi
(G_chns-dfdphi_chns)*phix
(G_chns-dfdphi_chns)*phiy
dfdphi_chns*phix
dfdphi_chns*phiy
hmax1_chns
1-Vf_chns
Vf_chns
rho1_chns+(rho2_chns-rho1_chns)*Vf_chns
eta1_chns+(eta2_chns-eta1_chns)*Vf_chns
rho_chns*U_chns*h/eta_chns
rho_chns*(u*ux+v*uy)+px-F_x_chns-(eta_chns+etaT_chns)*(2*uxx+uyy+vxy)-G_chns*phix
rho_chns*(u*vx+v*vy)+py-F_y_chns-(eta_chns+etaT_chns)*(vxx+uyx+2*vyy)-rho_chns*gy_chns-G_chns*phiy
rho_chns*u
rho_chns*v
eta_chns+etaT_chns
rho_chns
nojac(1/max(2*rho_chns*sqrt(emetric(u,v)),48*(eta_chns+etaT_chns)/h^2))
0.5*nojac(if(u^2+v^2
rho_chns*divU_chns
nojac(0.5*h/max(rho_chns*U_chns,6*Dlogk_chns/h))
nojac(0.5*h/max(rho_chns*U_chns,6*Dlogd_chns/h))
0
0
1
0
0
0
0
power
0
0
0
1
0
0
0.5
0.25
sc
0.35
1
0.1
0
Lagp2p1
Lagp2p1
0
0
Lagp2p1
1000
1.2
1e-005
2e-005
0.07000000000000001
1
50
0.72*hmax_chns
0
1
0
-9.810000000000001
U_chns
{0}
{1}
{"u","v","p","logk","logd","phi","psi"}
(eta_chns+etaT_chns)*(2*nx_chns*ux+ny_chns*(uy+vx))
-nx_chns*p+2*nx_chns*(eta_chns+etaT_chns)*ux+ny_chns*(eta_chns+etaT_chns)*(uy+vx)
(eta_chns+etaT_chns)*(nx_chns*(vx+uy)+2*ny_chns*vy)
-ny_chns*p+nx_chns*(eta_chns+etaT_chns)*(vx+uy)+2*ny_chns*(eta_chns+etaT_chns)*vy
-1+2*Vf0_chns
phix*nx_chns+phiy*ny_chns
phix-nx_chns*ndot_gradphi_chns
phiy-ny_chns*ndot_gradphi_chns
sqrt(Wx_chns^2+Wy_chns^2+eps)
epsilon_chns^2*tan(0.5*pi-theta_chns)*W_chns
0
0
0
0
1
0
0
0
0
0
exp(logk)
0.005
10
1
0.01
0.05
0.5*h
rho_chns*dw_chns*Cmu_chns^0.25*exp(0.5*logk)/eta_chns
0
0
0
1
0
1
0
0.5*pi
h
nx
ny
0
0
{0,0,0,0}
{1,2,2,3}
{"u","v","p","logk","logd","phi","psi"}
0
{0,0,0,0}
{1,1,1,1}
10
1.44
1.92
1
1.3
0.09
5.5
0.42
elcplscalar
{"1"}
{"hmax1_chns"}
{"1"}
{"1"}
{"R1"}
{0}
{0}
chns
weak
SI
COMSOL 3.5
a
0
603
$Name: $
$Date: 2008/12/03 17:02:19 $
1
geom
on
on
off
on
on
on
on
0.1
t
NavierStokes
cartesian
{"u","v","p","logk","logd","logw","phi","psi","nxw","nyw"}
{"x","y","z"}
chns
CHEM
{"shlag(2,'u')","shlag(2,'v')","shlag(1,'p')","shlag(2,'logk')","shlag(2,'logd')","shlag(2,'phi')","shlag(2,'psi')"}
2
off
_chns
Lagp2p1
time
Off
Off
k-epsilon
Off
Off
Off
pfnoncons
ref
off
{"lm1","lm2","lm3","lm4","lm5","lm6","lm7","lm8"}
ideal
0
0
{1,1,1,1}
uv
wall
ntotstress
totstress
cont
U0in
0
0
0
0
1
0
0
0
0
0
0.005
0.005
10
LTIT
open
1
0.01
0.05
dw
h/2
100
0
0
0
1
0
1
velocity
0
0
pi/2
h
{1,2,2,3}
1
1
1
0
0
1
0
0
0
0
power
0
0
0
1
1
0
0
0
0.5
0
pgc
0.25
0
sc
0.35
0
1
0.1
0
navierstokes
usesELM
0
Lagp2p1
0
Lagp2p1
0
0
0
0.25
1
gls
0.5
1
1
Lagp2p1
isoT
1000
1.2
1e-5
2e-5
7e-2
1
100
hmax_chns*0.5*1.2^2
0
1
{1}
10
1.44
1.92
1.0
1.3
0.09
5.5
0.42
{"x","y"}
{"ref"}
{"shlag(2,'u')","shlag(2,'v')","shlag(1,'p')","shlag(2,'logk')","shlag(2,'logd')","shlag(2,'phi')","shlag(2,'psi')"}
2
on
1
SI
{1}
{"u","v","p","logk","logd","phi","psi"}
sqrt(u^2+v^2)
vx-uy
ux+vy
etaT_chns/sigmad_chns+eta_chns
Cd1_chns*etaT_chns*exp(-logk)*pT_chns-Cd2_chns*rho_chns*exp(logd-logk)
exp(logd)
rho_chns*Cmu_chns*exp(2*logk-logd)
2*ux^2+2*vy^2+(uy+vx)^2
exp(logk)
etaT_chns/sigmak_chns+eta_chns
etaT_chns*exp(-logk)*pT_chns-rho_chns*exp(logd-logk)
-Dlogk_chns*(logkxx+logkyy)+rho_chns*(logkt+u*logkx+v*logky)-Dlogk_chns*(logkx^2+logky^2)-Slogk_chns
-Dlogd_chns*(logdxx+logdyy)+rho_chns*(logdt+u*logdx+v*logdy)-Dlogd_chns*(logdx^2+logdy^2)-Slogd_chns
3*sigma_chns*epsilon_chns/sqrt(8)
psi*lam_chns/epsilon_chns^2
min(max(0.5*(1+phi),0),1)
chi_chns*epsilon_chns^2
0.5*lam_chns*(phix^2+phiy^2)+0.25*lam_chns*(-1+phi^2)^2/epsilon_chns^2+dfdphi_chns*phi
(G_chns-dfdphi_chns)*phix
(G_chns-dfdphi_chns)*phiy
dfdphi_chns*phix
dfdphi_chns*phiy
hmax1_chns
1-Vf_chns
Vf_chns
rho1_chns+(rho2_chns-rho1_chns)*Vf_chns
eta1_chns+(eta2_chns-eta1_chns)*Vf_chns
rho_chns*U_chns*h/eta_chns
rho_chns*(u*ux+v*uy)+px-F_x_chns-(eta_chns+etaT_chns)*(2*uxx+uyy+vxy)-G_chns*phix
rho_chns*(u*vx+v*vy)+py-F_y_chns-(eta_chns+etaT_chns)*(vxx+uyx+2*vyy)-rho_chns*gy_chns-G_chns*phiy
rho_chns*u
rho_chns*v
eta_chns+etaT_chns
rho_chns
nojac(1/max(2*rho_chns*sqrt(emetric(u,v)),48*(eta_chns+etaT_chns)/h^2))
0.5*nojac(if(u^2+v^2
rho_chns*divU_chns
nojac(0.5*h/max(rho_chns*U_chns,6*Dlogk_chns/h))
nojac(0.5*h/max(rho_chns*U_chns,6*Dlogd_chns/h))
0
0
1
0
0
0
0
power
0
0
0
1
0
0
0.5
0.25
sc
0.35
1
0.1
0
Lagp2p1
Lagp2p1
0
0
Lagp2p1
1000
1.2
1e-005
2e-005
0.07000000000000001
1
100
0.72*hmax_chns
0
1
0
-9.810000000000001
U_chns
{0}
{1}
{"u","v","p","logk","logd","phi","psi"}
(eta_chns+etaT_chns)*(2*nx_chns*ux+ny_chns*(uy+vx))
-nx_chns*p+2*nx_chns*(eta_chns+etaT_chns)*ux+ny_chns*(eta_chns+etaT_chns)*(uy+vx)
(eta_chns+etaT_chns)*(nx_chns*(vx+uy)+2*ny_chns*vy)
-ny_chns*p+nx_chns*(eta_chns+etaT_chns)*(vx+uy)+2*ny_chns*(eta_chns+etaT_chns)*vy
-1+2*Vf0_chns
phix*nx_chns+phiy*ny_chns
phix-nx_chns*ndot_gradphi_chns
phiy-ny_chns*ndot_gradphi_chns
sqrt(Wx_chns^2+Wy_chns^2+eps)
epsilon_chns^2*tan(0.5*pi-theta_chns)*W_chns
0
0
0
0
1
0
0
0
0
0
exp(logk)
0.005
10
1
0.01
0.05
0.5*h
rho_chns*dw_chns*Cmu_chns^0.25*exp(0.5*logk)/eta_chns
0
0
0
1
0
1
0
0.5*pi
h
nx
ny
0
0
{0,0,0,0}
{1,2,2,3}
{"u","v","p","logk","logd","phi","psi"}
0
{0,0,0,0}
{1,1,1,1}
10
1.44
1.92
1
1.3
0.09
5.5
0.42
elcplscalar
{"1"}
{"hmax1_chns"}
{"1"}
{"1"}
{0}
{0}
chns
weak
SI
COMSOL 3.5
a
0
603
$Name: $
$Date: 2008/12/03 17:02:19 $
1
geom
on
on
off
on
on
on
on
0.1
t
NavierStokes
cartesian
{"u","v","p","logk","logd","logw","phi","psi","nxw","nyw"}
{"x","y","z"}
chns
CHEM
{"shlag(2,'phi')","shlag(2,'psi')"}
4
2
2
off
_chns
Lagp2p1
init_time
Off
Off
k-epsilon
Off
Off
Off
pfnoncons
ref
off
{"lm1","lm2","lm3","lm4","lm5","lm6","lm7","lm8"}
ideal
0
0
{1,1,1,1}
uv
wall
ntotstress
totstress
cont
U0in
0
0
0
0
1
0
0
0
0
0
0.005
0.005
10
LTIT
open
1
0.01
0.05
dw
h/2
100
0
0
0
1
0
1
velocity
0
0
pi/2
h
{1,2,2,3}
1
1
1
0
0
1
0
0
0
0
power
0
0
0
1
1
0
0
0
0.5
0
pgc
0.25
0
sc
0.35
0
1
0.1
0
navierstokes
usesELM
0
Lagp2p1
0
Lagp2p1
0
0
0
0.25
1
gls
0.5
1
1
Lagp2p1
isoT
1000
1.2
1e-5
2e-5
7e-2
1
100
hmax_chns*0.5*1.2^2
0
1
{1}
10
1.44
1.92
1.0
1.3
0.09
5.5
0.42
{"x","y"}
{"ref"}
{"shlag(2,'phi')","shlag(2,'psi')"}
4
2
2
on
1
SI
{1}
{"phi","psi"}
3*sigma_chns*epsilon_chns/sqrt(8)
psi*lam_chns/epsilon_chns^2
min(max(0.5*(1+phi),0),1)
chi_chns*epsilon_chns^2
0.5*lam_chns*(phix^2+phiy^2)+0.25*lam_chns*(-1+phi^2)^2/epsilon_chns^2+dfdphi_chns*phi
(G_chns-dfdphi_chns)*phix
(G_chns-dfdphi_chns)*phiy
dfdphi_chns*phix
dfdphi_chns*phiy
hmax1_chns
1-Vf_chns
Vf_chns
rho1_chns+(rho2_chns-rho1_chns)*Vf_chns
eta1_chns+(eta2_chns-eta1_chns)*Vf_chns
0
0
1
0
0
0
0
power
0
0
0
1
0
0
0.5
0.25
sc
0.35
1
0.1
0
Lagp2p1
Lagp2p1
0
0
Lagp2p1
1000
1.2
1e-005
2e-005
0.07000000000000001
1
100
0.72*hmax_chns
0
0
-9.810000000000001
{0}
{1}