Skip to content

Commit

Permalink
Add two more tests
Browse files Browse the repository at this point in the history
  • Loading branch information
@blaisb @oguevremont
blaisb authored and oguevremont committed Sep 24, 2024
1 parent d954dfd commit 5970820
Show file tree
Hide file tree
Showing 4 changed files with 522 additions and 0 deletions.
115 changes: 115 additions & 0 deletions applications_tests/lethe-fluid/dg_tracer_diffusion_mms_2d.output
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,115 @@
Running on 1 MPI rank(s)...
Number of active cells: 4
Number of degrees of freedom: 27
Volume of triangulation: 1
Number of tracer degrees of freedom: 16

*****************************
Steady iteration: 1/4
*****************************
---------------
Fluid Dynamics
---------------
Newton iteration: 0 - Residual: 0
-Tolerance of iterative solver is : 1e-14
-Iterative solver took : 0 steps to reach a residual norm of 0
alpha = 1 res = 0
||du||_L2 = 0 ||du||_Linfty = 0
||dp||_L2 = 0 ||dp||_Linfty = 0
-------
Tracer
-------
Newton iteration: 0 - Residual: 2.152
-Tolerance of iterative solver is : 2.152e-10
-Iterative solver took : 4 steps to reach a residual norm of 2.576e-13
alpha = 1 res = 2.573e-13 ||dx||_L2 = 2.382 ||dx||_Linfty = 1.185
L2 error velocity : 0
L2 error tracer : 0.00557821

*****************************
Steady iteration: 2/4
*****************************
Number of active cells: 16
Number of degrees of freedom: 75
Volume of triangulation: 1
Number of tracer degrees of freedom: 64
---------------
Fluid Dynamics
---------------
Newton iteration: 0 - Residual: 0
-Tolerance of iterative solver is : 1e-14
-Iterative solver took : 0 steps to reach a residual norm of 0
alpha = 1 res = 0
||du||_L2 = 0 ||du||_Linfty = 0
||dp||_L2 = 0 ||dp||_Linfty = 0
-------
Tracer
-------
Newton iteration: 0 - Residual: 1.193
-Tolerance of iterative solver is : 1.193e-10
-Iterative solver took : 13 steps to reach a residual norm of 7.414e-11
alpha = 1 res = 7.414e-11 ||dx||_L2 = 0.9338 ||dx||_Linfty = 0.1901
L2 error velocity : 0
L2 error tracer : 0.000580698

*****************************
Steady iteration: 3/4
*****************************
Number of active cells: 64
Number of degrees of freedom: 243
Volume of triangulation: 1
Number of tracer degrees of freedom: 256
---------------
Fluid Dynamics
---------------
Newton iteration: 0 - Residual: 0
-Tolerance of iterative solver is : 1e-14
-Iterative solver took : 0 steps to reach a residual norm of 0
alpha = 1 res = 0
||du||_L2 = 0 ||du||_Linfty = 0
||dp||_L2 = 0 ||dp||_Linfty = 0
-------
Tracer
-------
Newton iteration: 0 - Residual: 0.504
-Tolerance of iterative solver is : 5.04e-11
-Iterative solver took : 25 steps to reach a residual norm of 8.774e-12
alpha = 1 res = 8.774e-12 ||dx||_L2 = 0.5188 ||dx||_Linfty = 0.09607
L2 error velocity : 0
L2 error tracer : 4.0131e-05

*****************************
Steady iteration: 4/4
*****************************
Number of active cells: 256
Number of degrees of freedom: 867
Volume of triangulation: 1
Number of tracer degrees of freedom: 1024
---------------
Fluid Dynamics
---------------
Newton iteration: 0 - Residual: 0
-Tolerance of iterative solver is : 1e-14
-Iterative solver took : 0 steps to reach a residual norm of 0
alpha = 1 res = 0
||du||_L2 = 0 ||du||_Linfty = 0
||dp||_L2 = 0 ||dp||_Linfty = 0
-------
Tracer
-------
Newton iteration: 0 - Residual: 0.2278
-Tolerance of iterative solver is : 2.278e-11
-Iterative solver took : 40 steps to reach a residual norm of 2.141e-11
alpha = 1 res = 2.141e-11 ||dx||_L2 = 0.2654 ||dx||_Linfty = 0.02761
L2 error velocity : 0
L2 error tracer : 2.56353e-06
cells error_velocity error_pressure
4 0.000000e+00 - 0.000000e+00 -
16 0.000000e+00 nan 0.000000e+00 nan
64 0.000000e+00 nan 0.000000e+00 nan
256 0.000000e+00 nan 0.000000e+00 nan
cells error_tracer
4 5.578210e-03 -
16 5.806981e-04 3.26
64 4.013104e-05 3.85
256 2.563532e-06 3.97
177 changes: 177 additions & 0 deletions applications_tests/lethe-fluid/dg_tracer_diffusion_mms_2d.prm
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,177 @@
# Listing of Parameters
#----------------------

set dimension = 2

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

subsection simulation control
set method = steady
set number mesh adapt = 3
set output frequency = 0
end

#---------------------------------------------------
# FEM
#---------------------------------------------------

subsection FEM
set velocity order = 1
set pressure order = 1
set tracer order = 1
set tracer uses dg = true
end

#---------------------------------------------------
# Timer
#---------------------------------------------------

subsection timer
set type = none # <none|iteration|end>
end

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

subsection initial conditions
set type = nodal
subsection uvwp
set Function expression = 0; 0; 0
end
end

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

subsection physical properties
set number of fluids = 1
subsection fluid 0
set tracer diffusivity = 1
end
end

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

subsection mesh
set type = dealii
set grid type = hyper_rectangle
set grid arguments = 0, 0 : 1, 1.0 : false
set initial refinement = 1
end

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

subsection multiphysics
set tracer = true
end

#---------------------------------------------------
# Mesh Adaptation Control
#---------------------------------------------------

subsection mesh adaptation
set type = uniform
end

#---------------------------------------------------
# Analytical Solution
#---------------------------------------------------

subsection analytical solution
set enable = true
set verbosity = verbose
subsection uvwp
set Function expression = 0 ; 0 ; 0
end
subsection tracer
set Function expression = sin(pi*x)*sin(pi*y)
end
end

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

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

subsection boundary conditions tracer
set number = 1
subsection bc 0
set id = 0
set type = dirichlet
subsection dirichlet
set Function expression = 0
end
end
end

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

subsection source term
subsection tracer
set Function expression = 2*pi*pi*sin(pi*x)*sin(pi*y)
end
end

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

subsection non-linear solver
subsection tracer
set verbosity = verbose
set tolerance = 1e-10
set max iterations = 1
end
subsection fluid dynamics
set verbosity = verbose
set tolerance = 1e-10
set max iterations = 10
end
end

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

subsection linear solver
subsection fluid dynamics
set verbosity = verbose
set method = gmres
set relative residual = 1e-13
set minimum residual = 1e-14
set preconditioner = ilu
set ilu preconditioner fill = 0
set ilu preconditioner absolute tolerance = 1e-14
set ilu preconditioner relative tolerance = 1.00
end
subsection tracer
set verbosity = extra verbose
set method = gmres
set relative residual = 1e-10
set minimum residual = 1e-14
set preconditioner = ilu
set ilu preconditioner fill = 0
set ilu preconditioner absolute tolerance = 1e-14
set ilu preconditioner relative tolerance = 1.00
set max krylov vectors = 100
set max iters = 200

end
end
39 changes: 39 additions & 0 deletions applications_tests/lethe-fluid/dg_tracer_high_pe.output
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,39 @@
Running on 1 MPI rank(s)...
Number of active cells: 64
Number of degrees of freedom: 243
Volume of triangulation: 1
Number of tracer degrees of freedom: 256

*****************************
Steady iteration: 1/3
*****************************
L2 error velocity : 5.55112e-17
L2 error tracer : 0.0717243

*****************************
Steady iteration: 2/3
*****************************
Number of active cells: 256
Number of degrees of freedom: 867
Volume of triangulation: 1
Number of tracer degrees of freedom: 1024
L2 error velocity : 5.55112e-17
L2 error tracer : 0.011948

*****************************
Steady iteration: 3/3
*****************************
Number of active cells: 1024
Number of degrees of freedom: 3267
Volume of triangulation: 1
Number of tracer degrees of freedom: 4096
L2 error velocity : 5.55112e-17
L2 error tracer : 0.0013688
cells error_velocity error_pressure
64 5.551115e-17 - 0.000000e+00 -
256 5.551115e-17 0.00 0.000000e+00 nan
1024 5.551115e-17 0.00 0.000000e+00 nan
cells error_tracer
64 7.172431e-02 -
256 1.194799e-02 2.59
1024 1.368801e-03 3.13
Loading

0 comments on commit 5970820

Please sign in to comment.