I am trying to add phase transition in a subduction model using visco-plastic material model. I only add phase transition to one of the material (oceanic crust). I started with adding one transition pressure, temperature, and clapeyron slope.
set Define transition by depth instead of pressure = false
set Phase transition pressures= OC: 1.2e9
set Phase transition pressure widths= OC: 0.4e9
set Phase transition temperatures= OC: 803
set Phase transition Clapeyron slopes=OC: 4e6
set Densities = background: 3337.4, OC: 2950|3450, OML: 3337.4, UC: 2811.2, LC: 2929, CML: 3337.4
The transition occurs at the surface of the model, which doesn’t reach the transition pressure or temperature. I guess it is because of the clapeyron slope, or am I misunderstanding it? Is there a way to set a lower pressure/temp limit so that transition doesn’t happen below this limit? Or limit with another clapeyron slope.
I want to include the eclogite phase transition, like in this phase diagram.