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.67e08
Permittivity of Vacuum = 8.8542e12
Boltzmann Constant = 1.3807e23
Unit Charge = 1.602e19
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.0e5
Nonlinear System Convergence Tolerance = 1.0e7
Nonlinear System Max Iterations = 1
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 1.0e3
Nonlinear System Relaxation Factor = 0.5
Linear System Solver = Direct
Linear System Direct Method = Umfpack
End
Solver 1
Equation = NavierStokes
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.0e4
Nonlinear System Convergence Tolerance = 1.0e7
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 = KEpsilon
Procedure = "KESolver" "KESolver"
Exec Solver = Always
Stabilize = True
Bubbles = False
Lumped Mass Matrix = False
Optimize Bandwidth = True
Steady State Convergence Tolerance = 1.0e5
Nonlinear System Convergence Tolerance = 1.0e7
Nonlinear System Max Iterations = 20
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 1.0e3
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 = KEpsilon
End
Material 3
Name = "Material 3"
Viscosity Model = KEpsilon
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.0e4
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
Elmer+Paraview

 Posts: 24
 Joined: 28 Jul 2020, 11:24
 Antispam: Yes
Elmer+Paraview
 Attachments

 Mesh_7.zip
 (789.17 KiB) Downloaded 19 times
Re: Elmer+Paraview
Dear armaghantelgerafchi ,
1. Paraview is working well as you have put the following code in Simulation section.
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
1. Paraview is working well as you have put the following code in Simulation section.
Code: Select all
Simulation
...
Post File = case.vtu
End
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
Department of Materials Engineering, KU Leuven, Belgium