{ Fill in the following sections (removing comment marks ! if necessary), and delete those that are unused.} TITLE 'New Problem' { the problem identification } COORDINATES cartesian3 { coordinate system, 1D,2D,3D, etc } VARIABLES { system variables } ux uy uz vfx vfy vfz hp SELECT regrid=off deltat = 1e-4 ! initial timestep, choose 1e-6 or lower for PFM smoothinit=off prefer_stability=off prefer_speed=off errlim = 1.0e-3 threads = 6 nominmax notags painted DEFINITIONS mat lam_0s mue_0s ns0s cfc0s mue_fr ks0s jacobi = dx(ux) + dy(uy) + dz(uz) + dx(ux)*dy(uy) - (dx(uy)*dy(ux)) + dx(ux)*dz(uz) - (dx(uz)*dz(ux)) + dy(uy)*dz(uz) - (dy(uz)*dz(uy)) + dx(ux)*dy(uy)*dz(uz) - (dx(ux)*dy(uz)*dz(uy)) - (dx(uy)*dy(ux)*dz(uz)) + dx(uy)*dy(uz)*dz(ux) + dx(uz)*dy(ux)*dz(uy) - (dx(uz)*dy(uy)*dz(ux)) + 1 ns = ns0s*(1/(dx(ux) + dy(uy) + dz(uz) + dx(ux)*dy(uy) - (dx(uy)*dy(ux)) + dx(ux)*dz(uz) - (dx(uz)*dz(ux)) + dy(uy)*dz(uz) - (dy(uz)*dz(uy)) + dx(ux)*dy(uy)*dz(uz) - (dx(ux)*dy(uz)*dz(uy)) - (dx(uy)*dy(ux)*dz(uz)) + dx(uy)*dy(uz)*dz(ux) + dx(uz)*dy(ux)*dz(uy) - (dx(uz)*dy(uy)*dz(ux)) + 1)) nf = 1-ns ks = ks0s R = 8.314 temp_osm_deg = 25 temp_osm = 298.15 + temp_osm_deg cm_ = 0.15 alpha=1e3 nf0s=1-ns0s cfcm=(cfc0s*nf0s)/(jacobi-ns0s) vx=dt(ux) vy=dt(uy) vz=dt(uz) vv=vector(vx,vy,vz) vf=vector(vfx,vfy,vfz) !--- strain tensor (Isotropic) epsxx = dx(ux) epsyy = dy(uy) epszz = dz(uz) epsxy = 0.5*(dx(uy)+dy(ux)) epsyz = 0.5*(dy(uz)+dz(uy)) epsxz = 0.5*(dx(uz)+dz(ux)) epsyx = epsxy epszy = epsyz epszx = epsxz treps = epsxx+epsyy+epszz !--- isotropic stress tensor sigxx = 2*lam_0s*epsxx + mue_0s*treps sigyy = 2*lam_0s*epsyy + mue_0s*treps sigzz = 2*lam_0s*epszz + mue_0s*treps sigxy = 2*lam_0s*epsxy sigyz = 2*lam_0s*epsyz sigxz = 2*lam_0s*epsxz sigyx = 2*lam_0s*epsyx sigzx = 2*lam_0s*epszx sigzy =2*lam_0s*epszy !--- delta_pi delta_pi=R*temp_osm*(sqrt(4*cm_^2+cfcm^2)-(2*cm_))*0.001 !--- load P = ramp(t, 0, 2, 0.1) INITIAL VALUES ux = 0 uy = 0 uz = 0 vfx = 0 vfy = 0 vfz = 0 hp = 0 EQUATIONS ux: dx(sigxx) + dy(sigyx) + dz(sigzx) - dx(hp+delta_pi)= 0 uy: dx(sigxy) + dy(sigyy) + dz(sigzy) - dy(hp+delta_pi)= 0 uz: dx(sigxz) + dy(sigyz) + dz(sigzz) - dz(hp+delta_pi)= 0 vfx: - dx(hp) - nf*mue_fr/ks*(vfx-vx) = 0 vfy: - dy(hp) - nf*mue_fr/ks*(vfy-vy)= 0 vfz: - dz(hp) - nf*mue_fr/ks*(vfz-vz) = 0 hp: div(ns*vv+nf*vf) - alpha*del2(hp) = 0 extrusion z= -3.1,-3,-2,0,0.5,0.6 BOUNDARIES surface 1 value(ux)=0 value(uy)=0 value(uz)=0 value(hp)=0 surface 6 value(ux)=0 value(uy)=0 load(uz)=P value(hp)=0 Limited REGION 1 mat = 1 mue_fr = 3.8e-8 mue_0s = 96 lam_0s= 112.7 ns0s=0.2 ks0s=2.7e-5 cfc0s=0 Layer 1 Layer 5 START(-1,-1) value(hp)=0 LINE TO (-1,1) value(hp)=0 line TO (1,1) value(hp)=0 line TO (1,-1) value(hp)=0 Line to close Limited Region 2 mat = 2 mue_fr=6.9e-10 mue_0s = 0.95 lam_0s= 2.2 ns0s=0.35 ks0s=6.2e-12 cfc0s=0.1 Layer 2 Layer 3 Layer 4 START(-1,-1) value(hp)=0 LINE TO (-1,1) value(hp)=0 line TO (1,1) value(hp)=0 line TO (1,-1) value(hp)=0 Line to close Limited Region 3 mat = 3 mue_fr=6.9e-10 mue_0s=0.5 lam_0s = 0.75 ns0s=0.2 ks0s=3.5e-12 cfc0s=0.3 Layer 2 Layer 3 Layer 4 START(-0.5,-0.5) LINE TO (-0.5,0.5) line TO (0.5,0.5) line TO (0.5,-0.5) Line to close TIME 0 TO 1 PLOTS for cycle = 1 grid(x+ux,y+uy,z+uz) grid(x+ux,z+uz) on y=0 contour(hp) on y=0 contour(hp) on z=-2 contour(uz) on y=0 END