!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!Parameters !!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! $ bw=0.03 ! delta coeficient of level set $ E=1e-4 ! smooth coeficient originariamente E=0.05 $ rhoMin=1.0 ! air density kg/m^3 $ rhoPlus = 6500.0 ! zinc density kg/m^3 $ muMin = 0.00001 ! air dynamic viscosity [Pa/s] $ muPlus = 0.0029 ! zinc dynamic viscosity [Pa/s] $ nuPlus = muPlus/rhoPlus $ Up= 0.005 ! wall velocity $ g= 9.81 $ h = 0.5 ! $ h1 = (nuPlus*Up/g)^0.5 ! [m] !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Header Mesh DB "." "mesh" Include Path "" Results Directory "" End Simulation Max Output Level = 5 ! For verbosity Coordinate system = "Cartesian 2D" Coordinate mapping(3) = 1 2 3 Simulation Type = Transient ! Steady state Timestepping Method = String BDF BDF Order = Integer 2 Timestep sizes = Real 0.0001 Timestep Intervals = Integer 2 Timestep Function Real Procedure "LevelSet" "LevelSetTimestep" LevelSet Courant Number = Real 0.5 Steady State Max Iterations = Integer 1 ! 1 Output Intervals = Integer 1 ! Coordinate Scaling = 0.01 ! Output File = File "data.dat" Post File = File "case.vtu" End Body 1 Name = "Vessel" Equation = 1 Material = 1 Body Force = 1 Initial Condition = 1 End Body Force 1 Name = "BodyForce 1" Flow Bodyforce 2 = -9.816 End Initial Condition 1 ! Pressure = 0.0 Surface = Variable Coordinate 1 Real MATC "-tx+h*h1" ! preaviously 0.0575originaly 0.2 Velocity 1 = 0 Velocity 2 = Variable Coordinate 1 Real MATC "if (tx<=h*h1) -0.005*(tx/h1)^2+ h*(tx/h1) -0.005; else 0;" End ! -0.5*g*tx^2/nuPlus + g*h*tx/nuPlus -Up Equation 1 Active Solvers(5) = 1 2 3 4 5 End Material 1 ! Air has negative values of LS marker function Name = "Zinc-air" Density = Variable Surface Real MATC "(rhoPlus-rhoMin)*(0.5*(1+tanh(3.14*tx/E)))+rhoMin" Viscosity = Variable Surface Real MATC "(muPlus-muMin)*(0.5*(1+tanh(3.14*tx/E)))+muMin" LevelSet Bandwidth = Real $ bw Levelset Velocity 1 = Equals Velocity 1 Levelset Velocity 2 = Equals Velocity 2 End Solver 1 Equation = String "Level Set" Variable = String Surface Variable DOFs = Integer 1 Procedure = File "LevelSet" "LevelSetSolver" Steady State Convergence Tolerance = Real 1e-6 Stabilize = Logical True Nonlinear System Relaxation Factor = Real 1.0 Nonlinear System Max Iterations = Integer 1 Nonlinear System Convergence Tolerance = Real 1.0e-6 Linear System Solver = Iterative ! direct Linear System Iterative Method = BiCGstab ! umfpack Linear System Max Iterations = 500 Linear System Convergence Tolerance = 1.0e-8 BiCGstabl polynomial degree = 2 Linear System Preconditioning = ILU0 End Solver 2 Equation = String "Reinitialize" Variable = String Renorm Variable DOFs = Integer 1 Procedure = File "LevelSet" "LevelSetDistance" Extract Interval = Integer 2 Reinitialize Interval = Integer 10 End Solver 3 Equation = "Navier-Stokes" Linear System Solver = Iterative ! direct Linear System Iterative Method = BiCGstab ! umfpack Linear System Max Iterations = 800 Linear System Convergence Tolerance = 1.0e-5 BiCGstabl polynomial degree = 2 Linear System Preconditioning = ILU1 Stabilize = Logical True Nonlinear System Relaxation Factor = Real 0.5 Nonlinear System Max Iterations = Integer 10 Nonlinear System Convergence Tolerance = Real 1.0e-3 Nonlinear System Newton After Iterations = 10 Nonlinear System Newton After Tolerance = 0 End Solver 4 Exec Solver = always Procedure = File "SaveData" "SaveMaterials" Parameter 1 = String Density Parameter 2 = Viscosity End Solver 5 Equation = "result output" Procedure = "ResultOutputSolve" "ResultOutputSolver" No Fileindex = Logical True Binary Output = Logical True Single Precision = Logical True Vtu Format = Logical True End Boundary Condition 1 Target Boundaries = 3 Name = "wall" ! Velocity 1 = 0 ! Noslip wall BC = True ! Periodic BC = 1 ! Normal-Tangential Velocity = True Velocity 1 = 0 Velocity 2 = 0 End Boundary Condition 2 name = "moving wall" target boundaries = 1 ! Normal-Tangential Velocity = True Velocity 1 = 0 Velocity 2 = $ Up ! [m/s] End Boundary Condition 3 name = "Inlet" target boundaries = 4 ! Normal-Tangential Velocity = True Velocity 1 = 0 Velocity 2 = Variable Coordinate 1 Real MATC "if (tx<=h*h1) -0.005*(tx/h1)^2+ h*(tx/h1) -0.005; else 0;" ! Periodic BC = 4 !è il numero della boundary condition dove applica la periodicità ovvero la voglio anche al "top" ! Periodic BC Velocity 1 = Logical True ! Periodic BC Velocity 2 = Logical True End Boundary Condition 4 name = "Outlet" target boundaries = 2 ! Pressure 1 = 0 ! Normal-Tangential Velocity = True Velocity 1 = 0 Velocity 2 = Variable Coordinate 1 Real MATC "if (tx<=h*h1) -0.005*(tx/h1)^2+ h*(tx/h1) -0.005; else 0;" End