Ranque-Hilsch Vortex Tube Modelling on ELMER

Numerical methods and mathematical models of Elmer
Post Reply
serge_13
Posts: 24
Joined: 16 Oct 2015, 17:15
Antispam: Yes

Ranque-Hilsch Vortex Tube Modelling on ELMER

Post by serge_13 »

Dear friends!

I'd like to present you a new paper on Elmer's modelling:
"Model of anisotropic turbulence in finite-element package Elmer in calculation of a vortex flow parameters and thermal separation in the Ranque-Hilsch vortex tube"
It is saved here:
http://www.t-flow.ru/images/pdf_article ... _b_eng.pdf

Best regards !
Serge.
annier
Posts: 1168
Joined: 27 Aug 2013, 13:51
Antispam: Yes

Re: Ranque-Hilsch Vortex Tube Modelling on ELMER

Post by annier »

Serge,
Excellent work on modelling by you using elmer . Keep it up.

Yours Sincerely,
Anil Kunwar
Anil Kunwar
Faculty of Mechanical Engineering, Silesian University of Technology, Gliwice
serge_13
Posts: 24
Joined: 16 Oct 2015, 17:15
Antispam: Yes

Re: Ranque-Hilsch Vortex Tube Modelling on ELMER

Post by serge_13 »

Thank You !
serge_13
Posts: 24
Joined: 16 Oct 2015, 17:15
Antispam: Yes

Re: Ranque-Hilsch Vortex Tube Modelling on ELMER

Post by serge_13 »

I'm sorry to all.
There some mistakes are in the paper.
Of course, k-turb and mu-turb values in the system of PDEs (1-5) are the tensors. Main diagonals are not empty. For mu-turb diagonal elements are: mu_turb_xx, mu_turb_yy, mu_turb_zz. For k-turb the diagonal elements are: k_turb_xx, k_turb_yy, k_turb_zz.

The texts of paper will be updated soon.

Best regards,
Serge
annier
Posts: 1168
Joined: 27 Aug 2013, 13:51
Antispam: Yes

Re: Ranque-Hilsch Vortex Tube Modelling on ELMER

Post by annier »

Hi Serge,
That's excellent that this forum has been a great centre for learning and reviewing the updated modeling results.

Yours Sincerely,
Anil Kunwar
Anil Kunwar
Faculty of Mechanical Engineering, Silesian University of Technology, Gliwice
serge_13
Posts: 24
Joined: 16 Oct 2015, 17:15
Antispam: Yes

Re: Ranque-Hilsch Vortex Tube Modelling on ELMER

Post by serge_13 »

Good day to All!

I prepared final paper of Ranque-Hilsch tube investigation by Elmer's software. Here is the title:
"CFD-study of the adiabatic deceleration in the role of a heat source for the effect of temperature separation in Ranque-Hilsch vortex tube". I think the method will be helpful in your investigations.

You may download it from the URL:
http://www.t-flow.ru/images/pdf_article ... on_eng.pdf

Wbr.
Serge.
annier
Posts: 1168
Joined: 27 Aug 2013, 13:51
Antispam: Yes

Re: Ranque-Hilsch Vortex Tube Modelling on ELMER

Post by annier »

Hi Serge,
Great update.

Yours Sincerely,
Anil Kunwar
Anil Kunwar
Faculty of Mechanical Engineering, Silesian University of Technology, Gliwice
serge_13
Posts: 24
Joined: 16 Oct 2015, 17:15
Antispam: Yes

Re: Ranque-Hilsch Vortex Tube Modelling on ELMER

Post by serge_13 »

Dear friends !

You may to model Ranque-Hilsh Tube, Vortex Chambers and ... yourself.
The Elmer's modification with Anisotropic model of the turbulence is prepared by me. See the following reference:
“T-Flow Lab” Modification of Elmer FEM Package for Vortex and Complex Flows Modeling"
http://www.t-flow.ru/images/pdf_article ... ab_eng.pdf

Best regards,
Serge.
serge_13
Posts: 24
Joined: 16 Oct 2015, 17:15
Antispam: Yes

Re: Ranque-Hilsch Vortex Tube Modelling on ELMER

Post by serge_13 »

Dear friends. All old links are available again at https://www.t-flow.ru (remastered after hacking destroy) ...

Sorry, the ssl certificate will be updated soon ...

New 'CGNSSaveIterative' solver module was created to save transient data for Post Processing via cgns format. You can see temperature segregation (transient) process in Ranque-Hilsch tube. Starting conditions are: - constant temperature; - quasi-solid flow rotation in the tube. Time axis: - 31 step ; - start : 0.000245 s; - end : 0.001726 s.

Best regards,
Serge.
serge_13
Posts: 24
Joined: 16 Oct 2015, 17:15
Antispam: Yes

Re: Ranque-Hilsch Vortex Tube Modelling on ELMER

Post by serge_13 »

Dear friends some interesting results of 2019 with Ranque-Hilsch Vortex Tube (RHVT) with quadratic vortex channel is presented here. Sif file is included below. The non-stationary flow in RHVT was visualized in Paraview 5.7.
For reference see:
https://www.t-flow.ru/rhvt-quad-aero-all

Best regards,
Serge.

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

$function uy(c) {x=c(0); y=c(1); z=c(2); _uy=20*z/0.160}
$function uz(c) {x=c(0); y=c(1); z=c(2); _uz=-20*y/0.160}

$function press_1(c){x=c(0); y=c(1); z=c(2); arg_1=x/0.4095; arg_2=-(x-0.4315)/0.022; if(x < 0.4095) _press_1= 3.5e+5*arg_1; else _press_1= 3.5e+5*arg_2;}


Simulation
Max Output Level = 6
Coordinate System = Cartesian
Coordinate Mapping(3) = 1 2 3
Simulation Type = Transient
Steady State Max Iterations = 1
Timestepping Method = BDF

! Old ...
! Timestep Sizes(4) = 0.00001 0.00025 0.0001 0.00005
! Timestep Intervals(4) = 5 5 10 100

Timestep Sizes(3) = 0.00001 0.000025 0.0000025
Timestep Intervals(3) = 2 40 1000


! Timestep Sizes(1) = 0.0001
! Timestep Intervals(1) = 1000

! Output Intervals = 10
! Output File = myout
! Restart File = myout
! Restart Time = 4.10e-1

BDF Order = 1

! Post File = case.ep
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) = 11
Name = "Body 1"
Equation = 1
Material = 1
Initial condition = 1
End

Solver 2
Equation = Heat Equation
Procedure = "HeatSolve" "HeatSolver"
Variable = Temperature
Exec Solver = Always
Stabilize = False
Bubbles = False
Lumped Mass Matrix = False
Optimize Bandwidth = True
Steady State Convergence Tolerance = 1.0e-7
Nonlinear System Convergence Tolerance = 1.0e-7
Nonlinear System Max Iterations = 2
Nonlinear System Newton After Iterations = 20
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Relaxation Factor = 1.0
Linear System Solver = Iterative
Linear System Iterative Method = BiCGStab
Linear System Max Iterations = 300
Linear System Convergence Tolerance = 1.0e-10
Linear System Preconditioning = ILU1
Linear System ILUT Tolerance = 1.0e-3
Linear System Abort Not Converged = False
Linear System Residual Output = 1
Linear System Precondition Recompute = 1
Temperature Lower Limit = 250.0
Temperature Upper Limit = 450.0
End

Solver 1
Equation = Navier-Stokes
Procedure = "FlowSolve" "FlowSolver"
Variable = Flow Solution[Velocity:3 Pressure:1]
Exec Solver = Always
Stabilize = True
Bubbles = False
Lumped Mass Matrix = False
Optimize Bandwidth = True
! Div Discretization = True
Bouleev Anisotropic Viscosity = True
Steady State Convergence Tolerance = 1.0e-6
Nonlinear System Convergence Tolerance = 1.0e-6
Nonlinear System Max Iterations = 3
Nonlinear System Newton After Iterations = 15
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Relaxation Factor = 1
Linear System Solver = Iterative
Linear System Iterative Method = BiCGStab
Linear System Max Iterations = 500
Linear System Convergence Tolerance = 1.0e-7
Linear System Preconditioning = ILU1
Linear System ILUT Tolerance = 1.0e-3
Linear System Abort Not Converged = False
Linear System Residual Output = 1
Linear System Precondition Recompute = 1
End

Solver 3
Equation = Directed Scales
Exec Solver = Before Simulation
Procedure = "BouleevDirScales" "BouleevDirScales"
Variable = "Directed Scale"
Variable DOFs = 4
Nonlinear System Max Iterations = 100
Nonlinear System Convergence Tolerance = 1.0e-5
Display Model = False
Linear System Solver = Iterative
Steady State Convergence Tolerance = 1e-06
Ext Logging = False
Load from file = True
Save to file = True
End

Solver 4
Equation = Bouleev Anisotropic Viscosity
Procedure = "BouleevAnisotropic" "BouleevAnisotropic"
Variable = "AnisoTurbProps[TVis:3 TRe:1]"
Variable DOFs = 4
Nonlinear System Max Iterations = 1
Nonlinear System Convergence Tolerance = 1.0e-5
Display Model = False
Linear System Solver = Iterative
Steady State Convergence Tolerance = 1e-06
End

Solver 5
Equation = Divergence Solver
Procedure= "DivergenceSolver" "DivergenceSolver"
Target Variable = "Velocity"
Constant Bulk Matrix = True

Linear System Solver = Iterative
Linear System Iterative Method = cg
Linear System Preconditioning = None
Linear System Max Iterations = 500
Linear System Convergence Tolerance= 1.0e-10
End


Solver 6
Equation = CGNS Save Iterative

Procedure = "CGNSSaveSolveIterative" "CGNSSaveIterative"
Variable = "Dummy"

Nonlinear System Max Iterations = 100
Nonlinear System Convergence Tolerance = 1.0e-5

! Display Model = False
EveryNStep = 5
! CreateNewAfter = 100
EliminateSections = 1



Linear System Solver = Iterative
Steady State Convergence Tolerance = 1e-06
End


Equation 1
Name = "Equation 1"
Navier-Stokes = True
Heat Equation = True
Convection = Computed
Active Solvers(6) = 3 2 1 4 5 6
End

Material 1
Name = Air
Density = 1.293
Viscosity = 1.67e-5
Compressibility Model = "Perfect Gas"
Reference Pressure = 1.e+05
Reference Temperature = 293.0
Heat expansion Coefficient = 0.207e-3
Heat Conductivity = 0.02621
Sound speed = 330.0
Heat Capacity = 1.006e+03
End


Initial Condition 1
Name = "InitialCondition 1"

Velocity 1 = 0.0

Velocity 2 = Variable Coordinate
REAL MATC "uy(tx)"

Velocity 3 = Variable Coordinate
REAL MATC "uz(tx)"

! Velocity 2 = 0
! Velocity 3 = 0

Pressure = Variable Coordinate
REAL MATC "press_1(tx)"

Temperature = 293.0
TVis 1 = 1000
TVis 2 = 1000
TVis 3 = 1000
End

Boundary Condition 1
Target Boundaries(1) = 2
Name = "WALL"
Noslip wall BC = True
End

Boundary Condition 2
Target Boundaries(1) = 11
Name = "FLOW Input 1"
External Pressure = -2.5e+5

Normal-Tangential Velocity = True
Velocity 2 = 0
Velocity 3 = 0


Temperature = 293.0
End

Boundary Condition 3
Target Boundaries(1) = 1
Name = "FLOW Input 2"
External Pressure = -2.5e+5
Normal-Tangential Velocity = True
Velocity 2 = 0
Velocity 3 = 0
Temperature = 293.0
End


Boundary Condition 4
Target Boundaries(1) = 4
Name = "COLD Output"
Velocity 1 = 20.0
! Velocity 3 = 0
! Velocity 2 = 0
! External Pressure = 0
End

Boundary Condition 5
Target Boundaries(1) = 26
Name = "HOT Output"
! Normal-Tangential Velocity = Logical True
! Velocity 1 = Real 0.1
External Pressure = 0

End
Post Reply