Does the mass-density in DualSPHysics code have a limit?

PhoenixJiang

Hi,

Firstly, let me show the problem I was confronted.
In my simulation, the mass-density must be changed from 1000 Kilograms per cubic meter(kg/m3) to 2100 kg/m3.

And there are 2 modifing schemes:
1. Only Change the "Type.h" File
I changed the "Type.h" file, as follows:
change #define RHOPZERO 1000.f to #define RHOPZERO 2100.f
change #define OVERRHOPZERO 0.001f to #define OVERRHOPZERO 4.762e-4f

The source code was compiled successfully.
However, almost 90% particles ran out from the computational domain when I run this program at the first time step.
And the program failed.

2. Change Both the "Type.h" File And the "*_Def.xml" File
I changed the "Type.h" file, as follows:
change #define RHOPZERO 1000.f to #define RHOPZERO 2100.f
change #define OVERRHOPZERO 0.001f to #define OVERRHOPZERO 4.762e-4f
The source code was compiled successfully.
I also changed the "*_Def.xml" file, as follows:
change rhop0 value="1000.0" to rhop0 value="2100.0"

Also almost 90% particles ran out from the computational domain when I run this program at the first time step.
And the program failed.

So I doubt whether there is a limit for mass-density in DualSPHysics code. If there is a limit, what is the limit?
And how can I change the DualSPHysics code in order to run the program successfully in the case that 2100kg/m3(for mass-density)is necessary?


Thanks and regards!

Phoenix Jiang
2013.7.30

PhoenixJiang

Additionally, there is no problem both for the compilation and running when I change the mass-density from 1000kg/m3 to 946kg/m3(less than 1000 kg/m3).