bretherton.cc File Reference

#include <algorithm>
#include "generic.h"
#include "navier_stokes.h"
#include "fluid_interface.h"
#include "meshes/bretherton_spine_mesh.h"

Classes
class	BrethertonElement< ELEMENT >
class	oomph::FaceGeometry< BrethertonElement< SpineElement< QCrouzeixRaviartElement< 2 > > > >
	Face geometry of the Bretherton 2D Crouzeix_Raviart spine elements. More...
class	oomph::FaceGeometry< BrethertonElement< SpineElement< QTaylorHoodElement< 2 > > > >
	Face geometry of the Bretherton 2D Taylor Hood spine elements. More...
class	oomph::FaceGeometry< FaceGeometry< BrethertonElement< SpineElement< QCrouzeixRaviartElement< 2 > > > > >
class	oomph::FaceGeometry< FaceGeometry< BrethertonElement< SpineElement< QTaylorHoodElement< 2 > > > > >
class	BrethertonProblem< ELEMENT >
	Bretherton problem. More...

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

Functions
Vector< double >	Global_Physical_Variables::G (2)
	Direction of gravity. More...
void	Global_Physical_Variables::inflow (const Vector< double > &x, Vector< double > &veloc)
int	main (int argc, char *argv[])

Variables
double *	Global_Physical_Variables::H_lo_pt
	Pointer to film thickness at outflow on the lower wall. More...
double *	Global_Physical_Variables::H_up_pt
	Pointer to film thickness at outflow on the upper wall. More...
double *	Global_Physical_Variables::Y_lo_pt
	Pointer to y-position at inflow on the lower wall. More...
double *	Global_Physical_Variables::Y_up_pt
	Pointer to y-position at inflow on the upper wall. 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 a single step.

{
 
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 // Set physical parameters:
 
 //Set direction of gravity: Vertically downwards
 Global_Physical_Variables::G[0] = 0.0;
 Global_Physical_Variables::G[1] = -1.0;
 
 // Womersley number = Reynolds number (St = 1)
 Global_Physical_Variables::ReSt = 0.0;
 Global_Physical_Variables::Re = Global_Physical_Variables::ReSt;
 
 // The Capillary number
 Global_Physical_Variables::Ca = 0.05;
 
 // Re/Fr -- a measure of gravity...
 Global_Physical_Variables::ReInvFr = 0.0;
 
 //Set up the problem
 BrethertonProblem<BrethertonElement<SpineElement<QCrouzeixRaviartElement<2> > > > problem;
 
 // Self test:
 problem.self_test();
 
 // Number of steps: 
 unsigned nstep;
 if (CommandLineArgs::Argc>1)
  {
   // Validation run: Just one step
   nstep=2;
  }
 else
  {
   // Full run otherwise
   nstep=30;
  }
 
 // Run the parameter study: Perform nstep steps
 problem.parameter_study(nstep);
 
}

References oomph::CommandLineArgs::Argc, Global_Physical_Variables::Ca, Global_Physical_Variables::G, problem, Global_Physical_Variables::Re, Global_Physical_Variables::ReInvFr, Global_Physical_Variables::ReSt, and Flag_definition::setup().