flow_past_cylinder.cc File Reference

#include "generic.h"
#include "navier_stokes.h"
#include "assert.h"

Classes
class	GeneralEllipse
	My own Ellipse class. More...
class	RectangleWithHoleDomain
	Rectangular domain with circular whole. More...
class	RectangleWithHoleMesh< ELEMENT >
	Domain-based mesh for rectangular mesh with circular hole. More...
class	RefineableRectangleWithHoleMesh< ELEMENT >
class	FlowAroundCylinderProblem< ELEMENT >
	Flow around a cylinder in rectangular domain. More...

Namespaces
	Global_Parameters
	Namespace for global parameters.

Functions
int	main ()
	Driver. More...

Function Documentation

main()

int main

(

)

Driver.

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

{
 
 // Length and height of domain
 double length=2.2/0.1;
 double height=0.41/0.1;
 
 //Create a new ellipse object as the central cylinder
 GeneralEllipse* cylinder_pt = 
  new GeneralEllipse(0.2/0.1,0.2/0.1,0.05/0.1,0.05/0.1);
 
 // Create Problem
 FlowAroundCylinderProblem 
 <RefineableQCrouzeixRaviartElement<2> > problem(cylinder_pt,length,height);
 
 //Refine the problem a couple of times
 problem.refine_uniformly();
 problem.refine_uniformly();
 
 // Solve adaptively with up to max_adapt rounds of refinement
 unsigned max_adapt=5;
 problem.newton_solve(max_adapt);
 
 //Doc result
 ofstream outfile("soln.dat");
 problem.mesh_pt()->output(outfile,2);
 outfile.close();
 
}

References Global_Physical_Variables::height(), and problem.