channel_with_leaflet.cc File Reference

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

Classes
class	Leaflet
class	ChannelWithLeafletProblem< ELEMENT >
	Problem class. More...

Namespaces
	Global_Physical_Variables
	Global variables.

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

Function Documentation

main()

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

//////////////////////////////////////////////////////////////////// /////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////// Driver code – pass a command line argument if you want to run the code in validation mode where it only performs a few steps

Initialise timestep

Set max. number of adaptations (reduced for validation)

Output steady solution

Reduce the max number of adaptations for time-dependent simulation

{
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 // Set up doc info
 DocInfo doc_info; 
 doc_info.set_directory("RESLT");
 doc_info.number()=0;
 
 // Parameters for the leaflet
 //----------------------------
 
 // Height
 double h_leaflet = 0.5;
 
 
 // Tip deflection
 double d_x = 0.25;
 double d_y = -0.05;
 
 // x-positon of root
 double x_0 = 3.0; 
 
 // Period of the oscillation on the natural timescale of the flow
 double period = 20.0;
 
 
 //Parameters for the domain
 //-------------------------
 
 // Length of the mesh to right and left of the leaflet
 double l_left =2.0;
 double l_right= 3.0;
 
 // Total height of domain (unity because lengths have been scaled on it)
 double h_tot=1.0;
 
 // Initial number of element rows/columns in various mesh regions
 unsigned nleft=8;
 unsigned nright=12;
 unsigned ny1=2;
 unsigned ny2=2; 
 
 //Build the problem
 ChannelWithLeafletProblem<AlgebraicElement<RefineableQTaylorHoodElement<2> > >
  problem(l_left, l_right, h_leaflet,
          h_tot,nleft, nright,ny1,ny2,
          d_x, d_y, x_0,
          period);
 
 
 // Number of timesteps per period
 unsigned nsteps_per_period=40;
 
 // Number of periods
 unsigned nperiod=3; 
 
 // Number of timesteps (reduced for validation)
 unsigned nstep=nsteps_per_period*nperiod;
 if (CommandLineArgs::Argc>1)
  {
   nstep=3;
  }
 
 //Timestep: 
 double dt=period/double(nsteps_per_period);
 
 problem.initialise_dt(dt);
 
 
 unsigned max_adapt=5;
 if (CommandLineArgs::Argc>1)
  {
   max_adapt=2;
  }
 
 // Do steady solve first -- this also sets the history values
 // to those corresponding to an impulsive start from the
 // steady solution
 problem.steady_newton_solve(max_adapt);
 
 problem.doc_solution(doc_info);
 doc_info.number()++;
 
 
 max_adapt=1;
 
 // We don't want to re-assign the initial condition 
 bool first=false;
 
// Timestepping loop
 for (unsigned istep=0;istep<nstep;istep++)
  { 
   // Solve the problem
   problem.unsteady_newton_solve(dt,max_adapt,first);
   
   // Output the solution
   problem.doc_solution(doc_info);
   
   // Step number
   doc_info.number()++;
 
  }
 
 
}//end of main

References oomph::CommandLineArgs::Argc, oomph::DocInfo::number(), problem, oomph::DocInfo::set_directory(), and Flag_definition::setup().