��sr!com.femlab.server.ModelFileHeader�D��%LtagstLcom/femlab/util/FlStringList;Ltypesq~LvrsntLcom/femlab/util/FlVersion;xpwsrcom.femlab.util.FlVersion��%�/B = IbuildImajorLdatetLjava/lang/String;Lextq~Lnameq~Lrcsq~LreactionExtq~L reactionNameq~L scriptExtq~L scriptNameq~xpwtCOMSOL Script 1.3tt#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��{Gxp t modelinfotxfemtguitg7tg2tg1tg4tg10tfem13t mfileinfouq~ q~t femstructt guistructtdrawq~q~q~q~tgeomq~xsrcom.femlab.api.client.ModelInfo�^��%Ldescrq~LdocURLq~[imaget[Bxpwptpxuq~ t��clear 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 = 'TIME'; 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,''lm1'')','shlag(2,''lm2'')','shlag(1,''lm3'')','shlag(2,''u'')','shlag(2,''v'')','shlag(1,''p'')'}; appl.gporder = {30,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='ale'; clear weakconstr weakconstr.value = 'on'; weakconstr.dim = {'lm1','lm2','lm3'}; prop.weakconstr = weakconstr; prop.constrtype='non-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,1,1,1,1,1,1, ... 1,1,1,1,1,1]; appl.pnt = pnt; clear bnd bnd.name = {'','','','','','',''}; bnd.type = {'inlet','walltype','int','walltype','walltype','int','outlet'}; bnd.intype = {'p','p','uv','uv','uv','p','uv'}; bnd.outtype = {'p','p','p','p','p','p','p'}; bnd.walltype = {'noslip','noslip','noslip','noslip','lwall','noslip','noslip'}; bnd.opentype = {'ntotstress','ntotstress','ntotstress','ntotstress','ntotstress', ... 'ntotstress','ntotstress'}; bnd.stresstype = {'totstress','totstress','totstress','totstress','totstress','totstress', ... 'totstress'}; bnd.inttype = {'cont','cont','cont','cont','cont','cont','cont'}; bnd.velType = {'U0in','U0in','U0in','U0in','U0in','U0in','U0in'}; bnd.u0 = {'0','0','0','0','0','0','0'}; bnd.v0 = {'0','0','0','0','0','0','0'}; bnd.p0 = {'pressure*(t[1/s]*trange1+2*trange2)','10','0','0','0','pressure*(t[1/s]*trange1+trange2)', ... '0'}; bnd.f0 = {'0','0','0','0','0','0','0'}; bnd.Fbnd = {{'0';'0'},{'0';'0'},{'0';'0'},{'0';'0'},{'0';'0'},{'0';'0'},{'0'; ... '0'}}; bnd.U0in = {'.5','.5','1','1','1','.5','1'}; bnd.U0out = {'0','0','0','0','0','0','0'}; bnd.uvw = {'0','0','0','0','0','0','0'}; bnd.uwall = {'0','0','0','0','xt','0','0'}; bnd.vwall = {'0','0','0','0','yt','0','0'}; bnd.ww = {'0','0','0','0','0','0','0'}; bnd.weakconstr = {'1','1','1','1','1','1','1'}; bnd.wcshape = {[1;2;3],[1;2;3],[1;2;3],[1;2;3],[1;2;3],[1;2;3],[1;2;3]}; bnd.wcgporder = {{1;1;1},{1;1;1},{1;1;1},{1;1;1},{1;1;1},{1;1;1},{1;1;1}}; bnd.wcinit = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'},{'0'; ... '0';'0'},{'0';'0';'0'},{'0';'0';'0'}}; bnd.style = {{{'0'},{'0','0','0'},{'solid'}},{{'0'},{'0','0','0'},{'solid'}}, ... {{'0'},{'0','255','255'},{'solid'}},{{'0'},{'0','0','255'},{'solid'}},{{'0'}, ... {'255','0','255'},{'solid'}},{{'0'},{'0','0','0'},{'solid'}},{{'0'},{'255', ... '0','255'},{'solid'}}}; bnd.ind = [1,2,1,3,4,5,5,5,5,3,3,3,5,4,6,5,6,4,3,4,3,3,3,3,4,7,7,5,5,5, ... 5,5,5,5,5,5,5,5,5,5,5,5,5]; appl.bnd = bnd; clear equ equ.shape = {[4;5;6]}; equ.gporder = {{2;2;3}}; equ.cporder = {{1;1;2}}; equ.init = {{'0';'0';'0';'0';'0'}}; equ.usage = {1}; equ.rho = {'1.2'}; equ.eta = {'1.8e-5'}; equ.F_x = {'0'}; equ.F_y = {'0'}; equ.idon = {'0'}; equ.delid = {'0.5'}; 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.1'}; 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,1,1,1,1]; appl.equ = equ; appl.var = {'visc_vel_fact','10'}; xfem.appl{1} = appl; clear appl appl.mode.class = 'MovingMesh'; appl.mode.type = 'cartesian'; appl.dim = {}; appl.sdim = {'Xm','Ym','Zm'}; appl.name = 'ale'; appl.shape = {'shlag(2,''lm4'')','shlag(2,''lm5'')','shlag(2,''x'')','shlag(2,''y'')'}; appl.gporder = {30,4}; appl.cporder = 2; appl.sshape = 2; appl.border = 'on'; appl.assignsuffix = '_ale'; clear prop prop.elemdefault='Lag2'; prop.smoothing='winslow'; prop.analysis='transient'; prop.allowremesh='on'; prop.definesframe='ale'; prop.origrefframe='ref'; prop.frame='mesh'; clear weakconstr weakconstr.value = 'on'; weakconstr.dim = {'lm4','lm5'}; prop.weakconstr = weakconstr; prop.constrtype='non-ideal'; appl.prop = prop; clear bnd bnd.name = {'','','','','','','','','','','','','','','','','','','','','', ... '','','','','','','','',''}; bnd.constrcoord = {'global','global','global','global','global','global','global', ... 'global','global','global','global','global','global','global','global', ... 'global','global','global','global','global','global','global','global', ... 'global','global','global','global','global','global','global'}; bnd.type = {'def','def','def','def','def','def','def','def','def','def','def', ... 'def','def','def','def','def','def','def','def','def','def','def','def', ... 'def','def','def','def','def','def','def'}; bnd.deform = {{'0';'0'},{'0';'-m5*(s)'},{'0';'m5*(s)'},{'m4*(1-s)';'-(m2-m3)*s-m3'}, ... {'m4*(1-s)';'(m2-m3)*s+m3'},{'0';'-m2*(1-s)'},{'0';'m2*s'},{'0';'-m1*(s)'}, ... {'0';'-m1'},{'0';'m1*(1-s)'},{'0';'m1'},{'0';'-m1*s'},{'0';'-m1*(1-s)'}, ... {'0';'m1*s'},{'m17*s';'-(m16-m5)*s-m5'},{'m17*s';'(m16-m5)*s+m5'},{'(m9-m17)*s+m17'; ... '-(m8-m16)*s-m16'},{'(m9-m17)*s+m17';'(m8-m16)*s+m16'},{'(m15-m9)*s+m9'; ... '-(m14-m8)*s-m8'},{'(m15-m9)*s+m9';'(m14-m8)*s+m8'},{'(m4-m15)*s+m15';'-(m3-m14)*s-m14'}, ... {'(m4-m15)*s+m15';'(m3-m14)*s+m14'},{'m13*s';'-(m12-m2)*s-m2'},{'m13*s'; ... '(m12-m2)*s+m2'},{'(m7-m13)*s+m13';'-(m6-m12)*s-m12'},{'(m7-m13)*s+m13'; ... '(m6-m12)*s+m12'},{'(m11-m7)*s+m7';'-(m10-m6)*s-m6'},{'(m11-m7)*s+m7';'(m10-m6)*s+m6'}, ... {'m11*(1-s)';'-(m1-m10)*s-m10'},{'m11*(1-s)';'(m1-m10)*s+m10'}}; bnd.veldeform = {{'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.defflag = {{'1';'1'},{'1';'1'},{'1';'1'},{'1';'1'},{'1';'1'},{'1';'1'},{'1'; ... '1'},{'1';'1'},{'1';'1'},{'1';'1'},{'1';'1'},{'1';'1'},{'1';'1'},{'1';'1'}, ... {'1';'1'},{'1';'1'},{'1';'1'},{'1';'1'},{'1';'1'},{'1';'1'},{'1';'1'},{'1'; ... '1'},{'1';'1'},{'1';'1'},{'1';'1'},{'1';'1'},{'1';'1'},{'1';'1'},{'1';'1'}, ... {'1';'1'}}; bnd.veldefflag = {{'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.weakconstr = {'1','1','1','1','1','1','1','1','1','1','1','1','1','1','1','1', ... '1','1','1','1','1','1','1','1','1','1','1','1','1','1'}; bnd.wcshape = {[1;2],[1;2],[1;2],[1;2],[1;2],[1;2],[1;2],[1;2],[1;2],[1;2],[1; ... 2],[1;2],[1;2],[1;2],[1;2],[1;2],[1;2],[1;2],[1;2],[1;2],[1;2],[1;2],[1; ... 2],[1;2],[1;2],[1;2],[1;2],[1;2],[1;2],[1;2]}; bnd.wcgporder = {{1;1},{1;1},{1;1},{1;1},{1;1},{1;1},{1;1},{1;1},{1;1},{1;1},{1; ... 1},{1;1},{1;1},{1;1},{1;1},{1;1},{1;1},{1;1},{1;1},{1;1},{1;1},{1;1},{1; ... 1},{1;1},{1;1},{1;1},{1;1},{1;1},{1;1},{1;1}}; bnd.wcinit = {{'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.style = {{{'0'},{'0','0','255'},{'solid'}},{{'0'},{'0','255','255'},{'solid'}}, ... {{'0'},{'0','0','255'},{'solid'}},{{'0'},{'0','255','255'},{'solid'}},{{'0'}, ... {'0','255','0'},{'solid'}},{{'0'},{'0','255','0'},{'solid'}},{{'0'},{'0', ... '0','255'},{'solid'}},{{'0'},{'0','0','255'},{'solid'}},{{'0'},{'0','0', ... '0'},{'solid'}},{{'0'},{'0','255','0'},{'solid'}},{{'0'},{'255','0','255'}, ... {'solid'}},{{'0'},{'0','0','0'},{'solid'}},{{'0'},{'255','0','255'},{'solid'}}, ... {{'0'},{'0','255','255'},{'solid'}},{{'0'},{'255','0','255'},{'solid'}}, ... {{'0'},{'0','0','255'},{'solid'}},{{'0'},{'0','0','0'},{'solid'}},{{'0'}, ... {'0','0','0'},{'solid'}},{{'0'},{'0','0','255'},{'solid'}},{{'0'},{'255', ... '0','255'},{'solid'}},{{'0'},{'0','255','0'},{'solid'}},{{'0'},{'0','255', ... '255'},{'solid'}},{{'0'},{'255','0','255'},{'solid'}},{{'0'},{'0','0','255'}, ... {'solid'}},{{'0'},{'0','0','0'},{'solid'}},{{'0'},{'0','0','0'},{'solid'}}, ... {{'0'},{'0','0','255'},{'solid'}},{{'0'},{'255','0','255'},{'solid'}},{{'0'}, ... {'0','255','0'},{'solid'}},{{'0'},{'0','255','255'},{'solid'}}}; bnd.ind = [1,1,1,1,1,2,3,4,5,6,7,1,8,1,9,10,11,1,12,1,13,14,1,10,1,1,1, ... 15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30]; appl.bnd = bnd; clear equ equ.shape = {[3;4]}; equ.gporder = {{2;2}}; equ.cporder = {{1;1}}; equ.init = {{'xinit_ale';'yinit_ale'}}; equ.usage = {1}; equ.type = {'free'}; equ.physexpr = {{'0';'0'}}; equ.presexpr = {{'0';'0'}}; equ.dinit = {{'0';'0'}}; equ.style = {{{'0'},{'193','193','193'}}}; equ.ind = [1,1,1,1,1,1,1,1]; appl.equ = equ; xfem.appl{2} = appl; xfem.geom = flbinary('fem13','geom','rampsupport.mph'); xfem.sdim = {{'Xm','Ym'},{'X','Y'},{'x','y'}}; xfem.frame = {'mesh','ref','ale'}; xfem.shape = {'shlag(''order'',2,''basename'',''u'',''frame'',''ale'')','shlag(''order'',2,''basename'',''v'',''frame'',''ale'')','shlag(''order'',1,''basename'',''p'',''frame'',''ale'')','shlag(''order'',2,''basename'',''x'',''frame'',''ale'')','shlag(''order'',2,''basename'',''y'',''frame'',''ale'')','shlag(2,''X'')','shlag(2,''Y'')','shlag(1,''lm3'')','shlag(2,''lm4'')','shlag(2,''lm5'')','shlag(2,''lm1'')','shlag(2,''lm2'')'}; xfem.gporder = {4,2,30}; xfem.cporder = {2,1}; clear sshape sshape.dvolname = 'dvol'; sshape.frame = 'mesh'; sshape.sorder = 2; sshape.type = 'fixed'; xfem.sshape{1} = sshape; clear sshape sshape.dvolname = 'dvol_ale'; sshape.frame = 'ale'; sshape.sorder = 2; sshape.type = 'moving_abs'; sshape.refframe = 'mesh'; xfem.sshape{2} = sshape; clear sshape sshape.dvolname = 'dvol_ref'; sshape.frame = 'ref'; sshape.sorder = 2; sshape.type = 'moving_abs'; sshape.refframe = 'mesh'; xfem.sshape{3} = sshape; xfem.simplify = 'on'; xfem.border = 1; xfem.form = 'weak'; clear units; units.basesystem = 'SI'; xfem.units = units; clear equ equ.shape = {[1;2;3;4;5;6;7]}; equ.gporder = {{1;1;2;1;1}}; equ.cporder = {{1;1;2;1;1}}; equ.init = {{'0';'0';'0';'xinit_ale';'yinit_ale'}}; equ.dinit = {{'0';'0';'0';'0';'0'}}; equ.weak = {{'(-2*eta_ns*ux+p)*test(ux)-eta_ns*(uy+vx)*test(uy)-rho_ns*(u*ux+v*uy)*test(u)+gls_ns+shock_capt_ns'; ... '-eta_ns*(vx+uy)*test(vx)+(-2*eta_ns*vy+p)*test(vy)-rho_ns*(u*vx+v*vy)*test(v)'; ... '-divU_ns*test(p)';'-Xmx*test(xx)-Xmy*test(xy)';'-Ymx*test(yx)-Ymy*test(yy)'}}; equ.dweak = {{'rho_ns*ut*test(u)';'rho_ns*vt*test(v)';'0';'0';'0'}}; equ.constr = {{'0';'0';'0';'0';'0'}}; equ.constrf = {{'0';'0';'0';'0';'0'}}; equ.sshape = {[1;2;3]}; equ.sshapedim = {{2;2;2;2;2}}; equ.ind = [1,1,1,1,1,1,1,1]; equ.dim = {'u','v','p','x','y'}; equ.var = {'U_ns',{'sqrt(u^2+v^2)'}, ... 'V_ns',{'vx-uy'}, ... 'divU_ns',{'ux+vy'}, ... 'cellRe_ns',{'rho_ns*U_ns*h_ale/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_ale(u,v)),48*eta_ns/h_ale^2))'}, ... 'tauc_ns',{'0.5*nojac(if(u^2+v^2 r��L{��?��:��O��Fv��?BezierCurve��Mb��?��j!k��?BezierCurve��Mb��?��Mb`?��Mb��?BezierCurve�"�a8|��?�?BezierCurve�?��Mb`?�?BezierCurve��Mb`?��Mb��?{�G�zt?��Mb��?BezierCurve��Mb`?�?{�G�zt?�?BezierCurve{�G�zt?��Mb��?{�G�zt?��Mb��?BezierCurve{�G�zt?��Mb��?{�G�zt?�?BezierCurve�� p�k�� ?��:�xp�+҄j/݃�3��?N�2.p��E�躃�@Tf��?BezierCurve�� p�Uc�y��?��:�xp��v��y�3��?N�2.p�N� r��L{��?��:��O��Fv��?BezierCurveЈ-��O��j�C��?Ј-��O��Mb��?BezierCurveЈ-��O��Mb��?Ј-��O��?v��?BezierCurve͈-��O��Mb��?��Mb`?��Mb��?BezierCurve��Mb��?��j!k��?BezierCurve��Mb��?��Mb`?��Mb��?BezierCurve�@=��j!k��?��Mb`?��j!k��?BezierCurve�"�a8|��?�?BezierCurve�@=�"�a8|��?��Mb`?�"�a8|��?BezierCurve�?��Mb`?�?BezierCurve��Mb`?��Mb��?��Mb`?��j!k��?BezierCurve��Mb`?��Mb��?{�G�zt?��Mb��?BezierCurve��Mb`?��j!k��?��Mb`?��Mb��?BezierCurve��Mb`?��Mb��?��Mb`?�"�a8|��?BezierCurve��Mb`?��Mb��?{�G�zt?��Mb��?BezierCurve��Mb`?�"�a8|��?��Mb`?�?BezierCurve��Mb`?�?{�G�zt?�?BezierCurve{�G�zt?��Mb��?{�G�zt?��Mb��?BezierCurve{�G�zt?��Mb��?{�G�zt?�?BezierCurve�� p�k�� ?��:�xp�+҄j/݃�3��?N�2.p��E�躃�@Tf��?BezierCurve�� p�Uc�y��?��:�xp��v��y�3��?N�2.p�N�

map