From 4369a8c88c5f1ee335fcedb953f058b714729946 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Thu, 18 Oct 2007 20:26:19 +0000 Subject: [PATCH] Let gravity point outward instead of upward. Use instable temperature layering. git-svn-id: https://svn.dealii.org/trunk@15346 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-22/step-22.cc | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) diff --git a/deal.II/examples/step-22/step-22.cc b/deal.II/examples/step-22/step-22.cc index 861d753592..b65ec3f140 100644 --- a/deal.II/examples/step-22/step-22.cc +++ b/deal.II/examples/step-22/step-22.cc @@ -164,7 +164,7 @@ InitialValues::value (const Point &p, const unsigned int component) const { if (component == dim+1) - return (p.distance (Point(.05,-.95)) < .05 ? 1 : 0); + return (p.distance (Point()) < .6 ? 1 : 0); else return 0; } @@ -560,9 +560,12 @@ void BoussinesqFlowProblem::assemble_system () * fe_values.JxW(q); } + const Point gravity = fe_values.quadrature_point(q) / + fe_values.quadrature_point(q).norm(); + Assert (dim == 2, ExcInternalError()); local_rhs(i) += (Raleigh_number * - phi_i_u[1] * old_temperature)* + gravity * phi_i_u * old_temperature)* fe_values.JxW(q); } } -- 2.39.5