LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement Class Reference

#include <multi_domain_linearised_axisym_navier_stokes_elements.h>

+ Inheritance diagram for LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement:

Public Member Functions

 LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement ()
 Constructor: call the underlying constructors. More...
 
virtual void get_base_flow_u (const double &time, const unsigned &ipt, const Vector< double > &x, Vector< double > &result) const
 
virtual void get_base_flow_dudx (const double &time, const unsigned &ipt, const Vector< double > &x, DenseMatrix< double > &result) const
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 Compute the element's residual vector and the Jacobian matrix. More...
 
 LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement ()
 Constructor: call the underlying constructors. More...
 
virtual void get_base_flow_u (const double &time, const unsigned &ipt, const Vector< double > &x, Vector< double > &result) const
 
virtual void get_base_flow_dudx (const double &time, const unsigned &ipt, const Vector< double > &x, DenseMatrix< double > &result) const
 
virtual void get_base_flow_p (const double &time, const unsigned &ipt, const Vector< double > &x, double &result) const
 
virtual void get_base_flow_dudt (const double &time, const unsigned &ipt, const Vector< double > &x, Vector< double > &result) const
 
virtual void get_base_flow_duds (const double &time, const unsigned &ipt, const Vector< double > &x, DenseMatrix< double > &result) const
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 Compute the element's residual vector and the Jacobian matrix. More...
 
 LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement ()
 Constructor: call the underlying constructors. More...
 
virtual void get_base_flow_u (const double &time, const unsigned &ipt, const Vector< double > &x, Vector< double > &result) const
 
virtual void get_base_flow_dudx (const double &time, const unsigned &ipt, const Vector< double > &x, DenseMatrix< double > &result) const
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 Compute the element's residual vector and the Jacobian matrix. More...
 
- Public Member Functions inherited from oomph::LinearisedAxisymmetricQCrouzeixRaviartElement
 LinearisedAxisymmetricQCrouzeixRaviartElement ()
 
virtual unsigned required_nvalue (const unsigned &n) const
 Return number of values (pinned or dofs) required at local node n. More...
 
double p_linearised_axi_nst (const unsigned &i_internal, const unsigned &i) const
 
unsigned npres_linearised_axi_nst () const
 
void fix_pressure (const unsigned &p_dof, const double &pvalue)
 
int p_local_eqn (const unsigned &n, const unsigned &i)
 
void output (std::ostream &outfile)
 Redirect output to NavierStokesEquations output. More...
 
void output (std::ostream &outfile, const unsigned &n_plot)
 Redirect output to NavierStokesEquations output. More...
 
void output (FILE *file_pt)
 Redirect output to NavierStokesEquations output. More...
 
void output (FILE *file_pt, const unsigned &n_plot)
 Redirect output to NavierStokesEquations output. More...
 
 LinearisedAxisymmetricQCrouzeixRaviartElement ()
 
virtual unsigned required_nvalue (const unsigned &n) const
 Return number of values (pinned or dofs) required at local node n. More...
 
double p_lin_axi_nst (const unsigned &i_internal, const unsigned &i) const
 
unsigned npres_lin_axi_nst () const
 
void fix_pressure (const unsigned &p_dof, const double &pvalue)
 
void fix_cosine_component_of_pressure (const unsigned &p_dof, const double &pvalue)
 
void fix_sine_component_of_pressure (const unsigned &p_dof, const double &pvalue)
 
int p_local_eqn (const unsigned &n, const unsigned &i)
 
void output (std::ostream &outfile)
 Redirect output to NavierStokesEquations output. More...
 
void output (std::ostream &outfile, const unsigned &n_plot)
 Redirect output to NavierStokesEquations output. More...
 
void output (FILE *file_pt)
 Redirect output to NavierStokesEquations output. More...
 
void output (FILE *file_pt, const unsigned &n_plot)
 Redirect output to NavierStokesEquations output. More...
 
unsigned ndof_types () const
 
 LinearisedAxisymmetricQCrouzeixRaviartElement ()
 
virtual unsigned required_nvalue (const unsigned &n) const
 Return number of values (pinned or dofs) required at local node n. More...
 
double p_linearised_axi_nst (const unsigned &i_internal, const unsigned &i) const
 
unsigned npres_linearised_axi_nst () const
 
void fix_pressure (const unsigned &p_dof, const double &pvalue)
 
int p_local_eqn (const unsigned &n, const unsigned &i)
 
void output (std::ostream &outfile)
 Redirect output to NavierStokesEquations output. More...
 
void output (std::ostream &outfile, const unsigned &n_plot)
 Redirect output to NavierStokesEquations output. More...
 
void output (FILE *file_pt)
 Redirect output to NavierStokesEquations output. More...
 
void output (FILE *file_pt, const unsigned &n_plot)
 Redirect output to NavierStokesEquations output. More...
 
unsigned ndof_types () const
 
- Public Member Functions inherited from oomph::LinearisedAxisymmetricNavierStokesEquations
 LinearisedAxisymmetricNavierStokesEquations ()
 
const doublere () const
 Reynolds number. More...
 
const doublere_st () const
 Product of Reynolds and Strouhal number (=Womersley number) More...
 
const intazimuthal_mode_number () const
 Azimuthal mode number k in e^ik(theta) decomposition. More...
 
double *& re_pt ()
 Pointer to Reynolds number. More...
 
double *& re_st_pt ()
 Pointer to product of Reynolds and Strouhal number (=Womersley number) More...
 
int *& azimuthal_mode_number_pt ()
 Pointer to azimuthal mode number k in e^ik(theta) decomposition. More...
 
const doubleviscosity_ratio () const
 
double *& viscosity_ratio_pt ()
 Pointer to the viscosity ratio. More...
 
const doubledensity_ratio () const
 
double *& density_ratio_pt ()
 Pointer to the density ratio. More...
 
virtual unsigned u_index_linearised_axi_nst (const unsigned &i) const
 
double du_dt_linearised_axi_nst (const unsigned &n, const unsigned &i) const
 
void disable_ALE ()
 
void enable_ALE ()
 
virtual int p_index_linearised_axi_nst (const unsigned &i) const
 Which nodal value represents the pressure? More...
 
void strain_rate (const Vector< double > &s, DenseMatrix< double > &strain_rate) const
 
void output_veloc (std::ostream &outfile, const unsigned &nplot, const unsigned &t)
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Compute the element's residual Vector. More...
 
void fill_in_contribution_to_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
double interpolated_u_linearised_axi_nst (const Vector< double > &s, const unsigned &i) const
 
double interpolated_p_linearised_axi_nst (const Vector< double > &s, const unsigned &i) const
 
 LinearisedAxisymmetricNavierStokesEquations ()
 
const doublere () const
 Reynolds number. More...
 
double *& re_pt ()
 Pointer to Reynolds number. More...
 
const doublere_st () const
 Product of Reynolds and Strouhal number (=Womersley number) More...
 
double *& re_st_pt ()
 Pointer to product of Reynolds and Strouhal number (=Womersley number) More...
 
const doublere_invfr () const
 Global inverse Froude number. More...
 
double *& re_invfr_pt ()
 Pointer to global inverse Froude number. More...
 
const doublere_invro () const
 Global Reynolds number multiplied by inverse Rossby number. More...
 
double *& re_invro_pt ()
 Pointer to global Reynolds number multiplied by inverse Rossby number. More...
 
const Vector< double > & g () const
 Vector of gravitational components. More...
 
Vector< double > *& g_pt ()
 Pointer to Vector of gravitational components. More...
 
const intazimuthal_mode_number () const
 Azimuthal mode number k in e^ik(theta) decomposition. More...
 
int *& azimuthal_mode_number_pt ()
 Pointer to azimuthal mode number k in e^ik(theta) decomposition. More...
 
const doubleviscosity_ratio () const
 
double *& viscosity_ratio_pt ()
 Pointer to the viscosity ratio. More...
 
const doubledensity_ratio () const
 
double *& density_ratio_pt ()
 Pointer to the density ratio. More...
 
virtual int get_local_eqn_number_corresponding_to_geometric_dofs (const unsigned &n, const unsigned &i)=0
 
virtual unsigned xhat_index_lin_axi_nst (const unsigned &i) const
 
virtual unsigned u_index_lin_axi_nst (const unsigned &i) const
 
double du_dt_lin_axi_nst (const unsigned &n, const unsigned &i) const
 
double dnodal_position_perturbation_dt_lin_axi_nst (const unsigned &n, const unsigned &i) const
 
unsigned nexternal_H_data ()
 Return the number of external perturbed spine "heights" data. More...
 
void disable_ALE ()
 
void enable_ALE ()
 
virtual int p_index_lin_axi_nst (const unsigned &i) const
 Which nodal value represents the pressure? More...
 
void dkin_energy_dt (double &dkin_en_dt, double &kin_en) const
 Get integral of kinetic energy over element plus deriv w.r.t. time. More...
 
void strain_rate (const Vector< double > &s, DenseMatrix< double > &strain_rate) const
 
void output_veloc (std::ostream &outfile, const unsigned &nplot, const unsigned &t)
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Compute the element's residual Vector. More...
 
void fill_in_contribution_to_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
double interpolated_u_lin_axi_nst (const Vector< double > &s, const unsigned &i) const
 
double interpolated_p_lin_axi_nst (const Vector< double > &s, const unsigned &i) const
 
double interpolated_nodal_position_perturbation_lin_axi_nst (const Vector< double > &s, const unsigned &i) const
 
void assign_additional_local_eqn_numbers ()
 Define the local equation numbering schemes. More...
 
 LinearisedAxisymmetricNavierStokesEquations ()
 
const doublere () const
 Reynolds number. More...
 
const doublere_st () const
 Product of Reynolds and Strouhal number (=Womersley number) More...
 
const intazimuthal_mode_number () const
 Azimuthal mode number k in e^ik(theta) decomposition. More...
 
double *& re_pt ()
 Pointer to Reynolds number. More...
 
double *& re_st_pt ()
 Pointer to product of Reynolds and Strouhal number (=Womersley number) More...
 
int *& azimuthal_mode_number_pt ()
 Pointer to azimuthal mode number k in e^ik(theta) decomposition. More...
 
const doubleviscosity_ratio () const
 
double *& viscosity_ratio_pt ()
 Pointer to the viscosity ratio. More...
 
const doubledensity_ratio () const
 
double *& density_ratio_pt ()
 Pointer to the density ratio. More...
 
virtual unsigned u_index_linearised_axi_nst (const unsigned &i) const
 
double du_dt_linearised_axi_nst (const unsigned &n, const unsigned &i) const
 
void disable_ALE ()
 
void enable_ALE ()
 
virtual int p_index_linearised_axi_nst (const unsigned &i) const
 Which nodal value represents the pressure? More...
 
void strain_rate (const Vector< double > &s, DenseMatrix< double > &strain_rate) const
 
void output_veloc (std::ostream &outfile, const unsigned &nplot, const unsigned &t)
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Compute the element's residual Vector. More...
 
void fill_in_contribution_to_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
double interpolated_u_linearised_axi_nst (const Vector< double > &s, const unsigned &i) const
 
double interpolated_p_linearised_axi_nst (const Vector< double > &s, const unsigned &i) const
 
- Public Member Functions inherited from oomph::FiniteElement
void set_dimension (const unsigned &dim)
 
void set_nodal_dimension (const unsigned &nodal_dim)
 
void set_nnodal_position_type (const unsigned &nposition_type)
 Set the number of types required to interpolate the coordinate. More...
 
void set_n_node (const unsigned &n)
 
int nodal_local_eqn (const unsigned &n, const unsigned &i) const
 
double dJ_eulerian_at_knot (const unsigned &ipt, Shape &psi, DenseMatrix< double > &djacobian_dX) const
 
 FiniteElement ()
 Constructor. More...
 
virtual ~FiniteElement ()
 
 FiniteElement (const FiniteElement &)=delete
 Broken copy constructor. More...
 
virtual bool local_coord_is_valid (const Vector< double > &s)
 Broken assignment operator. More...
 
virtual void move_local_coord_back_into_element (Vector< double > &s) const
 
void get_centre_of_gravity_and_max_radius_in_terms_of_zeta (Vector< double > &cog, double &max_radius) const
 
virtual void local_coordinate_of_node (const unsigned &j, Vector< double > &s) const
 
virtual void local_fraction_of_node (const unsigned &j, Vector< double > &s_fraction)
 
virtual double local_one_d_fraction_of_node (const unsigned &n1d, const unsigned &i)
 
virtual void set_macro_elem_pt (MacroElement *macro_elem_pt)
 
MacroElementmacro_elem_pt ()
 Access function to pointer to macro element. More...
 
void get_x (const Vector< double > &s, Vector< double > &x) const
 
void get_x (const unsigned &t, const Vector< double > &s, Vector< double > &x)
 
virtual void get_x_from_macro_element (const Vector< double > &s, Vector< double > &x) const
 
virtual void get_x_from_macro_element (const unsigned &t, const Vector< double > &s, Vector< double > &x)
 
virtual void set_integration_scheme (Integral *const &integral_pt)
 Set the spatial integration scheme. More...
 
Integral *const & integral_pt () const
 Return the pointer to the integration scheme (const version) More...
 
virtual void shape (const Vector< double > &s, Shape &psi) const =0
 
virtual void shape_at_knot (const unsigned &ipt, Shape &psi) const
 
virtual void dshape_local (const Vector< double > &s, Shape &psi, DShape &dpsids) const
 
virtual void dshape_local_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsids) const
 
virtual void d2shape_local (const Vector< double > &s, Shape &psi, DShape &dpsids, DShape &d2psids) const
 
virtual void d2shape_local_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsids, DShape &d2psids) const
 
virtual double J_eulerian (const Vector< double > &s) const
 
virtual double J_eulerian_at_knot (const unsigned &ipt) const
 
void check_J_eulerian_at_knots (bool &passed) const
 
void check_jacobian (const double &jacobian) const
 
double dshape_eulerian (const Vector< double > &s, Shape &psi, DShape &dpsidx) const
 
virtual double dshape_eulerian_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsidx) const
 
virtual double dshape_eulerian_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsi, DenseMatrix< double > &djacobian_dX, RankFourTensor< double > &d_dpsidx_dX) const
 
double d2shape_eulerian (const Vector< double > &s, Shape &psi, DShape &dpsidx, DShape &d2psidx) const
 
virtual double d2shape_eulerian_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsidx, DShape &d2psidx) const
 
virtual void assign_nodal_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual void describe_nodal_local_dofs (std::ostream &out, const std::string &current_string) const
 
virtual void assign_all_generic_local_eqn_numbers (const bool &store_local_dof_pt)
 
Node *& node_pt (const unsigned &n)
 Return a pointer to the local node n. More...
 
Node *const & node_pt (const unsigned &n) const
 Return a pointer to the local node n (const version) More...
 
unsigned nnode () const
 Return the number of nodes. More...
 
virtual unsigned nnode_1d () const
 
double raw_nodal_position (const unsigned &n, const unsigned &i) const
 
double raw_nodal_position (const unsigned &t, const unsigned &n, const unsigned &i) const
 
double raw_dnodal_position_dt (const unsigned &n, const unsigned &i) const
 
double raw_dnodal_position_dt (const unsigned &n, const unsigned &j, const unsigned &i) const
 
double raw_nodal_position_gen (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double raw_nodal_position_gen (const unsigned &t, const unsigned &n, const unsigned &k, const unsigned &i) const
 
double raw_dnodal_position_gen_dt (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double raw_dnodal_position_gen_dt (const unsigned &j, const unsigned &n, const unsigned &k, const unsigned &i) const
 
double nodal_position (const unsigned &n, const unsigned &i) const
 
double nodal_position (const unsigned &t, const unsigned &n, const unsigned &i) const
 
double dnodal_position_dt (const unsigned &n, const unsigned &i) const
 Return the i-th component of nodal velocity: dx/dt at local node n. More...
 
double dnodal_position_dt (const unsigned &n, const unsigned &j, const unsigned &i) const
 
double nodal_position_gen (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double nodal_position_gen (const unsigned &t, const unsigned &n, const unsigned &k, const unsigned &i) const
 
double dnodal_position_gen_dt (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double dnodal_position_gen_dt (const unsigned &j, const unsigned &n, const unsigned &k, const unsigned &i) const
 
virtual void get_dresidual_dnodal_coordinates (RankThreeTensor< double > &dresidual_dnodal_coordinates)
 
unsigned nnodal_position_type () const
 
bool has_hanging_nodes () const
 
unsigned nodal_dimension () const
 Return the required Eulerian dimension of the nodes in this element. More...
 
virtual unsigned nvertex_node () const
 
virtual Nodevertex_node_pt (const unsigned &j) const
 
virtual Nodeconstruct_node (const unsigned &n)
 
virtual Nodeconstruct_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 
virtual Nodeconstruct_boundary_node (const unsigned &n)
 
virtual Nodeconstruct_boundary_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 
int get_node_number (Node *const &node_pt) const
 
virtual Nodeget_node_at_local_coordinate (const Vector< double > &s) const
 
double raw_nodal_value (const unsigned &n, const unsigned &i) const
 
double raw_nodal_value (const unsigned &t, const unsigned &n, const unsigned &i) const
 
double nodal_value (const unsigned &n, const unsigned &i) const
 
double nodal_value (const unsigned &t, const unsigned &n, const unsigned &i) const
 
unsigned dim () const
 
virtual ElementGeometry::ElementGeometry element_geometry () const
 Return the geometry type of the element (either Q or T usually). More...
 
virtual double interpolated_x (const Vector< double > &s, const unsigned &i) const
 Return FE interpolated coordinate x[i] at local coordinate s. More...
 
virtual double interpolated_x (const unsigned &t, const Vector< double > &s, const unsigned &i) const
 
virtual void interpolated_x (const Vector< double > &s, Vector< double > &x) const
 Return FE interpolated position x[] at local coordinate s as Vector. More...
 
virtual void interpolated_x (const unsigned &t, const Vector< double > &s, Vector< double > &x) const
 
virtual double interpolated_dxdt (const Vector< double > &s, const unsigned &i, const unsigned &t)
 
virtual void interpolated_dxdt (const Vector< double > &s, const unsigned &t, Vector< double > &dxdt)
 
unsigned ngeom_data () const
 
Datageom_data_pt (const unsigned &j)
 
void position (const Vector< double > &zeta, Vector< double > &r) const
 
void position (const unsigned &t, const Vector< double > &zeta, Vector< double > &r) const
 
void dposition_dt (const Vector< double > &zeta, const unsigned &t, Vector< double > &drdt)
 
virtual double zeta_nodal (const unsigned &n, const unsigned &k, const unsigned &i) const
 
void interpolated_zeta (const Vector< double > &s, Vector< double > &zeta) const
 
void locate_zeta (const Vector< double > &zeta, GeomObject *&geom_object_pt, Vector< double > &s, const bool &use_coordinate_as_initial_guess=false)
 
virtual void node_update ()
 
virtual void identify_field_data_for_interactions (std::set< std::pair< Data *, unsigned >> &paired_field_data)
 
virtual void identify_geometric_data (std::set< Data * > &geometric_data_pt)
 
virtual double s_min () const
 Min value of local coordinate. More...
 
virtual double s_max () const
 Max. value of local coordinate. More...
 
double size () const
 
virtual double compute_physical_size () const
 
virtual void point_output_data (const Vector< double > &s, Vector< double > &data)
 
void point_output (std::ostream &outfile, const Vector< double > &s)
 
virtual unsigned nplot_points_paraview (const unsigned &nplot) const
 
virtual unsigned nsub_elements_paraview (const unsigned &nplot) const
 
void output_paraview (std::ofstream &file_out, const unsigned &nplot) const
 
virtual void write_paraview_output_offset_information (std::ofstream &file_out, const unsigned &nplot, unsigned &counter) const
 
virtual void write_paraview_type (std::ofstream &file_out, const unsigned &nplot) const
 
virtual void write_paraview_offsets (std::ofstream &file_out, const unsigned &nplot, unsigned &offset_sum) const
 
virtual unsigned nscalar_paraview () const
 
virtual void scalar_value_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot) const
 
virtual void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt) const
 
virtual void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt) const
 
virtual std::string scalar_name_paraview (const unsigned &i) const
 
virtual void output (const unsigned &t, std::ostream &outfile, const unsigned &n_plot) const
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 Output an exact solution over the element. More...
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt)
 Output a time-dependent exact solution over the element. More...
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, const SolutionFunctorBase &exact_soln) const
 Output a time-dependent exact solution over the element. More...
 
virtual void get_s_plot (const unsigned &i, const unsigned &nplot, Vector< double > &s, const bool &shifted_to_interior=false) const
 
virtual std::string tecplot_zone_string (const unsigned &nplot) const
 
virtual void write_tecplot_zone_footer (std::ostream &outfile, const unsigned &nplot) const
 
virtual void write_tecplot_zone_footer (FILE *file_pt, const unsigned &nplot) const
 
virtual unsigned nplot_points (const unsigned &nplot) const
 
virtual void compute_error (FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 Calculate the norm of the error and that of the exact solution. More...
 
virtual void compute_error (FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 Calculate the norm of the error and that of the exact solution. More...
 
virtual void compute_error (FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_abs_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error)
 
void integrate_fct (FiniteElement::SteadyExactSolutionFctPt integrand_fct_pt, Vector< double > &integral)
 Integrate Vector-valued function over element. More...
 
void integrate_fct (FiniteElement::UnsteadyExactSolutionFctPt integrand_fct_pt, const double &time, Vector< double > &integral)
 Integrate Vector-valued time-dep function over element. More...
 
virtual void build_face_element (const int &face_index, FaceElement *face_element_pt)
 
virtual unsigned self_test ()
 
virtual unsigned get_bulk_node_number (const int &face_index, const unsigned &i) const
 
virtual int face_outer_unit_normal_sign (const int &face_index) const
 Get the sign of the outer unit normal on the face given by face_index. More...
 
virtual unsigned nnode_on_face () const
 
void face_node_number_error_check (const unsigned &i) const
 Range check for face node numbers. More...
 
virtual CoordinateMappingFctPt face_to_bulk_coordinate_fct_pt (const int &face_index) const
 
virtual BulkCoordinateDerivativesFctPt bulk_coordinate_derivatives_fct_pt (const int &face_index) const
 
- Public Member Functions inherited from oomph::GeneralisedElement
 GeneralisedElement ()
 Constructor: Initialise all pointers and all values to zero. More...
 
virtual ~GeneralisedElement ()
 Virtual destructor to clean up any memory allocated by the object. More...
 
 GeneralisedElement (const GeneralisedElement &)=delete
 Broken copy constructor. More...
 
void operator= (const GeneralisedElement &)=delete
 Broken assignment operator. More...
 
Data *& internal_data_pt (const unsigned &i)
 Return a pointer to i-th internal data object. More...
 
Data *const & internal_data_pt (const unsigned &i) const
 Return a pointer to i-th internal data object (const version) More...
 
Data *& external_data_pt (const unsigned &i)
 Return a pointer to i-th external data object. More...
 
Data *const & external_data_pt (const unsigned &i) const
 Return a pointer to i-th external data object (const version) More...
 
unsigned long eqn_number (const unsigned &ieqn_local) const
 
int local_eqn_number (const unsigned long &ieqn_global) const
 
unsigned add_external_data (Data *const &data_pt, const bool &fd=true)
 
bool external_data_fd (const unsigned &i) const
 
void exclude_external_data_fd (const unsigned &i)
 
void include_external_data_fd (const unsigned &i)
 
void flush_external_data ()
 Flush all external data. More...
 
void flush_external_data (Data *const &data_pt)
 Flush the object addressed by data_pt from the external data array. More...
 
unsigned ninternal_data () const
 Return the number of internal data objects. More...
 
unsigned nexternal_data () const
 Return the number of external data objects. More...
 
unsigned ndof () const
 Return the number of equations/dofs in the element. More...
 
void dof_vector (const unsigned &t, Vector< double > &dof)
 Return the vector of dof values at time level t. More...
 
void dof_pt_vector (Vector< double * > &dof_pt)
 Return the vector of pointers to dof values. More...
 
void set_internal_data_time_stepper (const unsigned &i, TimeStepper *const &time_stepper_pt, const bool &preserve_existing_data)
 
void assign_internal_eqn_numbers (unsigned long &global_number, Vector< double * > &Dof_pt)
 
void describe_dofs (std::ostream &out, const std::string &current_string) const
 
void add_internal_value_pt_to_map (std::map< unsigned, double * > &map_of_value_pt)
 
virtual void assign_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual void complete_setup_of_dependencies ()
 
virtual void get_residuals (Vector< double > &residuals)
 
virtual void get_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
virtual void get_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &mass_matrix)
 
virtual void get_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
virtual void get_dresiduals_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam)
 
virtual void get_djacobian_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
 
virtual void get_djacobian_and_dmass_matrix_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam, DenseMatrix< double > &dmass_matrix_dparam)
 
virtual void get_hessian_vector_products (Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
 
virtual void get_inner_products (Vector< std::pair< unsigned, unsigned >> const &history_index, Vector< double > &inner_product)
 
virtual void get_inner_product_vectors (Vector< unsigned > const &history_index, Vector< Vector< double >> &inner_product_vector)
 
virtual void compute_norm (Vector< double > &norm)
 
virtual void compute_norm (double &norm)
 
virtual void get_dof_numbers_for_unknowns (std::list< std::pair< unsigned long, unsigned >> &dof_lookup_list) const
 
- Public Member Functions inherited from oomph::GeomObject
 GeomObject ()
 Default constructor. More...
 
 GeomObject (const unsigned &ndim)
 
 GeomObject (const unsigned &nlagrangian, const unsigned &ndim)
 
 GeomObject (const unsigned &nlagrangian, const unsigned &ndim, TimeStepper *time_stepper_pt)
 
 GeomObject (const GeomObject &dummy)=delete
 Broken copy constructor. More...
 
void operator= (const GeomObject &)=delete
 Broken assignment operator. More...
 
virtual ~GeomObject ()
 (Empty) destructor More...
 
unsigned nlagrangian () const
 Access function to # of Lagrangian coordinates. More...
 
unsigned ndim () const
 Access function to # of Eulerian coordinates. More...
 
void set_nlagrangian_and_ndim (const unsigned &n_lagrangian, const unsigned &n_dim)
 Set # of Lagrangian and Eulerian coordinates. More...
 
TimeStepper *& time_stepper_pt ()
 
TimeSteppertime_stepper_pt () const
 
virtual void position (const double &t, const Vector< double > &zeta, Vector< double > &r) const
 
virtual void dposition (const Vector< double > &zeta, DenseMatrix< double > &drdzeta) const
 
virtual void d2position (const Vector< double > &zeta, RankThreeTensor< double > &ddrdzeta) const
 
virtual void d2position (const Vector< double > &zeta, Vector< double > &r, DenseMatrix< double > &drdzeta, RankThreeTensor< double > &ddrdzeta) const
 
- Public Member Functions inherited from oomph::ElementWithExternalElement
 ElementWithExternalElement ()
 
virtual ~ElementWithExternalElement ()
 The destructor, clean up any allocated memory. More...
 
 ElementWithExternalElement (const ElementWithExternalElement &)=delete
 Broken copy constructor. More...
 
void operator= (const ElementWithExternalElement &)=delete
 Broken assignment operator. More...
 
bool storage_has_been_allocated () const
 Helper function to check if storage has actually been allocated. More...
 
FiniteElement *& external_element_pt (const unsigned &interaction_index, const unsigned &ipt)
 
FiniteElement *const & external_element_pt (const unsigned &interaction_index, const unsigned &ipt) const
 Access function to source element, const version. More...
 
Vector< double > & external_element_local_coord (const unsigned &interaction_index, const unsigned &ipt)
 
Vector< double > const & external_element_local_coord (const unsigned &interaction_index, const unsigned &ipt) const
 Access function to get source element's coords, const version. More...
 
void output_external_elements (std::ostream &outfile, const unsigned &interaction_index)
 
void initialise_external_element_storage ()
 
void flush_all_external_element_storage ()
 Flush the storage for external elements. More...
 
void set_ninteraction (const unsigned &n_interaction)
 
unsigned ninteraction () const
 Return the number of interactions in the element. More...
 
virtual void identify_all_field_data_for_external_interaction (Vector< std::set< FiniteElement * >> const &external_elements_pt, std::set< std::pair< Data *, unsigned >> &paired_interaction_data)
 
virtual void identify_all_geometric_data_for_external_interaction (Vector< std::set< FiniteElement * >> const &external_elements_pt, std::set< Data * > &external_geometric_data_pt)
 
unsigned nexternal_interaction_field_data () const
 
Vector< Data * > external_interaction_field_data_pt () const
 
unsigned nexternal_interaction_geometric_data () const
 
bool add_external_geometric_data ()
 Are we including external geometric data in the element's Jacobian. More...
 
bool add_external_interaction_data ()
 Are we including external data in the element's Jacobian. More...
 
Vector< Data * > external_interaction_geometric_data_pt () const
 
void ignore_external_geometric_data ()
 
void ignore_external_interaction_data ()
 
void include_external_geometric_data ()
 
void include_external_interaction_data ()
 
bool external_geometric_data_is_included () const
 
void describe_local_dofs (std::ostream &out, const std::string &curr_string) const
 

Additional Inherited Members

- Public Types inherited from oomph::FiniteElement
typedef void(* SteadyExactSolutionFctPt) (const Vector< double > &, Vector< double > &)
 
typedef void(* UnsteadyExactSolutionFctPt) (const double &, const Vector< double > &, Vector< double > &)
 
- Public Attributes inherited from oomph::LinearisedAxisymmetricNavierStokesEquations
void(*&)(const double &time, const Vector< double > &x, Vector< double > &fbase_flow_u_fct_pt ()
 Access function for the base flow solution pointer. More...
 
void(*&)(const double &time, const Vector< double > &x, DenseMatrix< double > &fbase_flow_dudx_fct_pt ()
 
void(*&)(const double &time, const Vector< double > &x, Vector< double > &fbase_flow_dudt_fct_pt ()
 
void(*&)(const double &time, const Vector< double > &x, double &fbase_flow_p_fct_pt ()
 Access function for the base flow pressure pointer. More...
 
void(*&)(const double &time, const Vector< double > &x, RankFourTensor< double > &fbase_flow_d_dudx_dX_fct_pt ()
 
void(*&)(const double &time, const Vector< double > &x, Vector< double > &fbody_force_fct_pt ()
 Access function for the body-force pointer. More...
 
double(*&)(const double &time, const Vector< double > &xsource_fct_pt ()
 Access function for the source-function pointer. More...
 
void(*&)(const double &time, const Vector< double > &x, DenseMatrix< double > &fbase_flow_duds_fct_pt ()
 
- Static Public Attributes inherited from oomph::LinearisedAxisymmetricNavierStokesEquations
static Vector< doubleGamma
 Vector to decide whether the stress-divergence form is used or not. More...
 
- Static Public Attributes inherited from oomph::FiniteElement
static double Tolerance_for_singular_jacobian = 1.0e-16
 Tolerance below which the jacobian is considered singular. More...
 
static bool Accept_negative_jacobian = false
 
static bool Suppress_output_while_checking_for_inverted_elements
 
- Static Public Attributes inherited from oomph::GeneralisedElement
static bool Suppress_warning_about_repeated_internal_data
 
static bool Suppress_warning_about_repeated_external_data = true
 
static double Default_fd_jacobian_step = 1.0e-8
 
- Protected Member Functions inherited from oomph::LinearisedAxisymmetricQCrouzeixRaviartElement
double dshape_and_dtest_eulerian_linearised_axi_nst (const Vector< double > &s, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
double dshape_and_dtest_eulerian_at_knot_linearised_axi_nst (const unsigned &ipt, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
void pshape_linearised_axi_nst (const Vector< double > &s, Shape &psi) const
 Compute the pressure shape functions at local coordinate s. More...
 
void pshape_linearised_axi_nst (const Vector< double > &s, Shape &psi, Shape &test) const
 Compute the pressure shape and test functions at local coordinate s. More...
 
double dshape_and_dtest_eulerian_lin_axi_nst (const Vector< double > &s, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
double dshape_and_dtest_eulerian_at_knot_lin_axi_nst (const unsigned &ipt, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
double dshape_and_dtest_eulerian_and_dnodal_coordinates_at_knot_lin_axi_nst (const unsigned &ipt, Shape &psi, DShape &dpsidx, RankFourTensor< double > &d_dpsidx_dX, Shape &test, DShape &dtestdx, RankFourTensor< double > &d_dtestdx_dX, DenseMatrix< double > &djacobian_dX) const
 
void pshape_lin_axi_nst (const Vector< double > &s, Shape &psi) const
 Compute the pressure shape functions at local coordinate s. More...
 
void pshape_lin_axi_nst (const Vector< double > &s, Shape &psi, Shape &test) const
 Compute the pressure shape and test functions at local coordinate s. More...
 
double dshape_and_dtest_eulerian_linearised_axi_nst (const Vector< double > &s, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
double dshape_and_dtest_eulerian_at_knot_linearised_axi_nst (const unsigned &ipt, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
void pshape_linearised_axi_nst (const Vector< double > &s, Shape &psi) const
 Compute the pressure shape functions at local coordinate s. More...
 
void pshape_linearised_axi_nst (const Vector< double > &s, Shape &psi, Shape &test) const
 Compute the pressure shape and test functions at local coordinate s. More...
 
- Protected Member Functions inherited from oomph::LinearisedAxisymmetricNavierStokesEquations
virtual void fill_in_generic_residual_contribution_linearised_axi_nst (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix, unsigned flag)
 
virtual void get_base_flow_d_dudx_dX (const double &time, const unsigned &ipt, const Vector< double > &x, RankFourTensor< double > &result) const
 
void get_body_force_base_flow (const double &time, const unsigned &ipt, const Vector< double > &x, Vector< double > &result)
 
double get_source_base_flow (const double &time, const unsigned &ipt, const Vector< double > &x)
 
void get_body_force_gradient_base_flow (const double &time, const unsigned &ipt, const Vector< double > &x, DenseMatrix< double > &result)
 
void get_source_gradient_base_flow (const double &time, const unsigned &ipt, const Vector< double > &x, Vector< double > &result)
 
virtual void fill_in_generic_residual_contribution_lin_axi_nst (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix, unsigned flag)
 
virtual void fill_in_generic_residual_contribution_linearised_axi_nst (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix, unsigned flag)
 
- Protected Member Functions inherited from oomph::FiniteElement
virtual void assemble_local_to_eulerian_jacobian (const DShape &dpsids, DenseMatrix< double > &jacobian) const
 
virtual void assemble_local_to_eulerian_jacobian2 (const DShape &d2psids, DenseMatrix< double > &jacobian2) const
 
virtual void assemble_eulerian_base_vectors (const DShape &dpsids, DenseMatrix< double > &interpolated_G) const
 
template<unsigned DIM>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
virtual double invert_jacobian_mapping (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
virtual double local_to_eulerian_mapping (const DShape &dpsids, DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
double local_to_eulerian_mapping (const DShape &dpsids, DenseMatrix< double > &inverse_jacobian) const
 
virtual double local_to_eulerian_mapping_diagonal (const DShape &dpsids, DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
virtual void dJ_eulerian_dnodal_coordinates (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<unsigned DIM>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
virtual void d_dshape_eulerian_dnodal_coordinates (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<unsigned DIM>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
virtual void transform_derivatives (const DenseMatrix< double > &inverse_jacobian, DShape &dbasis) const
 
void transform_derivatives_diagonal (const DenseMatrix< double > &inverse_jacobian, DShape &dbasis) const
 
virtual void transform_second_derivatives (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<unsigned DIM>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<unsigned DIM>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
virtual void fill_in_jacobian_from_nodal_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void fill_in_jacobian_from_nodal_by_fd (DenseMatrix< double > &jacobian)
 
virtual void update_before_nodal_fd ()
 
virtual void reset_after_nodal_fd ()
 
virtual void update_in_nodal_fd (const unsigned &i)
 
virtual void reset_in_nodal_fd (const unsigned &i)
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 Zero-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 One-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 Two-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
- Protected Member Functions inherited from oomph::GeneralisedElement
unsigned add_internal_data (Data *const &data_pt, const bool &fd=true)
 
bool internal_data_fd (const unsigned &i) const
 
void exclude_internal_data_fd (const unsigned &i)
 
void include_internal_data_fd (const unsigned &i)
 
void clear_global_eqn_numbers ()
 
void add_global_eqn_numbers (std::deque< unsigned long > const &global_eqn_numbers, std::deque< double * > const &global_dof_pt)
 
int internal_local_eqn (const unsigned &i, const unsigned &j) const
 
int external_local_eqn (const unsigned &i, const unsigned &j)
 
void fill_in_jacobian_from_internal_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_internal_by_fd (DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_external_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_external_by_fd (DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
virtual void update_before_internal_fd ()
 
virtual void reset_after_internal_fd ()
 
virtual void update_in_internal_fd (const unsigned &i)
 
virtual void reset_in_internal_fd (const unsigned &i)
 
virtual void update_before_external_fd ()
 
virtual void reset_after_external_fd ()
 
virtual void update_in_external_fd (const unsigned &i)
 
virtual void reset_in_external_fd (const unsigned &i)
 
virtual void fill_in_contribution_to_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &mass_matrix)
 
virtual void fill_in_contribution_to_dresiduals_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam)
 
virtual void fill_in_contribution_to_djacobian_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
 
virtual void fill_in_contribution_to_djacobian_and_dmass_matrix_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam, DenseMatrix< double > &dmass_matrix_dparam)
 
virtual void fill_in_contribution_to_hessian_vector_products (Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
 
virtual void fill_in_contribution_to_inner_products (Vector< std::pair< unsigned, unsigned >> const &history_index, Vector< double > &inner_product)
 
virtual void fill_in_contribution_to_inner_product_vectors (Vector< unsigned > const &history_index, Vector< Vector< double >> &inner_product_vector)
 
- Protected Member Functions inherited from oomph::ElementWithExternalElement
void assign_internal_and_external_local_eqn_numbers (const bool &store_local_dof_pt)
 
void assign_external_interaction_data_local_eqn_numbers (const bool &store_local_dof_pt)
 
void fill_in_jacobian_from_external_interaction_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void fill_in_jacobian_from_external_interaction_by_fd (DenseMatrix< double > &jacobian)
 
void fill_in_jacobian_from_external_interaction_field_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void fill_in_jacobian_from_external_interaction_field_by_fd (DenseMatrix< double > &jacobian)
 
void fill_in_jacobian_from_external_interaction_geometric_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void fill_in_jacobian_from_external_interaction_geometric_by_fd (DenseMatrix< double > &jacobian)
 
virtual void update_before_external_interaction_field_fd ()
 
virtual void reset_after_external_interaction_field_fd ()
 
virtual void update_in_external_interaction_field_fd (const unsigned &i)
 
virtual void reset_in_external_interaction_field_fd (const unsigned &i)
 
virtual void update_before_external_interaction_geometric_fd ()
 
virtual void reset_after_external_interaction_geometric_fd ()
 
virtual void update_in_external_interaction_geometric_fd (const unsigned &i)
 
virtual void reset_in_external_interaction_geometric_fd (const unsigned &i)
 
- Protected Attributes inherited from oomph::LinearisedAxisymmetricQCrouzeixRaviartElement
Vector< unsignedP_linearised_axi_nst_internal_index
 
Vector< unsignedP_lin_axi_nst_internal_index
 
- Protected Attributes inherited from oomph::LinearisedAxisymmetricNavierStokesEquations
doubleViscosity_Ratio_pt
 
doubleDensity_Ratio_pt
 
doubleRe_pt
 Pointer to global Reynolds number. More...
 
doubleReSt_pt
 Pointer to global Reynolds number x Strouhal number (=Womersley) More...
 
intAzimuthal_Mode_Number_pt
 Pointer to azimuthal mode number k in e^ik(theta) decomposition. More...
 
void(* Base_flow_u_fct_pt )(const double &time, const Vector< double > &x, Vector< double > &result)
 Pointer to base flow solution (velocity components) function. More...
 
void(* Base_flow_dudx_fct_pt )(const double &time, const Vector< double > &x, DenseMatrix< double > &result)
 
bool ALE_is_disabled
 
doubleReInvFr_pt
 
doubleReInvRo_pt
 
Vector< double > * G_pt
 Pointer to global gravity Vector. More...
 
void(* Base_flow_dudt_fct_pt )(const double &time, const Vector< double > &x, Vector< double > &result)
 
void(* Base_flow_p_fct_pt )(const double &time, const Vector< double > &x, double &result)
 Pointer to base flow solution (pressure) function. More...
 
void(* Base_flow_d_dudx_dX_fct_pt )(const double &time, const Vector< double > &x, RankFourTensor< double > &result)
 
void(* Body_force_fct_pt )(const double &time, const Vector< double > &x, Vector< double > &result)
 Pointer to (base flow) body force function. More...
 
double(* Source_fct_pt )(const double &time, const Vector< double > &x)
 Pointer to (base flow) volumetric source function. More...
 
void(* Base_flow_duds_fct_pt )(const double &time, const Vector< double > &x, DenseMatrix< double > &result)
 
Vector< intExternal_H_local_eqn
 
Vector< unsignedExternal_node
 
- Protected Attributes inherited from oomph::FiniteElement
MacroElementMacro_elem_pt
 Pointer to the element's macro element (NULL by default) More...
 
- Protected Attributes inherited from oomph::GeomObject
unsigned NLagrangian
 Number of Lagrangian (intrinsic) coordinates. More...
 
unsigned Ndim
 Number of Eulerian coordinates. More...
 
TimeStepperGeom_object_time_stepper_pt
 
- Protected Attributes inherited from oomph::ElementWithExternalElement
bool Add_external_interaction_data
 Boolean flag to indicate whether to include the external data. More...
 
bool Add_external_geometric_data
 Boolean flag to indicate whether to include the external geometric data. More...
 
Data ** External_interaction_field_data_pt
 
Data ** External_interaction_geometric_data_pt
 
- Static Protected Attributes inherited from oomph::FiniteElement
static const unsigned Default_Initial_Nvalue = 0
 Default value for the number of values at a node. More...
 
static const double Node_location_tolerance = 1.0e-14
 
static const unsigned N2deriv [] = {0, 1, 3, 6}
 
- Static Protected Attributes inherited from oomph::GeneralisedElement
static DenseMatrix< doubleDummy_matrix
 
static std::deque< double * > Dof_pt_deque
 

Detailed Description

////////////////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////////////// Build a LinearisedAxisymmetricQCrouzeixRaviart element that inherits from ElementWithExternalElement so that it can "communicate" with an axisymmetric Navier-Stokes element that provides the base flow velocities and their derivatives w.r.t. global coordinates (r and z)

Constructor & Destructor Documentation

◆ LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement() [1/3]

LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement::LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement ( )
inline

Constructor: call the underlying constructors.

Do not include any external geometric data when computing the element's Jacobian.

190  :
193  {
194  // There are two interactions: the base flow velocities and their
195  // derivatives w.r.t. global coordinates
196  this->set_ninteraction(2);
197 
198  // Do not include any external interaction data when computing
199  // the element's Jacobian
200  ElementWithExternalElement::ignore_external_interaction_data();
201 
204  ElementWithExternalElement::ignore_external_geometric_data();
205  }
void set_ninteraction(const unsigned &n_interaction)
Definition: element_with_external_element.h:178
ElementWithExternalElement()
Definition: element_with_external_element.h:60
LinearisedAxisymmetricQCrouzeixRaviartElement()
Definition: demo_drivers/axisym_navier_stokes/counter_rotating_disks/linearised_axisym_navier_stokes_elements.h:524

References oomph::ElementWithExternalElement::ignore_external_geometric_data(), and oomph::ElementWithExternalElement::ignore_external_interaction_data().

◆ LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement() [2/3]

LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement::LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement ( )
inline

Constructor: call the underlying constructors.

Do not include any external geometric data when computing the element's Jacobian.

60  :
63  {
64  // There are three interactions: the base flow velocities, their
65  // derivatives w.r.t. global coordinates and the base flow pressure
66  this->set_ninteraction(5);
67 
68  // Do not include any external interaction data when computing
69  // the element's Jacobian
70  ElementWithExternalElement::ignore_external_interaction_data();
71 
74  ElementWithExternalElement::ignore_external_geometric_data();
75  }

References oomph::ElementWithExternalElement::ignore_external_geometric_data(), and oomph::ElementWithExternalElement::ignore_external_interaction_data().

◆ LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement() [3/3]

LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement::LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement ( )
inline

Constructor: call the underlying constructors.

Do not include any external geometric data when computing the element's Jacobian.

179  {
180  // There are two interactions: the base flow velocities and their
181  // derivatives w.r.t. global coordinates
182  this->set_ninteraction(2);
183 
184  // Do not include any external interaction data when computing
185  // the element's Jacobian
186  ElementWithExternalElement::ignore_external_interaction_data();
187 
190  ElementWithExternalElement::ignore_external_geometric_data();
191  }

References oomph::ElementWithExternalElement::ignore_external_geometric_data(), and oomph::ElementWithExternalElement::ignore_external_interaction_data().

Member Function Documentation

◆ fill_in_contribution_to_jacobian() [1/3]

void LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement::fill_in_contribution_to_jacobian ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian 
)
inlinevirtual

Compute the element's residual vector and the Jacobian matrix.

Reimplemented from oomph::LinearisedAxisymmetricNavierStokesEquations.

283  {
284  // Get the analytical contribution from the basic
285  // LinearisedAxisymmetricQCrouzeixRaviartElement
287  fill_in_contribution_to_jacobian(residuals,jacobian);
288 
289 //#ifdef USE_FD_FOR_DERIVATIVES_WRT_EXTERNAL_DATA
290 
291  // Get the off-diagonal terms by finite differencing
292  this->fill_in_jacobian_from_external_interaction_by_fd(residuals,jacobian);
293 
294 //#else
295 
296  // Get the off-diagonal terms analytically
297 // this->fill_in_off_diagonal_block_analytic(residuals,jacobian);
298 
299 //#endif
300  }
void fill_in_jacobian_from_external_interaction_by_fd(Vector< double > &residuals, DenseMatrix< double > &jacobian)
Definition: element_with_external_element.h:345
void fill_in_contribution_to_jacobian(Vector< double > &residuals, DenseMatrix< double > &jacobian)
Definition: gen_axisym_advection_diffusion_elements.h:536

References oomph::LinearisedAxisymmetricNavierStokesEquations::fill_in_contribution_to_jacobian().

◆ fill_in_contribution_to_jacobian() [2/3]

void LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement::fill_in_contribution_to_jacobian ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian 
)
inlinevirtual

Compute the element's residual vector and the Jacobian matrix.

Reimplemented from oomph::LinearisedAxisymmetricNavierStokesEquations.

253  {
254  // Get the analytical contribution from the basic patricklinearised element
256  fill_in_contribution_to_jacobian(residuals,jacobian);
257  }

References oomph::LinearisedAxisymmetricNavierStokesEquations::fill_in_contribution_to_jacobian().

◆ fill_in_contribution_to_jacobian() [3/3]

void LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement::fill_in_contribution_to_jacobian ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian 
)
inlinevirtual

Compute the element's residual vector and the Jacobian matrix.

Reimplemented from oomph::ElementWithExternalElement.

268  {
269  // Get the analytical contribution from the basic linearised element
271  fill_in_contribution_to_jacobian(residuals, jacobian);
272 
273  // Get the off-diagonal terms by finite differencing
274  this->fill_in_jacobian_from_external_interaction_by_fd(residuals, jacobian);
275  }

References oomph::LinearisedAxisymmetricNavierStokesEquations::fill_in_contribution_to_jacobian().

◆ get_base_flow_duds()

virtual void LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement::get_base_flow_duds ( const double time,
const unsigned ipt,
const Vector< double > &  x,
DenseMatrix< double > &  result 
) const
inlinevirtual

Overload get_base_flow_duds(...) to return the derivatives of the external element's velocity components w.r.t. local coordinates at the integration point

Reimplemented from oomph::LinearisedAxisymmetricNavierStokesEquations.

218  {
219  // Set interaction index to 4
220  const unsigned interaction = 4;
221 
222  // Get a pointer to the external element that computes the base flow.
223  // We know that it's an axisymmetric Navier-Stokes element.
224  const AxisymmetricQCrouzeixRaviartElement* base_flow_el_pt =
226  external_element_pt(interaction,ipt));
227 
228  // Provide storage for local coordinates in the external element
229  // which correspond to the integration point ipt
230  Vector<double> s_external(2);
231 
232  // Determine local coordinates in the external element which correspond
233  // to the integration point ipt
234  s_external = external_element_local_coord(interaction,ipt);
235 
236  // Loop over velocity components
237  for(unsigned i=0;i<3;i++)
238  {
239  // Loop over coordinate directions and get derivatives of velocity
240  // components from the external element
241  for(unsigned j=0;j<2;j++)
242  {
243  result(i,j)
244  = base_flow_el_pt->interpolated_duds_axi_nst(s_external,i,j);
245  }
246  }
247 
248  } // End of overloaded get_base_flow_duds function
int i
Definition: BiCGSTAB_step_by_step.cpp:9
double interpolated_duds_axi_nst(const Vector< double > &s, const unsigned &i, const unsigned &j) const
Definition: axisym_navier_stokes_elements.h:1023
Definition: axisym_navier_stokes_elements.h:1234
Vector< double > & external_element_local_coord(const unsigned &interaction_index, const unsigned &ipt)
Definition: element_with_external_element.h:136
FiniteElement *& external_element_pt(const unsigned &interaction_index, const unsigned &ipt)
Definition: element_with_external_element.h:107
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2

References i, oomph::AxisymmetricNavierStokesEquations::interpolated_duds_axi_nst(), and j.

◆ get_base_flow_dudt()

virtual void LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement::get_base_flow_dudt ( const double time,
const unsigned ipt,
const Vector< double > &  x,
Vector< double > &  result 
) const
inlinevirtual

Overload get_base_flow_dudt(...) to return the derivatives of the external element's velocity components w.r.t. time at the integration point

Reimplemented from oomph::LinearisedAxisymmetricNavierStokesEquations.

184  {
185  // Set interaction index to 3
186  const unsigned interaction = 3;
187 
188  // Get a pointer to the external element that computes the base flow.
189  // We know that it's an axisymmetric Navier-Stokes element.
190  const AxisymmetricQCrouzeixRaviartElement* base_flow_el_pt =
192  external_element_pt(interaction,ipt));
193 
194  // Provide storage for local coordinates in the external element
195  // which correspond to the integration point ipt
196  Vector<double> s_external(2);
197 
198  // Determine local coordinates in the external element which correspond
199  // to the integration point ipt
200  s_external = external_element_local_coord(interaction,ipt);
201 
202  // Loop over velocity components
203  for(unsigned i=0;i<3;i++)
204  {
205  // Get dudt components from the external element
206  result[i] = base_flow_el_pt->interpolated_dudt_axi_nst(s_external,i);
207  }
208 
209  } // End of overloaded get_base_flow_dudt function
double interpolated_dudt_axi_nst(const Vector< double > &s, const unsigned &i) const
Definition: axisym_navier_stokes_elements.h:1087

References i, and oomph::AxisymmetricNavierStokesEquations::interpolated_dudt_axi_nst().

◆ get_base_flow_dudx() [1/3]

virtual void LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement::get_base_flow_dudx ( const double time,
const unsigned ipt,
const Vector< double > &  x,
DenseMatrix< double > &  result 
) const
inlinevirtual

Overload get_base_flow_dudx(...) to return the derivatives of the external element's velocity components w.r.t. global coordinates at the integration point

Reimplemented from oomph::LinearisedAxisymmetricNavierStokesEquations.

246  {
247  // Set interaction index to 1
248  const unsigned interaction = 1;
249 
250  // Get a pointer to the external element that computes the base flow.
251  // We know that it's an axisymmetric Navier-Stokes element.
252  const AxisymmetricQCrouzeixRaviartElement* base_flow_el_pt =
254  external_element_pt(interaction,ipt));
255 
256  // Provide storage for local coordinates in the external element
257  // which correspond to the integration point ipt
258  Vector<double> s_external(2);
259 
260  // Determine local coordinates in the external element which correspond
261  // to the integration point ipt
262  s_external = external_element_local_coord(interaction,ipt);
263 
264  // Loop over velocity components
265  for(unsigned i=0;i<3;i++)
266  {
267  // Loop over coordinate directions and get derivatives of velocity
268  // components from the external element
269  for(unsigned j=0;j<2;j++)
270  {
271  result(i,j)
272  = base_flow_el_pt->interpolated_dudx_axi_nst(s_external,i,j);
273  }
274  }
275 
276  } // End of overloaded get_base_flow_dudx function
double interpolated_dudx_axi_nst(const Vector< double > &s, const unsigned &i, const unsigned &j) const
Definition: axisym_navier_stokes_elements.h:1055

References i, oomph::AxisymmetricNavierStokesEquations::interpolated_dudx_axi_nst(), and j.

◆ get_base_flow_dudx() [2/3]

virtual void LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement::get_base_flow_dudx ( const double time,
const unsigned ipt,
const Vector< double > &  x,
DenseMatrix< double > &  result 
) const
inlinevirtual

Overload get_base_flow_dudx(...) to return the derivatives of the external element's velocity components w.r.t. global coordinates at the integration point

Reimplemented from oomph::LinearisedAxisymmetricNavierStokesEquations.

116  {
117  // Set interaction index to 1
118  const unsigned interaction = 1;
119 
120  // Get a pointer to the external element that computes the base flow.
121  // We know that it's an axisymmetric Navier-Stokes element.
122  const AxisymmetricQCrouzeixRaviartElement* base_flow_el_pt =
124  external_element_pt(interaction,ipt));
125 
126  // Provide storage for local coordinates in the external element
127  // which correspond to the integration point ipt
128  Vector<double> s_external(2);
129 
130  // Determine local coordinates in the external element which correspond
131  // to the integration point ipt
132  s_external = external_element_local_coord(interaction,ipt);
133 
134  // Loop over velocity components
135  for(unsigned i=0;i<3;i++)
136  {
137  // Loop over coordinate directions and get derivatives of velocity
138  // components from the external element
139  for(unsigned j=0;j<2;j++)
140  {
141  result(i,j)
142  = base_flow_el_pt->interpolated_dudx_axi_nst(s_external,i,j);
143  }
144  }
145 
146  } // End of overloaded get_base_flow_dudx function

References i, oomph::AxisymmetricNavierStokesEquations::interpolated_dudx_axi_nst(), and j.

◆ get_base_flow_dudx() [3/3]

virtual void LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement::get_base_flow_dudx ( const double time,
const unsigned ipt,
const Vector< double > &  x,
DenseMatrix< double > &  result 
) const
inlinevirtual

Overload get_base_flow_dudx(...) to return the derivatives of the external element's velocity components w.r.t. global coordinates at the integration point

Reimplemented from oomph::LinearisedAxisymmetricNavierStokesEquations.

232  {
233  // Set interaction index to 1
234  const unsigned interaction = 1;
235 
236  // Get a pointer to the external element that computes the base flow.
237  // We know that it's an axisymmetric Navier-Stokes element.
238  const AxisymmetricQCrouzeixRaviartElement* base_flow_el_pt =
240  external_element_pt(interaction, ipt));
241 
242  // Provide storage for local coordinates in the external element
243  // which correspond to the integration point ipt
244  Vector<double> s_external(2);
245 
246  // Determine local coordinates in the external element which correspond
247  // to the integration point ipt
248  s_external = external_element_local_coord(interaction, ipt);
249 
250  // Loop over velocity components
251  for (unsigned i = 0; i < 3; i++)
252  {
253  // Loop over coordinate directions and get derivatives of velocity
254  // components from the external element
255  for (unsigned j = 0; j < 2; j++)
256  {
257  result(i, j) =
258  base_flow_el_pt->interpolated_dudx_axi_nst(s_external, i, j);
259  }
260  }
261 
262  } // End of overloaded get_base_flow_dudx function

References i, oomph::AxisymmetricNavierStokesEquations::interpolated_dudx_axi_nst(), and j.

◆ get_base_flow_p()

virtual void LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement::get_base_flow_p ( const double time,
const unsigned ipt,
const Vector< double > &  x,
double result 
) const
inlinevirtual

Overload get_base_flow_p(...) to return the external element's pressure at the integration point

Reimplemented from oomph::LinearisedAxisymmetricNavierStokesEquations.

154  {
155  // Set interaction index to 2
156  const unsigned interaction = 2;
157 
158  // Get a pointer to the external element that computes the base flow.
159  // We know that it's an axisymmetric Navier-Stokes element.
160  const AxisymmetricQCrouzeixRaviartElement* base_flow_el_pt =
162  external_element_pt(interaction,ipt));
163 
164  // Provide storage for local coordinates in the external element
165  // which correspond to the integration point ipt
166  Vector<double> s_external(2);
167 
168  // Determine local coordinates in the external element which correspond
169  // to the integration point ipt
170  s_external = external_element_local_coord(interaction,ipt);
171 
172  // Get the pressure interpolated from the external element
173  result = base_flow_el_pt->interpolated_p_axi_nst(s_external);
174 
175  } // End of overloaded get_base_flow_p function
double interpolated_p_axi_nst(const Vector< double > &s) const
Return FE interpolated pressure at local coordinate s.
Definition: axisym_navier_stokes_elements.h:1001

References oomph::AxisymmetricNavierStokesEquations::interpolated_p_axi_nst().

◆ get_base_flow_u() [1/3]

virtual void LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement::get_base_flow_u ( const double time,
const unsigned ipt,
const Vector< double > &  x,
Vector< double > &  result 
) const
inlinevirtual

Overload get_base_flow_u(...) to return the external element's velocity components at the integration point

Reimplemented from oomph::LinearisedAxisymmetricNavierStokesEquations.

213  {
214  // Set interaction index to 0
215  const unsigned interaction = 0;
216 
217  // Get a pointer to the external element that computes the base flow.
218  // We know that it's an axisymmetric Navier-Stokes element.
219  const AxisymmetricQCrouzeixRaviartElement* base_flow_el_pt =
221  external_element_pt(interaction,ipt));
222 
223  // Provide storage for local coordinates in the external element
224  // which correspond to the integration point ipt
225  Vector<double> s_external(2);
226 
227  // Determine local coordinates in the external element which correspond
228  // to the integration point ipt
229  s_external = external_element_local_coord(interaction,ipt);
230 
231  // Get the three velocity components interpolated from the external element
232  for(unsigned i=0;i<3;i++)
233  {
234  result[i] = base_flow_el_pt->interpolated_u_axi_nst(s_external,i);
235  }
236 
237  } // End of overloaded get_base_flow_u function
void interpolated_u_axi_nst(const Vector< double > &s, Vector< double > &veloc) const
Compute vector of FE interpolated velocity u at local coordinate s.
Definition: axisym_navier_stokes_elements.h:865

References i, and oomph::AxisymmetricNavierStokesEquations::interpolated_u_axi_nst().

◆ get_base_flow_u() [2/3]

virtual void LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement::get_base_flow_u ( const double time,
const unsigned ipt,
const Vector< double > &  x,
Vector< double > &  result 
) const
inlinevirtual

Overload get_base_flow_u(...) to return the external element's velocity components at the integration point

Reimplemented from oomph::LinearisedAxisymmetricNavierStokesEquations.

83  {
84  // Set interaction index to 0
85  const unsigned interaction = 0;
86 
87  // Get a pointer to the external element that computes the base flow.
88  // We know that it's an axisymmetric Navier-Stokes element.
89  const AxisymmetricQCrouzeixRaviartElement* base_flow_el_pt =
91  external_element_pt(interaction,ipt));
92 
93  // Provide storage for local coordinates in the external element
94  // which correspond to the integration point ipt
95  Vector<double> s_external(2);
96 
97  // Determine local coordinates in the external element which correspond
98  // to the integration point ipt
99  s_external = external_element_local_coord(interaction,ipt);
100 
101  // Get the three velocity components interpolated from the external element
102  for(unsigned i=0;i<3;i++)
103  {
104  result[i] = base_flow_el_pt->interpolated_u_axi_nst(s_external,i);
105  }
106 
107  } // End of overloaded get_base_flow_u function

References i, and oomph::AxisymmetricNavierStokesEquations::interpolated_u_axi_nst().

◆ get_base_flow_u() [3/3]

virtual void LinearisedAxisymmetricQCrouzeixRaviartMultiDomainElement::get_base_flow_u ( const double time,
const unsigned ipt,
const Vector< double > &  x,
Vector< double > &  result 
) const
inlinevirtual

Overload get_base_flow_u(...) to return the external element's velocity components at the integration point

Reimplemented from oomph::LinearisedAxisymmetricNavierStokesEquations.

199  {
200  // Set interaction index to 0
201  const unsigned interaction = 0;
202 
203  // Get a pointer to the external element that computes the base flow.
204  // We know that it's an axisymmetric Navier-Stokes element.
205  const AxisymmetricQCrouzeixRaviartElement* base_flow_el_pt =
207  external_element_pt(interaction, ipt));
208 
209  // Provide storage for local coordinates in the external element
210  // which correspond to the integration point ipt
211  Vector<double> s_external(2);
212 
213  // Determine local coordinates in the external element which correspond
214  // to the integration point ipt
215  s_external = external_element_local_coord(interaction, ipt);
216 
217  // Get the three velocity components interpolated from the external element
218  for (unsigned i = 0; i < 3; i++)
219  {
220  result[i] = base_flow_el_pt->interpolated_u_axi_nst(s_external, i);
221  }
222 
223  } // End of overloaded get_base_flow_u function

References i, and oomph::AxisymmetricNavierStokesEquations::interpolated_u_axi_nst().


The documentation for this class was generated from the following file: