multi_domain_axisym_heat_sphere.cc File Reference

#include "generic.h"
#include "half_rectangle_with_hole_mesh.h"
#include "multi_domain_axisym_boussinesq_elements.h"

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

Namespaces
	Global_Parameters
	Namespace for global parameters.

Functions
int	main ()
	Driver. More...

Function Documentation

main()

int main

(

)

Driver.

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

{
 //Multi_domain_functions::Doc_stats=true;
 //Multi_domain_functions::Doc_full_stats=true;
 
 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 
  <RefineableQAxisymCrouzeixRaviartBoussinesqElement,
   RefineableQAxisymAdvectionDiffusionBoussinesqElement> 
  problem(cylinder_pt,radius,length);
 
 //Refine the problem a couple of times
 //problem.refine_uniformly();
 //problem.refine_uniformly();
 
 
 // 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,"vel_soln_Re%g.dat", Global_Parameters::Re);
   //Doc result
   ofstream outfile(filename);
   problem.nst_mesh_pt()->output(outfile,5);
   outfile.close();
   
   sprintf(filename,"conc_soln_Re%g.dat", Global_Parameters::Re);
   outfile.open(filename);
   problem.adv_diff_mesh_pt()->output(outfile,5);
   outfile.close();
 
   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.