Elmer+Paraview

The graphical user interface of Elmer
Post Reply
armaghantelgerafchi
Posts: 24
Joined: 28 Jul 2020, 11:24
Antispam: Yes

Elmer+Paraview

Post by armaghantelgerafchi »

Hello
I have a problem with Paraview. I used the Elmer program for modeling the geometry of a cylinder that has an empty tube inside that. It seems that In Paraview , the velocity distribution inside the empty tube doesn't work. I expect to see higher velocity inside the tube that the bottom and top of tube. But, it seems that the elmer doesn't calculate velocity for empty tube inside cylinder. I would be pleased if you help me with this issue.
I attached mesh file and SIF file for your references.

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

Simulation
Max Output Level = 5
Coordinate System = Cartesian
Coordinate Mapping(3) = 1 2 3
Simulation Type = Steady state
Steady State Max Iterations = 1
Output Intervals = 1
Timestepping Method = BDF
BDF Order = 1
Solver Input File = case.sif
Post File = case.vtu
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

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

Body 3
Target Bodies(1) = 3
Name = "Body 3"
Equation = 2
Material = 3
Initial condition = 1
End

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

Body 5
Target Bodies(1) = 5
Name = "Body 5"
Equation = 2
Material = 3
Initial condition = 1
End

Body 6
Target Bodies(1) = 6
Name = "Body 6"
Equation = 2
Material = 3
Initial condition = 1
End

Solver 3
Equation = Heat Equation
Procedure = "HeatSolve" "HeatSolver"
Variable = Temperature
Exec Solver = Always
Stabilize = True
Bubbles = False
Lumped Mass Matrix = False
Optimize Bandwidth = True
Steady State Convergence Tolerance = 1.0e-5
Nonlinear System Convergence Tolerance = 1.0e-7
Nonlinear System Max Iterations = 1
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Relaxation Factor = 0.5
Linear System Solver = Direct
Linear System Direct Method = Umfpack
End

Solver 1
Equation = Navier-Stokes
Variable = Flow Solution[Velocity:3 Pressure:1]
Procedure = "FlowSolve" "FlowSolver"
Exec Solver = Always
Stabilize = True
Bubbles = False
Lumped Mass Matrix = False
Optimize Bandwidth = True
Steady State Convergence Tolerance = 1.0e-4
Nonlinear System Convergence Tolerance = 1.0e-7
Nonlinear System Max Iterations = 1
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 0.0
Nonlinear System Relaxation Factor = 0.5
Linear System Solver = Direct
Linear System Direct Method = Umfpack
End

Solver 2
Equation = K-Epsilon
Procedure = "KESolver" "KESolver"
Exec Solver = Always
Stabilize = True
Bubbles = False
Lumped Mass Matrix = False
Optimize Bandwidth = True
Steady State Convergence Tolerance = 1.0e-5
Nonlinear System Convergence Tolerance = 1.0e-7
Nonlinear System Max Iterations = 20
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Relaxation Factor = 1
Linear System Solver = Direct
Linear System Direct Method = Umfpack
End

Equation 1
Name = "Equation 1"
Convection = Computed
Active Solvers(1) = 3
End

Equation 2
Name = "Equation 2"
Convection = Computed
Active Solvers(3) = 3 1 2
End

Material 1
Name = "Material 1"
Heat Conductivity = 40
Density = 8000
End

Material 2
Name = "Material 2"
Density = 1800
Viscosity = 0.0009
Heat Conductivity = 88
Viscosity Model = K-Epsilon
End

Material 3
Name = "Material 3"
Viscosity Model = K-Epsilon
Density = 0.05
Viscosity = 0.0008
Heat Conductivity = 85
End

Initial Condition 1
Name = "InitialCondition 1"
Velocity 2 = 0
Kinetic Energy = 0.00457
Temperature = 1073
Velocity 1 = 0
Velocity 3 = 0
Kinetic Dissipation = 1.0e-4
End

Boundary Condition 1
Target Boundaries(1) = 7
Name = "BoundaryCondition 1"
Velocity 1 = 0
Velocity 2 = 0
Velocity 3 = 0
Temperature = 1073
End

Boundary Condition 2
Target Boundaries(1) = 10
Name = "BoundaryCondition 2"
Temperature = 1073
External Pressure = Variable Temperature
Real MATC " exp(-7750/(tx) + 12.79) + exp(-1.41)+tx+5400"
Pressure 1 = Variable Temperature
Real MATC " exp(-7750/(tx) + 12.79) + exp(-1.41)+tx+5400"
Pressure 2 = Variable Temperature
Real MATC " exp(-7750/(tx) + 12.79) + exp(-1.41)+tx+5400"
End

Boundary Condition 3
Target Boundaries(1) = 13
Name = "BoundaryCondition 3"
Pressure 2 = Variable Temperature
Real MATC " exp(-7750/(tx) + 12.79) + exp(-1.41)+tx"
Pressure 1 = Variable Temperature
Real MATC " exp(-7750/(tx) + 12.79) + exp(-1.41)+tx"
Temperature = 1073
Pressure 3 = Variable Temperature
Real MATC " exp(-7750/(tx) + 12.79) + exp(-1.41)+tx"
End

Boundary Condition 4
Target Boundaries(12) = 14 16 18 20 22 23 24 27 28 29 30 31
Name = "BoundaryCondition 4"
Noslip wall BC = True
Wall Law = True
Boundary Layer Thickness = 0.1
End

Boundary Condition 5
Target Boundaries(1) = 32
Name = "BoundaryCondition 5"
Temperature = 1072
Pressure 1 = Variable Temperature
Real MATC " exp(-7750/(tx) + 12.79) + exp(-1.41)+tx"
Pressure 2 = Variable Temperature
Real MATC " exp(-7750/(tx) + 12.79) + exp(-1.41)+tx"
Pressure 3 = Variable Temperature
Real MATC " exp(-7750/(tx) + 12.79) + exp(-1.41)+tx"
End
Attachments
Mesh_7.zip
(789.17 KiB) Downloaded 19 times
annier
Posts: 1141
Joined: 27 Aug 2013, 13:51
Antispam: Yes

Re: Elmer+Paraview

Post by annier »

Dear armaghantelgerafchi ,
1. Paraview is working well as you have put the following code in Simulation section.

Code: Select all

Simulation
...
Post File = case.vtu
End
2. You can try changing the boundary conditions that fit the actual physical conditions.
The boundary conditions pressure 1, pressure 2 and pressure 3 are the forces at the boundary.
The absolute pressure at the boundary is just mentioned as pressure. or, the normal pressure to the surface is referred to as External pressure.
Reference: Page 30 Elmer Models Manual (http://www.nic.funet.fi/pub/sci/physics ... Manual.pdf)


Yours Sincerely,
Anil Kunwar
Anil Kunwar
Department of Materials Engineering, KU Leuven, Belgium
Post Reply