Hi so i usually used to use ASPECT and i was wonder if i want to define a new rhyology in UW so i use directly the numbers of prefactor, activation volume and energy in rhyology papers and UW will scale the uniaxial strain to plain strain or I have to do the conversion then input? for example UW have built in lower crust of wet anorthite based on Rhybacki 2000 paper where as there is the updated 2006 paper so how can i update the rhyology material?
I guess this is the template: viscosity_dislocation = GEO.ViscousCreep(preExponentialFactor=1.0,
stressExponent=1.0,
activationVolume=0.,
activationEnergy=200 * u.kilojoules,
waterFugacity=0.0,
grainSize=0.0,
meltFraction=0.,
grainSizeExponent=0.,
waterFugacityExponent=0.,
meltFractionFactor=0.0,
f=1.0)
viscosity_diffusion = GEO.ViscousCreep(preExponentialFactor=1.0,
stressExponent=1.0,
activationVolume=0.,
activationEnergy=200 * u.kilojoules,
waterFugacity=0.0,
grainSize=0.0,
meltFraction=0.,
grainSizeExponent=0.,
waterFugacityExponent=0.,
meltFractionFactor=0.0,
f=1.0)
combined_viscosity = GEO.CompositeViscosity([viscosity_dislocation, viscosity_diffusion])
Hi so i usually used to use ASPECT and i was wonder if i want to define a new rhyology in UW so i use directly the numbers of prefactor, activation volume and energy in rhyology papers and UW will scale the uniaxial strain to plain strain or I have to do the conversion then input? for example UW have built in lower crust of wet anorthite based on Rhybacki 2000 paper where as there is the updated 2006 paper so how can i update the rhyology material?
I guess this is the template: viscosity_dislocation = GEO.ViscousCreep(preExponentialFactor=1.0,
stressExponent=1.0,
activationVolume=0.,
activationEnergy=200 * u.kilojoules,
waterFugacity=0.0,
grainSize=0.0,
meltFraction=0.,
grainSizeExponent=0.,
waterFugacityExponent=0.,
meltFractionFactor=0.0,
f=1.0)
viscosity_diffusion = GEO.ViscousCreep(preExponentialFactor=1.0,
stressExponent=1.0,
activationVolume=0.,
activationEnergy=200 * u.kilojoules,
waterFugacity=0.0,
grainSize=0.0,
meltFraction=0.,
grainSizeExponent=0.,
waterFugacityExponent=0.,
meltFractionFactor=0.0,
f=1.0)
combined_viscosity = GEO.CompositeViscosity([viscosity_dislocation, viscosity_diffusion])