underworldcode / underworld2

underworld2: A parallel, particle-in-cell, finite element code for Geodynamics.
http://www.underworldcode.org/
Other
168 stars 58 forks source link

A question about ProjectstressTensor and PressureField #629

Open tyszwh opened 2 years ago

tyszwh commented 2 years ago

Hi,Dear Underworld Team

I encountered some problems in calculating overpressure,

(1) In UWgeo, the output prostress Tensor should be the deviatoric stress,does it right? When reading the output from ParaView xdmf, Projectstress Tensor 0, and Projectstress Tensor 2 corresponds to sigmaxx, sigmayy, that is, the elements on the diagonal (For 2D case). (A previous issue explained this problem). However, It seems that Projectstress Tensor 0 + Projectstress Tensor_ 1 is closer to 0. Why?

(2) Whether the output pressure field is the dynamic pressure (1 / 3(tauxx + tauyy + tauzz)). If that is true, I can use this pressure field to calculate the total stress.

(3) How to quickly calculate the lithostatic pressure from the output results? It is very slow to cycle each depth to calculate the lithostatic pressure at each point of the model