unstructured_torus.cc File Reference

#include "generic.h"
#include "axisym_navier_stokes.h"
#include "meshes/triangle_mesh.h"

Classes
class	GeneralCircle
class	oomph::MyAxisymmetricFluidElement< ELEMENT >
	Overload Element to allow calculation of the flux. More...
class	oomph::FaceGeometry< MyAxisymmetricFluidElement< ELEMENT > >
class	UnstructuredTorusProblem< ELEMENT >
	Solve the Axisymmetric Navier Stokes equations in a torus. More...

Namespaces
	Global_Physical_Variables
	Global variables.
	oomph
	DRAIG: Change all instances of (SPATIAL_DIM) to (DIM-1).

Functions
int	main ()

Function Documentation

main()

int main

(

)

{
 
 //Construct and solve the problem
 //This maximum refinement level means that we fit (easily) into a 
 //1G machine, but it could probably go higher if you start to 
 //see refinement being overruled.
 //If you have this too high you get ridiculous refinement at early 
 //times that isn't really necessary.
 double max_error = 1.0e-3;
 double min_error = 1.0e-5;
 
 
{
  UnstructuredTorusProblem<MyAxisymmetricFluidElement<
   ProjectableAxisymmetricCrouzeixRaviartElement<
   AxisymmetricTCrouzeixRaviartElement> > >  problem(
    min_error,max_error);
  
  //Now timestep
  problem.solve_system(0.01,2,"RESLT_CR");
 }
 
 {
  UnstructuredTorusProblem<MyAxisymmetricFluidElement<
   ProjectableAxisymmetricTaylorHoodElement<
   AxisymmetricTTaylorHoodElement> > >  problem(
    min_error,max_error);
  
  //Now timestep
  problem.solve_system(0.01,2,"RESLT_TH");
 }
 
 
 }

References MeshRefinement::max_error, MeshRefinement::min_error, and problem.