pipe_variable_diff.cc File Reference

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

Classes
class	RefineableAdvectionDiffusionPipeProblem< ELEMENT >

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
void	Global_Physical_Variables::wind_function (const Vector< double > &x, Vector< double > &wind)
	Wind. More...
void	Global_Physical_Variables::swimming (const Vector< double > &x, Vector< double > &swim)
void	Global_Physical_Variables::diff_function (const Vector< double > &x, DenseMatrix< double > &D)
int	main ()
	Driver code for 2D AdvectionDiffusion problem. More...

Variables
double	Global_Physical_Variables::Peclet_St = 1.0
	Peclet number multiplied by Strouhal number. More...

Function Documentation

main()

int main

(

)

Driver code for 2D AdvectionDiffusion problem.

{
 
 //Set up the problem
 //------------------
 
 // Create the problem with axisymmetric nine-node refineable elements from the
 // RefineableQAdvectionDiffusionElement family.
 RefineableAdvectionDiffusionPipeProblem<
 RefineableQGeneralisedAxisymAdvectionDiffusionElement<3> > problem;
 
 //Refine uniformly once
 problem.refine_uniformly();
 //Set the initial distribution
 problem.set_initial_condition();
 problem.doc_solution();
 
 //Set the timestep
 double dt = 0.001;
 problem.assign_initial_values_impulsive(dt);
 
 //Now take a few timesteps
 bool first = true;
 unsigned max_adapt=2;
 for(unsigned t=0;t<3;t++)
  {
   problem.unsteady_newton_solve(dt,max_adapt,first);
   problem.doc_solution();
   if(t==0) {first=false;}
  }
 
} // end of main

References problem, and plotPSD::t.