surfactant.cc File Reference

#include "generic.h"
#include "advection_diffusion.h"
#include "navier_stokes.h"
#include "multi_physics.h"
#include "fluid_interface.h"
#include "meshes/single_layer_spine_mesh.h"

Classes
class	oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >
class	ConvectionProblem< NST_ELEMENT, AD_ELEMENT >

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

Functions
int	main (int argc, char **argv)
	Driver code for 2D Boussinesq convection problem. More...

Variables
double	Global_Physical_Variables::Scaled_Bond = 0.0
double	Global_Physical_Variables::Biot = 0.0
	Biot number. More...
double	Global_Physical_Variables::Marangoni = 125.0
	Marangoni number. More...
double	Global_Physical_Variables::Capillary = 0.0045
	Capillary number. More...
double	Global_Physical_Variables::Beta = 100.0
	Surface Elasticity number. More...
double	Global_Physical_Variables::Peclet_S = 1.0
	Surface Peclet number. More...
double	Global_Physical_Variables::Peclet_St_S = 100.0
	\shorT Sufrace Peclet number multiplied by Strouhal number More...
Vector< double >	Global_Physical_Variables::Direction_of_gravity (2)
	Gravity vector. More...

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

Driver code for 2D Boussinesq convection problem.

{
 ofstream trace("RESLT/trace.dat");
 
 // Set the direction of gravity
 Global_Physical_Variables::Direction_of_gravity[0] = 0.0;
 Global_Physical_Variables::Direction_of_gravity[1] = -1.0;
 
 //Construct our problem
 ConvectionProblem<SpineElement<BuoyantQCrouzeixRaviartElement<2> >,
  SpineLineMarangoniSurfactantFluidInterfaceElement<BuoyantQCrouzeixRaviartElement<2> > > 
problem;
 
 // Apply the boundary condition at time zero
 problem.set_boundary_conditions(0.0);
 
 //Perform a single steady Newton solve
 problem.steady_newton_solve();
 
 //Document the solution
 problem.doc_solution(trace);
 
 //Now release the interface for real fun
 //problem.unpin_surface();
 
 //Set the timestep
 double dt = 0.1;
 
 //Initialise the value of the timestep and set initial values 
 //of previous time levels assuming an impulsive start.
 problem.assign_initial_values_impulsive(dt);
 
 //Set the number of timesteps to our default value
 unsigned n_steps = 1000;
 
 //If we have a command line argument, perform fewer steps 
 //(used for self-test runs)
 if(argc > 1) {n_steps = 5;}
 
 //Perform n_steps timesteps
 for(unsigned i=0;i<n_steps;++i)
  {
   problem.unsteady_newton_solve(dt);
   problem.doc_solution(trace);
  }
 
} // end of main

References Global_Physical_Variables::Direction_of_gravity, i, and problem.