[cig-commits] commit: Add geologic extension
Mercurial
hg at geodynamics.org
Wed Jun 20 14:33:22 PDT 2012
changeset: 453:b035ceed4a7b
tag: tip
user: Walter Landry <wlandry at caltech.edu>
date: Wed Jun 20 14:33:04 2012 -0700
files: input/geologic/extension.json
description:
Add geologic extension
diff -r 83ebbcc4b06e -r b035ceed4a7b input/geologic/extension.json
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/input/geologic/extension.json Wed Jun 20 14:33:04 2012 -0700
@@ -0,0 +1,294 @@
+{
+ "EulerDeform":
+ {
+ "systems": [
+ {
+ "mesh": "v-mesh",
+ "p-mesh": "p-mesh",
+ "displacementField": "DisplacementField",
+ "remesher": "velocityRemesher",
+ "velocityField": "VelocityField",
+ "wrapTop": "True",
+ "staticBottom": "True"
+ }
+ ]
+ },
+ "components":
+ {
+ "buoyancyForceTerm":
+ {
+ "Type": "BuoyancyForceTerm",
+ "ForceVector": "mom_force",
+ "Swarm": "gaussSwarm",
+ "gravity": "gravity",
+ "TemperatureField": "TemperatureField"
+ },
+ "backgroundShape":
+ {
+ "Type": "EquationShape",
+ "equation": "1"
+ },
+ "upperCrustShape":
+ {
+ "Type": "EquationShape",
+ "equation": "y-(200000-15000)"
+ },
+ "magmaShape":
+ {
+ "Type": "EquationShape",
+ "equation": "1^2 - (((x-400000)/25000)^2 + ((y-(200000-35000))/5000)^2)"
+ },
+ "lowerCrustLayerShape":
+ {
+ "Type": "EquationShape",
+ "equation": "((200000-15000)-y)*(y-(200000-45000))"
+ },
+ "lowerCrustShape":
+ {
+ "Type": "Intersection",
+ "shapes":
+ [
+ "lowerCrustLayerShape",
+ "!magmaShape"
+ ]
+ },
+ "mantleShape":
+ {
+ "Type": "Intersection",
+ "shapes":
+ [
+ "backgroundShape",
+ "!upperCrustShape",
+ "!lowerCrustShape"
+ ]
+ },
+ "nonNewtonian":
+ {
+ "Type": "NonNewtonian",
+ "StrainRateInvariantField": "StrainRateInvariantField",
+ "TemperatureField": "TemperatureField",
+ "n": "3",
+ "T_0": "51717.1158914",
+ "A": "1.7578e7",
+ "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"
+ },
+ "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": [
+ "nonNewtonian",
+ "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",
+ "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",
+ "storeViscosity",
+ "storeStress"
+ ]
+ },
+ "magma":
+ {
+ "Type": "RheologyMaterial",
+ "Shape": "magmaShape",
+ "density": "3000.0",
+ "alpha": "1.0e-5",
+ "diffusivity": "1e15",
+ "heatingElements": [
+ {
+ "Q": "0.0",
+ "lambda": "1.0"
+ }
+ ],
+ "Rheology": [
+ "fixedViscosity",
+ "storeViscosity",
+ "storeStress"
+ ]
+ },
+ "stressBC":
+ {
+ "Type": "StressBC",
+ "ForceVector": "mom_force",
+ "wall": "bottom",
+ "normal_value": "9.81*(2700*15000 + 2800*20000 + 3300*165000)"
+ }
+ },
+ "velocityBCs" :{
+ "type": "CompositeVC",
+ "vcList": [
+ {
+ "type": "WallVC",
+ "wall": "left",
+ "variables": [
+ {
+ "name": "vx",
+ "value": "-0.01/(365*86400)*1e21"
+ }
+ ]
+ },
+ {
+ "type": "WallVC",
+ "wall": "right",
+ "variables": [
+ {
+ "name": "vx",
+ "value": "0.01/(365*86400)*1e21"
+ }
+ ]
+ },
+ {
+ "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": "bottom",
+ "variables": [
+ {
+ "name": "temperature",
+ "value": "1660.0"
+ }
+ ]
+ }
+ ]
+ },
+ "temperatureICs":
+ {
+ "type": "CompositeVC",
+ "vcList": [
+ {
+ "type": "AllNodesVC",
+ "variables":[
+ {
+ "name": "temperature",
+ "value": "y>100000 ? (200000-y)/(100000)*(1660-273) + 273 : 1660"
+ }
+ ]
+ }
+ ]
+ },
+
+ "FieldVariablesToCheckpoint": [
+ "StrainRateInvariantField",
+ "VelocityField",
+ "TemperatureField",
+ "PressureField"
+ ],
+ "enable-thermal": true,
+ "maxTimeSteps": "100",
+ "outputPath": "./extension",
+ "dim": "2",
+ "minX": "0",
+ "minY": "0",
+ "minZ": "0",
+ "maxX": "800000",
+ "maxY": "200000",
+ "maxZ": "0.3",
+ "nx": "64",
+ "ny": "16",
+ "nz": "4",
+ "particlesPerCell": "40",
+ "seed": "13",
+ "checkpointEvery": "1",
+ "gravity": "9.81",
+ "dtFactor": "0.25"
+}
\ No newline at end of file
More information about the CIG-COMMITS
mailing list