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>    &lt;struct name=&quot;stokesEqn&quot;&gt;</div><div>      &lt;param name=&quot;Type&quot;&gt;Stokes_SLE&lt;/param&gt;</div><div>      &lt;param name=&quot;SLE_Solver&quot;&gt;uzawa&lt;/param&gt;</div>
<div>      &lt;param name=&quot;Context&quot;&gt;context&lt;/param&gt;</div><div>      &lt;param name=&quot;StressTensorMatrix&quot;&gt;k_matrix&lt;/param&gt;</div><div>      &lt;param name=&quot;GradientMatrix&quot;&gt;g_matrix&lt;/param&gt;</div>
<div>      &lt;param name=&quot;DivergenceMatrix&quot;&gt;&lt;/param&gt;</div><div>      &lt;param name=&quot;CompressibilityMatrix&quot;&gt;c_matrix&lt;/param&gt; &lt;-- without this line the model run, bit somehow compressibility is not working</div>
<div>      &lt;param name=&quot;VelocityVector&quot;&gt;solutionVelocity&lt;/param&gt;</div><div>      &lt;param name=&quot;PressureVector&quot;&gt;solutionPressure&lt;/param&gt;</div><div>      &lt;param name=&quot;ForceVector&quot;&gt;mom_force&lt;/param&gt;</div>
<div>      &lt;param name=&quot;ContinuityForceVector&quot;&gt;cont_force&lt;/param&gt;</div><div>      &lt;param name=&quot;nonLinearMaxIterations&quot;&gt;nonLinearMaxIterations&lt;/param&gt;</div><div>      &lt;param name=&quot;nonLinearTolerance&quot;&gt;nonLinearTolerance&lt;/param&gt;</div>
<div>      &lt;param name=&quot;makeConvergenceFile&quot;&gt;false&lt;/param&gt;</div><div>    &lt;/struct&gt;</div><div>    &lt;struct name=&quot;c_matrix&quot;&gt;</div><div>      &lt;param name=&quot;Type&quot;&gt;StiffnessMatrix&lt;/param&gt;</div>
<div>      &lt;param name=&quot;RowVariable&quot;&gt;PressureField&lt;/param&gt;</div><div>      &lt;param name=&quot;ColumnVariable&quot;&gt;PressureField&lt;/param&gt;</div><div>      &lt;param name=&quot;RHS&quot;&gt;cont_force&lt;/param&gt;</div>
<div>      &lt;param name=&quot;allowZeroElementContributions&quot;&gt;True&lt;/param&gt;</div><div>    &lt;/struct&gt;</div><div>    &lt;struct name=&quot;mixedStabiliser&quot;&gt;</div><div>      &lt;param name=&quot;Type&quot;&gt;MixedStabiliserTerm&lt;/param&gt;</div>
<div>      &lt;param name=&quot;Swarm&quot;&gt;gaussSwarm&lt;/param&gt;</div><div>      &lt;param name=&quot;picSwarm&quot;&gt;picIntegrationPoints&lt;/param&gt;</div><div>      &lt;param name=&quot;storeVisc&quot;&gt;storeViscosity&lt;/param&gt;</div>
<div>      &lt;param name=&quot;StiffnessMatrix&quot;&gt;c_matrix&lt;/param&gt;</div><div>    &lt;/struct&gt;</div></div><div><br></div><div>and the Air layer:</div><div><br></div><div><div>    &lt;struct name=&quot;airShape&quot;&gt;</div>
<div>      &lt;param name=&quot;Type&quot;&gt;Box&lt;/param&gt;</div><div>      &lt;param name=&quot;startX&quot;&gt;minX&lt;/param&gt;</div><div>      &lt;param name=&quot;endX&quot;&gt;maxX&lt;/param&gt;</div><div>      &lt;param name=&quot;startY&quot;&gt;420000&lt;/param&gt;</div>
<div>      &lt;param name=&quot;endY&quot;&gt;maxAirY&lt;/param&gt;</div><div>      &lt;param name=&quot;startZ&quot;&gt;minZ&lt;/param&gt;</div><div>      &lt;param name=&quot;endZ&quot;&gt;maxZ&lt;/param&gt;</div><div>    &lt;/struct&gt;</div>
<div><span class="Apple-tab-span" style="white-space:pre">        </span>&lt;struct name=&quot;airViscosity&quot;&gt;</div><div>      &lt;param name=&quot;Type&quot;&gt;MaterialViscosity&lt;/param&gt;</div><div>      &lt;param name=&quot;eta0&quot;&gt;1e-4&lt;/param&gt;</div>
<div>    &lt;/struct&gt;</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>&lt;struct name=&quot;compresibility&quot;&gt;</div><div>      &lt;param name=&quot;Type&quot;&gt;Compressible&lt;/param&gt;</div>
<div>      &lt;param name=&quot;GeometryMesh&quot;&gt;v-mesh&lt;/param&gt;  &lt;-- I changed the mesh-linear to the v-mesh</div><div>      &lt;param name=&quot;StiffnessMatrix&quot;&gt;c_matrix&lt;/param&gt;</div><div>      &lt;param name=&quot;Swarm&quot;&gt;gaussSwarm&lt;/param&gt;</div>
<div>      &lt;param name=&quot;oneOnLambda&quot;&gt;10.0&lt;/param&gt;</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>&lt;/struct&gt;</div><div><br></div><div>    &lt;struct name=&quot;air&quot;&gt;</div>
<div>      &lt;param name=&quot;Type&quot;&gt;RheologyMaterial&lt;/param&gt;</div><div>      &lt;param name=&quot;Shape&quot;&gt;airShape&lt;/param&gt;</div><div>      &lt;param name=&quot;density&quot;&gt;0&lt;/param&gt;</div>
<div>      &lt;param name=&quot;diffusivity&quot;&gt;1e18&lt;/param&gt;</div><div>      &lt;list name=&quot;heatingElements&quot;&gt;</div><div>        &lt;struct&gt;</div><div>          &lt;param name=&quot;Q&quot;&gt;0&lt;/param&gt;</div>
<div>          &lt;param name=&quot;lambda&quot;&gt;0.0&lt;/param&gt;</div><div>        &lt;/struct&gt;</div><div>      &lt;/list&gt;</div><div>      &lt;list name=&quot;Rheology&quot;&gt;</div><div>        &lt;param&gt;airViscosity&lt;/param&gt;</div>
<div>        &lt;param&gt;storeViscosity&lt;/param&gt;</div><div>        &lt;param&gt;storeStress&lt;/param&gt;</div><div>      &lt;/list&gt;</div><div>      &lt;param name=&quot;Compressible&quot;&gt;compresibility&lt;/param&gt;</div>
<div>   &lt;/struct&gt;<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>