axisym_heat_sphere.cc File Reference

#include "generic.h"
#include "half_rectangle_with_hole_mesh.h"
#include "axisym_buoyant_navier_stokes.h"
#include "assert.h"

Classes
class	DragNusseltCalculationElement< ELEMENT >
class	FlowAroundHalfCylinderProblem< ELEMENT >
	Flow around a cylinder in rectangular domain. More...

Namespaces
	Global_Parameters
	Namespace for global parameters.
	StokesFlowExactWind

Functions
void	StokesFlowExactWind::get_wind (const Vector< double > &x, Vector< double > &wind)
	Wind that represents a constantly translating sphere. More...
int	main ()
	Driver. More...

Variables
double	Global_Parameters::Pr = 0.73
	Prandtl number. More...
double	Global_Parameters::Pe = Re*Pr
	Peclet number. More...
double	Global_Parameters::Ra = 0.0
	Rayleigh number. More...
Vector< double >	Global_Parameters::G
	Gravity. More...
double	Global_Parameters::Sphere_centre_z = 0.0
	Location of the centre of the sphere on the axis. More...

Function Documentation

main()

int main

(

)

Driver.

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

{
 Global_Parameters::G.resize(3);
 Global_Parameters::G[0] = 0.0;
 Global_Parameters::G[1] = -1.0;
 Global_Parameters::G[2] = 0.0;
 
 // radius and Length of domain
 double radius=50.0;
 double length=100.0;
 Global_Parameters::Sphere_centre_z = 50.0;
 
 //Create a new ellipse object as the central cylinder
 HalfEllipse* cylinder_pt = 
  new HalfEllipse(Global_Parameters::Sphere_centre_z,0.5,0.5);
 
 // Create Problem
 FlowAroundHalfCylinderProblem 
 <RefineableBuoyantQAxisymCrouzeixRaviartElement> 
  problem(cylinder_pt,radius,length);
 
 // Solve adaptively with up to max_adapt rounds of refinement
 unsigned max_adapt=1;
 
 // Output filename
 char filename[100];
 // Trace file
 std::ofstream trace("trace.dat");
 
 //Step up in Reynolds number (and therefore Peclet number)
 for(unsigned i=0;i<2;i++)
  {
   problem.newton_solve(max_adapt);
 
   //Open an output file
   sprintf(filename,"soln_Re%g.dat", Global_Parameters::Re);
   //Doc result
   ofstream outfile(filename);
   problem.mesh_pt()->output(outfile,5);
   outfile.close();
 
   //Compute drag and mass transfer
   Vector<double> drag(2); double nusselt = 0.0; double area = 0.0;
   problem.compute_drag_nusselt(drag,nusselt, area);
   
   trace << Global_Parameters::Re << " " << Global_Parameters::Pr
         << " " << Global_Parameters::Pe << " " 
         << drag[0] << " " << drag[1] << " " << 2.0*nusselt << std::endl;
 
   Global_Parameters::Re += 10.0;
   Global_Parameters::Pe = Global_Parameters::Pr*Global_Parameters::Re;
  }
 
 trace.close();
}

References MergeRestartFiles::filename, Global_Parameters::G, i, Global_Parameters::Pe, Global_Parameters::Pr, problem, UniformPSDSelfTest::radius, Global_Parameters::Re, and Global_Parameters::Sphere_centre_z.