3d_static_cap.cc File Reference

#include "generic.h"
#include "navier_stokes.h"
#include "fluid_interface.h"
#include "constitutive.h"
#include "solid.h"
#include "meshes/quarter_tube_mesh.h"

Classes
class	AxialSpineQuarterTubeMesh< ELEMENT, INTERFACE_ELEMENT >
class	FreeSurfaceRotationProblem< ELEMENT >
	Entry flow problem in quarter tube domain. More...

Namespaces
	Global_Physical_Variables
	Global variables.
	WallFunction

Functions
void	WallFunction::normal (const Vector< double > &x, Vector< double > &normal)
int	main (int argc, char *argv[])

Variables
Vector< double >	Global_Physical_Variables::G (3) = 1.0
	Direction of gravity. More...

Function Documentation

main()

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

Driver for 3D entry flow into a quarter tube. If there are any command line arguments, we regard this as a validation run and perform only a single adaptation

{
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 // Allow (up to) five rounds of fully automatic adapation in response to 
 //-----------------------------------------------------------------------
 // error estimate
 //---------------
 //unsigned max_adapt;
 double max_error_target,min_error_target;
 
 using namespace Global_Physical_Variables;
 
 G[0] = 0.0; G[1] = 0.0; G[2] = -1.0;
 
 
 // Set max number of adaptations in black-box Newton solver and
 // error targets for adaptation
 if (CommandLineArgs::Argc==1)
  {
   // Up to five adaptations
   //max_adapt=5;
 
   // Error targets for adaptive refinement
   max_error_target=0.005;
   min_error_target=0.0005;
  } 
 // Validation run: Only one adaptation. Relax error targets
 // to ensure that not all elements are refined so we're getting
 // some hanging nodes.
 else
  {
   // Validation run: Just one round of adaptation
   //max_adapt=1;
   
   // Error targets for adaptive refinement
   max_error_target=0.02;
   min_error_target=0.002;
  }
 // end max_adapt setup
 
 
 unsigned n_angles = 1;
 
 // Set up doc info
 DocInfo doc_info;
 
 // Do Taylor-Hood elements
 //------------------------
 {
  // Set output directory
  doc_info.set_directory("RESLT_TH_internal");
  
  // Step number
  doc_info.number()=0;
  
  // Build problem
  FreeSurfaceRotationProblem<
   Hijacked<SpineElement<RefineableQTaylorHoodElement<3> > > >
   problem(doc_info,min_error_target,max_error_target,1);
  
  cout << " Doing Taylor-Hood elements (internal hijack) " << std::endl;
 
  // Doc solution after solving
  problem.doc_solution();
  
  // Solve the problem  (DO NOT ADAPT)
  problem.newton_solve();
 
  //Decrease the contact angle
  for(unsigned i=0;i<n_angles;i++)
   {
    Angle -= 0.1;
    problem.newton_solve();
   }
   }
 
 
 // Do Taylor-Hood elements
 //------------------------
 {
  // Set output directory
  doc_info.set_directory("RESLT_TH_external");
  
  // Step number
  doc_info.number()=0;
  
  // Build problem
  FreeSurfaceRotationProblem<
   Hijacked<SpineElement<RefineableQTaylorHoodElement<3> > > >
   problem(doc_info,min_error_target,max_error_target,0);
  
  cout << " Doing Taylor-Hood elements (external hijack)" << std::endl;
 
  // Doc solution after solving
  problem.doc_solution();
  
  // Solve the problem  (DO NOT ADAPT)
  problem.newton_solve();
  
  for(unsigned i=0;i<n_angles;i++)
   {
    Angle -= 0.1;
    problem.newton_solve();
   }
   }
 
 
 // Do Crouzeix-Raviart elements
 //-----------------------------
 {
  // Set output directory
  doc_info.set_directory("RESLT_CR_internal");
  
  // Step number
  doc_info.number()=0;
  
  // Build problem
  FreeSurfaceRotationProblem<
   Hijacked<SpineElement<RefineableQCrouzeixRaviartElement<3> > > >
   problem(doc_info,min_error_target,max_error_target,2);
  
  cout << " Doing Crouzeix-Raviart elements (internal hijack) " << std::endl;
  
  // Solve the problem  DO NOT ADAPT
  problem.newton_solve();
 
   for(unsigned i=0;i<n_angles;i++)
   {
    Angle -= 0.1;
    problem.newton_solve();
   }
 
  }
 
 // Do Crouzeix-Raviart elements
 //-----------------------------
 {
  // Set output directory
  doc_info.set_directory("RESLT_CR_external");
  
  // Step number
  doc_info.number()=0;
  
  // Build problem
  FreeSurfaceRotationProblem<
   Hijacked< SpineElement<RefineableQCrouzeixRaviartElement<3> > > >
   problem(doc_info,min_error_target,max_error_target,0);
  
  cout << " Doing Crouzeix-Raviart elements (external hijack) " << std::endl;
  
  // Solve the problem  DO NOT ADAPT
  problem.newton_solve();
 
   for(unsigned i=0;i<n_angles;i++)
   {
    Angle -= 0.1;
    problem.newton_solve();
   }
 
   } 
 
 
} // end_of_main

References Global_Physical_Variables::Angle, oomph::CommandLineArgs::Argc, G, i, oomph::DocInfo::number(), problem, oomph::DocInfo::set_directory(), and Flag_definition::setup().