[cig-commits] commit: Add subduction
Mercurial
hg at geodynamics.org
Mon Jul 2 17:24:31 PDT 2012
changeset: 457:99a12f80251b
tag: tip
user: Walter Landry <wlandry at caltech.edu>
date: Mon Jul 02 17:24:12 2012 -0700
files: input/geologic/subduction.json
description:
Add subduction
diff -r 786b333eafa2 -r 99a12f80251b input/geologic/subduction.json
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/input/geologic/subduction.json Mon Jul 02 17:24:12 2012 -0700
@@ -0,0 +1,378 @@
+{
+ "EulerDeform":
+ {
+ "systems": [
+ {
+ "mesh": "v-mesh",
+ "p-mesh": "p-mesh",
+ "displacementField": "DisplacementField",
+ "remesher": "velocityRemesher",
+ "velocityField": "VelocityField",
+ "wrapTop": "True",
+ "staticBottom": "True",
+ "staticLeft": "True",
+ "staticRight": "True",
+ "left_equation": "0.01/(365*86400)*1e21 * t",
+ "right_equation": "3200000 - 0.01/(365*86400)*1e21 * t"
+ }
+ ]
+ },
+ "components":
+ {
+ "buoyancyForceTerm":
+ {
+ "Type": "BuoyancyForceTerm",
+ "ForceVector": "mom_force",
+ "Swarm": "gaussSwarm",
+ "gravity": "gravity",
+ "TemperatureField": "TemperatureField"
+ },
+ "backgroundShape":
+ {
+ "Type": "EquationShape",
+ "equation": "1"
+ },
+ "upperCrustShape":
+ {
+ "Type": "EquationShape",
+ "equation": "x<800000 ? y-(800000-15000) : -1"
+ },
+ "lowerCrustShape":
+ {
+ "Type": "EquationShape",
+ "equation": "x<800000 ? ((800000-15000)-y)*(y-(800000-45000)) : -1"
+ },
+ "leftLithosphericMantleShape":
+ {
+ "Type": "EquationShape",
+ "equation": "x < 800000 ? (y< 755000 ? (y > 700000 ? 1 : -1) : -1) : -1"
+ },
+ "rightLithosphericMantleShape":
+ {
+ "Type": "EquationShape",
+ "equation": "x < 950000 ? (x < 900000 ? (y > 700000 ? (y-700000) - (x-800000) : -1) : (y>700000 ? (x-900000) - (y-700000) : -1)) : (y > 700000 ? 1 : -1)"
+ },
+ // "damaged_lithospheric_mantle":
+ // {
+ // "Type": "EquationShape",
+ // "equation": "x < 900000 ? -1 : (y > 700000 ? 1 : -1)"
+ // },
+ "mantleShape":
+ {
+ "Type": "Intersection",
+ "shapes":
+ [
+ "backgroundShape",
+ "!upperCrustShape",
+ "!lowerCrustShape",
+ "!rightLithosphericMantleShape",
+ "!leftLithosphericMantleShape"
+ ]
+ },
+ "nonNewtonian_crust":
+ {
+ "Type": "NonNewtonian",
+ "StrainRateInvariantField": "StrainRateInvariantField",
+ "TemperatureField": "TemperatureField",
+ "n": "3",
+ "T_0": "51717.1158914",
+ "A": "1.7578e7",
+ "refStrainRate": "396372.399798",
+ "maxViscosity": "1e4",
+ "minViscosity": "1e0"
+ },
+ // Based on Karato & Wu 1993, Dry Olivine
+ "nonNewtonian_lithospheric_mantle":
+ {
+ "Type": "NonNewtonian",
+ "StrainRateInvariantField": "StrainRateInvariantField",
+ "TemperatureField": "TemperatureField",
+ "n": "3.5",
+ "T_0": "64947",
+ "A": "2.1e-16",
+ "refStrainRate": "396372.399798",
+ "maxViscosity": "1e4",
+ "minViscosity": "1e0"
+ },
+ "strainWeakening":
+ {
+ "Type": "StrainWeakening",
+ "TimeIntegrator": "timeIntegrator",
+ "MaterialPointsSwarm": "materialSwarm",
+ "initialSofteningStrain": "0.5",
+ "finalSofteningStrain": "1.5",
+ "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": "20e6",
+ "cohesionAfterSoftening": "1e6",
+ "frictionCoefficient": "0.258819045102",
+ "frictionCoefficientAfterSoftening": "0.0348994967025",
+ "minimumViscosity": "1.0e0"
+ },
+ "strainWeakening_damaged":
+ {
+ "Type": "StrainWeakening",
+ "TimeIntegrator": "timeIntegrator",
+ "MaterialPointsSwarm": "materialSwarm",
+ "initialSofteningStrain": "0.5",
+ "finalSofteningStrain": "1.5",
+ "initialDamageFraction": "0.0",
+ "initialDamageWavenumber": "0.5",
+ "initialDamageFactor": "0.5",
+ "healingRate": "0.0"
+ },
+ "yielding_damaged":
+ {
+ "Type": "DruckerPrager",
+ "PressureField": "PressureField",
+ "VelocityGradientsField": "VelocityGradientsField",
+ "MaterialPointsSwarm": "materialSwarm",
+ "Context": "context",
+ "StrainWeakening": "strainWeakening_damaged",
+ "StrainRateField": "StrainRateField",
+ "cohesion": "20e6",
+ "cohesionAfterSoftening": "1e6",
+ "frictionCoefficient": "0.258819045102",
+ "frictionCoefficientAfterSoftening": "0.0348994967025",
+ "minimumViscosity": "1.0e0"
+ },
+ "fixedViscosity":
+ {
+ "Type": "MaterialViscosity",
+ "eta0": "1.0e0"
+ },
+ "mantle":
+ {
+ "Type": "RheologyMaterial",
+ "Shape": "mantleShape",
+ "density": "3300.0",
+ "alpha": "1.0e-5",
+ "diffusivity": "1e15",
+ "heatingElements": [
+ {
+ "Q": "0.0",
+ "lambda": "1.0"
+ }
+ ],
+ "Rheology": [
+ "fixedViscosity",
+ "storeViscosity",
+ "storeStress"
+ ]
+ },
+ "rightLithosphericMantle":
+ {
+ "Type": "RheologyMaterial",
+ "Shape": "rightLithosphericMantleShape",
+ "density": "3300.0",
+ "alpha": "1.0e-5",
+ "diffusivity": "1e15",
+ "heatingElements": [
+ {
+ "Q": "0.0",
+ "lambda": "1.0"
+ }
+ ],
+ "Rheology": [
+ "nonNewtonian_lithospheric_mantle",
+ "yielding",
+ "storeViscosity",
+ "storeStress"
+ ]
+ },
+ "leftLithosphericMantle":
+ {
+ "Type": "RheologyMaterial",
+ "Shape": "leftLithosphericMantleShape",
+ "density": "3300.0",
+ "alpha": "1.0e-5",
+ "diffusivity": "1e15",
+ "heatingElements": [
+ {
+ "Q": "0.0",
+ "lambda": "1.0"
+ }
+ ],
+ "Rheology": [
+ "nonNewtonian_lithospheric_mantle",
+ "yielding",
+ "storeViscosity",
+ "storeStress"
+ ]
+ },
+ "upperCrust":
+ {
+ "Type": "RheologyMaterial",
+ "Shape": "upperCrustShape",
+ "density": "2700.0",
+ "alpha": "1.0e-5",
+ "diffusivity": "1e15",
+ "heatingElements": [
+ {
+ "Q": "0.0",
+ "lambda": "1.0"
+ }
+ ],
+ "Rheology": [
+ "nonNewtonian_crust",
+ "yielding",
+ "storeViscosity",
+ "storeStress"
+ ]
+ },
+ "lowerCrust":
+ {
+ "Type": "RheologyMaterial",
+ "Shape": "lowerCrustShape",
+ "density": "2800.0",
+ "alpha": "1.0e-5",
+ "diffusivity": "1e15",
+ "heatingElements": [
+ {
+ "Q": "0.0",
+ "lambda": "1.0"
+ }
+ ],
+ "Rheology": [
+ "nonNewtonian_crust",
+ "yielding",
+ "storeViscosity",
+ "storeStress"
+ ]
+ },
+ "stressBC":
+ {
+ "Type": "StressBC",
+ "ForceVector": "mom_force",
+ "wall": "bottom",
+ "normal_value": "9.81*(2700*15000 + 2800*20000 + 3300*765000)*0.98"
+ }
+ },
+ "velocityBCs" :{
+ "type": "CompositeVC",
+ "vcList": [
+ {
+ "type": "WallVC",
+ "wall": "left",
+ "variables": [
+ {
+ "name": "vx",
+ "value": "0.01/(365*86400)*1e21*(y<400000 ? 0 : (y<700000 ? (y-400000)/300000 : 1))"
+ }
+ ]
+ },
+ {
+ "type": "WallVC",
+ "wall": "right",
+ "variables": [
+ {
+ "name": "vx",
+ "value": "-0.01/(365*86400)*1e21*(y<400000 ? 0 : (y<700000 ? (y-400000)/300000 : 1))"
+ }
+ ]
+ },
+ {
+ "type": "WallVC",
+ "wall": "right",
+ "variables": [
+ {
+ "name": "vy",
+ "value": "0"
+ }
+ ]
+ },
+ {
+ "type": "WallVC",
+ "wall": "left",
+ "variables": [
+ {
+ "name": "vy",
+ "value": "0"
+ }
+ ]
+ }
+ ]
+ },
+ "temperatureBCs": {
+ "type": "CompositeVC",
+ "vcList": [
+ {
+ "type": "WallVC",
+ "wall": "top",
+ "variables": [
+ {
+ "name": "temperature",
+ "value": "273.0"
+ }
+ ]
+ },
+ {
+ "type": "WallVC",
+ "wall": "bottom",
+ "variables": [
+ {
+ "name": "temperature",
+ "value": "1660.0"
+ }
+ ]
+ }
+ ]
+ },
+ "temperatureICs":
+ {
+ "type": "CompositeVC",
+ "vcList": [
+ {
+ "type": "AllNodesVC",
+ "variables":[
+ {
+ "name": "temperature",
+ "value": "y>700000 ? (800000-y)/(100000)*(1660-273) + 273 : 1660"
+ }
+ ]
+ }
+ ]
+ },
+
+ "FieldVariablesToCheckpoint": [
+ "StrainRateInvariantField",
+ "VelocityField",
+ "TemperatureField",
+ "PressureField"
+ ],
+ "enable-thermal": true,
+ "maxTimeSteps": "100",
+ "outputPath": "./subduction",
+ "dim": "2",
+ "minX": "0",
+ "minY": "0",
+ "minZ": "0",
+ "maxX": "3200000",
+ "maxY": "800000",
+ "maxZ": "0.3",
+ "nx": "64",
+ "ny": "32",
+ "nz": "4",
+ "particlesPerCell": "40",
+ "seed": "13",
+ "checkpointEvery": "1",
+ "gravity": "9.81",
+ "dtFactor": "0.25"
+// ,
+// "journal.info": "True",
+// "journal.debug": "True",
+// "journal-level.info": "1",
+// "journal-level.debug": "1"
+}
\ No newline at end of file
More information about the CIG-COMMITS
mailing list