[CIG-LONG] ductile

[Walter Landry]

Hi Johnny,

Cc'ing the list.

Johnny Seales <johnny.seales at gmail.com> wrote:
> I am currently working with GALE on a modeling problem that incorporates a
> ductile lower crust. I have been unable to locate examples showing an input
> file incorporating this type of layer. Do you know where I may be able to
> find one? Thanks.

I have attached simple, non-realistic one.  Is that the sort of thing
that you want, or is temperature important?

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"
        },
        "mantleShape":
        {
            "Type": "EquationShape",
            "equation": "0.1-y"
        },
        "nonmantleShape":
        {
            "Type": "Intersection",
            "shapes":
            [
                "backgroundShape",
                "!mantleShape"
            ]
        },
        "crustViscosity":
        {
            "Type": "MaterialViscosity",
            "eta0": "10.0"
        },
        "strainWeakening":
        {
            "Type": "StrainWeakening",
            "TimeIntegrator": "timeIntegrator",
            "MaterialPointsSwarm": "materialSwarm",
            "softeningStrain": "0.1",
            "initialDamageFraction": "0.0",
            "initialDamageWavenumber": "0.5",
            "initialDamageFactor": "0.5",
            "healingRate": "0.0"
        },
        "yielding":
        {
            "Type": "DruckerPrager",
            "PressureField": "PressureField",
            "VelocityGradientsField": "VelocityGradientsField",
            "MaterialPointsSwarm": "materialSwarm",
            "Context": "context",
            "StrainWeakening": "strainWeakening",
            "StrainRateField": "StrainRateField",
            "cohesion": "1.0",
            "cohesionAfterSoftening": "0.0001",
            "frictionCoefficient": "0.0",
            "frictionCoefficientAfterSoftening": "0.0",
            "minimumViscosity": "1.0e-4"
        },
        "viscous":
        {
            "Type": "RheologyMaterial",
            "Shape": "nonmantleShape",
            "density": "1.0",
            "Rheology": [
                "crustViscosity",
                "yielding",
                "storeViscosity",
                "storeStress"
            ]
        },
        "mantleViscosity":
        {
            "Type": "MaterialViscosity",
            "eta0": "1.0"
        },
        "mantleViscous":
        {
            "Type": "RheologyMaterial",
            "Shape": "mantleShape",
            "density": "1.0",
            "Rheology": [
                "mantleViscosity",
                "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": [
        "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"
}