[CIG-LONG] Isostasy in GALE

[Robert Gray]

This is the input file.
It's modified version of deforming_bottom.json
Thanks,
Rob


{

    "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*3300.0*(0.036e6-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": "20",

    "outputPath": "./output",

    "dim": "2",

    "minX": "0",

    "minY": "0",

    "minZ": "0",

    "maxX": "1.0e6",

    "maxY": "0.035e6",

    "maxZ": "0.3",

    "nx": "16",

    "ny": "4",

    "nz": "4",

    "particlesPerCell": "40",

    "seed": "13",

    "checkpointEvery": "1",

    "gravity": "9.8"

}



On Tue, Feb 28, 2012 at 4:07 PM, Walter Landry <walter at geodynamics.org>wrote:

> Robert Gray <graywacke at gmail.com> wrote:
> > Hello all,
> > I'm trying to model crustal scale collision where the bottom boundary
> > deflects downwards to reflect isostatic deflection due to a crustal root,
> > as would be done in MICROFEM.  That said, the entire crustal layer
> deflects
> > downwards way too much.  Has anyone tried modelling anything similar?
> > Instead of velocity BC's on the bottom boundary I've put:
> >  ,"stressBC":
> >         {
> >             "Type": "StressBC",
> >             "ForceVector": "mom_force",
> >             "wall": "bottom",
> >             "normal_value": "-9.81*3300.0*(0.036e6-y)"
> >         }
> > at the end of the components section.
>
> Do you have a non-zero temperature?  Otherwise, that should work.  Can
> you send your complete input file?
>
> Cheers,
> Walter
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-long/attachments/20120228/54e8b462/attachment-0001.htm