[CIG-LONG] defining the model bounding geometry

Walter Landry walter at geodynamics.org
Mon Dec 19 16:51:24 PST 2011


Matty Mookerjee <matty.mookerjee at sonoma.edu> wrote:
> Walter,
> 
> Thanks for the reply.  I will try using the SurfaceAdaptor function
> this week.  I should have thought of that, but for some reason I
> only thought of SurfaceAdaptor as modeling topography not defining
> the shape of the entire model.  When I give that a shot, I will send
> you my input file if the frictional boundary conditions are still
> not working.  The file will be in the old format.  I am still not
> quite comfortable with the new file format.  For instance, I have
> not been able to get a program to run with Polygons.  I have input
> pretty much exactly what you have in your revised manual, but I
> haven't yet been able to get it to run.  If I can't figure this out,
> I might send you one of those input files as well.

Did you ever manage to figure that out?  I am attaching an example
that uses a polygon instead of a sphere for the multi-material
cookbook.

> I know that you are busy with AGU, so I won't be expecting too quick
> a reply.
> 
> Good luck with all your stuff, and I will send you those files later
> this week.

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": "PolygonShape",
            "vertices": [
            "asciidata",
            ["x","y"],
            1.0, .25,
            1.0707, .2207,
            1.1, .15,
            1.0707, .07929,
            1.0, .05,
            .92929, .07929,
            .9, .15,
            .92929, .2207
            ]
        },
        "nonsphereShape":
        {
            "Type": "Intersection",
            "shapes":
            [
                "backgroundShape",
                "!sphereShape"
            ]
        },
        "backgroundViscosity":
        {
            "Type": "MaterialViscosity",
            "eta0": "1.0"
        },
        "viscous":
        {
            "Type": "RheologyMaterial",
            "Shape": "nonsphereShape",
            "density": "1.0",
            "Rheology": [
                "backgroundViscosity",
                "storeViscosity",
                "storeStress"
            ]
        },
        "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": [
        "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"
}


More information about the CIG-LONG mailing list