Book a Demo!
CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutPoliciesSign UpSign In
ElmerCSC
GitHub Repository: ElmerCSC/elmerfem
Path: blob/devel/fem/tests/CavityLid2/case.sif
3206 views
! transient cavity lid case with Re=100

$ p0 = 0.0
$ v0 = 1.0
$ visc = 0.01
$ dens = 1.0


Header
  CHECK KEYWORDS Warn
  Mesh DB "." "square"
  Include Path ""
  Results Directory ""
End

Simulation
  Max Output Level = 5
  Coordinate System = Cartesian
  Simulation Type = Transient
  Steady State Max Iterations = 1

  Timestep Intervals = 20
  Timestep Sizes = 0.1
  Output Intervals = 0

  Timestepping method = bdf
  BDF Order = 2

  Post File = case.vtu
  Set Dirichlet BCs by BC Numbering = True
End

Constants
  Gravity(4) = 0 -1 0 9.82
  Stefan Boltzmann = 5.67e-08
  Permittivity of Vacuum = 8.8542e-12
  Boltzmann Constant = 1.3807e-23
  Unit Charge = 1.602e-19
End

Body 1
  Target Bodies(1) = 1
  Name = "Body 1"
  Equation = 1
  Material = 1
  Initial Condition = 1
End


Initial Condition 1 
  Velocity 1 = 1.0e-9
  Velocity 2 = 0.0
  Pressure = $ p0
End 

Solver 1
  Equation = Navier-Stokes

  Stabilization Method=Stabilized

  Steady State Convergence Tolerance = 1.0e-5

  Nonlinear System Convergence Tolerance = 1.0e-8
  Nonlinear System Max Iterations = 20
  Nonlinear System Newton After Iterations = 3
  Nonlinear System Relaxation Factor = 1
  Nonlinear System Newton After Tolerance = 1.0e-3

  Linear System Solver = Iterative
  Linear System Symmetric = True
  Linear System Iterative Method = BicgstabL
  BiCGStabl polynomial degree = 4
  Linear System Max Iterations = 500
  Linear System Convergence Tolerance = 1.0e-10
  Linear System Residual Output = 10
  Linear System Preconditioning = ILU0
End


Equation 1 :: Active Solvers=1

Material 1
  Density = $ dens
  Viscosity = $ visc
End

Boundary Condition 1
  Target Boundaries(3) = 1 2 4
  Name = "NoSlip"
  Velocity 1 = 0.0
  Velocity 2 = 0.0
End	

Boundary Condition 2
  Target Boundaries(1) = 3
  Name = "Top"
  Velocity 1 = $ v0
  Velocity 2 = 0.0
End

! Fix pressure at one node in the center of mesh
Boundary Condition 3 
  Target Coordinates(1,2) = 0.5 0.5
  Name = "FixPressure"
  Pressure = $ p0
End 




Solver 1 :: Reference Norm = Real  2.056714E-01
RUN