[CIG-LONG] Compressible layer

[Nicolas.Riel]

Hi everyone,

I have a little question about how to add a compressible layer with Gale
2.0.1.


I added these lines (like for Gale 1.6.1) in the input.xml file created in
the output directory.

    <struct name="stokesEqn">
      <param name="Type">Stokes_SLE</param>
      <param name="SLE_Solver">uzawa</param>
      <param name="Context">context</param>
      <param name="StressTensorMatrix">k_matrix</param>
      <param name="GradientMatrix">g_matrix</param>
      <param name="DivergenceMatrix"></param>
      <param name="CompressibilityMatrix">c_matrix</param> <-- without this
line the model run, bit somehow compressibility is not working
      <param name="VelocityVector">solutionVelocity</param>
      <param name="PressureVector">solutionPressure</param>
      <param name="ForceVector">mom_force</param>
      <param name="ContinuityForceVector">cont_force</param>
      <param name="nonLinearMaxIterations">nonLinearMaxIterations</param>
      <param name="nonLinearTolerance">nonLinearTolerance</param>
      <param name="makeConvergenceFile">false</param>
    </struct>
    <struct name="c_matrix">
      <param name="Type">StiffnessMatrix</param>
      <param name="RowVariable">PressureField</param>
      <param name="ColumnVariable">PressureField</param>
      <param name="RHS">cont_force</param>
      <param name="allowZeroElementContributions">True</param>
    </struct>
    <struct name="mixedStabiliser">
      <param name="Type">MixedStabiliserTerm</param>
      <param name="Swarm">gaussSwarm</param>
      <param name="picSwarm">picIntegrationPoints</param>
      <param name="storeVisc">storeViscosity</param>
      <param name="StiffnessMatrix">c_matrix</param>
    </struct>

and the Air layer:

    <struct name="airShape">
      <param name="Type">Box</param>
      <param name="startX">minX</param>
      <param name="endX">maxX</param>
      <param name="startY">420000</param>
      <param name="endY">maxAirY</param>
      <param name="startZ">minZ</param>
      <param name="endZ">maxZ</param>
    </struct>
<struct name="airViscosity">
      <param name="Type">MaterialViscosity</param>
      <param name="eta0">1e-4</param>
    </struct>
<struct name="compresibility">
      <param name="Type">Compressible</param>
      <param name="GeometryMesh">v-mesh</param>  <-- I changed the
mesh-linear to the v-mesh
      <param name="StiffnessMatrix">c_matrix</param>
      <param name="Swarm">gaussSwarm</param>
      <param name="oneOnLambda">10.0</param>
</struct>

    <struct name="air">
      <param name="Type">RheologyMaterial</param>
      <param name="Shape">airShape</param>
      <param name="density">0</param>
      <param name="diffusivity">1e18</param>
      <list name="heatingElements">
        <struct>
          <param name="Q">0</param>
          <param name="lambda">0.0</param>
        </struct>
      </list>
      <list name="Rheology">
        <param>airViscosity</param>
        <param>storeViscosity</param>
        <param>storeStress</param>
      </list>
      <param name="Compressible">compresibility</param>
   </struct>
However there is segmentation fault error...

Is anyone knows how to use compressible layers ?


Thank you in advance,

Nicolas Riel
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-long/attachments/20130325/b577aa32/attachment.htm