inclined_plane.cc File Reference

#include <iostream>
#include <fstream>
#include <string>
#include <cmath>
#include "generic.h"
#include "navier_stokes.h"
#include "solid.h"
#include "fluid_interface.h"
#include "meshes/simple_rectangular_quadmesh.h"

Classes
class	InclinedPlaneProblem< ELEMENT, INTERFACE_ELEMENT >
class	SpineInclinedPlaneMesh< ELEMENT >
	Create a spine mesh for the problem. More...
class	SpineInclinedPlaneProblem< ELEMENT, TIMESTEPPER >
class	ElasticInclinedPlaneMesh< ELEMENT >
	Create an Elastic mesh for the problem. More...
class	ElasticInclinedPlaneProblem< ELEMENT, TIMESTEPPER >

Namespaces
	Global_Physical_Variables
	Global variables.

Macros
#define	FLUID_ELEMENT   QTaylorHoodElement<2>

Functions
Vector< double >	Global_Physical_Variables::G (2, 0.0)
	The Vector direction of gravity, set in main() More...
void	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	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...
void	Global_Physical_Variables::hydrostatic_pressure_outlet (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::hydrostatic_pressure_inlet (const double &time, const Vector< double > &x, const Vector< double > &n, Vector< double > &traction)
	Function that prescribes hydrostatic pressure field at the inlet. More...
int	main (int argc, char **argv)

Variables
double	Global_Physical_Variables::N_wave = 3
	Set the number of waves desired in the domain. More...
double	Global_Physical_Variables::Inlet_Angle = 2.0*atan(1.0)
	The contact angle that is imposed at the inlet (pi) More...

Macro Definition Documentation

FLUID_ELEMENT

#define FLUID_ELEMENT   QTaylorHoodElement<2>

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

{
 using namespace Global_Physical_Variables;
 
 //Set the constitutive law for the mesh deformation
 Constitutive_law_pt  = new GeneralisedHookean(&Global_Physical_Variables::Nu);
 
#ifdef CR_ELEMENT
#define FLUID_ELEMENT QCrouzeixRaviartElement<2>
#else
#define FLUID_ELEMENT QTaylorHoodElement<2>
#endif 
 
 //Initialise physical parameters
 //Scale Reynolds number to be independent of alpha.
 Re = 4.0/sin(Alpha); 
 
 //Set the direction of gravity
 G[0] = sin(Alpha);
 G[1] = -cos(Alpha);
 
 //The wall normal to the inlet is in the negative x direction
 Wall_normal.resize(2);
 Wall_normal[0] = -1.0;
 Wall_normal[1] = 0.0;
 
 //Spine problem
 {
  //Create the problem
  SpineInclinedPlaneProblem<SpineElement<FLUID_ELEMENT >, BDF<2> > 
   problem(30,4,Global_Physical_Variables::Length);
  
  //Solve the steady problem
  problem.solve_steady();
  
  //Prepare the problem for timestepping 
  //(assume that it's been at the flat-film solution for all previous time)
  double dt = 0.1;
  problem.assign_initial_values_impulsive(dt);
  
  //Timestep it 
  problem.timestep(dt,2);
 } //End of spine problem
 
 
 //Elastic problem
 {
  //Create the problem
  ElasticInclinedPlaneProblem<
   PseudoSolidNodeUpdateElement<FLUID_ELEMENT,QPVDElement<2,3> >, BDF<2> > 
     problem(30,4,Global_Physical_Variables::Length);
  
  //Solve the steady problem
  problem.solve_steady();
  
  //Prepare the problem for timestepping 
  //(assume that it's been at the flat-film solution for all previous time)
  double dt = 0.1;
  problem.assign_initial_values_impulsive(dt);
 
  //Timestep it
  problem.timestep(dt,2);
 } //End of elastic problem
}

References TanhSolnForAdvectionDiffusion::Alpha, Constitutive::Constitutive_law_pt, cos(), G, Global_Physical_Variables::Length, Global_Physical_Variables::Nu, problem, GlobalPhysicalVariables::Re, sin(), and Global_Physical_Variables::Wall_normal.