3d_free_surface.cc File Reference

#include <iostream>
#include <fstream>
#include <cmath>
#include <typeinfo>
#include <algorithm>
#include <cstdio>
#include <complex>
#include "generic.h"
#include "navier_stokes.h"
#include "fluid_interface.h"
#include "constitutive.h"
#include "solid.h"
#include "meshes/single_layer_cubic_spine_mesh.h"

Classes
class	InterfaceProblem< ELEMENT, TIMESTEPPER >
	Single axisymmetric fluid interface problem in rectangular domain. More...
class	ElasticSimpleCubicMesh< ELEMENT >
	Elastic version of the simple cubic mesh. More...
class	ElasticInterfaceProblem< ELEMENT, TIMESTEPPER >
	Elastic version of Single fluid interface problem. More...

Namespaces
	Global_Physical_Variables
	Global variables.

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

Variables
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] = -1.0;
 
 // Womersley number = Reynolds number (St = 1)
 Global_Physical_Variables::ReSt = 5.0;
 Global_Physical_Variables::Re = Global_Physical_Variables::ReSt;
 
 // The Capillary number
 Global_Physical_Variables::Ca = 0.01;
 
 // Re/Fr 
 Global_Physical_Variables::ReInvFr = 0.0;
 
 //Axial lngth of domain
 double lx = 2.0;
 
 double ly = 2.0;
   
 // Number of elements in x- and y-direction
 unsigned nx=5;
 
 unsigned ny = 5;
 
 // Number of elements in layer
 unsigned nz=7;
 
 // Height of layer
 double h=1.0;
 
 {
  //Set up the spine test problem
  InterfaceProblem<SpineElement<QCrouzeixRaviartElement<3> >,BDF<2> >
   problem(nx,ny,nz,lx,ly,h);
  
  // Number of steps: 
  unsigned nstep;
  if(argc > 1)
   {
    // Validation run: Just two steps
    nstep=2;
   }
  else
   {
    // Full run otherwise
    nstep=160;
   }
  
  // Run the unsteady simulation
  problem.unsteady_run(nstep);
 }
 
 //Set the constituive law
   Global_Physical_Variables::Constitutive_law_pt = 
  new GeneralisedHookean(&Global_Physical_Variables::Nu);
 
 {
  //Set up the spine test problem
  ElasticInterfaceProblem<PseudoSolidNodeUpdateElement
                          <QCrouzeixRaviartElement<3>,
                           QPVDElementWithPressure<3> >,BDF<2> >
   problem(nx,ny,nz,lx,ly,h);
  
  // Number of steps: 
  unsigned nstep;
  if(argc > 1)
   {
    // Validation run: Just two steps
    nstep=2;
   }
  else
   {
    // Full run otherwise
    nstep=160;
   }
  
  // Run the unsteady simulation
  problem.unsteady_run(nstep);
 }
 
 
}

References Global_Physical_Variables::Ca, Global_Physical_Variables::Constitutive_law_pt, Global_Physical_Variables::G, Mesh_Parameters::lx, Mesh_Parameters::ly, Global_Physical_Variables::Nu, Mesh_Parameters::nx, Mesh_Parameters::ny, Mesh_Parameters::nz, problem, Global_Physical_Variables::Re, Global_Physical_Variables::ReInvFr, and Global_Physical_Variables::ReSt.