3d_rayleigh_instability_surfactant.cc File Reference

#include "generic.h"
#include "navier_stokes.h"
#include "fluid_interface.h"
#include "meshes/single_layer_cubic_spine_mesh.h"

Classes
class	ElementCmp
	Function-type-object to perform comparison of elements. More...
class	oomph::AnnularSpineMesh< ELEMENT >
class	InterfaceProblem< ELEMENT, TIMESTEPPER >
	Single axisymmetric fluid interface problem in rectangular domain. More...

Namespaces
	Global_Physical_Variables
	Global variables.
	oomph
	DRAIG: Change all instances of (SPATIAL_DIM) to (DIM-1).

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

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

Function Documentation

main()

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

Driver code for unsteady two-layer fluid problem. If there are any command line arguments, we regard this as a validation run and perform only two steps.

{
 
 // Set physical parameters:
 
 //Set direction of gravity: Vertically downwards
 Global_Physical_Variables::G[0] = 0.0;
 Global_Physical_Variables::G[1] = 0.0;
 Global_Physical_Variables::G[2] = 0.0;
 
 //Set the film thickness
 const double h = Global_Physical_Variables::Film_Thickness;
 
 //Axial lngth of domain
 double r_max = 1.0;
 double r_min = r_max - h;
 
 const double pi = MathematicalConstants::Pi;
 
 double alpha = Global_Physical_Variables::Alpha;
 
 double l_y = pi/alpha;
 
 double theta_max = 0.5*pi;
  
 // Number of elements in r and theta direction
 unsigned n_r = 4;
 
 unsigned n_theta = 4;
 
 // Number of elements in axial direction
 unsigned n_y = 20;
 
 
 {
  //Set up the spine test problem
  //The minimum values are all assumed to be zero.
  InterfaceProblem<SpineElement<QCrouzeixRaviartElement<3> >,BDF<2> >
   problem(n_r,n_y,n_theta,r_min,r_max,l_y,theta_max);
  
  // Number of steps: 
  unsigned nstep;
  if(argc > 1)
   {
    // Validation run: Just five steps
    nstep=5;
   }
  else
   {
    // Full run otherwise
    nstep=1000;
   }
  
  // Run the unsteady simulation
  problem.unsteady_run(nstep);
 }
} // End of main

References alpha, Global_Physical_Variables::Alpha, Global_Physical_Variables::Film_Thickness, Global_Physical_Variables::G, constants::pi, BiharmonicTestFunctions2::Pi, problem, BiharmonicTestFunctions2::r_max, and BiharmonicTestFunctions2::r_min.