Closed BenedicteCT closed 5 years ago
Hi @BenedicteCT ,
Was the first plot done right after Model inititialization, meaning the pressure is the lithostatic pressure?
After having run model init, the Pressure profile looks like the attached screenshot.
OK that's fine. But what about the first plot you posted. Was it just after the model init or after reading at least one solve?
The two first plots are generated through post processing (python script) of the data recorded on Lagrangian particles after the model has finished running.
@BenedicteCT Q1/P0 elements for Stokes (the default FEM element types in underworld) will produce a discontinuous pressure solutions across elements, so the jumps aren't surprising.
@rbeucher is the Model.pressureField projected? If not, projecting should help, although it'll be less accurate and high gradients could cause oscillations.
@julian you mean on the mesh? I will check but I'm pretty sure it's not.
@Julian you mean on the mesh? I will check but I'm pretty sure it's not.
Yes
Hi Benedicte,
So it's not really an issue but more a consequence of the type of elements we use by default in Underworld. The pressure is calculated on discrete elements Q0 that basically give us a value of the pressure for each cell and the transition from one element to one neighbour is slightly discontinuous (or very if the resolution is poor).
If I understand the plot you sent me you basically have a set of tracers with different depths...
I could use some tricks to make it look smoother but that won't make it more accurate...the temperature is better resolved which explains why your second plot is smoother.
You could increase the resolution and that should make things better.
Now Q0 elements are clearly not ideal, you could use Q1 elements which can resolve the pressure gradient inside elements. I am working on it now but I want to test things properly.
Romain
The Pressure recorded on Lagrangian Particles as (Lag_Grid.add_tracked_field(Model.pressureField, name="LagGridPress", units=u.megapascal, dataType="float") is discontinuous, showing steps in pressure (see screen shot1).
However, plotting against y (screen shot2) yields a smooth curve.
Is there any possibility to record a smoothed instead of discretised pressure in order to get a continuous PT trajectory ?