diff --git a/CHANGELOG.md b/CHANGELOG.md index d2ede6d0ca..cb70a9a481 100644 --- a/CHANGELOG.md +++ b/CHANGELOG.md @@ -3,6 +3,12 @@ All notable changes to the Lethe project will be documented in this file. The format is based on [Keep a Changelog](http://keepachangelog.com/). +## [Master] - 2024-11-28 + +### Changed + +- MINOR The squared term that was added in the correction direction vector of the VOF DCDD stabilization [#1103](https://github.com/chaos-polymtl/lethe/pull/1103) has been removed after observing regression in computed results of VOF examples. [#1390](https://github.com/chaos-polymtl/lethe/pull/1390) + ## [Master] - 2024-11-21 ### Changed diff --git a/applications_tests/lethe-fluid/gls_droplet_marangoni_effect.output b/applications_tests/lethe-fluid/gls_droplet_marangoni_effect.output index d23f8e0a39..5e172018b5 100644 --- a/applications_tests/lethe-fluid/gls_droplet_marangoni_effect.output +++ b/applications_tests/lethe-fluid/gls_droplet_marangoni_effect.output @@ -4,35 +4,35 @@ Running on 1 MPI rank(s)... Volume of triangulation: 0.0064 Number of thermal degrees of freedom: 4225 Number of VOF degrees of freedom: 4225 -Enstrophy : 0 -Kinetic energy : 0 +Enstrophy: 0 +Kinetic energy: 0 ******************************************************************************* -Transient iteration: 1 Time: 0.001 Time step: 0.001 CFL: 0 +Transient iteration: 1 Time: 0.001 Time step: 0.001 CFL: 0 ******************************************************************************* -Enstrophy : 0.42658 -Kinetic energy : 1.5988e-05 +Enstrophy: 0.42658 +Kinetic energy: 1.5988e-05 ******************************************************************************** Transient iteration: 2 Time: 0.002 Time step: 0.001 CFL: 0.0080812 ******************************************************************************** -Enstrophy : 97.119 -Kinetic energy : 0.0010102 +Enstrophy: 97.103 +Kinetic energy: 0.001011 ******************************************************************************* -Transient iteration: 3 Time: 0.003 Time step: 0.001 CFL: 0.17532 +Transient iteration: 3 Time: 0.003 Time step: 0.001 CFL: 0.17531 ******************************************************************************* -Enstrophy : 299.29 -Kinetic energy : 0.0039173 +Enstrophy: 299.47 +Kinetic energy: 0.0039228 ******************************************************************************* -Transient iteration: 4 Time: 0.004 Time step: 0.001 CFL: 0.30998 +Transient iteration: 4 Time: 0.004 Time step: 0.001 CFL: 0.31002 ******************************************************************************* -Enstrophy : 513.69 -Kinetic energy : 0.0079426 +Enstrophy: 513.9 +Kinetic energy: 0.0079503 ******************************************************************************* -Transient iteration: 5 Time: 0.005 Time step: 0.001 CFL: 0.40316 +Transient iteration: 5 Time: 0.005 Time step: 0.001 CFL: 0.40314 ******************************************************************************* -Enstrophy : 705.36 -Kinetic energy : 0.012344 +Enstrophy: 705.51 +Kinetic energy: 0.012354 diff --git a/applications_tests/lethe-fluid/gls_vof_1_isothermal_compressible_fluid.output b/applications_tests/lethe-fluid/gls_vof_1_isothermal_compressible_fluid.output index 9289d631e8..aface43f69 100644 --- a/applications_tests/lethe-fluid/gls_vof_1_isothermal_compressible_fluid.output +++ b/applications_tests/lethe-fluid/gls_vof_1_isothermal_compressible_fluid.output @@ -28,4 +28,4 @@ Transient iteration: 2 Time: 0.2 Time step: 0.1 CFL: 1.64181 VOF Mass Conservation ---------------------- time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -2.000000e-01 6.579242e-01 7.895070e-01 -2.229965e-02 4.765510e-02 3.420758e-01 3.420758e+02 1.863059e+01 -4.003717e+01 5.000000e-01 +2.000000e-01 6.579141e-01 7.894950e-01 -2.231310e-02 4.763344e-02 3.420859e-01 3.420859e+02 1.864223e+01 -4.002050e+01 5.000000e-01 diff --git a/applications_tests/lethe-fluid/gls_vof_2_isothermal_compressible_fluids.output b/applications_tests/lethe-fluid/gls_vof_2_isothermal_compressible_fluids.output index 3ad993b357..7f6eef16f0 100644 --- a/applications_tests/lethe-fluid/gls_vof_2_isothermal_compressible_fluids.output +++ b/applications_tests/lethe-fluid/gls_vof_2_isothermal_compressible_fluids.output @@ -28,4 +28,4 @@ Transient iteration: 2 Time: 0.2 Time step: 0.1 CFL: 1.53178 VOF Mass Conservation ---------------------- time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -2.000000e-01 6.579652e-01 7.895557e-01 -2.004270e-02 4.448404e-02 3.420348e-01 1.710176e+01 8.361416e-01 -1.863500e+00 5.000000e-01 +2.000000e-01 6.579587e-01 7.895478e-01 -2.005946e-02 4.446851e-02 3.420413e-01 1.710209e+01 8.365478e-01 -1.862997e+00 5.000000e-01 diff --git a/applications_tests/lethe-fluid/gls_vof_dirichlet_boundary_condition.output b/applications_tests/lethe-fluid/gls_vof_dirichlet_boundary_condition.output index 8633599240..770f87b559 100644 --- a/applications_tests/lethe-fluid/gls_vof_dirichlet_boundary_condition.output +++ b/applications_tests/lethe-fluid/gls_vof_dirichlet_boundary_condition.output @@ -18,44 +18,44 @@ Transient iteration: 1 Time: 0.05 Time step: 0.05 CFL: 2.23662 VOF Mass Conservation ---------------------- time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -5.000000e-02 2.023368e-02 2.023368e+00 -2.916883e-03 -2.022545e+00 9.797663e-01 2.939299e+02 3.285580e-03 -2.938668e+02 5.000000e-01 +5.000000e-02 2.022997e-02 2.022997e+00 -2.904127e-03 -2.022178e+00 9.797700e-01 2.939310e+02 3.247343e-03 -2.938679e+02 5.000000e-01 ******************************************************************************* -Transient iteration: 2 Time: 0.1 Time step: 0.05 CFL: 2.93874 +Transient iteration: 2 Time: 0.1 Time step: 0.05 CFL: 2.93873 ******************************************************************************* ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -1.000000e-01 3.991214e-02 3.991214e+00 -3.514015e-03 -3.993601e+00 9.600879e-01 2.880264e+02 3.315622e-03 -2.880185e+02 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +1.000000e-01 4.000801e-02 4.000801e+00 -3.499492e-03 -4.003198e+00 9.599920e-01 2.879976e+02 3.271545e-03 -2.879897e+02 5.000000e-01 ******************************************************************************* -Transient iteration: 3 Time: 0.15 Time step: 0.05 CFL: 2.87622 +Transient iteration: 3 Time: 0.15 Time step: 0.05 CFL: 2.87622 ******************************************************************************* ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -1.500000e-01 6.362990e-02 6.362990e+00 -3.677075e-03 -6.366046e+00 9.363701e-01 2.809110e+02 3.288925e-03 -2.809018e+02 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +1.500000e-01 6.365731e-02 6.365731e+00 -3.662590e-03 -6.368794e+00 9.363427e-01 2.809028e+02 3.244410e-03 -2.808936e+02 5.000000e-01 ******************************************************************************* -Transient iteration: 4 Time: 0.2 Time step: 0.05 CFL: 2.86088 +Transient iteration: 4 Time: 0.2 Time step: 0.05 CFL: 2.86088 ******************************************************************************* ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -2.000000e-01 8.856688e-02 8.856688e+00 -3.716378e-03 -8.860082e+00 9.114331e-01 2.734299e+02 2.660659e-03 -2.734197e+02 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +2.000000e-01 8.852284e-02 8.852284e+00 -3.702138e-03 -8.855684e+00 9.114772e-01 2.734431e+02 2.617015e-03 -2.734329e+02 5.000000e-01 ******************************************************************************* -Transient iteration: 5 Time: 0.25 Time step: 0.05 CFL: 2.85364 +Transient iteration: 5 Time: 0.25 Time step: 0.05 CFL: 2.85365 ******************************************************************************* ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -2.500000e-01 1.138431e-01 1.138431e+01 -3.712321e-03 -1.138794e+01 8.861569e-01 2.658471e+02 1.672484e-03 -2.658361e+02 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +2.500000e-01 1.137550e-01 1.137550e+01 -3.698586e-03 -1.137913e+01 8.862450e-01 2.658735e+02 1.630379e-03 -2.658626e+02 5.000000e-01 diff --git a/applications_tests/lethe-fluid/gls_vof_ht_adapt-dissipation.output b/applications_tests/lethe-fluid/gls_vof_ht_adapt-dissipation.output index de238004e4..ae4d418960 100644 --- a/applications_tests/lethe-fluid/gls_vof_ht_adapt-dissipation.output +++ b/applications_tests/lethe-fluid/gls_vof_ht_adapt-dissipation.output @@ -8,52 +8,52 @@ Running on 1 MPI rank(s)... ******************************************************************************* Transient iteration: 1 Time: 0.02 Time step: 0.02 CFL: 1.11831 ******************************************************************************* -L2 error velocity : 0.99151 +L2 error velocity: 0.99151 L2 error temperature : 4.95887e-13 L2 error phase : 0.00279114 ******************************************************************************* -Transient iteration: 2 Time: 0.04 Time step: 0.02 CFL: 1.25786 +Transient iteration: 2 Time: 0.04 Time step: 0.02 CFL: 1.25786 ******************************************************************************* -L2 error velocity : 1.09115 -L2 error temperature : 4.68382e-12 -L2 error phase : 0.00284754 +L2 error velocity: 1.091 +L2 error temperature : 4.73269e-12 +L2 error phase : 0.00286026 ******************************************************************************* -Transient iteration: 3 Time: 0.06 Time step: 0.02 CFL: 1.28682 +Transient iteration: 3 Time: 0.06 Time step: 0.02 CFL: 1.28679 ******************************************************************************* -L2 error velocity : 1.12524 -L2 error temperature : 2.14831e-11 -L2 error phase : 0.00292325 +L2 error velocity: 1.12496 +L2 error temperature : 2.16659e-11 +L2 error phase : 0.00295017 ******************************************************************************* -Transient iteration: 4 Time: 0.08 Time step: 0.02 CFL: 1.29512 +Transient iteration: 4 Time: 0.08 Time step: 0.02 CFL: 1.29508 ******************************************************************************* -L2 error velocity : 1.13541 -L2 error temperature : 6.41882e-11 -L2 error phase : 0.00301194 +L2 error velocity: 1.13507 +L2 error temperature : 6.43194e-11 +L2 error phase : 0.00305334 ******************************************************************************* -Transient iteration: 5 Time: 0.1 Time step: 0.02 CFL: 1.2975 -******************************************************************************* -L2 error velocity : 1.137 -L2 error temperature : 1.43343e-10 -L2 error phase : 0.0031096 - time error_velocity error_pressure -0.0200 9.915105e-01 12.1545 -0.0400 1.091154e+00 11.7510 -0.0600 1.125244e+00 11.6955 -0.0800 1.135406e+00 11.8167 -0.1000 1.136995e+00 11.9363 -cells error_temperature - 256 4.958868e-13 - 256 4.683819e-12 - 256 2.148313e-11 - 256 6.418823e-11 - 256 1.433433e-10 - time error_phase -0.0200 2.791137e-03 -0.0400 2.847543e-03 -0.0600 2.923250e-03 -0.0800 3.011940e-03 -0.1000 3.109602e-03 +Transient iteration: 5 Time: 0.1 Time step: 0.02 CFL: 1.29745 +******************************************************************************* +L2 error velocity: 1.13666 +L2 error temperature : 1.42948e-10 +L2 error phase : 0.00316573 + time error_velocity error_pressure +0.0200 9.915105e-01 12.1545 +0.0400 1.091004e+00 11.7602 +0.0600 1.124961e+00 11.7061 +0.0800 1.135072e+00 11.8248 +0.1000 1.136664e+00 11.9419 +cells error_temperature + 256 4.958868e-13 + 256 4.732688e-12 + 256 2.166589e-11 + 256 6.431939e-11 + 256 1.429483e-10 + time error_phase +0.0200 2.791137e-03 +0.0400 2.860260e-03 +0.0600 2.950171e-03 +0.0800 3.053344e-03 +0.1000 3.165733e-03 diff --git a/applications_tests/lethe-fluid/gls_vof_hydrostat_initial_refine.output b/applications_tests/lethe-fluid/gls_vof_hydrostat_initial_refine.output index 076246303f..cfa0452cbd 100644 --- a/applications_tests/lethe-fluid/gls_vof_hydrostat_initial_refine.output +++ b/applications_tests/lethe-fluid/gls_vof_hydrostat_initial_refine.output @@ -47,7 +47,7 @@ Running on 1 MPI rank(s)... ******************************************************************************* Transient iteration: 1 Time: 0.02 Time step: 0.02 CFL: 0 ******************************************************************************* -L2 error velocity : 1.09509e-07 +L2 error velocity: 1.09509e-07 L2 error phase : 0.000174446 ********************************************************************************** @@ -57,23 +57,23 @@ Transient iteration: 2 Time: 0.04 Time step: 0.02 CFL: 7.71695e-0 Number of degrees of freedom: 6948 Volume of triangulation: 1 Number of VOF degrees of freedom: 2316 -L2 error velocity : 1.09511e-07 -L2 error phase : 0.000174449 +L2 error velocity: 1.0951e-07 +L2 error phase : 0.000174451 -********************************************************************************** -Transient iteration: 3 Time: 0.06 Time step: 0.02 CFL: 7.71707e-06 -********************************************************************************** +******************************************************************************** +Transient iteration: 3 Time: 0.06 Time step: 0.02 CFL: 7.717e-06 +******************************************************************************** Number of active cells: 2032 Number of degrees of freedom: 6948 Volume of triangulation: 1 Number of VOF degrees of freedom: 2316 -L2 error velocity : 1.0951e-07 -L2 error phase : 0.000174452 - time error_velocity error_pressure -0.0200 1.095089e-07 0.0028 -0.0400 1.095107e-07 0.0028 -0.0600 1.095098e-07 0.0028 - time error_phase -0.0200 1.744461e-04 -0.0400 1.744490e-04 -0.0600 1.744518e-04 +L2 error velocity: 1.09508e-07 +L2 error phase : 0.000174457 + time error_velocity error_pressure +0.0200 1.095089e-07 0.0028 +0.0400 1.095098e-07 0.0028 +0.0600 1.095081e-07 0.0028 + time error_phase +0.0200 1.744461e-04 +0.0400 1.744515e-04 +0.0600 1.744569e-04 diff --git a/applications_tests/lethe-fluid/gls_vof_hydrostat_mesh-adapt.output b/applications_tests/lethe-fluid/gls_vof_hydrostat_mesh-adapt.output index 0a95e320fa..68632f363d 100644 --- a/applications_tests/lethe-fluid/gls_vof_hydrostat_mesh-adapt.output +++ b/applications_tests/lethe-fluid/gls_vof_hydrostat_mesh-adapt.output @@ -7,7 +7,7 @@ Running on 1 MPI rank(s)... ******************************************************************************* Transient iteration: 1 Time: 0.02 Time step: 0.02 CFL: 0 ******************************************************************************* -L2 error velocity : 1.98025e-05 +L2 error velocity: 1.98025e-05 L2 error phase : 0.0104167 ********************************************************************************** @@ -17,47 +17,47 @@ Transient iteration: 2 Time: 0.04 Time step: 0.02 CFL: 6.16399e-0 Number of degrees of freedom: 1596 Volume of triangulation: 1 Number of VOF degrees of freedom: 532 -L2 error velocity : 3.09853e-06 -L2 error phase : 0.00425559 +L2 error velocity: 3.0985e-06 +L2 error phase : 0.0042558 -********************************************************************************** -Transient iteration: 3 Time: 0.06 Time step: 0.02 CFL: 1.63314e-05 -********************************************************************************** +********************************************************************************* +Transient iteration: 3 Time: 0.06 Time step: 0.02 CFL: 1.6332e-05 +********************************************************************************* Number of active cells: 400 Number of degrees of freedom: 1437 Volume of triangulation: 1 Number of VOF degrees of freedom: 479 -L2 error velocity : 2.63133e-06 -L2 error phase : 0.00425566 +L2 error velocity: 2.63119e-06 +L2 error phase : 0.00425592 ********************************************************************************** -Transient iteration: 4 Time: 0.08 Time step: 0.02 CFL: 1.41975e-05 +Transient iteration: 4 Time: 0.08 Time step: 0.02 CFL: 1.41952e-05 ********************************************************************************** Number of active cells: 568 Number of degrees of freedom: 1950 Volume of triangulation: 1 Number of VOF degrees of freedom: 650 -L2 error velocity : 2.06792e-06 -L2 error phase : 0.00425573 +L2 error velocity: 2.06779e-06 +L2 error phase : 0.00425606 ********************************************************************************** -Transient iteration: 5 Time: 0.1 Time step: 0.02 CFL: 1.31728e-05 +Transient iteration: 5 Time: 0.1 Time step: 0.02 CFL: 1.31716e-05 ********************************************************************************** Number of active cells: 952 Number of degrees of freedom: 3300 Volume of triangulation: 1 Number of VOF degrees of freedom: 1100 -L2 error velocity : 1.55492e-06 -L2 error phase : 0.0042558 - time error_velocity error_pressure -0.0200 1.980252e-05 0.0684 -0.0400 3.098533e-06 0.0668 -0.0600 2.631326e-06 0.0669 -0.0800 2.067917e-06 0.0669 -0.1000 1.554917e-06 0.0669 - time error_phase -0.0200 1.041667e-02 -0.0400 4.255591e-03 -0.0600 4.255663e-03 -0.0800 4.255733e-03 -0.1000 4.255803e-03 +L2 error velocity: 1.55472e-06 +L2 error phase : 0.00425618 + time error_velocity error_pressure +0.0200 1.980252e-05 0.0684 +0.0400 3.098495e-06 0.0668 +0.0600 2.631190e-06 0.0669 +0.0800 2.067786e-06 0.0669 +0.1000 1.554718e-06 0.0669 + time error_phase +0.0200 1.041667e-02 +0.0400 4.255796e-03 +0.0600 4.255924e-03 +0.0800 4.256056e-03 +0.1000 4.256177e-03 diff --git a/applications_tests/lethe-fluid/gls_vof_periodic_boundary_condition.output b/applications_tests/lethe-fluid/gls_vof_periodic_boundary_condition.output index f675c584df..68a693454e 100644 --- a/applications_tests/lethe-fluid/gls_vof_periodic_boundary_condition.output +++ b/applications_tests/lethe-fluid/gls_vof_periodic_boundary_condition.output @@ -28,34 +28,34 @@ Transient iteration: 2 Time: 0.1 Time step: 0.05 CFL: 1.39102 VOF Mass Conservation ---------------------- time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -1.000000e-01 5.000000e-01 5.000000e+01 2.136042e-09 -4.905000e+01 5.000000e-01 1.500000e+02 -2.990494e-08 -1.471500e+02 5.000000e-01 +1.000000e-01 5.000000e-01 5.000000e+01 2.504256e-09 -4.905000e+01 5.000000e-01 1.500000e+02 -2.949111e-08 -1.471500e+02 5.000000e-01 ******************************************************************************* -Transient iteration: 3 Time: 0.15 Time step: 0.05 CFL: 2.78204 +Transient iteration: 3 Time: 0.15 Time step: 0.05 CFL: 2.78204 ******************************************************************************* ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -1.500000e-01 5.000000e-01 5.000000e+01 -4.189748e-09 -7.357500e+01 5.000000e-01 1.500000e+02 -5.910110e-08 -2.207250e+02 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +1.500000e-01 5.000000e-01 5.000000e+01 -3.430902e-09 -7.357500e+01 5.000000e-01 1.500000e+02 -5.777233e-08 -2.207250e+02 5.000000e-01 ******************************************************************************* -Transient iteration: 4 Time: 0.2 Time step: 0.05 CFL: 4.17307 +Transient iteration: 4 Time: 0.2 Time step: 0.05 CFL: 4.17307 ******************************************************************************* ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -2.000000e-01 5.000000e-01 5.000000e+01 3.434894e-09 -9.810001e+01 5.000000e-01 1.500000e+02 -9.171519e-08 -2.943000e+02 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +2.000000e-01 5.000000e-01 5.000000e+01 4.282132e-09 -9.809999e+01 5.000000e-01 1.500000e+02 -8.945293e-08 -2.943000e+02 5.000000e-01 ******************************************************************************* -Transient iteration: 5 Time: 0.25 Time step: 0.05 CFL: 5.56409 +Transient iteration: 5 Time: 0.25 Time step: 0.05 CFL: 5.56409 ******************************************************************************* ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -2.500000e-01 5.000000e-01 5.000000e+01 2.125296e-08 -1.226250e+02 5.000000e-01 1.500000e+02 -5.811637e-08 -3.678750e+02 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +2.500000e-01 4.999999e-01 4.999999e+01 2.209575e-08 -1.226250e+02 5.000001e-01 1.500000e+02 -5.519808e-08 -3.678750e+02 5.000000e-01 diff --git a/applications_tests/lethe-fluid/gls_vof_static_droplet_surface_tension_force_with_tanh_filter.output b/applications_tests/lethe-fluid/gls_vof_static_droplet_surface_tension_force_with_tanh_filter.output index 6f8fd8b835..68c75069a9 100644 --- a/applications_tests/lethe-fluid/gls_vof_static_droplet_surface_tension_force_with_tanh_filter.output +++ b/applications_tests/lethe-fluid/gls_vof_static_droplet_surface_tension_force_with_tanh_filter.output @@ -3,7 +3,7 @@ Running on 1 MPI rank(s)... Number of degrees of freedom: 3267 Volume of triangulation: 0.0064 Number of VOF degrees of freedom: 1089 -Filtered phase values: +Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 @@ -1095,9 +1095,9 @@ Filtered phase values: 2.06115e-09 ******************************************************************************* -Transient iteration: 1 Time: 0.001 Time step: 0.001 CFL: 0 +Transient iteration: 1 Time: 0.001 Time step: 0.001 CFL: 0 ******************************************************************************* -Filtered phase values: +Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 @@ -2187,18 +2187,13 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -L2 error velocity : 3.83829e-07 +L2 error velocity: 3.83829e-07 L2 error phase : 3.40945e-06 ********************************************************************************** -Transient iteration: 2 Time: 0.002 Time step: 0.001 CFL: 6.72415e-06 +Transient iteration: 2 Time: 0.002 Time step: 0.001 CFL: 6.72416e-06 ********************************************************************************** -Filtered phase values: -2.06115e-09 -2.06115e-09 -2.06115e-09 -2.06115e-09 -2.06115e-09 +Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 @@ -2277,6 +2272,9 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 +2.06116e-09 +2.06116e-09 +2.06114e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -2284,16 +2282,18 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06117e-09 +2.06116e-09 2.06115e-09 +2.06118e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 +2.06116e-09 2.06115e-09 -2.06117e-09 -2.06114e-09 +2.06118e-09 +2.06112e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -2302,40 +2302,40 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06116e-09 -2.06113e-09 +2.06117e-09 +2.0611e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 +2.06117e-09 +2.0611e-09 +2.06116e-09 2.06116e-09 -2.06113e-09 -2.06114e-09 -2.06114e-09 0.981292 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 -2.06115e-09 -2.06115e-09 2.06116e-09 -2.06111e-09 -2.0613e-09 -2.06115e-09 2.06115e-09 +2.06118e-09 +2.06106e-09 +2.06154e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06116e-09 -2.06111e-09 -2.0613e-09 -0.514685 +2.06115e-09 +2.06118e-09 +2.06106e-09 +2.06154e-09 +0.514665 1 -0.514685 +0.514665 1 1 2.06115e-09 @@ -2345,17 +2345,17 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06118e-09 -1.18953e-05 +2.06114e-09 +2.0612e-09 +1.1895e-05 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06118e-09 -1.18953e-05 +2.06114e-09 +2.0612e-09 +1.1895e-05 1 1 1 @@ -2372,7 +2372,7 @@ Filtered phase values: 2.06114e-09 2.06122e-09 2.06092e-09 -0.881063 +0.881062 2.06115e-09 2.06115e-09 2.06115e-09 @@ -2380,7 +2380,7 @@ Filtered phase values: 2.06114e-09 2.06122e-09 2.06092e-09 -0.881063 +0.881062 1 1 1 @@ -2396,17 +2396,17 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06116e-09 -2.06113e-09 -2.06121e-09 +2.06115e-09 +2.06115e-09 +2.06115e-09 0.999993 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 -2.06116e-09 -2.06113e-09 -2.06121e-09 +2.06115e-09 +2.06115e-09 +2.06115e-09 0.999993 1 1 @@ -2424,18 +2424,18 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06118e-09 -2.06105e-09 -2.93057e-07 +2.06114e-09 +2.06119e-09 +2.06102e-09 +2.93074e-07 1 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06118e-09 -2.06105e-09 -2.93057e-07 +2.06114e-09 +2.06119e-09 +2.06102e-09 +2.93074e-07 1 1 1 @@ -2454,19 +2454,19 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06115e-09 2.06116e-09 2.06114e-09 -2.06119e-09 -1.85633e-06 +2.06122e-09 +2.06091e-09 +1.85739e-06 1 2.06115e-09 2.06115e-09 -2.06115e-09 2.06116e-09 2.06114e-09 -2.06119e-09 -1.85633e-06 +2.06122e-09 +2.06091e-09 +1.85739e-06 1 1 1 @@ -2488,10 +2488,10 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06118e-09 -2.06105e-09 -2.93057e-07 +2.06114e-09 +2.06119e-09 +2.06102e-09 +2.93074e-07 1 1 1 @@ -2504,10 +2504,10 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06118e-09 -2.06105e-09 -2.93057e-07 +2.06114e-09 +2.06119e-09 +2.06102e-09 +2.93074e-07 1 1 1 @@ -2524,9 +2524,9 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06116e-09 -2.06113e-09 -2.06121e-09 +2.06115e-09 +2.06115e-09 +2.06115e-09 0.999993 1 1 @@ -2540,9 +2540,9 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06116e-09 -2.06113e-09 -2.06121e-09 +2.06115e-09 +2.06115e-09 +2.06115e-09 0.999993 1 1 @@ -2564,7 +2564,7 @@ Filtered phase values: 2.06114e-09 2.06122e-09 2.06092e-09 -0.881063 +0.881062 1 1 1 @@ -2580,7 +2580,7 @@ Filtered phase values: 2.06114e-09 2.06122e-09 2.06092e-09 -0.881063 +0.881062 1 1 1 @@ -2601,9 +2601,9 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06118e-09 -1.18953e-05 +2.06114e-09 +2.0612e-09 +1.1895e-05 1 1 1 @@ -2617,9 +2617,9 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06118e-09 -1.18953e-05 +2.06114e-09 +2.0612e-09 +1.1895e-05 1 1 1 @@ -2640,12 +2640,12 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06115e-09 -2.06115e-09 2.06116e-09 -2.06111e-09 -2.0613e-09 -0.514685 +2.06115e-09 +2.06118e-09 +2.06106e-09 +2.06154e-09 +0.514665 1 1 1 @@ -2656,12 +2656,12 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06115e-09 -2.06115e-09 2.06116e-09 -2.06111e-09 -2.0613e-09 -0.514685 +2.06115e-09 +2.06118e-09 +2.06106e-09 +2.06154e-09 +0.514665 1 1 1 @@ -2686,9 +2686,9 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 +2.06117e-09 +2.0611e-09 2.06116e-09 -2.06113e-09 -2.06114e-09 0.981292 1 1 @@ -2702,9 +2702,9 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 +2.06117e-09 +2.0611e-09 2.06116e-09 -2.06113e-09 -2.06114e-09 0.981292 1 1 @@ -2730,12 +2730,12 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 +2.06116e-09 2.06115e-09 -2.06115e-09 -2.06117e-09 -2.06114e-09 -2.06114e-09 -0.514685 +2.06118e-09 +2.06112e-09 +2.06116e-09 +0.514665 1 1 1 @@ -2746,12 +2746,12 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 +2.06116e-09 2.06115e-09 -2.06115e-09 -2.06117e-09 -2.06114e-09 -2.06114e-09 -0.514685 +2.06118e-09 +2.06112e-09 +2.06116e-09 +0.514665 1 1 1 @@ -2765,13 +2765,11 @@ Filtered phase values: 1 1 1 -0.514685 -2.06114e-09 -0.514685 -2.06114e-09 -2.06114e-09 -2.06115e-09 -2.06115e-09 +0.514665 +2.06116e-09 +0.514665 +2.06116e-09 +2.06112e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -2780,11 +2778,13 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06117e-09 -2.06113e-09 -2.0613e-09 -1.18954e-05 -0.881063 +2.06116e-09 +2.06114e-09 +2.06118e-09 +2.0611e-09 +2.06154e-09 +1.1895e-05 +0.881062 0.999993 1 1 @@ -2794,33 +2794,31 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06115e-09 -2.06115e-09 -2.06117e-09 -2.06113e-09 -2.0613e-09 -1.18954e-05 -0.881063 +2.06116e-09 +2.06114e-09 +2.06118e-09 +2.0611e-09 +2.06154e-09 +1.1895e-05 +0.881062 0.999993 1 1 1 0.999993 -0.881063 -1.18954e-05 +0.881062 +1.1895e-05 1 0.999993 -0.881063 -1.18954e-05 -2.0613e-09 -2.06113e-09 -2.0613e-09 -2.06113e-09 -2.06117e-09 -2.06117e-09 -2.06115e-09 -2.06115e-09 -2.06115e-09 +0.881062 +1.1895e-05 +2.06154e-09 +2.0611e-09 +2.06154e-09 +2.0611e-09 +2.06118e-09 +2.06118e-09 +2.06114e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -2831,14 +2829,14 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06111e-09 -2.06118e-09 -2.06092e-09 -2.06121e-09 -2.93056e-07 -1.85633e-06 2.06115e-09 +2.06117e-09 +2.06106e-09 +2.0612e-09 +2.06092e-09 2.06115e-09 +2.93074e-07 +1.85738e-06 2.06115e-09 2.06115e-09 2.06115e-09 @@ -2847,30 +2845,30 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06111e-09 -2.06118e-09 -2.06092e-09 -2.06121e-09 -2.93056e-07 -1.85633e-06 -2.93056e-07 -2.06121e-09 +2.06115e-09 +2.06117e-09 +2.06106e-09 +2.0612e-09 2.06092e-09 -2.06118e-09 -2.06111e-09 -2.06116e-09 2.06115e-09 +2.93074e-07 +1.85738e-06 +2.93074e-07 2.06115e-09 -2.93056e-07 -2.06121e-09 2.06092e-09 -2.06118e-09 -2.06111e-09 -2.06116e-09 -2.06115e-09 +2.0612e-09 +2.06106e-09 +2.06117e-09 2.06115e-09 +2.06116e-09 +2.93074e-07 2.06115e-09 +2.06092e-09 +2.0612e-09 +2.06106e-09 +2.06117e-09 2.06115e-09 +2.06116e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -2884,12 +2882,12 @@ Filtered phase values: 2.06115e-09 2.06116e-09 2.06115e-09 +2.06118e-09 +2.06114e-09 2.06122e-09 -2.06113e-09 -2.06105e-09 -2.06119e-09 -2.06115e-09 2.06115e-09 +2.06102e-09 +2.06091e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -2900,21 +2898,25 @@ Filtered phase values: 2.06115e-09 2.06116e-09 2.06115e-09 +2.06118e-09 +2.06114e-09 2.06122e-09 -2.06113e-09 -2.06105e-09 -2.06119e-09 -2.06105e-09 -2.06113e-09 +2.06115e-09 +2.06102e-09 +2.06091e-09 +2.06102e-09 +2.06115e-09 2.06122e-09 +2.06114e-09 +2.06118e-09 2.06115e-09 2.06116e-09 2.06115e-09 +2.06102e-09 2.06115e-09 -2.06115e-09 -2.06105e-09 -2.06113e-09 2.06122e-09 +2.06114e-09 +2.06118e-09 2.06115e-09 2.06116e-09 2.06115e-09 @@ -2934,13 +2936,11 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06115e-09 -2.06115e-09 2.06116e-09 2.06114e-09 -2.06116e-09 -2.06118e-09 -2.06114e-09 +2.06115e-09 +2.06119e-09 +2.06122e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -2954,25 +2954,21 @@ Filtered phase values: 2.06115e-09 2.06116e-09 2.06114e-09 -2.06116e-09 -2.06118e-09 -2.06114e-09 -2.06118e-09 -2.06116e-09 -2.06114e-09 -2.06116e-09 -2.06115e-09 -2.06115e-09 2.06115e-09 +2.06119e-09 +2.06122e-09 +2.06119e-09 2.06115e-09 -2.06118e-09 -2.06116e-09 2.06114e-09 2.06116e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 +2.06119e-09 +2.06115e-09 +2.06114e-09 +2.06116e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -2992,12 +2988,14 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06116e-09 2.06115e-09 2.06115e-09 2.06116e-09 2.06115e-09 +2.06116e-09 2.06115e-09 +2.06114e-09 +2.06114e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -3010,21 +3008,23 @@ Filtered phase values: 2.06115e-09 2.06116e-09 2.06115e-09 -2.06115e-09 2.06116e-09 2.06115e-09 +2.06114e-09 +2.06114e-09 +2.06114e-09 2.06115e-09 2.06116e-09 2.06115e-09 +2.06116e-09 2.06115e-09 2.06115e-09 2.06115e-09 -2.06115e-09 -2.06115e-09 +2.06114e-09 2.06115e-09 2.06116e-09 2.06115e-09 -2.06115e-09 +2.06116e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -3052,7 +3052,7 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 +2.06116e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -3068,7 +3068,7 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 +2.06116e-09 2.06116e-09 2.06115e-09 2.06115e-09 @@ -3283,13 +3283,13 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -L2 error velocity : 4.80362e-07 -L2 error phase : 3.40951e-06 +L2 error velocity: 4.80365e-07 +L2 error phase : 3.40956e-06 ********************************************************************************** Transient iteration: 3 Time: 0.003 Time step: 0.001 CFL: 8.34842e-06 ********************************************************************************** -Filtered phase values: +Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 @@ -3370,7 +3370,7 @@ Filtered phase values: 2.06115e-09 2.06116e-09 2.06116e-09 -2.06114e-09 +2.06113e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -3379,37 +3379,37 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06118e-09 +2.06114e-09 +2.0612e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06118e-09 -2.06112e-09 +2.06114e-09 +2.0612e-09 +2.06108e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 +2.06116e-09 2.06115e-09 +2.06118e-09 +2.06105e-09 2.06115e-09 -2.06117e-09 -2.06111e-09 2.06115e-09 2.06115e-09 2.06115e-09 +2.06116e-09 2.06115e-09 +2.06118e-09 +2.06105e-09 2.06115e-09 2.06115e-09 -2.06117e-09 -2.06111e-09 -2.06111e-09 -2.06111e-09 0.981289 2.06115e-09 2.06115e-09 @@ -3417,21 +3417,21 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06118e-09 -2.06107e-09 -2.06146e-09 +2.06114e-09 +2.06121e-09 +2.06096e-09 +2.06193e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06118e-09 -2.06107e-09 -2.06146e-09 -0.514702 +2.06114e-09 +2.06121e-09 +2.06096e-09 +2.06193e-09 +0.514661 1 -0.514702 +0.514661 1 1 2.06115e-09 @@ -3441,17 +3441,17 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 -2.06122e-09 -1.18961e-05 +2.06113e-09 +2.06124e-09 +1.18953e-05 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 -2.06122e-09 -1.18961e-05 +2.06113e-09 +2.06124e-09 +1.18953e-05 1 1 1 @@ -3467,16 +3467,16 @@ Filtered phase values: 2.06116e-09 2.06112e-09 2.06129e-09 -2.06065e-09 -0.881038 +2.06066e-09 +0.881037 2.06115e-09 2.06115e-09 2.06115e-09 2.06116e-09 2.06112e-09 2.06129e-09 -2.06065e-09 -0.881038 +2.06066e-09 +0.881037 1 1 1 @@ -3492,17 +3492,17 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06117e-09 +2.06115e-09 +2.06115e-09 2.06111e-09 -2.06127e-09 0.999993 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 -2.06117e-09 +2.06115e-09 +2.06115e-09 2.06111e-09 -2.06127e-09 0.999993 1 1 @@ -3520,18 +3520,18 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 -2.06121e-09 -2.06093e-09 -2.93143e-07 +2.06113e-09 +2.06124e-09 +2.06083e-09 +2.93205e-07 1 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 -2.06121e-09 -2.06093e-09 -2.93143e-07 +2.06113e-09 +2.06124e-09 +2.06083e-09 +2.93205e-07 1 1 1 @@ -3550,19 +3550,19 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06115e-09 2.06116e-09 -2.06113e-09 -2.06123e-09 -1.85659e-06 +2.06111e-09 +2.06131e-09 +2.06058e-09 +1.85904e-06 1 2.06115e-09 2.06115e-09 -2.06115e-09 2.06116e-09 -2.06113e-09 -2.06123e-09 -1.85659e-06 +2.06111e-09 +2.06131e-09 +2.06058e-09 +1.85904e-06 1 1 1 @@ -3584,10 +3584,10 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 -2.06121e-09 -2.06093e-09 -2.93143e-07 +2.06113e-09 +2.06124e-09 +2.06083e-09 +2.93205e-07 1 1 1 @@ -3600,10 +3600,10 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 -2.06121e-09 -2.06093e-09 -2.93143e-07 +2.06113e-09 +2.06124e-09 +2.06083e-09 +2.93205e-07 1 1 1 @@ -3620,9 +3620,9 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06117e-09 +2.06115e-09 +2.06115e-09 2.06111e-09 -2.06127e-09 0.999993 1 1 @@ -3636,9 +3636,9 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06117e-09 +2.06115e-09 +2.06115e-09 2.06111e-09 -2.06127e-09 0.999993 1 1 @@ -3659,8 +3659,8 @@ Filtered phase values: 2.06116e-09 2.06112e-09 2.06129e-09 -2.06065e-09 -0.881038 +2.06066e-09 +0.881037 1 1 1 @@ -3675,8 +3675,8 @@ Filtered phase values: 2.06116e-09 2.06112e-09 2.06129e-09 -2.06065e-09 -0.881038 +2.06066e-09 +0.881037 1 1 1 @@ -3697,9 +3697,9 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 -2.06122e-09 -1.18961e-05 +2.06113e-09 +2.06124e-09 +1.18953e-05 1 1 1 @@ -3713,9 +3713,9 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 -2.06122e-09 -1.18961e-05 +2.06113e-09 +2.06124e-09 +1.18953e-05 1 1 1 @@ -3737,11 +3737,11 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06118e-09 -2.06107e-09 -2.06146e-09 -0.514702 +2.06114e-09 +2.06121e-09 +2.06096e-09 +2.06193e-09 +0.514661 1 1 1 @@ -3753,11 +3753,11 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06118e-09 -2.06107e-09 -2.06146e-09 -0.514702 +2.06114e-09 +2.06121e-09 +2.06096e-09 +2.06193e-09 +0.514661 1 1 1 @@ -3780,11 +3780,11 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 +2.06116e-09 2.06115e-09 +2.06118e-09 +2.06105e-09 2.06115e-09 -2.06117e-09 -2.06111e-09 -2.06111e-09 0.981289 1 1 @@ -3796,11 +3796,11 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 +2.06116e-09 2.06115e-09 +2.06118e-09 +2.06105e-09 2.06115e-09 -2.06117e-09 -2.06111e-09 -2.06111e-09 0.981289 1 1 @@ -3827,11 +3827,11 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 +2.06114e-09 +2.0612e-09 +2.06108e-09 2.06115e-09 -2.06118e-09 -2.06112e-09 -2.06111e-09 -0.514702 +0.514661 1 1 1 @@ -3843,11 +3843,11 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 +2.06114e-09 +2.0612e-09 +2.06108e-09 2.06115e-09 -2.06118e-09 -2.06112e-09 -2.06111e-09 -0.514702 +0.514661 1 1 1 @@ -3861,11 +3861,11 @@ Filtered phase values: 1 1 1 -0.514702 -2.06111e-09 -0.514702 -2.06111e-09 -2.06112e-09 +0.514661 +2.06115e-09 +0.514661 +2.06115e-09 +2.06108e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -3875,12 +3875,12 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 -2.06118e-09 -2.06111e-09 -2.06146e-09 -1.18961e-05 -0.881038 +2.06113e-09 +2.0612e-09 +2.06105e-09 +2.06193e-09 +1.18953e-05 +0.881037 0.999993 1 1 @@ -3891,30 +3891,30 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 -2.06118e-09 -2.06111e-09 -2.06146e-09 -1.18961e-05 -0.881038 +2.06113e-09 +2.0612e-09 +2.06105e-09 +2.06193e-09 +1.18953e-05 +0.881037 0.999993 1 1 1 0.999993 -0.881038 -1.18961e-05 +0.881037 +1.18953e-05 1 0.999993 -0.881038 -1.18961e-05 -2.06146e-09 -2.06111e-09 -2.06146e-09 -2.06111e-09 -2.06118e-09 -2.06118e-09 -2.06114e-09 +0.881037 +1.18953e-05 +2.06193e-09 +2.06105e-09 +2.06193e-09 +2.06105e-09 +2.0612e-09 +2.0612e-09 +2.06113e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -3925,14 +3925,14 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06117e-09 -2.06107e-09 -2.06122e-09 -2.06065e-09 -2.06127e-09 -2.93142e-07 -1.85658e-06 +2.06114e-09 +2.06118e-09 +2.06096e-09 +2.06124e-09 +2.06066e-09 +2.06111e-09 +2.93205e-07 +1.85901e-06 2.06115e-09 2.06115e-09 2.06115e-09 @@ -3941,29 +3941,29 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06117e-09 -2.06107e-09 -2.06122e-09 -2.06065e-09 -2.06127e-09 -2.93142e-07 -1.85658e-06 -2.93142e-07 -2.06127e-09 -2.06065e-09 -2.06122e-09 -2.06107e-09 -2.06117e-09 -2.06115e-09 +2.06114e-09 +2.06118e-09 +2.06096e-09 +2.06124e-09 +2.06066e-09 +2.06111e-09 +2.93205e-07 +1.85901e-06 +2.93205e-07 +2.06111e-09 +2.06066e-09 +2.06124e-09 +2.06096e-09 +2.06118e-09 +2.06114e-09 2.06116e-09 -2.93142e-07 -2.06127e-09 -2.06065e-09 -2.06122e-09 -2.06107e-09 -2.06117e-09 -2.06115e-09 +2.93205e-07 +2.06111e-09 +2.06066e-09 +2.06124e-09 +2.06096e-09 +2.06118e-09 +2.06114e-09 2.06116e-09 2.06115e-09 2.06115e-09 @@ -3978,12 +3978,12 @@ Filtered phase values: 2.06115e-09 2.06116e-09 2.06115e-09 -2.06118e-09 -2.06114e-09 +2.06121e-09 +2.06113e-09 2.06129e-09 -2.06111e-09 -2.06093e-09 -2.06123e-09 +2.06115e-09 +2.06083e-09 +2.06059e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -3994,29 +3994,27 @@ Filtered phase values: 2.06115e-09 2.06116e-09 2.06115e-09 -2.06118e-09 -2.06114e-09 -2.06129e-09 -2.06111e-09 -2.06093e-09 -2.06123e-09 -2.06093e-09 -2.06111e-09 +2.06121e-09 +2.06113e-09 2.06129e-09 -2.06114e-09 -2.06118e-09 2.06115e-09 -2.06116e-09 +2.06083e-09 +2.06059e-09 +2.06083e-09 2.06115e-09 -2.06093e-09 -2.06111e-09 2.06129e-09 -2.06114e-09 -2.06118e-09 +2.06113e-09 +2.06121e-09 2.06115e-09 2.06116e-09 2.06115e-09 +2.06083e-09 +2.06115e-09 +2.06129e-09 +2.06113e-09 +2.06121e-09 2.06115e-09 +2.06116e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -4033,12 +4031,12 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 +2.06114e-09 +2.06116e-09 2.06112e-09 -2.06117e-09 -2.06121e-09 -2.06113e-09 -2.06115e-09 2.06115e-09 +2.06124e-09 +2.0613e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -4049,24 +4047,26 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 +2.06114e-09 +2.06116e-09 2.06112e-09 -2.06117e-09 -2.06121e-09 -2.06113e-09 -2.06121e-09 -2.06117e-09 +2.06115e-09 +2.06124e-09 +2.0613e-09 +2.06124e-09 +2.06115e-09 2.06112e-09 2.06116e-09 +2.06114e-09 +2.06116e-09 2.06115e-09 2.06115e-09 +2.06124e-09 2.06115e-09 -2.06115e-09 -2.06121e-09 -2.06117e-09 2.06112e-09 2.06116e-09 -2.06115e-09 -2.06115e-09 +2.06114e-09 +2.06116e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -4090,8 +4090,8 @@ Filtered phase values: 2.06115e-09 2.06116e-09 2.06115e-09 -2.06114e-09 -2.06116e-09 +2.06113e-09 +2.06111e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -4106,9 +4106,9 @@ Filtered phase values: 2.06115e-09 2.06116e-09 2.06115e-09 -2.06114e-09 -2.06116e-09 -2.06114e-09 +2.06113e-09 +2.06111e-09 +2.06113e-09 2.06115e-09 2.06116e-09 2.06115e-09 @@ -4116,7 +4116,7 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06114e-09 +2.06113e-09 2.06115e-09 2.06116e-09 2.06115e-09 @@ -4148,7 +4148,7 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 +2.06116e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -4164,7 +4164,7 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 +2.06116e-09 2.06116e-09 2.06115e-09 2.06115e-09 @@ -4379,13 +4379,13 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -L2 error velocity : 5.17376e-07 -L2 error phase : 3.40958e-06 +L2 error velocity: 5.17381e-07 +L2 error phase : 3.40971e-06 ********************************************************************************** -Transient iteration: 4 Time: 0.004 Time step: 0.001 CFL: 9.07931e-06 +Transient iteration: 4 Time: 0.004 Time step: 0.001 CFL: 9.07917e-06 ********************************************************************************** -Filtered phase values: +Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 @@ -4466,7 +4466,7 @@ Filtered phase values: 2.06115e-09 2.06116e-09 2.06116e-09 -2.06114e-09 +2.06112e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -4476,7 +4476,7 @@ Filtered phase values: 2.06115e-09 2.06116e-09 2.06114e-09 -2.06119e-09 +2.06122e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -4484,8 +4484,8 @@ Filtered phase values: 2.06115e-09 2.06116e-09 2.06114e-09 -2.06119e-09 -2.06111e-09 +2.06122e-09 +2.06105e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -4493,19 +4493,19 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06117e-09 -2.06109e-09 +2.06114e-09 +2.06119e-09 +2.06101e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06117e-09 -2.06109e-09 -2.06106e-09 -2.06106e-09 +2.06114e-09 +2.06119e-09 +2.06101e-09 +2.06112e-09 +2.06112e-09 0.981284 2.06115e-09 2.06115e-09 @@ -4513,21 +4513,21 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 -2.06119e-09 -2.06102e-09 -2.06162e-09 +2.06113e-09 +2.06123e-09 +2.06087e-09 +2.06229e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 -2.06119e-09 -2.06102e-09 -2.06162e-09 -0.514716 +2.06113e-09 +2.06123e-09 +2.06087e-09 +2.06229e-09 +0.514657 1 -0.514716 +0.514657 1 1 2.06115e-09 @@ -4537,17 +4537,17 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06112e-09 -2.06126e-09 -1.18968e-05 +2.06111e-09 +2.06129e-09 +1.18955e-05 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 2.06116e-09 -2.06112e-09 -2.06126e-09 -1.18968e-05 +2.06111e-09 +2.06129e-09 +1.18955e-05 1 1 1 @@ -4562,17 +4562,17 @@ Filtered phase values: 2.06115e-09 2.06117e-09 2.0611e-09 -2.06137e-09 -2.06038e-09 -0.881012 +2.06136e-09 +2.0604e-09 +0.881011 2.06115e-09 2.06115e-09 2.06115e-09 2.06117e-09 2.0611e-09 -2.06137e-09 -2.06038e-09 -0.881012 +2.06136e-09 +2.0604e-09 +0.881011 1 1 1 @@ -4588,17 +4588,17 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06117e-09 -2.06109e-09 -2.06134e-09 +2.06115e-09 +2.06116e-09 +2.06105e-09 0.999993 2.06115e-09 2.06115e-09 2.06115e-09 2.06115e-09 -2.06117e-09 -2.06109e-09 -2.06134e-09 +2.06115e-09 +2.06116e-09 +2.06105e-09 0.999993 1 1 @@ -4616,18 +4616,18 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06113e-09 -2.06124e-09 -2.06081e-09 -2.93232e-07 +2.06111e-09 +2.0613e-09 +2.06061e-09 +2.93357e-07 1 2.06115e-09 2.06115e-09 2.06116e-09 -2.06113e-09 -2.06124e-09 -2.06081e-09 -2.93232e-07 +2.06111e-09 +2.0613e-09 +2.06061e-09 +2.93357e-07 1 1 1 @@ -4644,21 +4644,21 @@ Filtered phase values: 1 2.06115e-09 2.06115e-09 -2.06115e-09 -2.06115e-09 -2.06115e-09 2.06116e-09 -2.06112e-09 -2.06127e-09 -1.85689e-06 -1 -2.06115e-09 -2.06115e-09 2.06115e-09 +2.06117e-09 +2.06109e-09 +2.0614e-09 +2.06022e-09 +1.86086e-06 +1 2.06116e-09 -2.06112e-09 -2.06127e-09 -1.85689e-06 +2.06115e-09 +2.06117e-09 +2.06109e-09 +2.0614e-09 +2.06022e-09 +1.86086e-06 1 1 1 @@ -4680,10 +4680,10 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06113e-09 -2.06124e-09 -2.06081e-09 -2.93232e-07 +2.06111e-09 +2.0613e-09 +2.06061e-09 +2.93357e-07 1 1 1 @@ -4696,10 +4696,10 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06113e-09 -2.06124e-09 -2.06081e-09 -2.93232e-07 +2.06111e-09 +2.0613e-09 +2.06061e-09 +2.93357e-07 1 1 1 @@ -4716,9 +4716,9 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06117e-09 -2.06109e-09 -2.06134e-09 +2.06115e-09 +2.06116e-09 +2.06105e-09 0.999993 1 1 @@ -4732,9 +4732,9 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06117e-09 -2.06109e-09 -2.06134e-09 +2.06115e-09 +2.06116e-09 +2.06105e-09 0.999993 1 1 @@ -4754,9 +4754,9 @@ Filtered phase values: 2.06115e-09 2.06117e-09 2.0611e-09 -2.06137e-09 -2.06038e-09 -0.881012 +2.06136e-09 +2.0604e-09 +0.881011 1 1 1 @@ -4770,9 +4770,9 @@ Filtered phase values: 2.06115e-09 2.06117e-09 2.0611e-09 -2.06137e-09 -2.06038e-09 -0.881012 +2.06136e-09 +2.0604e-09 +0.881011 1 1 1 @@ -4793,9 +4793,9 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06112e-09 -2.06126e-09 -1.18968e-05 +2.06111e-09 +2.06129e-09 +1.18955e-05 1 1 1 @@ -4809,9 +4809,9 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06112e-09 -2.06126e-09 -1.18968e-05 +2.06111e-09 +2.06129e-09 +1.18955e-05 1 1 1 @@ -4833,11 +4833,11 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 -2.06119e-09 -2.06102e-09 -2.06162e-09 -0.514716 +2.06113e-09 +2.06123e-09 +2.06087e-09 +2.06229e-09 +0.514657 1 1 1 @@ -4849,11 +4849,11 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 -2.06119e-09 -2.06102e-09 -2.06162e-09 -0.514716 +2.06113e-09 +2.06123e-09 +2.06087e-09 +2.06229e-09 +0.514657 1 1 1 @@ -4877,10 +4877,10 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06117e-09 -2.06109e-09 -2.06106e-09 +2.06114e-09 +2.06119e-09 +2.06101e-09 +2.06112e-09 0.981284 1 1 @@ -4893,10 +4893,10 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06117e-09 -2.06109e-09 -2.06106e-09 +2.06114e-09 +2.06119e-09 +2.06101e-09 +2.06112e-09 0.981284 1 1 @@ -4924,10 +4924,10 @@ Filtered phase values: 2.06115e-09 2.06116e-09 2.06114e-09 -2.06119e-09 -2.06111e-09 -2.06106e-09 -0.514716 +2.06122e-09 +2.06105e-09 +2.06112e-09 +0.514657 1 1 1 @@ -4940,10 +4940,10 @@ Filtered phase values: 2.06115e-09 2.06116e-09 2.06114e-09 -2.06119e-09 -2.06111e-09 -2.06106e-09 -0.514716 +2.06122e-09 +2.06105e-09 +2.06112e-09 +0.514657 1 1 1 @@ -4957,11 +4957,11 @@ Filtered phase values: 1 1 1 -0.514716 -2.06106e-09 -0.514716 -2.06106e-09 -2.06111e-09 +0.514657 +2.06112e-09 +0.514657 +2.06112e-09 +2.06105e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -4971,12 +4971,12 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 -2.06119e-09 -2.06109e-09 -2.06162e-09 -1.18968e-05 -0.881012 +2.06112e-09 +2.06122e-09 +2.06101e-09 +2.06229e-09 +1.18955e-05 +0.881011 0.999993 1 1 @@ -4987,30 +4987,30 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 -2.06119e-09 -2.06109e-09 -2.06162e-09 -1.18968e-05 -0.881012 +2.06112e-09 +2.06122e-09 +2.06101e-09 +2.06229e-09 +1.18955e-05 +0.881011 0.999993 1 1 1 0.999993 -0.881012 -1.18968e-05 +0.881011 +1.18955e-05 1 0.999993 -0.881012 -1.18968e-05 -2.06162e-09 -2.06109e-09 -2.06162e-09 -2.06109e-09 -2.06119e-09 -2.06119e-09 -2.06114e-09 +0.881011 +1.18955e-05 +2.06229e-09 +2.06101e-09 +2.06229e-09 +2.06101e-09 +2.06122e-09 +2.06122e-09 +2.06112e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -5022,13 +5022,13 @@ Filtered phase values: 2.06115e-09 2.06116e-09 2.06114e-09 -2.06117e-09 -2.06102e-09 -2.06126e-09 -2.06038e-09 -2.06134e-09 -2.93231e-07 -1.85688e-06 +2.06119e-09 +2.06087e-09 +2.06129e-09 +2.06041e-09 +2.06105e-09 +2.93358e-07 +1.86082e-06 2.06115e-09 2.06115e-09 2.06115e-09 @@ -5038,27 +5038,27 @@ Filtered phase values: 2.06115e-09 2.06116e-09 2.06114e-09 -2.06117e-09 -2.06102e-09 -2.06126e-09 -2.06038e-09 -2.06134e-09 -2.93231e-07 -1.85688e-06 -2.93231e-07 -2.06134e-09 -2.06038e-09 -2.06126e-09 -2.06102e-09 -2.06117e-09 +2.06119e-09 +2.06087e-09 +2.06129e-09 +2.06041e-09 +2.06105e-09 +2.93358e-07 +1.86082e-06 +2.93358e-07 +2.06105e-09 +2.06041e-09 +2.06129e-09 +2.06087e-09 +2.06119e-09 2.06114e-09 2.06116e-09 -2.93231e-07 -2.06134e-09 -2.06038e-09 -2.06126e-09 -2.06102e-09 -2.06117e-09 +2.93358e-07 +2.06105e-09 +2.06041e-09 +2.06129e-09 +2.06087e-09 +2.06119e-09 2.06114e-09 2.06116e-09 2.06115e-09 @@ -5073,13 +5073,13 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06119e-09 -2.06112e-09 -2.06137e-09 -2.06109e-09 -2.06082e-09 -2.06127e-09 +2.06114e-09 +2.06123e-09 +2.06111e-09 +2.06136e-09 +2.06116e-09 +2.06061e-09 +2.06023e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -5089,27 +5089,27 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 -2.06119e-09 -2.06112e-09 -2.06137e-09 -2.06109e-09 -2.06082e-09 -2.06127e-09 -2.06082e-09 -2.06109e-09 -2.06137e-09 -2.06112e-09 -2.06119e-09 -2.06115e-09 +2.06114e-09 +2.06123e-09 +2.06111e-09 +2.06136e-09 +2.06116e-09 +2.06061e-09 +2.06023e-09 +2.06061e-09 +2.06116e-09 +2.06136e-09 +2.06111e-09 +2.06123e-09 +2.06114e-09 2.06116e-09 2.06115e-09 -2.06082e-09 -2.06109e-09 -2.06137e-09 -2.06112e-09 -2.06119e-09 -2.06115e-09 +2.06061e-09 +2.06116e-09 +2.06136e-09 +2.06111e-09 +2.06123e-09 +2.06114e-09 2.06116e-09 2.06115e-09 2.06115e-09 @@ -5127,12 +5127,12 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 +2.06113e-09 2.06116e-09 2.0611e-09 -2.06117e-09 -2.06124e-09 -2.06112e-09 +2.06115e-09 +2.0613e-09 +2.0614e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -5143,25 +5143,25 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06114e-09 +2.06113e-09 2.06116e-09 2.0611e-09 -2.06117e-09 -2.06124e-09 -2.06112e-09 -2.06124e-09 -2.06117e-09 +2.06115e-09 +2.0613e-09 +2.0614e-09 +2.0613e-09 +2.06115e-09 2.0611e-09 2.06116e-09 -2.06114e-09 +2.06113e-09 2.06116e-09 2.06115e-09 2.06115e-09 -2.06124e-09 -2.06117e-09 +2.0613e-09 +2.06115e-09 2.0611e-09 2.06116e-09 -2.06114e-09 +2.06113e-09 2.06116e-09 2.06115e-09 2.06115e-09 @@ -5186,8 +5186,8 @@ Filtered phase values: 2.06115e-09 2.06117e-09 2.06115e-09 -2.06113e-09 -2.06116e-09 +2.06111e-09 +2.06109e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -5202,9 +5202,9 @@ Filtered phase values: 2.06115e-09 2.06117e-09 2.06115e-09 -2.06113e-09 -2.06116e-09 -2.06113e-09 +2.06111e-09 +2.06109e-09 +2.06111e-09 2.06115e-09 2.06117e-09 2.06115e-09 @@ -5212,7 +5212,7 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -2.06113e-09 +2.06111e-09 2.06115e-09 2.06117e-09 2.06115e-09 @@ -5244,7 +5244,7 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 +2.06117e-09 2.06115e-09 2.06115e-09 2.06115e-09 @@ -5260,7 +5260,7 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06116e-09 -2.06115e-09 +2.06117e-09 2.06116e-09 2.06115e-09 2.06115e-09 @@ -5475,15 +5475,15 @@ Filtered phase values: 2.06115e-09 2.06115e-09 2.06115e-09 -L2 error velocity : 5.36113e-07 -L2 error phase : 3.40966e-06 - time error_velocity error_pressure -0.0010 3.838286e-07 0.0124 -0.0020 4.803615e-07 0.0125 -0.0030 5.173757e-07 0.0125 -0.0040 5.361132e-07 0.0125 - time error_phase -0.0010 3.409448e-06 -0.0020 3.409511e-06 -0.0030 3.409583e-06 -0.0040 3.409655e-06 +L2 error velocity: 5.36117e-07 +L2 error phase : 3.40986e-06 + time error_velocity error_pressure +0.0010 3.838286e-07 0.0124 +0.0020 4.803652e-07 0.0125 +0.0030 5.173809e-07 0.0125 +0.0040 5.361174e-07 0.0125 + time error_phase +0.0010 3.409448e-06 +0.0020 3.409563e-06 +0.0030 3.409706e-06 +0.0040 3.409861e-06 diff --git a/applications_tests/lethe-fluid/gls_vof_tanh_filter.output b/applications_tests/lethe-fluid/gls_vof_tanh_filter.output index a4e8aba797..74df8c8730 100644 --- a/applications_tests/lethe-fluid/gls_vof_tanh_filter.output +++ b/applications_tests/lethe-fluid/gls_vof_tanh_filter.output @@ -146,69 +146,69 @@ Filtered phase values: Transient iteration: 2 Time: 0.1 Time step: 0.05 CFL: 0.0034067 ******************************************************************************** Filtered phase values: -0.999663 -0.999667 -0.999664 +0.99967 +0.999654 +0.999668 +0.999659 +0.000364357 0.999666 -0.000344323 -0.999665 -0.000341521 -0.999662 -0.00034173 -0.00033297 +0.000348685 +0.99966 +0.000341457 +0.000327928 0.999668 -0.000333707 -0.000333652 +0.000331845 +0.000333724 0.999664 -0.000339663 -0.000334198 -0.000336033 -0.999661 -0.000335821 -0.000335837 -0.00033568 +0.00033892 +0.000334395 +0.0003375 +0.99966 +0.000336358 +0.000335816 +0.000335623 0.999652 -0.000338936 -0.000334391 -0.000335625 -0.000335009 -0.000335115 -0.000335107 -0.000335185 -0.000335213 -0.000329252 -0.000335331 -0.000330554 -0.000335843 -0.000335218 -0.000335388 -0.000337003 -0.000335355 -0.000336647 -0.000335218 -0.000335386 +0.000339174 +0.000334328 +0.000335643 +0.00033428 +0.000334847 +0.000335117 +0.000335214 +0.000335204 +0.000329536 0.00033534 -0.000334908 -0.000335349 -0.000335003 -0.000335385 -0.000335341 -0.000335353 -0.000335469 -0.00033535 -0.000335444 -0.000335341 +0.000329975 +0.000335613 +0.00033528 +0.00033537 +0.000336925 0.000335353 +0.000336805 +0.00033528 +0.000335369 +0.000335345 +0.000334929 0.000335349 -0.000335316 +0.000334961 +0.000335369 +0.000335345 +0.000335352 +0.000335464 +0.00033535 +0.000335455 +0.000335345 +0.000335351 +0.00033535 +0.000335318 +0.00033535 +0.00033532 +0.000335352 0.00033535 -0.000335323 -0.000335353 -0.000335349 0.00033535 -0.000335367 +0.000335366 0.00033535 -0.000335363 +0.000335365 0.000335349 0.00033535 0.00033535 diff --git a/applications_tests/lethe-fluid/heat_transfer_vof_phase_change_constrain_solid_domain.output b/applications_tests/lethe-fluid/heat_transfer_vof_phase_change_constrain_solid_domain.output index fd44c4e08a..50f55f2c51 100644 --- a/applications_tests/lethe-fluid/heat_transfer_vof_phase_change_constrain_solid_domain.output +++ b/applications_tests/lethe-fluid/heat_transfer_vof_phase_change_constrain_solid_domain.output @@ -65,13 +65,13 @@ Total pressure drop: -0.0113477 Pa VOF Mass Conservation ---------------------- time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -3.000000e-02 5.912813e-01 5.912813e-01 5.080623e-07 -2.095655e-06 1.227187e-01 1.227187e-01 -1.074572e-07 -1.036538e-08 5.000000e-01 +3.000000e-02 5.912813e-01 5.912813e-01 5.080623e-07 -2.095656e-06 1.227187e-01 1.227187e-01 -1.074573e-07 -1.036449e-08 5.000000e-01 --------------- VOF Barycenter --------------- - time x_vof y_vof vx_vof vy_vof -3.000000e-02 5.000000e-01 6.152842e-02 -8.756382e-07 -8.446450e-08 + time x_vof y_vof vx_vof vy_vof +3.000000e-02 5.000000e-01 6.152842e-02 -8.756391e-07 -8.445725e-08 ********************************************************************************** Transient iteration: 4 Time: 0.04 Time step: 0.01 CFL: 0.000185963 @@ -82,11 +82,11 @@ Total pressure drop: -0.0104538 Pa ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -4.000000e-02 5.912813e-01 5.912813e-01 5.180054e-07 -2.109152e-06 1.227187e-01 1.227187e-01 -1.545717e-07 -2.010638e-08 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +4.000000e-02 5.912813e-01 5.912813e-01 5.180058e-07 -2.109154e-06 1.227187e-01 1.227187e-01 -1.545721e-07 -2.010420e-08 5.000000e-01 --------------- VOF Barycenter --------------- - time x_vof y_vof vx_vof vy_vof -4.000000e-02 5.000000e-01 6.152842e-02 -1.259561e-06 -1.638412e-07 + time x_vof y_vof vx_vof vy_vof +4.000000e-02 5.000000e-01 6.152843e-02 -1.259564e-06 -1.638235e-07 diff --git a/applications_tests/lethe-fluid/heat_transfer_vof_phase_change_constrain_solid_domain_with_plane.output b/applications_tests/lethe-fluid/heat_transfer_vof_phase_change_constrain_solid_domain_with_plane.output index 76f445bc8b..b0843834b6 100644 --- a/applications_tests/lethe-fluid/heat_transfer_vof_phase_change_constrain_solid_domain_with_plane.output +++ b/applications_tests/lethe-fluid/heat_transfer_vof_phase_change_constrain_solid_domain_with_plane.output @@ -47,46 +47,46 @@ Total pressure drop: -0.00497741 Pa VOF Mass Conservation ---------------------- time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -2.000000e-02 5.912813e-01 5.912813e-01 7.205532e-07 -2.095201e-06 1.227187e-01 1.227187e-01 -5.507064e-07 -1.576839e-07 5.000000e-01 +2.000000e-02 5.912813e-01 5.912813e-01 7.205516e-07 -2.095203e-06 1.227187e-01 1.227187e-01 -5.507048e-07 -1.576816e-07 5.000000e-01 --------------- VOF Barycenter --------------- - time x_vof y_vof vx_vof vy_vof -2.000000e-02 5.000001e-01 6.152840e-02 -4.487550e-06 -1.284921e-06 + time x_vof y_vof vx_vof vy_vof +2.000000e-02 5.000001e-01 6.152841e-02 -4.487537e-06 -1.284902e-06 ********************************************************************************** Transient iteration: 3 Time: 0.03 Time step: 0.01 CFL: 0.000176621 ********************************************************************************** -Pressure drop: -0.00508754 Pa -Total pressure drop: -0.00508754 Pa +Pressure drop: -0.00508755 Pa +Total pressure drop: -0.00508755 Pa ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -3.000000e-02 5.912813e-01 5.912813e-01 7.519397e-07 -2.126174e-06 1.227187e-01 1.227187e-01 -5.724598e-07 -1.644859e-07 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +3.000000e-02 5.912813e-01 5.912813e-01 7.519369e-07 -2.126180e-06 1.227187e-01 1.227187e-01 -5.724564e-07 -1.644793e-07 5.000000e-01 --------------- VOF Barycenter --------------- - time x_vof y_vof vx_vof vy_vof -3.000000e-02 5.000001e-01 6.152841e-02 -4.664813e-06 -1.340349e-06 + time x_vof y_vof vx_vof vy_vof +3.000000e-02 5.000001e-01 6.152842e-02 -4.664785e-06 -1.340296e-06 ********************************************************************************** Transient iteration: 4 Time: 0.04 Time step: 0.01 CFL: 0.000187224 ********************************************************************************** -Pressure drop: -0.00517123 Pa -Total pressure drop: -0.00517123 Pa +Pressure drop: -0.00517128 Pa +Total pressure drop: -0.00517128 Pa ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -4.000000e-02 5.912813e-01 5.912813e-01 7.667839e-07 -2.136658e-06 1.227187e-01 1.227187e-01 -5.803016e-07 -1.666419e-07 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +4.000000e-02 5.912813e-01 5.912813e-01 7.667819e-07 -2.136667e-06 1.227187e-01 1.227187e-01 -5.802978e-07 -1.666305e-07 5.000000e-01 --------------- VOF Barycenter --------------- - time x_vof y_vof vx_vof vy_vof -4.000000e-02 5.000001e-01 6.152842e-02 -4.728715e-06 -1.357918e-06 + time x_vof y_vof vx_vof vy_vof +4.000000e-02 5.000001e-01 6.152844e-02 -4.728684e-06 -1.357825e-06 diff --git a/applications_tests/lethe-fluid/heat_transfer_vof_weakly_compressible_flow.output b/applications_tests/lethe-fluid/heat_transfer_vof_weakly_compressible_flow.output index 406d77f06d..766b17eac2 100644 --- a/applications_tests/lethe-fluid/heat_transfer_vof_weakly_compressible_flow.output +++ b/applications_tests/lethe-fluid/heat_transfer_vof_weakly_compressible_flow.output @@ -16,17 +16,17 @@ Total pressure drop: -156.575 Pa ******************************************************************************* Transient iteration: 2 Time: 0.2 Time step: 0.1 CFL: 0.317772 ******************************************************************************* -Pressure drop: -161.55 Pa -Total pressure drop: -161.55 Pa +Pressure drop: -161.604 Pa +Total pressure drop: -161.604 Pa ******************************************************************************* -Transient iteration: 3 Time: 0.3 Time step: 0.1 CFL: 0.364515 +Transient iteration: 3 Time: 0.3 Time step: 0.1 CFL: 0.363124 ******************************************************************************* -Pressure drop: -162.151 Pa -Total pressure drop: -162.151 Pa +Pressure drop: -162.235 Pa +Total pressure drop: -162.235 Pa ******************************************************************************* -Transient iteration: 4 Time: 0.4 Time step: 0.1 CFL: 0.382727 +Transient iteration: 4 Time: 0.4 Time step: 0.1 CFL: 0.380725 ******************************************************************************* -Pressure drop: -162.279 Pa -Total pressure drop: -162.279 Pa +Pressure drop: -162.369 Pa +Total pressure drop: -162.369 Pa diff --git a/applications_tests/lethe-fluid/time_dependent_boundaries_vof.output b/applications_tests/lethe-fluid/time_dependent_boundaries_vof.output index 6c4973d526..55eab388ef 100644 --- a/applications_tests/lethe-fluid/time_dependent_boundaries_vof.output +++ b/applications_tests/lethe-fluid/time_dependent_boundaries_vof.output @@ -23,27 +23,27 @@ Transient iteration: 1 Time: 0.5 Time step: 0.5 CFL: 2.79578 ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -5.000000e-01 8.788982e-01 8.788982e-02 1.001574e-02 1.280383e-03 1.211018e-01 1.211018e+01 1.784765e+00 1.414586e-01 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +5.000000e-01 8.778856e-01 8.778856e-02 9.871119e-03 1.269337e-03 1.221144e-01 1.221144e+01 1.815053e+00 1.569605e-01 5.000000e-01 --------------- VOF Barycenter --------------- - time x_vof y_vof vx_vof vy_vof -5.000000e-01 4.827262e-01 5.083468e-01 1.473773e-01 1.168096e-02 + time x_vof y_vof vx_vof vy_vof +5.000000e-01 4.800128e-01 5.116058e-01 1.486354e-01 1.285355e-02 ******************************************************************************* -Transient iteration: 2 Time: 1 Time step: 0.5 CFL: 3.72386 +Transient iteration: 2 Time: 1 Time step: 0.5 CFL: 3.72386 ******************************************************************************* ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -1.000000e+00 8.266414e-01 8.266414e-02 8.650876e-03 -2.094834e-04 1.733586e-01 1.733586e+01 7.790878e+00 -2.329163e-01 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +1.000000e+00 8.268896e-01 8.268896e-02 8.653578e-03 -2.851727e-04 1.731104e-01 1.731104e+01 7.785316e+00 -2.389366e-01 5.000000e-01 --------------- VOF Barycenter --------------- - time x_vof y_vof vx_vof vy_vof -1.000000e+00 5.135102e-01 6.388068e-01 4.494082e-01 -1.343552e-02 + time x_vof y_vof vx_vof vy_vof +1.000000e+00 5.135601e-01 6.386556e-01 4.497314e-01 -1.380256e-02 diff --git a/applications_tests/lethe-fluid/vof-velocity-extrapolation.output b/applications_tests/lethe-fluid/vof-velocity-extrapolation.output index 8d3311cd18..acf236a20d 100644 --- a/applications_tests/lethe-fluid/vof-velocity-extrapolation.output +++ b/applications_tests/lethe-fluid/vof-velocity-extrapolation.output @@ -38,7 +38,7 @@ Transient iteration: 1 Time: 0.04 Time step: 0.04 CFL: 0 VOF Mass Conservation ---------------------- time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -4.000000e-02 1.998729e-01 1.998729e-04 -3.160239e-14 1.465966e-08 2.000127e+00 2.000127e+00 -2.931234e-11 -1.733959e-05 5.000000e-01 +4.000000e-02 1.998729e-01 1.998729e-04 -3.141666e-14 1.465966e-08 2.000127e+00 2.000127e+00 -2.941613e-11 -1.733959e-05 5.000000e-01 ******************************************************************************** Transient iteration: 2 Time: 0.106 Time step: 0.066 CFL: 0.0369272 @@ -47,195 +47,195 @@ Transient iteration: 2 Time: 0.106 Time step: 0.066 CFL: 0.0369272 ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -1.060000e-01 1.998956e-01 1.998956e-04 -2.432663e-13 3.812792e-08 2.000104e+00 2.000104e+00 -9.809253e-12 -4.094545e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +1.060000e-01 1.998954e-01 1.998954e-04 -2.494595e-13 3.809632e-08 2.000105e+00 2.000105e+00 -9.415353e-12 -4.091637e-05 5.000000e-01 ******************************************************************************* -Transient iteration: 3 Time: 0.216 Time step: 0.11 CFL: 0.145606 +Transient iteration: 3 Time: 0.216 Time step: 0.11 CFL: 0.145625 ******************************************************************************* ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -2.160000e-01 1.999246e-01 1.999246e-04 -1.097792e-13 6.931116e-08 2.000075e+00 2.000075e+00 -1.386806e-11 -7.210426e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +2.160000e-01 1.999235e-01 1.999235e-04 -1.199760e-13 6.879648e-08 2.000077e+00 2.000077e+00 -1.135472e-11 -7.161852e-05 5.000000e-01 ******************************************************************************* -Transient iteration: 4 Time: 0.319419 Time step: 0.103419 CFL: 0.425454 +Transient iteration: 4 Time: 0.319668 Time step: 0.103668 CFL: 0.42443 ******************************************************************************* ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -3.194190e-01 1.999461e-01 1.999461e-04 -1.827551e-13 8.720432e-08 2.000054e+00 2.000054e+00 -3.803980e-11 -9.007122e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +3.196684e-01 1.999441e-01 1.999441e-04 -7.846626e-14 8.547453e-08 2.000056e+00 2.000056e+00 -2.258717e-11 -8.814305e-05 5.000000e-01 ******************************************************************************** -Transient iteration: 5 Time: 0.40004 Time step: 0.0806207 CFL: 0.513114 +Transient iteration: 5 Time: 0.401429 Time step: 0.0817608 CFL: 0.507179 ******************************************************************************** ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -4.000397e-01 1.999615e-01 1.999615e-04 -9.479269e-14 9.162297e-08 2.000039e+00 2.000039e+00 -8.765327e-11 -9.414623e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +4.014292e-01 1.999579e-01 1.999579e-04 -6.338498e-14 9.027331e-08 2.000042e+00 2.000042e+00 -7.169848e-11 -9.275883e-05 5.000000e-01 ******************************************************************************** -Transient iteration: 6 Time: 0.472086 Time step: 0.0720464 CFL: 0.447604 +Transient iteration: 6 Time: 0.474477 Time step: 0.0730481 CFL: 0.44771 ******************************************************************************** ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -4.720862e-01 1.999731e-01 1.999731e-04 -2.116659e-13 9.008146e-08 2.000027e+00 2.000027e+00 -7.869427e-11 -9.237840e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +4.744773e-01 1.999686e-01 1.999686e-04 -1.698121e-13 8.913680e-08 2.000031e+00 2.000031e+00 -6.363754e-11 -9.136797e-05 5.000000e-01 ******************************************************************************** -Transient iteration: 7 Time: 0.539958 Time step: 0.0678721 CFL: 0.424601 +Transient iteration: 7 Time: 0.542935 Time step: 0.0684574 CFL: 0.426824 ******************************************************************************** ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -5.399583e-01 1.999826e-01 1.999826e-04 -2.509053e-13 8.506966e-08 2.000017e+00 2.000017e+00 -4.779690e-11 -8.744895e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +5.429348e-01 1.999770e-01 1.999770e-04 -2.264778e-13 8.420472e-08 2.000023e+00 2.000023e+00 -3.004115e-11 -8.665868e-05 5.000000e-01 ******************************************************************************** -Transient iteration: 8 Time: 0.605071 Time step: 0.0651127 CFL: 0.416952 +Transient iteration: 8 Time: 0.608582 Time step: 0.0656473 CFL: 0.417122 ******************************************************************************** ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -6.050710e-01 1.999920e-01 1.999920e-04 -2.905018e-13 7.709414e-08 2.000008e+00 2.000008e+00 -2.100524e-11 -7.958567e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +6.085821e-01 1.999857e-01 1.999857e-04 -2.624844e-13 7.635390e-08 2.000014e+00 2.000014e+00 -6.402991e-13 -7.886611e-05 5.000000e-01 ******************************************************************************** -Transient iteration: 9 Time: 0.668539 Time step: 0.0634678 CFL: 0.410367 +Transient iteration: 9 Time: 0.672579 Time step: 0.0639966 CFL: 0.410318 ******************************************************************************** ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -6.685388e-01 2.000016e-01 2.000016e-04 -3.130787e-13 6.596255e-08 1.999998e+00 1.999998e+00 -7.135564e-12 -6.848673e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +6.725787e-01 1.999945e-01 1.999945e-04 -2.437133e-13 6.531042e-08 2.000006e+00 2.000006e+00 1.039594e-10 -6.776133e-05 5.000000e-01 ******************************************************************************** -Transient iteration: 10 Time: 0.731211 Time step: 0.0626719 CFL: 0.40508 +Transient iteration: 10 Time: 0.735885 Time step: 0.0633065 CFL: 0.40436 ******************************************************************************** ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -7.312107e-01 2.000112e-01 2.000112e-04 -2.750099e-13 5.156979e-08 1.999989e+00 1.999989e+00 6.511611e-11 -5.392560e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +7.358852e-01 2.000031e-01 2.000031e-04 -2.076235e-13 5.092551e-08 1.999997e+00 1.999997e+00 1.893582e-10 -5.303707e-05 5.000000e-01 ******************************************************************************** -Transient iteration: 11 Time: 0.79412 Time step: 0.0629094 CFL: 0.39849 +Transient iteration: 11 Time: 0.799503 Time step: 0.0636178 CFL: 0.398043 ******************************************************************************** ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -7.941201e-01 2.000204e-01 2.000204e-04 -2.568808e-13 3.423597e-08 1.999980e+00 1.999980e+00 4.284320e-11 -3.614549e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +7.995030e-01 2.000103e-01 2.000103e-04 -1.981145e-13 3.369315e-08 1.999990e+00 1.999990e+00 1.892086e-10 -3.550299e-05 5.000000e-01 ******************************************************************************** -Transient iteration: 12 Time: 0.858447 Time step: 0.0643273 CFL: 0.391183 +Transient iteration: 12 Time: 0.864635 Time step: 0.0651318 CFL: 0.390702 ******************************************************************************** ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -8.584473e-01 2.000282e-01 2.000282e-04 -3.299305e-13 1.534834e-08 1.999972e+00 1.999972e+00 -1.635443e-11 -1.712386e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +8.646347e-01 2.000168e-01 2.000168e-04 -1.798912e-13 1.497962e-08 1.999983e+00 1.999983e+00 1.538475e-10 -1.689141e-05 5.000000e-01 ******************************************************************************** -Transient iteration: 13 Time: 0.924811 Time step: 0.0663637 CFL: 0.387725 +Transient iteration: 13 Time: 0.931897 Time step: 0.0672624 CFL: 0.387329 ******************************************************************************** ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -9.248111e-01 2.000363e-01 2.000363e-04 -3.060822e-13 -3.679444e-09 1.999964e+00 1.999964e+00 -6.483153e-12 1.734794e-06 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +9.318972e-01 2.000242e-01 2.000242e-04 -1.501923e-13 -3.712595e-09 1.999976e+00 1.999976e+00 1.579850e-10 1.606456e-06 5.000000e-01 ******************************************************************************** -Transient iteration: 14 Time: 0.993555 Time step: 0.0687442 CFL: 0.386149 +Transient iteration: 14 Time: 1.00169 Time step: 0.0697911 CFL: 0.385507 ******************************************************************************** ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -9.935553e-01 2.000454e-01 2.000454e-04 -3.181934e-13 -2.241224e-08 1.999955e+00 1.999955e+00 2.134781e-11 2.033038e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +1.001688e+00 2.000326e-01 2.000326e-04 -2.810896e-13 -2.182440e-08 1.999967e+00 1.999967e+00 1.733141e-10 1.985906e-05 5.000000e-01 ******************************************************************************** -Transient iteration: 15 Time: 1.06493 Time step: 0.0713699 CFL: 0.385284 +Transient iteration: 15 Time: 1.0744 Time step: 0.0727146 CFL: 0.383918 ******************************************************************************** ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -1.064925e+00 2.000552e-01 2.000552e-04 -3.110473e-13 -3.932277e-08 1.999945e+00 1.999945e+00 1.012938e-10 3.743795e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +1.074403e+00 2.000409e-01 2.000409e-04 -2.397834e-13 -3.842717e-08 1.999959e+00 1.999959e+00 2.485756e-10 3.666188e-05 5.000000e-01 ******************************************************************************** -Transient iteration: 16 Time: 1.14021 Time step: 0.0752821 CFL: 0.379213 +Transient iteration: 16 Time: 1.15149 Time step: 0.0770904 CFL: 0.377295 ******************************************************************************** ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -1.140207e+00 2.000645e-01 2.000645e-04 -2.573632e-13 -5.301125e-08 1.999935e+00 1.999935e+00 1.573958e-10 5.120977e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +1.151493e+00 2.000486e-01 2.000486e-04 -1.802561e-13 -5.176504e-08 1.999951e+00 1.999951e+00 3.045174e-10 5.000552e-05 5.000000e-01 ******************************************************************************** -Transient iteration: 17 Time: 1.22222 Time step: 0.0820177 CFL: 0.367151 +Transient iteration: 17 Time: 1.23601 Time step: 0.0845188 CFL: 0.364844 ******************************************************************************** ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -1.222225e+00 2.000741e-01 2.000741e-04 -2.561767e-13 -6.255171e-08 1.999926e+00 1.999926e+00 2.608132e-10 6.068088e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +1.236012e+00 2.000568e-01 2.000568e-04 -1.967321e-13 -6.087420e-08 1.999943e+00 1.999943e+00 3.707246e-10 5.896846e-05 5.000000e-01 ******************************************************************************** -Transient iteration: 18 Time: 1.31244 Time step: 0.0902194 CFL: 0.351978 +Transient iteration: 18 Time: 1.32898 Time step: 0.0929707 CFL: 0.348097 ******************************************************************************** ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -1.312444e+00 2.000851e-01 2.000851e-04 -2.376066e-13 -6.701049e-08 1.999915e+00 1.999915e+00 3.001884e-10 6.500908e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +1.328983e+00 2.000665e-01 2.000665e-04 -1.605929e-13 -6.453797e-08 1.999933e+00 1.999933e+00 4.040135e-10 6.248128e-05 5.000000e-01 -******************************************************************************** -Transient iteration: 19 Time: 1.41169 Time step: 0.0992414 CFL: 0.329891 -******************************************************************************** +******************************************************************************* +Transient iteration: 19 Time: 1.43125 Time step: 0.102268 CFL: 0.328105 +******************************************************************************* ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -1.411686e+00 2.000975e-01 2.000975e-04 -1.986559e-13 -6.580788e-08 1.999902e+00 1.999902e+00 3.284613e-10 6.366724e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +1.431250e+00 2.000778e-01 2.000778e-04 -1.130150e-13 -6.183007e-08 1.999922e+00 1.999922e+00 4.318061e-10 5.976764e-05 5.000000e-01 ******************************************************************************* -Transient iteration: 20 Time: 1.52085 Time step: 0.109165 CFL: 0.309595 +Transient iteration: 20 Time: 1.54374 Time step: 0.112495 CFL: 0.301612 ******************************************************************************* ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -1.520851e+00 2.001114e-01 2.001114e-04 -9.176963e-14 -5.713112e-08 1.999889e+00 1.999889e+00 3.629224e-10 5.499841e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +1.543745e+00 2.000899e-01 2.000899e-04 1.193763e-14 -5.135689e-08 1.999910e+00 1.999910e+00 4.632057e-10 4.930144e-05 5.000000e-01 ******************************************************************************* -Transient iteration: 21 Time: 1.64093 Time step: 0.120082 CFL: 0.262858 +Transient iteration: 21 Time: 1.66749 Time step: 0.123744 CFL: 0.245207 ******************************************************************************* ---------------------- VOF Mass Conservation ---------------------- - time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold -1.640933e+00 2.001256e-01 2.001256e-04 2.384551e-15 -3.763005e-08 1.999874e+00 1.999874e+00 4.077472e-10 3.542519e-05 5.000000e-01 + time surface_fluid_0 mass_per_length_fluid_0 momentum-x_fluid_0 momentum-y_fluid_0 surface_fluid_1 mass_per_length_fluid_1 momentum-x_fluid_1 momentum-y_fluid_1 sharpening_threshold +1.667489e+00 2.001026e-01 2.001026e-04 2.266145e-14 -3.069071e-08 1.999897e+00 1.999897e+00 5.092479e-10 2.862674e-05 5.000000e-01 diff --git a/doc/source/examples/multiphysics/3d-dam-break/3d-dam-break.rst b/doc/source/examples/multiphysics/3d-dam-break/3d-dam-break.rst index 86f4bd8215..0d533e51dc 100644 --- a/doc/source/examples/multiphysics/3d-dam-break/3d-dam-break.rst +++ b/doc/source/examples/multiphysics/3d-dam-break/3d-dam-break.rst @@ -90,7 +90,7 @@ Parameter File Simulation Control ~~~~~~~~~~~~~~~~~~ -Time integration is handled by a 1st order backward differentiation scheme (`bdf1`), for a :math:`6 \ \text{s}` simulation time with an initial time step of :math:`0.001 \ \text{s}`. Time-step adaptation is enabled using ``adapt=true`` +Time integration is handled by a 1st order backward differentiation scheme (`bdf1`), for a :math:`6 \ \text{s}` simulation time with an initial time-step of :math:`0.001 \ \text{s}`. Time-step adaptation is enabled using ``adapt=true`` and the max CFL is :math:`0.8`. .. note:: @@ -235,7 +235,11 @@ Running the Simulation We call the lethe-fluid by invoking: -``mpirun -np $number_of_CPU_cores lethe-fluid 3d-dam-break.prm`` + +.. code-block:: text + :class: copy-button + + mpirun -np $number_of_CPU_cores lethe-fluid 3d-dam-break.prm .. warning:: Make sure to compile Lethe in `Release` mode and run in parallel using ``mpirun``. This simulation took :math:`\approx` 17 hours on 64 processes (runned on the `Narval `_ cluster). diff --git a/doc/source/examples/multiphysics/air-bubble-compression/air-bubble-compression.rst b/doc/source/examples/multiphysics/air-bubble-compression/air-bubble-compression.rst index f0fdeb5b9c..2d1a936c78 100644 --- a/doc/source/examples/multiphysics/air-bubble-compression/air-bubble-compression.rst +++ b/doc/source/examples/multiphysics/air-bubble-compression/air-bubble-compression.rst @@ -5,6 +5,7 @@ Air Bubble Compression This example simulates the compression of an air bubble by surrounding liquid. The problem is inspired by the test case of Caltagirone *et al.* [#caltagirone2011]_ +**** -------- Features @@ -16,6 +17,7 @@ Features - Unsteady problem handled by an adaptive BDF2 time-stepping scheme - Usage of a python script for post-processing data +**** --------------------------- Files Used in This Example @@ -26,6 +28,7 @@ Both files mentioned below are located in the example's folder (``examples/multi - Parameter file: ``air-bubble-compression.prm`` - Postprocessing python script: ``air-bubble-compression-postprocessing.py`` +**** ----------------------- Description of the Case @@ -47,6 +50,8 @@ The initial configuration of this example is illustrated below. .. note:: In this example, gravity and surface tension forces are not considered. +**** + -------------- Parameter File -------------- @@ -54,8 +59,8 @@ Parameter File Simulation Control ~~~~~~~~~~~~~~~~~~ -Time integration is handled by a 2nd-order backward differentiation scheme (``bdf2``) with a variable time step. -The initial time step is set to :math:`0.1 \, \text{s}` and the simulation lasts :math:`2.5 \, \text{s}`. +Time integration is handled by a 2nd-order backward differentiation scheme (``bdf2``) with a variable time-step. +The initial time-step is set to :math:`0.1 \, \text{s}` and the simulation lasts :math:`2.5 \, \text{s}`. .. code-block:: text @@ -88,7 +93,7 @@ In the ``VOF`` subsection, the ``compressible``, the ``interface sharpening``, a The enabled ``compressible`` parameter allows interface compression by adding the term :math:`\phi (\nabla \cdot \mathbf{u})` to the VOF equation. The ``interface sharpening`` method and its parameters are explained in the :doc:`../dam-break/dam-break` example. The ``phase filtration`` filters the phase field used for the calculation of physical properties by stiffening the value of the phase fraction. -We refer the reader to :doc:`../../../../theory/multiphase/cfd/vof` theory guide for further explanation on the ``phase filtration``. +We refer the reader to :doc:`../../../../theory/multiphase/cfd/vof` theory guide for further explanation on ``phase filtration``. .. code-block:: text @@ -128,7 +133,7 @@ Boundary Conditions ~~~~~~~~~~~~~~~~~~~ On all four sides of the domain, water which is associated with the phase fraction :math:`\phi=1` is injected. -This is done in the simulation by setting the velocities of the fluid in the ``boundary conditions`` subsection and by selecting the correct fluid in the ``boundary conditions VOF`` subsection with a ``dirichlet`` boundary condition on the phase fraction as shown below. +This is done in the simulation by setting the velocities of the fluid in the ``boundary conditions`` subsection and by selecting the correct fluid in the ``boundary conditions VOF`` subsection using ``dirichlet`` boundary conditions on the phase fraction as shown below. Boundary Conditions - Fluid Dynamics ************************************ @@ -263,6 +268,7 @@ In the ``mesh adaptation`` subsection, adaptive mesh refinement is defined for t set initial refinement steps = 6 end +**** ----------------------- Running the Simulation @@ -280,6 +286,7 @@ to run the simulation using eight CPU cores. Feel free to use more. .. warning:: Make sure to compile lethe in `Release` mode and run in parallel using mpirun. This simulation takes :math:`\sim` 1.5 minute on 8 processes. +**** ------- Results @@ -309,7 +316,7 @@ The results can be post-processed by invoking the following command from the fol python3 air-bubble-compression-postprocessing.py . air-bubble-compression.prm .. important:: - You need to ensure that ``lethe_pyvista_tools`` is working on your machine. Click `here <../../../tools/postprocessing/postprocessing.html>`_ for details. + You need to ensure that ``lethe_pyvista_tools`` is working on your machine. Click :doc:`here <../../../tools/postprocessing/postprocessing_pyvista>` for details. The following figures present the comparison between the analytical results and the simulation results for the density and pressure evolutions evaluated at the center of the bubble. A pretty good agreement between the simulation and analytical results is observed. @@ -332,6 +339,7 @@ The following figures present the comparison between the analytical results and | | +-------------------------------------------------------------------------------------------------------------------+ +**** ---------- References diff --git a/doc/source/examples/multiphysics/capillary-wave/capillary-wave.rst b/doc/source/examples/multiphysics/capillary-wave/capillary-wave.rst index 3a6a5f61f3..c6113ff8bc 100644 --- a/doc/source/examples/multiphysics/capillary-wave/capillary-wave.rst +++ b/doc/source/examples/multiphysics/capillary-wave/capillary-wave.rst @@ -4,6 +4,7 @@ Capillary Wave This example simulates the damping of a small amplitude capillary wave for different time-steps allowing us to study the capillary time-step constraint. The problem is inspired by the test case of Denner *et al.* [#denner2022]_ +**** -------- Features @@ -15,6 +16,7 @@ Features - Bash scripts to write, launch, and postprocess multiple cases - Python scripts for postprocessing data +**** --------------------------- Files Used in This Example @@ -32,6 +34,7 @@ All files mentioned below are located in the example's folder (``examples/multip - Postprocessing Bash script: ``capillary-wave-time-step-sensitivity-postprocess.sh`` - Python script for computing quantities of interest: ``capillary-wave-calculation.py`` +**** ----------------------- Description of the Case @@ -88,6 +91,7 @@ with the shortest unambiguously resolved capillary wave having a wavelength of : Therefore, in order to get stable simulation results, :math:`\Delta t < \Delta t_\sigma` should be respected. In this example, different time-steps will be used to explore the stability limit of Lethe's current implementation. +**** -------------- Parameter File @@ -204,12 +208,13 @@ In the ``physical properties`` subsection, we define the fluids such that both f end end +**** ----------------------- Running the Simulation ----------------------- -We can call ``lethe-fluid`` for each time step value. For :math:`\Delta t \approx 0.95\Delta t_\sigma`, this can be done by invoking the following command: +We can call ``lethe-fluid`` for each time-step value. For :math:`\Delta t \approx 0.95\Delta t_\sigma`, this can be done by invoking the following command: .. code-block:: text :class: copy-button @@ -253,6 +258,7 @@ to run the simulation using four CPU cores. Feel free to use more CPU cores. This script runs the ``capillary-wave-calculation.py`` script before generating the different cases. Make sure that the information entered in the Python script corresponds to the ones you wish to simulate. +**** ------- Results @@ -361,6 +367,7 @@ By increasing the mesh resolution by an additional refinement, the :math:`\mathr | | +-----------------------------------------------------------------------------------------------------------------------------+ +**** ---------------- Acknowledgment @@ -368,6 +375,7 @@ Acknowledgment We would like to thank Prof. Fabian Denner for sharing his time and knowledge throughout the process of developing this example. +**** ---------- References diff --git a/doc/source/examples/multiphysics/capillary-wave/images/TSM_comparison_figure.png b/doc/source/examples/multiphysics/capillary-wave/images/TSM_comparison_figure.png index 1a010b637a..0942daac28 100644 Binary files a/doc/source/examples/multiphysics/capillary-wave/images/TSM_comparison_figure.png and b/doc/source/examples/multiphysics/capillary-wave/images/TSM_comparison_figure.png differ diff --git a/doc/source/examples/multiphysics/capillary-wave/images/figure-TSM-0_95.png b/doc/source/examples/multiphysics/capillary-wave/images/figure-TSM-0_95.png index 3eb4ee15a1..adb7f03be2 100644 Binary files a/doc/source/examples/multiphysics/capillary-wave/images/figure-TSM-0_95.png and b/doc/source/examples/multiphysics/capillary-wave/images/figure-TSM-0_95.png differ diff --git a/doc/source/examples/multiphysics/capillary-wave/images/initial-state.svg b/doc/source/examples/multiphysics/capillary-wave/images/initial-state.svg index c129b03e33..1d4a26974b 100644 --- a/doc/source/examples/multiphysics/capillary-wave/images/initial-state.svg +++ b/doc/source/examples/multiphysics/capillary-wave/images/initial-state.svg @@ -1,1619 +1,225 @@ - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + \ No newline at end of file diff --git a/doc/source/examples/multiphysics/dam-break/dam-break.rst b/doc/source/examples/multiphysics/dam-break/dam-break.rst index 634eb2df67..9a7e0df05b 100644 --- a/doc/source/examples/multiphysics/dam-break/dam-break.rst +++ b/doc/source/examples/multiphysics/dam-break/dam-break.rst @@ -4,6 +4,7 @@ Dam-Break This example simulates the dam break experiments of Martin and Moyce [#martin1952]_. +**** ---------------------------------- Features @@ -14,6 +15,7 @@ Features - Unsteady problem handled by an adaptive BDF1 time-stepping scheme - The use of a python script for post-processing data +**** --------------------------- Files Used in This Example @@ -23,6 +25,7 @@ Both files mentioned below are located in the example's folder (``examples/multi - Parameter file: ``dam-break-Martin-and-Moyce.prm`` - Postprocessing Python script: ``dam-break-2d.py`` +**** --------------------------- Description of the Case @@ -47,6 +50,7 @@ The following schematic describes the geometry and dimensions of the simulation All the four boundary conditions are ``slip``, and an external gravity field of :math:`-1` is applied in the :math:`y` direction. +**** -------------- Parameter File @@ -57,7 +61,7 @@ Simulation Control Time integration is handled by a 1st order backward differentiation scheme `(bdf1)`, for a :math:`4.1` s simulation time with an initial -time step of :math:`0.01` seconds. +time-step of :math:`0.01` seconds. .. note:: This example uses an adaptive time-stepping method, where the @@ -88,14 +92,10 @@ and off `(false)` the physics of interest. Here ``VOF`` is chosen. subsection multiphysics set VOF = true - end - - -Interface Sharpening Parameters -~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + end VOF -*** +~~~ If the ``interface sharpening`` is not enabled in :doc:`VOF <../../../parameters/cfd/volume_of_fluid>` subsection, the interface between phases will become blurry (due to diffusion). Furthermore, the ``phase filtration`` is enabled in this example. We refer the reader to the :doc:`../../../theory/multiphase/cfd/vof` documentation for more explanation on both methods. @@ -114,11 +114,8 @@ If the ``interface sharpening`` is not enabled in :doc:`VOF <../../../parameters end end -Fluid Phase Parameters -~~~~~~~~~~~~~~~~~~~~~~ - Initial Conditions -****************** +~~~~~~~~~~~~~~~~~~ In the ``initial conditions``, the initial velocity and initial position of the liquid phase are defined. The liquid phase is initially @@ -137,7 +134,7 @@ defined as rectangle of length :math:`= 3.5` and height :math:`= 7`. end Source Term -*********** +~~~~~~~~~~~ The ``source term`` subsection defines the gravitational acceleration: @@ -150,7 +147,7 @@ The ``source term`` subsection defines the gravitational acceleration: end Physical Properties -******************* +~~~~~~~~~~~~~~~~~~~ Two fluids are present in this simulation, hence in the ``physical properties`` subsection, their physical properties should be defined: @@ -214,6 +211,7 @@ is adapted to the initial condition for the phase. set initial refinement steps = 4 end +**** ---------------------- Running the Simulation @@ -234,6 +232,7 @@ to run the simulation using six CPU cores. Feel free to use more. Make sure to compile Lethe in `Release` mode and run in parallel using mpirun. This simulation takes :math:`\sim \, 3` minutes on :math:`6` processes. +**** ----------------------- Results and Discussion @@ -247,8 +246,8 @@ The red area corresponds to the liquid phase and the blue area corresponds to th :alt: time-shots :align: center -A python post-processing code `(dam-break-2d.py)` -is added to the example folder to post-process the results. +A python postprocessing code `(dam-break-2d.py)` +is added to the example folder to postprocess the results. Run .. code-block:: text @@ -256,10 +255,10 @@ Run python3 ./dam-break-2d.py ./output -to execute this post-processing code, where ``./output`` is the directory that contains the simulation results. -In post-processing, the maximum dimensionless lateral position of the liquid phase is tracked +to execute this postprocessing code, where ``./output`` is the directory that contains the simulation results. +In postprocessing, the maximum dimensionless lateral position of the liquid phase is tracked through time and compared with the experiments of Martin and Moyce (1952) [#martin1952]_. -The following figure shows the result of the post-processing, with a good agreement between the simulation and the experiment: +The following figure shows the result of the postprocessing, with a good agreement between the simulation and the experiment: .. image:: images/xmax-t.png :alt: xmax_t @@ -275,6 +274,7 @@ and refines the meshes on the interface. :alt: refinement :align: center +**** ---------------------------- References diff --git a/doc/source/examples/multiphysics/dam-break/images/dam-break-initial-configuration.svg b/doc/source/examples/multiphysics/dam-break/images/dam-break-initial-configuration.svg index 5c6ee9a76e..e7a480a8f1 100644 --- a/doc/source/examples/multiphysics/dam-break/images/dam-break-initial-configuration.svg +++ b/doc/source/examples/multiphysics/dam-break/images/dam-break-initial-configuration.svg @@ -1,589 +1,206 @@ - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + \ No newline at end of file diff --git a/doc/source/examples/multiphysics/dam-break/images/refinement.png b/doc/source/examples/multiphysics/dam-break/images/refinement.png index 44ec3d0153..3b8bbb0065 100644 Binary files a/doc/source/examples/multiphysics/dam-break/images/refinement.png and b/doc/source/examples/multiphysics/dam-break/images/refinement.png differ diff --git a/doc/source/examples/multiphysics/dam-break/images/time-series.png b/doc/source/examples/multiphysics/dam-break/images/time-series.png index cb073f0066..7796212ff7 100644 Binary files a/doc/source/examples/multiphysics/dam-break/images/time-series.png and b/doc/source/examples/multiphysics/dam-break/images/time-series.png differ diff --git a/doc/source/examples/multiphysics/rayleigh-taylor-instability/images/adaptive_mass.png b/doc/source/examples/multiphysics/rayleigh-taylor-instability/images/adaptive_mass.png index 018923cf3d..8c184d280b 100644 Binary files a/doc/source/examples/multiphysics/rayleigh-taylor-instability/images/adaptive_mass.png and b/doc/source/examples/multiphysics/rayleigh-taylor-instability/images/adaptive_mass.png differ diff --git a/doc/source/examples/multiphysics/rayleigh-taylor-instability/images/constant_mass.png b/doc/source/examples/multiphysics/rayleigh-taylor-instability/images/constant_mass.png index 32e76d29c7..e4c236d2eb 100644 Binary files a/doc/source/examples/multiphysics/rayleigh-taylor-instability/images/constant_mass.png and b/doc/source/examples/multiphysics/rayleigh-taylor-instability/images/constant_mass.png differ diff --git a/doc/source/examples/multiphysics/rayleigh-taylor-instability/images/rayleigh-taylor-instability-initial-state.svg b/doc/source/examples/multiphysics/rayleigh-taylor-instability/images/rayleigh-taylor-instability-initial-state.svg index e70e4a5e02..e43d144044 100644 --- a/doc/source/examples/multiphysics/rayleigh-taylor-instability/images/rayleigh-taylor-instability-initial-state.svg +++ b/doc/source/examples/multiphysics/rayleigh-taylor-instability/images/rayleigh-taylor-instability-initial-state.svg @@ -1,685 +1,239 @@ - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Fluid 1 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Fluid 0 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + \ No newline at end of file diff --git a/doc/source/examples/multiphysics/rayleigh-taylor-instability/rayleigh-taylor-instability.rst b/doc/source/examples/multiphysics/rayleigh-taylor-instability/rayleigh-taylor-instability.rst index 5894a2a0ad..7c587d225a 100644 --- a/doc/source/examples/multiphysics/rayleigh-taylor-instability/rayleigh-taylor-instability.rst +++ b/doc/source/examples/multiphysics/rayleigh-taylor-instability/rayleigh-taylor-instability.rst @@ -4,6 +4,7 @@ Rayleigh-Taylor Instability This example simulates the dynamic evolution of the single-mode Rayleigh-Taylor instability [#he1999]_ by density contrast. +**** -------- Features @@ -16,6 +17,7 @@ Features - Monitoring mass conservation - Interface sharpening +**** -------------------------- Files Used in This Example @@ -26,6 +28,7 @@ All files mentioned below are located in the example's folder (``examples/multip - Parameter file for constant sharpening: ``rayleigh-taylor-instability-constant-sharpening.prm`` - Postprocessing Python script: ``rayleigh-taylor_postprocess.py`` +**** ----------------------- Description of the Case @@ -56,8 +59,9 @@ which result in Reynolds and Atwood numbers equal to At = \frac{\rho_r - 1}{\rho_r + 1} = 0.5 -A perturbed interface defined as :math:`y = 2H + 0.1 H \cos{(2 \pi x / H)}` separates the fluids. At the top and bottom boundaries, a ``no-slip`` boundary condition is applied, while on the left and right walls, a ``periodic`` boundary condition is used. The temporal evolution of the interface is visually compared with the simulations of Garoosi and Hooman [#garoosi2022]_ at dimensionless times (:math:`t^* = t \sqrt{\mathbf{g} / H}`) of :math:`1.5`, :math:`2.5`, :math:`3.5`, :math:`4.0` and :math:`4.5`. The temporal evolution of the spike and the bubble positions are then compared to the results of He *et al.* [#he1999]_ The term "spike" refers to the lowest point of ``fluid 1`` and the term "bubble" refers to the highest point of ``fluid 0``. +A perturbed interface defined as :math:`y = 2H + 0.1 H \cos{(2 \pi x / H)}` separates the fluids. At the top and bottom boundaries, a ``no-slip`` boundary condition is applied, while on the left and right walls, ``slip`` boundary conditions are used. The temporal evolution of the interface is visually compared with the simulations of Garoosi and Hooman [#garoosi2022]_ at dimensionless times (:math:`t^* = t \sqrt{\mathbf{g} / H}`) of :math:`1.5`, :math:`2.5`, :math:`3.5`, :math:`4.0` and :math:`4.5`. The temporal evolution of the spike and the bubble positions are then compared to the results of He *et al.* [#he1999]_ The term "spike" refers to the lowest point of ``fluid 1`` and the term "bubble" refers to the highest point of ``fluid 0``. +**** -------------- Parameter File @@ -68,7 +72,7 @@ Simulation Control Time integration is handled by a 2nd order backward differentiation scheme (``bdf2``), for a :math:`0.75\, \text{s}` simulation time with an initial -time step of :math:`0.0002` seconds. Time-step adaptation is enabled using ``adapt = true`` +time-step of :math:`0.0002` seconds. Time-step adaptation is enabled using ``adapt = true`` and the max CFL is :math:`0.8`. .. note:: @@ -192,23 +196,25 @@ We set ``initial refinement steps = 4`` to adapt the mesh to the initial value o Boundary Conditions ~~~~~~~~~~~~~~~~~~~ -The boundary conditions applied on the left and right boundaries are ``periodic``, while a ``noslip`` boundary condition is used for the top and bottom walls. In the definition of a ``periodic`` boundary, we need to specify the ``periodic_id`` and the ``periodic_direction`` (in this example, :math:`0` which indicates the :math:`x` direction). +The boundary conditions applied on the left and right boundaries are ``slip``, while a ``noslip`` boundary condition is used for the top and bottom walls. .. code-block:: text subsection boundary conditions - set number = 3 + set number = 4 subsection bc 0 - set id = 0 - set type = periodic - set periodic_id = 1 - set periodic_direction = 0 + set id = 0 + set type = slip end subsection bc 1 + set id = 1 + set type = slip + end + subsection bc 2 set id = 2 set type = noslip end - subsection bc 2 + subsection bc 3 set id = 3 set type = noslip end @@ -274,6 +280,8 @@ In the ``post-processing`` subsection, the output of the mass of each fluid is e set calculate mass conservation = true end +**** + --------------------------- Running the Simulation --------------------------- @@ -293,6 +301,7 @@ to run the simulations using eight CPU cores. Feel free to use more. run in parallel using mpirun. On :math:`8` processes, this simulation takes :math:`\sim` :math:`2` minutes for the ``adaptive`` sharpening and :math:`\sim` :math:`4` minutes for ``constant`` sharpening. +**** ----------------------- Results and Discussion @@ -365,6 +374,8 @@ The following figures show the mass of ``fluid 1`` throughout the simulation wit | | | +---------------------------------------------+---------------------------------------------+ +**** + ----------- References ----------- diff --git a/doc/source/examples/multiphysics/sloshing-in-rectangular-tank/sloshing-in-rectangular-tank.rst b/doc/source/examples/multiphysics/sloshing-in-rectangular-tank/sloshing-in-rectangular-tank.rst index f5fadfe040..606c4f9eb0 100644 --- a/doc/source/examples/multiphysics/sloshing-in-rectangular-tank/sloshing-in-rectangular-tank.rst +++ b/doc/source/examples/multiphysics/sloshing-in-rectangular-tank/sloshing-in-rectangular-tank.rst @@ -4,6 +4,7 @@ Sloshing in a Rectangular Tank This example simulates the damping of a small amplitude wave for Reynolds number of (:math:`2`, :math:`20`, :math:`200` and :math:`2000`). The problem is inspired by the test case of Carrica *et al.* [#carrica2007]_ +**** -------- Features @@ -11,9 +12,10 @@ Features - Solver: ``lethe-fluid`` - Volume of fluid (VOF) -- Unsteady problem handled by an adaptive BDF2 time-stepping scheme +- Unsteady problem handled by an adaptive BDF2 time-stepping scheme - Usage of a python script for post-processing data +**** --------------------------- Files Used in This Example @@ -29,12 +31,13 @@ All files mentioned below are located in the example's folder (``examples/multip - Parameter file for :math:`Re=2000`: ``sloshing_2000/sloshing-in-rectangular-tank_Re2000.prm`` - Postprocessing Python script: ``sloshing_postprocessing.py`` +**** ----------------------- Description of the Case ----------------------- -Predicting the dynamics of free surface waves is essential for many industrial applications (e.g. transport of liquified natural gas). +Predicting the dynamics of free surface waves is essential for many industrial applications (e.g. transport of liquefied natural gas). Yet, simulating their dynamics is difficult, especially for high Reynolds number values. Indeed, in this case, the amplitude of the waves dampen very slowly. This leads to an oscillatory wave problem which is highly sensitive to the time integration scheme and the coupling between the VOF solver and the Navier-Stokes solver. @@ -56,6 +59,7 @@ In this problem, we simulate the damping of a small amplitude wave in a rectangu Four values of the Reynolds number are investigated: :math:`2`, :math:`20`, :math:`200` and :math:`2000`. +**** -------------- Parameter File @@ -64,7 +68,7 @@ Parameter File Simulation Control ~~~~~~~~~~~~~~~~~~ -The results for this problem are highly sensitive to the accuracy of the time-stepping scheme. For this reason, we use a 2nd order backward differentiation scheme (``bdf2``) with a variable time step. The ``adaptive time step scaling`` is set to :math:`1.025` to ensure that the time-step does not rise too quickly during wave oscillations. +The results for this problem are highly sensitive to the accuracy of the time-stepping scheme. For this reason, we use a 2nd order backward differentiation scheme (``bdf2``) with a variable time-step. The ``adaptive time step scaling`` is set to :math:`1.025` to ensure that the time-step does not rise too quickly during wave oscillations. .. code-block:: text @@ -154,6 +158,7 @@ The ``source term`` subsection is used to enable the gravitational acceleration end end +**** ----------------------- Running the Simulation @@ -174,12 +179,13 @@ to run the simulation using eight CPU cores. Feel free to use more. run in parallel using mpirun. This simulation takes :math:`\sim \, 8` minutes (:math:`Re=2`) to :math:`6` hours (:math:`Re=2000`) on :math:`8` processes. +**** ------- Results ------- -We compare the relative height of the free surface at :math:`x=0` with an analytical solution proposed by Wu *et al.* [#wu2001]_ For the Reynolds number of :math:`2`, :math:`20` and :math:`200`, data were directly extracted from Carrica *et al.* [#carrica2007]_, whereas for the Reynolds of :math:`2000`, the simplified analytical expression of Wu *et al.* [#wu2001]_ is used. The results for Reynolds number of :math:`2`, :math:`20`, :math:`200` and :math:`2000` can be post-processed by invoking the following command from the folder of the Reynolds number of interest (:math:`Re=20` in the example below): +We compare the relative height of the free surface at :math:`x=0` with an analytical solution proposed by Wu *et al.* [#wu2001]_ For the Reynolds number of :math:`2`, :math:`20` and :math:`200`, data were directly extracted from Carrica *et al.* [#carrica2007]_, whereas for the Reynolds of :math:`2000`, the simplified analytical expression of Wu *et al.* [#wu2001]_ is used. The results for Reynolds number of :math:`2`, :math:`20`, :math:`200` and :math:`2000` can be postprocessed by invoking the following command from the folder of the Reynolds number of interest (:math:`Re=20` in the example below): .. code-block:: text :class: copy-button @@ -212,6 +218,7 @@ The following table presents a comparison between the analytical results and the | | :align: center | +-------------+--------------------------------------------------------------------------------------------------------+ +**** ---------- References diff --git a/examples/multiphysics/3d-dam-break/3d-dam-break_postprocess.py b/examples/multiphysics/3d-dam-break/3d-dam-break_postprocess.py index 0528c0326e..fe3daf39e8 100644 --- a/examples/multiphysics/3d-dam-break/3d-dam-break_postprocess.py +++ b/examples/multiphysics/3d-dam-break/3d-dam-break_postprocess.py @@ -104,5 +104,5 @@ exec(f'ax{i-1}.legend(loc="best")') exec(f'for label in (ax{i-1}.get_xticklabels() + ax{i-1}.get_yticklabels()):\n' f'\t label.set_fontsize(fs2)') -fig0.savefig('./H1_to_H4_evolution_2.png') +fig0.savefig('./H1_to_H4_evolution.png') plt.show() diff --git a/examples/multiphysics/capillary-wave/capillary-wave-TSM-0_95.prm b/examples/multiphysics/capillary-wave/capillary-wave-TSM-0_95.prm index d3e792162f..317f53366e 100644 --- a/examples/multiphysics/capillary-wave/capillary-wave-TSM-0_95.prm +++ b/examples/multiphysics/capillary-wave/capillary-wave-TSM-0_95.prm @@ -62,16 +62,6 @@ subsection initial conditions end end -#--------------------------------------------------- -# Source term -#--------------------------------------------------- - -subsection source term - subsection fluid dynamics - set Function expression = 0; 0; 0 - end -end - #--------------------------------------------------- # Physical Properties #--------------------------------------------------- @@ -125,7 +115,6 @@ subsection mesh adaptation set initial refinement steps = 4 end - # -------------------------------------------------- # Boundary Conditions #--------------------------------------------------- diff --git a/examples/multiphysics/capillary-wave/capillary-wave.tpl b/examples/multiphysics/capillary-wave/capillary-wave.tpl index f7c56505fe..2388e1419a 100644 --- a/examples/multiphysics/capillary-wave/capillary-wave.tpl +++ b/examples/multiphysics/capillary-wave/capillary-wave.tpl @@ -1,3 +1,6 @@ +# SPDX-FileCopyrightText: Copyright (c) 2023-2024 The Lethe Authors +# SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception OR LGPL-2.1-or-later + # Listing of Parameters # --------------------- @@ -59,16 +62,6 @@ subsection initial conditions end end -#--------------------------------------------------- -# Source term -#--------------------------------------------------- - -subsection source term - subsection fluid dynamics - set Function expression = 0; 0; 0 - end -end - #--------------------------------------------------- # Physical Properties #--------------------------------------------------- @@ -122,7 +115,6 @@ subsection mesh adaptation set initial refinement steps = 4 end - # -------------------------------------------------- # Boundary Conditions #--------------------------------------------------- @@ -145,6 +137,10 @@ subsection boundary conditions end end +subsection boundary conditions VOF + set number = 4 +end + #--------------------------------------------------- # FEM #--------------------------------------------------- diff --git a/examples/multiphysics/rayleigh-taylor-instability/rayleigh-taylor-instability-adaptive-sharpening.prm b/examples/multiphysics/rayleigh-taylor-instability/rayleigh-taylor-instability-adaptive-sharpening.prm index 992273f6b2..7d5a2baa06 100644 --- a/examples/multiphysics/rayleigh-taylor-instability/rayleigh-taylor-instability-adaptive-sharpening.prm +++ b/examples/multiphysics/rayleigh-taylor-instability/rayleigh-taylor-instability-adaptive-sharpening.prm @@ -137,18 +137,20 @@ end #--------------------------------------------------- subsection boundary conditions - set number = 3 + set number = 4 subsection bc 0 - set id = 0 - set type = periodic - set periodic_id = 1 - set periodic_direction = 0 + set id = 0 + set type = slip end subsection bc 1 + set id = 1 + set type = slip + end + subsection bc 2 set id = 2 set type = noslip end - subsection bc 2 + subsection bc 3 set id = 3 set type = noslip end diff --git a/examples/multiphysics/rayleigh-taylor-instability/rayleigh-taylor-instability-constant-sharpening.prm b/examples/multiphysics/rayleigh-taylor-instability/rayleigh-taylor-instability-constant-sharpening.prm index 4fa1f216f0..fd42d3401f 100644 --- a/examples/multiphysics/rayleigh-taylor-instability/rayleigh-taylor-instability-constant-sharpening.prm +++ b/examples/multiphysics/rayleigh-taylor-instability/rayleigh-taylor-instability-constant-sharpening.prm @@ -133,18 +133,20 @@ end #--------------------------------------------------- subsection boundary conditions - set number = 3 + set number = 4 subsection bc 0 - set id = 0 - set type = periodic - set periodic_id = 1 - set periodic_direction = 0 + set id = 0 + set type = slip end subsection bc 1 + set id = 1 + set type = slip + end + subsection bc 2 set id = 2 set type = noslip end - subsection bc 2 + subsection bc 3 set id = 3 set type = noslip end diff --git a/examples/multiphysics/rayleigh-taylor-instability/rayleigh-taylor_postprocess.py b/examples/multiphysics/rayleigh-taylor-instability/rayleigh-taylor_postprocess.py index 38e477ab0e..cd1ffeba91 100644 --- a/examples/multiphysics/rayleigh-taylor-instability/rayleigh-taylor_postprocess.py +++ b/examples/multiphysics/rayleigh-taylor-instability/rayleigh-taylor_postprocess.py @@ -116,7 +116,7 @@ ax0.set_xlim([0, 4.5]) ax0.set_ylim([0.1, 0.8]) ax0.legend(loc="upper left") -plt.title("Spike and bubble evolution with {} interface sharpening".format(output_path[9:-1])) +plt.title("Spike and bubble evolution with {} interface sharpening\n".format(output_path[9:-1])) fig0.savefig('./spike_and_bubble_evolution_{}.png'.format(output_path[9:-1])) plt.show() @@ -159,11 +159,14 @@ def read_my_data(results_path): list_of_list_of_vars_name,list_of_list_of_vars=read_my_data(output_path + "/mass_conservation_information.dat") +relative_mass_variation_fluid_1 = (list_of_list_of_vars[0][6]/list_of_list_of_vars[0][6][0]-1)*100 -plt.plot(list_of_list_of_vars[0][0],list_of_list_of_vars[0][4]) -plt.title("Evolution of the mass of fluid 1 with {} interface sharpening".format(output_path[9:-1])) -plt.xlabel("Time (s)") -plt.ylabel("Mass of fluid 1 (kg)") -plt.ylim(37.4,37.7) -plt.savefig('./mass_of_fluid_1_{}.png'.format(output_path[9:-1])) +dpi_val = 300 +fig0 = plt.figure(figsize=(1920/dpi_val, 1440/dpi_val)) +plt.plot(list_of_list_of_vars[0][0],relative_mass_variation_fluid_1) +plt.title("Evolution of the mass of fluid 1 with {} interface sharpening\n".format(output_path[9:-1])) +plt.xlabel(r"Time $(t)$ [s]") +plt.ylabel(r"Relative mass variation percentage $\left(\frac{m_t-m_0}{m_0} \ 100\% \right)$") +plt.ylim(-0.2,0.4) +plt.savefig('./mass_of_fluid_1_{}.png'.format(output_path[9:-1]), dpi=dpi_val) plt.show() diff --git a/source/solvers/vof_assemblers.cc b/source/solvers/vof_assemblers.cc index 38b5fe755b..af5db222f3 100644 --- a/source/solvers/vof_assemblers.cc +++ b/source/solvers/vof_assemblers.cc @@ -364,8 +364,7 @@ VOFAssemblerDCDDStabilization::assemble_matrix( // Directionality tensor describing [rr - (r⋅s)^2*ss] const Tensor<2, dim> dir_tensor = outer_product(gradient_unit_vector, gradient_unit_vector) - - (Utilities::fixed_power<2>(gradient_unit_vector * - velocity_unit_vector) * + ((gradient_unit_vector * velocity_unit_vector) * outer_product(velocity_unit_vector, velocity_unit_vector)); // Compute the artificial viscosity of the shock capture @@ -436,8 +435,7 @@ VOFAssemblerDCDDStabilization::assemble_rhs( // Directionality tensor describing [rr - (r⋅s)^2*ss] const Tensor<2, dim> dir_tensor = outer_product(gradient_unit_vector, gradient_unit_vector) - - (Utilities::fixed_power<2>(gradient_unit_vector * - velocity_unit_vector) * + ((gradient_unit_vector * velocity_unit_vector) * outer_product(velocity_unit_vector, velocity_unit_vector)); // Compute the artificial viscosity of the shock capture