[CIG-LONG] Initial temperature conditions
Robert Gray
graywacke at gmail.com
Tue Feb 21 19:27:05 PST 2012
Sorry,
The attachment doesn't work.
Here's the input file:
{
"EulerDeform":
{
"systems": [
{
"mesh": "v-mesh",
"p-mesh": "p-mesh",
"T-mesh": "T-mesh",
"displacementField": "DisplacementField",
"remesher": "velocityRemesher",
"velocityField": "VelocityField",
"wrapTop": "True",
"staticRight":"True",
"staticRightTop":"True",
"staticLeft":"True",
"staticLeftTop":"True",
"staticBottom":"True"
}
]
},
"components":
{
"buoyancyForceTerm":
{
"Type": "BuoyancyForceTerm",
"ForceVector": "mom_force",
"Swarm": "gaussSwarm",
"gravity": "gravity",
"TemperatureField": "TemperatureField"
},
"backgroundShape":
{
"Type": "Box",
"startX": "0.0e6",
"endX": "2.0e6",
"startY": "0.0e6",
"endY": "0.6e6"
},
"WZShape":
{
"Type": "Box",
"startX": "0.94e6",
"endX": "1.06e6",
"startY": "0.48e6",
"endY": "0.6e6"
},
"MantleShape":
{
"Type": "Intersection",
"shapes":
[
"backgroundShape",
"!WZShape"
]
},
"MantlenonNewtonian":
{
"Type": "NonNewtonian",
"StrainRateInvariantField": "StrainRateInvariantField",
"TemperatureField": "TemperatureField",
"n": "3.5",
"T_0": "6.194370941e4",
"A": "303.85",
"refStrainRate": "1.0e8",
"minViscosity": "1.0e-4",
"maxViscosity": "1.0"
},
"MantlestrainWeakening":
{
"Type": "StrainWeakening",
"TimeIntegrator": "timeIntegrator",
"MaterialPointsSwarm": "materialSwarm",
"initialSofteningStrain": "0.5",
"finalSofteningStrain": "1.5",
"initialDamageFraction": "0.0",
"initialDamageWavenumber": "0.0",
"initialDamageFactor": "0.0",
"healingRate": "0.0"
},
"Mantleyielding":
{
"Type": "DruckerPrager",
"PressureField": "PressureField",
"VelocityGradientsField": "VelocityGradientsField",
"MaterialPointsSwarm": "materialSwarm",
"Context": "context",
"StrainWeakening": "MantlestrainWeakening",
"StrainRateField": "StrainRateField",
"cohesion": "1.0e7",
"cohesionAfterSoftening": "1.0e7",
"frictionCoefficient": "0.267949192",
"frictionCoefficientAfterSoftening": "0.03492",
"minimumViscosity": "1.0e-4"
},
"Mantleviscous":
{
"Type": "RheologyMaterial",
"Shape": "MantleShape",
"density": "3300.0",
"alpha": "2.8e-5",
"diffusivity": "1.0e17",
"heatingElements": [
{
"Q": "0.0",
"lambda": "0.0"
}
],
"Rheology": [
"MantlenonNewtonian",
"Mantleyielding",
"storeViscosity",
"storeStress"
]
},
"WZviscosity":
{
"Type": "MaterialViscosity",
"eta0": "1.0e-4"
},
"WZstrainWeakening":
{
"Type": "StrainWeakening",
"TimeIntegrator": "timeIntegrator",
"MaterialPointsSwarm": "materialSwarm",
"initialSofteningStrain": "0.5",
"finalSofteningStrain": "1.5",
"initialDamageFraction": "0.0",
"initialDamageWavenumber": "0.0",
"initialDamageFactor": "0.0",
"healingRate": "0.0"
},
"WZyielding":
{
"Type": "DruckerPrager",
"PressureField": "PressureField",
"VelocityGradientsField": "VelocityGradientsField",
"MaterialPointsSwarm": "materialSwarm",
"Context": "context",
"StrainWeakening": "WZstrainWeakening",
"StrainRateField": "StrainRateField",
"cohesion": "1.0e7",
"cohesionAfterSoftening": "1.0e7",
"frictionCoefficient": "0.08748864",
"frictionCoefficientAfterSoftening": "0.017455065",
"minimumViscosity": "1.0e-4"
},
"WZviscous":
{
"Type": "RheologyMaterial",
"Shape": "WZShape",
"density": "3300.0",
"alpha": "2.8e-5",
"diffusivity": "1.0e17",
"heatingElements": [
{
"Q": "0.0",
"lambda": "0.0"
}
],
"Rheology": [
"WZviscosity",
"WZyielding",
"storeViscosity",
"storeStress"
]
}
},
"velocityBCs" :{
"type": "CompositeVC",
"vcList": [
{
"type": "WallVC",
"wall": "left",
"variables": [
{
"name": "vx",
"value": "-2.0*6.072087827e12*step(0.47e6-y)"
}
]
},
{
"type": "WallVC",
"wall": "right",
"variables": [
{
"name": "vx",
"value": "-4.756468798e13*step(y-0.48e6)"
}
]
},
{
"type": "WallVC",
"wall": "bottom",
"variables": [
{
"name": "vx",
"value": "0.0"
},
{
"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": "1623.0"
}
]
}
]
},
"temperatureICs":
{
"type": "CompositeVC",
"vcList": [
{
"type": "AllNodesVC",
"variables":[
{
"name": "temperature",
// "value": "(-1350.0/0.6e6)*y + 1623.0"
"value": "(y>=0.48e6) ? -1.125e-2*y + 7023.0 :
1623.0"
}
]
}
]
},
"FieldVariablesToCheckpoint": [
"StrainRateInvariantField",
"VelocityField",
"TemperatureField",
"PressureField"
],
"enable-thermal": true,
"maxTimeSteps": "50",
"outputPath": "./output",
"dim": "2",
"minX": "0",
"minY": "0",
"minZ": "0",
"maxX": "2.0e6",
"maxY": "0.6e6",
"maxZ": "0.3e6",
"nx": "50",
"ny": "10",
"nz": "4",
"particlesPerCell": "40",
"seed": "13",
"checkpointEvery": "10",
"gravity": "9.8"
}
Rob
On Tue, Feb 21, 2012 at 1:36 PM, Robert Gray <graywacke at gmail.com> wrote:
> Hey all,
> Here's the input file that deals with my previous inquiry.
> Cheers,
> Rob
>
>
>
> On Tue, Feb 21, 2012 at 1:16 PM, Walter Landry <walter at geodynamics.org>wrote:
>
>> Robert Gray <graywacke at gmail.com> wrote:
>> > Hello everyone,
>> > I'm trying to set up an initial temperature profile in Gale that
>> increases
>> > steadily from 273 K to 1623 K in the top 120 km of the box, while
>> > everything below that is 1623 K. To do this I've used the following
>> > condition:
>> >
>> > "value": "(y>=0.48e6) ? -1.125e-2*y + 7023.0 : 1623.0"
>> >
>> > It works great for the first time step, but after that temperatures jump
>> > several orders of magnitude to unrealistic values. That said, the
>> > incremental increase :
>> >
>> > "value": "(-1350.0/0.6e6)*y + 1623.0"
>> >
>> > works fine but does not produce the initial temperature profile I need.
>>
>> Can you send the input file to the mailing list?
>>
>> Cheers,
>> Walter Landry
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-long/attachments/20120221/bc63d576/attachment-0001.htm
More information about the CIG-LONG
mailing list