thermal_fibre.cc File Reference

#include "generic.h"
#include "axisym_navier_stokes.h"
#include "fluid_interface.h"
#include "steady_axisym_advection_diffusion.h"
#include "meshes/horizontal_single_layer_spine_mesh.h"
#include "non_isothermal_axisym_nst_elements.h"

Classes
class	AxisymFreeSurfaceNozzleAdvDiffRobinProblem< ELEMENT >

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
Vector< double >	Global_Physical_Variables::Direction_of_gravity (3)
	Gravity vector. More...
void	Global_Physical_Variables::viscosity_ratio_function (double &temperature, double &result)
	Viscosity ratio function modeled following a Arrhenius fashion. More...
void	Global_Physical_Variables::prescribed_beta_on_fixed_r_boundary (const Vector< double > &x_vector, double &beta)
	Beta on a boundary on which r is fixed. More...
void	Global_Physical_Variables::prescribed_alpha_on_fixed_r_boundary (const Vector< double > &x_vect, double &alpha)
	Alfa on a boundary on which r is fixed. More...
int	main (int argc, char **argv)
	Driver code for 2D Boussinesq convection problem. More...

Variables
double	Global_Physical_Variables::Bi = 0.01
	Biot number. More...
double	Global_Physical_Variables::T_ext = 0.0
	Exterior temperature. More...
double	Global_Physical_Variables::T_inlet = 1.0
	Inlet temperature. More...
double	Global_Physical_Variables::eta = 0.0
	Eta factor (exponent of the exponecial function) More...

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

Driver code for 2D Boussinesq convection problem.

//////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////////////

{
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 // # of elements in r-direction
 unsigned n_r = 8; //20;
 
 // # of elements in z-direction
 unsigned n_z = 22;
 
 // Domain length in r-direction
 double l_r = 1.0;
 
 // Domain length in z-direction
 double h = 11.0;
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 //FiniteElement::Accept_negative_jacobian = true;
 
 // Set direction of gravity (vertically downwards)
 Global_Physical_Variables::Direction_of_gravity[0] = 0.0;
 Global_Physical_Variables::Direction_of_gravity[1] =-1.0;
 Global_Physical_Variables::Direction_of_gravity[2] = 0.0;
 
 //Construct our problem
 AxisymFreeSurfaceNozzleAdvDiffRobinProblem<SpineElement<NonIsothermalAxisymmetricQCrouzeixRaviartElement> > 
 problem(n_r,n_z,l_r,h);
 
 // Run the steady simulation
 problem.steady_run();
 
} // end of main

References Global_Physical_Variables::Direction_of_gravity, problem, and Flag_definition::setup().