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
Elmer+Paraview
-
- Posts: 26
- Joined: 28 Jul 2020, 11:24
- Antispam: Yes
Elmer+Paraview
- Attachments
-
- Mesh_7.zip
- (789.17 KiB) Downloaded 220 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
Faculty of Mechanical Engineering, Silesian University of Technology, Gliwice
Faculty of Mechanical Engineering, Silesian University of Technology, Gliwice