Public Member Functions | Private Member Functions | List of all members
oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM > Class Template Reference

#include <generalised_newtonian_refineable_navier_stokes_elements.h>

+ Inheritance diagram for oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >:

Public Member Functions

 PRefineableGeneralisedNewtonianQCrouzeixRaviartElement ()
 Constructor. More...
 
 ~PRefineableGeneralisedNewtonianQCrouzeixRaviartElement ()
 Destructor. More...
 
 PRefineableGeneralisedNewtonianQCrouzeixRaviartElement (const PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM > &dummy)=delete
 Broken copy constructor. More...
 
double p_nst (const unsigned &i) const
 Broken assignment operator. More...
 
double p_nst (const unsigned &t, const unsigned &i) const
 
unsigned npres_nst () const
 // Return number of pressure values More...
 
void fix_pressure (const unsigned &p_dof, const double &p_value)
 Pin p_dof-th pressure dof and set it to value specified by p_value. More...
 
unsigned required_nvalue (const unsigned &n) const
 Number of values (pinned or dofs) required at local node n. More...
 
unsigned ncont_interpolated_values () const
 Number of continuously interpolated values: DIM (velocities) More...
 
void rebuild_from_sons (Mesh *&mesh_pt)
 
unsigned nrecovery_order ()
 
unsigned nvertex_node () const
 Number of vertex nodes in the element. More...
 
Node * vertex_node_pt (const unsigned &j) const
 Pointer to the j-th vertex node in the element. More...
 
double dshape_and_dtest_eulerian_nst (const Vector< double > &s, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
double dshape_and_dtest_eulerian_at_knot_nst (const unsigned &ipt, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
void pshape_nst (const Vector< double > &s, Shape &psi) const
 Pressure shape functions at local coordinate s. More...
 
void pshape_nst (const Vector< double > &s, Shape &psi, Shape &test) const
 Pressure shape and test functions at local coordinte s. More...
 
void get_interpolated_values (const Vector< double > &s, Vector< double > &values)
 
void get_interpolated_values (const unsigned &t, const Vector< double > &s, Vector< double > &values)
 
void further_setup_hanging_nodes ()
 
void further_build ()
 
void further_build ()
 
void further_build ()
 
double dshape_and_dtest_eulerian_nst (const Vector< double > &s, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
double dshape_and_dtest_eulerian_at_knot_nst (const unsigned &ipt, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
double dshape_and_dtest_eulerian_nst (const Vector< double > &s, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
double dshape_and_dtest_eulerian_at_knot_nst (const unsigned &ipt, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
void pshape_nst (const Vector< double > &s, Shape &psi) const
 
void pshape_nst (const Vector< double > &s, Shape &psi, Shape &test) const
 Define the pressure shape and test functions. More...
 
void pshape_nst (const Vector< double > &s, Shape &psi) const
 
void pshape_nst (const Vector< double > &s, Shape &psi, Shape &test) const
 Define the pressure shape and test functions. More...
 
- Public Member Functions inherited from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >
 GeneralisedNewtonianQCrouzeixRaviartElement ()
 Constructor, there are DIM+1 internal values (for the pressure) More...
 
int p_local_eqn (const unsigned &n) const
 Return the local equation numbers for the pressure values. More...
 
void identify_load_data (std::set< std::pair< Data *, unsigned >> &paired_load_data)
 
void identify_pressure_data (std::set< std::pair< Data *, unsigned >> &paired_pressure_data)
 
void output (std::ostream &outfile)
 Redirect output to NavierStokesEquations output. More...
 
void output (std::ostream &outfile, const unsigned &nplot)
 Redirect output to NavierStokesEquations output. More...
 
void output (FILE *file_pt)
 Redirect output to NavierStokesEquations output. More...
 
void output (FILE *file_pt, const unsigned &nplot)
 Redirect output to NavierStokesEquations output. More...
 
void full_output (std::ostream &outfile)
 
void full_output (std::ostream &outfile, const unsigned &nplot)
 
unsigned ndof_types () const
 
void get_dof_numbers_for_unknowns (std::list< std::pair< unsigned long, unsigned >> &dof_lookup_list) const
 
- Public Member Functions inherited from oomph::GeneralisedNewtonianNavierStokesEquations< DIM >
 GeneralisedNewtonianNavierStokesEquations ()
 
const double & re () const
 Reynolds number. More...
 
const double & re_st () const
 Product of Reynolds and Strouhal number (=Womersley number) More...
 
double *& re_pt ()
 Pointer to Reynolds number. More...
 
double *& re_st_pt ()
 Pointer to product of Reynolds and Strouhal number (=Womersley number) More...
 
const double & viscosity_ratio () const
 
double *& viscosity_ratio_pt ()
 Pointer to Viscosity Ratio. More...
 
const double & density_ratio () const
 
double *& density_ratio_pt ()
 Pointer to Density ratio. More...
 
const double & re_invfr () const
 Global inverse Froude number. More...
 
double *& re_invfr_pt ()
 Pointer to global inverse Froude number. More...
 
const Vector< double > & g () const
 Vector of gravitational components. More...
 
Vector< double > *& g_pt ()
 Pointer to Vector of gravitational components. More...
 
NavierStokesBodyForceFctPt & body_force_fct_pt ()
 Access function for the body-force pointer. More...
 
NavierStokesBodyForceFctPt body_force_fct_pt () const
 Access function for the body-force pointer. Const version. More...
 
NavierStokesSourceFctPt & source_fct_pt ()
 Access function for the source-function pointer. More...
 
NavierStokesSourceFctPt source_fct_pt () const
 Access function for the source-function pointer. Const version. More...
 
GeneralisedNewtonianConstitutiveEquation< DIM > *& constitutive_eqn_pt ()
 Access function for the constitutive equation pointer. More...
 
void use_current_strainrate_to_compute_second_invariant ()
 Switch to fully implict evaluation of non-Newtonian const eqn. More...
 
void use_extrapolated_strainrate_to_compute_second_invariant ()
 Use extrapolation for non-Newtonian const eqn. More...
 
double u_nst (const unsigned &n, const unsigned &i) const
 
double u_nst (const unsigned &t, const unsigned &n, const unsigned &i) const
 
virtual unsigned u_index_nst (const unsigned &i) const
 
unsigned n_u_nst () const
 
double du_dt_nst (const unsigned &n, const unsigned &i) const
 
void disable_ALE ()
 
void enable_ALE ()
 
virtual int p_nodal_index_nst () const
 
double pressure_integral () const
 Integral of pressure over element. More...
 
void max_and_min_invariant_and_viscosity (double &min_invariant, double &max_invariant, double &min_viscosity, double &max_viscosity) const
 
double dissipation () const
 Return integral of dissipation over element. More...
 
double dissipation (const Vector< double > &s) const
 Return dissipation at local coordinate s. More...
 
void get_vorticity (const Vector< double > &s, Vector< double > &vorticity) const
 Compute the vorticity vector at local coordinate s. More...
 
double kin_energy () const
 Get integral of kinetic energy over element. More...
 
double d_kin_energy_dt () const
 Get integral of time derivative of kinetic energy over element. More...
 
void strain_rate (const Vector< double > &s, DenseMatrix< double > &strain_rate) const
 Strain-rate tensor: 1/2 (du_i/dx_j + du_j/dx_i) More...
 
void strain_rate (const unsigned &t, const Vector< double > &s, DenseMatrix< double > &strain_rate) const
 
virtual void extrapolated_strain_rate (const Vector< double > &s, DenseMatrix< double > &strain_rate) const
 
virtual void extrapolated_strain_rate (const unsigned &ipt, DenseMatrix< double > &strain_rate) const
 
void get_traction (const Vector< double > &s, const Vector< double > &N, Vector< double > &traction)
 
void get_traction (const Vector< double > &s, const Vector< double > &N, Vector< double > &traction_p, Vector< double > &traction_visc_n, Vector< double > &traction_visc_t)
 
void get_load (const Vector< double > &s, const Vector< double > &N, Vector< double > &load)
 
unsigned nscalar_paraview () const
 
void scalar_value_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot) const
 
std::string scalar_name_paraview (const unsigned &i) const
 
void full_output (std::ostream &outfile)
 
void full_output (std::ostream &outfile, const unsigned &nplot)
 
void output_veloc (std::ostream &outfile, const unsigned &nplot, const unsigned &t)
 
void output_vorticity (std::ostream &outfile, const unsigned &nplot)
 
void output_fct (std::ostream &outfile, const unsigned &nplot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 
void output_fct (std::ostream &outfile, const unsigned &nplot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt)
 
void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 
void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Compute the element's residual Vector. More...
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void fill_in_contribution_to_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
void fill_in_contribution_to_dresiduals_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam)
 Compute the element's residual Vector. More...
 
void fill_in_contribution_to_djacobian_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
 
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)
 
void pin_all_non_pressure_dofs (std::map< Data *, std::vector< int >> &eqn_number_backup)
 Pin all non-pressure dofs and backup eqn numbers. More...
 
void interpolated_u_nst (const Vector< double > &s, Vector< double > &veloc) const
 Compute vector of FE interpolated velocity u at local coordinate s. More...
 
double interpolated_u_nst (const Vector< double > &s, const unsigned &i) const
 Return FE interpolated velocity u[i] at local coordinate s. More...
 
double interpolated_u_nst (const unsigned &t, const Vector< double > &s, const unsigned &i) const
 
virtual double interpolated_p_nst (const Vector< double > &s) const
 Return FE interpolated pressure at local coordinate s. More...
 
double interpolated_p_nst (const unsigned &t, const Vector< double > &s) const
 Return FE interpolated pressure at local coordinate s at time level t. More...
 
void point_output_data (const Vector< double > &s, Vector< double > &data)
 
void get_vorticity (const Vector< double > &s, Vector< double > &vorticity) const
 
void get_vorticity (const Vector< double > &s, Vector< double > &vorticity) const
 Compute 3D vorticity vector at local coordinate s. More...
 
- Public Member Functions inherited from oomph::FSIFluidElement
 FSIFluidElement ()
 Constructor. More...
 
 FSIFluidElement (const FSIFluidElement &)=delete
 Broken copy constructor. More...
 
void operator= (const FSIFluidElement &)=delete
 Broken assignment operator. More...
 
- 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)
 
MacroElement * macro_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 describe_local_dofs (std::ostream &out, const std::string &current_string) const
 
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
 
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 Node * construct_node (const unsigned &n)
 
virtual Node * construct_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 
virtual Node * construct_boundary_node (const unsigned &n)
 
virtual Node * construct_boundary_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 
int get_node_number (Node *const &node_pt) const
 
virtual Node * get_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
 
Data * geom_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_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
 
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 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 void output (const unsigned &t, std::ostream &outfile, const unsigned &n_plot) const
 
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, 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)
 
- 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 ()
 
TimeStepper * time_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::GeneralisedNewtonianTemplateFreeNavierStokesEquationsBase
 GeneralisedNewtonianTemplateFreeNavierStokesEquationsBase ()
 Constructor (empty) More...
 
virtual ~GeneralisedNewtonianTemplateFreeNavierStokesEquationsBase ()
 Virtual destructor (empty) More...
 
- Public Member Functions inherited from oomph::NavierStokesElementWithDiagonalMassMatrices
 NavierStokesElementWithDiagonalMassMatrices ()
 Empty constructor. More...
 
virtual ~NavierStokesElementWithDiagonalMassMatrices ()
 Virtual destructor. More...
 
 NavierStokesElementWithDiagonalMassMatrices (const NavierStokesElementWithDiagonalMassMatrices &)=delete
 Broken copy constructor. More...
 
void operator= (const NavierStokesElementWithDiagonalMassMatrices &)=delete
 Broken assignment operator. More...
 
- Public Member Functions inherited from oomph::RefineableGeneralisedNewtonianNavierStokesEquations< DIM >
 RefineableGeneralisedNewtonianNavierStokesEquations ()
 Constructor. More...
 
virtual Node * pressure_node_pt (const unsigned &n_p)
 
void get_pressure_and_velocity_mass_matrix_diagonal (Vector< double > &press_mass_diag, Vector< double > &veloc_mass_diag, const unsigned &which_one=0)
 
unsigned num_Z2_flux_terms ()
 Number of 'flux' terms for Z2 error estimation. More...
 
void get_Z2_flux (const Vector< double > &s, Vector< double > &flux)
 
void dinterpolated_u_nst_ddata (const Vector< double > &s, const unsigned &i, Vector< double > &du_ddata, Vector< unsigned > &global_eqn_number)
 
- Public Member Functions inherited from oomph::RefineableElement
 RefineableElement ()
 
virtual ~RefineableElement ()
 Destructor, delete the allocated storage for the hanging equations. More...
 
 RefineableElement (const RefineableElement &)=delete
 Broken copy constructor. More...
 
void operator= (const RefineableElement &)=delete
 Broken assignment operator. More...
 
Tree * tree_pt ()
 Access function: Pointer to quadtree representation of this element. More...
 
void set_tree_pt (Tree *my_tree_pt)
 Set pointer to quadtree representation of this element. More...
 
virtual unsigned required_nsons () const
 
bool refinement_is_enabled ()
 Flag to indicate suppression of any refinement. More...
 
void disable_refinement ()
 Suppress of any refinement for this element. More...
 
void enable_refinement ()
 Emnable refinement for this element. More...
 
template<class ELEMENT >
void split (Vector< ELEMENT * > &son_pt) const
 
int local_hang_eqn (Node *const &node_pt, const unsigned &i)
 
virtual void build (Mesh *&mesh_pt, Vector< Node * > &new_node_pt, bool &was_already_built, std::ofstream &new_nodes_file)=0
 
void set_refinement_level (const int &refine_level)
 Set the refinement level. More...
 
unsigned refinement_level () const
 Return the Refinement level. More...
 
void select_for_refinement ()
 Select the element for refinement. More...
 
void deselect_for_refinement ()
 Deselect the element for refinement. More...
 
void select_sons_for_unrefinement ()
 Unrefinement will be performed by merging the four sons of this element. More...
 
void deselect_sons_for_unrefinement ()
 
bool to_be_refined ()
 Has the element been selected for refinement? More...
 
bool sons_to_be_unrefined ()
 Has the element been selected for unrefinement? More...
 
virtual void unbuild ()
 
virtual void deactivate_element ()
 
virtual bool nodes_built ()
 Return true if all the nodes have been built, false if not. More...
 
long number () const
 Element number (for debugging/plotting) More...
 
void set_number (const long &mynumber)
 Set element number (for debugging/plotting) More...
 
virtual Node * interpolating_node_pt (const unsigned &n, const int &value_id)
 
virtual double local_one_d_fraction_of_interpolating_node (const unsigned &n1d, const unsigned &i, const int &value_id)
 
virtual Node * get_interpolating_node_at_local_coordinate (const Vector< double > &s, const int &value_id)
 
virtual unsigned ninterpolating_node (const int &value_id)
 
virtual unsigned ninterpolating_node_1d (const int &value_id)
 
virtual void interpolating_basis (const Vector< double > &s, Shape &psi, const int &value_id) const
 
virtual void check_integrity (double &max_error)=0
 
void identify_field_data_for_interactions (std::set< std::pair< Data *, unsigned >> &paired_field_data)
 
void assign_nodal_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual RefineableElement * root_element_pt ()
 
virtual RefineableElement * father_element_pt () const
 Return a pointer to the father element. More...
 
void get_father_at_refinement_level (unsigned &refinement_level, RefineableElement *&father_at_reflevel_pt)
 
virtual void initial_setup (Tree *const &adopted_father_pt=0, const unsigned &initial_p_order=0)
 
virtual void pre_build (Mesh *&mesh_pt, Vector< Node * > &new_node_pt)
 Pre-build the element. More...
 
virtual void setup_hanging_nodes (Vector< std::ofstream * > &output_stream)
 
void get_dresidual_dnodal_coordinates (RankThreeTensor< double > &dresidual_dnodal_coordinates)
 
unsigned nshape_controlling_nodes ()
 
std::map< Node *, unsigned > shape_controlling_node_lookup ()
 
- Public Member Functions inherited from oomph::ElementWithZ2ErrorEstimator
 ElementWithZ2ErrorEstimator ()
 Default empty constructor. More...
 
 ElementWithZ2ErrorEstimator (const ElementWithZ2ErrorEstimator &)=delete
 Broken copy constructor. More...
 
void operator= (const ElementWithZ2ErrorEstimator &)=delete
 Broken assignment operator. More...
 
virtual unsigned ncompound_fluxes ()
 
virtual void compute_exact_Z2_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_flux_pt, double &error, double &norm)
 
virtual void get_Z2_compound_flux_indices (Vector< unsigned > &flux_index)
 
virtual double geometric_jacobian (const Vector< double > &x)
 
- Public Member Functions inherited from oomph::PRefineableQElement< DIM, 3 >
 PRefineableQElement ()
 Empty constuctor. More...
 

Private Member Functions

void unpin_elemental_pressure_dofs ()
 Unpin all internal pressure dofs. More...
 

Additional Inherited Members

- Public Types inherited from oomph::GeneralisedNewtonianNavierStokesEquations< DIM >
typedef void(* NavierStokesBodyForceFctPt) (const double &time, const Vector< double > &x, Vector< double > &body_force)
 
typedef double(* NavierStokesSourceFctPt) (const double &time, const Vector< double > &x)
 
typedef double(* NavierStokesPressureAdvDiffSourceFctPt) (const Vector< double > &x)
 
- Public Types inherited from oomph::FiniteElement
typedef void(* SteadyExactSolutionFctPt) (const Vector< double > &, Vector< double > &)
 
typedef void(* UnsteadyExactSolutionFctPt) (const double &, const Vector< double > &, Vector< double > &)
 
- Static Public Member Functions inherited from oomph::RefineableGeneralisedNewtonianNavierStokesEquations< DIM >
static void pin_redundant_nodal_pressures (const Vector< GeneralisedElement * > &element_pt)
 
static void unpin_all_pressure_dofs (const Vector< GeneralisedElement * > &element_pt)
 Unpin all pressure dofs in elements listed in vector. More...
 
- Static Public Member Functions inherited from oomph::RefineableElement
static double & max_integrity_tolerance ()
 Max. allowed discrepancy in element integrity check. More...
 
- Static Public Attributes inherited from oomph::GeneralisedNewtonianNavierStokesEquations< DIM >
static Vector< double > Gamma
 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::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >
double dshape_and_dtest_eulerian_at_knot_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
 
double dpshape_and_dptest_eulerian_nst (const Vector< double > &s, Shape &ppsi, DShape &dppsidx, Shape &ptest, DShape &dptestdx) const
 
double dshape_and_dtest_eulerian_at_knot_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
 
double dshape_and_dtest_eulerian_at_knot_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
 
double dpshape_and_dptest_eulerian_nst (const Vector< double > &s, Shape &ppsi, DShape &dppsidx, Shape &ptest, DShape &dptestdx) const
 
double dpshape_and_dptest_eulerian_nst (const Vector< double > &s, Shape &ppsi, DShape &dppsidx, Shape &ptest, DShape &dptestdx) const
 
- Protected Member Functions inherited from oomph::GeneralisedNewtonianNavierStokesEquations< DIM >
virtual void get_body_force_nst (const double &time, const unsigned &ipt, const Vector< double > &s, const Vector< double > &x, Vector< double > &result)
 
virtual void get_body_force_gradient_nst (const double &time, const unsigned &ipt, const Vector< double > &s, const Vector< double > &x, DenseMatrix< double > &d_body_force_dx)
 
virtual double get_source_nst (const double &time, const unsigned &ipt, const Vector< double > &x)
 
virtual void get_source_gradient_nst (const double &time, const unsigned &ipt, const Vector< double > &x, Vector< double > &gradient)
 
virtual void fill_in_generic_dresidual_contribution_nst (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam, DenseMatrix< double > &dmass_matrix_dparam, unsigned flag)
 
void fill_in_contribution_to_hessian_vector_products (Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
 
- Protected Member Functions inherited from oomph::FiniteElement
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 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
 
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)
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
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)
 
virtual void assign_internal_and_external_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual void assign_additional_local_eqn_numbers ()
 
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_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::RefineableGeneralisedNewtonianNavierStokesEquations< DIM >
virtual void pin_elemental_redundant_nodal_pressure_dofs ()
 
void fill_in_generic_residual_contribution_nst (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix, unsigned flag)
 
virtual void get_dresidual_dnodal_coordinates (RankThreeTensor< double > &dresidual_dnodal_coordinates)
 
- Protected Member Functions inherited from oomph::RefineableElement
void assemble_local_to_eulerian_jacobian (const DShape &dpsids, DenseMatrix< double > &jacobian) const
 
void assemble_local_to_eulerian_jacobian2 (const DShape &d2psids, DenseMatrix< double > &jacobian2) const
 
void assemble_eulerian_base_vectors (const DShape &dpsids, DenseMatrix< double > &interpolated_G) const
 
double local_to_eulerian_mapping_diagonal (const DShape &dpsids, DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
void assign_hanging_local_eqn_numbers (const bool &store_local_dof_pt)
 Assign the local equation numbers for hanging node variables. More...
 
virtual void fill_in_jacobian_from_nodal_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
- Static Protected Member Functions inherited from oomph::RefineableElement
static void check_value_id (const int &n_continuously_interpolated_values, const int &value_id)
 
- Protected Attributes inherited from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >
unsigned P_nst_internal_index
 
- Protected Attributes inherited from oomph::GeneralisedNewtonianNavierStokesEquations< DIM >
double * Viscosity_Ratio_pt
 
double * Density_Ratio_pt
 
double * Re_pt
 Pointer to global Reynolds number. More...
 
double * ReSt_pt
 Pointer to global Reynolds number x Strouhal number (=Womersley) More...
 
double * ReInvFr_pt
 
Vector< double > * G_pt
 Pointer to global gravity Vector. More...
 
NavierStokesBodyForceFctPt Body_force_fct_pt
 Pointer to body force function. More...
 
NavierStokesSourceFctPt Source_fct_pt
 Pointer to volumetric source function. More...
 
GeneralisedNewtonianConstitutiveEquation< DIM > * Constitutive_eqn_pt
 Pointer to the generalised Newtonian constitutive equation. More...
 
bool Use_extrapolated_strainrate_to_compute_second_invariant
 
bool ALE_is_disabled
 
- Protected Attributes inherited from oomph::FiniteElement
MacroElement * Macro_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...
 
TimeStepper * Geom_object_time_stepper_pt
 
- Protected Attributes inherited from oomph::RefineableElement
Tree * Tree_pt
 A pointer to a general tree object. More...
 
unsigned Refine_level
 Refinement level. More...
 
bool To_be_refined
 Flag for refinement. More...
 
bool Refinement_is_enabled
 Flag to indicate suppression of any refinement. More...
 
bool Sons_to_be_unrefined
 Flag for unrefinement. More...
 
long Number
 Global element number – for plotting/validation purposes. More...
 
- Static Protected Attributes inherited from oomph::GeneralisedNewtonianNavierStokesEquations< DIM >
static bool Pre_multiply_by_viscosity_ratio = false
 
- 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< double > Dummy_matrix
 
static std::deque< double * > Dof_pt_deque
 
- Static Protected Attributes inherited from oomph::RefineableElement
static double Max_integrity_tolerance = 1.0e-8
 Max. allowed discrepancy in element integrity check. More...
 

Detailed Description

template<unsigned DIM>
class oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >

p-refineable version of Crouzeix Raviart elements. Generic class definitions

Constructor & Destructor Documentation

◆ PRefineableGeneralisedNewtonianQCrouzeixRaviartElement() [1/2]

template<unsigned DIM>
oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement ( )
inline

Constructor.

1077  : RefineableElement(),
1078  RefineableGeneralisedNewtonianNavierStokesEquations<DIM>(),
1079  PRefineableQElement<DIM, 3>(),
1080  GeneralisedNewtonianQCrouzeixRaviartElement<DIM>()
1081  {
1082  // Set the p-order
1083  this->p_order() = 3;
1084 
1085  // Set integration scheme
1086  // (To avoid memory leaks in pre-build and p-refine where new
1087  // integration schemes are created)
1088  this->set_integration_scheme(new GaussLobattoLegendre<DIM, 3>);
1089 
1090  // Resize pressure storage
1091  // (Constructor for QCrouzeixRaviartElement sets up DIM+1 pressure values)
1092  if (this->internal_data_pt(this->P_nst_internal_index)->nvalue() <=
1093  this->npres_nst())
1094  {
1095  this->internal_data_pt(this->P_nst_internal_index)
1096  ->resize(this->npres_nst());
1097  }
1098  else
1099  {
1100  Data* new_data_pt = new Data(this->npres_nst());
1101  delete this->internal_data_pt(this->P_nst_internal_index);
1102  this->internal_data_pt(this->P_nst_internal_index) = new_data_pt;
1103  }
1104  }
oomph::Data::resize
virtual void resize(const unsigned &n_value)
Change (increase) the number of values that may be stored.
Definition: nodes.cc:1002
oomph::FiniteElement::set_integration_scheme
virtual void set_integration_scheme(Integral *const &integral_pt)
Set the spatial integration scheme.
Definition: elements.cc:3210
oomph::GeneralisedElement::internal_data_pt
Data *& internal_data_pt(const unsigned &i)
Return a pointer to i-th internal data object.
Definition: elements.h:622
oomph::GeneralisedNewtonianQCrouzeixRaviartElement::P_nst_internal_index
unsigned P_nst_internal_index
Definition: generalised_newtonian_navier_stokes_elements.h:1313
oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement::npres_nst
unsigned npres_nst() const
// Return number of pressure values
Definition: generalised_newtonian_refineable_navier_stokes_elements.h:1137
oomph::RefineableElement::RefineableElement
RefineableElement()
Definition: refineable_elements.h:188

References oomph::GeneralisedElement::internal_data_pt(), oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::npres_nst(), oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >::P_nst_internal_index, oomph::Data::resize(), and oomph::FiniteElement::set_integration_scheme().

◆ ~PRefineableGeneralisedNewtonianQCrouzeixRaviartElement()

template<unsigned DIM>
oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::~PRefineableGeneralisedNewtonianQCrouzeixRaviartElement ( )
inline

Destructor.

1108  {
1109  delete this->integral_pt();
1110  }
oomph::FiniteElement::integral_pt
Integral *const & integral_pt() const
Return the pointer to the integration scheme (const version)
Definition: elements.h:1963

References oomph::FiniteElement::integral_pt().

◆ PRefineableGeneralisedNewtonianQCrouzeixRaviartElement() [2/2]

template<unsigned DIM>
oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement ( const PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM > &  dummy)
delete

Broken copy constructor.

Member Function Documentation

◆ dshape_and_dtest_eulerian_at_knot_nst() [1/3]

template<unsigned DIM>
double oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::dshape_and_dtest_eulerian_at_knot_nst ( const unsigned &  ipt,
Shape &  psi,
DShape &  dpsidx,
Shape &  test,
DShape &  dtestdx 
) const
inlinevirtual

Velocity shape and test functions and their derivs w.r.t. to global coords at ipt-th integation point (taken from geometry) Return Jacobian of mapping between local and global coordinates.

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

◆ dshape_and_dtest_eulerian_at_knot_nst() [2/3]

double oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< 2 >::dshape_and_dtest_eulerian_at_knot_nst ( const unsigned &  ipt,
Shape &  psi,
DShape &  dpsidx,
Shape &  test,
DShape &  dtestdx 
) const
inlinevirtual

2D Derivatives of the shape functions and test functions w.r.t. to global (Eulerian) coordinates. Return Jacobian of mapping between local and global coordinates.

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

1825  {
1826  // Call the geometrical shape functions and derivatives
1827  double J = this->dshape_eulerian_at_knot(ipt, psi, dpsidx);
1828 
1829  // Loop over the test functions and derivatives and set them equal to the
1830  // shape functions
1831  for (unsigned i = 0; i < nnode_1d() * nnode_1d(); i++)
1832  {
1833  test[i] = psi[i];
1834  dtestdx(i, 0) = dpsidx(i, 0);
1835  dtestdx(i, 1) = dpsidx(i, 1);
1836  }
1837 
1838  // Return the jacobian
1839  return J;
1840  }
i
int i
Definition: BiCGSTAB_step_by_step.cpp:9
J
JacobiRotation< float > J
Definition: Jacobi_makeJacobi.cpp:3
oomph::FiniteElement::dshape_eulerian_at_knot
virtual double dshape_eulerian_at_knot(const unsigned &ipt, Shape &psi, DShape &dpsidx) const
Definition: elements.cc:3325
oomph::FiniteElement::nnode_1d
virtual unsigned nnode_1d() const
Definition: elements.h:2218
test
Definition: indexed_view.cpp:20

References i, and J.

◆ dshape_and_dtest_eulerian_at_knot_nst() [3/3]

double oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< 3 >::dshape_and_dtest_eulerian_at_knot_nst ( const unsigned &  ipt,
Shape &  psi,
DShape &  dpsidx,
Shape &  test,
DShape &  dtestdx 
) const
inlinevirtual

3D Derivatives of the shape functions and test functions w.r.t. to global (Eulerian) coordinates. Return Jacobian of mapping between local and global coordinates.

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

1886  {
1887  // Call the geometrical shape functions and derivatives
1888  double J = this->dshape_eulerian_at_knot(ipt, psi, dpsidx);
1889 
1890  // Loop over the test functions and derivatives and set them equal to the
1891  // shape functions
1892  for (unsigned i = 0; i < nnode_1d() * nnode_1d() * nnode_1d(); i++)
1893  {
1894  test[i] = psi[i];
1895  dtestdx(i, 0) = dpsidx(i, 0);
1896  dtestdx(i, 1) = dpsidx(i, 1);
1897  dtestdx(i, 2) = dpsidx(i, 2);
1898  }
1899 
1900  // Return the jacobian
1901  return J;
1902  }

References i, and J.

◆ dshape_and_dtest_eulerian_nst() [1/3]

template<unsigned DIM>
double oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::dshape_and_dtest_eulerian_nst ( const Vector< double > &  s,
Shape &  psi,
DShape &  dpsidx,
Shape &  test,
DShape &  dtestdx 
) const
inlinevirtual

Velocity shape and test functions and their derivs w.r.t. to global coords at local coordinate s (taken from geometry) Return Jacobian of mapping between local and global coordinates.

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

◆ dshape_and_dtest_eulerian_nst() [2/3]

double oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< 2 >::dshape_and_dtest_eulerian_nst ( const Vector< double > &  s,
Shape &  psi,
DShape &  dpsidx,
Shape &  test,
DShape &  dtestdx 
) const
inlinevirtual

2D Derivatives of the shape functions and test functions w.r.t. to global (Eulerian) coordinates. Return Jacobian of mapping between local and global coordinates.

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

1795  {
1796  // Call the geometrical shape functions and derivatives
1797  double J = this->dshape_eulerian(s, psi, dpsidx);
1798 
1799  // Loop over the test functions and derivatives and set them equal to the
1800  // shape functions
1801  for (unsigned i = 0; i < nnode_1d() * nnode_1d(); i++)
1802  {
1803  test[i] = psi[i];
1804  dtestdx(i, 0) = dpsidx(i, 0);
1805  dtestdx(i, 1) = dpsidx(i, 1);
1806  }
1807 
1808  // Return the jacobian
1809  return J;
1810  }
oomph::FiniteElement::dshape_eulerian
double dshape_eulerian(const Vector< double > &s, Shape &psi, DShape &dpsidx) const
Definition: elements.cc:3298

References i, J, and s.

◆ dshape_and_dtest_eulerian_nst() [3/3]

double oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< 3 >::dshape_and_dtest_eulerian_nst ( const Vector< double > &  s,
Shape &  psi,
DShape &  dpsidx,
Shape &  test,
DShape &  dtestdx 
) const
inlinevirtual

3D Derivatives of the shape functions and test functions w.r.t. to global (Eulerian) coordinates. Return Jacobian of mapping between local and global coordinates.

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

1855  {
1856  // Call the geometrical shape functions and derivatives
1857  double J = this->dshape_eulerian(s, psi, dpsidx);
1858 
1859  // Loop over the test functions and derivatives and set them equal to the
1860  // shape functions
1861  for (unsigned i = 0; i < nnode_1d() * nnode_1d() * nnode_1d(); i++)
1862  {
1863  test[i] = psi[i];
1864  dtestdx(i, 0) = dpsidx(i, 0);
1865  dtestdx(i, 1) = dpsidx(i, 1);
1866  dtestdx(i, 2) = dpsidx(i, 2);
1867  }
1868 
1869  // Return the jacobian
1870  return J;
1871  }

References i, J, and s.

◆ fix_pressure()

template<unsigned DIM>
void oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::fix_pressure ( const unsigned &  p_dof,
const double &  p_value 
)
inlinevirtual

Pin p_dof-th pressure dof and set it to value specified by p_value.

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

1144  {
1145  this->internal_data_pt(this->P_nst_internal_index)->pin(p_dof);
1146  this->internal_data_pt(this->P_nst_internal_index)
1147  ->set_value(p_dof, p_value);
1148  }
oomph::Data::pin
void pin(const unsigned &i)
Pin the i-th stored variable.
Definition: nodes.h:385
oomph::Data::set_value
void set_value(const unsigned &i, const double &value_)
Definition: nodes.h:271

References oomph::GeneralisedElement::internal_data_pt(), oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >::P_nst_internal_index, oomph::Data::pin(), and oomph::Data::set_value().

◆ further_build() [1/3]

void oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< 2 >::further_build ( )
virtual

2D Further build for Crouzeix_Raviart interpolates the internal pressure dofs from father element: Make sure pressure values and dp/ds agree between fathers and sons at the midpoints of the son elements.

Reimplemented from oomph::RefineableGeneralisedNewtonianNavierStokesEquations< DIM >.

1702  {
1703  if (this->tree_pt()->father_pt() != 0)
1704  {
1705  // Call the generic further build (if there is a father)
1706  RefineableGeneralisedNewtonianNavierStokesEquations<2>::further_build();
1707  }
1708  // Now do the PRefineableQElement further_build()
1709  PRefineableQElement<2, 3>::further_build();
1710 
1711  // Resize internal pressure storage
1712  if (this->internal_data_pt(this->P_nst_internal_index)->nvalue() <=
1713  this->npres_nst())
1714  {
1715  this->internal_data_pt(this->P_nst_internal_index)
1716  ->resize(this->npres_nst());
1717  }
1718  else
1719  {
1720  Data* new_data_pt = new Data(this->npres_nst());
1721  delete internal_data_pt(this->P_nst_internal_index);
1722  internal_data_pt(this->P_nst_internal_index) = new_data_pt;
1723  }
1724 
1725  if (this->tree_pt()->father_pt() != 0)
1726  {
1727  // Pointer to my father (in C-R element impersonation)
1728  PRefineableGeneralisedNewtonianQCrouzeixRaviartElement<2>*
1729  father_element_pt = dynamic_cast<
1730  PRefineableGeneralisedNewtonianQCrouzeixRaviartElement<2>*>(
1731  quadtree_pt()->father_pt()->object_pt());
1732 
1733  // If element has same p-order as father then do the projection problem
1734  // (called after h-refinement)
1735  if (father_element_pt->p_order() == this->p_order())
1736  {
1737  using namespace QuadTreeNames;
1738 
1739  // What type of son am I? Ask my quadtree representation...
1740  int son_type = quadtree_pt()->son_type();
1741 
1742  Vector<double> s_father(2);
1743 
1744  // Son midpoint is located at the following coordinates in father
1745  // element:
1746  switch (son_type)
1747  {
1748  case SW:
1749  // South west son
1750  s_father[0] = -0.5;
1751  s_father[1] = -0.5;
1752  break;
1753  case SE:
1754  // South east son
1755  s_father[0] = 0.5;
1756  s_father[1] = -0.5;
1757  break;
1758  case NE:
1759  // North east son
1760  s_father[0] = 0.5;
1761  s_father[1] = 0.5;
1762  break;
1763  case NW:
1764  // North west son
1765  s_father[0] = -0.5;
1766  s_father[1] = 0.5;
1767  break;
1768  default:
1769  throw OomphLibError("Invalid son type in",
1770  OOMPH_CURRENT_FUNCTION,
1771  OOMPH_EXCEPTION_LOCATION);
1772  break;
1773  }
1774 
1775  // Get pressure value in father element
1776  double press = father_element_pt->interpolated_p_nst(s_father);
1777 
1778  // Reset all pressures to zero
1779  for (unsigned p = 0; p < this->npres_nst(); p++)
1780  {
1781  internal_data_pt(this->P_nst_internal_index)->set_value(p, 0.0);
1782  }
1783 
1784  // Set pressure values from father (BENFLAG: projection problem hack)
1785  if (this->npres_nst() == 1)
1786  {
1787  internal_data_pt(this->P_nst_internal_index)->set_value(0, press);
1788  }
1789  else if (this->npres_nst() == 3)
1790  {
1791  internal_data_pt(this->P_nst_internal_index)->set_value(0, press);
1792  internal_data_pt(this->P_nst_internal_index)->set_value(1, press);
1793  internal_data_pt(this->P_nst_internal_index)->set_value(2, press);
1794  }
1795  else
1796  {
1797  internal_data_pt(this->P_nst_internal_index)->set_value(0, press);
1798  internal_data_pt(this->P_nst_internal_index)->set_value(1, press);
1799  internal_data_pt(this->P_nst_internal_index)->set_value(2, press);
1800  internal_data_pt(this->P_nst_internal_index)->set_value(3, press);
1801  }
1802  } // Otherwise this is called after p-refinement
1803  }
1804  }
p
float * p
Definition: Tutorial_Map_using.cpp:9
oomph::RefineableElement::tree_pt
Tree * tree_pt()
Access function: Pointer to quadtree representation of this element.
Definition: refineable_elements.h:211
oomph::RefineableElement::father_element_pt
virtual RefineableElement * father_element_pt() const
Return a pointer to the father element.
Definition: refineable_elements.h:539
oomph::RefineableGeneralisedNewtonianNavierStokesEquations::further_build
void further_build()
Further build, pass the pointers down to the sons.
Definition: generalised_newtonian_refineable_navier_stokes_elements.h:185
oomph::QuadTreeNames::SE
@ SE
Definition: quadtree.h:57
oomph::QuadTreeNames::NW
@ NW
Definition: quadtree.h:58
oomph::QuadTreeNames::NE
@ NE
Definition: quadtree.h:59
oomph::QuadTreeNames::SW
@ SW
Definition: quadtree.h:56
OOMPH_EXCEPTION_LOCATION
#define OOMPH_EXCEPTION_LOCATION
Definition: oomph_definitions.h:61
OOMPH_CURRENT_FUNCTION
#define OOMPH_CURRENT_FUNCTION
Definition: oomph_definitions.h:86

References oomph::RefineableGeneralisedNewtonianNavierStokesEquations< DIM >::further_build(), oomph::GeneralisedNewtonianNavierStokesEquations< DIM >::interpolated_p_nst(), oomph::QuadTreeNames::NE, oomph::QuadTreeNames::NW, OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, p, oomph::QuadTreeNames::SE, and oomph::QuadTreeNames::SW.

◆ further_build() [2/3]

void oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< 3 >::further_build ( )
virtual

2D Further build for Crouzeix_Raviart interpolates the internal pressure dofs from father element: Make sure pressure values and dp/ds agree between fathers and sons at the midpoints of the son elements.

Reimplemented from oomph::RefineableGeneralisedNewtonianNavierStokesEquations< DIM >.

1815  {
1816  if (this->tree_pt()->father_pt() != 0)
1817  {
1818  // Call the generic further build (if there is a father)
1819  RefineableGeneralisedNewtonianNavierStokesEquations<3>::further_build();
1820  }
1821  // Now do the PRefineableQElement further_build()
1822  PRefineableQElement<3, 3>::further_build();
1823 
1824  // Resize internal pressure storage
1825  if (this->internal_data_pt(this->P_nst_internal_index)->nvalue() <=
1826  this->npres_nst())
1827  {
1828  this->internal_data_pt(this->P_nst_internal_index)
1829  ->resize(this->npres_nst());
1830  }
1831  else
1832  {
1833  Data* new_data_pt = new Data(this->npres_nst());
1834  delete internal_data_pt(this->P_nst_internal_index);
1835  internal_data_pt(this->P_nst_internal_index) = new_data_pt;
1836  }
1837 
1838  if (this->tree_pt()->father_pt() != 0)
1839  {
1840  // Pointer to my father (in C-R element impersonation)
1841  PRefineableGeneralisedNewtonianQCrouzeixRaviartElement<3>*
1842  father_element_pt = dynamic_cast<
1843  PRefineableGeneralisedNewtonianQCrouzeixRaviartElement<3>*>(
1844  octree_pt()->father_pt()->object_pt());
1845 
1846  // If element has same p-order as father then do the projection problem
1847  // (called after h-refinement)
1848  if (father_element_pt->p_order() == this->p_order())
1849  {
1850  using namespace OcTreeNames;
1851 
1852  // What type of son am I? Ask my quadtree representation...
1853  int son_type = octree_pt()->son_type();
1854 
1855  Vector<double> s_father(3);
1856 
1857 
1858  // Son midpoint is located at the following coordinates in father
1859  // element:
1860  for (unsigned i = 0; i < 3; i++)
1861  {
1862  s_father[i] = 0.5 * OcTree::Direction_to_vector[son_type][i];
1863  }
1864 
1865  // Get pressure value in father element
1866  double press = father_element_pt->interpolated_p_nst(s_father);
1867 
1868  // Reset all pressures to zero
1869  for (unsigned p = 0; p < this->npres_nst(); p++)
1870  {
1871  internal_data_pt(this->P_nst_internal_index)->set_value(p, 0.0);
1872  }
1873 
1874  // Set pressure values from father (BENFLAG: projection problem hack)
1875  if (this->npres_nst() == 1)
1876  {
1877  internal_data_pt(this->P_nst_internal_index)->set_value(0, press);
1878  }
1879  else
1880  {
1881  internal_data_pt(this->P_nst_internal_index)->set_value(0, press);
1882  internal_data_pt(this->P_nst_internal_index)->set_value(1, press);
1883  internal_data_pt(this->P_nst_internal_index)->set_value(2, press);
1884  internal_data_pt(this->P_nst_internal_index)->set_value(3, press);
1885  internal_data_pt(this->P_nst_internal_index)->set_value(4, press);
1886  internal_data_pt(this->P_nst_internal_index)->set_value(5, press);
1887  internal_data_pt(this->P_nst_internal_index)->set_value(6, press);
1888  internal_data_pt(this->P_nst_internal_index)->set_value(7, press);
1889  }
1890  } // Otherwise this is called after p-refinement
1891  }
1892  }
oomph::OcTree::Direction_to_vector
static Vector< Vector< int > > Direction_to_vector
Definition: octree.h:353

References oomph::OcTree::Direction_to_vector, oomph::RefineableGeneralisedNewtonianNavierStokesEquations< DIM >::further_build(), i, oomph::GeneralisedNewtonianNavierStokesEquations< DIM >::interpolated_p_nst(), and p.

◆ further_build() [3/3]

template<unsigned DIM>
void oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::further_build ( )
virtual

Further build for Crouzeix_Raviart interpolates the internal pressure dofs from father element: Make sure pressure values and dp/ds agree between fathers and sons at the midpoints of the son elements. This must be specialised for each dimension.

Reimplemented from oomph::RefineableGeneralisedNewtonianNavierStokesEquations< DIM >.

◆ further_setup_hanging_nodes()

template<unsigned DIM>
void oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::further_setup_hanging_nodes ( )
inlinevirtual

Perform additional hanging node procedures for variables that are not interpolated by all nodes. Empty

Reimplemented from oomph::RefineableElement.

1289 {}

◆ get_interpolated_values() [1/2]

template<unsigned DIM>
void oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::get_interpolated_values ( const unsigned &  t,
const Vector< double > &  s,
Vector< double > &  values 
)
inlinevirtual

Get all function values [u,v..,p] at previous timestep t (t=0: present; t>0: previous timestep).

Note: Given the generality of the interface (this function is usually called from black-box documentation or interpolation routines), the values Vector sets its own size in here.

Note: No pressure history is kept, so pressure is always the current value.

Implements oomph::RefineableElement.

1258  {
1259  // Set size of Vector: u,v,p
1260  values.resize(DIM);
1261 
1262  // Initialise
1263  for (unsigned i = 0; i < DIM; i++)
1264  {
1265  values[i] = 0.0;
1266  }
1267 
1268  // Find out how many nodes there are
1269  unsigned n_node = this->nnode();
1270 
1271  // Shape functions
1272  Shape psif(n_node);
1273  this->shape(s, psif);
1274 
1275  // Calculate velocities: values[0],...
1276  for (unsigned i = 0; i < DIM; i++)
1277  {
1278  // Get the nodal index at which the i-th velocity component is stored
1279  unsigned u_nodal_index = this->u_index_nst(i);
1280  for (unsigned l = 0; l < n_node; l++)
1281  {
1282  values[i] += this->nodal_value(t, l, u_nodal_index) * psif[l];
1283  }
1284  }
1285  }
oomph::FiniteElement::nodal_value
double nodal_value(const unsigned &n, const unsigned &i) const
Definition: elements.h:2593
oomph::FiniteElement::shape
virtual void shape(const Vector< double > &s, Shape &psi) const =0
oomph::FiniteElement::nnode
unsigned nnode() const
Return the number of nodes.
Definition: elements.h:2210
oomph::GeneralisedNewtonianNavierStokesEquations::u_index_nst
virtual unsigned u_index_nst(const unsigned &i) const
Definition: generalised_newtonian_navier_stokes_elements.h:580
DIM
#define DIM
Definition: linearised_navier_stokes_elements.h:44

References DIM, i, oomph::FiniteElement::nnode(), oomph::FiniteElement::nodal_value(), oomph::FiniteElement::shape(), and oomph::GeneralisedNewtonianNavierStokesEquations< DIM >::u_index_nst().

◆ get_interpolated_values() [2/2]

template<unsigned DIM>
void oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::get_interpolated_values ( const Vector< double > &  s,
Vector< double > &  values 
)
inlinevirtual

Get the function value u in Vector. Note: Given the generality of the interface (this function is usually called from black-box documentation or interpolation routines), the values Vector sets its own size in here.

Reimplemented from oomph::RefineableElement.

1235  {
1236  // Set size of Vector: u,v,p and initialise to zero
1237  values.resize(DIM, 0.0);
1238 
1239  // Calculate velocities: values[0],...
1240  for (unsigned i = 0; i < DIM; i++)
1241  {
1242  values[i] = this->interpolated_u_nst(s, i);
1243  }
1244  }
oomph::GeneralisedNewtonianNavierStokesEquations::interpolated_u_nst
void interpolated_u_nst(const Vector< double > &s, Vector< double > &veloc) const
Compute vector of FE interpolated velocity u at local coordinate s.
Definition: generalised_newtonian_navier_stokes_elements.h:1093

References DIM, i, and oomph::GeneralisedNewtonianNavierStokesEquations< DIM >::interpolated_u_nst().

◆ ncont_interpolated_values()

template<unsigned DIM>
unsigned oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::ncont_interpolated_values ( ) const
inlinevirtual

Number of continuously interpolated values: DIM (velocities)

Implements oomph::RefineableElement.

1157  {
1158  return DIM;
1159  }

References DIM.

◆ npres_nst()

template<unsigned DIM>
unsigned oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::npres_nst ( ) const
inlinevirtual

// Return number of pressure values

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

1138  {
1139  return (this->p_order() - 2) * (this->p_order() - 2);
1140  }

Referenced by oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement(), oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::rebuild_from_sons(), and oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::unpin_elemental_pressure_dofs().

◆ nrecovery_order()

template<unsigned DIM>
unsigned oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::nrecovery_order ( )
inlinevirtual

Order of recovery shape functions for Z2 error estimation:

  • Same order as shape functions.
  • Constant recovery order, since recovery order of the first element is used for the whole mesh.

Implements oomph::ElementWithZ2ErrorEstimator.

1186  {
1187  return 3;
1188  }

◆ nvertex_node()

template<unsigned DIM>
unsigned oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::nvertex_node ( ) const
inlinevirtual

Number of vertex nodes in the element.

Implements oomph::ElementWithZ2ErrorEstimator.

1192  {
1193  return GeneralisedNewtonianQCrouzeixRaviartElement<DIM>::nvertex_node();
1194  }
oomph::FiniteElement::nvertex_node
virtual unsigned nvertex_node() const
Definition: elements.h:2491

References oomph::FiniteElement::nvertex_node().

◆ p_nst() [1/2]

template<unsigned DIM>
double oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::p_nst ( const unsigned &  i) const
inlinevirtual

Broken assignment operator.

Return the i-th pressure value (Discontinous pressure interpolation – no need to cater for hanging nodes).

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

1127  {
1128  return this->internal_data_pt(this->P_nst_internal_index)->value(i);
1129  }
oomph::Data::value
double value(const unsigned &i) const
Definition: nodes.h:293

References oomph::GeneralisedElement::internal_data_pt(), oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >::P_nst_internal_index, and oomph::Data::value().

◆ p_nst() [2/2]

template<unsigned DIM>
double oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::p_nst ( const unsigned &  t,
const unsigned &  i 
) const
inlinevirtual

Return the i-th pressure value (Discontinous pressure interpolation – no need to cater for hanging nodes).

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

1132  {
1133  return this->internal_data_pt(this->P_nst_internal_index)->value(t, i);
1134  }

References i, oomph::GeneralisedElement::internal_data_pt(), oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >::P_nst_internal_index, and oomph::Data::value().

◆ pshape_nst() [1/6]

template<unsigned DIM>
void oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::pshape_nst ( const Vector< double > &  s,
Shape &  psi 
) const
inlinevirtual

Pressure shape functions at local coordinate s.

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

◆ pshape_nst() [2/6]

void oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< 2 >::pshape_nst ( const Vector< double > &  s,
Shape &  psi 
) const
inlinevirtual

2D : Pressure shape functions

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

1911  {
1912  unsigned npres = this->npres_nst();
1913  if (npres == 1)
1914  {
1915  psi[0] = 1.0;
1916  }
1917  else
1918  {
1919  // Get number of pressure modes
1920  unsigned npres_1d = (int)std::sqrt((double)npres);
1921 
1922  // Local storage
1923  // Call the one-dimensional modal shape functions
1924  OneDimensionalModalShape psi1(npres_1d, s[0]);
1925  OneDimensionalModalShape psi2(npres_1d, s[1]);
1926 
1927  // Now let's loop over the nodal points in the element
1928  // s1 is the "x" coordinate, s2 the "y"
1929  for (unsigned i = 0; i < npres_1d; i++)
1930  {
1931  for (unsigned j = 0; j < npres_1d; j++)
1932  {
1933  // Multiply the two 1D functions together to get the 2D function
1934  psi[i * npres_1d + j] = psi2[i] * psi1[j];
1935  }
1936  }
1937  }
1938  }
sqrt
AnnoyingScalar sqrt(const AnnoyingScalar &x)
Definition: AnnoyingScalar.h:134
s
RealScalar s
Definition: level1_cplx_impl.h:130
int
return int(ret)+1
j
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2

References i, int(), j, s, and sqrt().

◆ pshape_nst() [3/6]

void oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< 3 >::pshape_nst ( const Vector< double > &  s,
Shape &  psi 
) const
inlinevirtual

3D : Pressure shape functions

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

1966  {
1967  unsigned npres = this->npres_nst();
1968  if (npres == 1)
1969  {
1970  psi[0] = 1.0;
1971  }
1972  else
1973  {
1974  // Get number of pressure modes
1975  unsigned npres_1d = (int)std::sqrt((double)npres);
1976 
1977  // Local storage
1978  // Call the one-dimensional modal shape functions
1979  OneDimensionalModalShape psi1(npres_1d, s[0]);
1980  OneDimensionalModalShape psi2(npres_1d, s[1]);
1981  OneDimensionalModalShape psi3(npres_1d, s[2]);
1982 
1983  // Now let's loop over the nodal points in the element
1984  // s1 is the "x" coordinate, s2 the "y"
1985  for (unsigned i = 0; i < npres_1d; i++)
1986  {
1987  for (unsigned j = 0; j < npres_1d; j++)
1988  {
1989  for (unsigned k = 0; k < npres_1d; k++)
1990  {
1991  // Multiply the two 1D functions together to get the 2D function
1992  psi[i * npres_1d * npres_1d + j * npres_1d + k] =
1993  psi3[i] * psi2[j] * psi1[k];
1994  }
1995  }
1996  }
1997  }
1998  }
k
char char char int int * k
Definition: level2_impl.h:374

References i, int(), j, k, s, and sqrt().

◆ pshape_nst() [4/6]

template<unsigned DIM>
void oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::pshape_nst ( const Vector< double > &  s,
Shape &  psi,
Shape &  test 
) const
inlinevirtual

Pressure shape and test functions at local coordinte s.

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

◆ pshape_nst() [5/6]

void oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< 2 >::pshape_nst ( const Vector< double > &  s,
Shape &  psi,
Shape &  test 
) const
inlinevirtual

Define the pressure shape and test functions.

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

1944  {
1945  // Call the pressure shape functions
1946  pshape_nst(s, psi);
1947 
1948  // Loop over the test functions and set them equal to the shape functions
1949  if (this->npres_nst() == 1)
1950  {
1951  test[0] = psi[0];
1952  }
1953  else
1954  {
1955  for (unsigned i = 0; i < this->npres_nst(); i++) test[i] = psi[i];
1956  }
1957  }
oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement::pshape_nst
void pshape_nst(const Vector< double > &s, Shape &psi) const
Pressure shape functions at local coordinate s.

References i, and s.

◆ pshape_nst() [6/6]

void oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< 3 >::pshape_nst ( const Vector< double > &  s,
Shape &  psi,
Shape &  test 
) const
inlinevirtual

Define the pressure shape and test functions.

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

2004  {
2005  // Call the pressure shape functions
2006  pshape_nst(s, psi);
2007 
2008  // Loop over the test functions and set them equal to the shape functions
2009  if (this->npres_nst() == 1)
2010  {
2011  test[0] = psi[0];
2012  }
2013  else
2014  {
2015  for (unsigned i = 0; i < this->npres_nst(); i++) test[i] = psi[i];
2016  }
2017  }

References i, and s.

◆ rebuild_from_sons()

template<unsigned DIM>
void oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::rebuild_from_sons ( Mesh *&  mesh_pt)
inlinevirtual

Rebuild from sons: Reconstruct pressure from the (merged) sons This must be specialised for each dimension.

Implements oomph::RefineableElement.

1164  {
1165  // Do p-refineable version
1166  PRefineableQElement<DIM, 3>::rebuild_from_sons(mesh_pt);
1167  // Do Crouzeix-Raviart version
1168  // Need to reconstruct pressure manually!
1169  for (unsigned p = 0; p < npres_nst(); p++)
1170  {
1171  // BENFLAG: Set to zero for now -- don't do projection problem yet
1172  this->internal_data_pt(this->P_nst_internal_index)->set_value(p, 0.0);
1173  }
1174  }

References oomph::GeneralisedElement::internal_data_pt(), oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::npres_nst(), p, oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >::P_nst_internal_index, and oomph::Data::set_value().

◆ required_nvalue()

template<unsigned DIM>
unsigned oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::required_nvalue ( const unsigned &  n) const
inlinevirtual

Number of values (pinned or dofs) required at local node n.

Number of values (pinned or dofs) required at node n.

Reimplemented from oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >.

1151  {
1152  return DIM;
1153  }

References DIM.

◆ unpin_elemental_pressure_dofs()

template<unsigned DIM>
void oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::unpin_elemental_pressure_dofs ( )
inlineprivatevirtual

Unpin all internal pressure dofs.

Implements oomph::RefineableGeneralisedNewtonianNavierStokesEquations< DIM >.

1064  {
1065  unsigned n_pres = this->npres_nst();
1066  n_pres = this->internal_data_pt(this->P_nst_internal_index)->nvalue();
1067  // loop over pressure dofs and unpin them
1068  for (unsigned l = 0; l < n_pres; l++)
1069  {
1070  this->internal_data_pt(this->P_nst_internal_index)->unpin(l);
1071  }
1072  }
oomph::Data::unpin
void unpin(const unsigned &i)
Unpin the i-th stored variable.
Definition: nodes.h:391
oomph::Data::nvalue
unsigned nvalue() const
Return number of values stored in data object (incl pinned ones).
Definition: nodes.h:483

References oomph::GeneralisedElement::internal_data_pt(), oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::npres_nst(), oomph::Data::nvalue(), oomph::GeneralisedNewtonianQCrouzeixRaviartElement< DIM >::P_nst_internal_index, and oomph::Data::unpin().

◆ vertex_node_pt()

template<unsigned DIM>
Node* oomph::PRefineableGeneralisedNewtonianQCrouzeixRaviartElement< DIM >::vertex_node_pt ( const unsigned &  j) const
inlinevirtual

Pointer to the j-th vertex node in the element.

Implements oomph::ElementWithZ2ErrorEstimator.

1198  {
1199  return GeneralisedNewtonianQCrouzeixRaviartElement<DIM>::vertex_node_pt(
1200  j);
1201  }
oomph::FiniteElement::vertex_node_pt
virtual Node * vertex_node_pt(const unsigned &j) const
Definition: elements.h:2500

References j, and oomph::FiniteElement::vertex_node_pt().

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