vmtk_fsi.cc File Reference

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

Classes
class	MySolidTetgenMesh< ELEMENT >
	Tetgen-based mesh upgraded to become a solid mesh. More...
class	MyFluidTetMesh< ELEMENT >
	Tetgen-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.

Enumerations
enum	{ Global_Parameters::Parallel_flow_lagrange_multiplier_id , Global_Parameters::FSI_interface_displacement_lagrange_multiplier_id }

Functions
int	main (int argc, char **argv)
	Demonstrate how to solve an unstructured 3D FSI problem. More...

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

Demonstrate how to solve an unstructured 3D FSI problem.

{
 
 // Label for output
 DocInfo doc_info;
 
 // Output directory
 doc_info.set_directory("RESLT");
 
 // Create generalised Hookean constitutive equations
 Global_Parameters::Constitutive_law_pt = 
  new GeneralisedHookean(&Global_Parameters::Nu);
 
 //Set up the problem
 UnstructuredFSIProblem<
 PseudoSolidNodeUpdateElement<TTaylorHoodElement<3>, TPVDElement<3,3> >,
  TPVDElement<3,3> > problem;
 
 //Output initial configuration
 problem.doc_solution(doc_info);
 doc_info.number()++;   
 
 // Parameter study
 unsigned nstep=2;
 
 // Increment in FSI parameter
 double q_increment=5.0e-2;
 
 for (unsigned istep=0;istep<nstep;istep++)
  {
   // Solve the problem
   problem.newton_solve();
   
   //Output solution
   problem.doc_solution(doc_info);
   doc_info.number()++;
 
   // Bump up FSI parameter
   Global_Parameters::Q+=q_increment;   
  }
 
} // end_of_main

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