axisym_navier_stokes/spin_up/spin_up.cc File Reference

#include "generic.h"
#include "navier_stokes.h"
#include "axisym_navier_stokes.h"
#include "meshes/rectangular_quadmesh.h"

Classes
class	RotatingCylinderProblem< ELEMENT, TIMESTEPPER >

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
int	main (int argc, char *argv[])
	Driver code for axisymmetric spin-up problem. More...

Function Documentation

main()

int main	(	int argc	,
		char *	argv[]
	)

Driver code for axisymmetric spin-up problem.

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

Maximum time

Duration of timestep

{
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 double t_max = 1.0;
 
 const double dt = 0.01;
 
 // If we are doing validation run, use smaller number of timesteps
 if(CommandLineArgs::Argc>1) { t_max = 0.02; }
 
 // Number of elements in radial (r) direction
 const unsigned n_r = 2;
 
 // Number of elements in axial (z) direction
 const unsigned n_z = 2;
 
 // Length in radial (r) direction
 const double l_r = 1.0;
 
 // Length in axial (z) direction
 const double l_z = 1.4;
 
 // -----------------------------------------
 // RefineableAxisymmetricQTaylorHoodElements
 // -----------------------------------------
 {
  cout << "Doing RefineableAxisymmetricQTaylorHoodElement" << std::endl;
 
  // Build the problem with RefineableAxisymmetricQTaylorHoodElements
  RotatingCylinderProblem
   <RefineableAxisymmetricQTaylorHoodElement, BDF<2> > 
   problem(n_r,n_z,l_r,l_z);
  
  // Solve the problem and output the solution
  problem.unsteady_run(t_max,dt,"RESLT_TH");
 }
 
 // ----------------------------------------------
 // RefineableAxisymmetricQCrouzeixRaviartElements
 // ----------------------------------------------
 {
  cout << "Doing RefineableAxisymmetricQCrouzeixRaviartElement" << std::endl;
 
  // Build the problem with RefineableAxisymmetricQCrouzeixRaviartElements
  RotatingCylinderProblem
   <RefineableAxisymmetricQCrouzeixRaviartElement, BDF<2> >
   problem(n_r,n_z,l_r,l_z);
 
  // Solve the problem and output the solution
  problem.unsteady_run(t_max,dt,"RESLT_CR");
 }
 
} // End of main

References oomph::CommandLineArgs::Argc, problem, and Flag_definition::setup().