mpi/distribution/circular_driven_cavity/circular_driven_cavity.cc File Reference

#include "generic.h"
#include "navier_stokes.h"
#include "meshes/quarter_circle_sector_mesh.h"

Classes
class	QuarterCircleDrivenCavityProblem< ELEMENT >

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
Vector< double >	Global_Physical_Variables::Gravity (2)
	Gravity vector. More...
void	Global_Physical_Variables::body_force (const double &time, const Vector< double > &x, Vector< double > &result)
	Functional body force. More...
void	Global_Physical_Variables::zero_body_force (const double &time, const Vector< double > &x, Vector< double > &result)
	Zero functional body force. More...
int	main (int argc, char **argv)
	Driver for QuarterCircleDrivenCavityProblem test problem. More...

Variables
double	Global_Physical_Variables::Re_invFr =100
	Reynolds/Froude number. More...

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

Driver for QuarterCircleDrivenCavityProblem test problem.

{
#ifdef OOMPH_HAS_MPI
 MPI_Helpers::init(argc,argv);
#endif
 
 // Set output directory and initialise count
 DocInfo doc_info;
 doc_info.set_directory("RESLT");
 DocInfo mesh_doc_info;
 mesh_doc_info.set_directory("RESLT_MESH");
 bool report_stats=false;
 DocInfo error_doc_info;
 error_doc_info.set_directory("RESLT_ERROR");
 
 // Set max. number of black-box adaptation
 unsigned max_adapt=3; 
 
 // Solve problem 1 with Taylor-Hood elements
 //--------------------------------------------
 {
  oomph_info << " " << std::endl
             << "-----------------------" << std::endl
             << "-- Running problem 1 --" << std::endl
             << "-----------------------" << std::endl;
 
  // Set up downwards-Gravity vector
  Global_Physical_Variables::Gravity[0] = 0.0;
  Global_Physical_Variables::Gravity[1] = -1.0;
  
  // Set up Gamma vector for stress-divergence form
  NavierStokesEquations<2>::Gamma[0]=1;
  NavierStokesEquations<2>::Gamma[1]=1;
 
  // Build problem with Gravity vector in stress divergence form, 
  // using zero body force function
  QuarterCircleDrivenCavityProblem<RefineableQTaylorHoodElement<2> > 
   problem(&Global_Physical_Variables::zero_body_force);
 
  // output error estimates during solver
  error_doc_info.number()=0;
 
  RefineableMeshBase* mmesh_pt = 
   dynamic_cast<RefineableMeshBase*>(problem.mesh_pt(0));
 
  mmesh_pt->doc_info_pt()=&error_doc_info;
 
  // Initial solve
  problem.newton_solve();
 
#ifdef OOMPH_HAS_MPI
  mesh_doc_info.number()=0;
  std::ifstream input_file;
  std::ofstream output_file;
  char filename[100];
 
  // Get the partition to be used from file
  unsigned n_partition=problem.mesh_pt()->nelement();
  Vector<unsigned> element_partition(n_partition);
  sprintf(filename,"circular_cavity_1_partition.dat");
  input_file.open(filename);
  std::string input_string;
  for (unsigned e=0;e<n_partition;e++)
   {
    getline(input_file,input_string,'\n');
    element_partition[e]=atoi(input_string.c_str());
   }
 
  // Distribute the problem and check halo schemes
  problem.distribute(element_partition,mesh_doc_info,report_stats);
  problem.check_halo_schemes(mesh_doc_info);
 
  oomph_info << "---- Now solve after distribute ----" << std::endl;
#endif // (OOMPH_HAS_MPI)
 
  // Solve
  problem.newton_solve(max_adapt);
 
  // Step number
  doc_info.number()=1;
   
  //Output solution
  problem.doc_solution(doc_info);
 
 } // end of problem 1
 
 
 
 // Solve problem 2 with Crouzeix Raviart elements
 //--------------------------------------------
 {
  oomph_info << " " << std::endl
             << "-----------------------" << std::endl
             << "-- Running problem 2 --" << std::endl
             << "-----------------------" << std::endl;
 
  // Set up zero-Gravity vector
  Global_Physical_Variables::Gravity[0] = 0.0;
  Global_Physical_Variables::Gravity[1] = 0.0;
 
  // Set up Gamma vector for simplified form
  NavierStokesEquations<2>::Gamma[0]=0;
  NavierStokesEquations<2>::Gamma[1]=0;
 
  // Build problem with body force function and simplified form,
  // using body force function
  QuarterCircleDrivenCavityProblem<RefineableQCrouzeixRaviartElement<2> >
   problem(&Global_Physical_Variables::body_force);
 
  // Initial solve
  problem.newton_solve();
 
#ifdef OOMPH_HAS_MPI
  std::ifstream input_file;
  std::ofstream output_file;
  char filename[100];
 
  // Get the partition to be used from file
  unsigned n_partition=problem.mesh_pt()->nelement();
  Vector<unsigned> element_partition(n_partition);
  sprintf(filename,"circular_cavity_2_partition.dat");
  input_file.open(filename);
  std::string input_string;
  for (unsigned e=0;e<n_partition;e++)
   {
    getline(input_file,input_string,'\n');
    element_partition[e]=atoi(input_string.c_str());
   }
 
  // Distribute problem and check halo schemes
  problem.distribute(element_partition,mesh_doc_info,report_stats);
  problem.check_halo_schemes(mesh_doc_info);
 
  oomph_info << "---- Now solve after distribute ----" << std::endl;
#endif // (OOMPH_HAS_MPI)
 
  // Solve the problem with automatic adaptation
  problem.newton_solve(max_adapt);
  
  // Step number
  doc_info.number()=2;
 
  // Output solution
  problem.doc_solution(doc_info);
 
 } // end of problem 2
 
#ifdef OOMPH_HAS_MPI
// Finalise MPI
 MPI_Helpers::finalize();
#endif
 
} // end_of_main

References Global_Physical_Variables::body_force(), oomph::RefineableMeshBase::doc_info_pt(), e(), MergeRestartFiles::filename, Global_Physical_Variables::Gravity, oomph::DocInfo::number(), oomph::oomph_info, problem, oomph::DocInfo::set_directory(), oomph::Global_string_for_annotation::string(), and Global_Physical_Variables::zero_body_force().