rayleigh_instability_soluble_surfactant.cc File Reference

#include "generic.h"
#include "axisym_navier_stokes.h"
#include "fluid_interface.h"
#include "solid.h"
#include "constitutive.h"
#include "meshes/rectangular_quadmesh.h"
#include "axisymmetric_advection_navier_stokes_elements.h"

Classes
class	ElasticAxisymmetricSolubleSurfactantTransportInterfaceElement< ELEMENT >
class	InterfaceProblem< ELEMENT, TIMESTEPPER >
	Single axisymmetric fluid interface problem in rectangular domain. More...

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
int	main (int argc, char *argv[])
	Driver code for single fluid axisymmetric horizontal interface problem. More...

Variables
double	Global_Physical_Variables::Alpha_absorption = 1.0
	Alpha for absorption kinetics. More...
double	Global_Physical_Variables::K = 0.01
	The ratio of adsorption-desorption times. More...
double	Global_Physical_Variables::Pe = 10000.0
	Peclet number. More...
double	Global_Physical_Variables::PeSt = Pe
double	Global_Physical_Variables::Pe_reference_scale =1.0
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 single fluid axisymmetric horizontal interface problem.

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

Maximum time

Duration of timestep

{
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 double t_max = 1000.0;
 
 const double dt = 0.1;
 
 // If we are doing validation run, use smaller number of timesteps
 if(CommandLineArgs::Argc>1) 
  { 
   t_max = 0.5; 
  }
 
 // Number of elements in radial (r) direction
 const unsigned n_r = 10;
   
 // Number of elements in axial (z) direction
 const unsigned n_z = 80;
 
 // Height of domain
 const double l_z = MathematicalConstants::Pi/Global_Physical_Variables::Alpha;
 
 // 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 up the spine test problem with AxisymmetricQCrouzeixRaviartElements,
 // using the BDF<2> timestepper
 InterfaceProblem<PseudoSolidNodeUpdateElement<AxisymmetricQAdvectionCrouzeixRaviartElement, QPVDElement<2,3> > ,BDF<2> >
  problem(n_r,n_z,l_z);
 
 // Run the unsteady simulation
 problem.unsteady_run(t_max,dt);
} // End of main

References Global_Physical_Variables::Alpha, oomph::CommandLineArgs::Argc, Global_Physical_Variables::G, oomph::MathematicalConstants::Pi, problem, and oomph::CommandLineArgs::setup().