mpi/distribution/three_d_entry_flow/three_d_entry_flow.cc File Reference

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

Classes
class	EntryFlowProblem< ELEMENT >
	Entry flow problem in quarter tube domain. More...

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
int	main (int argc, char *argv[])

Function Documentation

main()

int main	(	int argc	,
		char *	argv[]
	)

///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// Driver for 3D entry flow into a quarter tube. If there are any command line arguments, we regard this as a validation run and perform only a single adaptation

{
#ifdef OOMPH_HAS_MPI
 MPI_Helpers::init(argc,argv);
#endif
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 DocInfo mesh_doc_info;
 bool report_stats=true;
 mesh_doc_info.set_directory("RESLT_MESH");
 mesh_doc_info.number()=10;
 
 // Allow (up to) five rounds of fully automatic adapation in response to 
 //-----------------------------------------------------------------------
 // error estimate
 //---------------
 unsigned max_adapt;
 double max_error_target,min_error_target;
 
 // Set max number of adaptations in black-box Newton solver and
 // error targets for adaptation
 if (CommandLineArgs::Argc==1)
  {
   // Up to two adaptations
   max_adapt=2;
 
   // Error targets for adaptive refinement
   max_error_target=0.005;
   min_error_target=0.0005;
  } 
 // Validation run: Only one adaptation. Relax error targets
 // to ensure that not all elements are refined so we're getting
 // some hanging nodes.
 else
  {
   // Validation run: Just one round of adaptation
   max_adapt=1;
   
   // Error targets for adaptive refinement
   max_error_target=0.02;
   min_error_target=0.002;
  }
 // end max_adapt setup 
 
 // Set up doc info
 DocInfo doc_info;
 
 // Do Taylor-Hood elements
 //------------------------
 {
  // Set output directory
  doc_info.set_directory("RESLT_TH");
  
  // Step number
  doc_info.number()=0;
  
  // Build problem
  EntryFlowProblem<RefineableQTaylorHoodElement<3> > 
   problem(doc_info,min_error_target,max_error_target);
  
  cout << " Doing Taylor-Hood elements " << std::endl;
 
  // Initial refinement
  problem.refine_uniformly();
 
  // Initial solve
  problem.newton_solve();
 
#ifdef OOMPH_HAS_MPI
  // Distribute
  std::ifstream input_file;
  std::ofstream output_file;
  char filename[100];
 
  // Get partition from file
  unsigned n_partition=problem.mesh_pt()->nelement();
  Vector<unsigned> element_partition(n_partition,0);
  sprintf(filename,"three_d_entry_flow_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 and check halo schemes
  problem.distribute(element_partition,mesh_doc_info,report_stats);
  problem.check_halo_schemes(mesh_doc_info);
 
  oomph_info << "---- Now solve TH after distribute ----" << std::endl;
#endif
 
  // Solve the problem with adaptivity
  problem.newton_solve(max_adapt);
 
  // Doc solution after solving
  problem.doc_solution();
 }
 
 
 // Do Crouzeix-Raviart elements
 //-----------------------------
 {
  // Set output directory
  doc_info.set_directory("RESLT_CR");
  
  // Step number
  doc_info.number()=0;
  
  // Build problem
  EntryFlowProblem<RefineableQCrouzeixRaviartElement<3> >
   problem(doc_info,min_error_target,max_error_target);
  
  cout << " Doing Crouzeix-Raviart elements " << std::endl;
  
  // Initial refinement
  problem.refine_uniformly();
 
  // Initial solve
  problem.newton_solve();
 
#ifdef OOMPH_HAS_MPI
  // Distribute
  std::ifstream input_file;
  std::ofstream output_file;
  char filename[100];
 
  // Get partition from file
  unsigned n_partition=problem.mesh_pt()->nelement();
  Vector<unsigned> element_partition(n_partition,0);
  sprintf(filename,"three_d_entry_flow_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 and check halo schemes
  problem.distribute(element_partition,mesh_doc_info,report_stats);
  problem.check_halo_schemes(mesh_doc_info);
 
  oomph_info << "---- Now solve CR after distribute ----" << std::endl;
#endif
 
  // Solve the problem with adaptivity
  problem.newton_solve(max_adapt);
 
  // Doc solution after solving
  problem.doc_solution();  
 }
 
#ifdef OOMPH_HAS_MPI
 MPI_Helpers::finalize();
#endif
 
} // end_of_main

References oomph::CommandLineArgs::Argc, e(), MergeRestartFiles::filename, oomph::DocInfo::number(), oomph::oomph_info, problem, oomph::DocInfo::set_directory(), Flag_definition::setup(), and oomph::Global_string_for_annotation::string().