elastic_bretherton.cc File Reference

#include <iostream>
#include <fstream>
#include <cmath>
#include <typeinfo>
#include <algorithm>
#include <cstdio>
#include <complex>
#include <set>
#include "generic.h"
#include "navier_stokes.h"
#include "fluid_interface.h"
#include "beam.h"
#include "meshes/bretherton_spine_mesh.h"
#include "meshes/one_d_lagrangian_mesh.h"

Classes
class	ElementCmp
	Function-type-object to perform comparison of elements. More...
class	SpineGravityTractionElement< ELEMENT >
class	FixSpineHeightElement
	Point element that is used to set the bubble pressure. More...
class	FluxConstraint
	Generalised element used to set the pressure gradient. More...
class	AirwayReopeningProblem< ELEMENT >
	Bretherton problem. More...

Namespaces
	Global_Physical_Variables
	Global variables.
	No_Slip

Functions
Vector< double >	Global_Physical_Variables::G (2)
	Direction of gravity. More...
void	Global_Physical_Variables::hydrostatic_pressure (const double &time, const Vector< double > &x, const Vector< double > &n, Vector< double > &traction)
	Function that prescribes the hydrostatic pressure field at the outlet. More...
void	Global_Physical_Variables::spring_load (const Vector< double > &xi, const Vector< double > &x, const Vector< double > &N, Vector< double > &load)
	Load function for the wall elements. More...
void	Global_Physical_Variables::spring_load_lower (const Vector< double > &xi, const Vector< double > &x, const Vector< double > &N, Vector< double > &load)
	Load function for the wall elements. More...
void	No_Slip::no_slip_condition_first (Node *node_pt)
	Function that is used to set and update the no-slip boundary condition. More...
void	No_Slip::no_slip_condition_second (Node *node_pt)
int	main (int argc, char *argv[])

Variables
double	Global_Physical_Variables::ReCa =0.0
	Reynolds divided by Capillary number. More...
double	Global_Physical_Variables::Bo = 0.0
	Bond number. More...
double	Global_Physical_Variables::Pext = 0.0
	External Pressure. More...
double	Global_Physical_Variables::Gamma = 1.0e-7
double	Global_Physical_Variables::T = 0.0
	Period of oscillations. More...
double	Global_Physical_Variables::Kstiff = 0.5*50.0e-7
	Spring stiffness. More...
double	Global_Physical_Variables::Kseparation = H
	Huge stiffness for hitting. More...
double	Global_Physical_Variables::Ktable = 0.0001
	Huge stiffness for the table. More...
double	Global_Physical_Variables::Tube_width = 10.0
	Natural width of the open tube. More...
double	Global_Physical_Variables::Table_position = 10.78
	Position of the table. More...
double	Global_Physical_Variables::Q =1.0e-2
	Ratio of scales. More...
double	Global_Physical_Variables::Rest_length_linear =1.0
	Rest length of the linear springs. More...
bool	Global_Physical_Variables::Non_linear_springs = false
	Do we include the non-linear springs. More...
GeomObject *	Global_Physical_Variables::Upper_wall_pt
	Pointer to the upper wall. More...
GeomObject *	Global_Physical_Variables::Lower_wall_pt
	Pointer to the lower wall. More...
map< double, pair< GeomObject *, Vector< double > > >	Global_Physical_Variables::upper_map
	Upper wall map. More...
map< double, pair< GeomObject *, Vector< double > > >	Global_Physical_Variables::lower_map
	Lower wall map. 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:
 using namespace Global_Physical_Variables;
 
 //Set direction of gravity: Vertically downwards
 G[0] = 0.0;
 G[1] = -1.0;
 
 // Womersley number = Reynolds number (St = 1)
 ReSt = 0.0;
 Re = Global_Physical_Variables::ReSt;
 
 // The Capillary number
 Ca = 0.05;
 
 // The Bond number
 Bo = 0.0;
 
 // Re/Fr -- a measure of gravity...
 ReInvFr = Bo/Ca;
 
 // The ratio of scales
 Q = Ca*Gamma;
 
 // The value of sigma0
 T = 100.0*Gamma/H;
 
 //Set up the problem
 AirwayReopeningProblem<Hijacked<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=200;
  }
 
 // Run the parameter study: Perform nstep steps
 problem.parameter_study(nstep,false);
}

References oomph::CommandLineArgs::Argc, Global_Physical_Variables::Bo, Global_Physical_Variables::Ca, G, Global_Physical_Variables::Gamma, H, problem, Global_Physical_Variables::Q, GlobalPhysicalVariables::Re, GlobalPhysicalVariables::ReInvFr, GlobalPhysicalVariables::ReSt, Global_Physical_Variables::ReSt, and Flag_definition::setup().