unstructured_two_d_fluid.cc File Reference

#include "generic.h"
#include "navier_stokes.h"
#include "solid.h"
#include "constitutive.h"
#include "meshes/triangle_mesh.h"

Classes
class	ElasticTriangleMesh< ELEMENT >
	Triangle-based mesh upgraded to become a (pseudo-) solid mesh. More...
class	UnstructuredFluidProblem< ELEMENT >
	Unstructured Navier-Stokes ALE Problem. More...

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
int	main ()
	Driver for unstructured fluid test problem. More...

Function Documentation

main()

int main

(

)

Driver for unstructured fluid test problem.

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

{
 // Label for output
 DocInfo doc_info;
 
 //Set the constitutive law for the pseudo-elasticity
 Global_Physical_Variables::Constitutive_law_pt = 
  new GeneralisedHookean(&Global_Physical_Variables::Nu);
 
 //Taylor--Hood formulation
 {
  // Set output directory
  doc_info.set_directory("RESLT_TH");
  
  // Step number
  doc_info.number()=0;
  
  // Build the problem with TTaylorHoodElements
  UnstructuredFluidProblem<PseudoSolidNodeUpdateElement<
   TTaylorHoodElement<2>, 
   TPVDElement<2,3> > > problem;
  
  // Output boundaries 
  problem.fluid_mesh_pt()->output_boundaries("RESLT_TH/boundaries.dat");
  
  // Output the initial guess for the solution
  problem.doc_solution(doc_info);
  doc_info.number()++;
  
  // Parameter study
  double re_increment=5.0;
  unsigned nstep=2; // 10;
  for (unsigned i=0;i<nstep;i++)
   {
    // Solve the problem
    problem.newton_solve();
    
    // Output the solution
    problem.doc_solution(doc_info);
    doc_info.number()++;
    
    // Bump up Re
    Global_Physical_Variables::Re+=re_increment;
   }  //end_of_parameter_study
 }
 
 
 
 //Crouzeix--Raviart formulation
 {
 // Set output directory
 doc_info.set_directory("RESLT_CR");
 
 // Step number
 doc_info.number()=0;
 
 // Reset Reynolds number
 Global_Physical_Variables::Re = 0.0;
 
 // Build the problem with TCrouzeixRaviartElements
 UnstructuredFluidProblem<PseudoSolidNodeUpdateElement<
  TCrouzeixRaviartElement<2>, 
  TPVDBubbleEnrichedElement<2,3> > > problem;
 
 // Output boundaries 
  problem.fluid_mesh_pt()->output_boundaries("RESLT_CR/boundaries.dat");
  
  // Output the initial guess for the solution
  problem.doc_solution(doc_info);
  doc_info.number()++;
  
  // Parameter study
  double re_increment=5.0;
  unsigned nstep=2; // 10;
  for (unsigned i=0;i<nstep;i++)
   {
    // Solve the problem
    problem.newton_solve();
    
    // Output the solution
    problem.doc_solution(doc_info);
    doc_info.number()++;
    
    // Bump up Re
    Global_Physical_Variables::Re+=re_increment;
   }
 }
 
 
} // end_of_main

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