unstructured_two_d_fsi.cc File Reference

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

Classes
class	MySolidTriangleMesh< ELEMENT >
	Triangle-based mesh upgraded to become a solid mesh. More...
class	FluidTriangleMesh< ELEMENT >
	Triangle-based mesh upgraded to become a pseudo-solid mesh. More...
class	UnstructuredFSIProblem< FLUID_ELEMENT, SOLID_ELEMENT >
	Unstructured 3D FSI problem. More...

Namespaces
	Global_Parameters
	Namespace for global parameters.

Functions
void	Global_Parameters::gravity (const double &time, const Vector< double > &xi, Vector< double > &b)
	Non-dimensional gravity as body force. More...
bool	Global_Parameters::is_on_fsi_boundary (Node *nod_pt)
	Boolean to identify if node is on fsi boundary. More...
int	main ()
	Driver for unstructured fsi problem. More...

Function Documentation

main()

int main

(

)

Driver for unstructured fsi problem.

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

{
 // Label for output
 DocInfo doc_info;
 
 // Set output directory
 doc_info.set_directory("RESLT");
 
 //Create the constitutive law
 Global_Parameters::Constitutive_law_pt = new GeneralisedHookean(
  &Global_Parameters::Nu);
 
 // Build the problem with triangular Taylor Hood for fluid and solid
 UnstructuredFSIProblem<
 PseudoSolidNodeUpdateElement<TTaylorHoodElement<2>, TPVDElement<2,3> >, 
  TPVDElement<2,3> > problem;
 
 // Output boundaries 
 problem.fluid_mesh_pt()->output_boundaries("RESLT/fluid_boundaries.dat");
 problem.solid_mesh_pt()->output_boundaries("RESLT/solid_boundaries.dat");
 
 // Output the initial guess for the solution
 problem.doc_solution(doc_info);
 doc_info.number()++;
 
 // Parameter study
 Global_Parameters::Gravity=2.0e-4;
 double q_increment=1.0e-6;
 
 // Solve the problem at zero Re and Q
 problem.newton_solve();
 
 // Output the solution
 problem.doc_solution(doc_info);
 doc_info.number()++;
 
 // Bump up Re
 Global_Parameters::Re=10.0;
 
 // Now do proper parameter study with increase in Q
 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 Q
   Global_Parameters::Q+=q_increment;
  }
 
 
} // end_of_main

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