Hi everyone,<div><br></div><div>I have a little question about how to add a compressible layer with Gale 2.0.1.</div><div><br></div><div><br></div><div>I added these lines (like for Gale 1.6.1) in the input.xml file created in the output directory.</div>
<div><br></div><div><div> <struct name="stokesEqn"></div><div> <param name="Type">Stokes_SLE</param></div><div> <param name="SLE_Solver">uzawa</param></div>
<div> <param name="Context">context</param></div><div> <param name="StressTensorMatrix">k_matrix</param></div><div> <param name="GradientMatrix">g_matrix</param></div>
<div> <param name="DivergenceMatrix"></param></div><div> <param name="CompressibilityMatrix">c_matrix</param> <-- without this line the model run, bit somehow compressibility is not working</div>
<div> <param name="VelocityVector">solutionVelocity</param></div><div> <param name="PressureVector">solutionPressure</param></div><div> <param name="ForceVector">mom_force</param></div>
<div> <param name="ContinuityForceVector">cont_force</param></div><div> <param name="nonLinearMaxIterations">nonLinearMaxIterations</param></div><div> <param name="nonLinearTolerance">nonLinearTolerance</param></div>
<div> <param name="makeConvergenceFile">false</param></div><div> </struct></div><div> <struct name="c_matrix"></div><div> <param name="Type">StiffnessMatrix</param></div>
<div> <param name="RowVariable">PressureField</param></div><div> <param name="ColumnVariable">PressureField</param></div><div> <param name="RHS">cont_force</param></div>
<div> <param name="allowZeroElementContributions">True</param></div><div> </struct></div><div> <struct name="mixedStabiliser"></div><div> <param name="Type">MixedStabiliserTerm</param></div>
<div> <param name="Swarm">gaussSwarm</param></div><div> <param name="picSwarm">picIntegrationPoints</param></div><div> <param name="storeVisc">storeViscosity</param></div>
<div> <param name="StiffnessMatrix">c_matrix</param></div><div> </struct></div></div><div><br></div><div>and the Air layer:</div><div><br></div><div><div> <struct name="airShape"></div>
<div> <param name="Type">Box</param></div><div> <param name="startX">minX</param></div><div> <param name="endX">maxX</param></div><div> <param name="startY">420000</param></div>
<div> <param name="endY">maxAirY</param></div><div> <param name="startZ">minZ</param></div><div> <param name="endZ">maxZ</param></div><div> </struct></div>
<div><span class="Apple-tab-span" style="white-space:pre">        </span><struct name="airViscosity"></div><div> <param name="Type">MaterialViscosity</param></div><div> <param name="eta0">1e-4</param></div>
<div> </struct></div><div><span class="Apple-tab-span" style="white-space:pre">        </span><struct name="compresibility"></div><div> <param name="Type">Compressible</param></div>
<div> <param name="GeometryMesh">v-mesh</param> <-- I changed the mesh-linear to the v-mesh</div><div> <param name="StiffnessMatrix">c_matrix</param></div><div> <param name="Swarm">gaussSwarm</param></div>
<div> <param name="oneOnLambda">10.0</param></div><div><span class="Apple-tab-span" style="white-space:pre">        </span></struct></div><div><br></div><div> <struct name="air"></div>
<div> <param name="Type">RheologyMaterial</param></div><div> <param name="Shape">airShape</param></div><div> <param name="density">0</param></div>
<div> <param name="diffusivity">1e18</param></div><div> <list name="heatingElements"></div><div> <struct></div><div> <param name="Q">0</param></div>
<div> <param name="lambda">0.0</param></div><div> </struct></div><div> </list></div><div> <list name="Rheology"></div><div> <param>airViscosity</param></div>
<div> <param>storeViscosity</param></div><div> <param>storeStress</param></div><div> </list></div><div> <param name="Compressible">compresibility</param></div>
<div> </struct><span class="Apple-tab-span" style="white-space:pre">        </span></div><div><span class="Apple-tab-span" style="white-space:pre">        </span></div></div><div>However there is segmentation fault error...</div>
<div><br></div><div>Is anyone knows how to use compressible layers ?</div><div><br></div><div><br></div><div>Thank you in advance,</div><div><br></div><div>Nicolas Riel</div><div><br></div><div><br></div>