�� 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 ~ L scriptExtq ~ L scriptNameq ~ xpw t COMSOL Script 1.3t t #COMSOL Reaction Engineering Lab 1.5q ~ t COMSOL 3.5q ~ w �t $Name: $t $Date: 2008/09/19 16:09:48 $xur [Ljava.lang.String;��V�� {G xp t modelinfot xfemt guit g20t fem1t mfileinfouq ~ q ~ t femstructt guistructt drawt geomq ~ xsr com.femlab.api.client.ModelInfo� ^���% L descrq ~ L docURLq ~ [ imaget [Bxpw pt pxuq ~ teyclear xfem clear vrsn vrsn.name = 'COMSOL 3.5'; vrsn.ext = ''; vrsn.major = 0; vrsn.build = 494; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2008/09/19 16:09:48 $'; xfem.version = vrsn; xfem.id = 1; xfem.geomdata = 'geom'; xfem.eqvars = 'on'; xfem.cplbndeq = 'on'; xfem.cplbndsh = 'off'; xfem.drawvalid = 'on'; xfem.geomvalid = 'on'; xfem.solvalid = 'off'; xfem.linshape = 'on'; xfem.linshapetol = 0.1; xfem.meshtime = 't'; clear appl appl.mode.class = 'FlConvCond'; appl.mode.type = 'cartesian'; appl.dim = {'T'}; appl.sdim = {'x','y','z'}; appl.name = 'cc'; appl.shape = {'shlag(2,''T'')'}; appl.gporder = 4; appl.cporder = 2; appl.sshape = 2; appl.border = 'off'; appl.assignsuffix = '_cc'; clear prop prop.elemdefault='Lag2'; prop.analysis='time'; prop.frame='ref'; clear weakconstr weakconstr.value = 'off'; weakconstr.dim = {'lm1'}; prop.weakconstr = weakconstr; prop.constrtype='ideal'; appl.prop = prop; clear bnd bnd.name = {'','','',''}; bnd.q0 = {'0','0','qw','0'}; bnd.T0 = {'273.15','273.15','273.15','273.15'}; bnd.kbnd = {'0','0','0','0'}; bnd.d = {'1','1','1','1'}; bnd.type = {'q','q0','q','cont'}; bnd.relExpr = {{},{},{},{}}; bnd.style = {{{'0'},{'0','255','255'},{'solid'}},{{'0'},{'0','0','0'},{'solid'}}, ... {{'0'},{'0','255','0'},{'solid'}},{{'0'},{'0','0','0'},{'solid'}}}; bnd.ind = [1,2,3,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,3,1,4,4,4,4,4,4]; appl.bnd = bnd; clear equ equ.shape = {[1]}; equ.gporder = {{1}}; equ.cporder = {{1}}; equ.init = {{'273.15'}}; equ.usage = {1}; equ.k = {'mat1_k(T[1/K])[W/(m*K)]'}; equ.ktensor = {{'0.025','0';'0','0.025'}}; equ.ktype = {'iso'}; equ.Dts = {'1'}; equ.rho = {'mat1_rho(T[1/K])[kg/m^3]'}; equ.C = {'mat1_Cp(T[1/K])[J/(kg*K)]'}; equ.Q = {'0'}; equ.u = {'u'}; equ.v = {'v'}; equ.usesELM = {'0'}; equ.usedElement = {'Lag2'}; equ.usesGMG = {'0'}; equ.finestMCelem = {'Lag2'}; equ.stabElmType = {'usesELM'}; equ.idon = {'0'}; equ.delid = {'0.5'}; equ.sdon = {'1'}; equ.sdtype = {'gls'}; equ.delsd = {'0.25'}; equ.cdtypesystem = {'scalar'}; equ.cdon = {'0'}; equ.glim = {'0.01[K]/helem_T_cc'}; equ.gamma = {'1'}; equ.relExpr = {{}}; equ.style = {{{'0'},{'193','193','193'}}}; equ.ind = [1,1,1,1]; appl.equ = equ; xfem.appl{1} = appl; clear appl appl.mode.class = 'FlNavierStokes'; appl.mode.type = 'cartesian'; appl.dim = {'u','v','p','nxw','nyw'}; appl.sdim = {'x','y','z'}; appl.name = 'ns'; appl.shape = {'shlag(2,''u'')','shlag(2,''v'')','shlag(1,''p'')'}; appl.gporder = {4,2}; appl.cporder = {2,1}; appl.sshape = 2; appl.border = 'off'; appl.assignsuffix = '_ns'; clear prop prop.elemdefault='Lagp2p1'; prop.analysis='time'; prop.cornersmoothing='Off'; prop.frame='ref'; clear weakconstr weakconstr.value = 'off'; weakconstr.dim = {'lm2','lm3','lm4'}; prop.weakconstr = weakconstr; prop.constrtype='ideal'; appl.prop = prop; clear pnt pnt.p0 = {'0'}; pnt.pnton = {'0'}; pnt.style = {{{'0'},{'0','0','0'}}}; pnt.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]; appl.pnt = pnt; clear bnd bnd.name = {'','',''}; bnd.type = {'walltype','walltype','int'}; bnd.intype = {'uv','uv','uv'}; bnd.outtype = {'p','p','p'}; bnd.walltype = {'noslip','slip','noslip'}; bnd.opentype = {'ntotstress','ntotstress','ntotstress'}; bnd.stresstype = {'totstress','totstress','totstress'}; bnd.inttype = {'cont','cont','cont'}; bnd.velType = {'U0in','U0in','U0in'}; bnd.u0 = {'0','0','0'}; bnd.v0 = {'0','0','0'}; bnd.p0 = {'0','0','0'}; bnd.f0 = {'0','0','0'}; bnd.Fbnd = {{'0';'0'},{'0';'0'},{'0';'0'}}; bnd.U0in = {'1','1','1'}; bnd.U0out = {'0','0','0'}; bnd.uvw = {'0','0','0'}; bnd.uwall = {'0','0','0'}; bnd.vwall = {'0','0','0'}; bnd.ww = {'0','0','0'}; bnd.style = {{{'0'},{'0','0','0'},{'solid'}},{{'0'},{'0','255','0'},{'solid'}}, ... {{'0'},{'255','0','255'},{'solid'}}}; bnd.ind = [1,1,2,1,2,2,2,2,2,2,2,2,2,2,2,2,2,2,1,2,1,3,3,3,3,3,3]; appl.bnd = bnd; clear equ equ.shape = {[1;2;3]}; equ.gporder = {{1;1;2}}; equ.cporder = {{1;1;2}}; equ.init = {{'0';'0';'0';'0';'0'}}; equ.usage = {1}; equ.rho = {'mat1_rho(T[1/K])[kg/m^3]'}; equ.eta = {'mat1_eta(T[1/K])[Pa*s]'}; equ.F_x = {'0'}; equ.F_y = {'alpha*g*rho_cc*(T-T_m)*(T>=T_m)'}; equ.idon = {'1'}; equ.delid = {'0.7'}; equ.sdon = {'1'}; equ.sdtype = {'gls'}; equ.delsd = {'0.25'}; equ.cdon = {'1'}; equ.cdtype = {'sc'}; equ.delcd = {'0.35'}; equ.pson = {'0'}; equ.delps = {'1'}; equ.ck = {'0.3'}; equ.impreson = {'0'}; equ.artstabtype = {'navierstokes'}; equ.stabElmType = {'usesELM'}; equ.usesELM = {'0'}; equ.usedElement = {'Lagp2p1'}; equ.usesGMG = {'0'}; equ.finestMCelem = {'Lagp2p1'}; equ.gls = {'0'}; equ.g = {'0'}; equ.gijgij = {'0'}; equ.style = {{{'0'},{'193','193','193'}}}; equ.ind = [1,1,1,1]; appl.equ = equ; appl.var = {'visc_vel_fact','10'}; xfem.appl{2} = appl; clear appl appl.mode.class = 'FlPDEWBoundary'; appl.mode.type = 'cartesian'; appl.dim = {'lm','lm_t'}; appl.sdim = {'x','y','z'}; appl.name = 'wb'; appl.shape = {'shlag(1,''lm'')','shlag(2,''lm'')'}; appl.gporder = {2,4}; appl.cporder = {1,2}; appl.sshape = 1; appl.border = 'off'; appl.assignsuffix = '_wb'; clear prop prop.elemdefault='Lag2'; prop.wave='off'; prop.frame='ref'; clear weakconstr weakconstr.value = 'off'; weakconstr.dim = {'lm5','lm6'}; prop.weakconstr = weakconstr; appl.prop = prop; clear pnt pnt.weak = {{'0'}}; pnt.dweak = {{'0'}}; pnt.constr = {{'0'}}; pnt.constrf = {{'0'}}; pnt.constrtype = {'ideal'}; pnt.style = {{{'0'},{'0','0','0'}}}; pnt.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]; appl.pnt = pnt; clear bnd bnd.shape = {[1],[1],[2]}; bnd.gporder = {{1},{1},{2}}; bnd.cporder = {{1},{1},{2}}; bnd.init = {{'0';'0'},{'0';'0'},{'0';'0'}}; bnd.usage = {0,1,0}; bnd.name = {'','',''}; bnd.dinit = {{'0';'0'},{'0';'0'},{'0';'0'}}; bnd.weak = {{'0'},{'lm_test*(eta_ns*uy-d_gamma2*Tx)+u_test*lm'},{'0'}}; bnd.dweak = {{'0'},{'0'},{'0'}}; bnd.constr = {{'0'},{'0'},{'0'}}; bnd.constrf = {{'0'},{'0'},{'0'}}; bnd.constrtype = {'ideal','ideal','ideal'}; bnd.style = {{{'0'},{'0','0','0'},{'solid'}},{{'0'},{'0','0','255'},{'solid'}}, ... {{'0'},{'0','255','0'},{'solid'}}}; bnd.ind = [1,1,2,1,2,2,2,2,2,2,2,2,2,2,2,2,2,2,1,2,1,3,3,3,3,3,3]; appl.bnd = bnd; xfem.appl{3} = appl; xfem.geom = flbinary('fem1','geom','cc_iso_wb(A)_mesh4(T)_xxcoarse2_low_a_4.mph'); xfem.sdim = {'x','y'}; xfem.frame = {'ref'}; xfem.shape = {'shlag(2,''T'')','shlag(2,''u'')','shlag(2,''v'')','shlag(1,''p'')','shlag(1,''lm'')'}; xfem.gporder = {4,2}; xfem.cporder = {2,1}; xfem.sshape = 2; xfem.simplify = 'on'; xfem.border = 1; xfem.form = 'general'; clear units; units.basesystem = 'SI'; xfem.units = units; clear equ equ.shape = {[1;2;3;4]}; equ.gporder = {{1;1;1;2}}; equ.cporder = {{1;1;1;2}}; equ.init = {{'273.15';'0';'0';'0'}}; equ.dinit = {{'0';'0';'0';'0'}}; equ.weak = {{'-test(nojac(beta_T_x_cc)*Tx+nojac(beta_T_y_cc)*Ty-nojac(kmean_T_cc)*Txx-nojac(kmean_T_cc)*Tyy)*taue_T_cc*res_T_cc'; ... 'gls_ns+shock_capt_ns';'0';'0'}}; equ.dweak = {{'0';'0';'0';'0'}}; equ.constr = {{'0';'0';'0';'0'}}; equ.constrf = {{'0';'0';'0';'0'}}; equ.c = {{{'-diff(-k_T_cc*Tx,Tx)','-diff(-k_T_cc*Tx,Ty)';'-diff(-k_T_cc*Ty,Tx)', ... '-diff(-k_T_cc*Ty,Ty)'},{'-diff(-k_T_cc*Tx,ux)','-diff(-k_T_cc*Tx,uy)'; ... '-diff(-k_T_cc*Ty,ux)','-diff(-k_T_cc*Ty,uy)'},{'-diff(-k_T_cc*Tx,vx)', ... '-diff(-k_T_cc*Tx,vy)';'-diff(-k_T_cc*Ty,vx)','-diff(-k_T_cc*Ty,vy)'},{'-diff(-k_T_cc*Tx,px)', ... '-diff(-k_T_cc*Tx,py)';'-diff(-k_T_cc*Ty,px)','-diff(-k_T_cc*Ty,py)'};{'-diff(-2*eta_ns*ux+p,Tx)', ... '-diff(-2*eta_ns*ux+p,Ty)';'-diff(-eta_ns*(uy+vx),Tx)','-diff(-eta_ns*(uy+vx),Ty)'}, ... {'-diff(-2*eta_ns*ux+p,ux)','-diff(-2*eta_ns*ux+p,uy)';'-diff(-eta_ns*(uy+vx),ux)', ... '-diff(-eta_ns*(uy+vx),uy)'},{'-diff(-2*eta_ns*ux+p,vx)','-diff(-2*eta_ns*ux+p,vy)'; ... '-diff(-eta_ns*(uy+vx),vx)','-diff(-eta_ns*(uy+vx),vy)'},{'-diff(-2*eta_ns*ux+p,px)', ... '-diff(-2*eta_ns*ux+p,py)';'-diff(-eta_ns*(uy+vx),px)','-diff(-eta_ns*(uy+vx),py)'}; ... {'-diff(-eta_ns*(vx+uy),Tx)','-diff(-eta_ns*(vx+uy),Ty)';'-diff(-2*eta_ns*vy+p,Tx)', ... '-diff(-2*eta_ns*vy+p,Ty)'},{'-diff(-eta_ns*(vx+uy),ux)','-diff(-eta_ns*(vx+uy),uy)'; ... '-diff(-2*eta_ns*vy+p,ux)','-diff(-2*eta_ns*vy+p,uy)'},{'-diff(-eta_ns*(vx+uy),vx)', ... '-diff(-eta_ns*(vx+uy),vy)';'-diff(-2*eta_ns*vy+p,vx)','-diff(-2*eta_ns*vy+p,vy)'}, ... {'-diff(-eta_ns*(vx+uy),px)','-diff(-eta_ns*(vx+uy),py)';'-diff(-2*eta_ns*vy+p,px)', ... '-diff(-2*eta_ns*vy+p,py)'};{'0'},{'0'},{'0'},{'0'}}}; equ.a = {{'-diff(rho_cc*C_cc*(-u_cc*Tx-v_cc*Ty),T)','-diff(rho_cc*C_cc*(-u_cc*Tx-v_cc*Ty),u)', ... '-diff(rho_cc*C_cc*(-u_cc*Tx-v_cc*Ty),v)','-diff(rho_cc*C_cc*(-u_cc*Tx-v_cc*Ty),p)'; ... '-diff(-rho_ns*(u*ux+v*uy),T)','-diff(-rho_ns*(u*ux+v*uy),u)','-diff(-rho_ns*(u*ux+v*uy),v)', ... '-diff(-rho_ns*(u*ux+v*uy),p)';'-diff(F_y_ns-rho_ns*(u*vx+v*vy),T)','-diff(F_y_ns-rho_ns*(u*vx+v*vy),u)', ... '-diff(F_y_ns-rho_ns*(u*vx+v*vy),v)','-diff(F_y_ns-rho_ns*(u*vx+v*vy),p)'; ... '-diff(-divU_ns,T)','-diff(-divU_ns,u)','-diff(-divU_ns,v)','-diff(-divU_ns,p)'}}; equ.f = {{'rho_cc*C_cc*(-u_cc*Tx-v_cc*Ty)';'-rho_ns*(u*ux+v*uy)';'F_y_ns-rho_ns*(u*vx+v*vy)'; ... '-divU_ns'}}; equ.ea = {{'0','0','0','0';'0','0','0','0';'0','0','0','0';'0','0','0', ... '0'}}; equ.da = {{'rho_cc*C_cc','0','0','0';'0','rho_ns','0','0';'0','0','rho_ns', ... '0';'0','0','0','0'}}; equ.al = {{{'-diff(-k_T_cc*Tx,T)';'-diff(-k_T_cc*Ty,T)'},{'-diff(-k_T_cc*Tx,u)'; ... '-diff(-k_T_cc*Ty,u)'},{'-diff(-k_T_cc*Tx,v)';'-diff(-k_T_cc*Ty,v)'},{'-diff(-k_T_cc*Tx,p)'; ... '-diff(-k_T_cc*Ty,p)'};{'-diff(-2*eta_ns*ux+p,T)';'-diff(-eta_ns*(uy+vx),T)'}, ... {'-diff(-2*eta_ns*ux+p,u)';'-diff(-eta_ns*(uy+vx),u)'},{'-diff(-2*eta_ns*ux+p,v)'; ... '-diff(-eta_ns*(uy+vx),v)'},{'-diff(-2*eta_ns*ux+p,p)';'-diff(-eta_ns*(uy+vx),p)'}; ... {'-diff(-eta_ns*(vx+uy),T)';'-diff(-2*eta_ns*vy+p,T)'},{'-diff(-eta_ns*(vx+uy),u)'; ... '-diff(-2*eta_ns*vy+p,u)'},{'-diff(-eta_ns*(vx+uy),v)';'-diff(-2*eta_ns*vy+p,v)'}, ... {'-diff(-eta_ns*(vx+uy),p)';'-diff(-2*eta_ns*vy+p,p)'};{'0';'0'},{'0';'0'}, ... {'0';'0'},{'0';'0'}}}; equ.be = {{{'-diff(rho_cc*C_cc*(-u_cc*Tx-v_cc*Ty),Tx)';'-diff(rho_cc*C_cc*(-u_cc*Tx-v_cc*Ty),Ty)'}, ... {'-diff(rho_cc*C_cc*(-u_cc*Tx-v_cc*Ty),ux)';'-diff(rho_cc*C_cc*(-u_cc*Tx-v_cc*Ty),uy)'}, ... {'-diff(rho_cc*C_cc*(-u_cc*Tx-v_cc*Ty),vx)';'-diff(rho_cc*C_cc*(-u_cc*Tx-v_cc*Ty),vy)'}, ... {'-diff(rho_cc*C_cc*(-u_cc*Tx-v_cc*Ty),px)';'-diff(rho_cc*C_cc*(-u_cc*Tx-v_cc*Ty),py)'}; ... {'-diff(-rho_ns*(u*ux+v*uy),Tx)';'-diff(-rho_ns*(u*ux+v*uy),Ty)'},{'-diff(-rho_ns*(u*ux+v*uy),ux)'; ... '-diff(-rho_ns*(u*ux+v*uy),uy)'},{'-diff(-rho_ns*(u*ux+v*uy),vx)';'-diff(-rho_ns*(u*ux+v*uy),vy)'}, ... {'-diff(-rho_ns*(u*ux+v*uy),px)';'-diff(-rho_ns*(u*ux+v*uy),py)'};{'-diff(F_y_ns-rho_ns*(u*vx+v*vy),Tx)'; ... '-diff(F_y_ns-rho_ns*(u*vx+v*vy),Ty)'},{'-diff(F_y_ns-rho_ns*(u*vx+v*vy),ux)'; ... '-diff(F_y_ns-rho_ns*(u*vx+v*vy),uy)'},{'-diff(F_y_ns-rho_ns*(u*vx+v*vy),vx)'; ... '-diff(F_y_ns-rho_ns*(u*vx+v*vy),vy)'},{'-diff(F_y_ns-rho_ns*(u*vx+v*vy),px)'; ... '-diff(F_y_ns-rho_ns*(u*vx+v*vy),py)'};{'-diff(-divU_ns,Tx)';'-diff(-divU_ns,Ty)'}, ... {'-diff(-divU_ns,ux)';'-diff(-divU_ns,uy)'},{'-diff(-divU_ns,vx)';'-diff(-divU_ns,vy)'}, ... {'-diff(-divU_ns,px)';'-diff(-divU_ns,py)'}}}; equ.ga = {{{'-k_T_cc*Tx';'-k_T_cc*Ty'};{'-2*eta_ns*ux+p';'-eta_ns*(uy+vx)'}; ... {'-eta_ns*(vx+uy)';'-2*eta_ns*vy+p'};{'0';'0'}}}; equ.sshape = {[1]}; equ.sshapedim = {{1;1;1;1}}; equ.ind = [1,1,1,1]; equ.dim = {'T','u','v','p'}; equ.var = {'grad_T_x_cc',{'Tx'}, ... 'dflux_T_x_cc',{'-kxx_T_cc*Tx-kxy_T_cc*Ty'}, ... 'cflux_T_x_cc',{'rho_T_cc*C_T_cc*T*u_T_cc'}, ... 'tflux_T_x_cc',{'dflux_T_x_cc+cflux_T_x_cc'}, ... 'grad_T_y_cc',{'Ty'}, ... 'dflux_T_y_cc',{'-kyx_T_cc*Tx-kyy_T_cc*Ty'}, ... 'cflux_T_y_cc',{'rho_T_cc*C_T_cc*T*v_T_cc'}, ... 'tflux_T_y_cc',{'dflux_T_y_cc+cflux_T_y_cc'}, ... 'beta_T_x_cc',{'rho_T_cc*C_T_cc*u_T_cc'}, ... 'beta_T_y_cc',{'rho_T_cc*C_T_cc*v_T_cc'}, ... 'grad_T_cc',{'sqrt(grad_T_x_cc^2+grad_T_y_cc^2)'}, ... 'dflux_T_cc',{'sqrt(dflux_T_x_cc^2+dflux_T_y_cc^2)'}, ... 'cflux_T_cc',{'sqrt(cflux_T_x_cc^2+cflux_T_y_cc^2)'}, ... 'tflux_T_cc',{'sqrt(tflux_T_x_cc^2+tflux_T_y_cc^2)'}, ... 'cellPe_T_cc',{'h*sqrt(beta_T_x_cc^2+beta_T_y_cc^2)/Dm_T_cc'}, ... 'kmean_T_cc',{'k_cc'}, ... 'taue_T_cc',{'nojac(1/max(2*C_T_cc*rho_T_cc*sqrt(emetric(u_cc,v_cc))/gamma_cc,48*kmean_T_cc/h^2))'}, ... 'helem_T_cc',{'h'}, ... 'Dm_T_cc',{'rho_T_cc^2*C_T_cc^2*(kxx_T_cc*u_T_cc^2+kxy_T_cc*u_T_cc*v_T_cc+kyx_T_cc*v_T_cc*u_T_cc+kyy_T_cc*v_T_cc^2)/((rho_T_cc*C_T_cc*u_T_cc)^2+(rho_T_cc*C_T_cc*v_T_cc)^2+eps)'}, ... 'res_T_cc',{'-kxx_T_cc*Txx-kxy_T_cc*Txy-kyx_T_cc*Tyx-kyy_T_cc*Tyy+Tx*rho_T_cc*C_T_cc*u_T_cc+Ty*rho_T_cc*C_T_cc*v_T_cc-Q_T_cc+rho_cc*C_cc*Tt'}, ... 'da_T_cc',{'Dts_T_cc*rho_T_cc*C_T_cc'}, ... 'U_ns',{'sqrt(u^2+v^2)'}, ... 'V_ns',{'vx-uy'}, ... 'divU_ns',{'ux+vy'}, ... 'cellRe_ns',{'rho_ns*U_ns*h/eta_ns'}, ... 'res_u_ns',{'rho_ns*(ut+u*ux+v*uy)+px-F_x_ns-eta_ns*(2*uxx+uyy+vxy)'}, ... 'res_v_ns',{'rho_ns*(vt+u*vx+v*vy)+py-F_y_ns-eta_ns*(vxx+uyx+2*vyy)'}, ... 'beta_x_ns',{'rho_ns*u'}, ... 'beta_y_ns',{'rho_ns*v'}, ... 'Dm_ns',{'eta_ns'}, ... 'da_ns',{'rho_ns'}, ... 'taum_ns',{'nojac(1/max(2*rho_ns*sqrt(emetric(u,v)),48*eta_ns/h^2))'}, ... 'tauc_ns',{'0.5*nojac(if(u^2+v^2
=T_m)'}, ... 'delid_ns',{'0.7'}, ... 'delsd_ns',{'0.25'}, ... 'cdtype_ns',{'sc'}, ... 'delcd_ns',{'0.35'}, ... 'delps_ns',{'1'}, ... 'ck_ns',{'0.3'}, ... 'impreson_ns',{'0'}, ... 'usedElement_ns',{'Lagp2p1'}, ... 'finestMCelem_ns',{'Lagp2p1'}, ... 'gls_ns',{'-test(nojac(rho_ns)*(nojac(u)*ux+nojac(v)*uy)+px-nojac(eta_ns)*(2*uxx+uyy+vxy))*taum_ns*res_u_ns-test(nojac(rho_ns)*(nojac(u)*vx+nojac(v)*vy)+py-nojac(eta_ns)*(vxx+uyx+2*vyy))*taum_ns*res_v_ns-test(divU_ns)*tauc_ns*res_p_ns'}, ... 'g_ns',{'0'}, ... 'gijgij_ns',{'(diff(x,xi1)^2+diff(x,xi2)^2)^2+2*(diff(y,xi1)*diff(x,xi1)+diff(y,xi2)*diff(x,xi2))^2+(diff(y,xi1)^2+diff(y,xi2)^2)^2'}, ... 'g11_ns',{'diff(x,xi1)^2+diff(x,xi2)^2'}, ... 'g21_ns',{'diff(y,xi1)*diff(x,xi1)+diff(y,xi2)*diff(x,xi2)'}, ... 'g22_ns',{'diff(y,xi1)^2+diff(y,xi2)^2'}, ... 'shock_capt_ns',{'-nojac(0.5*min(ck_ns*sqrt((tauc_ns*res_p_ns/(nojac(rhoint_g1_ns/vol_g1_ns)*nojac(sqrt(u2int_g1_ns/vol_g1_ns))))^2+(taum_ns*res_u_ns)^2+(taum_ns*res_v_ns)^2),U_ns)*min(gijgij_ns^(-0.25),rho_ns*U_ns/(48*eta_ns)))*rho_ns*(g11_ns*(ux*test(ux)+vx*test(vx))+g21_ns*(uy*test(ux)+vy*test(vx))+g21_ns*(ux*test(uy)+vx*test(vy))+g22_ns*(uy*test(uy)+vy*test(vy)))'}}; equ.expr = {}; equ.bnd.weak = {{'0';'0';'0';'0'}}; equ.bnd.gporder = {{1;1;1;2}}; equ.bnd.ind = [1,1,1,1]; equ.bnd.var = {}; equ.bnd.expr = {}; equ.lock = [0,0,0,0]; equ.mlock = {[0,0,0,0]}; xfem.equ = equ; clear bnd bnd.weak = {{'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'lm_test*(eta_ns*uy-d_gamma2*Tx)+u_test*lm'}, ... {'0';'0';'0';'0';'0'}}; bnd.dweak = {{'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0'},{'0';'0';'0';'0'; ... '0'}}; bnd.constr = {{'0';'-u';'-v';'0';'0'},{'0';'-u*nx_ns-v*ny_ns';'0';'0';'0'}, ... {'0';'0';'0';'0';'0'}}; bnd.constrf = {{'0';'test(-u)';'test(-v)';'0';'0'},{'0';'test(-u*nx_ns-v*ny_ns)'; ... '0';'0';'0'},{'0';'0';'0';'0';'0'}}; bnd.q = {{'0','0','0','0','0';'0','0','0','0','0';'0','0','0','0','0'; ... '0','0','0','0','0';'0','0','0','0','0'},{'-diff(q0_cc,T)','-diff(q0_cc,u)', ... '-diff(q0_cc,v)','-diff(q0_cc,p)','-diff(q0_cc,lm)';'0','0','0','0','0'; ... '0','0','0','0','0';'0','0','0','0','0';'0','0','0','0','0'},{'0','0','0', ... '0','0';'0','0','0','0','0';'0','0','0','0','0';'0','0','0','0','0';'0', ... '0','0','0','0'}}; bnd.h = {{'0','0','0','0','0';'0','0','0','0','0';'0','0','0','0','0'; ... '0','0','0','0','0';'0','0','0','0','0'},{'0','0','0','0','0';'0','0','0', ... '0','0';'0','0','0','0','0';'0','0','0','0','0';'0','0','0','0','0'},{'0', ... '0','0','0','0';'0','0','0','0','0';'0','0','0','0','0';'0','0','0','0', ... '0';'0','0','0','0','0'}}; bnd.g = {{'0';'0';'0';'0';'0'},{'q0_cc';'0';'0';'0';'0'},{'0';'0';'0'; ... '0';'0'}}; bnd.r = {{'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0'},{'0';'0';'0';'0'; ... '0'}}; bnd.shape = {[1;2;3;4],[1;2;3;4;5],[1;2;3;4]}; bnd.sshape = {[1],[1],[1]}; bnd.sshapedim = {{1;1;1;1;1},{1;1;1;1;1},{1;1;1;1;1}}; bnd.gporder = {{1;1;1;2;2},{1;1;1;2;2},{1;1;1;2;1}}; bnd.cporder = {{1;1;1;2;2},{1;1;1;2;2},{1;1;1;2;1}}; bnd.init = {{'';'';'';'';'0'},{'';'';'';'';'0'},{'';'';'';'';'0'}}; bnd.dinit = {{'';'';'';'';'0'},{'';'';'';'';'0'},{'';'';'';'';'0'}}; bnd.ind = [1,1,2,1,2,2,2,2,2,2,2,2,2,2,2,2,2,2,1,2,1,3,3,3,3,3,3]; bnd.dim = {'T','u','v','p','lm'}; bnd.var = {'ndflux_T_cc',{'nx_cc*dflux_T_x_cc+ny_cc*dflux_T_y_cc','nx_cc*dflux_T_x_cc+ny_cc*dflux_T_y_cc', ... 'nx_cc*dflux_T_x_cc+ny_cc*dflux_T_y_cc'}, ... 'ncflux_T_cc',{'nx_cc*cflux_T_x_cc+ny_cc*cflux_T_y_cc','nx_cc*cflux_T_x_cc+ny_cc*cflux_T_y_cc', ... 'nx_cc*cflux_T_x_cc+ny_cc*cflux_T_y_cc'}, ... 'ntflux_T_cc',{'nx_cc*tflux_T_x_cc+ny_cc*tflux_T_y_cc','nx_cc*tflux_T_x_cc+ny_cc*tflux_T_y_cc', ... 'nx_cc*tflux_T_x_cc+ny_cc*tflux_T_y_cc'}, ... 'K_x_ns',{'eta_ns*(2*nx_ns*ux+ny_ns*(uy+vx))','eta_ns*(2*nx_ns*ux+ny_ns*(uy+vx))', ... 'eta_ns*(2*nx_ns*ux+ny_ns*(uy+vx))'}, ... 'T_x_ns',{'-nx_ns*p+2*nx_ns*eta_ns*ux+ny_ns*eta_ns*(uy+vx)','-nx_ns*p+2*nx_ns*eta_ns*ux+ny_ns*eta_ns*(uy+vx)', ... '-nx_ns*p+2*nx_ns*eta_ns*ux+ny_ns*eta_ns*(uy+vx)'}, ... 'K_y_ns',{'eta_ns*(nx_ns*(vx+uy)+2*ny_ns*vy)','eta_ns*(nx_ns*(vx+uy)+2*ny_ns*vy)', ... 'eta_ns*(nx_ns*(vx+uy)+2*ny_ns*vy)'}, ... 'T_y_ns',{'-ny_ns*p+nx_ns*eta_ns*(vx+uy)+2*ny_ns*eta_ns*vy','-ny_ns*p+nx_ns*eta_ns*(vx+uy)+2*ny_ns*eta_ns*vy', ... '-ny_ns*p+nx_ns*eta_ns*(vx+uy)+2*ny_ns*eta_ns*vy'},'q0_cc',{'0','qw','0'}, ... 'T0_cc',{'273.15','273.15','273.15'}, ... 'kbnd_cc',{'0','0','0'}, ... 'd_cc',{'1','1','1'}, ... 'nx_cc',{'nx','nx','nx'}, ... 'ny_cc',{'ny','ny','ny'}, ... 'N_T_cc',{'0','qw','0'}, ... 'c0_T_cc',{'273.15','273.15','273.15'}, ... 'u0_ns',{'0','0','0'}, ... 'v0_ns',{'0','0','0'}, ... 'p0_ns',{'0','0','0'}, ... 'f0_ns',{'0','0','0'}, ... 'U0in_ns',{'1','1','1'}, ... 'U0out_ns',{'0','0','0'}, ... 'uvw_ns',{'0','0','0'}, ... 'uwall_ns',{'0','0','0'}, ... 'vwall_ns',{'0','0','0'}, ... 'ww_ns',{'0','0','0'}, ... 'nx_ns',{'nx','nx','nx'}, ... 'ny_ns',{'ny','ny','ny'}, ... 'Fbndx_ns',{'0','0','0'}, ... 'Fbndy_ns',{'0','0','0'}}; bnd.expr = {}; bnd.lock = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; bnd.mlock = {[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]}; xfem.bnd = bnd; clear pnt pnt.weak = {{'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0'}}; pnt.dweak = {{'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0'}}; pnt.constr = {{'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0'}}; pnt.constrf = {{'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0'}}; pnt.shape = {[1;2;3;4],[1;2;3;4;5]}; pnt.sshape = {[1],[1]}; pnt.sshapedim = {{1;1;1;1;1},{1;1;1;1;1}}; pnt.init = {{'';'';'';'';''},{'';'';'';'';''}}; pnt.dinit = {{'';'';'';'';''},{'';'';'';'';''}}; pnt.ind = [1,2,1,2,2,2,2,2,2,2,2,1,1,1,2,2,2,2,2,2,1,2,1,2]; pnt.dim = {'T','u','v','p','lm'}; pnt.var = {'p0_ns',{'0','0'}}; pnt.expr = {}; pnt.lock = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; pnt.mlock = {[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]}; xfem.pnt = pnt; xfem.var = {'visc_vel_fact_ns','10'}; xfem.expr = {'qw','2*P*Acon/(rf^2*pi)*exp(-2*x^2/rf^2)*flc2hs(t[1/s]-0.05,0.05)', ... 'd_gamma','-A-R*gs*log(1+K*a)[1/kg]-(K*a/(1+K*a))*gs*(delta_H-HM)/T', ... 'd_gamma2','diff(gamma,T)', ... 'K','k*exp(-delta_H/(R*T)[1/kg])', ... 'gamma','gamma_m-A*(T-T_m)-R*T*gs*(log(1+k*a*exp(-delta_H/(R*T)[1/kg])))[1/kg]'}; clear descr descr.expr= {'qw','Wärmequelle'}; xfem.descr = descr; clear elemmph clear elem elem.elem = 'elcplscalar'; elem.g = {'1'}; src = cell(1,1); clear equ equ.expr = {{'nojac(rho_ns)'},{'nojac(1)'},{ ... 'nojac(max((+u^2+v^2),(eta_ns/rho_ns/h)^2))'},{ ... 'nojac(-1/max(emetric(u,v),(visc_vel_fact_ns*eta_ns/rho_ns)^2*emetric2))'}}; equ.ipoints = {{'0'},{'0'},{'0'},{'0'}}; equ.frame = {{'ref'},{'ref'},{'ref'},{'ref'}}; equ.ind = {{'1','2','3','4'}}; src{1} = {{},{},equ}; elem.src = src; geomdim = cell(1,1); geomdim{1} = {}; elem.geomdim = geomdim; elem.var = {'rhoint_g1_ns','vol_g1_ns','u2int_g1_ns','mdt2_cfl_g1_ns'}; elem.global = {'1','2','3','4'}; elem.maxvars = {'4'}; elemmph{1} = elem; clear elem elem.elem = 'elconst'; elem.var = {'dt_cfl_g1_ns','nojac(sqrt(-mdt2_cfl_g1_ns))'}; elemmph{2} = elem; xfem.elemmph = elemmph; clear draw draw.p.objs = {}; draw.p.name = {}; draw.c.objs = {}; draw.c.name = {}; draw.s.objs = {flbinary('g20','draw','cc_iso_wb(A)_mesh4(T)_xxcoarse2_low_a_4.mph')}; draw.s.name = {'CO1'}; xfem.draw = draw; xfem.const = {'P','500[W]','Acon','0.36','rf','1.5[mm]','k_const','50[W/(m*K)]','rho_const','7870[kg/m^3]','visc_const','0.0055[Pa*s]','cp_const','600[J/(kg*K)]','g','9.81[m/s^2]','alpha','12e-6[1/K]','T_m','1806[K]','A','4.3e-4[N/(m*K)]','R','8.3145[kJ/(mol*K)]','delta_H','-166.2e3[kJ/(kg*mol)]','k','0.00318','gs','1.3e-8[kg*mol/m^2]','HM','-6.2e3[kJ/(kg*mol)]','a','0.0010','gamma_m','1.943[N/m]'}; xfem.globalexpr = {}; clear fcns xfem.functions = {}; descr = xfem.descr; descr.const= {'g','Erdbeschleunigung','a','Grenzflächenaktivität','HM','partielle molare enthalpie','rf','Fokusierradius','Acon','Absorptionskoeffizient','k','entropiekonstante','delta_H','molare adsorptionsenthalpie','cp_const','konstante spez. Wärmekapazität','P','Abgestrahle Leistung','alpha','Temperaturausdehnungskoeffizient','R','universelle Gaskonstante','k_const','konstante Wärmeleitfähigkeit','rho_const','konstante Diche','gamma_m','Oberflächenspannung bei Raumtemperatur','visc_const','konstante Viskosität','gs','Sättigungswert'}; xfem.descr = descr; xfem.xmcases = [0]; xfem.mcases = [0]; flbinary clear; xfem.rulingmode = 'cc'; xfem.solform = 'weak'; clear lib lib.mat{1}.name='SR235JG'; lib.mat{1}.varname='mat1'; lib.mat{1}.variables.eta='eta(T[1/K])[Pa*s]'; lib.mat{1}.variables.C='Cp(T[1/K])[J/(kg*K)]'; lib.mat{1}.variables.rho='rho(T[1/K])[kg/m^3]'; lib.mat{1}.variables.k='k(T[1/K])[W/(m*K)]'; clear fcns fcns{1}.type='inline'; fcns{1}.name='rho(T)'; fcns{1}.expr='7870-855*flc2hs(T-1806,50)'; fcns{1}.dexpr={'diff(7870-855*flc2hs(T-1806,50),T)'}; fcns{2}.type='inline'; fcns{2}.name='eta(T)'; fcns{2}.expr='500-499.9945*flc2hs(T-1806,50)'; fcns{2}.dexpr={'diff(500-499.9945*flc2hs(T-1806,50),T)'}; fcns{3}.type='interp'; fcns{3}.name='k'; fcns{3}.method='cubic'; fcns{3}.extmethod='extrap'; fcns{3}.x={'293','373','473','573','673','773','873','5000'}; fcns{3}.data={'57','57','54','50','47','42','37','37'}; fcns{4}.type='interp'; fcns{4}.name='Cp'; fcns{4}.method='cubic'; fcns{4}.extmethod='extrap'; fcns{4}.x={'293','373','473','573','673','773','873','5000'}; fcns{4}.data={'430','470','500','540','580','660','780','780'}; lib.mat{1}.functions = fcns; lib.mat{2}.name='SR235JG_1'; lib.mat{2}.varname='mat2'; lib.mat{2}.variables.eta='eta(T[1/K])[Pa*s]'; lib.mat{2}.variables.C='Cp(T[1/K])[J/(kg*K)]'; lib.mat{2}.variables.rho='rho(T[1/K])[kg/m^3]'; lib.mat{2}.variables.k='k(T[1/K])[W/(m*K)]'; clear fcns fcns{1}.type='inline'; fcns{1}.name='rho(T)'; fcns{1}.expr='7870-855*flc2hs(T-1806,50)'; fcns{1}.dexpr={'diff(7870-855*flc2hs(T-1806,50),T)'}; fcns{2}.type='inline'; fcns{2}.name='eta(T)'; fcns{2}.expr='500-499.9945*flc2hs(T-1806,50)'; fcns{2}.dexpr={'diff(500-499.9945*flc2hs(T-1806,50),T)'}; fcns{3}.type='interp'; fcns{3}.name='k'; fcns{3}.method='cubic'; fcns{3}.extmethod='extrap'; fcns{3}.x={'293','373','473','573','673','773','873','5000'}; fcns{3}.data={'57','57','54','50','47','42','37','37'}; fcns{4}.type='interp'; fcns{4}.name='Cp'; fcns{4}.method='cubic'; fcns{4}.extmethod='extrap'; fcns{4}.x={'293','373','473','573','673','773','873','5000'}; fcns{4}.data={'430','470','500','540','580','660','780','780'}; lib.mat{2}.functions = fcns; lib.mat{3}.name='SR235JG_2'; lib.mat{3}.varname='mat3'; lib.mat{3}.variables.eta='eta(T[1/K])[Pa*s]'; lib.mat{3}.variables.C='Cp(T[1/K])[J/(kg*K)]'; lib.mat{3}.variables.rho='rho(T[1/K])[kg/m^3]'; lib.mat{3}.variables.k='k(T[1/K])[W/(m*K)]'; clear fcns fcns{1}.type='inline'; fcns{1}.name='rho(T)'; fcns{1}.expr='7870-855*flc2hs(T-1806,50)'; fcns{1}.dexpr={'diff(7870-855*flc2hs(T-1806,50),T)'}; fcns{2}.type='inline'; fcns{2}.name='eta(T)'; fcns{2}.expr='500-499.9945*flc2hs(T-1806,50)'; fcns{2}.dexpr={'diff(500-499.9945*flc2hs(T-1806,50),T)'}; fcns{3}.type='interp'; fcns{3}.name='k'; fcns{3}.method='cubic'; fcns{3}.extmethod='extrap'; fcns{3}.x={'293','373','473','573','673','773','873','5000'}; fcns{3}.data={'57','57','54','50','47','42','37','37'}; fcns{4}.type='interp'; fcns{4}.name='Cp'; fcns{4}.method='cubic'; fcns{4}.extmethod='extrap'; fcns{4}.x={'293','373','473','573','673','773','873','5000'}; fcns{4}.data={'430','470','500','540','580','660','780','780'}; lib.mat{3}.functions = fcns; xfem.lib = lib; clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; uq ~ t�`gui.solvemodel.toutcomp='off'; gui.solvemodel.currsolver='time'; gui.solvemodel.solveroption='plotcallback'; gui.solvemodel.postsolver='stationary'; gui.solvemodel.nonlin='auto'; gui.solvemodel.ntol='1.0E-6'; gui.solvemodel.maxiter='25'; gui.solvemodel.segterm='tol'; gui.solvemodel.maxsegiter='100'; gui.solvemodel.segiter='1'; gui.solvemodel.manualdamp='off'; gui.solvemodel.damping='on'; gui.solvemodel.hnlin='off'; gui.solvemodel.initstep='1.0'; gui.solvemodel.minstep='1.0E-4'; gui.solvemodel.rstep='10.0'; gui.solvemodel.useaugsolver='off'; gui.solvemodel.autoaugcomp='on'; gui.solvemodel.augcomp=''; gui.solvemodel.augtol='0.0010'; gui.solvemodel.augmaxiter='25'; gui.solvemodel.augsolver='lumped'; gui.solvemodel.nlsolver='automatic'; gui.solvemodel.timenonlin='auto'; gui.solvemodel.useratelimit='on'; gui.solvemodel.timentolfact='1'; gui.solvemodel.timemaxiter='4'; gui.solvemodel.timesegterm='tol'; gui.solvemodel.timemaxsegiter='10'; gui.solvemodel.timesegiter='1'; gui.solvemodel.timemanualdamp='off'; gui.solvemodel.timedtech='const'; gui.solvemodel.timedamp='1.0'; gui.solvemodel.timejtech='once'; gui.solvemodel.timeinitstep='1.0'; gui.solvemodel.timeminstep='1.0E-2'; gui.solvemodel.timerstep='10.0'; gui.solvemodel.atol='0.0010'; gui.solvemodel.rtol='0.01'; gui.solvemodel.tlist='0:0.02:1'; gui.solvemodel.tout='tlist'; gui.solvemodel.tsteps='free'; gui.solvemodel.odesolver='genalpha'; gui.solvemodel.timestep='0.01'; gui.solvemodel.incrdelay='off'; gui.solvemodel.incrdelaysteps='15'; gui.solvemodel.manualreassem='off'; gui.solvemodel.emassconst='on'; gui.solvemodel.massconst='on'; gui.solvemodel.loadconst='on'; gui.solvemodel.constrconst='on'; gui.solvemodel.jacobianconst='on'; gui.solvemodel.constrjacobianconst='on'; gui.solvemodel.manualstep='off'; gui.solvemodel.maxstepauto='on'; gui.solvemodel.initialstepauto='on'; gui.solvemodel.initialstep='0.0010'; gui.solvemodel.maxorder='5'; gui.solvemodel.minorder='1'; gui.solvemodel.maxstep='0.01'; gui.solvemodel.rhoinf='0.75'; gui.solvemodel.predictor='linear'; gui.solvemodel.timeusestopcond='off'; gui.solvemodel.timestopcond=''; gui.solvemodel.paramusestopcond='off'; gui.solvemodel.masssingular='maybe'; gui.solvemodel.consistent='bweuler'; gui.solvemodel.estrat='0'; gui.solvemodel.complex='off'; gui.solvemodel.neigs='6'; gui.solvemodel.shift='0'; gui.solvemodel.maxeigit='300'; gui.solvemodel.etol='0.0'; gui.solvemodel.krylovdim='0'; gui.solvemodel.eigname='lambda'; gui.solvemodel.eigref='0'; gui.solvemodel.pname=''; gui.solvemodel.plist=''; gui.solvemodel.pdistrib='off'; gui.solvemodel.porder='1'; gui.solvemodel.manualparam='off'; gui.solvemodel.pinitstep='0.0'; gui.solvemodel.pminstep='0.0'; gui.solvemodel.pmaxstep='0.0'; gui.solvemodel.autooldcomp='on'; gui.solvemodel.oldcomp=''; gui.solvemodel.outform='auto'; gui.solvemodel.symmetric='auto'; gui.solvemodel.symmhermit='auto'; gui.solvemodel.method='eliminate'; gui.solvemodel.nullfun='auto'; gui.solvemodel.blocksize='1000'; gui.solvemodel.blocksizeauto='on'; gui.solvemodel.uscale='auto'; gui.solvemodel.manscale=''; gui.solvemodel.rowscale='on'; gui.solvemodel.conjugate='off'; gui.solvemodel.complexfun='off'; gui.solvemodel.matherr='on'; gui.solvemodel.solfile='off'; gui.solvemodel.adaptgeom='currgeom'; gui.solvemodel.eefun='l2'; gui.solvemodel.eefunc=''; gui.solvemodel.maxt='10000000'; gui.solvemodel.rmethod='longest'; gui.solvemodel.resmethod='weak'; gui.solvemodel.resorderauto='on'; gui.solvemodel.resorder='0'; gui.solvemodel.l2scale='1'; gui.solvemodel.l2staborder='2'; gui.solvemodel.eigselect='1'; gui.solvemodel.tpfun='fltpft'; gui.solvemodel.ngen='2'; gui.solvemodel.tpmult='1.7'; gui.solvemodel.tpworst='0.5'; gui.solvemodel.tpfract='0.5'; gui.solvemodel.autosolver='on'; gui.solvemodel.varcomp=''; gui.solvemodel.oldvarcomp=''; gui.solvemodel.manualhessupd='off'; gui.solvemodel.manuallimitexpr='off'; gui.solvemodel.designsolver='sensitivity'; gui.solvemodel.sensmethod='adjoint'; gui.solvemodel.sensfunc=''; gui.solvemodel.sensfuncauto='on'; gui.solvemodel.qpsolver='cholesky'; gui.solvemodel.gradient='analytic'; gui.solvemodel.limitexpr=''; gui.solvemodel.nsolvemax='500'; gui.solvemodel.hessupd='10'; gui.solvemodel.opttol='1.0e-6'; gui.solvemodel.feastol='1.0e-6'; gui.solvemodel.majfeastol='1.0e-6'; gui.solvemodel.funcprec='1.0e-6'; gui.solvemodel.callblevel=''; gui.solvemodel.callblevelshow=''; gui.solvemodel.callbfreq=''; gui.solvemodel.callbackrough='0'; gui.solvemodel.callbclose='off'; gui.solvemodel.solcomp='T,v,u,lm,p'; gui.solvemodel.outcomp='v,T,u,lm,p'; gui.solvemodel.reacf='on'; gui.solvemodel.inittype='init_expr_currsol_radio'; gui.solvemodel.initsolnum='Automatic'; gui.solvemodel.inittime='0'; gui.solvemodel.utype='u_init_radio'; gui.solvemodel.usolnum='Automatic'; gui.solvemodel.utime='0'; gui.solvemodel.scriptcommands=''; gui.solvemodel.usescript='off'; gui.solvemodel.autoscript='off'; gui.solvemodel.sameaxis='off'; gui.solvemodel.linsolvernode.currlinsolver='umfpack'; gui.solvemodel.linsolvernode.type='linsolver'; gui.solvemodel.linsolvernode.droptol='0.0'; gui.solvemodel.linsolvernode.thresh='0.1'; gui.solvemodel.linsolvernode.umfalloc='0.7'; gui.solvemodel.linsolvernode.preorder='nd'; gui.solvemodel.linsolvernode.preroworder='on'; gui.solvemodel.linsolvernode.pivotstrategy='off'; gui.solvemodel.linsolvernode.pardreorder='nd'; gui.solvemodel.linsolvernode.pardrreorder='on'; gui.solvemodel.linsolvernode.pivotperturb='1.0E-8'; gui.solvemodel.linsolvernode.errorchk='on'; gui.solvemodel.linsolvernode.errorchkd='off'; gui.solvemodel.linsolvernode.termination='tol'; gui.solvemodel.linsolvernode.iter='2'; gui.solvemodel.linsolvernode.itol='1.0E-6'; gui.solvemodel.linsolvernode.rhob='400.0'; gui.solvemodel.linsolvernode.maxlinit='10000'; gui.solvemodel.linsolvernode.prefuntype='left'; gui.solvemodel.linsolvernode.prefuntype2='right'; gui.solvemodel.linsolvernode.iluiter='1'; gui.solvemodel.linsolvernode.itrestart='50'; gui.solvemodel.linsolvernode.seconditer='1'; gui.solvemodel.linsolvernode.relax='1.0'; gui.solvemodel.linsolvernode.amgauto='3'; gui.solvemodel.linsolvernode.mglevels='6'; gui.solvemodel.linsolvernode.mgcycle='v'; gui.solvemodel.linsolvernode.maxcoarsedof='5000'; gui.solvemodel.linsolvernode.oocmemory='512.0'; gui.solvemodel.linsolvernode.oocfilename=''; gui.solvemodel.linsolvernode.modified='off'; gui.solvemodel.linsolvernode.fillratio='2.0'; gui.solvemodel.linsolvernode.respectpattern='on'; gui.solvemodel.linsolvernode.droptype='droptol'; gui.solvemodel.linsolvernode.vankavars=''; gui.solvemodel.linsolvernode.vankasolv='gmres'; gui.solvemodel.linsolvernode.vankatol='0.02'; gui.solvemodel.linsolvernode.vankarestart='100'; gui.solvemodel.linsolvernode.vankarelax='0.8'; gui.solvemodel.linsolvernode.vankablocked='on'; gui.solvemodel.linsolvernode.sorblocked='on'; gui.solvemodel.linsolvernode.sorvecdof=''; gui.solvemodel.linsolvernode.mgauto='shape'; gui.solvemodel.linsolvernode.rmethod='regular'; gui.solvemodel.linsolvernode.coarseassem='on'; gui.solvemodel.linsolvernode.meshscale='2'; gui.solvemodel.linsolvernode.mgautolevels='2'; gui.solvemodel.linsolvernode.mgkeep='off'; gui.solvemodel.linsolvernode.mggeom='Geom1'; gui.solvemodel.linsolvernode.mcase0='on'; gui.solvemodel.linsolvernode.mgassem0='on'; gui.solvemodel.solversegmodel.seggrps{1}.segcomp='T'; gui.solvemodel.solversegmodel.seggrps{1}.ntol='1e-3'; gui.solvemodel.solversegmodel.seggrps{1}.timentol='1'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.currlinsolver='umfpack'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.type='linsolver'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.droptol='0.0'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.thresh='0.1'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.umfalloc='0.7'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.preorder='nd'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.preroworder='on'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.pivotstrategy='off'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.pardreorder='nd'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.pardrreorder='on'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.pivotperturb='1.0E-8'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.errorchk='on'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.errorchkd='off'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.termination='tol'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.iter='2'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.itol='1.0E-6'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.rhob='400.0'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.maxlinit='10000'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.prefuntype='left'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.prefuntype2='right'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.iluiter='1'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.itrestart='50'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.seconditer='1'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.relax='1.0'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.amgauto='3'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.mglevels='6'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.mgcycle='v'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.maxcoarsedof='5000'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.oocmemory='512.0'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.oocfilename=''; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.modified='off'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.fillratio='2.0'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.respectpattern='on'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.droptype='droptol'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.vankavars=''; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.vankasolv='gmres'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.vankatol='0.02'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.vankarestart='100'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.vankarelax='0.8'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.vankablocked='on'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.sorblocked='on'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.sorvecdof=''; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.mgauto='shape'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.rmethod='regular'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.coarseassem='on'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.meshscale='2'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.mgautolevels='2'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.mgkeep='off'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.mggeom='Geom1'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.mcase0='on'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.mgassem0='on'; gui.solvemodel.solversegmodel.seggrps{2}.segcomp='u v p'; gui.solvemodel.solversegmodel.seggrps{2}.ntol='1e-3'; gui.solvemodel.solversegmodel.seggrps{2}.timentol='1'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.currlinsolver='umfpack'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.type='linsolver'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.droptol='0.0'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.thresh='0.1'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.umfalloc='0.7'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.preorder='nd'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.preroworder='on'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.pivotstrategy='off'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.pardreorder='nd'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.pardrreorder='on'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.pivotperturb='1.0E-8'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.errorchk='on'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.errorchkd='off'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.termination='tol'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.iter='2'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.itol='1.0E-6'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.rhob='400.0'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.maxlinit='10000'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.prefuntype='left'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.prefuntype2='right'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.iluiter='1'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.itrestart='50'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.seconditer='1'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.relax='1.0'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.amgauto='3'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.mglevels='6'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.mgcycle='v'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.maxcoarsedof='5000'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.oocmemory='512.0'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.oocfilename=''; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.modified='off'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.fillratio='2.0'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.respectpattern='on'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.droptype='droptol'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.vankavars=''; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.vankasolv='gmres'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.vankatol='0.02'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.vankarestart='100'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.vankarelax='0.8'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.vankablocked='on'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.sorblocked='on'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.sorvecdof=''; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.mgauto='shape'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.rmethod='regular'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.coarseassem='on'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.meshscale='2'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.mgautolevels='2'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.mgkeep='off'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.mggeom='Geom1'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.mcase0='on'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.mgassem0='on'; gui.solvemodel.solversegmodel.seggrps{3}.segcomp='u2'; gui.solvemodel.solversegmodel.seggrps{3}.ntol='1e-3'; gui.solvemodel.solversegmodel.seggrps{3}.timentol='1'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.currlinsolver='umfpack'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.type='linsolver'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.droptol='0.0'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.thresh='0.1'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.umfalloc='0.7'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.preorder='nd'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.preroworder='on'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.pivotstrategy='off'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.pardreorder='nd'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.pardrreorder='on'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.pivotperturb='1.0E-8'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.errorchk='on'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.errorchkd='off'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.termination='tol'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.iter='2'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.itol='1.0E-6'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.rhob='400.0'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.maxlinit='10000'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.prefuntype='left'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.prefuntype2='right'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.iluiter='1'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.itrestart='50'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.seconditer='1'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.relax='1.0'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.amgauto='3'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.mglevels='6'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.mgcycle='v'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.maxcoarsedof='5000'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.oocmemory='512.0'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.oocfilename=''; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.modified='off'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.fillratio='2.0'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.respectpattern='on'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.droptype='droptol'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.vankavars=''; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.vankasolv='gmres'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.vankatol='0.02'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.vankarestart='100'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.vankarelax='0.8'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.vankablocked='on'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.sorblocked='on'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.sorvecdof=''; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.mgauto='shape'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.rmethod='regular'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.coarseassem='on'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.meshscale='2'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.mgautolevels='2'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.mgkeep='off'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.mggeom='Geom1'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.mcase0='on'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.mgassem0='on'; gui.solvemodel.solversegmodel.seggrps{4}.segcomp='lm'; gui.solvemodel.solversegmodel.seggrps{4}.ntol='1e-3'; gui.solvemodel.solversegmodel.seggrps{4}.timentol='1'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.currlinsolver='umfpack'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.type='linsolver'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.droptol='0.0'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.thresh='0.1'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.umfalloc='0.7'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.preorder='nd'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.preroworder='on'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.pivotstrategy='off'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.pardreorder='nd'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.pardrreorder='on'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.pivotperturb='1.0E-8'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.errorchk='on'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.errorchkd='off'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.termination='tol'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.iter='2'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.itol='1.0E-6'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.rhob='400.0'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.maxlinit='10000'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.prefuntype='left'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.prefuntype2='right'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.iluiter='1'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.itrestart='50'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.seconditer='1'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.relax='1.0'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.amgauto='3'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.mglevels='6'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.mgcycle='v'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.maxcoarsedof='5000'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.oocmemory='512.0'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.oocfilename=''; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.modified='off'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.fillratio='2.0'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.respectpattern='on'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.droptype='droptol'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.vankavars=''; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.vankasolv='gmres'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.vankatol='0.02'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.vankarestart='100'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.vankarelax='0.8'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.vankablocked='on'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.sorblocked='on'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.sorvecdof=''; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.mgauto='shape'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.rmethod='regular'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.coarseassem='on'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.meshscale='2'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.mgautolevels='2'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.mgkeep='off'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.mggeom='Geom1'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.mcase0='on'; gui.solvemodel.solversegmodel.seggrps{4}.linsolvernode.mgassem0='on'; gui.solvemodel.solversegmodel.segsteps{1}.segorder='1'; gui.solvemodel.solversegmodel.segsteps{1}.subterm='iter'; gui.solvemodel.solversegmodel.segsteps{1}.subdamp='0.5'; gui.solvemodel.solversegmodel.segsteps{1}.timesubdamp='1'; gui.solvemodel.solversegmodel.segsteps{1}.subiter='1'; gui.solvemodel.solversegmodel.segsteps{1}.maxsubiter='20'; gui.solvemodel.solversegmodel.segsteps{1}.timemaxsubiter='10'; gui.solvemodel.solversegmodel.segsteps{1}.subntol='1.0E-2'; gui.solvemodel.solversegmodel.segsteps{1}.subntolfact='1'; gui.solvemodel.solversegmodel.segsteps{1}.subdtech='const'; gui.solvemodel.solversegmodel.segsteps{1}.submandamp='off'; gui.solvemodel.solversegmodel.segsteps{1}.subinitstep='1.0'; gui.solvemodel.solversegmodel.segsteps{1}.subminstep='1.0E-4'; gui.solvemodel.solversegmodel.segsteps{1}.timesubminstep='1.0E-2'; gui.solvemodel.solversegmodel.segsteps{1}.subrstep='10.0'; gui.solvemodel.solversegmodel.segsteps{1}.timesubjtech='minimal'; gui.solvemodel.solversegmodel.segsteps{1}.subjtech='onevery'; gui.solvemodel.solversegmodel.segsteps{2}.segorder='2'; gui.solvemodel.solversegmodel.segsteps{2}.subterm='iter'; gui.solvemodel.solversegmodel.segsteps{2}.subdamp='0.5'; gui.solvemodel.solversegmodel.segsteps{2}.timesubdamp='1'; gui.solvemodel.solversegmodel.segsteps{2}.subiter='1'; gui.solvemodel.solversegmodel.segsteps{2}.maxsubiter='20'; gui.solvemodel.solversegmodel.segsteps{2}.timemaxsubiter='10'; gui.solvemodel.solversegmodel.segsteps{2}.subntol='1.0E-2'; gui.solvemodel.solversegmodel.segsteps{2}.subntolfact='1'; gui.solvemodel.solversegmodel.segsteps{2}.subdtech='const'; gui.solvemodel.solversegmodel.segsteps{2}.submandamp='off'; gui.solvemodel.solversegmodel.segsteps{2}.subinitstep='1.0'; gui.solvemodel.solversegmodel.segsteps{2}.subminstep='1.0E-4'; gui.solvemodel.solversegmodel.segsteps{2}.timesubminstep='1.0E-2'; gui.solvemodel.solversegmodel.segsteps{2}.subrstep='10.0'; gui.solvemodel.solversegmodel.segsteps{2}.timesubjtech='minimal'; gui.solvemodel.solversegmodel.segsteps{2}.subjtech='onevery'; gui.solvemodel.solversegmodel.segsteps{3}.segorder='3'; gui.solvemodel.solversegmodel.segsteps{3}.subterm='iter'; gui.solvemodel.solversegmodel.segsteps{3}.subdamp='0.5'; gui.solvemodel.solversegmodel.segsteps{3}.timesubdamp='1'; gui.solvemodel.solversegmodel.segsteps{3}.subiter='1'; gui.solvemodel.solversegmodel.segsteps{3}.maxsubiter='20'; gui.solvemodel.solversegmodel.segsteps{3}.timemaxsubiter='10'; gui.solvemodel.solversegmodel.segsteps{3}.subntol='1.0E-2'; gui.solvemodel.solversegmodel.segsteps{3}.subntolfact='1'; gui.solvemodel.solversegmodel.segsteps{3}.subdtech='const'; gui.solvemodel.solversegmodel.segsteps{3}.submandamp='off'; gui.solvemodel.solversegmodel.segsteps{3}.subinitstep='1.0'; gui.solvemodel.solversegmodel.segsteps{3}.subminstep='1.0E-4'; gui.solvemodel.solversegmodel.segsteps{3}.timesubminstep='1.0E-2'; gui.solvemodel.solversegmodel.segsteps{3}.subrstep='10.0'; gui.solvemodel.solversegmodel.segsteps{3}.timesubjtech='minimal'; gui.solvemodel.solversegmodel.segsteps{3}.subjtech='onevery'; gui.solvemodel.solversegmodel.segsteps{4}.segorder='4'; gui.solvemodel.solversegmodel.segsteps{4}.subterm='iter'; gui.solvemodel.solversegmodel.segsteps{4}.subdamp='0.5'; gui.solvemodel.solversegmodel.segsteps{4}.timesubdamp='1'; gui.solvemodel.solversegmodel.segsteps{4}.subiter='1'; gui.solvemodel.solversegmodel.segsteps{4}.maxsubiter='20'; gui.solvemodel.solversegmodel.segsteps{4}.timemaxsubiter='10'; gui.solvemodel.solversegmodel.segsteps{4}.subntol='1.0E-2'; gui.solvemodel.solversegmodel.segsteps{4}.subntolfact='1'; gui.solvemodel.solversegmodel.segsteps{4}.subdtech='const'; gui.solvemodel.solversegmodel.segsteps{4}.submandamp='off'; gui.solvemodel.solversegmodel.segsteps{4}.subinitstep='1.0'; gui.solvemodel.solversegmodel.segsteps{4}.subminstep='1.0E-4'; gui.solvemodel.solversegmodel.segsteps{4}.timesubminstep='1.0E-2'; gui.solvemodel.solversegmodel.segsteps{4}.subrstep='10.0'; gui.solvemodel.solversegmodel.segsteps{4}.timesubjtech='minimal'; gui.solvemodel.solversegmodel.segsteps{4}.subjtech='onevery'; gui.solvemodel.solversegmodel.manualsteps='off'; gui.solvemodel.solversegmodel.llimitdof=''; gui.solvemodel.solversegmodel.llimitval=''; gui.solvemodel.paramsweep.pname='a'; gui.solvemodel.paramsweep.plist='0.0001 0.00025 0.0005 0.001 0.0025 0.005 0.01 0.025 0.05 0.1 0.15 0.2 0.35 0'; gui.solvemodel.paramsweep.pdistrib='off'; gui.solvemodel.paramsweep.savefiles='on'; gui.solvemodel.paramsweep.varnames=''; gui.solvemodel.paramsweep.logfile=''; gui.registry.general_currentmodel='Geom1'; gui.registry.general_currmeshcase='0'; gui.registry.general_savedonserver='off'; gui.registry.general_savedchanges='off'; gui.registry.general_rulingmode=''; gui.registry.general_incompletemfilehistory='off'; gui.registry.saved_license='1022031'; gui.registry.saved_version='COMSOL 3.5.0.494'; gui.registry.info_modelname=''; gui.registry.info_author=''; gui.registry.info_company=''; gui.registry.info_department=''; gui.registry.info_reference=''; gui.registry.info_url=''; gui.registry.info_saveddate='1315487791492'; gui.registry.info_creationdate='1294398651117'; gui.registry.info_modelresult=''; gui.registry.spice_netlist=''; gui.registry.spice_forceac='off'; gui.reportregistry.report_contents=''; gui.reportregistry.report_outputformat='html'; gui.reportregistry.report_filename=''; gui.reportregistry.report_autoopen='off'; gui.reportregistry.report_paperformat='a4'; gui.reportregistry.report_includedefaults='off'; gui.reportregistry.report_template='full'; gui.reportregistry.report_showemptysections='off'; gui.flmodel{1}.modelname='Geom1'; gui.flmodel{1}.currmode='subdomain'; gui.flmodel{1}.currappl='0'; gui.flmodel{1}.axis.xmin='-0.004715756839344084'; gui.flmodel{1}.axis.xmax='0.004715756839344084'; gui.flmodel{1}.axis.ymin='-0.0011173167008345225'; gui.flmodel{1}.axis.ymax='0.004117316726911555'; gui.flmodel{1}.axis.zmin='-1.0'; gui.flmodel{1}.axis.zmax='1.0'; gui.flmodel{1}.axis.xspacing='5.0E-4'; gui.flmodel{1}.axis.yspacing='5.0E-4'; gui.flmodel{1}.axis.zspacing='0.01'; gui.flmodel{1}.axis.extrax=''; gui.flmodel{1}.axis.extray=''; gui.flmodel{1}.axis.extraz=''; gui.flmodel{1}.camera.xmin='-0.4400000208988786'; gui.flmodel{1}.camera.xmax='0.4400000208988786'; gui.flmodel{1}.camera.ymin='-0.26173167138730385'; gui.flmodel{1}.camera.ymax='0.26173167138730385'; gui.flmodel{1}.camera.camposx='0.0'; gui.flmodel{1}.camera.camposy='0.1500000013038516'; gui.flmodel{1}.camera.camposz='4.400000208988786'; gui.flmodel{1}.camera.camtargetx='0.0'; gui.flmodel{1}.camera.camtargety='0.1500000013038516'; gui.flmodel{1}.camera.camtargetz='0.0'; gui.flmodel{1}.camera.camupx='0.0'; gui.flmodel{1}.camera.camupy='1.0'; gui.flmodel{1}.camera.camupz='0.0'; gui.flmodel{1}.lightmodel.headlight.type='point'; gui.flmodel{1}.lightmodel.headlight.name='headlight'; gui.flmodel{1}.lightmodel.headlight.enable='on'; gui.flmodel{1}.lightmodel.headlight.colorr='255'; gui.flmodel{1}.lightmodel.headlight.colorg='255'; gui.flmodel{1}.lightmodel.headlight.colorb='255'; gui.flmodel{1}.lightmodel.headlight.xpos='0.0'; gui.flmodel{1}.lightmodel.headlight.ypos='0.15'; gui.flmodel{1}.lightmodel.headlight.zpos='4.4'; gui.flmodel{1}.lightmodel.scenelight{1}.type='directional'; gui.flmodel{1}.lightmodel.scenelight{1}.name='light 1'; gui.flmodel{1}.lightmodel.scenelight{1}.enable='on'; gui.flmodel{1}.lightmodel.scenelight{1}.colorr='255'; gui.flmodel{1}.lightmodel.scenelight{1}.colorg='255'; gui.flmodel{1}.lightmodel.scenelight{1}.colorb='255'; gui.flmodel{1}.lightmodel.scenelight{1}.xdir='1.0'; gui.flmodel{1}.lightmodel.scenelight{1}.ydir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{1}.zdir='1.0'; gui.flmodel{1}.lightmodel.scenelight{2}.type='directional'; gui.flmodel{1}.lightmodel.scenelight{2}.name='light 2'; gui.flmodel{1}.lightmodel.scenelight{2}.enable='on'; gui.flmodel{1}.lightmodel.scenelight{2}.colorr='255'; gui.flmodel{1}.lightmodel.scenelight{2}.colorg='255'; gui.flmodel{1}.lightmodel.scenelight{2}.colorb='255'; gui.flmodel{1}.lightmodel.scenelight{2}.xdir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{2}.ydir='1.0'; gui.flmodel{1}.lightmodel.scenelight{2}.zdir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{3}.type='directional'; gui.flmodel{1}.lightmodel.scenelight{3}.name='light 3'; gui.flmodel{1}.lightmodel.scenelight{3}.enable='on'; gui.flmodel{1}.lightmodel.scenelight{3}.colorr='255'; gui.flmodel{1}.lightmodel.scenelight{3}.colorg='255'; gui.flmodel{1}.lightmodel.scenelight{3}.colorb='255'; gui.flmodel{1}.lightmodel.scenelight{3}.xdir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{3}.ydir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{3}.zdir='1.0'; gui.flmodel{1}.lightmodel.scenelight{4}.type='directional'; gui.flmodel{1}.lightmodel.scenelight{4}.name='light 4'; gui.flmodel{1}.lightmodel.scenelight{4}.enable='on'; gui.flmodel{1}.lightmodel.scenelight{4}.colorr='255'; gui.flmodel{1}.lightmodel.scenelight{4}.colorg='255'; gui.flmodel{1}.lightmodel.scenelight{4}.colorb='255'; gui.flmodel{1}.lightmodel.scenelight{4}.xdir='1.0'; gui.flmodel{1}.lightmodel.scenelight{4}.ydir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{4}.zdir='-1.0'; gui.flmodel{1}.registry.axis_visible='on'; gui.flmodel{1}.registry.axis_auto='on'; gui.flmodel{1}.registry.axis_autoy='on'; gui.flmodel{1}.registry.axis_autoz='on'; gui.flmodel{1}.registry.axis_box='off'; gui.flmodel{1}.registry.axis_equal='on'; gui.flmodel{1}.registry.axis_csys='on'; gui.flmodel{1}.registry.grid_visible='on'; gui.flmodel{1}.registry.grid_auto='on'; gui.flmodel{1}.registry.grid_autoz='on'; gui.flmodel{1}.registry.grid_labels='on'; gui.flmodel{1}.registry.labels_object='on'; gui.flmodel{1}.registry.labels_vertex='off'; gui.flmodel{1}.registry.labels_edge='off'; gui.flmodel{1}.registry.labels_face='off'; gui.flmodel{1}.registry.labels_subdomain='off'; gui.flmodel{1}.registry.symbols_vertexscale='1.0'; gui.flmodel{1}.registry.symbols_edgescale='1.0'; gui.flmodel{1}.registry.symbols_facescale='1.0'; gui.flmodel{1}.registry.select_draw2d='size'; gui.flmodel{1}.registry.select_adj='cycle'; gui.flmodel{1}.registry.light_headlight='off'; gui.flmodel{1}.registry.light_scenelight='off'; gui.flmodel{1}.registry.light_shininess='0.5'; gui.flmodel{1}.registry.camera_mouse='orbit'; gui.flmodel{1}.registry.camera_camconstr='none'; gui.flmodel{1}.registry.camera_mouseconstr='none'; gui.flmodel{1}.registry.camera_perspective='off'; gui.flmodel{1}.registry.camera_moveasbox='off'; gui.flmodel{1}.registry.draw_assembly='off'; gui.flmodel{1}.registry.draw_dialog='off'; gui.flmodel{1}.registry.draw_keepborders='on'; gui.flmodel{1}.registry.draw_keepedges='off'; gui.flmodel{1}.registry.draw_multi='off'; gui.flmodel{1}.registry.draw_snap2grid='on'; gui.flmodel{1}.registry.draw_snap2vtx='on'; gui.flmodel{1}.registry.draw_solid='on'; gui.flmodel{1}.registry.draw_workplane_coordsys='on'; gui.flmodel{1}.registry.draw_workplane_showgeom='on'; gui.flmodel{1}.registry.draw_repair='on'; gui.flmodel{1}.registry.draw_repairtol='1.0E-6'; gui.flmodel{1}.registry.draw_projection='intersection'; gui.flmodel{1}.registry.transparency_value='1.0'; gui.flmodel{1}.registry.mesh_geomdetail='normal'; gui.flmodel{1}.registry.mesh_showquality='on'; gui.flmodel{1}.registry.post_cameraview='2'; gui.flmodel{1}.registry.graphics_scale='100.0'; gui.flmodel{1}.registry.render_mesh='off'; gui.flmodel{1}.registry.render_bndarrow='on'; gui.flmodel{1}.registry.render_vertex='off'; gui.flmodel{1}.registry.render_edge='on'; gui.flmodel{1}.registry.render_face='off'; gui.flmodel{1}.registry.highlight_vertex='off'; gui.flmodel{1}.registry.highlight_edge='on'; gui.flmodel{1}.registry.highlight_face='on'; gui.flmodel{1}.meshparam.hauto='5'; gui.flmodel{1}.meshparam.usehauto='on'; gui.flmodel{1}.meshparam.hmax=''; gui.flmodel{1}.meshparam.hmaxfact='1'; gui.flmodel{1}.meshparam.hcurve='0.3'; gui.flmodel{1}.meshparam.hgrad='1.3'; gui.flmodel{1}.meshparam.hcutoff='0.001'; gui.flmodel{1}.meshparam.hnarrow='1'; gui.flmodel{1}.meshparam.hpnt='10'; gui.flmodel{1}.meshparam.xscale='1.0'; gui.flmodel{1}.meshparam.yscale='1.0'; gui.flmodel{1}.meshparam.jiggle='on'; gui.flmodel{1}.meshparam.mcase='0'; gui.flmodel{1}.meshparam.rmethod='regular'; gui.flmodel{1}.meshparam.hmaxvtx={'','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hgradvtx={'','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hmaxedg={'','','','','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hcutoffedg={'','','','','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hcurveedg={'','','','','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hgradedg={'','','','','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hgradsub={'1.2','1.2','1.15','1.1'}; gui.flmodel{1}.meshparam.methodsub={'triaf','triaf','triaf','triaf'}; gui.flmodel{1}.meshparam.hmaxsub={'2e-4','6e-5','5e-5','4e-5'}; gui.flmodel{1}.mapmeshparam.edgelem={9,'100',10,'100',11,'100',13,'100',14,'100'}; gui.flmodel{1}.mapmeshparam.nedges='27'; gui.flmodel{1}.mapmeshparam.nelemtype=[2,2,2,2,2]; gui.flmodel{1}.mapmeshparam.symmetric=[0,0,0,0,0]; gui.flmodel{1}.mapmeshparam.reverse=[0,0,0,0,0]; gui.flmodel{1}.mapmeshparam.ratio={'1','1','1','1','1'}; gui.flmodel{1}.mapmeshparam.hauto='5'; gui.flmodel{1}.postmodel.postplot.triplot='on'; gui.flmodel{1}.postmodel.postplot.tridata={'T'}; gui.flmodel{1}.postmodel.postplot.trirangeauto='off'; gui.flmodel{1}.postmodel.postplot.trirangemin='500'; gui.flmodel{1}.postmodel.postplot.trirangemax='2500'; gui.flmodel{1}.postmodel.postplot.tricont='on'; gui.flmodel{1}.postmodel.postplot.trirecover='off'; gui.flmodel{1}.postmodel.postplot.triunit='K'; gui.flmodel{1}.postmodel.postplot.triheightdata={'T'}; gui.flmodel{1}.postmodel.postplot.triheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.triheightunit='K'; gui.flmodel{1}.postmodel.postplot.trimap='jet'; gui.flmodel{1}.postmodel.postplot.trimapreverse='off'; gui.flmodel{1}.postmodel.postplot.tribar='on'; gui.flmodel{1}.postmodel.postplot.triusemap='on'; gui.flmodel{1}.postmodel.postplot.tricolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.tricoloring='interp'; gui.flmodel{1}.postmodel.postplot.trifill='fill'; gui.flmodel{1}.postmodel.postplot.contplot='on'; gui.flmodel{1}.postmodel.postplot.contdata={'T'}; gui.flmodel{1}.postmodel.postplot.contcont='on'; gui.flmodel{1}.postmodel.postplot.contrecover='off'; gui.flmodel{1}.postmodel.postplot.contunit='K'; gui.flmodel{1}.postmodel.postplot.contheightdata={'T'}; gui.flmodel{1}.postmodel.postplot.contheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.contheightunit='K'; gui.flmodel{1}.postmodel.postplot.contcolordata={'T'}; gui.flmodel{1}.postmodel.postplot.contcolorrangeauto='on'; gui.flmodel{1}.postmodel.postplot.contcolorrangemin=''; gui.flmodel{1}.postmodel.postplot.contcolorrangemax=''; gui.flmodel{1}.postmodel.postplot.contcolordatacheck='off'; gui.flmodel{1}.postmodel.postplot.contcolorunit='K'; gui.flmodel{1}.postmodel.postplot.contmap='jet'; gui.flmodel{1}.postmodel.postplot.contmapreverse='off'; gui.flmodel{1}.postmodel.postplot.contbar='on'; gui.flmodel{1}.postmodel.postplot.contusemap='off'; gui.flmodel{1}.postmodel.postplot.contcolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.contlevels='20'; gui.flmodel{1}.postmodel.postplot.contvectorlevels='1806'; gui.flmodel{1}.postmodel.postplot.contisvector='on'; gui.flmodel{1}.postmodel.postplot.contlabel='off'; gui.flmodel{1}.postmodel.postplot.contfill='off'; gui.flmodel{1}.postmodel.postplot.linplot='off'; gui.flmodel{1}.postmodel.postplot.lindata={'T'}; gui.flmodel{1}.postmodel.postplot.linrangeauto='on'; gui.flmodel{1}.postmodel.postplot.linrangemin=''; gui.flmodel{1}.postmodel.postplot.linrangemax=''; gui.flmodel{1}.postmodel.postplot.lincont='on'; gui.flmodel{1}.postmodel.postplot.linrecover='off'; gui.flmodel{1}.postmodel.postplot.linunit='K'; gui.flmodel{1}.postmodel.postplot.linheightdata={'T'}; gui.flmodel{1}.postmodel.postplot.linheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.linheightunit='K'; gui.flmodel{1}.postmodel.postplot.linmap='jet'; gui.flmodel{1}.postmodel.postplot.linmapreverse='off'; gui.flmodel{1}.postmodel.postplot.linbar='on'; gui.flmodel{1}.postmodel.postplot.linusemap='on'; gui.flmodel{1}.postmodel.postplot.lincolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.lincoloring='interp'; gui.flmodel{1}.postmodel.postplot.arrowplot='on'; gui.flmodel{1}.postmodel.postplot.arrowploton='sub'; gui.flmodel{1}.postmodel.postplot.arrowdata={'u','v'}; gui.flmodel{1}.postmodel.postplot.arrowrecover='off'; gui.flmodel{1}.postmodel.postplot.arrowunit='m/s'; gui.flmodel{1}.postmodel.postplot.arrowbnddata={'',''}; gui.flmodel{1}.postmodel.postplot.arrowbndrecover='off'; gui.flmodel{1}.postmodel.postplot.arrowheightdata={'T'}; gui.flmodel{1}.postmodel.postplot.arrowheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.arrowheightunit='K'; gui.flmodel{1}.postmodel.postplot.arrowxspacing='100'; gui.flmodel{1}.postmodel.postplot.arrowxvectorspacing=''; gui.flmodel{1}.postmodel.postplot.arrowxisvector='off'; gui.flmodel{1}.postmodel.postplot.arrowyspacing='100'; gui.flmodel{1}.postmodel.postplot.arrowyvectorspacing=''; gui.flmodel{1}.postmodel.postplot.arrowyisvector='off'; gui.flmodel{1}.postmodel.postplot.arrowtype='arrow'; gui.flmodel{1}.postmodel.postplot.arrowlength='proportional'; gui.flmodel{1}.postmodel.postplot.arrowcolor='0,0,0'; gui.flmodel{1}.postmodel.postplot.arrowautoscale='on'; gui.flmodel{1}.postmodel.postplot.arrowscale='1'; gui.flmodel{1}.postmodel.postplot.princplot='off'; gui.flmodel{1}.postmodel.postplot.princdata={'','','','','','','','','','','',''}; gui.flmodel{1}.postmodel.postplot.princrecover='off'; gui.flmodel{1}.postmodel.postplot.princheightdata={'T'}; gui.flmodel{1}.postmodel.postplot.princheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.princheightunit='K'; gui.flmodel{1}.postmodel.postplot.princxspacing='8'; gui.flmodel{1}.postmodel.postplot.princxvectorspacing=''; gui.flmodel{1}.postmodel.postplot.princxisvector='off'; gui.flmodel{1}.postmodel.postplot.princyspacing='8'; gui.flmodel{1}.postmodel.postplot.princyvectorspacing=''; gui.flmodel{1}.postmodel.postplot.princyisvector='off'; gui.flmodel{1}.postmodel.postplot.princtype='arrow'; gui.flmodel{1}.postmodel.postplot.princlength='proportional'; gui.flmodel{1}.postmodel.postplot.princcolor='0,153,0'; gui.flmodel{1}.postmodel.postplot.princautoscale='on'; gui.flmodel{1}.postmodel.postplot.princscale='1'; gui.flmodel{1}.postmodel.postplot.flowplot='off'; gui.flmodel{1}.postmodel.postplot.flowdata={'tflux_T_x_cc','tflux_T_y_cc'}; gui.flmodel{1}.postmodel.postplot.flowunit='W/m^2'; gui.flmodel{1}.postmodel.postplot.flowuseexpression='off'; gui.flmodel{1}.postmodel.postplot.flowcolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.flowcolordata={'T'}; gui.flmodel{1}.postmodel.postplot.flowcolorunit='K'; gui.flmodel{1}.postmodel.postplot.flowmap='jet'; gui.flmodel{1}.postmodel.postplot.flowmapreverse='off'; gui.flmodel{1}.postmodel.postplot.flowbar='on'; gui.flmodel{1}.postmodel.postplot.flowheightdata={'T'}; gui.flmodel{1}.postmodel.postplot.flowheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.flowheightunit='K'; gui.flmodel{1}.postmodel.postplot.flowlines='20'; gui.flmodel{1}.postmodel.postplot.flowstart='sub'; gui.flmodel{1}.postmodel.postplot.flowstartx=''; gui.flmodel{1}.postmodel.postplot.flowstarty=''; gui.flmodel{1}.postmodel.postplot.flowisstartvector='off'; gui.flmodel{1}.postmodel.postplot.flowtol='0.001'; gui.flmodel{1}.postmodel.postplot.flowstattol='0.01'; gui.flmodel{1}.postmodel.postplot.flowlooptol='0.01'; gui.flmodel{1}.postmodel.postplot.flowmaxtime='Inf'; gui.flmodel{1}.postmodel.postplot.flowmaxsteps='5000'; gui.flmodel{1}.postmodel.postplot.flowback='on'; gui.flmodel{1}.postmodel.postplot.flownormal='off'; gui.flmodel{1}.postmodel.postplot.flowdistuniform='0.05'; gui.flmodel{1}.postmodel.postplot.flowlinesvel='20'; gui.flmodel{1}.postmodel.postplot.flowseedmanual='off'; gui.flmodel{1}.postmodel.postplot.flowseed1=''; gui.flmodel{1}.postmodel.postplot.flowseed2=''; gui.flmodel{1}.postmodel.postplot.flowinitref='1'; gui.flmodel{1}.postmodel.postplot.flowignoredist='0.5'; gui.flmodel{1}.postmodel.postplot.flowsat='1.3'; gui.flmodel{1}.postmodel.postplot.flowdistend='0.5'; gui.flmodel{1}.postmodel.postplot.flowdens='none'; gui.flmodel{1}.postmodel.postplot.partplot='off'; gui.flmodel{1}.postmodel.postplot.partmasstype='mass'; gui.flmodel{1}.postmodel.postplot.partplotas='lines'; gui.flmodel{1}.postmodel.postplot.predefforces=''; gui.flmodel{1}.postmodel.postplot.partforceparam0='radius,partr,1e-4'; gui.flmodel{1}.postmodel.postplot.partmass='4*pi/3*1e-9'; gui.flmodel{1}.postmodel.postplot.partforce={'','',''}; gui.flmodel{1}.postmodel.postplot.part_massless_flowdata={'tflux_T_x_cc','tflux_T_y_cc'}; gui.flmodel{1}.postmodel.postplot.part_massless_flowunit='W/m^2'; gui.flmodel{1}.postmodel.postplot.parttstartauto='on'; gui.flmodel{1}.postmodel.postplot.parttstart=''; gui.flmodel{1}.postmodel.postplot.partvelstart={'0','0','0'}; gui.flmodel{1}.postmodel.postplot.partstartptssel='partstart_explicit'; gui.flmodel{1}.postmodel.postplot.partstartdl=''; gui.flmodel{1}.postmodel.postplot.partstartedim1levels='10'; gui.flmodel{1}.postmodel.postplot.partstartedim1vectorlevels=''; gui.flmodel{1}.postmodel.postplot.partstartedim1isvector='off'; gui.flmodel{1}.postmodel.postplot.explicitcoord={'0','0'}; gui.flmodel{1}.postmodel.postplot.partuseexpression='off'; gui.flmodel{1}.postmodel.postplot.partcolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.partcolordata={'T'}; gui.flmodel{1}.postmodel.postplot.partcolorunit='K'; gui.flmodel{1}.postmodel.postplot.partmap='jet'; gui.flmodel{1}.postmodel.postplot.partmapreverse='off'; gui.flmodel{1}.postmodel.postplot.partbar='on'; gui.flmodel{1}.postmodel.postplot.partpointcolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.partpointautoscale='on'; gui.flmodel{1}.postmodel.postplot.partpointscale='1'; gui.flmodel{1}.postmodel.postplot.partdroptype='once'; gui.flmodel{1}.postmodel.postplot.partdroptimes=''; gui.flmodel{1}.postmodel.postplot.partdropfreq=''; gui.flmodel{1}.postmodel.postplot.partbnd='stick'; gui.flmodel{1}.postmodel.postplot.partmasslessrtol='0.001'; gui.flmodel{1}.postmodel.postplot.partmasslessatolmanual='off'; gui.flmodel{1}.postmodel.postplot.partmasslessatol={''}; gui.flmodel{1}.postmodel.postplot.partmasslessatolexpanded={''}; gui.flmodel{1}.postmodel.postplot.partmasslessstepsizemanual='off'; gui.flmodel{1}.postmodel.postplot.partmasslesstendauto='on'; gui.flmodel{1}.postmodel.postplot.partmasslessmaxstepsauto='on'; gui.flmodel{1}.postmodel.postplot.partmasslessedgetol='0.001'; gui.flmodel{1}.postmodel.postplot.partmasslesstvar='partt'; gui.flmodel{1}.postmodel.postplot.partmasslessstatic='off'; gui.flmodel{1}.postmodel.postplot.partmasslessres='5'; gui.flmodel{1}.postmodel.postplot.partrtol='0.001'; gui.flmodel{1}.postmodel.postplot.partatolmanual='off'; gui.flmodel{1}.postmodel.postplot.partatol={'',''}; gui.flmodel{1}.postmodel.postplot.partatolexpanded={'',''}; gui.flmodel{1}.postmodel.postplot.partstepsizemanual='off'; gui.flmodel{1}.postmodel.postplot.parttendauto='on'; gui.flmodel{1}.postmodel.postplot.partmaxstepsauto='on'; gui.flmodel{1}.postmodel.postplot.partedgetol='0.001'; gui.flmodel{1}.postmodel.postplot.partvelvar={'partu','partv','partw'}; gui.flmodel{1}.postmodel.postplot.parttvar='partt'; gui.flmodel{1}.postmodel.postplot.partstatic='off'; gui.flmodel{1}.postmodel.postplot.partres='5'; gui.flmodel{1}.postmodel.postplot.maxminplot='off'; gui.flmodel{1}.postmodel.postplot.maxminsubdata={'T'}; gui.flmodel{1}.postmodel.postplot.maxminsubrecover='off'; gui.flmodel{1}.postmodel.postplot.maxminsubdatacheck='on'; gui.flmodel{1}.postmodel.postplot.maxminsubunit='K'; gui.flmodel{1}.postmodel.postplot.maxminbnddata={'T'}; gui.flmodel{1}.postmodel.postplot.maxminbndrecover='off'; gui.flmodel{1}.postmodel.postplot.maxminbnddatacheck='off'; gui.flmodel{1}.postmodel.postplot.maxminbndunit='K'; gui.flmodel{1}.postmodel.postplot.geom='on'; gui.flmodel{1}.postmodel.postplot.roughplot='off'; gui.flmodel{1}.postmodel.postplot.autorefine='on'; gui.flmodel{1}.postmodel.postplot.refine='2'; gui.flmodel{1}.postmodel.postplot.geomnum={'Geom1'}; gui.flmodel{1}.postmodel.postplot.phase='0'; gui.flmodel{1}.postmodel.postplot.solnum='0'; gui.flmodel{1}.postmodel.postplot.selectvia='stored'; gui.flmodel{1}.postmodel.postplot.autotitle='on'; gui.flmodel{1}.postmodel.postplot.customtitle=''; gui.flmodel{1}.postmodel.postplot.smoothinternal='on'; gui.flmodel{1}.postmodel.postplot.useellogic='off'; gui.flmodel{1}.postmodel.postplot.ellogic=''; gui.flmodel{1}.postmodel.postplot.ellogictype='all'; gui.flmodel{1}.postmodel.postplot.complexfun='on'; gui.flmodel{1}.postmodel.postplot.matherr='off'; gui.flmodel{1}.postmodel.postplot.deformplot='off'; gui.flmodel{1}.postmodel.postplot.deformsub='on'; gui.flmodel{1}.postmodel.postplot.deformbnd='on'; gui.flmodel{1}.postmodel.postplot.deformsubdata={'tflux_T_x_cc','tflux_T_y_cc'}; gui.flmodel{1}.postmodel.postplot.deformsubunit='W/m^2'; gui.flmodel{1}.postmodel.postplot.deformbnddata={'',''}; gui.flmodel{1}.postmodel.postplot.deformautoscale='on'; gui.flmodel{1}.postmodel.postplot.deformscale='1'; gui.flmodel{1}.postmodel.postplot.animate_solnum='0'; gui.flmodel{1}.postmodel.postplot.animate_selectvia='stored'; gui.flmodel{1}.postmodel.postplot.outputtype='moviefile'; gui.flmodel{1}.postmodel.postplot.filetype='AVI'; gui.flmodel{1}.postmodel.postplot.width='640'; gui.flmodel{1}.postmodel.postplot.height='480'; gui.flmodel{1}.postmodel.postplot.fps='10'; gui.flmodel{1}.postmodel.postplot.statfunctype='full'; gui.flmodel{1}.postmodel.postplot.statnframes='11'; gui.flmodel{1}.postmodel.postplot.reverse='off'; gui.flmodel{1}.postmodel.postplot.movieinmatlab='off'; gui.flmodel{1}.postmodel.postplot.copyaxis='off'; gui.flmodel{1}.postmodel.intdata{1}.intdata={'T'}; gui.flmodel{1}.postmodel.intdata{1}.intrecover='off'; gui.flmodel{1}.postmodel.intdata{1}.intunit='K'; gui.flmodel{1}.postmodel.intdata{1}.phase='0'; gui.flmodel{1}.postmodel.intdata{1}.solnum='0'; gui.flmodel{1}.postmodel.intdata{1}.selectvia='stored'; gui.flmodel{1}.postmodel.intdata{2}.multiplyexpr='off'; gui.flmodel{1}.postmodel.intdata{2}.method='auto'; gui.flmodel{1}.postmodel.intdata{2}.order='4'; gui.flmodel{1}.postmodel.intdata{2}.autoorder='on'; gui.flmodel{1}.postmodel.intdata{2}.intdata={'T'}; gui.flmodel{1}.postmodel.intdata{2}.intrecover='off'; gui.flmodel{1}.postmodel.intdata{2}.intunit='m*K'; gui.flmodel{1}.postmodel.intdata{2}.phase='0'; gui.flmodel{1}.postmodel.intdata{2}.solnum='0'; gui.flmodel{1}.postmodel.intdata{2}.selectvia='stored'; gui.flmodel{1}.postmodel.intdata{3}.multiplyexpr='off'; gui.flmodel{1}.postmodel.intdata{3}.method='auto'; gui.flmodel{1}.postmodel.intdata{3}.order='4'; gui.flmodel{1}.postmodel.intdata{3}.autoorder='on'; gui.flmodel{1}.postmodel.intdata{3}.intdata={'T'}; gui.flmodel{1}.postmodel.intdata{3}.intrecover='off'; gui.flmodel{1}.postmodel.intdata{3}.intunit='m^2*K'; gui.flmodel{1}.postmodel.intdata{3}.phase='0'; gui.flmodel{1}.postmodel.intdata{3}.solnum='0'; gui.flmodel{1}.postmodel.intdata{3}.selectvia='stored'; gui.flmodel{1}.postmodel.domainplot.colordata={'T'}; gui.flmodel{1}.postmodel.domainplot.colorrangeauto='on'; gui.flmodel{1}.postmodel.domainplot.colorrangemin=''; gui.flmodel{1}.postmodel.domainplot.colorrangemax=''; gui.flmodel{1}.postmodel.domainplot.colorcont='on'; gui.flmodel{1}.postmodel.domainplot.colorrecover='off'; gui.flmodel{1}.postmodel.domainplot.colorunit='K'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacemap='jet'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacemapreverse='off'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacebar='on'; gui.flmodel{1}.postmodel.domainplot.surfacesurfaceusemap='on'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacecolor='255,0,0'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacecoloring='interp'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacefill='fill'; gui.flmodel{1}.postmodel.domainplot.extrusion='off'; gui.flmodel{1}.postmodel.domainplot.lineyaxisdata={'T'}; gui.flmodel{1}.postmodel.domainplot.lineyaxiscont='on'; gui.flmodel{1}.postmodel.domainplot.lineyaxisrecover='off'; gui.flmodel{1}.postmodel.domainplot.lineyaxisunit='K'; gui.flmodel{1}.postmodel.domainplot.linexaxisxaxistype='arc'; gui.flmodel{1}.postmodel.domainplot.linexaxisuseexpr='off'; gui.flmodel{1}.postmodel.domainplot.linexaxisdata={'T'}; gui.flmodel{1}.postmodel.domainplot.linexaxisunit='K'; gui.flmodel{1}.postmodel.domainplot.linelinestyle='solid'; gui.flmodel{1}.postmodel.domainplot.linelinecolor='cyclecolor'; gui.flmodel{1}.postmodel.domainplot.linelinemarker='none'; gui.flmodel{1}.postmodel.domainplot.linelegend='off'; gui.flmodel{1}.postmodel.domainplot.linelinelabels='off'; gui.flmodel{1}.postmodel.domainplot.linecolor='255,0,0'; gui.flmodel{1}.postmodel.domainplot.linesurfacemap='jet'; gui.flmodel{1}.postmodel.domainplot.linesurfacemapreverse='off'; gui.flmodel{1}.postmodel.domainplot.linesurfacebar='on'; gui.flmodel{1}.postmodel.domainplot.linesurfaceusemap='on'; gui.flmodel{1}.postmodel.domainplot.linesurfacecolor='255,0,0'; gui.flmodel{1}.postmodel.domainplot.linesurfacecoloring='interp'; gui.flmodel{1}.postmodel.domainplot.linesurfacefill='fill'; gui.flmodel{1}.postmodel.domainplot.pointyaxisdata={'T'}; gui.flmodel{1}.postmodel.domainplot.pointyaxisrecover='off'; gui.flmodel{1}.postmodel.domainplot.pointyaxisunit='K'; gui.flmodel{1}.postmodel.domainplot.pointxxaxistype=''; gui.flmodel{1}.postmodel.domainplot.pointxuseexpr='off'; gui.flmodel{1}.postmodel.domainplot.pointxdata={'T'}; gui.flmodel{1}.postmodel.domainplot.pointxunit='K'; gui.flmodel{1}.postmodel.domainplot.pointlinestyle='solid'; gui.flmodel{1}.postmodel.domainplot.pointlinecolor='cyclecolor'; gui.flmodel{1}.postmodel.domainplot.pointlinemarker='none'; gui.flmodel{1}.postmodel.domainplot.pointlegend='off'; gui.flmodel{1}.postmodel.domainplot.pointlinelabels='off'; gui.flmodel{1}.postmodel.domainplot.pointcolor='255,0,0'; gui.flmodel{1}.postmodel.domainplot.crossdispcolor='255,0,0'; gui.flmodel{1}.postmodel.domainplot.phase='0'; gui.flmodel{1}.postmodel.domainplot.solnum='0'; gui.flmodel{1}.postmodel.domainplot.selectvia='stored'; gui.flmodel{1}.postmodel.domainplot.autotitle='on'; gui.flmodel{1}.postmodel.domainplot.customtitle=''; gui.flmodel{1}.postmodel.domainplot.autolabelx='on'; gui.flmodel{1}.postmodel.domainplot.customlabelx=''; gui.flmodel{1}.postmodel.domainplot.autolabely='on'; gui.flmodel{1}.postmodel.domainplot.customlabely=''; gui.flmodel{1}.postmodel.domainplot.axistype={'lin','lin'}; gui.flmodel{1}.postmodel.domainplot.smoothinternal='on'; gui.flmodel{1}.postmodel.domainplot.autorefine='on'; gui.flmodel{1}.postmodel.domainplot.refine='1'; gui.flmodel{1}.postmodel.domainplot.plottypeind='0'; gui.flmodel{1}.postmodel.crossplot.extrusion='off'; gui.flmodel{1}.postmodel.crossplot.lineyaxisdata={'T'}; gui.flmodel{1}.postmodel.crossplot.lineyaxisrecover='off'; gui.flmodel{1}.postmodel.crossplot.lineyaxisunit='K'; gui.flmodel{1}.postmodel.crossplot.linexaxisxaxistype='arc'; gui.flmodel{1}.postmodel.crossplot.linexaxisuseexpr='off'; gui.flmodel{1}.postmodel.crossplot.linexaxisdata={'T'}; gui.flmodel{1}.postmodel.crossplot.linexaxisunit='K'; gui.flmodel{1}.postmodel.crossplot.linelinestyle='solid'; gui.flmodel{1}.postmodel.crossplot.linelinecolor='cyclecolor'; gui.flmodel{1}.postmodel.crossplot.linelinemarker='none'; gui.flmodel{1}.postmodel.crossplot.linelegend='off'; gui.flmodel{1}.postmodel.crossplot.linelinelabels='off'; gui.flmodel{1}.postmodel.crossplot.linecolor='255,0,0'; gui.flmodel{1}.postmodel.crossplot.linesurfacemap='jet'; gui.flmodel{1}.postmodel.crossplot.linesurfacemapreverse='off'; gui.flmodel{1}.postmodel.crossplot.linesurfacebar='on'; gui.flmodel{1}.postmodel.crossplot.linesurfaceusemap='on'; gui.flmodel{1}.postmodel.crossplot.linesurfacecolor='255,0,0'; gui.flmodel{1}.postmodel.crossplot.linesurfacecoloring='interp'; gui.flmodel{1}.postmodel.crossplot.linesurfacefill='fill'; gui.flmodel{1}.postmodel.crossplot.lineresolution='200'; gui.flmodel{1}.postmodel.crossplot.linecoord={'0','1','0','0'}; gui.flmodel{1}.postmodel.crossplot.linelevels='5'; gui.flmodel{1}.postmodel.crossplot.linevectorlevels=''; gui.flmodel{1}.postmodel.crossplot.lineisvector='off'; gui.flmodel{1}.postmodel.crossplot.lineactive='off'; gui.flmodel{1}.postmodel.crossplot.pointyaxisdata={'T'}; gui.flmodel{1}.postmodel.crossplot.pointyaxisrecover='off'; gui.flmodel{1}.postmodel.crossplot.pointyaxisunit='K'; gui.flmodel{1}.postmodel.crossplot.pointxxaxistype=''; gui.flmodel{1}.postmodel.crossplot.pointxuseexpr='off'; gui.flmodel{1}.postmodel.crossplot.pointxdata={'T'}; gui.flmodel{1}.postmodel.crossplot.pointxunit='m^2*K'; gui.flmodel{1}.postmodel.crossplot.pointlinestyle='solid'; gui.flmodel{1}.postmodel.crossplot.pointlinecolor='cyclecolor'; gui.flmodel{1}.postmodel.crossplot.pointlinemarker='none'; gui.flmodel{1}.postmodel.crossplot.pointlegend='off'; gui.flmodel{1}.postmodel.crossplot.pointlinelabels='off'; gui.flmodel{1}.postmodel.crossplot.pointcolor='255,0,0'; gui.flmodel{1}.postmodel.crossplot.pointcoord={'0','0'}; gui.flmodel{1}.postmodel.crossplot.crossdispcolor='255,0,0'; gui.flmodel{1}.postmodel.crossplot.phase='0'; gui.flmodel{1}.postmodel.crossplot.solnum='0'; gui.flmodel{1}.postmodel.crossplot.selectvia='stored'; gui.flmodel{1}.postmodel.crossplot.autotitle='on'; gui.flmodel{1}.postmodel.crossplot.customtitle=''; gui.flmodel{1}.postmodel.crossplot.autolabelx='on'; gui.flmodel{1}.postmodel.crossplot.customlabelx=''; gui.flmodel{1}.postmodel.crossplot.autolabely='on'; gui.flmodel{1}.postmodel.crossplot.customlabely=''; gui.flmodel{1}.postmodel.crossplot.axistype={'lin','lin'}; gui.flmodel{1}.postmodel.crossplot.smoothinternal='on'; gui.flmodel{1}.postmodel.crossplot.plottypeind='0'; gui.flmodel{1}.postmodel.dataexport.pntdata={'T'}; gui.flmodel{1}.postmodel.dataexport.pntrecover='off'; gui.flmodel{1}.postmodel.dataexport.pntunit='K'; gui.flmodel{1}.postmodel.dataexport.pntlocation='element'; gui.flmodel{1}.postmodel.dataexport.pntlagorder='2'; gui.flmodel{1}.postmodel.dataexport.bnddata={'T'}; gui.flmodel{1}.postmodel.dataexport.bndcont='off'; gui.flmodel{1}.postmodel.dataexport.bndrecover='off'; gui.flmodel{1}.postmodel.dataexport.bndunit='K'; gui.flmodel{1}.postmodel.dataexport.bndlocation='element'; gui.flmodel{1}.postmodel.dataexport.bndlagorder='2'; gui.flmodel{1}.postmodel.dataexport.subdata={'T'}; gui.flmodel{1}.postmodel.dataexport.subcont='off'; gui.flmodel{1}.postmodel.dataexport.subrecover='off'; gui.flmodel{1}.postmodel.dataexport.subunit='K'; gui.flmodel{1}.postmodel.dataexport.subxspacing='10'; gui.flmodel{1}.postmodel.dataexport.subxvectorspacing=''; gui.flmodel{1}.postmodel.dataexport.subxisvector='off'; gui.flmodel{1}.postmodel.dataexport.subyspacing='10'; gui.flmodel{1}.postmodel.dataexport.subyvectorspacing=''; gui.flmodel{1}.postmodel.dataexport.subyisvector='off'; gui.flmodel{1}.postmodel.dataexport.sublocation='element'; gui.flmodel{1}.postmodel.dataexport.sublagorder='2'; gui.flmodel{1}.postmodel.dataexport.phase='0'; gui.flmodel{1}.postmodel.dataexport.solnum='0'; gui.flmodel{1}.postmodel.dataexport.selectvia='stored'; gui.flmodel{1}.postmodel.dataexport.smoothinternal='on'; gui.flmodel{1}.postmodel.dataexport.exportformat='ptd'; gui.flmodel{1}.postmodel.dataexport.struct='off'; gui.flmodel{1}.postmodel.dataexport.exportedim='2'; gui.flmodel{1}.postmodel.dataexport.plotexportformat='coorddata'; gui.flmodel{1}.postmodel.dataexport.plotstruct='off'; gui.flmodel{1}.postmodel.datadisplay.fullprecision='off'; gui.flmodel{1}.postmodel.datadisplay.smoothinternal='on'; gui.flmodel{1}.postmodel.datadisplay.phase='0'; gui.flmodel{1}.postmodel.datadisplay.solnum='0'; gui.flmodel{1}.postmodel.datadisplay.selectvia='stored'; gui.flmodel{1}.postmodel.datadisplay.interpdata={'T'}; gui.flmodel{1}.postmodel.datadisplay.interprecover='off'; gui.flmodel{1}.postmodel.datadisplay.interpunit='K'; gui.flmodel{1}.postmodel.datadisplay.coord={'0','0'}; gui.flmodel{1}.postmodel.globalplot.globalyaxisexprs={}; gui.flmodel{1}.postmodel.globalplot.globalyaxisexprsdisp={}; gui.flmodel{1}.postmodel.globalplot.linexaxisxaxistype=''; gui.flmodel{1}.postmodel.globalplot.linexaxisuseexpr='off'; gui.flmodel{1}.postmodel.globalplot.globallinestyle='solid'; gui.flmodel{1}.postmodel.globalplot.globallinecolor='cyclecolor'; gui.flmodel{1}.postmodel.globalplot.globallinemarker='none'; gui.flmodel{1}.postmodel.globalplot.globallegend='off'; gui.flmodel{1}.postmodel.globalplot.globallinelabels='off'; gui.flmodel{1}.postmodel.globalplot.globalcolor='255,0,0'; gui.flmodel{1}.postmodel.globalplot.autotitle='on'; gui.flmodel{1}.postmodel.globalplot.customtitle=''; gui.flmodel{1}.postmodel.globalplot.autolabelx='on'; gui.flmodel{1}.postmodel.globalplot.customlabelx=''; gui.flmodel{1}.postmodel.globalplot.autolabely='on'; gui.flmodel{1}.postmodel.globalplot.customlabely=''; gui.flmodel{1}.postmodel.globalplot.axistype={'lin','lin'}; gui.flmodel{1}.postmodel.globalplot.phase='0'; gui.flmodt}el{1}.postmodel.globalplot.solnum='0'; gui.flmodel{1}.postmodel.globalplot.selectvia='stored'; gui.flmodel{1}.geommodel.suppressed{1}=[]; gui.flmodel{1}.geommodel.suppressed{2}=[]; gui.flmodel{1}.geommodel.suppressed{3}=[]; gui.flmodel{1}.workplane.type='0'; gui.flmodel{1}.workplane.wrkpln='0,1,0,0,0,1,0,0,0'; gui.flmodel{1}.workplane.localsys='0,1,0,0,0,0,1,0,0,0,0,1'; gui.flmodel{1}.workplane.model2d='Geom1'; gui.flmodel{1}.workplane.quicktype='10'; gui.flmodel{1}.workplane.parameter='0'; gui.flmodel{1}.workplane.zdir='up'; gui.flmodel{1}.meshmodel.meshplot.subplot='on'; gui.flmodel{1}.meshmodel.meshplot.meshsubuseexpression='on'; gui.flmodel{1}.meshmodel.meshplot.meshsubcolor='128,128,128'; gui.flmodel{1}.meshmodel.meshplot.meshsubbordercheck='on'; gui.flmodel{1}.meshmodel.meshplot.meshwiresubplot='on'; gui.flmodel{1}.meshmodel.meshplot.meshwiresubcolor='192,192,192'; gui.flmodel{1}.meshmodel.meshplot.bndplot='on'; gui.flmodel{1}.meshmodel.meshplot.meshbndcolor='0,0,0'; gui.flmodel{1}.meshmodel.meshplot.useellogic='off'; gui.flmodel{1}.meshmodel.meshplot.ellogic=''; gui.flmodel{1}.meshmodel.meshplot.ellogictype='all'; gui.flmodel{1}.meshmodel.meshplot.meshkeepfraction='1'; gui.flmodel{1}.meshmodel.meshplot.meshkeeptype='random'; gui.flmodel{1}.meshmodel.meshplot.meshmap='jet'; gui.flmodel{1}.meshmodel.meshplot.meshmapreverse='off'; gui.flmodel{1}.meshmodel.meshplot.meshbar='on'; 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? ����Mbp� ���� �����Mbp��~j�t�h?���� ��~j�t�h� ���� ��~j�t�h��~j�t�h?���� ��J�4a��~j�t�h?���� ��~j�t�X��~j�t�h?���� �Ǻ���V��~j�t�h?���� �����MbP��~j�t�h?���� ��~j�t�H��~j�t�h?���� �Ǻ���F��~j�t�h?���� �����Mb@��~j�t�h?���� � ��H�}M?���� � -C ��6Z?���� � HP�s�b?���� �����Mb@?�~j�t�h?���� �Ǻ���F?�~j�t�h?���� �����MbP?�~j�t�h?���� �Ǻ���V?�~j�t�h?���� ��~j�t�X?�~j�t�h?���� ��J�4a?�~j�t�h?���� ��~j�t�h? ���� ��~j�t�h?�~j�t�h?���� �����Mbp? ���� �����Mbp?�~j�t�h?���� � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � B e z i e r C u r v e ����Mbp� �?����Mbp��~j�t�h? �? B e z i e r C u r v e ����Mbp� �?�~j�t�h� �? B e z i e r C u r v e ����Mbp��~j�t�h? �?�~j�t�h��~j�t�h? �? B e z i e r C u r v e �~j�t�h� �?�~j�t�h? �? B e z i e r C u r v e �~j�t�h��~j�t�h? �?�J�4a��~j�t�h? �? B e z i e r C u r v e �J�4a��~j�t�h? �?�~j�t�X��~j�t�h? �? B e z i e r C u r v e �~j�t�X��~j�t�h? �?Ⱥ���V��~j�t�h? �? B e z i e r C u r v e Ⱥ���V��~j�t�h? �?����MbP��~j�t�h? �? B e z i e r C u r v e ����MbP��~j�t�h? �?�~j�t�H��~j�t�h? �? B e z i e r C u r v e �~j�t�H��~j�t�h? �?Ǻ���F��~j�t�h? �? B e z i e r C u r v e Ǻ���F��~j�t�h? �?����Mb@��~j�t�h? �? B e z i e r C u r v e ����Mb@��~j�t�h? �?����Mb@?�~j�t�h? �? B e z i e r C u r v e ����Mb@?�~j�t�h? �?Ǻ���F?�~j�t�h? �? B e z i e r C u r v e Ǻ���F?�~j�t�h? �?����MbP?�~j�t�h? �? B e z i e r C u r v e ����MbP?�~j�t�h? �?Ǻ���V?�~j�t�h? �? B e z i e r C u r v e Ǻ���V?�~j�t�h? �?�~j�t�X?�~j�t�h? �? B e z i e r C u r v e �~j�t�X?�~j�t�h? �?�J�4a?�~j�t�h? �? B e z i e r C u r v e �J�4a?�~j�t�h? �?�~j�t�h?�~j�t�h? �? B e z i e r C u r v e �~j�t�h? �?����Mbp? �? B e z i e r C u r v e �~j�t�h?�~j�t�h? �?����Mbp?�~j�t�h? �? B e z i e r C u r v e ����Mbp? �?����Mbp?�~j�t�h? �? B e z i e r C u r v e �J�4a��~j�t�h? �?�J�4a���H�}M?�;f���? ��H�}M? �? B e z i e r C u r v e Ǻ���V��~j�t�h? �?Ǻ���V�-C ��6Z?�;f���? -C ��6Z? �? B e z i e r C u r v e Ǻ���F��~j�t�h? �?Ǻ���F�HP�s�b?�;f���? HP�s�b? �? B e z i e r C u r v e ��H�}M? �?�J�4a?��H�}M?�;f���?�J�4a?�~j�t�h? �? B e z i e r C u r v e -C ��6Z? �?Ǻ���V?-C ��6Z?�;f���?Ǻ���V?�~j�t�h? �? B e z i e r C u r v e HP�s�b? �?Ǻ���F?HP�s�b?�;f���?Ǻ���F?�~j�t�h? �? A s s o c A t t r i b V e c t o r I n t z xw xsq ~ #w q ~ w uq ~ ' G e o m 2 �����|�=-C ��6? ����Mbp� ���� �����Mbp��~j�t�h?���� ��~j�t�h� ���� ��~j�t�h��~j�t�h?���� ��J�4a��~j�t�h?���� ��~j�t�X��~j�t�h?���� �Ǻ���V��~j�t�h?���� �����MbP��~j�t�h?���� ��~j�t�H��~j�t�h?���� �Ǻ���F��~j�t�h?���� �����Mb@��~j�t�h?���� � ��H�}M?���� � -C ��6Z?���� � HP�s�b?���� �����Mb@?�~j�t�h?���� �Ǻ���F?�~j�t�h?���� �����MbP?�~j�t�h?���� �Ǻ���V?�~j�t�h?���� ��~j�t�X?�~j�t�h?���� ��J�4a?�~j�t�h?���� ��~j�t�h? ���� ��~j�t�h?�~j�t�h?���� �����Mbp? ���� �����Mbp?�~j�t�h?���� � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � �? � B e z i e r C u r v e ����Mbp� �?����Mbp��~j�t�h? �? B e z i e r C u r v e ����Mbp� �?�~j�t�h� �? B e z i e r C u r v e ����Mbp��~j�t�h? �?�~j�t�h��~j�t�h? �? B e z i e r C u r v e �~j�t�h� �?�~j�t�h? �? B e z i e r C u r v e �~j�t�h��~j�t�h? �?�J�4a��~j�t�h? �? B e z i e r C u r v e �J�4a��~j�t�h? �?�~j�t�X��~j�t�h? �? B e z i e r C u r v e �~j�t�X��~j�t�h? �?Ⱥ���V��~j�t�h? �? B e z i e r C u r v e Ⱥ���V��~j�t�h? �?����MbP��~j�t�h? �? B e z i e r C u r v e ����MbP��~j�t�h? �?�~j�t�H��~j�t�h? �? B e z i e r C u r v e �~j�t�H��~j�t�h? �?Ǻ���F��~j�t�h? �? B e z i e r C u r v e Ǻ���F��~j�t�h? �?����Mb@��~j�t�h? �? B e z i e r C u r v e ����Mb@��~j�t�h? �?����Mb@?�~j�t�h? �? B e z i e r C u r v e ����Mb@?�~j�t�h? �?Ǻ���F?�~j�t�h? �? B e z i e r C u r v e Ǻ���F?�~j�t�h? �?����MbP?�~j�t�h? �? B e z i e r C u r v e ����MbP?�~j�t�h? �?Ǻ���V?�~j�t�h? �? B e z i e r C u r v e Ǻ���V?�~j�t�h? �?�~j�t�X?�~j�t�h? �? B e z i e r C u r v e �~j�t�X?�~j�t�h? �?�J�4a?�~j�t�h? �? B e z i e r C u r v e �J�4a?�~j�t�h? �?�~j�t�h?�~j�t�h? �? B e z i e r C u r v e �~j�t�h? �?����Mbp? �? B e z i e r C u r v e �~j�t�h?�~j�t�h? �?����Mbp?�~j�t�h? �? B e z i e r C u r v e ����Mbp? �?����Mbp?�~j�t�h? �? B e z i e r C u r v e �J�4a��~j�t�h? �?�J�4a���H�}M?�;f���? ��H�}M? �? B e z i e r C u r v e Ǻ���V��~j�t�h? �?Ǻ���V�-C ��6Z?�;f���? -C ��6Z? �? B e z i e r C u r v e Ǻ���F��~j�t�h? �?Ǻ���F�HP�s�b?�;f���? HP�s�b? �? B e z i e r C u r v e ��H�}M? �?�J�4a?��H�}M?�;f���?�J�4a?�~j�t�h? �? B e z i e r C u r v e -C ��6Z? �?Ǻ���V?-C ��6Z?�;f���?Ǻ���V?�~j�t�h? �? B e z i e r C u r v e HP�s�b? �?Ǻ���F?HP�s�b?�;f���?Ǻ���F?�~j�t�h? �? A s s o c A t t r i b V e c t o r I n t # xw xsr com.femlab.api.client.MFileInfo��3$�$ L femNameq ~ [ historyt [Ljava/lang/String;[ mfileTagsAndTypest [[Ljava/lang/String;[ resetHistoryq ~ ,[ storedNamesq ~ ,L versionq ~ xpw sq ~ w q ~ q ~ q ~ q ~ t COMSOL 3.5t w �t $Name: $t $Date: 2008/09/19 16:09:48 $xuq ~ t4�% COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) flclear fem % COMSOL version clear vrsn vrsn.name = 'COMSOL 3.5'; vrsn.ext = ''; vrsn.major = 0; vrsn.build = 494; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2008/09/19 16:09:48 $'; fem.version = vrsn; % Geometry g1=rect2('0.006','0.003','base','corner','pos',{'-0.003','0'},'rot','0'); g2=ellip2(2.5E-4,2.5E-4,'base','corner','pos',[0,0.0030]); g3=ellip2(5.0E-4,5.0E-4,'base','corner','pos',[-5.0E-4,0.0025]); g4=geomcomp({g1,g3},'ns',{'R1','E1'},'sf','R1+R1*E1','edge','none'); g5=geomcomp({g4},'ns',{'CO1'},'sf','CO1','edge','none'); g6=ellip2(0.0010,5.0E-4,'base','corner','pos',[-0.0010,0.0020]); g7=ellip2(5.0E-4,7.5E-4,'base','corner','pos',[0,0.0020]); g8=ellip2(0.0010,7.5E-4,'base','corner','pos',[-0.0010,0.0020]); g9=ellip2(0.0010,0.0010,'base','corner','pos',[-0.0010,0.0020]); g10=geomcomp({g5,g9},'ns',{'CO2','E1'},'sf','CO2+CO2*E1','edge','none'); g11=geomcomp({g10},'ns',{'CO1'},'sf','CO1','edge','none'); g12=ellip2(5.0E-4,5.0E-4,'base','center','pos',[-0.0015,0.0015]); g13=ellip2(0.00125,0.00125,'base','corner','pos',[-0.0015,0.0015]); g14=ellip2(0.0015,0.0015,'base','corner','pos',[-0.0015,0.0015]); g15=geomcomp({g11,g14},'ns',{'CO2','E1'},'sf','CO2+CO2*E1','edge','none'); g16=geomcomp({g15},'ns',{'CO1'},'sf','CO1','edge','none'); % Analyzed geometry clear s s.objs={g16}; s.name={'CO2'}; s.tags={'g16'}; fem.draw=struct('s',s); fem.geom=geomcsg(fem); % Constants fem.const = {'P','500[W]', ... 'Acon','0.36', ... 'rf','1.5[mm]', ... 'k_const','50[W/(m*K)]', ... 'rho_const','7870[kg/m^3]', ... 'visc_const','0.0055[Pa*s]', ... 'cp_const','600[J/(kg*K)]'}; % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) % Constants fem.const = {'P','500[W]', ... 'Acon','0.36', ... 'rf','1.5[mm]', ... 'k_const','50[W/(m*K)]', ... 'rho_const','7870[kg/m^3]', ... 'visc_const','0.0055[Pa*s]', ... 'cp_const','600[J/(kg*K)]', ... 'g','9.81[m/s^2]', ... 'alpha','12e-6[1/K]', ... 'T_m','1806[K]', ... 'A','4.3e-4[N/(m*K)]', ... 'R','8.3145[kJ/(mol*K)]', ... 'delta_H','-166.2e3[kJ/(kg*mol)]', ... 'k','0.00318', ... 'gs','1.3e-8[kg*mol/m^2]', ... 'HM','-6.2e3[kJ/(kg*mol)]', ... 'a','0.05'}; % Geometry parr={point2(0,0.0030)}; g1=geomcoerce('point',parr); % Analyzed geometry clear p s p.objs={g1}; p.name={'PT1'}; p.tags={'g1'}; s.objs={g16}; s.name={'CO2'}; s.tags={'g16'}; fem.draw=struct('p',p,'s',s); fem.geom=geomcsg(fem); carr={curve2([0,0],[0.0030,0],[1,1])}; g2=geomcoerce('curve',carr); % Analyzed geometry clear p c s p.objs={g1}; p.name={'PT1'}; p.tags={'g1'}; c.objs={g2}; c.name={'B1'}; c.tags={'g2'}; s.objs={g16}; s.name={'CO2'}; s.tags={'g16'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); % Analyzed geometry clear s s.objs={g16}; s.name={'CO2'}; s.tags={'g16'}; fem.draw=struct('s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) % Constants fem.const = {'P','500[W]', ... 'Acon','0.36', ... 'rf','1.5[mm]', ... 'k_const','50[W/(m*K)]', ... 'rho_const','7870[kg/m^3]', ... 'visc_const','0.0055[Pa*s]', ... 'cp_const','600[J/(kg*K)]', ... 'g','9.81[m/s^2]', ... 'alpha','12e-6[1/K]', ... 'T_m','1806[K]', ... 'A','4.3e-4[N/(m*K)]', ... 'R','8.3145[kJ/(mol*K)]', ... 'delta_H','-166.2e3[kJ/(kg*mol)]', ... 'k','0.00318', ... 'gs','1.3e-8[kg*mol/m^2]', ... 'HM','-6.2e3[kJ/(kg*mol)]', ... 'a','0.05', ... 'gamma_m','1.943[N/m]'}; % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) % Geometry g1=rect2(0.0010,0.0030,'base','corner','pos',[-0.0040,0]); g2=rect2(0.0010,0.0030,'base','corner','pos',[0.0030,0]); g3=geomcomp({g16,g1,g2},'ns',{'CO2','R1','R2'},'sf','CO2+R1+R2','edge','none'); % Analyzed geometry clear s s.objs={g3}; s.name={'CO1'}; s.tags={'g3'}; fem.draw=struct('s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) % Constants fem.const = {'P','500[W]', ... 'Acon','0.36', ... 'rf','1.5[mm]', ... 'k_const','50[W/(m*K)]', ... 'rho_const','7870[kg/m^3]', ... 'visc_const','0.0055[Pa*s]', ... 'cp_const','600[J/(kg*K)]', ... 'g','9.81[m/s^2]', ... 'alpha','12e-6[1/K]', ... 'T_m','1806[K]', ... 'A','4.3e-4[N/(m*K)]', ... 'R','8.3145[kJ/(mol*K)]', ... 'delta_H','-166.2e3[kJ/(kg*mol)]', ... 'k','0.00318', ... 'gs','1.3e-8[kg*mol/m^2]', ... 'HM','-6.2e3[kJ/(kg*mol)]', ... 'a','0.005', ... 'gamma_m','1.943[N/m]'}; % Geometry gg=geomedit(g3); gg{19}=beziercurve2([-7.5E-4,-7.5E-4,0],[0.0030,0.00225,0.00225],[1,0.7071067811865482,1]); g1=geomedit(g3,gg); gg=geomedit(g1); gg([23])={[]}; g2=geomedit(g1,gg); gg=geomedit(g2); gg([20])={[]}; g4=geomedit(g2,gg); gg=geomedit(g4); gg([21])={[]}; g5=geomedit(g4,gg); gg=geomedit(g5); gg([19])={[]}; g6=geomedit(g5,gg); gg=geomedit(g6); gg([19])={[]}; g7=geomedit(g6,gg); gg=geomedit(g7); gg([18])={[]}; g8=geomedit(g7,gg); g9=geomdel(g8); g10=ellip2(7.0E-4,7.0E-4,'base','center','pos',[0,0.0030]); g11=circ2('0.006','base','center','pos',{'0','0.03'},'rot','0'); g12=circ2('0.0007','base','center','pos',{'0','0.003'},'rot','0'); g13=geomcomp({g9,g12},'ns',{'CO2','C1'},'sf','CO2+CO2*C1','edge','none'); g14=ellip2(0.0010,5.0E-4,'base','center','pos',[0,0.0030]); g15=circ2('0.0014','base','center','pos',{'0','0.003'},'rot','0'); g16=geomcomp({g13,g15},'ns',{'CO1','C1'},'sf','CO1+C1*CO1','edge','none'); g17=geomcomp({g16},'ns',{'CO2'},'sf','CO2','edge','none'); g18=circ2('0.0021','base','center','pos',{'0','0.003'},'rot','0'); g19=geomcomp({g17,g18},'ns',{'CO1','C1'},'sf','CO1+CO1*C1','edge','none'); g20=geomcomp({g19},'ns',{'CO2'},'sf','CO2','edge','none'); % Analyzed geometry clear s s.objs={g20}; s.name={'CO1'}; s.tags={'g20'}; fem.draw=struct('s',s); fem.geom=geomcsg(fem); % Constants fem.const = {'P','500[W]', ... 'Acon','0.36', ... 'rf','1.5[mm]', ... 'k_const','50[W/(m*K)]', ... 'rho_const','7870[kg/m^3]', ... 'visc_const','0.0055[Pa*s]', ... 'cp_const','600[J/(kg*K)]', ... 'g','9.81[m/s^2]', ... 'alpha','12e-6[1/K]', ... 'T_m','1806[K]', ... 'A','4.3e-4[N/(m*K)]', ... 'R','8.3145[kJ/(mol*K)]', ... 'delta_H','-166.2e3[kJ/(kg*mol)]', ... 'k','0.00318', ... 'gs','1.3e-8[kg*mol/m^2]', ... 'HM','-6.2e3[kJ/(kg*mol)]', ... 'a','0.15', ... 'gamma_m','1.943[N/m]'}; % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) % Constants fem.const = {'P','500[W]', ... 'Acon','0.36', ... 'rf','1.5[mm]', ... 'k_const','50[W/(m*K)]', ... 'rho_const','7870[kg/m^3]', ... 'visc_const','0.0055[Pa*s]', ... 'cp_const','600[J/(kg*K)]', ... 'g','9.81[m/s^2]', ... 'alpha','12e-6[1/K]', ... 'T_m','1806[K]', ... 'A','4.3e-4[N/(m*K)]', ... 'R','8.3145[kJ/(mol*K)]', ... 'delta_H','-166.2e3[kJ/(kg*mol)]', ... 'k','0.00318', ... 'gs','1.3e-8[kg*mol/m^2]', ... 'HM','-6.2e3[kJ/(kg*mol)]', ... 'a','0.0010', ... 'gamma_m','1.943[N/m]'}; % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) % (Default values are not included) % Application mode 1 clear appl appl.mode.class = 'FlConvCond'; appl.assignsuffix = '_cc'; clear bnd bnd.q0 = {0,0,'qw',0}; bnd.type = {'q','q0','q','cont'}; bnd.ind = [1,2,3,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,3,1,4,4,4,4,4,4]; appl.bnd = bnd; clear equ equ.C = 'mat1_Cp(T[1/K])[J/(kg*K)]'; equ.rho = 'mat1_rho(T[1/K])[kg/m^3]'; equ.k = 'mat1_k(T[1/K])[W/(m*K)]'; equ.v = 'v'; equ.u = 'u'; equ.ind = [1,1,1,1]; appl.equ = equ; fem.appl{1} = appl; % Application mode 2 clear appl appl.mode.class = 'FlNavierStokes'; appl.gporder = {4,2}; appl.cporder = {2,1}; appl.assignsuffix = '_ns'; clear prop clear weakconstr weakconstr.value = 'off'; weakconstr.dim = {'lm2','lm3','lm4'}; prop.weakconstr = weakconstr; appl.prop = prop; clear bnd bnd.type = {'walltype','walltype','int'}; bnd.walltype = {'noslip','slip','noslip'}; bnd.ind = [1,1,2,1,2,2,2,2,2,2,2,2,2,2,2,2,2,2,1,2,1,3,3,3,3,3,3]; appl.bnd = bnd; clear equ equ.eta = 'mat1_eta(T[1/K])[Pa*s]'; equ.gporder = {{1;1;2}}; equ.idon = 1; equ.F_y = 'alpha*g*rho_cc*(T-T_m)*(T>=T_m)'; equ.rho = 'mat1_rho(T[1/K])[kg/m^3]'; equ.cporder = {{1;1;2}}; equ.delid = 0.7; equ.ck = 0.3; equ.ind = [1,1,1,1]; appl.equ = equ; fem.appl{2} = appl; % Application mode 3 clear appl appl.mode.class = 'FlPDEWBoundary'; appl.dim = {'lm','lm_t'}; appl.shape = {'shlag(1,''lm'')','shlag(2,''lm'')'}; appl.gporder = {2,4}; appl.cporder = {1,2}; appl.assignsuffix = '_wb'; clear prop clear weakconstr weakconstr.value = 'off'; weakconstr.dim = {'lm5','lm6'}; prop.weakconstr = weakconstr; appl.prop = prop; clear bnd bnd.cporder = {1,1,2}; bnd.shape = {1,1,2}; bnd.gporder = {1,1,2}; bnd.usage = {0,1,0}; bnd.weak = {0,'lm_test*(eta_ns*uy-d_gamma2*Tx)+u_test*lm',0}; bnd.ind = [1,1,2,1,2,2,2,2,2,2,2,2,2,2,2,2,2,2,1,2,1,3,3,3,3,3,3]; appl.bnd = bnd; fem.appl{3} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Scalar expressions fem.expr = {'qw','2*P*Acon/(rf^2*pi)*exp(-2*x^2/rf^2)*flc2hs(t[1/s]-0.05,0.05)', ... 'd_gamma','-A-R*gs*log(1+K*a)[1/kg]-(K*a/(1+K*a))*gs*(delta_H-HM)/T', ... 'd_gamma2','diff(gamma,T)', ... 'K','k*exp(-delta_H/(R*T)[1/kg])', ... 'gamma','gamma_m-A*(T-T_m)-R*T*gs*(log(1+k*a*exp(-delta_H/(R*T)[1/kg])))[1/kg]'}; % Descriptions clear descr descr.expr= {'qw','Wärmequelle'}; fem.descr = descr; % Descriptions descr = fem.descr; descr.const= {'g','Erdbeschleunigung','a','Grenzflächenaktivität','HM','partielle molare enthalpie','rf','Fokusierradius','Acon','Absorptionskoeffizient','k','entropiekonstante','delta_H','molare adsorptionsenthalpie','cp_const','konstante spez. Wärmekapazität','P','Abgestrahle Leistung','alpha','Temperaturausdehnungskoeffizient','R','universelle Gaskonstante','k_const','konstante Wärmeleitfähigkeit','rho_const','konstante Diche','gamma_m','Oberflächenspannung bei Raumtemperatur','visc_const','konstante Viskosität','gs','Sättigungswert'}; fem.descr = descr; % Library materials clear lib lib.mat{1}.name='SR235JG'; lib.mat{1}.varname='mat1'; lib.mat{1}.variables.eta='eta(T[1/K])[Pa*s]'; lib.mat{1}.variables.C='Cp(T[1/K])[J/(kg*K)]'; lib.mat{1}.variables.rho='rho(T[1/K])[kg/m^3]'; lib.mat{1}.variables.k='k(T[1/K])[W/(m*K)]'; clear fcns fcns{1}.type='inline'; fcns{1}.name='rho(T)'; fcns{1}.expr='7870-855*flc2hs(T-1806,50)'; fcns{1}.dexpr={'diff(7870-855*flc2hs(T-1806,50),T)'}; fcns{2}.type='inline'; fcns{2}.name='eta(T)'; fcns{2}.expr='500-499.9945*flc2hs(T-1806,50)'; fcns{2}.dexpr={'diff(500-499.9945*flc2hs(T-1806,50),T)'}; fcns{3}.type='interp'; fcns{3}.name='k'; fcns{3}.method='cubic'; fcns{3}.extmethod='extrap'; fcns{3}.x={'293','373','473','573','673','773','873','5000'}; fcns{3}.data={'57','57','54','50','47','42','37','37'}; fcns{4}.type='interp'; fcns{4}.name='Cp'; fcns{4}.method='cubic'; fcns{4}.extmethod='extrap'; fcns{4}.x={'293','373','473','573','673','773','873','5000'}; fcns{4}.data={'430','470','500','540','580','660','780','780'}; lib.mat{1}.functions = fcns; lib.mat{2}.name='SR235JG_1'; lib.mat{2}.varname='mat2'; lib.mat{2}.variables.eta='eta(T[1/K])[Pa*s]'; lib.mat{2}.variables.C='Cp(T[1/K])[J/(kg*K)]'; lib.mat{2}.variables.rho='rho(T[1/K])[kg/m^3]'; lib.mat{2}.variables.k='k(T[1/K])[W/(m*K)]'; clear fcns fcns{1}.type='inline'; fcns{1}.name='rho(T)'; fcns{1}.expr='7870-855*flc2hs(T-1806,50)'; fcns{1}.dexpr={'diff(7870-855*flc2hs(T-1806,50),T)'}; fcns{2}.type='inline'; fcns{2}.name='eta(T)'; fcns{2}.expr='500-499.9945*flc2hs(T-1806,50)'; fcns{2}.dexpr={'diff(500-499.9945*flc2hs(T-1806,50),T)'}; fcns{3}.type='interp'; fcns{3}.name='k'; fcns{3}.method='cubic'; fcns{3}.extmethod='extrap'; fcns{3}.x={'293','373','473','573','673','773','873','5000'}; fcns{3}.data={'57','57','54','50','47','42','37','37'}; fcns{4}.type='interp'; fcns{4}.name='Cp'; fcns{4}.method='cubic'; fcns{4}.extmethod='extrap'; fcns{4}.x={'293','373','473','573','673','773','873','5000'}; fcns{4}.data={'430','470','500','540','580','660','780','780'}; lib.mat{2}.functions = fcns; lib.mat{3}.name='SR235JG_2'; lib.mat{3}.varname='mat3'; lib.mat{3}.variables.eta='eta(T[1/K])[Pa*s]'; lib.mat{3}.variables.C='Cp(T[1/K])[J/(kg*K)]'; lib.mat{3}.variables.rho='rho(T[1/K])[kg/m^3]'; lib.mat{3}.variables.k='k(T[1/K])[W/(m*K)]'; clear fcns fcns{1}.type='inline'; fcns{1}.name='rho(T)'; fcns{1}.expr='7870-855*flc2hs(T-1806,50)'; fcns{1}.dexpr={'diff(7870-855*flc2hs(T-1806,50),T)'}; fcns{2}.type='inline'; fcns{2}.name='eta(T)'; fcns{2}.expr='500-499.9945*flc2hs(T-1806,50)'; fcns{2}.dexpr={'diff(500-499.9945*flc2hs(T-1806,50),T)'}; fcns{3}.type='interp'; fcns{3}.name='k'; fcns{3}.method='cubic'; fcns{3}.extmethod='extrap'; fcns{3}.x={'293','373','473','573','673','773','873','5000'}; fcns{3}.data={'57','57','54','50','47','42','37','37'}; fcns{4}.type='interp'; fcns{4}.name='Cp'; fcns{4}.method='cubic'; fcns{4}.extmethod='extrap'; fcns{4}.x={'293','373','473','573','673','773','873','5000'}; fcns{4}.data={'430','470','500','540','580','660','780','780'}; lib.mat{3}.functions = fcns; fem.lib = lib; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; fem.ode=ode; % Multiphysics fem=multiphysics(fem); uq ~ t &% COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) flclear fem % COMSOL version clear vrsn vrsn.name = 'COMSOL 3.5'; vrsn.ext = ''; vrsn.major = 0; vrsn.build = 494; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2008/09/19 16:09:48 $'; fem.version = vrsn; % Geometry g1=rect2('0.006','0.003','base','corner','pos',{'-0.003','0'},'rot','0'); g2=ellip2(2.5E-4,2.5E-4,'base','corner','pos',[0,0.0030]); g3=ellip2(5.0E-4,5.0E-4,'base','corner','pos',[-5.0E-4,0.0025]); g4=geomcomp({g1,g3},'ns',{'R1','E1'},'sf','R1+R1*E1','edge','none'); g5=geomcomp({g4},'ns',{'CO1'},'sf','CO1','edge','none'); g6=ellip2(0.0010,5.0E-4,'base','corner','pos',[-0.0010,0.0020]); g7=ellip2(5.0E-4,7.5E-4,'base','corner','pos',[0,0.0020]); g8=ellip2(0.0010,7.5E-4,'base','corner','pos',[-0.0010,0.0020]); g9=ellip2(0.0010,0.0010,'base','corner','pos',[-0.0010,0.0020]); g10=geomcomp({g5,g9},'ns',{'CO2','E1'},'sf','CO2+CO2*E1','edge','none'); g11=geomcomp({g10},'ns',{'CO1'},'sf','CO1','edge','none'); g12=ellip2(5.0E-4,5.0E-4,'base','center','pos',[-0.0015,0.0015]); g13=ellip2(0.00125,0.00125,'base','corner','pos',[-0.0015,0.0015]); g14=ellip2(0.0015,0.0015,'base','corner','pos',[-0.0015,0.0015]); g15=geomcomp({g11,g14},'ns',{'CO2','E1'},'sf','CO2+CO2*E1','edge','none'); g16=geomcomp({g15},'ns',{'CO1'},'sf','CO1','edge','none'); % Analyzed geometry clear s s.objs={g16}; s.name={'CO2'}; s.tags={'g16'}; fem.draw=struct('s',s); fem.geom=geomcsg(fem); % Constants fem.const = {'P','500[W]', ... 'Acon','0.36', ... 'rf','1.5[mm]', ... 'k_const','50[W/(m*K)]', ... 'rho_const','7870[kg/m^3]', ... 'visc_const','0.0055[Pa*s]', ... 'cp_const','600[J/(kg*K)]'}; % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) % Constants fem.const = {'P','500[W]', ... 'Acon','0.36', ... 'rf','1.5[mm]', ... 'k_const','50[W/(m*K)]', ... 'rho_const','7870[kg/m^3]', ... 'visc_const','0.0055[Pa*s]', ... 'cp_const','600[J/(kg*K)]', ... 'g','9.81[m/s^2]', ... 'alpha','12e-6[1/K]', ... 'T_m','1806[K]', ... 'A','4.3e-4[N/(m*K)]', ... 'R','8.3145[kJ/(mol*K)]', ... 'delta_H','-166.2e3[kJ/(kg*mol)]', ... 'k','0.00318', ... 'gs','1.3e-8[kg*mol/m^2]', ... 'HM','-6.2e3[kJ/(kg*mol)]', ... 'a','0.05'}; % Geometry parr={point2(0,0.0030)}; g1=geomcoerce('point',parr); % Analyzed geometry clear p s p.objs={g1}; p.name={'PT1'}; p.tags={'g1'}; s.objs={g16}; s.name={'CO2'}; s.tags={'g16'}; fem.draw=struct('p',p,'s',s); fem.geom=geomcsg(fem); carr={curve2([0,0],[0.0030,0],[1,1])}; g2=geomcoerce('curve',carr); % Analyzed geometry clear p c s p.objs={g1}; p.name={'PT1'}; p.tags={'g1'}; c.objs={g2}; c.name={'B1'}; c.tags={'g2'}; s.objs={g16}; s.name={'CO2'}; s.tags={'g16'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); % Analyzed geometry clear s s.objs={g16}; s.name={'CO2'}; s.tags={'g16'}; fem.draw=struct('s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) % Constants fem.const = {'P','500[W]', ... 'Acon','0.36', ... 'rf','1.5[mm]', ... 'k_const','50[W/(m*K)]', ... 'rho_const','7870[kg/m^3]', ... 'visc_const','0.0055[Pa*s]', ... 'cp_const','600[J/(kg*K)]', ... 'g','9.81[m/s^2]', ... 'alpha','12e-6[1/K]', ... 'T_m','1806[K]', ... 'A','4.3e-4[N/(m*K)]', ... 'R','8.3145[kJ/(mol*K)]', ... 'delta_H','-166.2e3[kJ/(kg*mol)]', ... 'k','0.00318', ... 'gs','1.3e-8[kg*mol/m^2]', ... 'HM','-6.2e3[kJ/(kg*mol)]', ... 'a','0.05', ... 'gamma_m','1.943[N/m]'}; % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) % Geometry g1=rect2(0.0010,0.0030,'base','corner','pos',[-0.0040,0]); g2=rect2(0.0010,0.0030,'base','corner','pos',[0.0030,0]); g3=geomcomp({g16,g1,g2},'ns',{'CO2','R1','R2'},'sf','CO2+R1+R2','edge','none'); % Analyzed geometry clear s s.objs={g3}; s.name={'CO1'}; s.tags={'g3'}; fem.draw=struct('s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) % Constants fem.const = {'P','500[W]', ... 'Acon','0.36', ... 'rf','1.5[mm]', ... 'k_const','50[W/(m*K)]', ... 'rho_const','7870[kg/m^3]', ... 'visc_const','0.0055[Pa*s]', ... 'cp_const','600[J/(kg*K)]', ... 'g','9.81[m/s^2]', ... 'alpha','12e-6[1/K]', ... 'T_m','1806[K]', ... 'A','4.3e-4[N/(m*K)]', ... 'R','8.3145[kJ/(mol*K)]', ... 'delta_H','-166.2e3[kJ/(kg*mol)]', ... 'k','0.00318', ... 'gs','1.3e-8[kg*mol/m^2]', ... 'HM','-6.2e3[kJ/(kg*mol)]', ... 'a','0.005', ... 'gamma_m','1.943[N/m]'}; % Geometry gg=geomedit(g3); gg{19}=beziercurve2([-7.5E-4,-7.5E-4,0],[0.0030,0.00225,0.00225],[1,0.7071067811865482,1]); g1=geomedit(g3,gg); gg=geomedit(g1); gg([23])={[]}; g2=geomedit(g1,gg); gg=geomedit(g2); gg([20])={[]}; g4=geomedit(g2,gg); gg=geomedit(g4); gg([21])={[]}; g5=geomedit(g4,gg); gg=geomedit(g5); gg([19])={[]}; g6=geomedit(g5,gg); gg=geomedit(g6); gg([19])={[]}; g7=geomedit(g6,gg); gg=geomedit(g7); gg([18])={[]}; g8=geomedit(g7,gg); g9=geomdel(g8); g10=ellip2(7.0E-4,7.0E-4,'base','center','pos',[0,0.0030]); g11=circ2('0.006','base','center','pos',{'0','0.03'},'rot','0'); g12=circ2('0.0007','base','center','pos',{'0','0.003'},'rot','0'); g13=geomcomp({g9,g12},'ns',{'CO2','C1'},'sf','CO2+CO2*C1','edge','none'); g14=ellip2(0.0010,5.0E-4,'base','center','pos',[0,0.0030]); g15=circ2('0.0014','base','center','pos',{'0','0.003'},'rot','0'); g16=geomcomp({g13,g15},'ns',{'CO1','C1'},'sf','CO1+C1*CO1','edge','none'); g17=geomcomp({g16},'ns',{'CO2'},'sf','CO2','edge','none'); g18=circ2('0.0021','base','center','pos',{'0','0.003'},'rot','0'); g19=geomcomp({g17,g18},'ns',{'CO1','C1'},'sf','CO1+CO1*C1','edge','none'); g20=geomcomp({g19},'ns',{'CO2'},'sf','CO2','edge','none'); % Analyzed geometry clear s s.objs={g20}; s.name={'CO1'}; s.tags={'g20'}; fem.draw=struct('s',s); fem.geom=geomcsg(fem); % Constants fem.const = {'P','500[W]', ... 'Acon','0.36', ... 'rf','1.5[mm]', ... 'k_const','50[W/(m*K)]', ... 'rho_const','7870[kg/m^3]', ... 'visc_const','0.0055[Pa*s]', ... 'cp_const','600[J/(kg*K)]', ... 'g','9.81[m/s^2]', ... 'alpha','12e-6[1/K]', ... 'T_m','1806[K]', ... 'A','4.3e-4[N/(m*K)]', ... 'R','8.3145[kJ/(mol*K)]', ... 'delta_H','-166.2e3[kJ/(kg*mol)]', ... 'k','0.00318', ... 'gs','1.3e-8[kg*mol/m^2]', ... 'HM','-6.2e3[kJ/(kg*mol)]', ... 'a','0.15', ... 'gamma_m','1.943[N/m]'}; % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) % Constants fem.const = {'P','500[W]', ... 'Acon','0.36', ... 'rf','1.5[mm]', ... 'k_const','50[W/(m*K)]', ... 'rho_const','7870[kg/m^3]', ... 'visc_const','0.0055[Pa*s]', ... 'cp_const','600[J/(kg*K)]', ... 'g','9.81[m/s^2]', ... 'alpha','12e-6[1/K]', ... 'T_m','1806[K]', ... 'A','4.3e-4[N/(m*K)]', ... 'R','8.3145[kJ/(mol*K)]', ... 'delta_H','-166.2e3[kJ/(kg*mol)]', ... 'k','0.00318', ... 'gs','1.3e-8[kg*mol/m^2]', ... 'HM','-6.2e3[kJ/(kg*mol)]', ... 'a','0.0010', ... 'gamma_m','1.943[N/m]'}; tclear mfile clear vrsn vrsn.name = 'COMSOL 3.5'; vrsn.ext = ''; vrsn.major = 0; vrsn.build = 494; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2008/09/19 16:09:48 $'; mfile.version=vrsn; mfile.fem='fem'; mfile.stored={'fem0','fem1'}; mfile.tags={}; mfile.types={}; x