two_d_adv_diff_SUPG.cc File Reference

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

Classes
class	SUPGAdvectionDiffusionProblem< ELEMENT >

Namespaces
	GlobalPhysicalParameters
	Namespace for Womersley problem.

Functions
void	GlobalPhysicalParameters::get_boundary_values (const Vector< double > &x, Vector< double > &u)
	Some "solution" for assignment of boundary values. More...
void	GlobalPhysicalParameters::source_function (const Vector< double > &x_vect, double &source)
	Zero source function. More...
void	GlobalPhysicalParameters::wind_function (const Vector< double > &x, Vector< double > &wind)
	Wind. More...
int	main ()
	Driver code for 2D AdvectionDiffusion problem. More...

Variables
double	GlobalPhysicalParameters::Peclet =200.0
	Peclet number. More...
double	GlobalPhysicalParameters::Alpha =50.0
	Parameter for steepness of step in boundary values. More...
double	GlobalPhysicalParameters::TanPhi =1.0
	Parameter for angle of step in boundary values: 45 degrees. More...

Function Documentation

main()

int main

(

)

Driver code for 2D AdvectionDiffusion problem.

{
 
 //Set up the problem with stabilisation
 {
  bool use_stabilisation=true;
  
  // Create the problem with 2D nine-node elements from the
  // QAdvectionDiffusionElement family. Pass pointer to 
  // source and wind function. 
  SUPGAdvectionDiffusionProblem<QSUPGAdvectionDiffusionElement<2,3> > 
   problem(&GlobalPhysicalParameters::source_function,
           &GlobalPhysicalParameters::wind_function,
           use_stabilisation);
  
  // Solve the problem
  problem.newton_solve();
  
  //Output the solution
  problem.doc_solution();
  
 }
 
 
 
 //Set up the problem without stabilisation
 {
  
  bool use_stabilisation=false;
  
  // Create the problem with 2D nine-node elements from the
  // QAdvectionDiffusionElement family. Pass pointer to 
  // source and wind function. 
  SUPGAdvectionDiffusionProblem<QSUPGAdvectionDiffusionElement<2,3> > 
   problem(&GlobalPhysicalParameters::source_function,
           &GlobalPhysicalParameters::wind_function,
           use_stabilisation);
  
  // Solve the problem
  problem.newton_solve();
  
  //Output the solution
  problem.doc_solution();
  
 }
 
 
 
} // end of main

References problem, GlobalPhysicalParameters::source_function(), and GlobalPhysicalParameters::wind_function().