jeffery_orbit.cc File Reference

#include "generic.h"
#include "navier_stokes.h"
#include "solid.h"
#include "constitutive.h"
#include "rigid_body.h"
#include "meshes/triangle_mesh.h"
#include "my_taylor_hood_elements.h"
#include <algorithm>

Classes
class	GeneralEllipse
	My own Ellipse class. More...
class	UnstructuredImmersedEllipseProblem< ELEMENT >

Namespaces
	Problem_Parameter
	Namespace for Problem Parameter.
	Jeffery_Solution

Functions
double	Jeffery_Solution::null (const double &t)
	Null function. More...
double	Jeffery_Solution::angle (const double &t)
	Angular position as a function of time t. More...
double	Jeffery_Solution::velocity (const double &t)
	Angular velocity as function of time t. More...
double	Jeffery_Solution::acceleration (const double &t)
	Angular acceleration as a function of time t (should always be zero) More...
int	main (int argc, char **argv)
	Driver code for immersed ellipse problem. More...

Variables
double	Problem_Parameter::Re =1.0
	Reynolds number. More...
double	Problem_Parameter::St = 1.0
	Strouhal number. More...
double	Problem_Parameter::Density_ratio = 1.0
	Density ratio (Solid density / Fluid density) More...
double	Problem_Parameter::A = 0.25
	Initial axis of the elliptical solid in x-direction. More...
double	Problem_Parameter::B = 0.5
double	Problem_Parameter::Nu =0.3
	Pseudo-solid (mesh) Poisson ratio. More...
double	Problem_Parameter::Lambda_sq =0.0
ConstitutiveLaw *	Problem_Parameter::Constitutive_law_pt
	Constitutive law used to determine the mesh deformation. More...

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

Driver code for immersed ellipse problem.

{
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 // Define possible command line arguments and parse the ones that
 // were actually specified
 
 // Validation?
 CommandLineArgs::specify_command_line_flag("--validation");
 
 // Parse command line
 CommandLineArgs::parse_and_assign(); 
 
 // Doc what has actually been specified on the command line
 CommandLineArgs::doc_specified_flags();
  
 // Create problem in initial configuration
 UnstructuredImmersedEllipseProblem<
 ProjectableTaylorHoodElement<MyTaylorHoodElement> > problem;  
 
 //Initially ensure that the nodal positions are consistent with 
 //their weak imposition
 problem.solve_for_consistent_nodal_positions();
 
 // Initialise timestepper
 double dt=0.05;
 problem.initialise_dt(dt);
 
 // Perform impulsive start
 problem.assign_initial_values_impulsive();
 
 // Output initial conditions
 problem.doc_solution();
 
 // Solve problem a few times on given mesh
 unsigned nstep=3;
 for (unsigned i=0;i<nstep;i++)
  {
   // Solve the problem
   problem.unsteady_newton_solve(dt);    
   problem.doc_solution();
  }
 
 // Now do a couple of adaptations
 unsigned ncycle=200;
 if (CommandLineArgs::command_line_flag_has_been_set("--validation"))
  {
   ncycle=1;
   oomph_info << "Only doing one cycle during validation\n";
  }
 
 for (unsigned j=0;j<ncycle;j++)
  {       
   // Adapt the problem
   problem.adapt();
 
   //Solve problem a few times
   for (unsigned i=0;i<nstep;i++)
    {     
     // Solve the problem
     problem.unsteady_newton_solve(dt);
     problem.doc_solution();
    }
  }
 
} //end of main

References oomph::CommandLineArgs::command_line_flag_has_been_set(), oomph::CommandLineArgs::doc_specified_flags(), i, j, oomph::oomph_info, oomph::CommandLineArgs::parse_and_assign(), problem, Flag_definition::setup(), and oomph::CommandLineArgs::specify_command_line_flag().