oomph::Global_Physical_Variables Namespace Reference

Namespace for the physical parameters in the problem. More...

Functions
void	hydrostatic_pressure_outlet_upper (const double &time, const Vector< double > &x, const Vector< double > &n, Vector< double > &traction)
void	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	hydrostatic_pressure_outlet_lower (const double &time, const Vector< double > &x, const Vector< double > &n, Vector< double > &traction)
void	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	wall_unit_normal_inlet_fct (const Vector< double > &x, Vector< double > &normal)
	Function that specifies the wall unit normal at the inlet. More...
void	wall_unit_normal_outlet_fct (const Vector< double > &x, Vector< double > &normal)
	Function that specified the wall unit normal at the outlet. More...

Variables
double	L = 28.0
	Geometry. More...
double	R = 1.0
	Fluid property Ratios. More...
double	M = 0.5
double	H0 = 0.2
double	Ha = H0/5.0
double	Re = 0.0
	Hydrodynamic Parameters. More...
double	Ca = 0.001
double	ReInvFr = 0.0
double	Ma = 10.0
	Surfactant Parameters. More...
double	Beta_s =Ca*Ma
	Surface Elasticity number (Capillary number x Marangoni number) More...
double	Pe_s = 10.0
	Surface Peclet number. More...
double	Pe_b = 10.0
	Bulk Peclet number. More...
double	Pe_m = 10.0
	Micelle Pelect number. More...
double	Biot = 0.1
	Solubility Parameters. More...
double	K_b = 3.0
	The ratio of adsorption-desorption times. More...
double	Beta_b = 1.0
double	K_m = 1.0
double	N = 10.0
	Power of the concentration in bulk -> micelle flux expression. More...
double	Delta_P = 1.0
	The imposed pressure gradient. More...
Vector< double >	Tau (2, 1.0)
Vector< double >	D (2, 1.0)
Vector< double >	Direction_of_gravity (2)
	Gravity vector, will be set in the main code. More...
double	Nu = 0.1
	Pseudo-solid Poisson ratio. More...
double	Inlet_Angle = 2.0*atan(1.0)
Vector< double >	Wall_normal
	Direction of the wall normal vector (at the inlet) More...

Detailed Description

Namespace for the physical parameters in the problem.

Function Documentation

hydrostatic_pressure_inlet_lower()

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.

 {
  traction[0] = Delta_P +
   -ReInvFr*Direction_of_gravity[1]*(R*(1.0 - H0) + H0 - x[1]);
   traction[1] = 0.0;
 }

References Delta_P, Direction_of_gravity, H0, R, ReInvFr, and plotDoE::x.

hydrostatic_pressure_inlet_upper()

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.

 {
   traction[0] = Delta_P + -ReInvFr*R*Direction_of_gravity[1]*(1.0 - x[1]);
   traction[1] = 0.0;
 }

References Delta_P, Direction_of_gravity, R, ReInvFr, and plotDoE::x.

hydrostatic_pressure_outlet_lower()

void oomph::Global_Physical_Variables::hydrostatic_pressure_outlet_lower	(	const double &	time,
		const Vector< double > &	x,
		const Vector< double > &	n,
		Vector< double > &	traction
	)

Function that prescribes the hydrostatic pressure field at the outlet Must match pressure in lower fluid — This may be tricky if the interface is not pinned (i.e. we'll need to read out the interfacial position on the boundary). For now assume it's at H0.

 {
  traction[0] = ReInvFr*Direction_of_gravity[1]*(R*(1.0 - H0) + H0 - x[1]);
  traction[1] = 0.0;
 }

References Direction_of_gravity, H0, R, ReInvFr, and plotDoE::x.

hydrostatic_pressure_outlet_upper()

void oomph::Global_Physical_Variables::hydrostatic_pressure_outlet_upper	(	const double &	time,
		const Vector< double > &	x,
		const Vector< double > &	n,
		Vector< double > &	traction
	)

This next set of functions is only used if we do NOT have periodic conditions Function that prescribes the hydrostatic pressure field at the outlet Let's fix things so that the pressure at the top of the channel is zero.

 {
  traction[0] = ReInvFr*R*Direction_of_gravity[1]*(1.0 - x[1]);
  traction[1] = 0.0;
 }

References Direction_of_gravity, R, ReInvFr, and plotDoE::x.

wall_unit_normal_inlet_fct()

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.

 {
  normal=Wall_normal;
 }

References WallFunction::normal(), and Wall_normal.

wall_unit_normal_outlet_fct()

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.

 {
  //Set the normal
  normal = Wall_normal;
  //and flip the sign
  unsigned n_dim = normal.size();
  for(unsigned i=0;i<n_dim;++i) {normal[i] *= -1.0;}
 }

References i, WallFunction::normal(), and Wall_normal.

Variable Documentation

Beta_b

double oomph::Global_Physical_Variables::Beta_b = 1.0

Referenced by oomph::SolubleSurfactantTransportInterfaceElement::add_additional_residual_contributions_interface(), SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >::create_interface_elements(), and SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >::doc_solution().

Beta_s

double oomph::Global_Physical_Variables::Beta_s =Ca*Ma

Surface Elasticity number (Capillary number x Marangoni number)

Referenced by SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >::create_interface_elements().

Biot

double oomph::Global_Physical_Variables::Biot = 0.1

Solubility Parameters.

Biot number

Referenced by SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >::create_interface_elements().

Ca

double oomph::Global_Physical_Variables::Ca = 0.001

We do not need the Weber number, because we have non-dimensionalised pressure on the viscous scale and so multiplying the normal stress balanced by the Reynolds number gives a term in the Capillary number only (Ca Re = We). Capillary number (of which the results are independent for a pinned surface)

Referenced by SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >::create_interface_elements(), and InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem().

D

Vector< double > oomph::Global_Physical_Variables::D	(	2	,
		1.	0
	)

Diffusivity (should be 1/Pe_b, 1/Pe_m), which will be set in the main code

Delta_P

double oomph::Global_Physical_Variables::Delta_P = 1.0

The imposed pressure gradient.

Referenced by hydrostatic_pressure_inlet_lower(), and hydrostatic_pressure_inlet_upper().

Direction_of_gravity

Vector< double > oomph::Global_Physical_Variables::Direction_of_gravity

(

2

)

Gravity vector, will be set in the main code.

Referenced by ConvectionProblem< NST_ELEMENT, AD_ELEMENT >::ConvectionProblem(), DDConvectionProblem< NST_ELEMENT, AD_ELEMENT >::DDConvectionProblem(), hydrostatic_pressure_inlet_lower(), hydrostatic_pressure_inlet_upper(), hydrostatic_pressure_outlet_lower(), hydrostatic_pressure_outlet_upper(), main(), RefineableConvectionProblem< NST_ELEMENT, AD_ELEMENT >::RefineableConvectionProblem(), and RefineableDDConvectionProblem< NST_ELEMENT, AD_ELEMENT >::RefineableDDConvectionProblem().

H0

double oomph::Global_Physical_Variables::H0 = 0.2

Referenced by SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >::deform_interface(), SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >::doc_solution(), hydrostatic_pressure_inlet_lower(), hydrostatic_pressure_outlet_lower(), and SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >::SurfactantProblem().

Ha

double oomph::Global_Physical_Variables::Ha = H0/5.0

Referenced by main().

Inlet_Angle

double oomph::Global_Physical_Variables::Inlet_Angle = 2.0*atan(1.0)

K_b

double oomph::Global_Physical_Variables::K_b = 3.0

The ratio of adsorption-desorption times.

Referenced by SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >::create_interface_elements(), and SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >::set_boundary_conditions().

K_m

double oomph::Global_Physical_Variables::K_m = 1.0

Referenced by SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >::SurfactantProblem().

L

double oomph::Global_Physical_Variables::L = 28.0

Geometry.

Referenced by SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >::deform_interface(), and RefineablePorousChannelProblem< ELEMENT >::RefineablePorousChannelProblem().

M

double oomph::Global_Physical_Variables::M = 0.5

Ratio of viscosity in upper fluid to viscosity in lower fluid. Reynolds number etc. is based on viscosity in lower fluid.

Referenced by SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >::create_interface_elements(), and SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >::SurfactantProblem().

Ma

double oomph::Global_Physical_Variables::Ma = 10.0

Surfactant Parameters.

Marangoni number

Referenced by SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >::create_interface_elements(), oomph::SpineLineMarangoniFluidInterfaceElement< ELEMENT >::sigma(), oomph::FixedVolumeSpineLineMarangoniFluidInterfaceElement< ELEMENT >::sigma(), and oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::sigma().

N

double oomph::Global_Physical_Variables::N = 10.0

Power of the concentration in bulk -> micelle flux expression.

Nu

double oomph::Global_Physical_Variables::Nu = 0.1

Pseudo-solid Poisson ratio.

Referenced by main().

Pe_b

double oomph::Global_Physical_Variables::Pe_b = 10.0

Bulk Peclet number.

Referenced by main().

Pe_m

double oomph::Global_Physical_Variables::Pe_m = 10.0

Micelle Pelect number.

Referenced by main().

Pe_s

double oomph::Global_Physical_Variables::Pe_s = 10.0

Surface Peclet number.

Referenced by oomph::SpineAxisymmetricMarangoniSurfactantFluidInterfaceElement< ELEMENT >::add_additional_residual_contributions_interface(), oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::add_additional_residual_contributions_interface(), oomph::SurfactantTransportInterfaceElement::add_additional_residual_contributions_interface(), and SurfactantProblem< ELEMENT, INTERFACE_ELEMENT >::create_interface_elements().

R

double oomph::Global_Physical_Variables::R = 1.0

Fluid property Ratios.

Ratio of density in upper fluid to density in lower fluid. Reynolds number etc. is based on density in lower fluid.

Re

double oomph::Global_Physical_Variables::Re = 0.0

Hydrodynamic Parameters.

Reynolds number

Referenced by oomph::StreamfunctionProblem::header(), InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem(), main(), oomph::StreamfunctionProblem::my_output(), and RefineablePorousChannelProblem< ELEMENT >::RefineablePorousChannelProblem().

ReInvFr

double oomph::Global_Physical_Variables::ReInvFr = 0.0

In our non-dimensionalisation, we have a Reynolds number/ Froude number squared in (a material parameter that doesn't involve the velocity scale).

Referenced by hydrostatic_pressure_inlet_lower(), hydrostatic_pressure_inlet_upper(), hydrostatic_pressure_outlet_lower(), hydrostatic_pressure_outlet_upper(), and InterfaceProblem< ELEMENT, TIMESTEPPER >::InterfaceProblem().

Tau

Vector< double > oomph::Global_Physical_Variables::Tau	(	2	,
		1.	0
	)

Timescales for transport equations (identically one from our non-dimensionalisation)

Wall_normal

Vector< double > oomph::Global_Physical_Variables::Wall_normal

Direction of the wall normal vector (at the inlet)

Referenced by wall_unit_normal_inlet_fct(), and wall_unit_normal_outlet_fct().