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

  • That is error is not related to DualSPHysics....

    GenCase is the problem since it can not find the input XML file


    Regards

  • Hello there

    I tried to correct my xml file to fill a container by one inlet using DualSPHysics_NonNewtonian.

    I reached to the message:

    Please advise me what course of action should be done to have inlet/outlet option in multiphase formulation.


    ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

        <geometry>

          <definition dp="0.05" units_comment="metres (m)">

            <pointmin x="-2" y="-2" z="-1" />

            <pointmax x="4" y="4" z="5" />

          </definition>

          <commands>

            <mainlist>

              <setshapemode>dp | bound</setshapemode>

              <setdrawmode mode="full" />

              <setmkbound mk="0" />

              <drawbox>

                <boxfill>all^top</boxfill>

                <point x="-1.5" y="-1.5" z="0" />

                <size x="3" y="3" z="2" />

              </drawbox>

              <setmkvoid />

              <drawbox>

                <boxfill>bottom</boxfill>

                <point x="-0.4" y="-0.4" z="0" />

                <size x="0.8" y="0.8" z="1" />

              </drawbox>

              <setmkfluid mk="4" />

              <drawbox>

                <boxfill>solid</boxfill>

                <point x="-1.5" y="-1.4" z="0.2" />

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

              </drawbox>

              <shapeout file="" />

            </mainlist>

          </commands>

        </geometry>

      </casedef>

      <execution>

        <special>

          <inout>

    <memoryresize size0="10" size="10" comment="Initial memory resize (size0) and the following memory resizes according to initial inlet/outlet particles (default=2 and 4)"/>

            <useboxlimit value="true" comment="In/out process is only applied to InOut zones delimited by BoxLimit (default=true)">

              <freecentre x="0" y="0" z="0.3" comment="Centre of zone where InOut is not applied (default=centre of simulation domain)" units_comment="m" />

            </useboxlimit>

            <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="2" comment="Treatment of fluid entering the zone. 0:No changes, 1:Convert fluid (necessary for outlet), 2:Remove fluid (recommended for inlet)" />

              <layers value="6" comment="Number of inlet/outlet particle layers" />

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

                <particles mkfluid="4" direction="right" />

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

    <nnphases> %Defines non-newtonian phases parameters

            <phase mkfluid="4">

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

              <visco value="0.9" comment="Kinematic viscosity (or consistency index) value for phase (m2/s)" />

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

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

              <HBP_n value="5" 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>

        </special>

        <parameters>

          <parameter key="SavePosDouble" value="0" comment="Saves particle position using double precision (default=0)" />

          <parameter key="StepAlgorithm" value="2" comment="Step Algorithm 1:Verlet, 2:Symplectic (default=1)" />

          <parameter key="VerletSteps" value="40" comment="Verlet only: Number of steps to apply Euler timestepping (default=40)" />

          <parameter key="Kernel" value="2" comment="Interaction Kernel 1:Cubic Spline, 2:Wendland (default=2)" />

    %Choice of reology treatment, velocity gradient calculation and viscosity treatment  

    <parameter key="RheologyTreatment" value="2" comment="Reology formulation 1:Single-phase classic, 2: Single and multi-phase" />       

    <parameter key="VelocityGradientType" value="1" comment="Velocity gradient formulation 1:FDA, 2:SPH" />

          <parameter key="ViscoTreatment" value="2" comment="Viscosity formulation 1:Artificial, 2:Laminar+SPS, 3:Constitutive eq." />     

    %Wall boundary viscosity or/and artificial viscosity if ViscoTreatment is 1:Artificial 

    <parameter key="Visco" value="0.05" comment="Viscosity value" /> % Note alpha can depend on the resolution when using artificial viscosity

          <parameter key="ViscoBoundFactor" value="1" comment="Multiply viscosity value with boundary (default=1)" />

  • I think that the message is clear enough...

    Inlet and outlet conditions are not implemented yet to be used with multiphase...

    I do not understand your question

  • Thanks for your response.

    Is it possible to modify NNewtonian_DualSPHysics source files using inlet-outlet source files in classic DualSPHysicsv_5.0 to implement Inlet-outlet option in NNewtonian_DualSPHysics?

  • We have not carried out that before, but it is planned for the next release (v5.2).

    However since you have already all source files you can give a try on your own before.

    Regards

  • Hope to see that soon.

    Thanks

Sign In or Register to comment.