3d_multimesh_ref_b_convect.cc File Reference

#include "my_boussinesq_elements.h"
#include "multi_physics.h"
#include "meshes/simple_cubic_mesh.h"

Classes
class	ElementCmp
	Function-type-object to perform comparison of elements. More...
class	RefineableConvectionProblem< NST_ELEMENT, AD_ELEMENT >

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
Vector< double >	Global_Physical_Variables::Direction_of_gravity (3)
	Gravity vector. More...
int	main (int argc, char **argv)

Variables
double	Global_Physical_Variables::Reynolds =10.0
	Reynolds number. More...

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

Driver code for 3D Boussinesq convection problem with adaptivity.

{
 
 
//Uncomment this to turn this into a parallel code
//#ifdef OOMPH_HAS_MPI
// // Set up MPI_Helpers
// MPI_Helpers::init(argc,argv);
// #endif
 
 // Set the direction of gravity
 Global_Physical_Variables::Direction_of_gravity[0] = 0.0;
 Global_Physical_Variables::Direction_of_gravity[1] = -1.0;
 Global_Physical_Variables::Direction_of_gravity[2] = 0.0;
 
 // Create the problem with 3D twenty-seven-node refineable elements.
 RefineableConvectionProblem<
  RefineableQCrouzeixRaviartElementWithExternalElement<3>,
  RefineableQAdvectionDiffusionElementWithExternalElement<3> > problem;
 
 // Distribute the problem (including set sources)
// #ifdef OOMPH_HAS_MPI
//  DocInfo mesh_doc_info;
//  bool report_stats=true;
//  problem.distribute(mesh_doc_info,report_stats);
// #endif
  
 //Solve the problem with (up to) two levels of adaptation
 //problem.switch_to_iterative_linear_solver();
 problem.newton_solve();
 //Document the solution
 problem.doc_solution();
 
 
//#ifdef OOMPH_HAS_MPI 
// // finalize MPI
// MPI_Helpers::finalize();
//#endif
 
} // end of main

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