Elmer+Paraview
Posted: 28 Aug 2020, 13:27
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
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