full_tube.cc File Reference

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

Classes
class	MyCylinder
class	SteadyTubeProblem< ELEMENT >
	Entry flow problem in tapered tube domain. More...

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
int	main (int argc, char *argv[])
	Driver for 3D entry flow into a straight tube. More...

Function Documentation

main()

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

Driver for 3D entry flow into a straight tube.

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

{
 //The purpose of this code is to test that the position of nodes on 
 //curvilinear boundaries are correctly updated after nodes that were 
 //hanging become non-hanging.
 unsigned max_adapt=1;
 double max_error_target=0.02 , min_error_target=0.002;
 
 
 // Set up doc info
 DocInfo doc_info;
 
 // Do Taylor-Hood elements
 //------------------------
 {
  // Set output directory
  doc_info.set_directory("RESLT_TH");
  
  // Step number
  doc_info.number()=0;
  
  // Build problem
  SteadyTubeProblem<RefineableQTaylorHoodElement<3> > 
   problem(doc_info,min_error_target,max_error_target);
  
  cout << " Doing Taylor-Hood elements " << std::endl;
 
  for(unsigned i=0;i<2;i++)
   {
    // Solve the problem 
    problem.newton_solve(max_adapt);
    // Doc solution after solving
    problem.doc_solution();
 
    Global_Physical_Variables::Re += 25.0;
    ++doc_info.number();
   }
 }
 
 exit(1);
 // Do Crouzeix-Raviart elements
 //------------------------
 {
  // Set output directory
  doc_info.set_directory("RESLT_CR");
  
  // Step number
  doc_info.number()=0;
  
  // Build problem
  SteadyTubeProblem<RefineableQCrouzeixRaviartElement<3> > 
   problem(doc_info,min_error_target,max_error_target);
  
  cout << " Doing Crouzeix-Raviart elements " << std::endl;
 
  for(unsigned i=0;i<2;i++)
   {
    // Solve the problem 
    problem.newton_solve(max_adapt);
    // Doc solution after solving
    problem.doc_solution();
 
    Global_Physical_Variables::Re += 25.0;
    ++doc_info.number();
   }
 }
} // end_of_main

References i, oomph::DocInfo::number(), problem, Global_Physical_Variables::Re, and oomph::DocInfo::set_directory().