[CIG-LONG] About Gale output

[Walter Landry]

feng lin <rimbong at gmail.com> wrote:
> Dear Walter,

I am Cc'ing the list in case anyone else has this question.

> I'm a phd Student major in geological science in University of Missouri. I'm
> now trying to use Gale to do build a model to study extension process.
> Gale works well. However, I couldn't get the stress tensor at the nodes. What I
> can get is only the Strain rate Invariant on the nodes or Stress Tensor on the
> particles in each element.
> Do you think I can get the stress tensor directory from the code? If I could,
> can you tell me where should I make the change?

You can get it by adding a StressField component.  As in

        "stressField":
        {
            "Type": "StressField",
            "StrainRateField":"StrainRateField",
            "Context":"context",
            "Swarm":"picIntegrationPoints",
            "Mesh":"v-mesh",
            "ConstitutiveMatrix":"constitutiveMatrix"
        },


You also need to add "stressField" to "FieldVariablesToCheckpoint".  I
am attaching a modified version of input/cookbook/multi_material.json
that outputs the stress.

Cheers,
Walter Landry
-------------- next part --------------
{
    "EulerDeform":
    {
        "systems": [
            {
                "mesh": "v-mesh",
                "p-mesh": "p-mesh",
                "remesher": "velocityRemesher",
                "velocityField": "VelocityField",
                "wrapTop": "True"
            }
        ]
    },
    "components":
    {
        "buoyancyForceTerm":
        {
            "Type": "BuoyancyForceTerm",
            "ForceVector": "mom_force",
            "Swarm": "gaussSwarm",
            "gravity": "gravity"
        },
        "backgroundShape":
        {
            "Type": "EquationShape",
            "equation": "1"
        },
        "sphereShape":
        {
            "Type": "EquationShape",
            "equation": ".1^2 - ((x-1)^2 + (y-.15)^2)"
        },
        "nonsphereShape":
        {
            "Type": "Intersection",
            "shapes":
            [
                "backgroundShape",
                "!sphereShape"
            ]
        },
        "backgroundViscosity":
        {
            "Type": "MaterialViscosity",
            "eta0": "1.0"
        },
        "viscous":
        {
            "Type": "RheologyMaterial",
            "Shape": "nonsphereShape",
            "density": "1.0",
            "Rheology": [
                "backgroundViscosity",
                "storeViscosity",
                "storeStress"
            ]
        },
        "stressField":
        {
            "Type": "StressField",
            "StrainRateField":"StrainRateField",
            "Context":"context",
            "Swarm":"picIntegrationPoints",
            "Mesh":"v-mesh",
            "ConstitutiveMatrix":"constitutiveMatrix"
        },

        "sphereViscosity":
        {
            "Type": "MaterialViscosity",
            "eta0": "10.0"
        },
        "sphereViscous":
        {
            "Type": "RheologyMaterial",
            "Shape": "sphereShape",
            "density": "1.0",
            "Rheology": [
                "sphereViscosity",
                "storeViscosity",
                "storeStress"
            ]
        }
    },
    "velocityBCs" :{
        "type": "CompositeVC",
        "vcList": [
            {
                "type": "WallVC",
                "wall": "left",
                "variables": [
                    {
                        "name": "vx",
                        "value": "0"
                    }
                ]
            },
            {
                "type": "WallVC",
                "wall": "right",
                "variables": [
                    {
                        "name": "vx",
                        "value": "1.0"
                    }
                ]
            },
            {
                "type": "WallVC",
                "wall": "bottom",
                "variables": [
                    {
                        "name": "vy",
                        "value": "0"
                    }
                ]
            }
        ]
    },

    "FieldVariablesToCheckpoint": [
        "stressField",
        "StrainRateInvariantField",
        "VelocityField",
        "PressureField"
    ],
    "maxTimeSteps": "10",
    "outputPath": "./output",
    "dim": "2",
    "minX": "0",
    "minY": "0",
    "minZ": "0",
    "maxX": "2",
    "maxY": "0.35",
    "maxZ": "0.3",
    "nx": "16",
    "ny": "4",
    "nz": "4",
    "particlesPerCell": "40",
    "seed": "13",
    "checkpointEvery": "1",
    "gravity": "1.0"
}