]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Let gravity point outward instead of upward. Use instable temperature layering.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Thu, 18 Oct 2007 20:26:19 +0000 (20:26 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Thu, 18 Oct 2007 20:26:19 +0000 (20:26 +0000)
git-svn-id: https://svn.dealii.org/trunk@15346 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-22/step-22.cc

index 861d753592b45b7b33bc5bfcfdefb2088ea1fc7a..b65ec3f140ad18886a95e4606240a31df97dd92a 100644 (file)
@@ -164,7 +164,7 @@ InitialValues<dim>::value (const Point<dim>  &p,
                            const unsigned int component) const 
 {
   if (component == dim+1)
-    return (p.distance (Point<dim>(.05,-.95)) < .05 ? 1 : 0);
+    return (p.distance (Point<dim>()) < .6 ? 1 : 0);
   else
     return 0;
 }
@@ -560,9 +560,12 @@ void BoussinesqFlowProblem<dim>::assemble_system ()
                                       * fe_values.JxW(q);     
                }
 
+             const Point<dim> 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);
           }
        }

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.