From: bangerth Date: Thu, 18 Oct 2007 20:26:19 +0000 (+0000) Subject: Let gravity point outward instead of upward. Use instable temperature layering. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=f7c1742b5b444990359a1e66d5118a9b32f0203b;p=dealii-svn.git 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 --- 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); } }