three_d_breth.cc File Reference

#include <algorithm>
#include <iostream>
#include <cmath>
#include <map>
#include "generic.h"
#include "navier_stokes.h"
#include "fluid_interface.h"
#include "merge_meshes.h"
#include "st_mesh.h"
#include "extra_elements.h"

Classes
class	ThreeDimBethertonProblem< ELEMENT >
	Brethreton problem in square tube domain. More...

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
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...
int	main (int argc, char *argv[])

Variables
double	Global_Physical_Variables::Height = 1.0
double	Global_Physical_Variables::Length_tip = 0.50
double	Global_Physical_Variables::Length_liq = 1.0
double	Global_Physical_Variables::Length_can = 1.5
unsigned	Global_Physical_Variables::Ncan = 2
unsigned	Global_Physical_Variables::Ntip = 2
unsigned	Global_Physical_Variables::Nliq = 2
double	Global_Physical_Variables::Rat_press = 2.87
Vector< double >	Global_Physical_Variables::G (3) = 1.0
	Direction of gravity. More...

Function Documentation

main()

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

Driver for RectangularDrivenCavity test problem – test drive with two different types of element.

{
 // Doc info object
 DocInfo doc_info;
 
 // Set output directory
 char name_dir[20];
 sprintf(name_dir,"RESLT");
 
 doc_info.set_directory(name_dir);
 doc_info.number()=0;
 
 // Direction of gravity
 Global_Physical_Variables::G[0] = 0.0;
 Global_Physical_Variables::G[1] = 0.0;
 Global_Physical_Variables::G[2] = -1.0;
 
 //Construct the problem
 ThreeDimBethertonProblem<SpineElement<QCrouzeixRaviartElement<3> >  >
  problem;
 
 //Solve and document without gravity
 problem.newton_solve();
 problem.doc_solution(doc_info);
 
 //Crank up the Bond number
 for(unsigned i=0;i<2;i++) 
  {
   doc_info.number()++;
   using namespace Global_Physical_Variables;
   Bo += 0.1;
   ReInvFr = Bo/Ca;
   problem.newton_solve();
   problem.doc_solution(doc_info);
  }
 
}

References Global_Physical_Variables::Bo, Global_Physical_Variables::Ca, Global_Physical_Variables::G, i, oomph::DocInfo::number(), problem, GlobalPhysicalVariables::ReInvFr, and oomph::DocInfo::set_directory().