marangoni_convection_box.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	ComplexGreater< T >
	Driver for a multi-physics problem that couples the Navier–Stokes. More...
class	oomph::SpineVolumeConstraintPointElement< ELEMENT >
class	oomph::FixedVolumeSpineLineMarangoniFluidInterfaceElement< 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
void	Global_Physical_Variables::wall_unit_normal_left_fct (const Vector< double > &x, Vector< double > &normal)
	Function that specifies the wall unit normal. More...
void	Global_Physical_Variables::wall_unit_normal_right_fct (const Vector< double > &x, Vector< double > &normal)
	Function that specifies the wall unit normal. More...
int	main (int argc, char **argv)
	Driver code for 2D Boussinesq convection problem. More...

Variables
const double	Global_Physical_Variables::Pi =MathematicalConstants::Pi
	Set the value of Pi. More...
double	Global_Physical_Variables::Volume = 1.2
	The volume of the domain. More...
double	Global_Physical_Variables::Angle = 0.5*Pi
	The contact angle. 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.

{
 // 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> >,
  FixedVolumeSpineLineMarangoniFluidInterfaceElement<SpineElement<BuoyantQCrouzeixRaviartElement<2> > > > problem;
 
 //Set "interesting" Marangoni and Capillary numbers 
 Global_Physical_Variables::Marangoni = 200.0;
 Global_Physical_Variables::Capillary = 2.75e-3;
 
 // 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();
 
 //Decrease the contact angle
 for(unsigned i=0;i<3;i++)
  {
   Global_Physical_Variables::Angle -= 5*Global_Physical_Variables::Pi/180;
   problem.steady_newton_solve();
   problem.doc_solution();
  }
 
 //Now we timestep the system
 //Note that we are not taking any contact line dynamics into
 //account, we are merely evolving between steady solutions
 {
  problem.switch_boundary_conditions();
 
  //Set the timestep
  double dt = 0.5;
  
  //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 = 200;
  
  //If we have a command line argument, perform fewer steps 
  //(used for self-test runs)
  if(argc > 1) {n_steps = 2;}
  
  //Perform n_steps timesteps
  for(unsigned i=0;i<n_steps;++i)
   {
    problem.unsteady_newton_solve(dt);
    problem.doc_solution();
   }
 }
 
} // end of main

References Global_Physical_Variables::Angle, Global_Physical_Variables::Capillary, oomph::Global_Physical_Variables::Direction_of_gravity, i, Global_Physical_Variables::Marangoni, Global_Physical_Variables::Pi, and problem.