circular_driven_cavity_driver.cc File Reference

#include "generic.h"
#include "navier_stokes.h"
#include "meshes/quarter_circle_sector_mesh.h"

Classes
class	DrivenCavityProblem< ELEMENT, MESH >

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
int	main ()
	Driver for DrivenCavity test problem with mesh adaptation. More...

Function Documentation

main()

int main

(

)

Driver for DrivenCavity test problem with mesh adaptation.

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

{
 
 
 typedef AlgebraicElement<RefineableQCrouzeixRaviartElement<2> > ELEMENT0;
 typedef AlgebraicElement<RefineableQTaylorHoodElement<2> > ELEMENT1;
 typedef RefineableQCrouzeixRaviartElement<2> ELEMENT2;
 typedef RefineableQTaylorHoodElement<2> ELEMENT3;
 
 typedef AlgebraicRefineableQuarterCircleSectorMesh<ELEMENT0> MESH0;
 typedef AlgebraicRefineableQuarterCircleSectorMesh<ELEMENT1> MESH1;
 typedef RefineableQuarterCircleSectorMesh<ELEMENT2> MESH2;
 typedef RefineableQuarterCircleSectorMesh<ELEMENT3> MESH3;
 
// Counter for steps
 unsigned next_istep=0; 
  
 {
  
  cout << " " << std::endl;
  cout << "==========================================================" << std::endl;
  cout << "Algebraic Taylor Hood " << std::endl;
  cout << "==========================================================" << std::endl;
  cout << " " << std::endl;
  DrivenCavityProblem<ELEMENT1,MESH1>* problem_pt= new
   DrivenCavityProblem<ELEMENT1,MESH1>;
  
  problem_pt->run(next_istep);
 }
 
 
 {
  cout << " " << std::endl;
  cout << "==========================================================" << std::endl;
  cout << "Algebraic Crouzeix Raviart " << std::endl;
  cout << "==========================================================" << std::endl;
  cout << " " << std::endl;
  DrivenCavityProblem<ELEMENT0,MESH0>* problem_pt= new
   DrivenCavityProblem<ELEMENT0,MESH0>;
  
  problem_pt->run(next_istep);
 }
 
 {
  cout << " " << std::endl;
  cout << "==========================================================" << std::endl;
  cout << "Taylor Hood " << std::endl;
  cout << "==========================================================" << std::endl;
  cout << " " << std::endl;
  DrivenCavityProblem<ELEMENT3,MESH3>* problem_pt= new
   DrivenCavityProblem<ELEMENT3,MESH3>;
  
  problem_pt->run(next_istep);
 }
 
 {
  cout << " " << std::endl;
  cout << "==========================================================" << std::endl;
  cout << "Crouzeix Raviart " << std::endl;
  cout << "==========================================================" << std::endl;
  cout << " " << std::endl;
  DrivenCavityProblem<ELEMENT2,MESH2>* problem_pt= new
   DrivenCavityProblem<ELEMENT2,MESH2>;
  
  problem_pt->run(next_istep);
 }
 
 
}

References DrivenCavityProblem< ELEMENT, MESH >::run().