Classes | Namespaces | Functions | Variables
gen_axisym_advection_diffusion_elements.h File Reference

Go to the source code of this file.

Classes

class  QGeneralisedAxisymAdvectionDiffusionElement< NNODE_1D >
 QGeneralisedAxisymAdvectionDiffusionElement elements are linear/quadrilateral/brick-shaped Advection Diffusion elements with isoparametric interpolation for the function. More...
 
class  FaceGeometry< QGeneralisedAxisymAdvectionDiffusionElement< NNODE_1D > >
 Face geometry for the QGeneralisedAxisymAdvectionDiffusionElement elements: The spatial dimension of the face elements is one lower than that of the bulk element but they have the same number of points along their 1D edges. More...
 

Namespaces

 oomph
 

Functions

virtual unsigned oomph::u_index_cons_axisym_adv_diff () const
 A class for all elements that solve the Advection Diffusion equations in conservative form using isoparametric elements in a cylindrical polar coordinate system. {$$} ( Pe {$w$}({$x$}) u. More...
 
double oomph::du_dt_cons_axisym_adv_diff (const unsigned &n) const
 du/dt at local node n. More...
 
void oomph::disable_ALE ()
 Disable ALE, i.e. assert the mesh is not moving – you do this at your own risk! More...
 
void oomph::enable_ALE ()
 (Re-)enable ALE, i.e. take possible mesh motion into account when evaluating the time-derivative. Note: By default, ALE is enabled, at the expense of possibly creating unnecessary work in problems where the mesh is, in fact, stationary. More...
 
void oomph::output (std::ostream &outfile)
 Output with default number of plot points. More...
 
void oomph::output (std::ostream &outfile, const unsigned &nplot)
 Output FE representation of soln: r,z,u at nplot^2 plot points. More...
 
void oomph::output (FILE *file_pt)
 C_style output with default number of plot points. More...
 
void oomph::output (FILE *file_pt, const unsigned &n_plot)
 C-style output FE representation of soln: r,z,u at n_plot^2 plot points. More...
 
void oomph::output_fct (std::ostream &outfile, const unsigned &nplot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 Output exact soln: r,z,u_exact at nplot^2 plot points. More...
 
virtual void oomph::output_fct (std::ostream &outfile, const unsigned &nplot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt)
 Output exact soln: r,z,,u_exact at nplot^2 plot points (dummy time-dependent version to keep intel compiler happy) More...
 
void oomph::compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 Get error against and norm of exact solution. More...
 
void oomph::compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 Dummy, time dependent error checker. More...
 
double oomph::integrate_u ()
 Integrate the concentration over the element. More...
 
GeneralisedAxisymAdvectionDiffusionSourceFctPt & oomph::source_fct_pt ()
 Access function: Pointer to source function. More...
 
GeneralisedAxisymAdvectionDiffusionWindFctPt & oomph::wind_fct_pt ()
 Access function: Pointer to wind function. More...
 
GeneralisedAxisymAdvectionDiffusionWindFctPt & oomph::conserved_wind_fct_pt ()
 Access function: Pointer to additional (conservative) wind function. More...
 
GeneralisedAxisymAdvectionDiffusionDiffFctPt & oomph::diff_fct_pt ()
 Access function: Pointer to diffusion function. More...
 
const double & oomph::pe () const
 Peclet number. More...
 
double *& oomph::pe_pt ()
 Pointer to Peclet number. More...
 
const double & oomph::pe_st () const
 Peclet number multiplied by Strouhal number. More...
 
double *& oomph::pe_st_pt ()
 Pointer to Peclet number multipled by Strouha number. More...
 
virtual void oomph::get_source_cons_axisym_adv_diff (const unsigned &ipt, const Vector< double > &x, double &source) const
 Get source term at (Eulerian) position x. This function is virtual to allow overloading in multi-physics problems where the strength of the source function might be determined by another system of equations. More...
 
virtual void oomph::get_wind_cons_axisym_adv_diff (const unsigned &ipt, const Vector< double > &s, const Vector< double > &x, Vector< double > &wind) const
 Get wind at (Eulerian) position x and/or local coordinate s. This function is virtual to allow overloading in multi-physics problems where the wind function might be determined by another system of equations. More...
 
virtual void oomph::get_conserved_wind_cons_axisym_adv_diff (const unsigned &ipt, const Vector< double > &s, const Vector< double > &x, Vector< double > &wind) const
 Get additional (conservative) wind at (Eulerian) position x and/or local coordinate s. This function is virtual to allow overloading in multi-physics problems where the wind function might be determined by another system of equations. More...
 
virtual void oomph::get_diff_cons_axisym_adv_diff (const unsigned &ipt, const Vector< double > &s, const Vector< double > &x, DenseMatrix< double > &D) const
 Get diffusivity tensor at (Eulerian) position x and/or local coordinate s. This function is virtual to allow overloading in multi-physics problems where the wind function might be determined by another system of equations. More...
 
void oomph::get_flux (const Vector< double > &s, Vector< double > &flux) const
 Get flux: $\mbox{flux}[i] = \mbox{d}u / \mbox{d}x_i $. More...
 
void oomph::get_total_flux (const Vector< double > &s, Vector< double > &total_flux) const
 Get flux: $\mbox{flux}[i] = \mbox{d}u / \mbox{d}x_i $. More...
 
void oomph::fill_in_contribution_to_residuals (Vector< double > &residuals)
 Add the element's contribution to its residual vector (wrapper) More...
 
void oomph::fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 Add the element's contribution to its residual vector and the element Jacobian matrix (wrapper) More...
 
void oomph::fill_in_contribution_to_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
double oomph::interpolated_u_cons_axisym_adv_diff (const Vector< double > &s) const
 Return FE representation of function value u(s) at local coordinate s. More...
 
unsigned oomph::self_test ()
 Self-test: Return 0 for OK. More...
 
virtual double oomph::dshape_and_dtest_eulerian_cons_axisym_adv_diff (const Vector< double > &s, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const =0
 Shape/test functions and derivs w.r.t. to global coords at local coord. s; return Jacobian of mapping. More...
 
virtual double oomph::dshape_and_dtest_eulerian_at_knot_cons_axisym_adv_diff (const unsigned &ipt, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const =0
 Shape/test functions and derivs w.r.t. to global coords at integration point ipt; return Jacobian of mapping. More...
 
virtual void oomph::fill_in_generic_residual_contribution_cons_axisym_adv_diff (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix, unsigned flag)
 Add the element's contribution to its residual vector only (if flag=and/or element Jacobian matrix. More...
 

Variables

double * oomph::Pe_pt
 Pointer to global Peclet number. More...
 
double * oomph::PeSt_pt
 Pointer to global Peclet number multiplied by Strouhal number. More...
 
GeneralisedAxisymAdvectionDiffusionSourceFctPt oomph::Source_fct_pt
 Pointer to source function: More...
 
GeneralisedAxisymAdvectionDiffusionWindFctPt oomph::Wind_fct_pt
 Pointer to wind function: More...
 
GeneralisedAxisymAdvectionDiffusionWindFctPt oomph::Conserved_wind_fct_pt
 Pointer to additional (conservative) wind function: More...
 
GeneralisedAxisymAdvectionDiffusionDiffFctPt oomph::Diff_fct_pt
 Pointer to diffusivity funciton. More...
 
bool oomph::ALE_is_disabled
 Boolean flag to indicate if ALE formulation is disabled when time-derivatives are computed. Only set to false if you're sure that the mesh is stationary. More...
 
static double oomph::Default_peclet_number
 Static default value for the Peclet number. More...