[CIG-LONG] Isostasy in GALE

[Walter Landry]

Robert Gray <graywacke at gmail.com> wrote:
> This is the input file.
> It's modified version of deforming_bottom.json

The line

>             "normal_value": "-9.81*3300.0*(0.036e6-y)"

has the wrong sign (it should be positive).  Also, the density here is
3300, but is 2700 elsewhere.  Those two things are the main cause of
your problems.  There are other minor differences.  Gravity is 9.81
here and 9.8 elsewhere.  The height of the box is 0.036e6 here and
0.035e6 elsewhere.

With all of those changes, I get something more like the deforming
bottom.  I have attached a fixed version of the input file.

Cheers,
Walter Landry
-------------- next part --------------
{

    "EulerDeform":

    {

        "systems": [

            {

                "mesh": "v-mesh",

                "p-mesh": "p-mesh",

                "remesher": "velocityRemesher",

                "velocityField": "VelocityField",

                "wrapTop": "True",

                "staticRight":"True",

                "staticRightTop":"True",

                "staticRightBottom":"True",

                "staticLeft":"True",

                "staticLeftTop":"True",

                "staticLeftBottom":"True",

                "wrapBottom": "True"

//                "staticBottom":"True"

            }

        ]

    },

    "components":

    {

        "buoyancyForceTerm":

        {

            "Type": "BuoyancyForceTerm",

            "ForceVector": "mom_force",

            "Swarm": "gaussSwarm",

            "gravity": "gravity"

        },

        "backgroundShape":

        {

            "Type": "EquationShape",

            "equation": "1"

        },

        "backgroundViscosity":

        {

            "Type": "MaterialViscosity",

            "eta0": "1.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": "backgroundShape",

            // "density": "1.0",

            "density": "2700.0",

            "Rheology": [

                "backgroundViscosity",

//                "yielding",

                "storeViscosity",

                "storeStress"

            ]

        }

        ,"stressBC":

        {

            "Type": "StressBC",

            "ForceVector": "mom_force",

            "wall": "bottom",

            "normal_value": "9.81*2700.0*(0.035e6-y)"

//            "normal_value": "-2.3*(0.35-y)"

        }

    },

    "velocityBCs" :{

        "type": "CompositeVC",

        "vcList": [

            {

                "type": "WallVC",

                "wall": "left",

                "variables": [

                    {

                        "name": "vx",

                        "value": "0"

                    },

                    {

                        "name": "vy",

                        "value": "0"

                    }

                ]

            },

            {

                "type": "WallVC",

                "wall": "right",

                "variables": [

                    {

                        "name": "vx",

                        "value": "-4.756468798e13"

//                        "value": "-1.0"

                    },

                    {

                        "name": "vy",

                        "value": "0"

                    }

                ]

            }

        ]

    },


    "FieldVariablesToCheckpoint": [

        "StrainRateInvariantField",

        "VelocityField",

        "PressureField"

    ],

    "maxTimeSteps": "0",

    "outputPath": "./output",

    "dim": "2",

    "minX": "0",

    "minY": "0",

    "minZ": "0",

    "maxX": "1.0e6",

    "maxY": "0.035e6",

    "maxZ": "0.3",

    "nx": "16",

    "ny": "16",

    "nz": "4",

    "particlesPerCell": "40",

    "seed": "13",

    "checkpointEvery": "1",

    "gravity": "9.81"

}