DesignSPHysics
DesignSPHysics copied to clipboard
Please explain why my Inlet is not working. Here is my code.
trafficstars
<geometry>
<definition dp="0.003" comment="Initial inter-particle distance" units_comment="metres (m)">
<pointref x="0.0" y="0.0" z="0.0" />
<pointmin x="0.0" y="0.15" z="0.0" />
<pointmax x="0.5" y="0.15" z="0.5" />
</definition>
<commands>
<mainlist>
<setshapemode>actual | dp | bound</setshapemode>
<move x="0.135" y="0.135" z="0.3" />
<setmkfluid mk="0"/>
<fillbox x="0.015" y="0.015" z="0.01" objname="WaterFlow">
<modefill>void</modefill>
<point x="0" y="0" z="0" />
<size x="0.03" y="0.03" z="0.03" />
</fillbox>
<matrixreset />
<setmkbound mk="0"/>
<setdrawmode mode="face"/>
<drawbox objname="Domain">
<boxfill>front | back | bottom | left | right</boxfill>
<point x="0.0" y="0.0" z="0.0" />
<size x="0.3" y="0.3" z="0.42" />
</drawbox>
<setmkfluid mk="1"/>
<fillbox x="0.15" y="0.15" z="0.1" objname="Water">
<modefill>void</modefill>
<point x="0" y="0" z="0" />
<size x="0.3" y="0.3" z="0.2" />
</fillbox>
<setmkfluid mk="2"/>
<setdrawmode mode="face"/>
<drawbox objname="Inlet">
<boxfill>top</boxfill>
<point x="0.135" y="0.135" z="0.3" />
<size x="0.03" y="0.03" z="0.02" />
</drawbox>
<shapeout file="" />
</mainlist>
</commands>
</geometry>
</casedef>
<execution>
<special>
<initialize>
<fluidvelocity mkfluid="2">
<direction x="0.0" y="0.0" z="-2.0" />
<velocity v="-2.0" comment="Uniform profile velocity" units_comment="m/s" />
</fluidvelocity>
</initialize>
<inout>
<memoryresize size0="2.0" size="4.0" comment="Initial memory resize (size0) and the following memory resizes according to initial inlet/outlet particles (default=2 and 4)"/>
<determlimit value="1e+3" comment="Use 1e-3 for first_order or 1e+3 for zeroth_order (default=1e+3)" />
<extrapolatemode value="1" comment="Calculation mode for rhop and velocity extrapolation from ghost nodes 1:fast-single, 2:single, 3:double (default=1)" />
<inoutzone>
<refilling value="0" comment="Refilling mode. 0:Simple full, 1:Simple below zsurf, 2:Advanced for reverse flows (very slow) (default=1)" />
<inputtreatment value="0" comment="Treatment of fluid entering the zone. 0:No changes, 1:Convert fluid (necessary for outlet), 2:Remove fluid (recommended for inlet)" />
<layers value="4" comment="Number of inlet/outlet particle layers" />
<zone2d comment="Input zone for 2D simulations">
<particles mkfluid="2" direction="bottom" />
</zone2d>
<imposevelocity mode="0" comment="Imposed velocity 0:fixed value, 1:variable value, 2:Extrapolated value (default=0)">
<velocity v="2.0" comment="Uniform velocity" units_comment="m/s" />
</imposevelocity>
<imposerhop mode="0" comment="Outlet rhop 0:Imposed fixed value, 1:Hydrostatic, 2:Extrapolated from ghost nodes (default=0)" />
<imposezsurf mode="0" comment="Inlet Z-surface 0:Imposed fixed value, 1:Imposed variable value, 2:Calculated from fluid domain (default=0)">
<zbottom value="0.0" comment="Bottom level of water (used for Hydrostatic option)" units_comment="m" />
<zsurf value="0.2" comment="Characteristic inlet Z-surface (used for Hydrostatic option)" units_comment="m" />
</imposezsurf>
</inoutzone>
</inout>
</special>
I did this in DesignSPH. I cannot seem to get the inlet to work to create a impinging water jet into the free surface. Many Thanks for your help.