multimesh_ref_dd_convection.cc File Reference

#include "db_nst_external_elements.h"
#include "meshes/rectangular_quadmesh.h"

Classes
class	RefineableDDConvectionProblem< NST_ELEMENT, AD_ELEMENT >

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
int	main (int argc, char **argv)

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

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

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.
 RefineableDDConvectionProblem<
  RefineableQCrouzeixRaviartElementWithTwoExternalElement<2>,
  RefineableQAdvectionDiffusionElementWithExternalElement<2> > problem;
 
 
 
 // Apply the boundary condition at time zero
 problem.set_boundary_conditions(0.0);
 
 //Perform a single steady Newton solve
 problem.steady_newton_solve();
 
 //Document the solution
 problem.doc_solution();
 
  //Start a trace file
 ofstream trace("RESLT_ref_multimesh/trace.dat");
 //Local variables for the kinetic energy and its rate of change
 double E=0.0, Edot = 0.0;
 
 //Output to the trace file
 problem.get_kinetic_energy(E,Edot);
 trace << problem.time_pt()->time() << " " 
       << problem.nst_mesh_pt()->boundary_node_pt(1,8)->value(1) << " " 
       << E << " " << Edot << std::endl;
 
 //Set the timestep
 double dt = 0.01;
 
 //Initialise the value of the timestep and set an impulsive start
 problem.assign_initial_values_impulsive(dt);
 
 //Set the number of timesteps to our default value
 unsigned n_steps = 2000;
 
 //If we have a command line argument, perform fewer steps 
 //(used for self-test runs)
 if(argc > 1) {n_steps = 5;}
 
 bool first=true;
 
 unsigned n_refine = 0;
 
 //Perform n_steps timesteps
 for(unsigned i=0;i<n_steps;++i)
  {
   problem.unsteady_newton_solve(dt,n_refine,first);
   problem.doc_solution();
 
   first = false;
 
   //Refine the last case for the demo
   if((argc > 1) && (i==3)) {n_refine=1;}
   //Only start to refine after 700 steps once the
   //limit cycle has essentially settled down
   if(i > 700) {n_refine=1;}
 
 //Output to the trace file
   problem.get_kinetic_energy(E,Edot);
   trace << problem.time_pt()->time() << " " 
         << problem.nst_mesh_pt()->boundary_node_pt(1,8)->value(1) << " "
         << E << " " << Edot << std::endl;
  }
 
 trace.close();
 
 
} // end of main

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