diff --git a/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/3d-turbulent-taylor-couette.rst b/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/3d-turbulent-taylor-couette.rst index 8217b731e9..9433f049c6 100644 --- a/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/3d-turbulent-taylor-couette.rst +++ b/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/3d-turbulent-taylor-couette.rst @@ -9,7 +9,7 @@ Features --------- - Solvers: ``lethe-fluid-matrix-free`` (with Q2-Q2 or Q3-Q3) -- Transient problem using ``bdf3`` time integrator +- Transient problem using ``bdf2`` time integrator - Demonstrates the implementation of initial conditions for velocity and pressure - Demonstrates the postprocessing of the enstrophy @@ -237,16 +237,17 @@ To monitor the kinetic energy and the enstrophy, we set calculation to ``true`` Simulation Control ~~~~~~~~~~~~~~~~~~ -The ``simulation control`` subsection controls the flow of the simulation. To maximize the temporal accuracy of the simulation, we use a third order ``bdf3`` scheme. Results are written every 2 time-steps. To ensure a more adequate visualization of the high-order elements, we set ``subdivision = 2``. This will allow Paraview to render the high-order solutions with more fidelity. +The ``simulation control`` subsection controls the flow of the simulation. To maximize the temporal accuracy of the simulation, we use a second-order ``bdf2`` scheme. Results are written every 10 time-steps. To ensure a more adequate visualization of the high-order elements, we set ``subdivision = 2``. This will allow Paraview to render the high-order solutions with more fidelity. .. code-block:: text subsection simulation control - set method = bdf3 - set time step = 0.0125 - set number mesh adapt = 0 + set method = bdf2 + set time step = 0.01 + set adapt = true + set max cfl = 1 set time end = 60 - set output frequency = 2 + set output frequency = 10 set subdivision = 2 end diff --git a/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/images/enstrophy_comparison_655k_Q2.png b/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/images/enstrophy_comparison_655k_Q2.png index 4dc2d38a3b..1205fe375d 100644 Binary files a/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/images/enstrophy_comparison_655k_Q2.png and b/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/images/enstrophy_comparison_655k_Q2.png differ diff --git a/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/images/enstrophy_comparison_655k_Q3.png b/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/images/enstrophy_comparison_655k_Q3.png index 0406dbf215..6a717cf448 100644 Binary files a/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/images/enstrophy_comparison_655k_Q3.png and b/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/images/enstrophy_comparison_655k_Q3.png differ diff --git a/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/images/enstrophy_comparison_82k_Q2.png b/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/images/enstrophy_comparison_82k_Q2.png index 3cc9b5bcf2..c7e7636d3d 100644 Binary files a/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/images/enstrophy_comparison_82k_Q2.png and b/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/images/enstrophy_comparison_82k_Q2.png differ diff --git a/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/images/enstrophy_comparison_82k_Q3.png b/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/images/enstrophy_comparison_82k_Q3.png index 42fad81a78..4d048007c7 100644 Binary files a/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/images/enstrophy_comparison_82k_Q3.png and b/doc/source/examples/incompressible-flow/3d-turbulent-taylor-couette/images/enstrophy_comparison_82k_Q3.png differ diff --git a/examples/incompressible-flow/3d-turbulent-taylor-couette/tc-matrix-free.prm b/examples/incompressible-flow/3d-turbulent-taylor-couette/tc-matrix-free.prm index 02877386bc..64bf18790b 100644 --- a/examples/incompressible-flow/3d-turbulent-taylor-couette/tc-matrix-free.prm +++ b/examples/incompressible-flow/3d-turbulent-taylor-couette/tc-matrix-free.prm @@ -4,21 +4,23 @@ set dimension = 3 #--------------------------------------------------- -# Simulation Control +# Simulation Control #--------------------------------------------------- subsection simulation control - set method = bdf3 - set output name = taylor-couette - set output path = ./output/ - set time end = 60 - set time step = 0.0125 - set output frequency = 2 - set subdivision = 2 + set method = bdf2 + set output name = taylor-couette + set output path = ./output/ + set time end = 60 + set adapt = true + set max cfl = 1 + set time step = 0.01 + set output frequency = 10 + set subdivision = 2 end #------------------------ --------------------------- -# Physical Properties +# Physical Properties #--------------------------------------------------- subsection physical properties @@ -28,7 +30,7 @@ subsection physical properties end #--------------------------------------------------- -# Initial conditions +# Initial conditions #--------------------------------------------------- subsection initial conditions @@ -46,20 +48,19 @@ end #--------------------------------------------------- subsection mesh - set type = dealii + set type = dealii set grid type = cylinder_shell set grid arguments = 3.14159265359 : 0.5 : 1.0 : 5 : 4 : true set initial refinement = 4 end - #--------------------------------------------------- # Forces #--------------------------------------------------- subsection forces - set verbosity = quiet # Output force and torques in log - set calculate torque = true # Enable torque calculation + set verbosity = quiet # Output force and torques in log + set calculate torque = true # Enable torque calculation end #--------------------------------------------------- @@ -67,8 +68,8 @@ end #--------------------------------------------------- subsection FEM - set velocity order = 2 #3 for Q3 - set pressure order = 2 #3 for Q3 + set velocity order = 2 #3 for Q3 + set pressure order = 2 #3 for Q3 end #--------------------------------------------------- @@ -76,23 +77,21 @@ end #--------------------------------------------------- subsection post-processing - # Kinetic energy calculation - set calculate kinetic energy = true + set calculate kinetic energy = true # Enstrophy calculation - set calculate enstrophy = true + set calculate enstrophy = true end - #--------------------------------------------------- # Boundary Conditions #--------------------------------------------------- subsection boundary conditions - set number = 3 + set number = 3 set fix pressure constant = true - subsection bc 0 + subsection bc 0 set type = function subsection u set Function expression = -y @@ -104,10 +103,10 @@ subsection boundary conditions set Function expression = 0 end end - subsection bc 1 + subsection bc 1 set type = noslip end - subsection bc 2 + subsection bc 2 set type = periodic set id = 2 set periodic_id = 3 @@ -115,9 +114,18 @@ subsection boundary conditions end end +#--------------------------------------------------- +# Timer +#--------------------------------------------------- + +subsection timer + set type = iteration +end + #--------------------------------------------------- # Non-Linear Solver Control #--------------------------------------------------- + subsection non-linear solver subsection fluid dynamics set tolerance = 1e-3