fsi_chan_precond_driver.cc File Reference

#include "generic.h"
#include "navier_stokes.h"
#include "beam.h"
#include "multi_physics.h"
#include "meshes/one_d_lagrangian_mesh.h"
#include "meshes/collapsible_channel_mesh.h"
#include "fsi_chan_problem.h"

Classes
class	PreconditionedFSICollapsibleChannelProblem< ELEMENT >
	Problem class. More...

Namespaces
	Flags
	Extend namespace for flags.

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

Variables
unsigned	Flags::Solver_flag =0
unsigned	Flags::Solver_sub_flag =0

Function Documentation

main()

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

////////////////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////////////// Driver code for a collapsible channel problem with FSI. Presence of command line arguments indicates validation run with coarse resolution and small number of steps.

Resolution factor

Use displacement control (1) or not (0)

Min. y coordinate for parameter study with displacement control

Steady (1) or unsteady (0) run

Number of steps

Reynolds number

Womersley number

FSI parameter

{ 
//#ifdef OOMPH_HAS_MPI
//  MPI_Helpers::init(argc,argv);
//#endif
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 if (CommandLineArgs::Argc==1)
  {
   std::cout << "Using default settings for flags" << std::endl;
  }
 else if ((CommandLineArgs::Argc==10)||(CommandLineArgs::Argc==12))
  {
   Flags::Resolution_factor=atoi(argv[1]);
 
   Flags::Use_displ_control=atoi(argv[2]);
 
   Global_Physical_Variables::Yprescr_min=double(atof(argv[3]));
 
   Flags::Steady_flag=atoi(argv[4]);
 
   Flags::Nsteps=atoi(argv[5]);
 
   // Solver flag [0: direct; 1: exact; 2: simple; 3: Schur/Superlu; 
   // 4: Schur/Hypre]
   Flags::Solver_flag=atoi(argv[6]);
   
   // Solver sub flag [0: diag; 1: retain fluid on solid; 
   // 2: retain solid on fluid]
   Flags::Solver_sub_flag=atoi(argv[7]);
 
   Global_Physical_Variables::Re=double(atof(argv[8]));
 
   Global_Physical_Variables::ReSt=Global_Physical_Variables::Re;
 
   Global_Physical_Variables::Q=double(atof(argv[9]));
 
   // Restart?
   if (CommandLineArgs::Argc==12)
    {
     // Name of restart file
     Flags::Restart_file_name=argv[10];
     
     // Jump in pressure 
     Global_Physical_Variables::P_step=double(atof(argv[11]));
    }
  }
 else
  {
   std::cout 
    << "\n\n\n\n\n"
    << "Wrong number of command line args: Specify none, six or eight:"
    << std::endl
    << "- resolution factor [1,2,...]" << std::endl
    << "- use_displ_control [0/1]" << std::endl
    << "- min. y-coordinate of control point when using displ control" 
    <<    std::endl
    << "- steady_flag [0/1]" << std::endl
    << "- number of steps " << std::endl
    << "- solver flag [0: direct; 1: exact; 2: simple; 3: Schur/Superlu; 4: Schur/Hypre]" << std::endl
    << "- solver sub flag [0: diag; 1: retain fluid on solid; 2: retain solid on fluid]" << std::endl
    << "- Reynolds number" << std::endl
    << "- FSI parameter Q" << std::endl
    << "- restart file name [optional] " << std::endl
    << "- jump in pressure P_step [optional] " << std::endl
    << "You specified " << CommandLineArgs::Argc-1 << " command line arg[s]" 
    << "\n\n\n\n\n" << std::endl;
   abort();
  }
 Flags::doc_flags();
 
 
 // Number of elements in the domain
 unsigned nup=4*Flags::Resolution_factor;
 unsigned ncollapsible=20*Flags::Resolution_factor;
 unsigned ndown=40*Flags::Resolution_factor;
 unsigned ny=4*Flags::Resolution_factor;
  
 
 // Length of the domain
 double lup=1.0;
 double lcollapsible=5.0;
 double ldown=10.0;
 double ly=1.0;
 
 // Use displacement control?
 bool displ_control=false;
 if (Flags::Use_displ_control==1) displ_control=true;
 
 // Steady run?
 bool steady_flag=false;
 if (Flags::Steady_flag==1) steady_flag=true;
 
 // Build the problem with QTaylorHoodElements
 PreconditionedFSICollapsibleChannelProblem
  <AlgebraicElement<QTaylorHoodElement<2> > > 
 problem(nup, ncollapsible, ndown, ny, 
         lup, lcollapsible, ldown, ly, displ_control,
         steady_flag,Flags::Solver_flag,Flags::Solver_sub_flag);
  
 if (Flags::Steady_flag)
  {
   problem.steady_run();
  }
 else
  {
   problem.unsteady_run();
  }
 
//#ifdef OOMPH_HAS_MPI
// MPI_Helpers::finalize();
//#endif
 
 
}//end of main

References oomph::CommandLineArgs::Argc, Flags::doc_flags(), Mesh_Parameters::ly, Flags::Nsteps, Mesh_Parameters::ny, Global_Physical_Variables::P_step, problem, Global_Physical_Variables::Q, Global_Physical_Variables::Re, Flags::Resolution_factor, Global_Physical_Variables::ReSt, Flags::Restart_file_name, Flag_definition::setup(), Flags::Solver_flag, Flags::Solver_sub_flag, Flags::Steady_flag, Flags::Use_displ_control, and Global_Physical_Variables::Yprescr_min.