From af59c0e18f0edeedfd1f8d3da4308f32a792e657 Mon Sep 17 00:00:00 2001 From: bangerth Date: Thu, 29 Jan 2009 15:26:17 +0000 Subject: [PATCH] Bring step-32 back in line with step-31 (at least for the moment, so that we can compare more easily): revert to artificial viscosity with beta=0.015, alpha=1; use gravity pointing upward; introduce the factor of 2 in the Stokes matrix that makes the difference between the symmetric gradient and the gradient. git-svn-id: https://svn.dealii.org/trunk@18310 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-32/step-32.cc | 13 +++++-------- 1 file changed, 5 insertions(+), 8 deletions(-) diff --git a/deal.II/examples/step-32/step-32.cc b/deal.II/examples/step-32/step-32.cc index 7c2e42c43a..2598c73158 100644 --- a/deal.II/examples/step-32/step-32.cc +++ b/deal.II/examples/step-32/step-32.cc @@ -558,8 +558,8 @@ compute_viscosity (const std::vector &old_temperature, const double global_T_variation, const double cell_diameter) { - const double beta = 0.04 * dim; - const double alpha = 2; + const double beta = 0.015 * dim; + const double alpha = 1; if (global_u_infty == 0) return 5e-3 * cell_diameter; @@ -950,15 +950,14 @@ void BoussinesqFlowProblem::assemble_stokes_system () if (rebuild_stokes_matrix) for (unsigned int i=0; i gravity = stokes_fe_values.quadrature_point(q) / - stokes_fe_values.quadrature_point(q).norm(); + const Point gravity = ( (dim == 2) ? (Point (0,1)) : + (Point (0,0,1)) ); for (unsigned int i=0; i::assemble_temperature_system () const double global_u_infty = get_maximal_velocity(); const std::pair global_T_range = get_extrapolated_temperature_range(); - const double global_Omega_diameter = GridTools::diameter (triangulation); const TrilinosWrappers::BlockVector localized_stokes_solution (stokes_solution); @@ -1648,7 +1646,6 @@ void BoussinesqFlowProblem::run () old_temperature_solution = temperature_solution; } while (time <= 100); - } -- 2.39.5