No particles with mkfluid 1 (fluid with a prism shape)

Shadi

Dear DualSPHysics users,

I am trying to set a model with a non-Newtonian fluid on a ramp which its initial shape is a prism.

First I used drawprism and then fill the prism as following:

 <setmkbound mk="1" />

          <drawprism >

            <point x="50" y="12" z="170" />

            <point x="75" y="12" z="125" />

            <point x="100" y="12" z="170" />

            <point x="50" y="-12" z="170" />

            <point x="75" y="-12" z="125" />

            <point x="100" y="-12" z="170" />

             

          </drawprism>

           <shapeout file="Mud" /> 

 <!-- Mud-->   

           

           <setmkfluid mk="1" />

          <fillbox x="75" y="0" z="130">

            <modefill>void</modefill>

            <point x="-500" y="-200" z="125" />

            <size x="400" y="600" z="25" />

          </fillbox>

          <_pointsmkout file="fluid2" />

          <shapeout file="Mud" />       

           

But I faced with the error that there is no particles with mkfluid=1

I would be glad if you could help me to solve it.

Thank you.

TPouzol

Hello,

Your fillbox has an error in X.

point -500, size 400 would give you a max x at -100

and your seed is 75 so no particles generated

Regards

Shadi

https://forums.dual.sphysics.org/discussion/comment/5839#Comment_5839

Hello TPouzol,

Thanks for your response.

I changed it to:

<setmkbound mk="1" />

     <drawprism >

      <point x="50" y="12" z="170" />

      <point x="75" y="12" z="125" />

      <point x="100" y="12" z="170" />

      <point x="50" y="-12" z="170" />

      <point x="75" y="-12" z="125" />

      <point x="100" y="-12" z="170" />

        </drawprism>

  <shapeout file="Mud" /> 

 <setmkfluid mk="1" />

     <fillbox x="75" y="0" z="130">

      <modefill>void</modefill>

      <point x="-500" y="-200" z="125" />

      <size x="700" y="600" z="25" />

     </fillbox>

     <_pointsmkout file="fluid2" />

     <shapeout file="Mud" />   

But I faced with the same error.

TPouzol

Hello,

Prism should be in mkfluid also, and according to the prism and your drawing the fillbox is weird. Something like :

<setmkfluid mk="1" />

<drawprism >

      <point x="50" y="12" z="170" />

      <point x="75" y="12" z="125" />

      <point x="100" y="12" z="170" />

      <point x="50" y="-12" z="170" />

      <point x="75" y="-12" z="125" />

      <point x="100" y="-12" z="170" />

        </drawprism>

<fillbox x="75" y="0" z="130">

      <modefill>void</modefill>

      <point x="50" y="-12" z="125" />

      <size x="50" y="24" z="45" />

     </fillbox>

     <_pointsmkout file="fluid2" />

     <shapeout file="Mud" />

Also what is the resolution ?

And I think there is an option to autofill the prism (check documentation). I never use draw function but only import stl files can't help you here

Regards

Shadi

https://forums.dual.sphysics.org/discussion/comment/5842#Comment_5842

Hello,

Thank you for your suggestion. I could solve it and fluid particles could be stored. The resolution is dp=0.8 .

However, some particles penetrate from the boundary.

I changed the StepAlgorithm from Symplectic to verlet and still the particles peneterate throgh the boundary.

I am seeking any ideas or suggestions you may have to help resolve this matter.

Regards,

TPouzol

Hello,

Nice !

You should ensure that your boundaries are multi layered !

Regards

Shadi

https://forums.dual.sphysics.org/discussion/comment/5845#Comment_5845

Hello,

Thanks for your suggestion :)

I am trying to apply try mDBC.

Best regards,