spine_two_layer_interface.cc File Reference

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

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

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
Vector< double >	Global_Physical_Variables::G (2)
	Direction of gravity. More...
int	main (int argc, char *argv[])
	Driver code for two-dimensional two fluid interface problem. More...

Function Documentation

main()

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

Driver code for two-dimensional two fluid interface problem.

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

Maximum time

Duration of timestep

{
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 // Compute the Womersley number
 Global_Physical_Variables::ReSt =
  Global_Physical_Variables::Re*Global_Physical_Variables::St;
 
 double t_max = 0.6;
 
 const double dt = 0.0025;
 
 // If we are doing validation run, use smaller number of timesteps
 if(CommandLineArgs::Argc>1) { t_max = 0.005; }
 
 // Number of elements in x direction
 const unsigned n_x = 12;
   
 // Number of elements in y direction in lower fluid (fluid 1)
 const unsigned n_y1 = 12;
   
 // Number of elements in y direction in upper fluid (fluid 2)
 const unsigned n_y2 = 12;
 
 // Width of domain
 const double l_x = 1.0;
 
 // Height of lower fluid layer
 const double h1 = 1.0;
   
 // Height of upper fluid layer
 const double h2 = 1.0;
 
 // Set direction of gravity (vertically downwards)
 Global_Physical_Variables::G[0] = 0.0;
 Global_Physical_Variables::G[1] = -1.0;
 
 // Set up the spine test problem with QCrouzeixRaviartElements,
 // using the BDF<2> timestepper
 InterfaceProblem<SpineElement<QCrouzeixRaviartElement<2> >,BDF<2> >
  problem(n_x,n_y1,n_y2,l_x,h1,h2);
   
 // Run the unsteady simulation
 problem.unsteady_run(t_max,dt);
 
} // End of main

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