DualSPHysics / DesignSPHysics

Module for FreeCAD intended to use with DualSPHysics fluid simulator.
Other
111 stars 40 forks source link

Please explain why my Inlet is not working. Here is my code. #210

Closed SohailAusgur closed 1 week ago

SohailAusgur commented 5 months ago
    <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.

imestevez commented 1 week ago

Please, post your question into the forum: https://github.com/DualSPHysics/DualSPHysics/discussions