<?xml version="1.0" encoding="UTF-8" ?>
<case app="GenCase4 v4.0.068 (27-04-2018)" date="08-05-2018 11:05:08">
    <casedef>
        <constantsdef>
            <gravity x="0" y="0" z="-9.81" comment="Gravitational acceleration" units_comment="m/s^2" />
            <rhop0 value="1000" comment="Reference density of the fluid" units_comment="kg/m^3" />
            <hswl value="1" auto="false" comment="Maximum still water level to calculate speedofsound using coefsound" units_comment="metres (m)" />
			<!-- <hswl value="0" auto="true" comment="Maximum still water level to calculate speedofsound using coefsound" units_comment="metres (m)" /> -->
            <gamma value="7" comment="Polytropic constant for water used in the state equation" />
            <speedsystem value="0" auto="true" comment="Maximum system speed (by default the dam-break propagation is used)" />
            <coefsound value="20" comment="Coefficient to multiply speedsystem" />
            <speedsound value="0" auto="true" comment="Speed of sound to use in the simulation (by default speedofsound=coefsound*speedsystem)" />
            <coefh value="1.0" comment="Coefficient to calculate the smoothing length (h=coefh*sqrt(3*dp^2) in 3D)" />
            <cflnumber value="0.2" comment="Coefficient to multiply dt" />
        </constantsdef>
        <mkconfig boundcount="240" fluidcount="9" />
        <geometry>
            <definition dp="0.002" units_comment="metres (m)">
                <pointmin x="-2" y="0.2885" z="-1.0" />
                <pointmax x="4" y="0.2885" z="2.5" />
            </definition>
            <commands>
                <mainlist>
                    <setshapemode>actual | bound | fluid</setshapemode>
                    <setdrawmode mode="full" />
                    <setmkbound mk="0" />
                    <drawbeach mask="32"> %|
                        <point x="0" y="0.577" z="1.2" />
                        <point x="1.2" y="0.577" z="0" />
                        <point x="3.2" y="0.577" z="0" />
						<point x="3.2" y="0.577" z="1.2" />
                    </drawbeach>
                    <shapeout file="beach" reset="true" />					
					<setmkbound mk="10" />
					<drawprism>
					<point x="0.1495" y="0" z="1.0525" />
					<point x="0.5731" y="0" z="0.6289" />
					<point x="0.5731" y="0" z="0.7986" />
					<point x="0.2344" y="0" z="1.1373" />
					<point x="0.1495" y="2" z="1.0525" />
					<point x="0.5731" y="2" z="0.6289" />
					<point x="0.5731" y="2" z="0.7986" />
					<point x="0.2344" y="2" z="1.1373" />			
					</drawprism>    
					<setmkbound mk="10" />
                    <fillbox x="0.3613" y="0.2885" z="0.9255">
                        <modefill>void</modefill>
                        <point x="-2" y="-1" z="-1.0" />
                        <size x="6" y="6" z="2.5" />
                    </fillbox>	
                    <shapeout file="block" reset="true" />	
                    <setmkfluid mk="1" />
                    <fillbox x="2" y="0.2885" z="0.12">
                        <modefill>void</modefill>
                        <point x="-2" y="-1" z="-1.0" />
                        <size x="6" y="6" z="1.24" />
                    </fillbox>	
                    <shapeout file="" reset="true" />						
                </mainlist>
            </commands>
        </geometry>
        <floatings>
            <floating mkbound="10"  property="pvc" >
                <massbody value="60.14" /> 	
            </floating>							
        </floatings>
        <properties>
            <propertyfile file="Floating_Materials.xml" path="materials" />
            <property name="userdef01" Restitution_Coefficient_User="0.65" comment="User redefinition for Restitution Coefficient (-)" />
            <links>
                <link mkbound="0" property="steel + userdef01" comment="Property for the tank" />
            </links>
        </properties>
    </casedef>
    <execution>
		<special>
            <chrono>
                <_savedata value="0.01" comment="Saves CSV with data exchange for each time interval (0=all steps)" />
                <schemescale value="1" comment="Scale used to create the initial scheme of Chrono objects (default=1)" />
				<collision activate="true">
				    <distancedp value="0.5" comment="Allowed collision overlap according Dp (default=0.5)" />
					<contactmethod value="0" comment="Contact method type. 0:NSC (Non Smooth Contacts), 1:SMC (SMooth Contacts). (default=0)" />
					<ompthreads value="0" comment="Number of threads by host for parallel execution. 0:Muli Core, 1:Single Core (default=1)" />					
                </collision>
                <bodyfloating id="block" mkbound="10" modelfile="AutoActual" />
                <bodyfixed id="tank" mkbound="0" modelfile="AutoActual" modelnormal="invert" />
            </chrono>
        </special>
        <parameters>
            <parameter key="SavePosDouble" value="1" comment="Saves particle position using double precision (default=0)" />
            <parameter key="StepAlgorithm" value="2" comment="Step Algorithm 1:Verlet, 2:Symplectic (def=1)" />
            <parameter key="VerletSteps" value="40" comment="Verlet only: Number of steps to apply Eulerian equations (def=40)" />
            <parameter key="Kernel" value="2" comment="Interaction Kernel 1:Cubic Spline, 2:Wendland (def=1)" />
            <parameter key="ViscoTreatment" value="1" comment="Viscosity Formulation 1:Artificial, 2:Laminar+SPS (def=1)" />
            <parameter key="Visco" value="0.05" comment="Viscosity value" />
            <parameter key="ViscoBoundFactor" value="1" comment="Multiply viscosity value with boundary (default=1)" />
			<parameter key="DensityDT" value="2" comment="Density Diffusion Term 0:None, 1:Molteni, 2:Fourtakas, 3:Fourtakas(full) (default=0)" />
			<parameter key="DensityDTvalue" value="0.1" comment="DDT value (default=0.1)" />
            <parameter key="Shifting" value="0" comment="Shifting mode 0:None, 1:Ignore bound, 2:Ignore fixed, 3:Full (default=0)" />
            <parameter key="ShiftCoef" value="-2" comment="Coefficient for shifting computation (default=-2)" />
            <parameter key="ShiftTFS" value="0" comment="Threshold to detect free surface. Typically 1.5 for 2D and 2.75 for 3D (default=0)" />
            <parameter key="RigidAlgorithm" value="3" comment="Rigid Algorithm 0:collision-free, 1:SPH, 2:DEM, 3:Chrono (default=1)" />
            <parameter key="FtPause" value="0.0" comment="Time to freeze the floatings at simulation start (warmup) (def=0)" />
            <parameter key="DtIni" value="0" comment="Initial time step. Use 0 to defult use (default=h/speedsound)" units_comment="seconds" />
            <parameter key="DtMin" value="0" comment="Minimum time step. Use 0 to defult use (default=coefdtmin*h/speedsound)" units_comment="seconds" />
            <parameter key="DtAllParticles" value="0" comment="Velocity of particles used to calculate DT. 1:All, 0:Only fluid/floating (default=0)" />
            <parameter key="TimeMax" value="2.0" comment="Time of simulation" />
            <parameter key="TimeOut" value="0.01" comment="Time between output files" />
            <parameter key="PartsOutMax" value="1" comment="%/100 of fluid particles allowed to be excluded from domain (default=1)" units_comment="decimal" />
            <parameter key="RhopOutMin" value="700" comment="Minimum rhop valid (default=700)" units_comment="kg/m^3" />
            <parameter key="RhopOutMax" value="1800" comment="Maximum rhop valid (default=1300)" units_comment="kg/m^3" />
            <simulationdomain comment="Defines domain of simulation (default=Uses minimun and maximum position of the generated particles)">
                <posmin x="default" y="default" z="default" comment="e.g.: x=0.5, y=default-1, z=default-10%" />
                <posmax x="default+0.02" y="default" z="default + 20%" />
            </simulationdomain>
        </parameters>
    </execution>
</case>
