How to fill the water in the DPI

Kyushu_panda

Dear all

I am new to DualSPHysics, and I want to create a model describe a wedge slide along a slope into the water, just like the landslide tsunami test in SPHysics. I filled water using a seed point with box limitation just follow the way in CaseFloating, but when I executed the case, it told that

*** Exception (JSphGpuSingle::VisuConfig)
Text: Constant 'b' can not be zero.
'b' is zero when fluid height is zero (or fluid particles were not created)

I am not sure what does "b" means, and I suppose the fluid particles were not created because the GenCase also showed
FileShapes> _Dp.vtk <---- EMPTY FILE NOT SAVED
FileShapes> _Actual.vtk <---- EMPTY FILE NOT SAVED

I create simple river banks like a reservoir using drawprism
setmkbound B="255" G="51" R="51" mk="5" name="Reservior"
drawprism mask="2"
point x="0.5" y="0.0" z="0.0"
point x="1.5" y="0.0" z="0.0"
point x="1.5" y="1.0" z="0.0"
point x="0.5" y="1.0" z="0.0"
point x="0.0" y="0.0" z="0.5"
point x="2.0" y="0.0" z="0.5"
point x="2.0" y="1.0" z="0.5"
point x="0.0" y="1.0" z="0.5"
drawprism>
shapeout file="Reservior"

The fillbox defined in the XML file like this
setmkfluid mk="0" />
fillbox x="0.9" y="0.5" z="0.1"
modefill>void
point x="0" y="0.0" z="0.0"
size x="2" y="1.0" z="0.3"
fillbox
shapeout reset="true" file=""

For display, I have to delete the "<" and "\>"
Anyone can help?
Thanks and best regards!

George

Hi Kyushu_panda or tiger or lion,

Constant b is the c^2rho/gamma from the equation of state,
now you probably have not created any fluid in the domain, and dualsphysics is complaining... quite rightfully so! P=rho*g*h, if your h=0....

I am not going to go through your xml but have a look here:
http://dual.sphysics.org/vanilla/discussion/326/aboutfillvoidpoint,
or the guide please,

also is your gravity 0 in the z direction? if yes then your know what is wrong I suppose!

Cheers tiger, panda or any other animal in this planet
George

Kyushu_panda

Dear George

Thanks for your quick response.
Panda is just a nickname since they look cute.

Now I add a fluid geometry in the model and fill the region with the fillbox.
When calculating, the particles are out of the domain (Reservoir) step by step.

And I visualize the .VTK files in paraview found the boundary particles of the "reservoir" are missing where fluid particles exist. That is why the fluid particle directly move downwards and out of the limitation. I do not know why it occurred and how to fix that.
I have also met the same problem in SPHysics that the fluid particles move out of the boundary.
Would you please give some tips in the model creation.

If I want to make a wedge move down a slope, I should give a motion to the wedge and the wedge can not be set as floating. While I set a uniform rectilinear motion to the wedge, the wedge didn't move at all.

Sorry I have too many questions for you but I really need your help.

Thank you
Best regards!

Wang

Alex

In case of further problems, you can just send your XML to our dualsphysics@gmail.com account, explain the problem and we will help you.

Regards

Alex

Kyushu_panda

Hi Alex

Thanks for your kind help!
Best Regards.

Wang