Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Add compressibility to heat transfer #1051

Merged
merged 3 commits into from
Mar 2, 2024
Merged

Add compressibility to heat transfer #1051

Show file tree
Hide file tree
Changes from 2 commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
3812c64
Add compressibility to heat transfer
AmishgaAlphonius Feb 29, 2024
39fe148
Add test for VOF also
AmishgaAlphonius Feb 29, 2024
8bf3b78
Apply code review suggestion
AmishgaAlphonius Mar 1, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
32 changes: 32 additions & 0 deletions applications_tests/lethe-fluid/heat_transfer_vof_weakly_compressible_flow.output
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,32 @@
Running on 1 MPI rank(s)...
Number of active cells: 256
Number of degrees of freedom: 867
Volume of triangulation: 0.0625
Number of thermal degrees of freedom: 289
Number of VOF degrees of freedom: 289
Pressure drop: 0 Pa
Total pressure drop: 0 Pa

*******************************************************************************
Transient iteration: 1 Time: 0.1 Time step: 0.1 CFL: 0
*******************************************************************************
Pressure drop: -156.575 Pa
Total pressure drop: -156.575 Pa

*******************************************************************************
Transient iteration: 2 Time: 0.2 Time step: 0.1 CFL: 0.317772
*******************************************************************************
Pressure drop: -161.604 Pa
Total pressure drop: -161.604 Pa

*******************************************************************************
Transient iteration: 3 Time: 0.3 Time step: 0.1 CFL: 0.363124
*******************************************************************************
Pressure drop: -162.235 Pa
Total pressure drop: -162.235 Pa

*******************************************************************************
Transient iteration: 4 Time: 0.4 Time step: 0.1 CFL: 0.380725
*******************************************************************************
Pressure drop: -162.369 Pa
Total pressure drop: -162.369 Pa
205 changes: 205 additions & 0 deletions applications_tests/lethe-fluid/heat_transfer_vof_weakly_compressible_flow.prm
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,205 @@
# Listing of Parameters
#----------------------

set dimension = 2

#---------------------------------------------------
# Simulation Control
#---------------------------------------------------

subsection simulation control
set method = bdf1
set time end = 0.4
set time step = 0.1
set output frequency = 0
end

#---------------------------------------------------
# Multiphysics
#---------------------------------------------------

subsection multiphysics
set heat transfer = true
set VOF = true
end

#---------------------------------------------------
# Initial condition
#---------------------------------------------------

subsection initial conditions
subsection temperature
set Function expression = (40-25)*(0.25-x)*4+25
end
subsection VOF
set Function expression = if(y>=0.125, 0, 1)
end
end

#---------------------------------------------------
# Source term
#---------------------------------------------------

subsection source term
subsection xyz
set Function expression = 5 ; 0 ; 0
end
end

#---------------------------------------------------
# Physical Properties
#---------------------------------------------------

subsection physical properties
set number of fluids = 2
subsection fluid 0
set thermal conductivity = 0.5
set kinematic viscosity = 0.01
set thermal expansion = 0.0001
set specific heat = 100
set density model = isothermal_ideal_gas
subsection isothermal_ideal_gas
set density_ref = 10
set R = 287.05
set T = 298.15
end
end
subsection fluid 0
set thermal conductivity = 0.5
set kinematic viscosity = 0.01
set thermal expansion = 0.0001
set specific heat = 200
set density model = isothermal_ideal_gas
subsection isothermal_ideal_gas
set density_ref = 15
set R = 287.05
set T = 298.15
end
end
end

#---------------------------------------------------
# Mesh
#---------------------------------------------------

subsection mesh
set type = dealii
set grid type = hyper_cube
set grid arguments = 0 : 0.25 : true
set initial refinement = 4
end

#---------------------------------------------------
# Boundary Conditions
#---------------------------------------------------

subsection boundary conditions
set number = 4
subsection bc 0
set id = 0
set type = noslip
end
subsection bc 1
set id = 1
set type = noslip
end
subsection bc 2
set id = 2
set type = noslip
end
subsection bc 3
set id = 3
set type = noslip
end
end

subsection boundary conditions heat transfer
set number = 2
subsection bc 0
set id = 0
set type = temperature
subsection value
set Function expression = 40
end
end
subsection bc 1
set id = 1
set type = temperature
subsection value
set Function expression = 25
end
end
end

#---------------------------------------------------
# Post-processing
#---------------------------------------------------

subsection post-processing
set verbosity = verbose
set calculate pressure drop = true
set calculate mass conservation = false
end

#---------------------------------------------------
# Non-Linear Solver Control
#---------------------------------------------------

subsection non-linear solver
subsection heat transfer
set tolerance = 1e-3
set max iterations = 20
set verbosity = quiet
end
subsection fluid dynamics
set tolerance = 1e-4
set max iterations = 5
set verbosity = quiet
end
subsection VOF
set tolerance = 1e-9
set max iterations = 20
set verbosity = quiet
end
end

#---------------------------------------------------
# Linear Solver Control
#---------------------------------------------------

subsection linear solver
subsection fluid dynamics
set verbosity = quiet
set method = gmres
set max iters = 1000
set relative residual = 1e-3
set minimum residual = 1e-5
set preconditioner = ilu
set ilu preconditioner fill = 0
set ilu preconditioner absolute tolerance = 1e-12
set ilu preconditioner relative tolerance = 1.00
set max krylov vectors = 1000
end
subsection heat transfer
set verbosity = quiet
set method = gmres
set max iters = 1000
set relative residual = 1e-3
set minimum residual = 1e-5
set preconditioner = ilu
set ilu preconditioner fill = 0
set ilu preconditioner absolute tolerance = 1e-12
set ilu preconditioner relative tolerance = 1.00
end
subsection VOF
set verbosity = quiet
set method = gmres
set max iters = 1000
set relative residual = 1e-3
set minimum residual = 1e-10
set preconditioner = ilu
set ilu preconditioner fill = 0
set ilu preconditioner absolute tolerance = 1e-12
set ilu preconditioner relative tolerance = 1.00
end
end
31 changes: 31 additions & 0 deletions applications_tests/lethe-fluid/heat_transfer_weakly_compressible_flow.output
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,31 @@
Running on 1 MPI rank(s)...
Number of active cells: 256
Number of degrees of freedom: 867
Volume of triangulation: 0.0625
Number of thermal degrees of freedom: 289
Pressure drop: 0 Pa
Total pressure drop: 0 Pa

*******************************************************************************
Transient iteration: 1 Time: 0.1 Time step: 0.1 CFL: 0
*******************************************************************************
Pressure drop: -123.246 Pa
Total pressure drop: -123.246 Pa

**********************************************************************************
Transient iteration: 2 Time: 0.2 Time step: 0.1 CFL: 6.88457e-05
**********************************************************************************
Pressure drop: -123.485 Pa
Total pressure drop: -123.485 Pa

**********************************************************************************
Transient iteration: 3 Time: 0.3 Time step: 0.1 CFL: 4.78548e-05
**********************************************************************************
Pressure drop: -123.636 Pa
Total pressure drop: -123.636 Pa

**********************************************************************************
Transient iteration: 4 Time: 0.4 Time step: 0.1 CFL: 3.46342e-05
**********************************************************************************
Pressure drop: -123.757 Pa
Total pressure drop: -123.757 Pa
Loading
Loading