Pressure calculation in Mphase source code always keep zero

Kyu_MAKyu_MA
in DualSPHysics v5.0

Hi, every one.

I am trying to modify the DualSPHysics code in the multiphase Nnewtonia fluid scheme, and I found the Pressure calculated by the state equation always keen zero:

O%BZ])S@_QUF(8_JHNHGGYO.png


image.png

I have compared it to the Newtonia scheme, and the Pressure can be obtained regularly:

image.png

and also I have checked the state equation code:

image.png

In the NNewtonia scheme, the rhop_r0 always keeps a value of 1, As a result, the pressure is always zero,

I was really confused about this issue cause the momentum and continue equation can not be obtained if this happens.

Comments

  • PrincelxfPrincelxf

    Hello, may I ask which file the code shown in your screenshot is located in? I would like to take a look.

Sign In or Register to comment.

© Vanilla Keystone Theme 2024