multi_physics/boussinesq_convection/multi_domain_ref_b_convection.cc File Reference

#include "generic.h"
#include "advection_diffusion.h"
#include "navier_stokes.h"
#include "multi_physics.h"
#include "meshes/rectangular_quadmesh.h"

Classes
class	RefineableConvectionProblem< NST_ELEMENT, AD_ELEMENT >

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
int	main ()

Variables
Vector< double >	Global_Physical_Variables::Direction_of_gravity (2)
	Gravity vector. More...

Function Documentation

main()

int main

(

)

Driver code for 2D Boussinesq convection problem with adaptivity.

{
 
 
 // Set the direction of gravity
 Global_Physical_Variables::Direction_of_gravity[0] = 0.0;
 Global_Physical_Variables::Direction_of_gravity[1] = -1.0;
 
 // Create the problem with 2D nine-node refineable elements.
 RefineableConvectionProblem<
 RefineableNavierStokesBoussinesqElement<
 RefineableQCrouzeixRaviartElement<2>,
  RefineableQAdvectionDiffusionElement<2,3> > ,
  RefineableAdvectionDiffusionBoussinesqElement<
  RefineableQAdvectionDiffusionElement<2,3>,
  RefineableQCrouzeixRaviartElement<2> >
  > problem;
 
 // Apply a perturbation to the upper boundary condition to
 // force the solution into the symmetry-broken state.
 problem.enable_imperfection();
 
 //Solve the problem with (up to) two levels of adaptation
 problem.newton_solve(2);
 
 //Document the solution
 problem.doc_solution();
 
 // Make the boundary conditions perfect and solve again. 
 // Now the slightly perturbed symmetry broken state computed
 // above is used as the initial condition and the Newton solver
 // converges to the symmetry broken solution, even without
 // the perturbation
 problem.disable_imperfection();
 problem.newton_solve(2);
 problem.doc_solution();
 
} // end of main

References oomph::Global_Physical_Variables::Direction_of_gravity, and problem.