spherical_couette.cc File Reference

#include "generic.h"
#include "spherical_navier_stokes.h"
#include "meshes/rectangular_quadmesh.h"

Classes
class	TorqueCalculationElement< ELEMENT >
class	RefineableSphericalCouetteProblem< ELEMENT >
	Driven cavity problem in rectangular domain - time dependent version. More...

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
int	main ()
	Driver for RefineableSphericalCouette test problem using a. More...

Function Documentation

main()

int main

(

)

Driver for RefineableSphericalCouette test problem using a.

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

{
 // Label for output
 DocInfo doc_info;
 
 // Set output directory
 doc_info.set_directory("RESLT");
 
 // ---------------
 // end of Set up doc info
 
 
 // Doing QCrouzeixRaviartElements
 {
  // Build the problem with QCrouzeixRaviartElements
  RefineableSphericalCouetteProblem<
   RefineableQSphericalCrouzeixRaviartElement > 
   problem;
  cout << "Doing QCrouzeixRaviartElement" << std::endl;
  
  // Open a trace file
  ofstream trace_file;
  char filename[100];   
  sprintf(filename,"%s/trace.dat",doc_info.directory().c_str());
  trace_file.open(filename);
  
  // Over-ride the maximum and minimum permitted errors
  problem.mesh_pt()->max_permitted_error() = 1.0e-3; //Default = 1.0e-3
  problem.mesh_pt()->min_permitted_error() = 1.0e-5; //Default = 1.0e-5
  
  // Over-ride the maximum and minimum permitted refinement levels
  problem.mesh_pt()->max_refinement_level() = 4;//maximum_ref_level;
  problem.mesh_pt()->min_refinement_level() = 1;//minimum_ref_level;
  
  //Set the boundary conditions
  problem.set_boundary_conditions();
  
  //Set the maximum adaptation
  unsigned max_adapt = 10;
 
  for(unsigned i=0;i<5;i++)
   {
    // Solve the problem
    problem.steady_newton_solve(max_adapt);
 
    //Scale the torque to be consistent with the non-dimensionalisation
    //in Dennis, Singh & Ingham (1980) and many others
    double torque = 3.0*problem.compute_torque()/4.0;
    std::cout << "Torque is " 
              << torque
              << "\n";
    
    problem.doc_solution(doc_info,trace_file);
    
    trace_file << Global_Physical_Variables::Re << " " 
               << " " << torque
               << std::endl;
 
      Global_Physical_Variables::Re += 25.0;
   }
  
  // Close trace file
  trace_file.close();
  
 } // end of QCrouzeixRaviartElements
 
 
} // end_of_main

References oomph::DocInfo::directory(), MergeRestartFiles::filename, i, problem, Global_Physical_Variables::Re, and oomph::DocInfo::set_directory().