fibre.cc File Reference

#include "generic.h"
#include "navier_stokes.h"
#include "axisym_navier_stokes.h"
#include "fluid_interface.h"
#include "meshes/horizontal_single_layer_spine_mesh.h"

Classes
class	MeltSpinningProblem< ELEMENT >
	Melt spinning axisymmetric fluid interface problem in rectangular domain. More...

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
int	main (int argc, char *argv[])
	Driver for isothermal fibre spinning axisymmetric fluid problem. More...

Variables
double	Global_Physical_Variables::Re = 0.1
	Reynolds number. More...
double	Global_Physical_Variables::ReInvFr = 0.01
	Product of Reynolds number and inverse of Froude number. More...
double	Global_Physical_Variables::Ca = 10.0
	Capillary number. More...
double	Global_Physical_Variables::P_ext = 0.0
	External pressure. More...
Vector< double >	Global_Physical_Variables::G (3) = 1.0
	Direction of gravity. More...

Function Documentation

main()

int main	(	int argc	,
		char *	argv[]
	)

Driver for isothermal fibre spinning axisymmetric fluid problem.

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

{
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 // Number of elements in radial (r) direction
 const unsigned n_r = 20;
   
 // Number of elements in axial (z) direction
 const unsigned n_z = 22;
 
 // Width of domain
 const double l_r = 1.0;
 
 // Height of fluid layer
 const double h = 11.0;
 
 // Accept negative jacobian 
 
 //FiniteElement::Accept_negative_jacobian = true;
 
 // Set direction of gravity (vertically downwards)
 Global_Physical_Variables::G[0] = 0.0;
 Global_Physical_Variables::G[1] =-1.0;
 Global_Physical_Variables::G[2] = 0.0;
 
 // Set up the spine test problem with AxisymmetricQCrouzeixRaviartElements,
 MeltSpinningProblem<SpineElement<AxisymmetricQCrouzeixRaviartElement > >
  problem(n_r,n_z,l_r,h);
 
 // Run the steady simulation
  problem.steady_run();
 
} // End of main

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