refineable_two_layer_soluble_surfactant.cc File Reference

#include "generic.h"
#include "double_buoyant_navier_stokes_elements.h"
#include "fluid_interface.h"
#include "soluble_surfactant_transport_equations.h"
#include "constitutive.h"
#include "solid.h"
#include "meshes/rectangular_quadmesh.h"

Classes
class	DependentPositionPointElement
class	ElasticRefineableTwoLayerMesh< ELEMENT >
class	SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >

Namespaces
	oomph
	DRAIG: Change all instances of (SPATIAL_DIM) to (DIM-1).
	oomph::Control_Parameters
	oomph::Global_Physical_Variables
	Namespace for the physical parameters in the problem.

Functions
void	oomph::Global_Physical_Variables::hydrostatic_pressure_outlet_upper (const double &time, const Vector< double > &x, const Vector< double > &n, Vector< double > &traction)
void	oomph::Global_Physical_Variables::hydrostatic_pressure_inlet_upper (const double &time, const Vector< double > &x, const Vector< double > &n, Vector< double > &traction)
	Function that prescribes hydrostatic pressure field at the inlet. More...
void	oomph::Global_Physical_Variables::hydrostatic_pressure_outlet_lower (const double &time, const Vector< double > &x, const Vector< double > &n, Vector< double > &traction)
void	oomph::Global_Physical_Variables::hydrostatic_pressure_inlet_lower (const double &time, const Vector< double > &x, const Vector< double > &n, Vector< double > &traction)
	Function that prescribes hydrostatic pressure field at the inlet. More...
void	oomph::Global_Physical_Variables::wall_unit_normal_inlet_fct (const Vector< double > &x, Vector< double > &normal)
	Function that specifies the wall unit normal at the inlet. More...
void	oomph::Global_Physical_Variables::wall_unit_normal_outlet_fct (const Vector< double > &x, Vector< double > &normal)
	Function that specified the wall unit normal at the outlet. More...
int	main (int argc, char **argv)
	Driver code for 2D Boussinesq convection problem. More...

Variables
bool	oomph::Control_Parameters::Periodic_BCs = true
double	oomph::Global_Physical_Variables::L = 28.0
	Geometry. More...
double	oomph::Global_Physical_Variables::R = 1.0
	Fluid property Ratios. More...
double	oomph::Global_Physical_Variables::M = 0.5
double	oomph::Global_Physical_Variables::H0 = 0.2
double	oomph::Global_Physical_Variables::Ha = H0/5.0
double	oomph::Global_Physical_Variables::Re = 0.0
	Hydrodynamic Parameters. More...
double	oomph::Global_Physical_Variables::Ca = 0.001
double	oomph::Global_Physical_Variables::ReInvFr = 0.0
double	oomph::Global_Physical_Variables::Ma = 10.0
	Surfactant Parameters. More...
double	oomph::Global_Physical_Variables::Beta_s =Ca*Ma
	Surface Elasticity number (Capillary number x Marangoni number) More...
double	oomph::Global_Physical_Variables::Pe_s = 10.0
	Surface Peclet number. More...
double	oomph::Global_Physical_Variables::Pe_b = 10.0
	Bulk Peclet number. More...
double	oomph::Global_Physical_Variables::Pe_m = 10.0
	Micelle Pelect number. More...
double	oomph::Global_Physical_Variables::Biot = 0.1
	Solubility Parameters. More...
double	oomph::Global_Physical_Variables::K_b = 3.0
	The ratio of adsorption-desorption times. More...
double	oomph::Global_Physical_Variables::Beta_b = 1.0
double	oomph::Global_Physical_Variables::K_m = 1.0
double	oomph::Global_Physical_Variables::N = 10.0
	Power of the concentration in bulk -> micelle flux expression. More...
double	oomph::Global_Physical_Variables::Delta_P = 1.0
	The imposed pressure gradient. More...
Vector< double >	oomph::Global_Physical_Variables::Tau (2, 1.0)
Vector< double >	oomph::Global_Physical_Variables::D (2, 1.0)
Vector< double >	oomph::Global_Physical_Variables::Direction_of_gravity (2)
	Gravity vector, will be set in the main code. More...
double	oomph::Global_Physical_Variables::Nu = 0.1
	Pseudo-solid Poisson ratio. More...
double	oomph::Global_Physical_Variables::Inlet_Angle = 2.0*atan(1.0)
Vector< double >	oomph::Global_Physical_Variables::Wall_normal
	Direction of the wall normal vector (at the inlet) 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;
 
 //Set the diffusivities (inverse peclect numbers)
 Global_Physical_Variables::D[0] = 1.0/Global_Physical_Variables::Pe_b;
 Global_Physical_Variables::D[1] = 1.0/Global_Physical_Variables::Pe_m;
 
 Global_Physical_Variables::Wall_normal.resize(2);
 Global_Physical_Variables::Wall_normal[0] = -1.0;
 Global_Physical_Variables::Wall_normal[1] = 0.0;
 
 
 //Construct our problem
 SurfactantProblem<Hijacked<
                    RefineablePseudoSolidNodeUpdateElement<RefineableDoubleBuoyantQCrouzeixRaviartElement<2>, RefineableQPVDElement<2,3> > >,
           ElasticLineSolubleSurfactantTransportInterfaceElement<
                    RefineablePseudoSolidNodeUpdateElement<RefineableDoubleBuoyantQCrouzeixRaviartElement<2>, RefineableQPVDElement<2,3> > > > 
                   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.deform_interface(Global_Physical_Variables::Ha,1);
 problem.doc_solution(trace);
 problem.assign_initial_values_impulsive(dt);
 
 //Set the number of timesteps to our default value
 unsigned n_steps = 1000;
 
 //Set the freqency of refinement
 unsigned refine_after_n_steps = 5;
 
 //If we have a command line argument, perform fewer steps 
 //(used for self-test runs)
 if(argc > 1) {n_steps = 5; refine_after_n_steps=3;}
 
 //Perform n_steps timesteps
 for(unsigned i=0;i<n_steps;++i)
  {
   double dt_next = dt;
   //Spatial adapt every 5 steps with a fixed timestep
   if((i>0) && (i%refine_after_n_steps==0))
    {
      problem.unsteady_newton_solve(dt,1,false);
    }
    else
    {
     //Unsteady timestep (for fixed mesh)
      dt_next = problem.adaptive_unsteady_newton_solve(dt,1.0e-5);
    }
   dt = dt_next;
   problem.doc_solution(trace);
  }
 
} // end of main

References Global_Physical_Variables::D, Global_Physical_Variables::Direction_of_gravity, oomph::Global_Physical_Variables::Ha, i, oomph::Global_Physical_Variables::Pe_b, oomph::Global_Physical_Variables::Pe_m, problem, and Global_Physical_Variables::Wall_normal.