Inlet to fill a container with Non-Newtonian fluid

Hi there

I tried to set two inlets to fill a container by DualSPHysics_NonNewtonian.

This is a part of the XML file:

    <special>

    <nnphases> %Defines non-newtonian phases parameters

        <phase mkfluid="0">

          <rhop value="1000" comment="Density of the phase" />           

          <visco value="1e-6" comment="Kinematic viscosity (or consistency index) value for phase (m2/s)" />

          <tau_yield value="0" comment="User defined maximum specific yield stress of phase (Pa m3/kg)" />

          <HBP_m value="0" comment="Use 0 to reduce Newtonian liquid, order of 10 for power law and order of 100 for Bingham (sec) " />

          <HBP_n value="0.95" comment="Use 1 to reduce to Newtonian, <1 for shear thinning >1 for shear thickenning " />           

<phasetype value="0" comment="Non-Newtonian=0 only option in v5.0" />

        </phase>          

      </nnphases>

      <inout resizetime="0.5">

        <determlimit value="1e+3" comment="Use 1e-3 for first_order or 1e+3 for zeroth_order (default=1e+3)" />

        <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" />

          <zone3d comment="Input zone for 3-D simulations">

            <box>

              <point x="0.3" y="2.8" z="0.5" />

              <size x="0.5" y="0" z="0.4" />

              <direction x="0" y="-1" z="0" />

            </box>

          </zone3d>

          <imposevelocity mode="0" comment="Imposed velocity 0:fixed value, 1:variable value, 2:Extrapolated value (default=0)">

            <velocity v="2" comment="Uniform velocity" units_comment="m/s" />

          </imposevelocity>

        </inoutzone>

        <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" />

          <zone3d comment="Input zone for 3-D simulations">

            <circle>

              <point x="2.2" y="2.8" z="0.7" />

              <radius v="0.2" />

              <direction x="0" y="-1" z="0" />

            </circle>

          </zone3d>

          <imposevelocity mode="0" comment="Imposed velocity 0:fixed value, 1:variable value, 2:Extrapolated value (default=0)">

            <velocity v="2" comment="Uniform velocity" units_comment="m/s" />

          </imposevelocity>

        </inoutzone>

      </inout>

    </special>

If I run the case by DualSPHysics_NonNewtonian it faces with the following error:

If I run the case by DualSPHysics_v5.0 or _v4.0 it cannot load the [Multi-phase configuration]

Please advise me in this regard.

Comments

Sign In or Register to comment.