!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!Parameters !!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! $ bw=0.03 ! 1e-5 ! 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.5 ! wall velocity $ g= 9.81 $ h = 1 $ Up1=0.05 $ h1 = (nuPlus*Up/g)^0.5 ! [m] $ sigma= 1.6e7 $ epsilon= 1e-3 $ y1= epsilon/h1 $ B= 0.5 ! tesla $ Ha= B*sqrt(sigma*Up/(rhoPlus*g)) $ relax = 1e-2 $ a = 1.25e-4 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Header Mesh DB "." "." 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.01 Timestep Intervals = Integer 10000 Timestep Function Real Procedure "LevelSet" "LevelSetTimestep" LevelSet Courant Number = Real 0.5 Steady State Max Iterations = Integer 1 ! 1 Output Intervals = Integer 10 ! Coordinate Scaling = 0.01 ! Post File = File "data.ep" Solver Input File = case.sif Post 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 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.5*g*tx^2)/nuPlus - h*h1*tx*g/nuPlus +Up1; else 0;" End Equation 1 Active Solvers(4) = 1 2 3 4 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 3 Equation = Navier-Stokes Variable = Flow Solution[Velocity:2 Pressure:1] Procedure = "FlowSolve" "FlowSolver" Stabilize = True Optimize Bandwidth = True Steady State Convergence Tolerance = 1.0e-5 Nonlinear System Convergence Tolerance = 1.0e-4 Nonlinear System Max Iterations = 1 Nonlinear System Newton After Iterations = 100 Nonlinear System Newton After Tolerance = 0.0 Nonlinear System Relaxation Factor = $relax Linear System Solver = direct Linear System Direct Method = umfpack 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 Nonlinear System Newton After Iterations = 100 Nonlinear System Newton After Tolerance = 0.0 Linear System Solver = direct Linear System Direct Method = umfpack 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 4 Exec Solver = always Procedure = File "SaveData" "SaveMaterials" Parameter 1 = String Density Parameter 2 = Viscosity End Boundary Condition 1 Target Boundaries (2) = 3 5 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 = $ Up1 ! [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.5*(tx/h1)^2+ h*(tx/h1) -0.5; 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.5*(tx/h1)^2+ h*(tx/h1) -0.5; else 0;" End ! Solver 1 :: Reference Norm = 5.29208635E-07 ! Solver 2 :: Reference Norm = 6.00803071E-01