spherical_shell_convection.cc File Reference

#include "generic.h"
#include "spherical_advection_diffusion.h"
#include "spherical_buoyant_navier_stokes.h"
#include "meshes/rectangular_quadmesh.h"

Classes
class	RefineableSphereConvectionProblem< ELEMENT >
	Convection in a spherical shell. More...

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
int	main (int argc, char **argv)
	Driver for RefineableSphereConvection test problem using a. More...

Variables
double	Global_Physical_Variables::Ra = 10.0
	Rayleigh number. More...

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

Driver for RefineableSphereConvection test problem using a.

///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// /////////////////////////////////////////////////////////////////////

{
 //Set gravity
 Global_Physical_Variables::G.resize(3,0.0);
 Global_Physical_Variables::G[0] = -1.0;
 
 // Label for output
 DocInfo doc_info;
 
 // Set output directory
 doc_info.set_directory("RESLT_adapt");
 
 // ---------------
 // end of Set up doc info
 
 // Build the problem with QCrouzeixRaviartElements
 RefineableSphereConvectionProblem<
 RefineableBuoyantQSphericalCrouzeixRaviartElement> 
  problem;
 
 // Over-ride the maximum and minimum permitted errors
 problem.mesh_pt()->max_permitted_error() = 1.0e-3; //Default = 1.0e-3
 problem.mesh_pt()->min_permitted_error() = 1.0e-5; //Default = 1.0e-5
 
 // Over-ride the maximum and minimum permitted refinement levels
 problem.mesh_pt()->max_refinement_level() = 4;//maximum_ref_level;
 problem.mesh_pt()->min_refinement_level() = 1;//minimum_ref_level;
 
  //Set the boundary conditions
 problem.set_boundary_conditions(0.0);
 
 //Set the maximum adaptation
 unsigned max_adapt = 2;
  
 // Solve the problem
 problem.steady_newton_solve(max_adapt);
 problem.doc_solution(doc_info);
 
 //Crank up the Rayleigh number and kick it
 Global_Physical_Variables::Ra = 2500.0;
 
 //Set the timestep
 
 //Set the timestep
 double dt = 1.0;
 
 //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 = 100; //This is enough to get to steady state
 
 //If we have a command line argument, perform fewer steps
 if(argc > 1) {n_steps = 5;}
 
 //Perform n_steps timesteps
 for(unsigned i=0;i<n_steps;++i)
  {
   problem.unsteady_newton_solve(dt);
   problem.doc_solution(doc_info);
  }
 
 
 
} // end_of_main

References Global_Physical_Variables::G, i, problem, Global_Physical_Variables::Ra, and oomph::DocInfo::set_directory().