oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM > Class Template Reference

#include <double_buoyant_navier_stokes_elements.h>

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

Public Member Functions

 RefineableDoubleBuoyantQCrouzeixRaviartElement ()
 
void unfix_pressure (const unsigned &p_dof)
 UnPin p_dof-th pressure dof. More...
 
unsigned required_nvalue (const unsigned &n) const
 
const doublekm () const
 Access function for the transfer constant. More...
 
double *& km_pt ()
 Access function for the pointer to transfer constant. More...
 
const doublen () const
 Access function for the number of monomers in the micelle. More...
 
double *& n_pt ()
 Access function for the pointer to the number of monomers in the micelle. More...
 
void disable_ALE ()
 Final override for disable ALE. More...
 
void enable_ALE ()
 Final override for enable ALE. More...
 
void get_reaction_adv_diff_react (const unsigned &ipt, const Vector< double > &C, Vector< double > &R) const
 
void get_reaction_deriv_adv_diff_react (const unsigned &ipt, const Vector< double > &C, DenseMatrix< double > &dRdC) const
 
void compute_norm (double &norm)
 
void output (std::ostream &outfile)
 Overload the standard output function with the broken default. More...
 
void output (std::ostream &outfile, const unsigned &nplot)
 
void output (FILE *file_pt)
 C-style output function: Broken default. More...
 
void output (FILE *file_pt, const unsigned &n_plot)
 C-style output function: Broken default. More...
 
void output_fct (std::ostream &outfile, const unsigned &Nplot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 Output function for an exact solution: Broken default. More...
 
void output_fct (std::ostream &outfile, const unsigned &Nplot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt)
 
unsigned c_index_adv_diff_react (const unsigned &i) const
 
unsigned nvertex_node () const
 geometric element. More...
 
Nodevertex_node_pt (const unsigned &j) const
 
unsigned ncont_interpolated_values () const
 
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 rebuild_from_sons (Mesh *&mesh_pt)
 
void further_build ()
 
unsigned nrecovery_order ()
 The recovery order is that of the NavierStokes elements. More...
 
unsigned num_Z2_flux_terms ()
 
void get_Z2_flux (const Vector< double > &s, Vector< double > &flux)
 
unsigned ncompound_fluxes ()
 
void get_Z2_compound_flux_indices (Vector< unsigned > &flux_index)
 
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 get_wind_adv_diff_react (const unsigned &ipt, const Vector< double > &s, const Vector< double > &x, Vector< double > &wind) const
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 
void fill_in_off_diagonal_jacobian_blocks_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
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 integrated_C_and_M (double &int_C, double &int_M)
 
- Public Member Functions inherited from oomph::RefineableQAdvectionDiffusionReactionElement< 2, DIM, 3 >
 RefineableQAdvectionDiffusionReactionElement ()
 Empty Constructor: More...
 
 RefineableQAdvectionDiffusionReactionElement (const RefineableQAdvectionDiffusionReactionElement< NREAGENT, DIM, NNODE_1D > &dummy)=delete
 Broken copy constructor. More...
 
void operator= (const RefineableQAdvectionDiffusionReactionElement< NREAGENT, DIM, NNODE_1D > &)=delete
 Broken assignment operator. More...
 
unsigned ncont_interpolated_values () const
 Number of continuously interpolated values: NREAGENT. More...
 
unsigned nvertex_node () const
 Number of vertex nodes in the element. More...
 
Nodevertex_node_pt (const unsigned &j) const
 Pointer to the j-th vertex node in the element. More...
 
void rebuild_from_sons (Mesh *&mesh_pt)
 Rebuild from sons: empty. More...
 
unsigned nrecovery_order ()
 
void further_setup_hanging_nodes ()
 
- Public Member Functions inherited from oomph::QAdvectionDiffusionReactionElement< NREAGENT, DIM, NNODE_1D >
 QAdvectionDiffusionReactionElement ()
 
 QAdvectionDiffusionReactionElement (const QAdvectionDiffusionReactionElement< NREAGENT, DIM, NNODE_1D > &dummy)=delete
 Broken copy constructor. More...
 
void operator= (const QAdvectionDiffusionReactionElement< NREAGENT, DIM, NNODE_1D > &)=delete
 Broken assignment operator. More...
 
unsigned required_nvalue (const unsigned &n) const
 
void output (std::ostream &outfile)
 
void output (std::ostream &outfile, const unsigned &n_plot)
 
void output (FILE *file_pt)
 
void output (FILE *file_pt, const unsigned &n_plot)
 
void output_fct (std::ostream &outfile, const unsigned &n_plot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 
void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt)
 
- Public Member Functions inherited from oomph::AdvectionDiffusionReactionEquations< NREAGENT, DIM >
 AdvectionDiffusionReactionEquations ()
 
 AdvectionDiffusionReactionEquations (const AdvectionDiffusionReactionEquations &dummy)=delete
 Broken copy constructor. More...
 
void operator= (const AdvectionDiffusionReactionEquations &)=delete
 Broken assignment operator. More...
 
double dc_dt_adv_diff_react (const unsigned &n, const unsigned &r) const
 
void dc_dt_adv_diff_react (const unsigned &n, Vector< double > &dc_dt) const
 
void disable_ALE ()
 
void enable_ALE ()
 
void compute_norm (double &norm)
 Compute norm of the solution: sum of squares of L2 norms for reagents. More...
 
void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 Get error against and norm of exact solution. More...
 
void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 Dummy, time dependent error checker. More...
 
AdvectionDiffusionReactionSourceFctPtsource_fct_pt ()
 Access function: Pointer to source function. More...
 
AdvectionDiffusionReactionSourceFctPt source_fct_pt () const
 Access function: Pointer to source function. Const version. More...
 
AdvectionDiffusionReactionWindFctPtwind_fct_pt ()
 Access function: Pointer to wind function. More...
 
AdvectionDiffusionReactionWindFctPt wind_fct_pt () const
 Access function: Pointer to reaction function. Const version. More...
 
AdvectionDiffusionReactionReactionFctPtreaction_fct_pt ()
 Access function: Pointer to reaction function. More...
 
AdvectionDiffusionReactionReactionFctPt reaction_fct_pt () const
 Access function: Pointer to reaction function. Const version. More...
 
AdvectionDiffusionReactionReactionDerivFctPtreaction_deriv_fct_pt ()
 Access function: Pointer to reaction derivatives function. More...
 
AdvectionDiffusionReactionReactionDerivFctPt reaction_deriv_fct_pt () const
 Access function: Pointer to reaction function. Const version. More...
 
const Vector< double > & diff () const
 Vector of diffusion coefficients. More...
 
Vector< double > *& diff_pt ()
 Pointer to vector of diffusion coefficients. More...
 
const Vector< double > & tau () const
 Vector of dimensionless timescales. More...
 
Vector< double > *& tau_pt ()
 Pointer to vector of dimensionless timescales. More...
 
virtual void get_source_adv_diff_react (const unsigned &ipt, const Vector< double > &x, Vector< double > &source) const
 
void get_flux (const Vector< double > &s, Vector< double > &flux) const
 Get flux: \(\mbox{flux}[DIM r + i] = \mbox{d}C_{r} / \mbox{d}x_i \). More...
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Add the element's contribution to its residual vector (wrapper) 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)
 
double interpolated_c_adv_diff_react (const Vector< double > &s, const unsigned &i) const
 Return FE representation of function value c_i(s) at local coordinate s. More...
 
unsigned self_test ()
 Self-test: Return 0 for OK. More...
 
void integrate_reagents (Vector< double > &C) const
 Return the integrated reagent concentrations. 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)
 
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 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 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_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 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, 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 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)
 
- 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::RefineableAdvectionDiffusionReactionEquations< NREAGENT, DIM >
 RefineableAdvectionDiffusionReactionEquations ()
 Empty Constructor. More...
 
 RefineableAdvectionDiffusionReactionEquations (const RefineableAdvectionDiffusionReactionEquations< NREAGENT, DIM > &dummy)=delete
 Broken copy constructor. More...
 
void operator= (const RefineableAdvectionDiffusionReactionEquations< NREAGENT, DIM > &)=delete
 Broken assignment operator. More...
 
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 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_build ()
 
- 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...
 
Treetree_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 Nodeinterpolating_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 Nodeget_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 RefineableElementroot_element_pt ()
 
virtual RefineableElementfather_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)
 
unsigned nshape_controlling_nodes ()
 
std::map< Node *, unsignedshape_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 void compute_exact_Z2_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_flux_pt, double &error, double &norm)
 
virtual double geometric_jacobian (const Vector< double > &x)
 
- Public Member Functions inherited from oomph::RefineableQElement< DIM >
 RefineableQElement ()
 Empty constuctor. More...
 
- Public Member Functions inherited from oomph::RefineableQCrouzeixRaviartElement< DIM >
 RefineableQCrouzeixRaviartElement ()
 Constructor. More...
 
 RefineableQCrouzeixRaviartElement (const RefineableQCrouzeixRaviartElement< DIM > &dummy)=delete
 Broken copy constructor. More...
 
unsigned ncont_interpolated_values () const
 Broken assignment operator. More...
 
void rebuild_from_sons (Mesh *&mesh_pt)
 
unsigned nrecovery_order ()
 
unsigned nvertex_node () const
 Number of vertex nodes in the element. More...
 
Nodevertex_node_pt (const unsigned &j) const
 Pointer to the j-th vertex node in the element. 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 build_fp_press_adv_diff_robin_bc_element (const unsigned &face_index)
 
void identify_load_data (std::set< std::pair< Data *, unsigned >> &paired_load_data)
 
void rebuild_from_sons (Mesh *&mesh_pt)
 2D Rebuild from sons: Reconstruct pressure from the (merged) sons More...
 
void rebuild_from_sons (Mesh *&mesh_pt)
 3D Rebuild from sons: Reconstruct pressure from the (merged) sons More...
 
void further_build ()
 
void further_build ()
 
- Public Member Functions inherited from oomph::QCrouzeixRaviartElement< DIM >
 QCrouzeixRaviartElement ()
 Constructor, there are DIM+1 internal values (for the pressure) More...
 
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...
 
double p_nst (const unsigned &i) const
 
double p_nst (const unsigned &t, const unsigned &i) const
 
unsigned npres_nst () const
 Return number of pressure values. More...
 
double dpshape_and_dptest_eulerian_nst (const Vector< double > &s, Shape &ppsi, DShape &dppsidx, Shape &ptest, DShape &dptestdx) const
 
int p_local_eqn (const unsigned &n) const
 Return the local equation numbers for the 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...
 
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
 
void pshape_nst (const Vector< double > &s, Shape &psi) const
 
double dpshape_and_dptest_eulerian_nst (const Vector< double > &s, Shape &ppsi, DShape &dppsidx, Shape &ptest, DShape &dptestdx) const
 
void pshape_nst (const Vector< double > &s, Shape &psi) const
 
double dpshape_and_dptest_eulerian_nst (const Vector< double > &s, Shape &ppsi, DShape &dppsidx, Shape &ptest, DShape &dptestdx) const
 
- Public Member Functions inherited from oomph::NavierStokesEquations< DIM >
 NavierStokesEquations ()
 
const doublere () const
 Reynolds number. More...
 
const doublere_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 doubleviscosity_ratio () const
 
double *& viscosity_ratio_pt ()
 Pointer to Viscosity Ratio. More...
 
const doubledensity_ratio () const
 
double *& density_ratio_pt ()
 Pointer to Density ratio. More...
 
const doublere_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...
 
NavierStokesBodyForceFctPtbody_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...
 
NavierStokesSourceFctPtsource_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...
 
NavierStokesPressureAdvDiffSourceFctPtsource_fct_for_pressure_adv_diff ()
 
NavierStokesPressureAdvDiffSourceFctPt source_fct_for_pressure_adv_diff () const
 
intpinned_fp_pressure_eqn ()
 
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...
 
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...
 
void get_vorticity (const Vector< double > &s, double &vorticity) const
 Compute the scalar vorticity at local coordinate s (2D) 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 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
 
void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt) 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_norm (double &norm)
 Compute norm of solution: square of the L2 norm of the velocities. More...
 
void compute_norm (Vector< double > &norm)
 Compute the vector norm of the FEM solution. More...
 
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 compute_error (FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 
void compute_error (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 fill_in_pressure_advection_diffusion_residuals (Vector< double > &residuals)
 
void fill_in_pressure_advection_diffusion_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
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 output_pressure_advection_diffusion_robin_elements (std::ostream &outfile)
 
void delete_pressure_advection_diffusion_robin_elements ()
 
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...
 
double interpolated_dudx_nst (const Vector< double > &s, const unsigned &i, const unsigned &j) const
 
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, 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::TemplateFreeNavierStokesEquationsBase
 TemplateFreeNavierStokesEquationsBase ()
 Constructor (empty) More...
 
virtual ~TemplateFreeNavierStokesEquationsBase ()
 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::RefineableNavierStokesEquations< DIM >
 RefineableNavierStokesEquations ()
 Constructor. More...
 
virtual Nodepressure_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)
 

Private Member Functions

double Default_Physical_Constant_Value
 Set the default physical value to be one. More...
 

Private Attributes

doubleKm_pt
 Pointer to private data. The value of Km. More...
 
doubleN_pt
 

Static Private Attributes

static double Default_Physical_Constant_Value
 The static default value of the Rayleigh number. More...
 

Additional Inherited Members

- Public Types inherited from oomph::AdvectionDiffusionReactionEquations< NREAGENT, DIM >
typedef void(* AdvectionDiffusionReactionSourceFctPt) (const Vector< double > &x, Vector< double > &f)
 
typedef void(* AdvectionDiffusionReactionReactionFctPt) (const Vector< double > &c, Vector< double > &R)
 Function pointer to reaction terms. More...
 
typedef void(* AdvectionDiffusionReactionReactionDerivFctPt) (const Vector< double > &c, DenseMatrix< double > &dRdC)
 Function pointer to derivative of reaction terms. More...
 
typedef void(* AdvectionDiffusionReactionWindFctPt) (const double &time, const Vector< double > &x, Vector< double > &wind)
 
- 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 Types inherited from oomph::NavierStokesEquations< 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)
 
- Static Public Member Functions inherited from oomph::RefineableElement
static doublemax_integrity_tolerance ()
 Max. allowed discrepancy in element integrity check. More...
 
- Static Public Member Functions inherited from oomph::RefineableNavierStokesEquations< 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 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
 
- Static Public Attributes inherited from oomph::NavierStokesEquations< DIM >
static Vector< doubleGamma
 Vector to decide whether the stress-divergence form is used or not. More...
 
- Protected Member Functions inherited from oomph::QAdvectionDiffusionReactionElement< NREAGENT, DIM, NNODE_1D >
double dshape_and_dtest_eulerian_adv_diff_react (const Vector< double > &s, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
double dshape_and_dtest_eulerian_at_knot_adv_diff_react (const unsigned &ipt, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
- 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)
 
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::RefineableAdvectionDiffusionReactionEquations< NREAGENT, DIM >
void fill_in_generic_residual_contribution_adv_diff_react (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix, unsigned flag)
 
- 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)
 
- Protected Member Functions inherited from oomph::QCrouzeixRaviartElement< DIM >
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_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 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
 
- Protected Member Functions inherited from oomph::NavierStokesEquations< 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::RefineableNavierStokesEquations< 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)
 
void fill_in_generic_pressure_advection_diffusion_contribution_nst (Vector< double > &residuals, DenseMatrix< double > &jacobian, unsigned flag)
 
virtual void get_dresidual_dnodal_coordinates (RankThreeTensor< double > &dresidual_dnodal_coordinates)
 
- 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::AdvectionDiffusionReactionEquations< NREAGENT, DIM >
Vector< double > * Diff_pt
 Pointer to global diffusion coefficients. More...
 
Vector< double > * Tau_pt
 Pointer to global timescales. More...
 
AdvectionDiffusionReactionSourceFctPt Source_fct_pt
 Pointer to source function: More...
 
AdvectionDiffusionReactionWindFctPt Wind_fct_pt
 Pointer to wind function: More...
 
AdvectionDiffusionReactionReactionFctPt Reaction_fct_pt
 Pointer to reaction function. More...
 
AdvectionDiffusionReactionReactionDerivFctPt Reaction_deriv_fct_pt
 Pointer to reaction derivatives. More...
 
bool ALE_is_disabled
 
- 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::RefineableElement
TreeTree_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...
 
- Protected Attributes inherited from oomph::QCrouzeixRaviartElement< DIM >
unsigned P_nst_internal_index
 
- Protected Attributes inherited from oomph::NavierStokesEquations< DIM >
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...
 
doubleReInvFr_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...
 
NavierStokesPressureAdvDiffSourceFctPt Press_adv_diff_source_fct_pt
 
bool ALE_is_disabled
 
Vector< FpPressureAdvDiffRobinBCElementBase * > Pressure_advection_diffusion_robin_element_pt
 
int Pinned_fp_pressure_eqn
 
- Static Protected Attributes inherited from oomph::AdvectionDiffusionReactionEquations< NREAGENT, DIM >
static const unsigned N_reagent
 
- 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
 
- 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::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >

A RefinebleELement class that solves the coupled Navier–Stokes and energy equations by coupling two pre-existing classes. The QAdvectionDiffusionReactionElement with bi-quadratic interpolation for the scalar variables (temperature and concentration) and QCrouzeixRaviartElement which solves the Navier–Stokes equations using bi-quadratic interpolation for the velocities and a discontinuous bi-linear interpolation for the pressure. Note that we are free to choose the order in which we store the variables at the nodes. In this case we choose to store the variables in the order fluid velocities followed by bulk concentration followed by micelle concentration. We must, therefore, overload the function AdvectionDiffusionReactionEquations<2,DIM>::u_index_adv_diff_react() to indicate that the concentraion is stored at the DIM-th position not the 0-th. We do not need to overload the corresponding function in the NavierStokesEquations<DIM> class because the velocities are stored first.

Constructor & Destructor Documentation

◆ RefineableDoubleBuoyantQCrouzeixRaviartElement()

Constructor: call the underlying constructors and initialise the pointer to the Rayleigh number to point to the default value of 0.0.

603  :
604  RefineableQAdvectionDiffusionReactionElement<2,DIM,3>(),
605  RefineableQCrouzeixRaviartElement<DIM>()
606  {
609  }
static double Default_Physical_Constant_Value
The static default value of the Rayleigh number.
Definition: two_layer_soluble_surfactant/double_buoyant_navier_stokes_elements.h:596
double * N_pt
Definition: two_layer_soluble_surfactant/double_buoyant_navier_stokes_elements.h:593
double * Km_pt
Pointer to private data. The value of Km.
Definition: two_layer_soluble_surfactant/double_buoyant_navier_stokes_elements.h:590

References oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::Default_Physical_Constant_Value, oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::Km_pt, and oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::N_pt.

Member Function Documentation

◆ c_index_adv_diff_react()

template<unsigned DIM>
unsigned oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::c_index_adv_diff_react ( const unsigned i) const
inlinevirtual

Overload the index at which the temperature and solute concentration variables are stored.

Reimplemented from oomph::AdvectionDiffusionReactionEquations< NREAGENT, DIM >.

763  {return DIM+i;}
int i
Definition: BiCGSTAB_step_by_step.cpp:9
#define DIM
Definition: linearised_navier_stokes_elements.h:44

References DIM, and i.

Referenced by oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::fill_in_off_diagonal_jacobian_blocks_by_fd(), and oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::integrated_C_and_M().

◆ compute_error() [1/2]

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::compute_error ( std::ostream &  outfile,
FiniteElement::SteadyExactSolutionFctPt  exact_soln_pt,
double error,
double norm 
)
inlinevirtual

Validate against exact solution. Solution is provided via function pointer. Plot at a given number of plot points and compute L2 error and L2 norm of velocity solution over element Call the broken default

Reimplemented from oomph::FiniteElement.

972  {FiniteElement::compute_error(outfile,exact_soln_pt,error,norm);}
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.
Definition: elements.h:3198
int error
Definition: calibrate.py:297

References oomph::FiniteElement::compute_error(), and calibrate::error.

◆ compute_error() [2/2]

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::compute_error ( std::ostream &  outfile,
FiniteElement::UnsteadyExactSolutionFctPt  exact_soln_pt,
const double time,
double error,
double norm 
)
inlinevirtual

Validate against exact solution at given time Solution is provided via function pointer. Plot at a given number of plot points and compute L2 error and L2 norm of velocity solution over element Call the broken default

Reimplemented from oomph::FiniteElement.

961  {FiniteElement::compute_error(outfile,exact_soln_pt,
962  time,error,norm);}

References oomph::FiniteElement::compute_error(), and calibrate::error.

◆ compute_norm()

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::compute_norm ( double norm)
inlinevirtual

Compute norm of solution – broken virtual can be overloaded by element writer to implement whatever norm is desired for the specific element

Reimplemented from oomph::GeneralisedElement.

void compute_norm(double &norm)
Compute norm of solution: square of the L2 norm of the velocities.
Definition: navier_stokes_elements.cc:186

References oomph::NavierStokesEquations< DIM >::compute_norm().

◆ Default_Physical_Constant_Value()

double oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< 2 >::Default_Physical_Constant_Value
private

Set the default physical value to be one.

◆ disable_ALE()

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::disable_ALE ( )
inlinevirtual

Final override for disable ALE.

Reimplemented from oomph::FiniteElement.

641  {
642  //Disable ALE in both sets of equations
645  }
void disable_ALE()
Definition: advection_diffusion_reaction_elements.h:180
void disable_ALE()
Definition: navier_stokes_elements.h:909

References oomph::AdvectionDiffusionReactionEquations< NREAGENT, DIM >::disable_ALE(), and oomph::NavierStokesEquations< DIM >::disable_ALE().

◆ enable_ALE()

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::enable_ALE ( )
inlinevirtual

Final override for enable ALE.

Reimplemented from oomph::FiniteElement.

649  {
650  //Enable ALE in both sets of equations
653  }
void enable_ALE()
Definition: advection_diffusion_reaction_elements.h:189
void enable_ALE()
Definition: navier_stokes_elements.h:918

References oomph::AdvectionDiffusionReactionEquations< NREAGENT, DIM >::enable_ALE(), and oomph::NavierStokesEquations< DIM >::enable_ALE().

◆ fill_in_contribution_to_jacobian()

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::fill_in_contribution_to_jacobian ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian 
)
inlinevirtual

Compute the element's residual Vector and the Jacobian matrix. Use finite-differencing only for the off-diagonal blocks.

Reimplemented from oomph::GeneralisedElement.

1279  {
1280  //Calculate the Navier-Stokes contributions (diagonal block and residuals)
1282  fill_in_contribution_to_jacobian(residuals,jacobian);
1283 
1284  //Calculate the advection-diffusion contributions
1285  //(diagonal block and residuals)
1287  fill_in_contribution_to_jacobian(residuals,jacobian);
1288 
1289  //Add in the off-diagonal blocks
1290  fill_in_off_diagonal_jacobian_blocks_by_fd(residuals,jacobian);
1291 
1292  } //End of jacobian calculation
void fill_in_contribution_to_jacobian(Vector< double > &residuals, DenseMatrix< double > &jacobian)
Definition: advection_diffusion_reaction_elements.h:533
void fill_in_contribution_to_jacobian(Vector< double > &residuals, DenseMatrix< double > &jacobian)
Definition: navier_stokes_elements.h:1273
void fill_in_off_diagonal_jacobian_blocks_by_fd(Vector< double > &residuals, DenseMatrix< double > &jacobian)
Definition: two_layer_soluble_surfactant/double_buoyant_navier_stokes_elements.h:1019

References oomph::AdvectionDiffusionReactionEquations< NREAGENT, DIM >::fill_in_contribution_to_jacobian(), oomph::NavierStokesEquations< DIM >::fill_in_contribution_to_jacobian(), and oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::fill_in_off_diagonal_jacobian_blocks_by_fd().

◆ fill_in_contribution_to_jacobian_and_mass_matrix()

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::fill_in_contribution_to_jacobian_and_mass_matrix ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian,
DenseMatrix< double > &  mass_matrix 
)
inlinevirtual

Add the element's contribution to its residuals vector, jacobian matrix and mass matrix

Reimplemented from oomph::GeneralisedElement.

1302  {
1305  residuals,jacobian,mass_matrix);
1306 
1309  residuals,jacobian,mass_matrix);
1310 
1311  fill_in_off_diagonal_jacobian_blocks_by_fd(residuals,jacobian);
1312  }
void fill_in_contribution_to_jacobian_and_mass_matrix(Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
Definition: advection_diffusion_reaction_elements.h:544
void fill_in_contribution_to_jacobian_and_mass_matrix(Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
Definition: navier_stokes_elements.h:1283

References oomph::AdvectionDiffusionReactionEquations< NREAGENT, DIM >::fill_in_contribution_to_jacobian_and_mass_matrix(), oomph::NavierStokesEquations< DIM >::fill_in_contribution_to_jacobian_and_mass_matrix(), and oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::fill_in_off_diagonal_jacobian_blocks_by_fd().

◆ fill_in_contribution_to_residuals()

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::fill_in_contribution_to_residuals ( Vector< double > &  residuals)
inlinevirtual

Calculate the element's contribution to the residual vector. Recall that fill_in_* functions MUST NOT initialise the entries in the vector to zero. This allows us to call the fill_in_* functions of the constituent single-physics elements sequentially, without wiping out any previously computed entries.

Reimplemented from oomph::GeneralisedElement.

992  {
993  //Fill in the residuals of the Navier-Stokes equations
996 
997  //Fill in the residuals of the advection-diffusion eqautions
1000  }
void fill_in_contribution_to_residuals(Vector< double > &residuals)
Add the element's contribution to its residual vector (wrapper)
Definition: advection_diffusion_reaction_elements.h:519
void fill_in_contribution_to_residuals(Vector< double > &residuals)
Compute the element's residual Vector.
Definition: navier_stokes_elements.h:1260

References oomph::AdvectionDiffusionReactionEquations< NREAGENT, DIM >::fill_in_contribution_to_residuals(), and oomph::NavierStokesEquations< DIM >::fill_in_contribution_to_residuals().

◆ fill_in_off_diagonal_jacobian_blocks_by_fd()

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::fill_in_off_diagonal_jacobian_blocks_by_fd ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian 
)
inline

Helper function to get the off-diagonal blocks of the Jacobian matrix by finite differences

1021  {
1022  //Local storage for the index in the nodes at which the
1023  //Navier-Stokes velocities are stored (we know that this should be 0,1,2)
1024  unsigned u_nodal_nst[DIM];
1025  for(unsigned i=0;i<DIM;i++)
1026  {u_nodal_nst[i] = this->u_index_nst(i);}
1027 
1028  //Local storage for the index at which the temperature and
1029  //solute are stored
1030  unsigned C_nodal_adv_diff_react[2];
1031  for(unsigned r=0;r<2;r++)
1032  {C_nodal_adv_diff_react[r] = this->c_index_adv_diff_react(r);}
1033 
1034  //Find the total number of unknowns in the elements
1035  unsigned n_dof = this->ndof();
1036 
1037  //Temporary storage for residuals
1038  Vector<double> newres(n_dof);
1039 
1040  //Storage for local unknown and local equation
1041  int local_unknown =0, local_eqn = 0;
1042 
1043  //Set the finite difference step
1044  double fd_step = FiniteElement::Default_fd_jacobian_step;
1045 
1046  //Find the number of nodes
1047  unsigned n_node = this->nnode();
1048 
1049  //Calculate the contribution of the Navier--Stokes velocities to the
1050  //advection-diffusion equations
1051 
1052  //Loop over the nodes
1053  for(unsigned n=0;n<n_node;n++)
1054  {
1055  //Cache a local pointer to the node
1056  Node* const local_node_pt = this->node_pt(n);
1057 
1058  //There are DIM values of the velocities
1059  for(unsigned i=0;i<DIM;i++)
1060  {
1061  //Need to check for hanging nodes (if not hanging do the usual)
1062  if(local_node_pt->is_hanging(u_nodal_nst[i]) == false)
1063  {
1064  //Get the local velocity equation number
1065  local_unknown = this->nodal_local_eqn(n,u_nodal_nst[i]);
1066 
1067  //If it's not pinned
1068  if(local_unknown >= 0)
1069  {
1070  //Get a pointer to the velocity value
1071  double *value_pt = local_node_pt->value_pt(u_nodal_nst[i]);
1072 
1073  //Save the old value
1074  double old_var = *value_pt;
1075 
1076  //Increment the value
1077  *value_pt += fd_step;
1078 
1079  //Get the altered advection-diffusion residuals.
1080  //Do this using fill_in because get_residuals has never been
1081  //overloaded, and will actually compute all residuals which
1082  //is slightly inefficient.
1083  for(unsigned l=0;l<n_dof;l++) {newres[l] = 0.0;}
1086 
1087  //Now fill in the Advection-Diffusion-Reaction sub-block
1088  //of the jacobian
1089  for(unsigned l=0;l<n_node;l++)
1090  {
1091  for(unsigned r=0;r<2;r++)
1092  {
1093  //Local equation for temperature or solute
1094  local_eqn = this->nodal_local_eqn(l,C_nodal_adv_diff_react[r]);
1095 
1096  //If it's not a boundary condition
1097  if(local_eqn >= 0)
1098  {
1099  double sum = (newres[local_eqn] - residuals[local_eqn])/fd_step;
1100  jacobian(local_eqn,local_unknown) = sum;
1101  }
1102  }
1103  }
1104 
1105  //Reset the nodal data
1106  *value_pt = old_var;
1107  }
1108  }
1109  //Otherwise the value is hanging
1110  else
1111  {
1112  //Get the local hanging infor
1113  HangInfo *hang_info_pt = local_node_pt->hanging_pt(u_nodal_nst[i]);
1114  //Loop over the master nodes
1115  const unsigned n_master = hang_info_pt->nmaster();
1116  for(unsigned m=0;m<n_master;m++)
1117  {
1118  //Get the pointer to the master node
1119  Node* const master_node_pt = hang_info_pt->master_node_pt(m);
1120 
1121  //Get the number of the unknown
1122  local_unknown = this->local_hang_eqn(master_node_pt,u_nodal_nst[i]);
1123  //If the variable is free
1124  if(local_unknown >= 0)
1125  {
1126  //Get a pointer to the nodal value stored at the master node
1127  double* const value_pt = master_node_pt->value_pt(u_nodal_nst[i]);
1128  //Save the old value
1129  double old_var = *value_pt;
1130 
1131  //Increment the value
1132  *value_pt += fd_step;
1133 
1134  //Get the altered advection-diffusion residuals.
1135  //Do this using fill_in because get_residuals has never been
1136  //overloaded, and will actually compute all residuals which
1137  //is slightly inefficient.
1138  for(unsigned l=0;l<n_dof;l++) {newres[l] = 0.0;}
1141 
1142  //Now fill in the Advection-Diffusion-Reaction sub-block
1143  //of the jacobian
1144  for(unsigned l=0;l<n_node;l++)
1145  {
1146  for(unsigned r=0;r<2;r++)
1147  {
1148  //Local equation for temperature or solute
1149  local_eqn = this->nodal_local_eqn(l,C_nodal_adv_diff_react[r]);
1150 
1151  //If it's not a boundary condition
1152  if(local_eqn >= 0)
1153  {
1154  double sum = (newres[local_eqn] - residuals[local_eqn])/fd_step;
1155  jacobian(local_eqn,local_unknown) = sum;
1156  }
1157  }
1158  }
1159 
1160  //Reset the nodal data
1161  *value_pt = old_var;
1162  }
1163  } //End of loop over master nodes
1164  } //End of hanging case
1165  }
1166 
1167  //Calculate the contribution of the temperature to the Navier--Stokes
1168  //equations
1169  for(unsigned r=0;r<2;r++)
1170  {
1171  //Need to check for hanging nodes (if not hanging do the usual)
1172  if(local_node_pt->is_hanging(C_nodal_adv_diff_react[r]) == false)
1173  {
1174  //Get the local equation number
1175  local_unknown = this->nodal_local_eqn(n,C_nodal_adv_diff_react[r]);
1176 
1177  //If it's not pinned
1178  if(local_unknown >= 0)
1179  {
1180  //Get a pointer to the concentration value
1181  double *value_pt =
1182  local_node_pt->value_pt(C_nodal_adv_diff_react[r]);
1183 
1184  //Save the old value
1185  double old_var = *value_pt;
1186 
1187  //Increment the value (Really need access)
1188  *value_pt += fd_step;
1189 
1190  //Get the altered Navier--Stokes residuals
1191  //Do this using fill_in because get_residuals has never been
1192  //overloaded, and will actually compute all residuals which
1193  //is slightly inefficient.
1194  for(unsigned l=0;l<n_dof;l++) {newres[l] = 0.0;}
1196 
1197  //Now fill in the Navier-Stokes sub-block
1198  for(unsigned l=0;l<n_node;l++)
1199  {
1200  //Loop over the fluid velocities
1201  for(unsigned j=0;j<DIM;j++)
1202  {
1203  //Local fluid equation
1204  local_eqn = this->nodal_local_eqn(l,u_nodal_nst[j]);
1205  if(local_eqn >= 0)
1206  {
1207  double sum = (newres[local_eqn] - residuals[local_eqn])/fd_step;
1208  jacobian(local_eqn,local_unknown) = sum;
1209  }
1210  }
1211  }
1212 
1213  //Reset the nodal data
1214  *value_pt = old_var;
1215  }
1216  } //End of not hanging case
1217  else
1218  {
1219  //Get the local hanging infor
1220  HangInfo *hang_info_pt = local_node_pt->hanging_pt(C_nodal_adv_diff_react[r]);
1221  //Loop over the master nodes
1222  const unsigned n_master = hang_info_pt->nmaster();
1223  for(unsigned m=0;m<n_master;m++)
1224  {
1225  //Get the pointer to the master node
1226  Node* const master_node_pt = hang_info_pt->master_node_pt(m);
1227 
1228  //Get the number of the unknown
1229  local_unknown = this->local_hang_eqn(master_node_pt,C_nodal_adv_diff_react[r]);
1230  //If the variable is free
1231  if(local_unknown >= 0)
1232  {
1233  //Get a pointer to the nodal value stored at the master node
1234  double* const value_pt = master_node_pt->value_pt(C_nodal_adv_diff_react[r]);
1235 
1236  //Save the old value
1237  double old_var = *value_pt;
1238 
1239  //Increment the value (Really need access)
1240  *value_pt += fd_step;
1241 
1242  //Get the altered Navier--Stokes residuals
1243  //Do this using fill_in because get_residuals has never been
1244  //overloaded, and will actually compute all residuals which
1245  //is slightly inefficient.
1246  for(unsigned l=0;l<n_dof;l++) {newres[l] = 0.0;}
1248 
1249  //Now fill in the Navier-Stokes sub-block
1250  for(unsigned l=0;l<n_node;l++)
1251  {
1252  //Loop over the fluid velocities
1253  for(unsigned j=0;j<DIM;j++)
1254  {
1255  //Local fluid equation
1256  local_eqn = this->nodal_local_eqn(l,u_nodal_nst[j]);
1257  if(local_eqn >= 0)
1258  {
1259  double sum = (newres[local_eqn] - residuals[local_eqn])/fd_step;
1260  jacobian(local_eqn,local_unknown) = sum;
1261  }
1262  }
1263  }
1264 
1265  //Reset the nodal data
1266  *value_pt = old_var;
1267  }
1268  } //End of loop over master nodes
1269  } //End of hanging case
1270  } //End of loop over reagents
1271 
1272  }//End of loop over nodes
1273  }
Node *& node_pt(const unsigned &n)
Return a pointer to the local node n.
Definition: elements.h:2175
int nodal_local_eqn(const unsigned &n, const unsigned &i) const
Definition: elements.h:1432
unsigned nnode() const
Return the number of nodes.
Definition: elements.h:2210
static double Default_fd_jacobian_step
Definition: elements.h:1198
unsigned ndof() const
Return the number of equations/dofs in the element.
Definition: elements.h:835
virtual unsigned u_index_nst(const unsigned &i) const
Definition: navier_stokes_elements.h:866
const double & n() const
Access function for the number of monomers in the micelle.
Definition: two_layer_soluble_surfactant/double_buoyant_navier_stokes_elements.h:634
unsigned c_index_adv_diff_react(const unsigned &i) const
Definition: two_layer_soluble_surfactant/double_buoyant_navier_stokes_elements.h:762
int local_hang_eqn(Node *const &node_pt, const unsigned &i)
Definition: refineable_elements.h:278
int * m
Definition: level2_cplx_impl.h:294
r
Definition: UniformPSDSelfTest.py:20
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2

References oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::c_index_adv_diff_react(), oomph::GeneralisedElement::Default_fd_jacobian_step, DIM, oomph::AdvectionDiffusionReactionEquations< NREAGENT, DIM >::fill_in_contribution_to_residuals(), oomph::NavierStokesEquations< DIM >::fill_in_contribution_to_residuals(), oomph::Node::hanging_pt(), i, oomph::Node::is_hanging(), j, oomph::RefineableElement::local_hang_eqn(), m, oomph::HangInfo::master_node_pt(), oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::n(), oomph::GeneralisedElement::ndof(), oomph::HangInfo::nmaster(), oomph::FiniteElement::nnode(), oomph::FiniteElement::nodal_local_eqn(), oomph::FiniteElement::node_pt(), UniformPSDSelfTest::r, oomph::NavierStokesEquations< DIM >::u_index_nst(), and oomph::Data::value_pt().

Referenced by oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::fill_in_contribution_to_jacobian(), and oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::fill_in_contribution_to_jacobian_and_mass_matrix().

◆ further_build()

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::further_build ( )
inlinevirtual

Call the underlying single-physics element's further_build() functions and make sure that the pointer to the Rayleigh number is passed to the sons

Reimplemented from oomph::RefineableElement.

866  {
869 
870  //Cast the pointer to the father element to the specific
871  //element type
872  RefineableDoubleBuoyantQCrouzeixRaviartElement<DIM>* cast_father_element_pt
873  = dynamic_cast<RefineableDoubleBuoyantQCrouzeixRaviartElement<DIM>*>(
874  this->father_element_pt());
875 
876  //Set the pointer to the physical variables to be the same as
877  //the father
878  this->Km_pt = cast_father_element_pt->km_pt();
879  this->N_pt = cast_father_element_pt->n_pt();
880  } //end of further build
void further_build()
Definition: refineable_advection_diffusion_reaction_elements.h:166
virtual RefineableElement * father_element_pt() const
Return a pointer to the father element.
Definition: refineable_elements.h:539

References oomph::RefineableElement::father_element_pt(), oomph::RefineableAdvectionDiffusionReactionEquations< NREAGENT, DIM >::further_build(), oomph::RefineableQCrouzeixRaviartElement< DIM >::further_build(), oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::km_pt(), and oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::n_pt().

◆ further_setup_hanging_nodes()

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

The additional hanging node information must be set up for both single-physics elements.

Reimplemented from oomph::RefineableElement.

843  {
847  }
void further_setup_hanging_nodes()
Definition: refineable_advection_diffusion_reaction_elements.h:268
void further_setup_hanging_nodes()
Definition: refineable_navier_stokes_elements.h:1234

References oomph::RefineableQAdvectionDiffusionReactionElement< NREAGENT, DIM, NNODE_1D >::further_setup_hanging_nodes(), and oomph::RefineableQCrouzeixRaviartElement< DIM >::further_setup_hanging_nodes().

◆ get_interpolated_values() [1/2]

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

Get all continuously interpolated values at the local coordinate s at time level t (t=0: present; t>0: previous). We choose to put the fluid velocities first, followed by the temperature

Implements oomph::RefineableElement.

816  {
817  //Storage for the fluid velocities
818  Vector<double> nst_values;
819 
820  //Get the fluid velocities from the fluid element
822  get_interpolated_values(t,s,nst_values);
823 
824  //Storage for the temperature
825  Vector<double> advection_values;
826 
827  //Get the temperature from the advection-diffusion element
829  get_interpolated_values(s,advection_values);
830 
831  //Add the fluid velocities to the values vector
832  for(unsigned i=0;i<DIM;i++) {values.push_back(nst_values[i]);}
833 
834  //Add the concentration to the end
835  for(unsigned i=0;i<DIM;i++) {values.push_back(advection_values[i]);}
836 
837  } // end of get_interpolated_values
void get_interpolated_values(const Vector< double > &s, Vector< double > &values)
Definition: refineable_advection_diffusion_reaction_elements.h:97
void get_interpolated_values(const Vector< double > &s, Vector< double > &values)
Definition: refineable_navier_stokes_elements.h:1178
RealScalar s
Definition: level1_cplx_impl.h:130
t
Definition: plotPSD.py:36

References DIM, oomph::RefineableAdvectionDiffusionReactionEquations< NREAGENT, DIM >::get_interpolated_values(), oomph::RefineableQCrouzeixRaviartElement< DIM >::get_interpolated_values(), i, s, and plotPSD::t.

◆ get_interpolated_values() [2/2]

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

Get the continuously interpolated values at the local coordinate s. We choose to put the fluid velocities first, followed by the temperature.

Reimplemented from oomph::RefineableElement.

786  {
787  //Storage for the fluid velocities
788  Vector<double> nst_values;
789 
790  //Get the fluid velocities from the fluid element
792  get_interpolated_values(s,nst_values);
793 
794  //Storage for the temperature
795  Vector<double> advection_values;
796 
797  //Get the temperature from the advection-diffusion element
799  get_interpolated_values(s,advection_values);
800 
801  //Add the fluid velocities to the values vector
802  for(unsigned i=0;i<DIM;i++) {values.push_back(nst_values[i]);}
803 
804  //Add the concentration values to the end
805  for(unsigned i=0;i<2;i++) {values.push_back(advection_values[i]);}
806  }

References DIM, oomph::RefineableAdvectionDiffusionReactionEquations< NREAGENT, DIM >::get_interpolated_values(), oomph::RefineableQCrouzeixRaviartElement< DIM >::get_interpolated_values(), i, and s.

◆ get_reaction_adv_diff_react()

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::get_reaction_adv_diff_react ( const unsigned ipt,
const Vector< double > &  C,
Vector< double > &  R 
) const
inlinevirtual

Overload the reaction terms to couple the concentration and micelle terms

Reimplemented from oomph::AdvectionDiffusionReactionEquations< NREAGENT, DIM >.

660  {
661  //Compute the flux between equations (equation (2.22))
662  const double J_m = this->km()*(pow(C[0],this->n()) - C[1]);
663  //Return the reaction terms
664  R[0] = J_m;
665  R[1] = -J_m;
666  }
@ R
Definition: StatisticsVector.h:21
const double & km() const
Access function for the transfer constant.
Definition: two_layer_soluble_surfactant/double_buoyant_navier_stokes_elements.h:628
EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC bfloat16 pow(const bfloat16 &a, const bfloat16 &b)
Definition: BFloat16.h:625

References oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::km(), oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::n(), Eigen::bfloat16_impl::pow(), and R.

◆ get_reaction_deriv_adv_diff_react()

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::get_reaction_deriv_adv_diff_react ( const unsigned ipt,
const Vector< double > &  C,
DenseMatrix< double > &  dRdC 
) const
inlinevirtual

Get the derivatives of the reaction terms with respect to the concentration variables. If no explicit function pointer is set, these will be calculated by finite differences

Reimplemented from oomph::AdvectionDiffusionReactionEquations< NREAGENT, DIM >.

673  {
674  const double Km = this->km(); const double N = this->n();
675  //Fill in the derivative terms
676  dRdC(0,0) = Km*N*pow(C[0],N-1); dRdC(0,1) = -Km;
677  dRdC(1,0) = -dRdC(0,0); dRdC(1,1) = -dRdC(0,1);
678  }
Definition: matrices.h:74
@ N
Definition: constructor.cpp:22

References oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::km(), N, oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::n(), and Eigen::bfloat16_impl::pow().

◆ get_wind_adv_diff_react()

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::get_wind_adv_diff_react ( const unsigned ipt,
const Vector< double > &  s,
const Vector< double > &  x,
Vector< double > &  wind 
) const
inlinevirtual

Overload the wind function in the advection-diffusion equations. This provides the coupling from the Navier–Stokes equations to the advection-diffusion equations because the wind is the fluid velocity.

Reimplemented from oomph::AdvectionDiffusionReactionEquations< NREAGENT, DIM >.

980  {
981  //The wind function is simply the velocity at the points
982  this->interpolated_u_nst(s,wind);
983  }
void interpolated_u_nst(const Vector< double > &s, Vector< double > &veloc) const
Compute vector of FE interpolated velocity u at local coordinate s.
Definition: navier_stokes_elements.h:1505

References oomph::NavierStokesEquations< DIM >::interpolated_u_nst().

◆ get_Z2_compound_flux_indices()

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::get_Z2_compound_flux_indices ( Vector< unsigned > &  flux_index)
inlinevirtual

Fill in which flux components are associated with the fluid measure and which are associated with the temperature measure

Reimplemented from oomph::ElementWithZ2ErrorEstimator.

933  {
934  //Find the number of fluid fluxes
935  unsigned n_fluid_flux =
937  //Find the number of concentration fluxes
938  unsigned n_conc_flux =
940 
941  //The fluid fluxes are first
942  //The values of the flux_index vector are zero on entry, so we
943  //could omit this line
944  for(unsigned i=0;i<n_fluid_flux;i++) {flux_index[i] = 0;}
945 
946  //Set the concentration fluxes (the last set of fluxes
947  for(unsigned i=0;i<n_conc_flux;i++) {flux_index[n_fluid_flux + i] = 1;}
948 
949  } //end of get_Z2_compound_flux_indices
unsigned num_Z2_flux_terms()
Number of 'flux' terms for Z2 error estimation.
Definition: refineable_advection_diffusion_reaction_elements.h:80
unsigned num_Z2_flux_terms()
Number of 'flux' terms for Z2 error estimation.
Definition: refineable_navier_stokes_elements.h:395

References i, oomph::RefineableAdvectionDiffusionReactionEquations< NREAGENT, DIM >::num_Z2_flux_terms(), and oomph::RefineableNavierStokesEquations< DIM >::num_Z2_flux_terms().

◆ get_Z2_flux()

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::get_Z2_flux ( const Vector< double > &  s,
Vector< double > &  flux 
)
inlinevirtual

Get the Z2 flux by concatenating the fluxes from the fluid and the advection diffusion elements.

Implements oomph::ElementWithZ2ErrorEstimator.

900  {
901  //Find the number of fluid fluxes
902  unsigned n_fluid_flux =
904 
905  //Fill in the first flux entries as the velocity entries
907 
908  //Find the number of concentration fluxes
909  unsigned n_conc_flux =
911  Vector<double> conc_flux(n_conc_flux);
912 
913  //Get the temperature flux
915  get_Z2_flux(s,conc_flux);
916 
917  //Add the concentration flux to the end of the flux vector
918  for(unsigned i=0;i<n_conc_flux;i++)
919  {
920  flux[n_fluid_flux+i] = conc_flux[i];
921  }
922 
923  } //end of get_Z2_flux
void get_Z2_flux(const Vector< double > &s, Vector< double > &flux)
Definition: refineable_advection_diffusion_reaction_elements.h:87
void get_Z2_flux(const Vector< double > &s, Vector< double > &flux)
Definition: refineable_navier_stokes_elements.h:403
void flux(const double &time, const Vector< double > &x, double &flux)
Get flux applied along boundary x=0.
Definition: pretend_melt.cc:59

References ProblemParameters::flux(), oomph::RefineableAdvectionDiffusionReactionEquations< NREAGENT, DIM >::get_Z2_flux(), oomph::RefineableNavierStokesEquations< DIM >::get_Z2_flux(), i, oomph::RefineableAdvectionDiffusionReactionEquations< NREAGENT, DIM >::num_Z2_flux_terms(), oomph::RefineableNavierStokesEquations< DIM >::num_Z2_flux_terms(), and s.

◆ integrated_C_and_M()

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::integrated_C_and_M ( double int_C,
double int_M 
)
inline
1315  {
1316  //Vector to store the integrals
1317  Vector<double> sum(2,0.0);
1318 
1319  //Find the number of nodes
1320  const unsigned n_node = this->nnode();
1321  //Storage for the shape functions
1322  Shape psi(n_node);
1323 
1324  unsigned C_index[2];
1325  for(unsigned r=0;r<2;++r)
1326  {C_index[r] = this->c_index_adv_diff_react(r);}
1327 
1328  //Loop over the integration points
1329  const unsigned n_intpt = this->integral_pt()->nweight();
1330  for(unsigned ipt=0;ipt<n_intpt;ipt++)
1331  {
1332  //Get the integral weight
1333  double w = this->integral_pt()->weight(ipt);
1334  //Get the value of the Jacobian of the mapping to global coordinates
1335  double J = this->J_eulerian_at_knot(ipt);
1336  double W = w*J;
1337  //Get the shape function at the know
1338  this->shape_at_knot(ipt,psi);
1339 
1340  //Get the interpolated values
1341  Vector<double> interpolated_C(2,0.0);
1342  for(unsigned l=0;l<n_node;++l)
1343  {
1344  const double psi_ = psi(l);
1345  for(unsigned r=0;r<2;++r)
1346  {
1347  interpolated_C[r] += this->nodal_value(l,C_index[r])*psi_;
1348  }
1349  }
1350 
1351  for(unsigned r=0;r<2;++r)
1352  {
1353  sum[r] += interpolated_C[r]*W;
1354  }
1355  } //End of integration loop
1356 
1357  //Return the values
1358  int_C = sum[0]; int_M = sum[1];
1359  }
JacobiRotation< float > J
Definition: Jacobi_makeJacobi.cpp:3
RowVector3d w
Definition: Matrix_resize_int.cpp:3
double nodal_value(const unsigned &n, const unsigned &i) const
Definition: elements.h:2593
Integral *const & integral_pt() const
Return the pointer to the integration scheme (const version)
Definition: elements.h:1963
virtual double J_eulerian_at_knot(const unsigned &ipt) const
Definition: elements.cc:4168
virtual void shape_at_knot(const unsigned &ipt, Shape &psi) const
Definition: elements.cc:3220
virtual unsigned nweight() const =0
Return the number of integration points of the scheme.
virtual double weight(const unsigned &i) const =0
Return weight of i-th integration point.
@ W
Definition: quadtree.h:63

References oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::c_index_adv_diff_react(), oomph::FiniteElement::integral_pt(), J, oomph::FiniteElement::J_eulerian_at_knot(), oomph::FiniteElement::nnode(), oomph::FiniteElement::nodal_value(), oomph::Integral::nweight(), UniformPSDSelfTest::r, oomph::FiniteElement::shape_at_knot(), w, oomph::QuadTreeNames::W, and oomph::Integral::weight().

◆ km()

◆ km_pt()

template<unsigned DIM>
double* & oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::km_pt ( )
inline

Access function for the pointer to transfer constant.

631 {return Km_pt;}

References oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::Km_pt.

Referenced by oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::further_build().

◆ n()

◆ n_pt()

template<unsigned DIM>
double* & oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::n_pt ( )
inline

Access function for the pointer to the number of monomers in the micelle.

637 {return N_pt;}

References oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::N_pt.

Referenced by oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::further_build().

◆ ncompound_fluxes()

template<unsigned DIM>
unsigned oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::ncompound_fluxes ( )
inlinevirtual

The number of compound fluxes is two (one for the fluid and one for the concentrations) ALH: May want 3 fluxes here

Reimplemented from oomph::ElementWithZ2ErrorEstimator.

928 {return 2;}

◆ ncont_interpolated_values()

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

The total number of continously interpolated values is DIM+2 (DIM fluid velocities and TWO concentrations).

Implements oomph::RefineableElement.

779  {return DIM+2;}

References DIM.

◆ nrecovery_order()

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

The recovery order is that of the NavierStokes elements.

Implements oomph::ElementWithZ2ErrorEstimator.

unsigned nrecovery_order()
Definition: refineable_navier_stokes_elements.h:1157

References oomph::RefineableQCrouzeixRaviartElement< DIM >::nrecovery_order().

◆ num_Z2_flux_terms()

template<unsigned DIM>
unsigned oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::num_Z2_flux_terms ( )
inlinevirtual

The number of Z2 flux terms is the same as that in the fluid element plus that in the advection-diffusion element

Implements oomph::ElementWithZ2ErrorEstimator.

◆ nvertex_node()

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

geometric element.

Number of vertex nodes in the element is obtained from the

Implements oomph::ElementWithZ2ErrorEstimator.

◆ output() [1/4]

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::output ( FILE *  file_pt)
inlinevirtual

C-style output function: Broken default.

Reimplemented from oomph::FiniteElement.

735  {FiniteElement::output(file_pt);}
virtual void output(std::ostream &outfile)
Definition: elements.h:3050

References oomph::FiniteElement::output().

◆ output() [2/4]

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::output ( FILE *  file_pt,
const unsigned n_plot 
)
inlinevirtual

C-style output function: Broken default.

Reimplemented from oomph::FiniteElement.

739  {FiniteElement::output(file_pt,n_plot);}

References oomph::FiniteElement::output().

◆ output() [3/4]

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::output ( std::ostream &  outfile)
inlinevirtual

Overload the standard output function with the broken default.

Reimplemented from oomph::FiniteElement.

686 {FiniteElement::output(outfile);}

References oomph::FiniteElement::output().

◆ output() [4/4]

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::output ( std::ostream &  outfile,
const unsigned nplot 
)
inlinevirtual

Output function:
Output x, y, u, v, p, theta at Nplot^DIM plot points

Reimplemented from oomph::FiniteElement.

692  {
693  //vector of local coordinates
694  Vector<double> s(DIM);
695 
696  // Tecplot header info
697  outfile << this->tecplot_zone_string(nplot);
698 
699  // Loop over plot points
700  unsigned num_plot_points=this->nplot_points(nplot);
701  for (unsigned iplot=0;iplot<num_plot_points;iplot++)
702  {
703  // Get local coordinates of plot point
704  this->get_s_plot(iplot,nplot,s);
705 
706  // Output the position of the plot point
707  for(unsigned i=0;i<DIM;i++)
708  {outfile << this->interpolated_x(s,i) << " ";}
709 
710  // Output the fluid velocities at the plot point
711  for(unsigned i=0;i<DIM;i++)
712  {outfile << this->interpolated_u_nst(s,i) << " ";}
713 
714  // Output the fluid pressure at the plot point
715  outfile << this->interpolated_p_nst(s) << " ";
716 
717  //Output the temperature and the solute concentration
718  for(unsigned i=0;i<2;i++)
719  {
720  outfile << this->interpolated_c_adv_diff_react(s,i) << " ";
721  }
722  outfile << "\n";
723  }
724  outfile << std::endl;
725 
726  // Write tecplot footer (e.g. FE connectivity lists)
727  this->write_tecplot_zone_footer(outfile,nplot);
728  } //End of output function
double interpolated_c_adv_diff_react(const Vector< double > &s, const unsigned &i) const
Return FE representation of function value c_i(s) at local coordinate s.
Definition: advection_diffusion_reaction_elements.h:556
virtual std::string tecplot_zone_string(const unsigned &nplot) const
Definition: elements.h:3161
virtual double interpolated_x(const Vector< double > &s, const unsigned &i) const
Return FE interpolated coordinate x[i] at local coordinate s.
Definition: elements.cc:3962
virtual void get_s_plot(const unsigned &i, const unsigned &nplot, Vector< double > &s, const bool &shifted_to_interior=false) const
Definition: elements.h:3148
virtual unsigned nplot_points(const unsigned &nplot) const
Definition: elements.h:3186
virtual void write_tecplot_zone_footer(std::ostream &outfile, const unsigned &nplot) const
Definition: elements.h:3174
virtual double interpolated_p_nst(const Vector< double > &s) const
Return FE interpolated pressure at local coordinate s.
Definition: navier_stokes_elements.h:1639

References DIM, oomph::FiniteElement::get_s_plot(), i, oomph::AdvectionDiffusionReactionEquations< NREAGENT, DIM >::interpolated_c_adv_diff_react(), oomph::NavierStokesEquations< DIM >::interpolated_p_nst(), oomph::NavierStokesEquations< DIM >::interpolated_u_nst(), oomph::FiniteElement::interpolated_x(), oomph::FiniteElement::nplot_points(), s, oomph::FiniteElement::tecplot_zone_string(), and oomph::FiniteElement::write_tecplot_zone_footer().

◆ output_fct() [1/2]

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::output_fct ( std::ostream &  outfile,
const unsigned Nplot,
const double time,
FiniteElement::UnsteadyExactSolutionFctPt  exact_soln_pt 
)
inlinevirtual

Output function for a time-dependent exact solution: Broken default.

Reimplemented from oomph::AdvectionDiffusionReactionEquations< NREAGENT, DIM >.

754  {
756  output_fct(outfile,Nplot,time,exact_soln_pt);
757  }
virtual void output_fct(std::ostream &outfile, const unsigned &n_plot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
Output an exact solution over the element.
Definition: elements.h:3104

References oomph::FiniteElement::output_fct().

◆ output_fct() [2/2]

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::output_fct ( std::ostream &  outfile,
const unsigned Nplot,
FiniteElement::SteadyExactSolutionFctPt  exact_soln_pt 
)
inlinevirtual

Output function for an exact solution: Broken default.

Reimplemented from oomph::FiniteElement.

745  {FiniteElement::output_fct(outfile,Nplot,exact_soln_pt);}

References oomph::FiniteElement::output_fct().

◆ rebuild_from_sons()

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

Call the rebuild_from_sons functions for each of the constituent multi-physics elements.

Implements oomph::RefineableElement.

854  {
856  rebuild_from_sons(mesh_pt);
858  }
void rebuild_from_sons(Mesh *&mesh_pt)
Rebuild from sons: empty.
Definition: refineable_advection_diffusion_reaction_elements.h:257

References oomph::RefineableQAdvectionDiffusionReactionElement< NREAGENT, DIM, NNODE_1D >::rebuild_from_sons(), and oomph::RefineableQCrouzeixRaviartElement< DIM >::rebuild_from_sons().

◆ required_nvalue()

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

The required number of values stored at the nodes is the sum of the required values of the two single-physics elements. Note that this step is generic for any multi-physics element of this type.

Reimplemented from oomph::QCrouzeixRaviartElement< DIM >.

622  {return (
unsigned required_nvalue(const unsigned &n) const
Definition: advection_diffusion_reaction_elements.h:686
virtual unsigned required_nvalue(const unsigned &n) const
Number of values (pinned or dofs) required at local node n.
Definition: navier_stokes_elements.cc:2581

References oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::n().

◆ unfix_pressure()

template<unsigned DIM>
void oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::unfix_pressure ( const unsigned p_dof)
inline

UnPin p_dof-th pressure dof.

613  {
614  this->internal_data_pt(this->P_nst_internal_index)->unpin(p_dof);
615  }
void unpin(const unsigned &i)
Unpin the i-th stored variable.
Definition: nodes.h:391
Data *& internal_data_pt(const unsigned &i)
Return a pointer to i-th internal data object.
Definition: elements.h:622
unsigned P_nst_internal_index
Definition: navier_stokes_elements.h:1757

References oomph::GeneralisedElement::internal_data_pt(), oomph::QCrouzeixRaviartElement< DIM >::P_nst_internal_index, and oomph::Data::unpin().

◆ vertex_node_pt()

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

Pointer to the j-th vertex node in the element, Call the geometric element's function.

Implements oomph::ElementWithZ2ErrorEstimator.

References j.

Member Data Documentation

◆ Default_Physical_Constant_Value

template<unsigned DIM>
double oomph::RefineableDoubleBuoyantQCrouzeixRaviartElement< DIM >::Default_Physical_Constant_Value
staticprivate

◆ Km_pt

◆ N_pt


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