ray_circ_cavity_adapt.cc File Reference

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

Classes
class	MyEllipse
class	RayleighProblem< ELEMENT, TIMESTEPPER >

Namespaces
	Global_Parameters
	Namespace for global parameters.

Functions
int	main (int argc, char *argv[])
	Driver code for Rayleigh channel problem. More...

Function Documentation

main()

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

Driver code for Rayleigh channel problem.

Convert command line arguments (if any) into flags:

{
 
 if (argc==1)
  {
   cout << "No command line arguments specified -- using defaults." 
        << std::endl;
  }
 else if (argc>=2)
  {
   cout << "Command line arguments specified:" << std::endl;
   // Flag for long run
   Global_Parameters::Long_run_flag=0;
  }
 
 cout << "Long run flag: " 
      <<  Global_Parameters::Long_run_flag << std::endl;
 
 
 // Set physical parameters:
 
 // Womersley number = Reynolds number (St = 1)
 Global_Parameters::ReSt = 10.0;
 Global_Parameters::Re = Global_Parameters::ReSt;
 
 //Horizontal length of domain
 double lx = 1.0;
 
 //Vertical length of domain
 double ly = 1.0;
 
 // Number of elements in x-direction
 unsigned nx=5;
 
 // Number of elements in y-direction
 unsigned ny=10;
 
 // Loop over run with/without ALE
 for (unsigned ale_flag=0;ale_flag<2;ale_flag++)
  {
 
   bool use_ale=true;
   if (ale_flag==1) use_ale=false;
 
 
   // Solve with Crouzeix-Raviart elements
   {
    // Set up doc info
    DocInfo doc_info;
    doc_info.number()=0;
    if (ale_flag==0)
     {
      doc_info.set_directory("RESLT2_CR");
     }
    else
     {
      doc_info.set_directory("RESLT2_CR_ALE");
     }
 
    //Set up problem
    RayleighProblem<RefineableQCrouzeixRaviartElement<2>,BDF<2> > 
     problem(nx,ny,lx,ly,use_ale);
    
    // Run the unsteady simulation
    problem.unsteady_run(doc_info);
   }
   
   // Solve with Taylor-Hood elements
   {
    // Set up doc info
    DocInfo doc_info;
    doc_info.number()=0;
    if (ale_flag==0)
     {
      doc_info.set_directory("RESLT2_TH");
     }
    else
     {
      doc_info.set_directory("RESLT2_TH_ALE");
     }
    
    //Set up problem
    RayleighProblem<RefineableQTaylorHoodElement<2>,BDF<2> >
     problem(nx,ny,lx,ly,use_ale);
    
    // Run the unsteady simulation
    problem.unsteady_run(doc_info);
   }
  }
 
} // end of main

References Global_Parameters::Long_run_flag, Mesh_Parameters::lx, Mesh_Parameters::ly, oomph::DocInfo::number(), Mesh_Parameters::nx, Mesh_Parameters::ny, problem, Global_Parameters::Re, Global_Parameters::ReSt, and oomph::DocInfo::set_directory().