multi_domain_b_convect_sph.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=5.0;
 double length=10.0;
 Global_Parameters::Sphere_centre_z = 5.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_Ra%g.dat", Global_Parameters::Re,
           Global_Parameters::Ra);
   //Doc result
   ofstream outfile(filename);
   problem.nst_mesh_pt()->output(outfile,5);
   outfile.close();
   
   sprintf(filename,"conc_soln_Re%g_Ra%g.dat", Global_Parameters::Re,
           Global_Parameters::Ra);
   outfile.open(filename);
   problem.adv_diff_mesh_pt()->output(outfile,5);
   outfile.close();
 
   Global_Parameters::Ra += 40.0;
  }
 
 //trace.close();
}

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