Public Member Functions | Protected Member Functions | Protected Attributes | Private Attributes | List of all members
SpineGravityTractionElement< ELEMENT > Class Template Reference

#include <extra_elements.h>

+ Inheritance diagram for SpineGravityTractionElement< ELEMENT >:

Public Member Functions

 SpineGravityTractionElement (FiniteElement *const &element_pt, const int &face_index)
 
void hijack_all_nodes ()
 
void set_delta_p_pt (Data *const &delta_p_pt)
 Access function to the pointer to pressure gradient data. More...
 
const double delta_p () const
 Return the value of the pressure gradient. More...
 
unsigned nexternal_u_data ()
 Return the number of external velocity data. More...
 
const double & viscosity_ratio () const
 
double *& viscosity_ratio_pt ()
 Pointer to Viscosity Ratio. More...
 
const double & density_ratio () const
 
double *& density_ratio_pt ()
 Pointer to Density ratio. More...
 
double *& re_invfr_pt ()
 Pointer to Reynolds number divided by Froude number. More...
 
const double & re_invfr ()
 Reynolds number divided by Froude number. More...
 
const Vector< double > & g () const
 Vector of gravitational components. More...
 
Vector< double > *& g_pt ()
 Pointer to Vector of gravitational components. More...
 
double u (const unsigned &l, const unsigned &i)
 Access function for the velocity. N. B. HEAVY ASSUMPTIONS HERE. More...
 
double u (const unsigned &t, const unsigned &l, const unsigned &i) const
 
double du_dt (const unsigned &l, const unsigned &i) const
 i-th component of du/dt at local node l. More...
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Add the contribution to the residuals. More...
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 This function returns the residuals and the jacobian. More...
 
void add_generic_residual_contribution (Vector< double > &residuals, DenseMatrix< double > &jacobian, unsigned flag)
 
double get_flux ()
 Calculate the flux. More...
 
void assign_additional_local_eqn_numbers ()
 Define the local equation numbering schemes. More...
 
void output (ostream &outfile)
 Overload the output function. More...
 
void output (ostream &outfile, const unsigned &Np)
 Output function: x,y,[z],u,v,[w],p in tecplot format. More...
 
 SpineGravityTractionElement (FiniteElement *element_pt, int face_index)
 
void node_update ()
 
unsigned nexternal_u_data ()
 Return the number of external velocity data. More...
 
const double & viscosity_ratio () const
 
double *& viscosity_ratio_pt ()
 Pointer to Viscosity Ratio. More...
 
const double & density_ratio () const
 
double *& density_ratio_pt ()
 Pointer to Density ratio. More...
 
double *& re_invfr_pt ()
 Pointer to Reynolds number divided by Froude number. More...
 
const double & re_invfr () const
 Return the value of the Re/Fr number. More...
 
const Vector< double > & g () const
 Vector of gravitational components. More...
 
Vector< double > *& g_pt ()
 Pointer to Vector of gravitational components. More...
 
double flow ()
 Calculate the flow across the element. More...
 
double u (const unsigned &l, const unsigned &i)
 Access function for the velocity. N. B. HEAVY ASSUMPTIONS HERE. More...
 
double u (const unsigned &t, const unsigned &l, const unsigned &i) const
 
double du_dt (const unsigned &l, const unsigned &i) const
 i-th component of du/dt at local node l. More...
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Add the contribution to the residuals. More...
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 This function returns the residuals and the jacobian. More...
 
void add_generic_residual_contribution (Vector< double > &residuals, DenseMatrix< double > &jacobian, unsigned flag)
 
void assign_additional_local_eqn_numbers ()
 Define the local equation numbering schemes. More...
 
void output (std::ostream &outfile)
 Overload the output function. More...
 
void output (std::ostream &outfile, const unsigned &Np)
 Output function: x,y,[z],u,v,[w],p in tecplot format. More...
 
- Public Member Functions inherited from oomph::SpineElement< FaceGeometry< ELEMENT > >
 SpineElement ()
 Constructor, call the constructor of the base element. More...
 
 SpineElement (FiniteElement *const &element_pt, const int &face_index)
 Constructor used for spine face elements. More...
 
 ~SpineElement ()
 Destructor, clean up the storage allocated to the local equation numbers. More...
 
int spine_local_eqn (const unsigned &n)
 
- Public Member Functions inherited from oomph::ElementWithSpecificMovingNodes< ELEMENT, NODE_TYPE >
void describe_local_dofs (std::ostream &out, const std::string &current_string) const
 
 ElementWithSpecificMovingNodes ()
 Constructor, call the constructor of the base element. More...
 
 ElementWithSpecificMovingNodes (FiniteElement *const &element_pt, const int &face_index)
 Constructor used for face elements. More...
 
 ~ElementWithSpecificMovingNodes ()
 Empty Destructor,. More...
 
void describe_local_dofs (std::ostream &out, std::string &curr_str)
 Unique final overrider for describe_dofs. More...
 
Node * construct_node (const unsigned &n)
 
Node * construct_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 Overloaded node allocation for unsteady problems. More...
 
Node * construct_boundary_node (const unsigned &n)
 Overload the node assignment routine to assign boundary nodes. More...
 
Node * construct_boundary_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 Overloaded boundary node allocation for unsteady problems. More...
 
void assign_all_generic_local_eqn_numbers (const bool &store_local_dof_pt)
 
void get_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 Compute the element's residuals vector and jacobian matrix. More...
 
void get_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 Compute the element's residuals vector and jacobian matrix. More...
 
- Public Member Functions inherited from oomph::ElementWithMovingNodes
 ElementWithMovingNodes ()
 Constructor. More...
 
 ElementWithMovingNodes (const ElementWithMovingNodes &)=delete
 Broken copy constructor. More...
 
void operator= (const ElementWithMovingNodes &)=delete
 Broken assignment operator. More...
 
virtual ~ElementWithMovingNodes ()
 Virtual destructor (clean up and allocated memory) More...
 
unsigned ngeom_dof () const
 Number of geometric dofs. More...
 
int geometric_data_local_eqn (const unsigned &n, const unsigned &i)
 
void assemble_set_of_all_geometric_data (std::set< Data * > &unique_geom_data_pt)
 Return a set of all geometric data associated with the element. More...
 
void identify_geometric_data (std::set< Data * > &geometric_data_pt)
 
bool are_dresidual_dnodal_coordinates_always_evaluated_by_fd () const
 Return whether shape derivatives are evaluated by fd. More...
 
void enable_always_evaluate_dresidual_dnodal_coordinates_by_fd ()
 
void disable_always_evaluate_dresidual_dnodal_coordinates_by_fd ()
 
void evaluate_shape_derivs_by_direct_fd ()
 Evaluate shape derivatives by direct finite differencing. More...
 
void evaluate_shape_derivs_by_chain_rule (const bool &i_know_what_i_am_doing=false)
 
void evaluate_shape_derivs_by_fastest_method (const bool &i_know_what_i_am_doing=false)
 
int & method_for_shape_derivs ()
 Access to method (enumerated flag) for determination of shape derivs. More...
 
void enable_bypass_fill_in_jacobian_from_geometric_data ()
 Bypass the call to fill_in_jacobian_from_geometric_data. More...
 
void disable_bypass_fill_in_jacobian_from_geometric_data ()
 Do not bypass the call to fill_in_jacobian_from_geometric_data. More...
 
bool is_fill_in_jacobian_from_geometric_data_bypassed () const
 
unsigned ngeom_data () const
 
- Public Member Functions inherited from oomph::FiniteElement
void set_dimension (const unsigned &dim)
 
void set_nodal_dimension (const unsigned &nodal_dim)
 
void set_nnodal_position_type (const unsigned &nposition_type)
 Set the number of types required to interpolate the coordinate. More...
 
void set_n_node (const unsigned &n)
 
int nodal_local_eqn (const unsigned &n, const unsigned &i) const
 
double dJ_eulerian_at_knot (const unsigned &ipt, Shape &psi, DenseMatrix< double > &djacobian_dX) const
 
 FiniteElement ()
 Constructor. More...
 
virtual ~FiniteElement ()
 
 FiniteElement (const FiniteElement &)=delete
 Broken copy constructor. More...
 
virtual bool local_coord_is_valid (const Vector< double > &s)
 Broken assignment operator. More...
 
virtual void move_local_coord_back_into_element (Vector< double > &s) const
 
void get_centre_of_gravity_and_max_radius_in_terms_of_zeta (Vector< double > &cog, double &max_radius) const
 
virtual void local_coordinate_of_node (const unsigned &j, Vector< double > &s) const
 
virtual void local_fraction_of_node (const unsigned &j, Vector< double > &s_fraction)
 
virtual double local_one_d_fraction_of_node (const unsigned &n1d, const unsigned &i)
 
virtual void set_macro_elem_pt (MacroElement *macro_elem_pt)
 
MacroElement * macro_elem_pt ()
 Access function to pointer to macro element. More...
 
void get_x (const Vector< double > &s, Vector< double > &x) const
 
void get_x (const unsigned &t, const Vector< double > &s, Vector< double > &x)
 
virtual void get_x_from_macro_element (const Vector< double > &s, Vector< double > &x) const
 
virtual void get_x_from_macro_element (const unsigned &t, const Vector< double > &s, Vector< double > &x)
 
virtual void set_integration_scheme (Integral *const &integral_pt)
 Set the spatial integration scheme. More...
 
Integral *const & integral_pt () const
 Return the pointer to the integration scheme (const version) More...
 
virtual void shape (const Vector< double > &s, Shape &psi) const =0
 
virtual void shape_at_knot (const unsigned &ipt, Shape &psi) const
 
virtual void dshape_local (const Vector< double > &s, Shape &psi, DShape &dpsids) const
 
virtual void dshape_local_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsids) const
 
virtual void d2shape_local (const Vector< double > &s, Shape &psi, DShape &dpsids, DShape &d2psids) const
 
virtual void d2shape_local_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsids, DShape &d2psids) const
 
void check_J_eulerian_at_knots (bool &passed) const
 
void check_jacobian (const double &jacobian) const
 
double dshape_eulerian (const Vector< double > &s, Shape &psi, DShape &dpsidx) const
 
virtual double dshape_eulerian_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsidx) const
 
virtual double dshape_eulerian_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsi, DenseMatrix< double > &djacobian_dX, RankFourTensor< double > &d_dpsidx_dX) const
 
double d2shape_eulerian (const Vector< double > &s, Shape &psi, DShape &dpsidx, DShape &d2psidx) const
 
virtual double d2shape_eulerian_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsidx, DShape &d2psidx) const
 
virtual void assign_nodal_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual void describe_nodal_local_dofs (std::ostream &out, const std::string &current_string) const
 
Node *& node_pt (const unsigned &n)
 Return a pointer to the local node n. More...
 
Node *const & node_pt (const unsigned &n) const
 Return a pointer to the local node n (const version) More...
 
unsigned nnode () const
 Return the number of nodes. More...
 
virtual unsigned nnode_1d () const
 
double raw_nodal_position (const unsigned &n, const unsigned &i) const
 
double raw_nodal_position (const unsigned &t, const unsigned &n, const unsigned &i) const
 
double raw_dnodal_position_dt (const unsigned &n, const unsigned &i) const
 
double raw_dnodal_position_dt (const unsigned &n, const unsigned &j, const unsigned &i) const
 
double raw_nodal_position_gen (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double raw_nodal_position_gen (const unsigned &t, const unsigned &n, const unsigned &k, const unsigned &i) const
 
double raw_dnodal_position_gen_dt (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double raw_dnodal_position_gen_dt (const unsigned &j, const unsigned &n, const unsigned &k, const unsigned &i) const
 
double nodal_position (const unsigned &n, const unsigned &i) const
 
double nodal_position (const unsigned &t, const unsigned &n, const unsigned &i) const
 
double dnodal_position_dt (const unsigned &n, const unsigned &i) const
 Return the i-th component of nodal velocity: dx/dt at local node n. More...
 
double dnodal_position_dt (const unsigned &n, const unsigned &j, const unsigned &i) const
 
double nodal_position_gen (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double nodal_position_gen (const unsigned &t, const unsigned &n, const unsigned &k, const unsigned &i) const
 
double dnodal_position_gen_dt (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double dnodal_position_gen_dt (const unsigned &j, const unsigned &n, const unsigned &k, const unsigned &i) const
 
virtual void get_dresidual_dnodal_coordinates (RankThreeTensor< double > &dresidual_dnodal_coordinates)
 
virtual void disable_ALE ()
 
virtual void enable_ALE ()
 
virtual unsigned required_nvalue (const unsigned &n) 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 unsigned nvertex_node () const
 
virtual Node * vertex_node_pt (const unsigned &j) const
 
int get_node_number (Node *const &node_pt) const
 
virtual Node * get_node_at_local_coordinate (const Vector< double > &s) const
 
double raw_nodal_value (const unsigned &n, const unsigned &i) const
 
double raw_nodal_value (const unsigned &t, const unsigned &n, const unsigned &i) const
 
double nodal_value (const unsigned &n, const unsigned &i) const
 
double nodal_value (const unsigned &t, const unsigned &n, const unsigned &i) const
 
unsigned dim () const
 
virtual ElementGeometry::ElementGeometry element_geometry () const
 Return the geometry type of the element (either Q or T usually). More...
 
unsigned ngeom_data () const
 
Data * geom_data_pt (const unsigned &j)
 
void position (const Vector< double > &zeta, Vector< double > &r) const
 
void position (const unsigned &t, const Vector< double > &zeta, Vector< double > &r) const
 
void dposition_dt (const Vector< double > &zeta, const unsigned &t, Vector< double > &drdt)
 
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 identify_field_data_for_interactions (std::set< std::pair< Data *, unsigned >> &paired_field_data)
 
virtual double s_min () const
 Min value of local coordinate. More...
 
virtual double s_max () const
 Max. value of local coordinate. More...
 
double size () const
 
virtual double compute_physical_size () const
 
virtual void point_output_data (const Vector< double > &s, Vector< double > &data)
 
void point_output (std::ostream &outfile, const Vector< double > &s)
 
virtual unsigned nplot_points_paraview (const unsigned &nplot) const
 
virtual unsigned nsub_elements_paraview (const unsigned &nplot) const
 
void output_paraview (std::ofstream &file_out, const unsigned &nplot) const
 
virtual void write_paraview_output_offset_information (std::ofstream &file_out, const unsigned &nplot, unsigned &counter) const
 
virtual void write_paraview_type (std::ofstream &file_out, const unsigned &nplot) const
 
virtual void write_paraview_offsets (std::ofstream &file_out, const unsigned &nplot, unsigned &offset_sum) const
 
virtual unsigned nscalar_paraview () const
 
virtual void scalar_value_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot) const
 
virtual void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt) const
 
virtual void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt) const
 
virtual std::string scalar_name_paraview (const unsigned &i) const
 
virtual void output (const unsigned &t, std::ostream &outfile, const unsigned &n_plot) const
 
virtual void output (FILE *file_pt)
 
virtual void output (FILE *file_pt, const unsigned &n_plot)
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 Output an exact solution over the element. More...
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt)
 Output a time-dependent exact solution over the element. More...
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, const SolutionFunctorBase &exact_soln) const
 Output a time-dependent exact solution over the element. More...
 
virtual void get_s_plot (const unsigned &i, const unsigned &nplot, Vector< double > &s, const bool &shifted_to_interior=false) const
 
virtual std::string tecplot_zone_string (const unsigned &nplot) const
 
virtual void write_tecplot_zone_footer (std::ostream &outfile, const unsigned &nplot) const
 
virtual void write_tecplot_zone_footer (FILE *file_pt, const unsigned &nplot) const
 
virtual unsigned nplot_points (const unsigned &nplot) const
 
virtual void compute_error (FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 Calculate the norm of the error and that of the exact solution. More...
 
virtual void compute_error (FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 Calculate the norm of the error and that of the exact solution. More...
 
virtual void compute_error (FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_abs_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error)
 
void integrate_fct (FiniteElement::SteadyExactSolutionFctPt integrand_fct_pt, Vector< double > &integral)
 Integrate Vector-valued function over element. More...
 
void integrate_fct (FiniteElement::UnsteadyExactSolutionFctPt integrand_fct_pt, const double &time, Vector< double > &integral)
 Integrate Vector-valued time-dep function over element. More...
 
virtual void build_face_element (const int &face_index, FaceElement *face_element_pt)
 
virtual unsigned self_test ()
 
virtual unsigned get_bulk_node_number (const int &face_index, const unsigned &i) const
 
virtual int face_outer_unit_normal_sign (const int &face_index) const
 Get the sign of the outer unit normal on the face given by face_index. More...
 
virtual unsigned nnode_on_face () const
 
void face_node_number_error_check (const unsigned &i) const
 Range check for face node numbers. More...
 
virtual CoordinateMappingFctPt face_to_bulk_coordinate_fct_pt (const int &face_index) const
 
virtual BulkCoordinateDerivativesFctPt bulk_coordinate_derivatives_fct_pt (const int &face_index) const
 
- Public Member Functions inherited from oomph::GeneralisedElement
 GeneralisedElement ()
 Constructor: Initialise all pointers and all values to zero. More...
 
virtual ~GeneralisedElement ()
 Virtual destructor to clean up any memory allocated by the object. More...
 
 GeneralisedElement (const GeneralisedElement &)=delete
 Broken copy constructor. More...
 
void operator= (const GeneralisedElement &)=delete
 Broken assignment operator. More...
 
Data *& internal_data_pt (const unsigned &i)
 Return a pointer to i-th internal data object. More...
 
Data *const & internal_data_pt (const unsigned &i) const
 Return a pointer to i-th internal data object (const version) More...
 
Data *& external_data_pt (const unsigned &i)
 Return a pointer to i-th external data object. More...
 
Data *const & external_data_pt (const unsigned &i) const
 Return a pointer to i-th external data object (const version) More...
 
unsigned long eqn_number (const unsigned &ieqn_local) const
 
int local_eqn_number (const unsigned long &ieqn_global) const
 
unsigned add_external_data (Data *const &data_pt, const bool &fd=true)
 
bool external_data_fd (const unsigned &i) const
 
void exclude_external_data_fd (const unsigned &i)
 
void include_external_data_fd (const unsigned &i)
 
void flush_external_data ()
 Flush all external data. More...
 
void flush_external_data (Data *const &data_pt)
 Flush the object addressed by data_pt from the external data array. More...
 
unsigned ninternal_data () const
 Return the number of internal data objects. More...
 
unsigned nexternal_data () const
 Return the number of external data objects. More...
 
unsigned ndof () const
 Return the number of equations/dofs in the element. More...
 
void dof_vector (const unsigned &t, Vector< double > &dof)
 Return the vector of dof values at time level t. More...
 
void dof_pt_vector (Vector< double * > &dof_pt)
 Return the vector of pointers to dof values. More...
 
void set_internal_data_time_stepper (const unsigned &i, TimeStepper *const &time_stepper_pt, const bool &preserve_existing_data)
 
void assign_internal_eqn_numbers (unsigned long &global_number, Vector< double * > &Dof_pt)
 
void describe_dofs (std::ostream &out, const std::string &current_string) const
 
void add_internal_value_pt_to_map (std::map< unsigned, double * > &map_of_value_pt)
 
virtual void assign_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual void get_residuals (Vector< double > &residuals)
 
virtual void get_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &mass_matrix)
 
virtual void get_dresiduals_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam)
 
virtual void get_djacobian_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
 
virtual void get_djacobian_and_dmass_matrix_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam, DenseMatrix< double > &dmass_matrix_dparam)
 
virtual void get_hessian_vector_products (Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
 
virtual void get_inner_products (Vector< std::pair< unsigned, unsigned >> const &history_index, Vector< double > &inner_product)
 
virtual void get_inner_product_vectors (Vector< unsigned > const &history_index, Vector< Vector< double >> &inner_product_vector)
 
virtual void compute_norm (Vector< double > &norm)
 
virtual void compute_norm (double &norm)
 
virtual unsigned ndof_types () const
 
virtual void get_dof_numbers_for_unknowns (std::list< std::pair< unsigned long, unsigned >> &dof_lookup_list) const
 
- Public Member Functions inherited from oomph::GeomObject
 GeomObject ()
 Default constructor. More...
 
 GeomObject (const unsigned &ndim)
 
 GeomObject (const unsigned &nlagrangian, const unsigned &ndim)
 
 GeomObject (const unsigned &nlagrangian, const unsigned &ndim, TimeStepper *time_stepper_pt)
 
 GeomObject (const GeomObject &dummy)=delete
 Broken copy constructor. More...
 
void operator= (const GeomObject &)=delete
 Broken assignment operator. More...
 
virtual ~GeomObject ()
 (Empty) destructor More...
 
unsigned nlagrangian () const
 Access function to # of Lagrangian coordinates. More...
 
unsigned ndim () const
 Access function to # of Eulerian coordinates. More...
 
void set_nlagrangian_and_ndim (const unsigned &n_lagrangian, const unsigned &n_dim)
 Set # of Lagrangian and Eulerian coordinates. More...
 
TimeStepper *& time_stepper_pt ()
 
TimeStepper * time_stepper_pt () const
 
virtual void position (const double &t, const Vector< double > &zeta, Vector< double > &r) const
 
virtual void dposition (const Vector< double > &zeta, DenseMatrix< double > &drdzeta) const
 
virtual void d2position (const Vector< double > &zeta, RankThreeTensor< double > &ddrdzeta) const
 
virtual void d2position (const Vector< double > &zeta, Vector< double > &r, DenseMatrix< double > &drdzeta, RankThreeTensor< double > &ddrdzeta) const
 
- Public Member Functions inherited from oomph::SpineFiniteElement
 SpineFiniteElement ()
 Empty constructor. More...
 
virtual ~SpineFiniteElement ()
 Emtpty virtual destructor. More...
 
- Public Member Functions inherited from oomph::FaceElement
 FaceElement ()
 Constructor: Initialise all appropriate member data. More...
 
virtual ~FaceElement ()
 Empty virtual destructor. More...
 
 FaceElement (const FaceElement &)=delete
 Broken copy constructor. More...
 
const unsigned & boundary_number_in_bulk_mesh () const
 Broken assignment operator. More...
 
void set_boundary_number_in_bulk_mesh (const unsigned &b)
 Set function for the boundary number in bulk mesh. More...
 
double zeta_nodal (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double J_eulerian (const Vector< double > &s) const
 
double J_eulerian_at_knot (const unsigned &ipt) const
 
void check_J_eulerian_at_knots (bool &passed) const
 
double interpolated_x (const Vector< double > &s, const unsigned &i) const
 
double interpolated_x (const unsigned &t, const Vector< double > &s, const unsigned &i) const
 
void interpolated_x (const Vector< double > &s, Vector< double > &x) const
 
void interpolated_x (const unsigned &t, const Vector< double > &s, Vector< double > &x) const
 
double interpolated_dxdt (const Vector< double > &s, const unsigned &i, const unsigned &t)
 
void interpolated_dxdt (const Vector< double > &s, const unsigned &t, Vector< double > &dxdt)
 
int & normal_sign ()
 
int normal_sign () const
 
int & face_index ()
 
int face_index () const
 
const Vector< double > * tangent_direction_pt () const
 Public access function for the tangent direction pointer. More...
 
void set_tangent_direction (Vector< double > *tangent_direction_pt)
 Set the tangent direction vector. More...
 
void turn_on_warning_for_discontinuous_tangent ()
 
void turn_off_warning_for_discontinuous_tangent ()
 
void continuous_tangent_and_outer_unit_normal (const Vector< double > &s, Vector< Vector< double >> &tang_vec, Vector< double > &unit_normal) const
 
void continuous_tangent_and_outer_unit_normal (const unsigned &ipt, Vector< Vector< double >> &tang_vec, Vector< double > &unit_normal) const
 
void outer_unit_normal (const Vector< double > &s, Vector< double > &unit_normal) const
 Compute outer unit normal at the specified local coordinate. More...
 
void outer_unit_normal (const unsigned &ipt, Vector< double > &unit_normal) const
 Compute outer unit normal at ipt-th integration point. More...
 
FiniteElement *& bulk_element_pt ()
 Pointer to higher-dimensional "bulk" element. More...
 
FiniteElement * bulk_element_pt () const
 Pointer to higher-dimensional "bulk" element (const version) More...
 
CoordinateMappingFctPt & face_to_bulk_coordinate_fct_pt ()
 
CoordinateMappingFctPt face_to_bulk_coordinate_fct_pt () const
 
BulkCoordinateDerivativesFctPt & bulk_coordinate_derivatives_fct_pt ()
 
BulkCoordinateDerivativesFctPt bulk_coordinate_derivatives_fct_pt () const
 
Vector< double > local_coordinate_in_bulk (const Vector< double > &s) const
 
void get_local_coordinate_in_bulk (const Vector< double > &s, Vector< double > &s_bulk) const
 
void get_ds_bulk_ds_face (const Vector< double > &s, DenseMatrix< double > &dsbulk_dsface, unsigned &interior_direction) const
 
unsigned & bulk_position_type (const unsigned &i)
 
const unsigned & bulk_position_type (const unsigned &i) const
 
void bulk_node_number_resize (const unsigned &i)
 Resize the storage for the bulk node numbers. More...
 
unsigned & bulk_node_number (const unsigned &n)
 
const unsigned & bulk_node_number (const unsigned &n) const
 
void bulk_position_type_resize (const unsigned &i)
 Resize the storage for bulk_position_type to i entries. More...
 
unsigned & nbulk_value (const unsigned &n)
 
unsigned nbulk_value (const unsigned &n) const
 
void nbulk_value_resize (const unsigned &i)
 
void resize_nodes (Vector< unsigned > &nadditional_data_values)
 
void output_zeta (std::ostream &outfile, const unsigned &nplot)
 Output boundary coordinate zeta. More...
 

Protected Member Functions

int invca_local_eqn ()
 
- Protected Member Functions inherited from oomph::ElementWithMovingNodes
virtual void get_dnodal_coordinates_dgeom_dofs (RankThreeTensor< double > &dnodal_coordinates_dgeom_dofs)
 
void fill_in_jacobian_from_geometric_data (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void fill_in_jacobian_from_geometric_data (DenseMatrix< double > &jacobian)
 
- Protected Member Functions inherited from oomph::FiniteElement
virtual void assemble_local_to_eulerian_jacobian (const DShape &dpsids, DenseMatrix< double > &jacobian) const
 
virtual void assemble_local_to_eulerian_jacobian2 (const DShape &d2psids, DenseMatrix< double > &jacobian2) const
 
virtual void assemble_eulerian_base_vectors (const DShape &dpsids, DenseMatrix< double > &interpolated_G) const
 
template<unsigned DIM>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
virtual double invert_jacobian_mapping (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
virtual double local_to_eulerian_mapping (const DShape &dpsids, DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
double local_to_eulerian_mapping (const DShape &dpsids, DenseMatrix< double > &inverse_jacobian) const
 
virtual double local_to_eulerian_mapping_diagonal (const DShape &dpsids, DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
virtual void dJ_eulerian_dnodal_coordinates (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<unsigned DIM>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
virtual void d_dshape_eulerian_dnodal_coordinates (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<unsigned DIM>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
virtual void transform_derivatives (const DenseMatrix< double > &inverse_jacobian, DShape &dbasis) const
 
void transform_derivatives_diagonal (const DenseMatrix< double > &inverse_jacobian, DShape &dbasis) const
 
virtual void transform_second_derivatives (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<unsigned DIM>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<unsigned DIM>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
virtual void fill_in_jacobian_from_nodal_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void fill_in_jacobian_from_nodal_by_fd (DenseMatrix< double > &jacobian)
 
virtual void update_before_nodal_fd ()
 
virtual void reset_after_nodal_fd ()
 
virtual void update_in_nodal_fd (const unsigned &i)
 
virtual void reset_in_nodal_fd (const unsigned &i)
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 Zero-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 One-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 Two-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
- Protected Member Functions inherited from oomph::GeneralisedElement
unsigned add_internal_data (Data *const &data_pt, const bool &fd=true)
 
bool internal_data_fd (const unsigned &i) const
 
void exclude_internal_data_fd (const unsigned &i)
 
void include_internal_data_fd (const unsigned &i)
 
void clear_global_eqn_numbers ()
 
void add_global_eqn_numbers (std::deque< unsigned long > const &global_eqn_numbers, std::deque< double * > const &global_dof_pt)
 
virtual void assign_internal_and_external_local_eqn_numbers (const bool &store_local_dof_pt)
 
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_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
virtual void fill_in_contribution_to_dresiduals_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam)
 
virtual void fill_in_contribution_to_djacobian_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
 
virtual void fill_in_contribution_to_djacobian_and_dmass_matrix_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam, DenseMatrix< double > &dmass_matrix_dparam)
 
virtual void fill_in_contribution_to_hessian_vector_products (Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
 
virtual void fill_in_contribution_to_inner_products (Vector< std::pair< unsigned, unsigned >> const &history_index, Vector< double > &inner_product)
 
virtual void fill_in_contribution_to_inner_product_vectors (Vector< unsigned > const &history_index, Vector< Vector< double >> &inner_product_vector)
 
- Protected Member Functions inherited from oomph::FaceElement
void add_additional_values (const Vector< unsigned > &nadditional_values, const unsigned &id)
 

Protected Attributes

DenseMatrix< int > U_local_eqn
 
int Delta_P_local_eqn
 
unsigned External_Delta_P_index
 
DenseMatrix< int > External_u_local_eqn
 
Vector< unsigned > External_node
 
unsigned External_data_number_of_invca
 
Data * invca_data_pt
 Pointer to the Data item that stores the capillary number. More...
 
double * bond_pt
 
- Protected Attributes inherited from oomph::ElementWithMovingNodes
Vector< Data * > Geom_data_pt
 
- Protected Attributes inherited from oomph::FiniteElement
MacroElement * Macro_elem_pt
 Pointer to the element's macro element (NULL by default) More...
 
- Protected Attributes inherited from oomph::GeomObject
unsigned NLagrangian
 Number of Lagrangian (intrinsic) coordinates. More...
 
unsigned Ndim
 Number of Eulerian coordinates. More...
 
TimeStepper * Geom_object_time_stepper_pt
 
- Protected Attributes inherited from oomph::FaceElement
unsigned Boundary_number_in_bulk_mesh
 The boundary number in the bulk mesh to which this element is attached. More...
 
FiniteElement * Bulk_element_pt
 Pointer to the associated higher-dimensional "bulk" element. More...
 
Vector< unsigned > Bulk_node_number
 
Vector< unsigned > Nbulk_value
 
Vector< double > * Tangent_direction_pt
 

Private Attributes

unsigned Dim
 The highest dimension of the problem. More...
 
double * ReInvFr_pt
 Pointer to the global Reynold number divided by the Froude number. More...
 
Vector< double > * G_pt
 Pointer to global gravity Vector. More...
 
double * Viscosity_Ratio_pt
 
double * Density_Ratio_pt
 
Data * Delta_P_pt
 

Additional Inherited Members

- Public Types inherited from oomph::ElementWithMovingNodes
enum  { Shape_derivs_by_chain_rule , Shape_derivs_by_direct_fd , Shape_derivs_by_fastest_method }
 Public enumeration to choose method for computing shape derivatives. More...
 
- Public Types inherited from oomph::FiniteElement
typedef void(* SteadyExactSolutionFctPt) (const Vector< double > &, Vector< double > &)
 
typedef void(* UnsteadyExactSolutionFctPt) (const double &, const Vector< double > &, Vector< double > &)
 
- Static Public 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 Protected Attributes inherited from oomph::FiniteElement
static const unsigned Default_Initial_Nvalue = 0
 Default value for the number of values at a node. More...
 
static const double Node_location_tolerance = 1.0e-14
 
static const unsigned N2deriv [] = {0, 1, 3, 6}
 
- Static Protected Attributes inherited from oomph::GeneralisedElement
static DenseMatrix< double > Dummy_matrix
 
static std::deque< double * > Dof_pt_deque
 

Detailed Description

template<class ELEMENT>
class SpineGravityTractionElement< ELEMENT >

Special Face element used to calculate the additional inlet velocities and tractions when gravity has a component in the axial direction

Constructor & Destructor Documentation

◆ SpineGravityTractionElement() [1/2]

template<class ELEMENT >
SpineGravityTractionElement< ELEMENT >::SpineGravityTractionElement ( FiniteElement *const &  element_pt,
const int &  face_index 
)
inline

Constructor, which takes a "bulk" element and the value of the index and its limit

393  :
394  SpineElement<FaceGeometry<ELEMENT> >(), FaceElement()
395  {
396  //Attach the geometrical information to the element. N.B. This function
397  //also assigns nbulk_data from the required_nvalue of the bulk element
398  element_pt->build_face_element(face_index,this);
399 
400  //Set the dimension from the dimension of the first node
401  Dim = node_pt(0)->ndim();
402 
403  //Set the Physical values from the bulk elemenet
404  ELEMENT* cast_element_pt = dynamic_cast<ELEMENT*>(element_pt);
405  this->ReInvFr_pt = cast_element_pt->re_invfr_pt();
406  this->G_pt = cast_element_pt->g_pt();
407  this->Viscosity_Ratio_pt = cast_element_pt->viscosity_ratio_pt();
408  this->Density_Ratio_pt = cast_element_pt->density_ratio_pt();
409 
410  //Initialise the external pressure gradient pointer to zero
411  Delta_P_pt=0;
412 
413  //Hijack the nodes in the bulk element in the axial coordinate
414  unsigned n_node = this->nnode();
415  for(unsigned m=0;m<n_node;m++)
416  {
417  delete cast_element_pt->hijack_nodal_value(bulk_node_number(m),0);
418  }
419 
420  //The other nodes of the bulk element must be external data, because
421  //they can affect the derivatives that are used in this element
422  //Loop over the nodes of the parent element
423  unsigned n_node_parent = cast_element_pt->nnode();
424  for(unsigned n=0;n<n_node_parent;n++)
425  {
426  bool external=true;
427  //Loop over the face nodes
428  for(unsigned m=0;m<n_node;m++)
429  {
430  //If the parent's node is one of the face nodes continue
431  if(n == this->bulk_node_number(m)) {external=false;}
432  }
433 
434  //If it's external data add it, but do not finite difference
435  if(external)
436  {
437  this->add_external_data(cast_element_pt->node_pt(n),false);
438  External_node.push_back(n);
439  }
440  }
441 
442  //Now add the spines of the bulk elemnt as external data, which
443  //we will finite difference
444 
445  //Set of unique geometric data that is used to update the bulk,
446  //but is not used to update the face
447  std::set<Data*> unique_additional_geom_data;
448  //Get all the geometric data for this (bulk) element
449  cast_element_pt->assemble_set_of_all_geometric_data(
450  unique_additional_geom_data);
451 
452  //Now assemble the set of geometric data for the face element
453  std::set<Data*> unique_face_geom_data_pt;
454  this->assemble_set_of_all_geometric_data(unique_face_geom_data_pt);
455  //Erase the face geometric data from the additional data
456  for(std::set<Data*>::iterator it=unique_face_geom_data_pt.begin();
457  it!=unique_face_geom_data_pt.end();++it)
458  {unique_additional_geom_data.erase(*it);}
459 
460  //Finally add all unique additional data as external data
461  for(std::set<Data*>::iterator it = unique_additional_geom_data.begin();
462  it!= unique_additional_geom_data.end();++it)
463  {
464  this->add_external_data(*it);
465  }
466 
467  }
n
const unsigned n
Definition: CG3DPackingUnitTest.cpp:11
SpineGravityTractionElement::Delta_P_pt
Data * Delta_P_pt
Definition: elastic_bretherton.cc:364
SpineGravityTractionElement::Viscosity_Ratio_pt
double * Viscosity_Ratio_pt
Definition: elastic_bretherton.cc:356
SpineGravityTractionElement::Dim
unsigned Dim
The highest dimension of the problem.
Definition: elastic_bretherton.cc:346
SpineGravityTractionElement::ReInvFr_pt
double * ReInvFr_pt
Pointer to the global Reynold number divided by the Froude number.
Definition: elastic_bretherton.cc:349
SpineGravityTractionElement::External_node
Vector< unsigned > External_node
Definition: elastic_bretherton.cc:386
SpineGravityTractionElement::G_pt
Vector< double > * G_pt
Pointer to global gravity Vector.
Definition: elastic_bretherton.cc:352
SpineGravityTractionElement::Density_Ratio_pt
double * Density_Ratio_pt
Definition: elastic_bretherton.cc:360
oomph::ElementWithMovingNodes::assemble_set_of_all_geometric_data
void assemble_set_of_all_geometric_data(std::set< Data * > &unique_geom_data_pt)
Return a set of all geometric data associated with the element.
Definition: element_with_moving_nodes.cc:43
oomph::FaceElement::face_index
int & face_index()
Definition: elements.h:4626
oomph::FaceElement::FaceElement
FaceElement()
Constructor: Initialise all appropriate member data.
Definition: elements.h:4446
oomph::FaceElement::bulk_node_number
unsigned & bulk_node_number(const unsigned &n)
Definition: elements.h:4825
oomph::FiniteElement::node_pt
Node *& node_pt(const unsigned &n)
Return a pointer to the local node n.
Definition: elements.h:2175
oomph::FiniteElement::build_face_element
virtual void build_face_element(const int &face_index, FaceElement *face_element_pt)
Definition: elements.cc:5132
oomph::FiniteElement::nnode
unsigned nnode() const
Return the number of nodes.
Definition: elements.h:2210
oomph::GeneralisedElement::add_external_data
unsigned add_external_data(Data *const &data_pt, const bool &fd=true)
Definition: elements.cc:307
oomph::Node::ndim
unsigned ndim() const
Return (Eulerian) spatial dimension of the node.
Definition: nodes.h:1054
oomph::SpineElement
Definition: spines.h:477
m
int * m
Definition: level2_cplx_impl.h:294

References oomph::FiniteElement::build_face_element(), Global_Variables::Dim, m, and n.

◆ SpineGravityTractionElement() [2/2]

template<class ELEMENT >
SpineGravityTractionElement< ELEMENT >::SpineGravityTractionElement ( FiniteElement *  element_pt,
int  face_index 
)
inline

Constructor, which takes a "bulk" element and the value of the index and its limit

304  :
305  SpineElement<FaceGeometry<ELEMENT> >(), FaceElement()
306  {
307 
308  //Attach the geometrical information to the element. N.B. This function
309  //also assigns nbulk_data from the required_nvalue of the bulk element
310  element_pt->build_face_element(face_index,this);
311  //Set the dimension from the dimension of the first node
312  Dim = node_pt(0)->ndim();
313 
314  //Initializa pointer
315  invca_data_pt = 0;
316 
317  //Set the Physical values from the bulk elemenet
318  ELEMENT* cast_element_pt = dynamic_cast<ELEMENT*>(element_pt);
319  this->ReInvFr_pt = cast_element_pt->re_invfr_pt();
320  this->G_pt = cast_element_pt->g_pt();
321  this->Viscosity_Ratio_pt = cast_element_pt->viscosity_ratio_pt();
322  this->Density_Ratio_pt = cast_element_pt->density_ratio_pt();
323 
324 
325  //Hijack the nodes in the bulk element in the axial coordinate
326  unsigned n_node = this->nnode();
327 
328  for(unsigned m=0;m<n_node;m++)
329  {
330  delete cast_element_pt->hijack_nodal_value(bulk_node_number(m),1);
331  }
332 
333  //The other nodes of the bulk element must be external data, because
334  //they can affect the derivatives that are used in this element
335  //Loop over the nodes of the parent element
336  unsigned n_node_parent = cast_element_pt->nnode();
337  for(unsigned n=0;n<n_node_parent;n++)
338  {
339  bool external=true;
340  //Loop over the face nodes
341  for(unsigned m=0;m<n_node;m++)
342  {
343  //If the parent's node is one of the face nodes continue
344  if(n == this->bulk_node_number(m)) {external=false;}
345  }
346 
347  //If it's external data add it
348  if(external)
349  {
350  this->add_external_data(cast_element_pt->node_pt(n));
351  External_node.push_back(n);
352  }
353  }
354 
355  //Now add the spines of the bulk elemnt as external data, which
356  //we will finite difference
357 
358  //Set of unique geometric data that is used to update the bulk,
359  //but is not used to update the face
360  std::set<Data*> unique_additional_geom_data;
361  //Get all the geometric data for this (bulk) element
362  cast_element_pt->assemble_set_of_all_geometric_data(
363  unique_additional_geom_data);
364 
365  //Now assemble the set of geometric data for the face element
366  std::set<Data*> unique_face_geom_data_pt;
367  this->assemble_set_of_all_geometric_data(unique_face_geom_data_pt);
368  //Erase the face geometric data from the additional data
369  for(std::set<Data*>::iterator it=unique_face_geom_data_pt.begin();
370  it!=unique_face_geom_data_pt.end();++it)
371  {unique_additional_geom_data.erase(*it);}
372 
373  //Finally add all unique additional data as geometric data
374  /*for(std::set<Data*>::iterator it = unique_additional_geom_data.begin();
375  it!= unique_additional_geom_data.end();++it)
376  {
377  this->add_external_data(*it);
378  }*/
379  this->identify_geometric_data(unique_additional_geom_data);
380  }
SpineGravityTractionElement::invca_data_pt
Data * invca_data_pt
Pointer to the Data item that stores the capillary number.
Definition: extra_elements.h:283
oomph::FiniteElement::identify_geometric_data
virtual void identify_geometric_data(std::set< Data * > &geometric_data_pt)
Definition: elements.h:2789

References oomph::FiniteElement::build_face_element(), Global_Variables::Dim, m, and n.

Member Function Documentation

◆ add_generic_residual_contribution() [1/2]

template<class ELEMENT >
void SpineGravityTractionElement< ELEMENT >::add_generic_residual_contribution ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian,
unsigned  flag 
)
inline

This function returns the residuals for the Navier–Stokes equations; flag=1(or 0): do (or don't) compute the Jacobian as well.

586  {
587  //Find out how many nodes there are
588  unsigned n_node = nnode();
589 
590  //Set the value of n_intpt
591  unsigned n_intpt = integral_pt()->nweight();
592 
593  //Set the Vector to hold local coordinates
594  Vector<double> s(Dim-1);
595 
596  //Set the Vector to hold the local coordinates of the parent element
597  Vector<double> s_parent(Dim);
598 
599  //Get a pointer to the parent element
600  ELEMENT* bulk_el_pt = dynamic_cast<ELEMENT*>(bulk_element_pt());
601 
602  //Find the number of nodes in the parent element
603  unsigned n_node_parent = bulk_el_pt->nnode();
604  //Set up memory for shape functions and their derivatives
605  Shape psif_parent(n_node_parent), testf_parent(n_node_parent);
606  DShape dpsifdx_parent(n_node_parent,Dim),
607  dtestfdx_parent(n_node_parent,Dim);
608 
609  //Storage for the local equation number
610  int local_eqn=0, local_unknown=0;
611 
612  //Get the Physical Variable
613  double ReInvFr = re_invfr()*density_ratio();
614  double Viscosity_Ratio = viscosity_ratio();
615  Vector<double> G = g();
616  double Delta_p = delta_p();
617 
618  unsigned n_external_u = nexternal_u_data();
619 
620  //Loop over the integration points
621  for(unsigned ipt=0;ipt<n_intpt;ipt++)
622  {
623  //Assign values of s
624  for(unsigned i=0;i<(Dim-1);i++) {s[i] = integral_pt()->knot(ipt,i);}
625 
626  //Get the values of s_parent
627  this->get_local_coordinate_in_bulk(s,s_parent);
628 
629  //Get the integral weight
630  double w = integral_pt()->weight(ipt);
631 
632  //Find the shape functions and derivatives of the parent
633  (void)bulk_el_pt->
634  dshape_eulerian(s_parent,psif_parent,dpsifdx_parent);
635 
636  //Get the local jacobian for the FaceElement
637  double J = J_eulerian(s);
638 
639  //Premultiply the weights and the Jacobian
640  double W = w*J;
641 
642  //Need only to find the velocity derivatives
643  DenseMatrix<double> interpolated_dudx(Dim,Dim,0.0);
644 
645  Vector<double> interpolated_u(Dim,0.0);
646 
647  //Calculate velocities and derivatives
648  for(unsigned l=0;l<n_node_parent;l++)
649  {
650  //Loop over velocity components
651  for(unsigned i=0;i<Dim;i++)
652  {
653  for(unsigned j=0;j<Dim;j++)
654  {
655  interpolated_dudx(i,j) +=
656  bulk_el_pt->u_nst(l,i)*dpsifdx_parent(l,j);
657  }
658  dtestfdx_parent(l,i) = dpsifdx_parent(l,i);
659  interpolated_u[i] += bulk_el_pt->u_nst(l,i)*psif_parent(l);
660  }
661  //Set the test functions to be the same as the shape functions
662  testf_parent[l] = psif_parent[l];
663  }
664 
665  //Storage for the outer unit normal
666  Vector<double> normal(Dim);
667  outer_unit_normal(s,normal);
668 
669  //Loop over the test functions
670  for(unsigned l=0;l<n_node;l++)
671  {
672  //Do the first (axial) velocity component -- Poisson equation
673  {
674  unsigned i=0;
675 
676  local_eqn = U_local_eqn(l,i);
677  /*IF it's not a boundary condition*/
678  if(local_eqn >= 0)
679  {
680  //Add the pressure gradient
681  residuals[local_eqn] +=
682  Delta_p*testf_parent[bulk_node_number(l)]*W;
683 
684  //Add the gravitational body force term
685  residuals[local_eqn] +=
686  ReInvFr*testf_parent[bulk_node_number(l)]*G[i]*W;
687 
688  //Add in the Poisson terms (Only in y(1) direction)
689  residuals[local_eqn] -= Viscosity_Ratio*
690  (interpolated_dudx(i,1) /*+ interpolated_dudx(1,i)*/)
691  *dtestfdx_parent(bulk_node_number(l),1)*W;
692 
693  //Now add the jacobian terms
694  if(flag)
695  {
696  //Loop over all nodes again
697  for(unsigned l2=0;l2<n_node;l2++)
698  {
699  {
700  unsigned i2=0;
701  local_unknown = U_local_eqn(l2,i2);
702  if(local_unknown >= 0)
703  {
704  jacobian(local_eqn,local_unknown) -=
705  Viscosity_Ratio*(dpsifdx_parent(bulk_node_number(l2),1)*
706  dtestfdx_parent(bulk_node_number(l),1))*W;
707  }
708  }
709 
710  /*{
711  unsigned i2=1;
712  local_unknown = U_local_eqn(l2,i2);
713  if(local_unknown >= 0)
714  {
715  jacobian(local_eqn,local_unknown) -=
716  Viscosity_Ratio*(dpsifdx_parent(bulk_node_number(l2),i)*
717  dtestfdx_parent(bulk_node_number(l),1))*W;
718  }
719  }*/
720  }
721 
722  //Loop over external data
723  for(unsigned l2=0;l2<n_external_u;l2++) {
724  {
725  unsigned i2=0;
726  local_unknown = External_u_local_eqn(l2,i2);
727  if(local_unknown >= 0)
728  {
729  jacobian(local_eqn,local_unknown) -=
730  Viscosity_Ratio*(dpsifdx_parent(External_node[l2],1)*
731  dtestfdx_parent(bulk_node_number(l),1))*W;
732  }
733  }
734 
735  /*{
736  unsigned i2=1;
737  local_unknown = External_u_local_eqn(l2,i2);
738  if(local_unknown >= 0)
739  {
740  jacobian(local_eqn,local_unknown) -=
741  Viscosity_Ratio*(dpsifdx_parent(External_node[l2],i)*
742  dtestfdx_parent(bulk_node_number(l),1))*W;
743  }
744  }*/
745  }
746 
747  //Add in the pressure gradient term
748  local_unknown = Delta_P_local_eqn;
749  if(local_unknown >= 0)
750  {
751  jacobian(local_eqn,local_unknown) +=
752  testf_parent[bulk_node_number(l)]*W;
753  }
754  } //End of if flag
755  }
756 
757  }
758 
759  //Now do the second (traction) component
760  /* {
761  unsigned i=1;
762 
763  local_eqn = U_local_eqn(l,i);
764  //If it's not a boundary condition
765  if(local_eqn >= 0)
766  {
767  for(unsigned k=0;k<Dim;k++)
768  {
769  residuals[local_eqn] += Viscosity_Ratio*
770  (interpolated_dudx(i,k) + interpolated_dudx(k,i))*normal[k]
771  *testf_parent(bulk_node_number(l))*W;
772  }
773 
774  //Now add the jacobian terms
775  if(flag)
776  {
777  //Loop over all nodes again
778  for(unsigned l2=0;l2<n_node;l2++)
779  {
780  {
781  unsigned i2=0;
782  local_unknown = U_local_eqn(l2,i2);
783  if(local_unknown >= 0)
784  {
785  jacobian(local_eqn,local_unknown) +=
786  Viscosity_Ratio*dpsifdx_parent(bulk_node_number(l2),i)*
787  normal[i2]*testf_parent(bulk_node_number(l))*W;
788  }
789  }
790 
791  {
792  unsigned i2=1;
793  local_unknown = U_local_eqn(l2,i2);
794  if(local_unknown >= 0)
795  {
796  for(unsigned k=0;k<Dim;k++)
797  {
798  jacobian(local_eqn,local_unknown) +=
799  Viscosity_Ratio*dpsifdx_parent(bulk_node_number(l2),k)*
800  normal[k]*testf_parent(bulk_node_number(l))*W;
801  }
802 
803  jacobian(local_eqn,local_unknown) +=
804  Viscosity_Ratio*dpsifdx_parent(bulk_node_number(l2),i)*
805  normal[i2]*testf_parent(bulk_node_number(l))*W;
806  }
807  }
808  }
809 
810 
811  //Loop over all external data
812  for(unsigned l2=0;l2<n_external_u;l2++)
813  {
814  {
815  unsigned i2=0;
816  local_unknown = External_u_local_eqn(l2,i2);
817  if(local_unknown >= 0)
818  {
819  jacobian(local_eqn,local_unknown) +=
820  Viscosity_Ratio*dpsifdx_parent(External_node[l2],i)*
821  normal[i2]*testf_parent(bulk_node_number(l))*W;
822  }
823  }
824 
825  {
826  unsigned i2=1;
827  local_unknown = External_u_local_eqn(l2,i2);
828  if(local_unknown >= 0)
829  {
830  for(unsigned k=0;k<Dim;k++)
831  {
832  jacobian(local_eqn,local_unknown) +=
833  Viscosity_Ratio*dpsifdx_parent(External_node[l2],k)*
834  normal[k]*testf_parent(bulk_node_number(l))*W;
835  }
836 
837  jacobian(local_eqn,local_unknown) +=
838  Viscosity_Ratio*dpsifdx_parent(External_node[l2],i)*
839  normal[i2]*testf_parent(bulk_node_number(l))*W;
840  }
841  }
842  }
843  }
844  }
845  }*/
846  }
847 
848  //Add the contribution to the pressure gradient constraint,
849  //which states that the flux must be 2.0
850  local_eqn = Delta_P_local_eqn;
851  if(local_eqn >= 0)
852  {
853  residuals[local_eqn] += interpolated_u[0]*W;
854 
855  //Add the jacobian terms
856  if(flag)
857  {
858  //Loop over the nodes again
859  for(unsigned l2=0;l2<n_node;l2++)
860  {
861  //Only the x-velocities at the nodes on the boundary
862  //affect the residuals
863  local_unknown = U_local_eqn(l2,0);
864  if(local_unknown >= 0)
865  {
866  //Add the appropriate jacobian term
867  jacobian(local_eqn,local_unknown) +=
868  psif_parent(bulk_node_number(l2))*W;
869  }
870  }
871  }
872  }
873  }
874  }
i
int i
Definition: BiCGSTAB_step_by_step.cpp:9
G
JacobiRotation< float > G
Definition: Jacobi_makeGivens.cpp:2
J
JacobiRotation< float > J
Definition: Jacobi_makeJacobi.cpp:3
w
RowVector3d w
Definition: Matrix_resize_int.cpp:3
SpineGravityTractionElement::re_invfr
const double & re_invfr()
Reynolds number divided by Froude number.
Definition: elastic_bretherton.cc:516
SpineGravityTractionElement::delta_p
const double delta_p() const
Return the value of the pressure gradient.
Definition: elastic_bretherton.cc:489
SpineGravityTractionElement::U_local_eqn
DenseMatrix< int > U_local_eqn
Definition: elastic_bretherton.cc:370
SpineGravityTractionElement::viscosity_ratio
const double & viscosity_ratio() const
Definition: elastic_bretherton.cc:500
SpineGravityTractionElement::g
const Vector< double > & g() const
Vector of gravitational components.
Definition: elastic_bretherton.cc:519
SpineGravityTractionElement::density_ratio
const double & density_ratio() const
Definition: elastic_bretherton.cc:507
SpineGravityTractionElement::External_u_local_eqn
DenseMatrix< int > External_u_local_eqn
Definition: elastic_bretherton.cc:382
SpineGravityTractionElement::nexternal_u_data
unsigned nexternal_u_data()
Return the number of external velocity data.
Definition: elastic_bretherton.cc:496
SpineGravityTractionElement::Delta_P_local_eqn
int Delta_P_local_eqn
Definition: elastic_bretherton.cc:374
oomph::DShape
Definition: shape.h:278
oomph::FaceElement::outer_unit_normal
void outer_unit_normal(const Vector< double > &s, Vector< double > &unit_normal) const
Compute outer unit normal at the specified local coordinate.
Definition: elements.cc:6006
oomph::FaceElement::bulk_element_pt
FiniteElement *& bulk_element_pt()
Pointer to higher-dimensional "bulk" element.
Definition: elements.h:4735
oomph::FaceElement::J_eulerian
double J_eulerian(const Vector< double > &s) const
Definition: elements.cc:5242
oomph::FaceElement::get_local_coordinate_in_bulk
void get_local_coordinate_in_bulk(const Vector< double > &s, Vector< double > &s_bulk) const
Definition: elements.cc:6384
oomph::FiniteElement::integral_pt
Integral *const & integral_pt() const
Return the pointer to the integration scheme (const version)
Definition: elements.h:1963
oomph::FiniteElement::dshape_eulerian
double dshape_eulerian(const Vector< double > &s, Shape &psi, DShape &dpsidx) const
Definition: elements.cc:3298
oomph::Integral::knot
virtual double knot(const unsigned &i, const unsigned &j) const =0
Return local coordinate s[j] of i-th integration point.
oomph::Integral::nweight
virtual unsigned nweight() const =0
Return the number of integration points of the scheme.
oomph::Integral::weight
virtual double weight(const unsigned &i) const =0
Return weight of i-th integration point.
oomph::Shape
Definition: shape.h:76
s
RealScalar s
Definition: level1_cplx_impl.h:130
GlobalPhysicalVariables::ReInvFr
double ReInvFr
Product of Rynolds number and inverse of Froude number.
Definition: axisym_navier_stokes/counter_rotating_disks/counter_rotating_disks.cc:67
Global_Physical_Variables::Viscosity_Ratio
double Viscosity_Ratio
Definition: elastic_two_layer_interface_axisym.cc:76
WallFunction::normal
void normal(const Vector< double > &x, Vector< double > &normal)
Definition: free_surface_rotation.cc:65
oomph::QuadTreeNames::W
@ W
Definition: quadtree.h:63
j
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2

References Global_Variables::Dim, G, i, J, j, WallFunction::normal(), GlobalPhysicalVariables::ReInvFr, s, Global_Physical_Variables::Viscosity_Ratio, w, and oomph::QuadTreeNames::W.

◆ add_generic_residual_contribution() [2/2]

template<class ELEMENT >
void SpineGravityTractionElement< ELEMENT >::add_generic_residual_contribution ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian,
unsigned  flag 
)
inline

This function returns the residuals for the Navier–Stokes equations; flag=1(or 0): do (or don't) compute the Jacobian as well.

538  {
539  //Find out how many nodes there are
540  unsigned n_node = nnode();
541 
542  //Set the value of n_intpt
543  unsigned n_intpt = integral_pt()->nweight();
544 
545  //Set the Vector to hold local coordinates
546  Vector<double> s(Dim-1);
547 
548  //Set the Vector to hold the local coordinates of the parent element
549  Vector<double> s_parent(Dim);
550 
551  //Get a pointer to the parent element
552  ELEMENT* bulk_el_pt = dynamic_cast<ELEMENT*>(bulk_element_pt());
553 
554  //Find the number of nodes in the parent element
555  unsigned n_node_parent = bulk_el_pt->nnode();
556  //Set up memory for shape functions and their derivatives
557  Shape psif_parent(n_node_parent), testf_parent(n_node_parent);
558  DShape dpsifdx_parent(n_node_parent,Dim),
559  dtestfdx_parent(n_node_parent,Dim);
560 
561  //Storage for the local equation number
562  int local_eqn=0, local_unknown=0;
563 
564  //Get the Physical Variable
565  double ReInvFr = re_invfr()*density_ratio();
566  double Viscosity_Ratio = viscosity_ratio();
567  Vector<double> G = g();
568 
569 
570  unsigned n_external = this->nexternal_u_data();
571 
572  //Loop over the integration points
573  for(unsigned ipt=0;ipt<n_intpt;ipt++)
574  {
575  //Assign values of s and s_parent
576  for(unsigned i=0;i<(Dim-1);i++) {s[i] = integral_pt()->knot(ipt,i);}
577 
578  //Get the local coordinate in the bulk
579  this->get_local_coordinate_in_bulk(s,s_parent);
580 
581  //Get the integral weight
582  double w = integral_pt()->weight(ipt);
583 
584  //Find the shape functions and derivatives of the parent
585  (void)bulk_el_pt->
586  dshape_eulerian(s_parent,psif_parent,dpsifdx_parent);
587 
588  //Get the local jacobian for the FaceElement
589  double J = J_eulerian(s);
590 
591  //Premultiply the weights and the Jacobian
592  double W = w*J;
593 
594  //Need only to find the velocity derivatives
595  DenseMatrix<double> interpolated_dudx(Dim,Dim,0.0);
596 
597  Vector<double> interpolated_u(Dim,0.0);
598 
599  //Calculate velocities and derivatives
600  for(unsigned l=0;l<n_node_parent;l++)
601  {
602  //Loop over velocity components
603  for(unsigned i=0;i<Dim;i++)
604  {
605  for(unsigned j=0;j<Dim;j++)
606  {
607  interpolated_dudx(i,j) +=
608  bulk_el_pt->u_nst(l,i)*dpsifdx_parent(l,j);
609  }
610  dtestfdx_parent(l,i) = dpsifdx_parent(l,i);
611  interpolated_u[i] += bulk_el_pt->u_nst(l,i)*psif_parent(l);
612  }
613  //Set the test functions to be the same as the shape functions
614  testf_parent[l] = psif_parent[l];
615  }
616 
617  //Storage for the outer unit normal
618  Vector<double> normal(Dim);
619  outer_unit_normal(s,normal);
620 
621  //Loop over the test functions
622  for(unsigned l=0;l<n_node;l++)
623  {
624  //Do the first (axial) velocity component -- Poisson equation
625  {
626  unsigned i=1;
627 
628  local_eqn = U_local_eqn(l,i);
629  /*IF it's not a boundary condition*/
630  if(local_eqn >= 0)
631  {
632  //Add the gravitational body force term
633  residuals[local_eqn] +=
634  ReInvFr*testf_parent[bulk_node_number(l)]*G[i]*W;
635 
636  //Add in the equation term (equation constant traction in the axial direction d(Tyy+p)/dy =0)
637  residuals[local_eqn] -= Viscosity_Ratio*W* (
638  ( interpolated_dudx(i,0) /*+ interpolated_dudx(0,i)*/ )
639  *dtestfdx_parent(bulk_node_number(l),0) +
640  ( interpolated_dudx(i,2) /*+ interpolated_dudx(2,i) */)
641  *dtestfdx_parent(bulk_node_number(l),2) );
642 
643  //Now add the jacobian terms
644  if(flag)
645  {
646 
647  //Loop over all nodes again
648  for(unsigned l2=0;l2<n_node;l2++)
649  {
650 
651  /* {
652  unsigned i2=0;
653  local_unknown = U_local_eqn(l2,i2);
654  if(local_unknown >= 0)
655  {
656  jacobian(local_eqn,local_unknown) -=
657  Viscosity_Ratio*(dpsifdx_parent(bulk_node_number(l2),i)*
658  dtestfdx_parent(bulk_node_number(l),0))*W;
659  }
660  } */
661 
662  //Poisson equation term
663  {
664  unsigned i2=1;
665  local_unknown = U_local_eqn(l2,i2);
666  if(local_unknown >= 0)
667  {
668  jacobian(local_eqn,local_unknown) -=
669  Viscosity_Ratio*(dpsifdx_parent(bulk_node_number(l2),0)*
670  dtestfdx_parent(bulk_node_number(l),0) +
671  dpsifdx_parent(bulk_node_number(l2),2)*
672  dtestfdx_parent(bulk_node_number(l),2) )*W;
673  }
674  }
675 
676  /* {
677  unsigned i2=2;
678  local_unknown = U_local_eqn(l2,i2);
679  if(local_unknown >= 0)
680  {
681  jacobian(local_eqn,local_unknown) -=
682  Viscosity_Ratio*(dpsifdx_parent(bulk_node_number(l2),i)*
683  dtestfdx_parent(bulk_node_number(l),2))*W;
684  }
685  }*/
686 
687 
688  }
689 
690  //Loop over external data
691  for(unsigned l2=0;l2<n_external;l2++)
692  {
693  /* {
694  unsigned i2=0;
695  local_unknown = External_u_local_eqn(l2,i2);
696  if(local_unknown >= 0)
697  {
698  jacobian(local_eqn,local_unknown) -=
699  Viscosity_Ratio*(dpsifdx_parent(External_node[l2],i)*
700  dtestfdx_parent(bulk_node_number(l),0))*W;
701  }
702  }*/
703 
704  //Poisson equation term
705  {
706  unsigned i2=1;
707  local_unknown = External_u_local_eqn(l2,i2);
708  if(local_unknown >= 0)
709  {
710  jacobian(local_eqn,local_unknown) -=
711  Viscosity_Ratio*(dpsifdx_parent(External_node[l2],0)*
712  dtestfdx_parent(bulk_node_number(l),0) +
713  dpsifdx_parent(External_node[l2],2)*
714  dtestfdx_parent(bulk_node_number(l),2) )*W;
715  }
716  }
717 
718  /*{
719  unsigned i2=2;
720  local_unknown = External_u_local_eqn(l2,i2);
721  if(local_unknown >= 0)
722  {
723  jacobian(local_eqn,local_unknown) -=
724  Viscosity_Ratio*(dpsifdx_parent(External_node[l2],i)*
725  dtestfdx_parent(bulk_node_number(l),2))*W;
726  }
727  }*/
728  }
729  } //End of if flag
730  }
731 
732  }
733 
734  //Now do the other (traction) components
735  /* for(unsigned i =0;i<3;i+=2)
736  {
737  local_eqn = U_local_eqn(l,i);
738  //If it's not a boundary condition
739  if(local_eqn >= 0)
740  {
741  for(unsigned k=0;k<Dim;k++)
742  {
743  residuals[local_eqn] += Viscosity_Ratio*
744  (interpolated_dudx(i,k) + interpolated_dudx(k,i))*normal[k]
745  *testf_parent(bulk_node_number(l))*W;
746  }
747 
748  //Now add the jacobian terms
749  if(flag)
750  {
751  //Loop over all nodes again
752  for(unsigned l2=0;l2<n_node;l2++)
753  {
754  for(unsigned i2=0;i2<Dim;i2++)
755  {
756  if(i2!=i)
757  {
758  // unsigned i2=0;
759  local_unknown = U_local_eqn(l2,i2);
760  if(local_unknown >= 0)
761  {
762  jacobian(local_eqn,local_unknown) +=
763  Viscosity_Ratio*dpsifdx_parent(bulk_node_number(l2),i)*
764  normal[i2]*testf_parent(bulk_node_number(l))*W;
765  }
766  }
767  else
768  {
769  // unsigned i2=1;
770  local_unknown = U_local_eqn(l2,i2);
771  if(local_unknown >= 0)
772  {
773  for(unsigned k=0;k<Dim;k++)
774  {
775  jacobian(local_eqn,local_unknown) +=
776  Viscosity_Ratio*dpsifdx_parent(bulk_node_number(l2),k)*
777  normal[k]*testf_parent(bulk_node_number(l))*W;
778  }
779 
780  jacobian(local_eqn,local_unknown) +=
781  Viscosity_Ratio*dpsifdx_parent(bulk_node_number(l2),i)*
782  normal[i2]*testf_parent(bulk_node_number(l))*W;
783  }
784  }
785  }
786  }
787 
788 
789  //Loop over all external data
790  for(unsigned l2=0;l2<n_external;l2++)
791  {
792  for(unsigned i2=0;i2<Dim;i2++)
793  {
794  if (i2!=i)
795  {
796  //unsigned i2=0;
797  local_unknown = External_u_local_eqn(l2,i2);
798  if(local_unknown >= 0)
799  {
800  jacobian(local_eqn,local_unknown) +=
801  Viscosity_Ratio*dpsifdx_parent(External_node[l2],i)*
802  normal[i2]*testf_parent(bulk_node_number(l))*W;
803  }
804  }
805  else
806  {
807  //unsigned i2=1;
808  local_unknown = External_u_local_eqn(l2,i2);
809  if(local_unknown >= 0)
810  {
811  for(unsigned k=0;k<Dim;k++)
812  {
813  jacobian(local_eqn,local_unknown) +=
814  Viscosity_Ratio*dpsifdx_parent(External_node[l2],k)*
815  normal[k]*testf_parent(bulk_node_number(l))*W;
816  }
817 
818  jacobian(local_eqn,local_unknown) +=
819  Viscosity_Ratio*dpsifdx_parent(External_node[l2],i)*
820  normal[i2]*testf_parent(bulk_node_number(l))*W;
821  }
822  }
823  }
824  }
825  }
826  }
827  } */
828  }
829  }
830 
831  }

References Global_Variables::Dim, G, i, J, j, WallFunction::normal(), GlobalPhysicalVariables::ReInvFr, s, Global_Physical_Variables::Viscosity_Ratio, w, and oomph::QuadTreeNames::W.

◆ assign_additional_local_eqn_numbers() [1/2]

template<class ELEMENT >
void SpineGravityTractionElement< ELEMENT >::assign_additional_local_eqn_numbers ( )
inlinevirtual

Define the local equation numbering schemes.

Reimplemented from oomph::GeneralisedElement.

936  {
937  //Get number of nodes
938  unsigned n_node = nnode();
939  //Resize the equation counters
940  U_local_eqn.resize(n_node,Dim);
941 
942  //Loop over the nodes
943  for(unsigned i=0;i<n_node;i++)
944  {
945  //Loop over the nodal values
946  for(unsigned j=0;j<Dim;j++)
947  {
948  U_local_eqn(i,j) = this->nodal_local_eqn(i,j);
949  }
950  }
951 
952  //Now sort out the external degrees of freedom
953  //Find the number of external degrees of freedom
954  unsigned n_external_u = nexternal_u_data();
955  //Resize the external degree of freedom counter
956  External_u_local_eqn.resize(n_external_u,Dim);
957  //Loop over the external degrees of freedom
958  for(unsigned e=0;e<n_external_u;e++)
959  {
960  //Loop over the velocity values
961  for(unsigned i=0;i<Dim;i++)
962  {
963  External_u_local_eqn(e,i) = external_local_eqn(e,i);
964  }
965  }
966 
967  //Add the final additional data object for the external pressure
968  //gradient
969  if(Delta_P_pt!=0)
970  {
971  //This is the final external data object
972  Delta_P_local_eqn = external_local_eqn(External_Delta_P_index,0);
973  }
974  else
975  {
976  Delta_P_local_eqn = -1;
977  }
978  }
e
Array< double, 1, 3 > e(1./3., 0.5, 2.)
SpineGravityTractionElement::External_Delta_P_index
unsigned External_Delta_P_index
Definition: elastic_bretherton.cc:378
oomph::DenseMatrix::resize
void resize(const unsigned long &n)
Definition: matrices.h:498
oomph::FiniteElement::nodal_local_eqn
int nodal_local_eqn(const unsigned &n, const unsigned &i) const
Definition: elements.h:1432
oomph::GeneralisedElement::external_local_eqn
int external_local_eqn(const unsigned &i, const unsigned &j)
Definition: elements.h:311

References Global_Variables::Dim, e(), i, j, and oomph::DenseMatrix< T >::resize().

◆ assign_additional_local_eqn_numbers() [2/2]

template<class ELEMENT >
void SpineGravityTractionElement< ELEMENT >::assign_additional_local_eqn_numbers ( )
inlinevirtual

Define the local equation numbering schemes.

Reimplemented from oomph::GeneralisedElement.

836  {
837  //Get number of nodes
838  unsigned n_node = nnode();
839  //Resize the equation counters
840  U_local_eqn.resize(n_node,Dim);
841 
842  //Loop over the nodes
843  for(unsigned i=0;i<n_node;i++)
844  {
845  //Loop over the nodal values
846  for(unsigned j=0;j<Dim;j++)
847  {
848  U_local_eqn(i,j) = this->nodal_local_eqn(i,j);
849  }
850  }
851 
852  //Now sort out the external degrees of freedom
853  //Find the number of external degrees of freedom
854  //ACHTUNG !!!!!!
855 // unsigned n_external = nexternal_data();
856  unsigned n_external = External_node.size();
857 
858  //Resize the external degree of freedom counter
859  External_u_local_eqn.resize(n_external,Dim);
860  //Loop over the external degrees of freedom
861  for(unsigned e=0;e<n_external;e++)
862  {
863  //Loop over the velocity values
864  for(unsigned i=0;i<Dim;i++)
865  {
866  External_u_local_eqn(e,i) = external_local_eqn(e,i);
867  }
868  }
869  }

References Global_Variables::Dim, e(), i, j, and oomph::DenseMatrix< T >::resize().

◆ delta_p()

template<class ELEMENT >
const double SpineGravityTractionElement< ELEMENT >::delta_p ( ) const
inline

Return the value of the pressure gradient.

490  {
491  if(Delta_P_pt==0) {return 0.0;}
492  else {return Delta_P_pt->value(0);}
493  }
oomph::Data::value
double value(const unsigned &i) const
Definition: nodes.h:293

References oomph::Data::value().

◆ density_ratio() [1/2]

template<class ELEMENT >
const double& SpineGravityTractionElement< ELEMENT >::density_ratio ( ) const
inline

Density ratio for element: Element's density relative to the viscosity used in the definition of the Reynolds number

507 {return *Density_Ratio_pt;}

◆ density_ratio() [2/2]

template<class ELEMENT >
const double& SpineGravityTractionElement< ELEMENT >::density_ratio ( ) const
inline

Density ratio for element: Element's density relative to the viscosity used in the definition of the Reynolds number

403 {return *Density_Ratio_pt;}

◆ density_ratio_pt() [1/2]

template<class ELEMENT >
double* & SpineGravityTractionElement< ELEMENT >::density_ratio_pt ( )
inline

Pointer to Density ratio.

510 {return Density_Ratio_pt;}

◆ density_ratio_pt() [2/2]

template<class ELEMENT >
double* & SpineGravityTractionElement< ELEMENT >::density_ratio_pt ( )
inline

Pointer to Density ratio.

406 {return Density_Ratio_pt;}

◆ du_dt() [1/2]

template<class ELEMENT >
double SpineGravityTractionElement< ELEMENT >::du_dt ( const unsigned &  l,
const unsigned &  i 
) const
inline

i-th component of du/dt at local node l.

536  {
537  // Get the data's timestepper
538  TimeStepper* time_stepper_pt=node_pt(l)->time_stepper_pt();
539 
540  // Number of timsteps (past & present)
541  unsigned n_time = time_stepper_pt->ntstorage();
542 
543  double dudt=0.0;
544 
545  //Loop over the timesteps
546  if (time_stepper_pt->type()!="Steady")
547  {
548  for(unsigned t=0;t<n_time;t++)
549  {
550  dudt+=time_stepper_pt->weight(1,t)*u(t,l,i);
551  }
552  }
553 
554  return dudt;
555  }
SpineGravityTractionElement::u
double u(const unsigned &l, const unsigned &i)
Access function for the velocity. N. B. HEAVY ASSUMPTIONS HERE.
Definition: elastic_bretherton.cc:525
oomph::Data::time_stepper_pt
TimeStepper *& time_stepper_pt()
Return the pointer to the timestepper.
Definition: nodes.h:238
oomph::GeomObject::time_stepper_pt
TimeStepper *& time_stepper_pt()
Definition: geom_objects.h:192
oomph::TimeStepper
Definition: timesteppers.h:231
oomph::TimeStepper::ntstorage
unsigned ntstorage() const
Definition: timesteppers.h:601
oomph::TimeStepper::weight
virtual double weight(const unsigned &i, const unsigned &j) const
Access function for j-th weight for the i-th derivative.
Definition: timesteppers.h:594
oomph::TimeStepper::type
std::string type() const
Definition: timesteppers.h:490
plotPSD.t
t
Definition: plotPSD.py:36

References i, oomph::TimeStepper::ntstorage(), plotPSD::t, oomph::TimeStepper::type(), and oomph::TimeStepper::weight().

◆ du_dt() [2/2]

template<class ELEMENT >
double SpineGravityTractionElement< ELEMENT >::du_dt ( const unsigned &  l,
const unsigned &  i 
) const
inline

i-th component of du/dt at local node l.

488  {
489  // Get the data's timestepper
490  TimeStepper* time_stepper_pt=node_pt(l)->time_stepper_pt(); // cgj: previously Node_pt[l] -- but (interpreting as FiniteElement::Node_pt as clang does) this is private!
491 
492  // Number of timsteps (past & present)
493  unsigned n_time = time_stepper_pt->ntstorage();
494 
495  double dudt=0.0;
496 
497  //Loop over the timesteps
498  if (time_stepper_pt->type()!="Steady")
499  {
500  for(unsigned t=0;t<n_time;t++)
501  {
502  dudt+=time_stepper_pt->weight(1,t)*u(t,l,i);
503  }
504  }
505 
506  return dudt;
507  }

References i, oomph::TimeStepper::ntstorage(), plotPSD::t, oomph::TimeStepper::type(), and oomph::TimeStepper::weight().

◆ fill_in_contribution_to_jacobian() [1/2]

template<class ELEMENT >
void SpineGravityTractionElement< ELEMENT >::fill_in_contribution_to_jacobian ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian 
)
inlinevirtual

This function returns the residuals and the jacobian.

Reimplemented from oomph::GeneralisedElement.

569  {
570  //Call the generic routine with the flag set to 1
571  add_generic_residual_contribution(residuals,jacobian,1);
572  //Add in external data contributions
573  this->fill_in_jacobian_from_external_by_fd(jacobian);
574  //Add the spine contributions
575  this->fill_in_jacobian_from_geometric_data(jacobian);
576  }
SpineGravityTractionElement::add_generic_residual_contribution
void add_generic_residual_contribution(Vector< double > &residuals, DenseMatrix< double > &jacobian, unsigned flag)
Definition: elastic_bretherton.cc:583
oomph::ElementWithMovingNodes::fill_in_jacobian_from_geometric_data
void fill_in_jacobian_from_geometric_data(Vector< double > &residuals, DenseMatrix< double > &jacobian)
Definition: element_with_moving_nodes.cc:323
oomph::GeneralisedElement::fill_in_jacobian_from_external_by_fd
void fill_in_jacobian_from_external_by_fd(Vector< double > &residuals, DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
Definition: elements.cc:1199

◆ fill_in_contribution_to_jacobian() [2/2]

template<class ELEMENT >
void SpineGravityTractionElement< ELEMENT >::fill_in_contribution_to_jacobian ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian 
)
inlinevirtual

This function returns the residuals and the jacobian.

Reimplemented from oomph::GeneralisedElement.

521  {
522  //Call the generic routine with the flag set to 1
523  add_generic_residual_contribution(residuals,jacobian,1);
524  //Add the external data
525  this->fill_in_jacobian_from_external_by_fd(jacobian);
526  //Add the spine contributions
527  this->fill_in_jacobian_from_geometric_data(jacobian);
528  }

◆ fill_in_contribution_to_residuals() [1/2]

template<class ELEMENT >
void SpineGravityTractionElement< ELEMENT >::fill_in_contribution_to_residuals ( Vector< double > &  residuals)
inlinevirtual

Add the contribution to the residuals.

Reimplemented from oomph::GeneralisedElement.

559  {
560  //Create a dummy matrix
561  DenseMatrix<double> dummy(1);
562  //Call the generic residuals function with flag set to 0
563  add_generic_residual_contribution(residuals,dummy,0);
564  }

◆ fill_in_contribution_to_residuals() [2/2]

template<class ELEMENT >
void SpineGravityTractionElement< ELEMENT >::fill_in_contribution_to_residuals ( Vector< double > &  residuals)
inlinevirtual

Add the contribution to the residuals.

Reimplemented from oomph::GeneralisedElement.

511  {
512  //Create a dummy matrix
513  DenseMatrix<double> dummy(1);
514  //Call the generic residuals function with flag set to 0
515  add_generic_residual_contribution(residuals,dummy,0);
516  }

◆ flow()

template<class ELEMENT >
double SpineGravityTractionElement< ELEMENT >::flow ( )
inline

Calculate the flow across the element.

422  {
423  //Storage for the normal flux
424  double flux = 0.0;
425 
426  //Set the value of n_intpt
427  const unsigned n_intpt = integral_pt()->nweight();
428 
429  //Set the Vector to hold local coordinates
430  Vector<double> s(Dim-1);
431 
432  //Set the Vector to hold the local coordinates of the parent element
433  Vector<double> s_parent(Dim);
434 
435  //Get a pointer to the parent element
436  ELEMENT* bulk_el_pt = dynamic_cast<ELEMENT*>(bulk_element_pt());
437 
438  //Loop over the integration points
439  for(unsigned ipt=0;ipt<n_intpt;ipt++)
440  {
441  //Assign values of s and s_parent
442  for(unsigned i=0;i<(Dim-1);i++) {s[i] = integral_pt()->knot(ipt,i);}
443 
444  //Get the local coordinate in the bulk
445  this->get_local_coordinate_in_bulk(s,s_parent);
446 
447  //Get the integral weight
448  double w = integral_pt()->weight(ipt);
449 
450  //Get the local jacobian for the FaceElement
451  double J = J_eulerian(s);
452 
453  //Premultiply the weights and the Jacobian
454  double W = w*J;
455 
456  //Get the velocity from the parent
457  Vector<double> interpolated_u(Dim,0.0);
458  bulk_el_pt->interpolated_u_nst(s_parent,interpolated_u);
459 
460  //Get the outer unit normal
461  Vector<double> normal(Dim,0);
462  outer_unit_normal(s,normal);
463 
464  //Find the normal flux
465  double normal_flux = 0.0;
466  for(unsigned i=0;i<Dim;i++) {normal_flux += interpolated_u[i]*normal[i];}
467 
468  //Add to the flux
469  flux += normal_flux*W;
470 
471  } //End of loop over integration points
472  //Return the flux
473  return flux;
474  }
ProblemParameters::flux
void flux(const double &time, const Vector< double > &x, double &flux)
Get flux applied along boundary x=0.
Definition: pretend_melt.cc:59

References Global_Variables::Dim, ProblemParameters::flux(), i, J, WallFunction::normal(), s, w, and oomph::QuadTreeNames::W.

◆ g() [1/2]

template<class ELEMENT >
const Vector<double>& SpineGravityTractionElement< ELEMENT >::g ( ) const
inline

Vector of gravitational components.

519 {return *G_pt;}

◆ g() [2/2]

template<class ELEMENT >
const Vector<double>& SpineGravityTractionElement< ELEMENT >::g ( ) const
inline

Vector of gravitational components.

415 {return *G_pt;}

◆ g_pt() [1/2]

template<class ELEMENT >
Vector<double>* & SpineGravityTractionElement< ELEMENT >::g_pt ( )
inline

Pointer to Vector of gravitational components.

522 {return G_pt;}

Referenced by AirwayReopeningProblem< ELEMENT >::AirwayReopeningProblem(), and ThreeDimBethertonProblem< ELEMENT >::ThreeDimBethertonProblem().

◆ g_pt() [2/2]

template<class ELEMENT >
Vector<double>* & SpineGravityTractionElement< ELEMENT >::g_pt ( )
inline

Pointer to Vector of gravitational components.

418 {return G_pt;}

◆ get_flux()

template<class ELEMENT >
double SpineGravityTractionElement< ELEMENT >::get_flux ( )
inline

Calculate the flux.

879  {
880  //Initialise flux to zero
881  double flux = 0.0;
882 
883  //Find out how many nodes there are
884  unsigned n_node = nnode();
885 
886  //Set the value of n_intpt
887  unsigned n_intpt = integral_pt()->nweight();
888 
889  //Set the Vector to hold local coordinates
890  Vector<double> s(Dim-1);
891 
892  //Set up memory for shape functions and their derivatives
893  Shape psif(n_node);
894 
895  //Loop over the integration points
896  for(unsigned ipt=0;ipt<n_intpt;ipt++)
897  {
898  //Assign values of s and s_parent
899  for(unsigned i=0;i<(Dim-1);i++)
900  {s[i] = integral_pt()->knot(ipt,i);}
901 
902  //Get the integral weight
903  double w = integral_pt()->weight(ipt);
904 
905  //Find the shape functions and derivatives of the parent
906  shape(s,psif);
907 
908  //Get the local jacobian for the FaceElement
909  double J = J_eulerian(s);
910 
911  //Premultiply the weights and the Jacobian
912  double W = w*J;
913 
914  Vector<double> interpolated_u(Dim,0.0);
915 
916  //Calculate velocities and derivatives
917  for(unsigned l=0;l<n_node;l++)
918  {
919  //Loop over velocity components
920  for(unsigned i=0;i<Dim;i++)
921  {
922  interpolated_u[i] += u(l,i)*psif(l);
923  }
924  }
925 
926 
927  flux += interpolated_u[0]*W;
928  }
929 
930  return flux;
931  }
oomph::FiniteElement::shape
virtual void shape(const Vector< double > &s, Shape &psi) const =0

References Global_Variables::Dim, ProblemParameters::flux(), i, J, s, oomph::OneDimLagrange::shape(), w, and oomph::QuadTreeNames::W.

◆ hijack_all_nodes()

template<class ELEMENT >
void SpineGravityTractionElement< ELEMENT >::hijack_all_nodes ( )
inline
470  {
471  ELEMENT* cast_element_pt = dynamic_cast<ELEMENT*>(bulk_element_pt());
472  //Hijack the nodes in the bulk element in the axial coordinate
473  unsigned n_node = this->nnode();
474  for(unsigned m=0;m<n_node;m++)
475  {
476  delete cast_element_pt->hijack_nodal_value(bulk_node_number(m),0);
477  }
478  }

References m.

◆ invca_local_eqn()

template<class ELEMENT >
int SpineGravityTractionElement< ELEMENT >::invca_local_eqn ( )
inlineprotected
291 {
292  return external_local_eqn(External_data_number_of_invca,0);
293 
294 }
SpineGravityTractionElement::External_data_number_of_invca
unsigned External_data_number_of_invca
Definition: extra_elements.h:280

◆ nexternal_u_data() [1/2]

template<class ELEMENT >
unsigned SpineGravityTractionElement< ELEMENT >::nexternal_u_data ( )
inline

Return the number of external velocity data.

496 {return External_node.size();}

◆ nexternal_u_data() [2/2]

template<class ELEMENT >
unsigned SpineGravityTractionElement< ELEMENT >::nexternal_u_data ( )
inline

Return the number of external velocity data.

392 {return External_node.size();}

◆ node_update()

template<class ELEMENT >
void SpineGravityTractionElement< ELEMENT >::node_update ( )
inlinevirtual

Update the positions of all nodes in the element using each node update function. The default implementation may be overloaded so that more efficient versions can be written

Loop over all nodes in the element and update their positions using each node's (algebraic) update function

Reimplemented from oomph::FiniteElement.

386  {
387  this->bulk_element_pt()->node_update();
388  }
oomph::FiniteElement::node_update
virtual void node_update()
Definition: elements.cc:5072

◆ output() [1/4]

template<class ELEMENT >
void SpineGravityTractionElement< ELEMENT >::output ( ostream &  outfile)
inlinevirtual

Overload the output function.

Reimplemented from oomph::FiniteElement.

981 { }

◆ output() [2/4]

template<class ELEMENT >
void SpineGravityTractionElement< ELEMENT >::output ( ostream &  outfile,
const unsigned &  Np 
)
inlinevirtual

Output function: x,y,[z],u,v,[w],p in tecplot format.

Reimplemented from oomph::FiniteElement.

984 { }

◆ output() [3/4]

template<class ELEMENT >
void SpineGravityTractionElement< ELEMENT >::output ( std::ostream &  outfile)
inlinevirtual

Overload the output function.

Reimplemented from oomph::FiniteElement.

872 { }

◆ output() [4/4]

template<class ELEMENT >
void SpineGravityTractionElement< ELEMENT >::output ( std::ostream &  outfile,
const unsigned &  Np 
)
inlinevirtual

Output function: x,y,[z],u,v,[w],p in tecplot format.

Reimplemented from oomph::FiniteElement.

876  {
877  SpineElement<FaceGeometry<ELEMENT> >::output(outfile,Np);
878  }
SpineGravityTractionElement::output
void output(ostream &outfile)
Overload the output function.
Definition: elastic_bretherton.cc:981

References oomph::output().

◆ re_invfr() [1/2]

template<class ELEMENT >
const double& SpineGravityTractionElement< ELEMENT >::re_invfr ( )
inline

Reynolds number divided by Froude number.

516 {return *ReInvFr_pt;}

◆ re_invfr() [2/2]

template<class ELEMENT >
const double& SpineGravityTractionElement< ELEMENT >::re_invfr ( ) const
inline

Return the value of the Re/Fr number.

412 {return *ReInvFr_pt;}

◆ re_invfr_pt() [1/2]

template<class ELEMENT >
double* & SpineGravityTractionElement< ELEMENT >::re_invfr_pt ( )
inline

Pointer to Reynolds number divided by Froude number.

513 {return ReInvFr_pt;}

Referenced by AirwayReopeningProblem< ELEMENT >::AirwayReopeningProblem(), and ThreeDimBethertonProblem< ELEMENT >::ThreeDimBethertonProblem().

◆ re_invfr_pt() [2/2]

template<class ELEMENT >
double* & SpineGravityTractionElement< ELEMENT >::re_invfr_pt ( )
inline

Pointer to Reynolds number divided by Froude number.

409 {return ReInvFr_pt;}

◆ set_delta_p_pt()

template<class ELEMENT >
void SpineGravityTractionElement< ELEMENT >::set_delta_p_pt ( Data *const &  delta_p_pt)
inline

Access function to the pointer to pressure gradient data.

482  {
483  Delta_P_pt = delta_p_pt;
484  //Add the external data and do not finite difference it!
485  External_Delta_P_index = add_external_data(delta_p_pt,false);
486  }

Referenced by AirwayReopeningProblem< ELEMENT >::construct_poisson_outlet_elements().

◆ u() [1/4]

template<class ELEMENT >
double SpineGravityTractionElement< ELEMENT >::u ( const unsigned &  l,
const unsigned &  i 
)
inline

Access function for the velocity. N. B. HEAVY ASSUMPTIONS HERE.

526  {return nodal_value(l,i);}
oomph::FiniteElement::nodal_value
double nodal_value(const unsigned &n, const unsigned &i) const
Definition: elements.h:2593

References i.

◆ u() [2/4]

template<class ELEMENT >
double SpineGravityTractionElement< ELEMENT >::u ( const unsigned &  l,
const unsigned &  i 
)
inline

Access function for the velocity. N. B. HEAVY ASSUMPTIONS HERE.

478  {return this->nodal_value(l,i);}

References i.

◆ u() [3/4]

template<class ELEMENT >
double SpineGravityTractionElement< ELEMENT >::u ( const unsigned &  t,
const unsigned &  l,
const unsigned &  i 
) const
inline

/// Velocity i at local node l at timestep t (t=0: present; t>0: previous). SIMILAR HEAVY ASSUMPTIONS

532  {return nodal_value(t,l,i);}

References i, and plotPSD::t.

◆ u() [4/4]

template<class ELEMENT >
double SpineGravityTractionElement< ELEMENT >::u ( const unsigned &  t,
const unsigned &  l,
const unsigned &  i 
) const
inline

/// Velocity i at local node l at timestep t (t=0: present; t>0: previous). SIMILAR HEAVY ASSUMPTIONS

484  {return this->nodal_value(t,l,i);}

References i, and plotPSD::t.

◆ viscosity_ratio() [1/2]

template<class ELEMENT >
const double& SpineGravityTractionElement< ELEMENT >::viscosity_ratio ( ) const
inline

Viscosity ratio for element: Element's viscosity relative to the viscosity used in the definition of the Reynolds number

500 {return *Viscosity_Ratio_pt;}

◆ viscosity_ratio() [2/2]

template<class ELEMENT >
const double& SpineGravityTractionElement< ELEMENT >::viscosity_ratio ( ) const
inline

Viscosity ratio for element: Element's viscosity relative to the viscosity used in the definition of the Reynolds number

396 {return *Viscosity_Ratio_pt;}

◆ viscosity_ratio_pt() [1/2]

template<class ELEMENT >
double* & SpineGravityTractionElement< ELEMENT >::viscosity_ratio_pt ( )
inline

Pointer to Viscosity Ratio.

503 {return Viscosity_Ratio_pt;}

◆ viscosity_ratio_pt() [2/2]

template<class ELEMENT >
double* & SpineGravityTractionElement< ELEMENT >::viscosity_ratio_pt ( )
inline

Pointer to Viscosity Ratio.

399 {return Viscosity_Ratio_pt;}

Member Data Documentation

◆ bond_pt

template<class ELEMENT >
double* SpineGravityTractionElement< ELEMENT >::bond_pt
protected

◆ Delta_P_local_eqn

template<class ELEMENT >
int SpineGravityTractionElement< ELEMENT >::Delta_P_local_eqn
protected

The local equation number for the external data associated with the unknown pressure gradient

◆ Delta_P_pt

template<class ELEMENT >
Data* SpineGravityTractionElement< ELEMENT >::Delta_P_pt
private

Pointer to an External Data object that represents the an unknown pressure gradient

◆ Density_Ratio_pt

template<class ELEMENT >
double * SpineGravityTractionElement< ELEMENT >::Density_Ratio_pt
private

Pointer to the density ratio (relative to the density used in the definition of the Reynolds number)

◆ Dim

template<class ELEMENT >
unsigned SpineGravityTractionElement< ELEMENT >::Dim
private

The highest dimension of the problem.

◆ External_data_number_of_invca

template<class ELEMENT >
unsigned SpineGravityTractionElement< ELEMENT >::External_data_number_of_invca
protected

The Data that contains the capillary number is stored as external Data for the element. Which external Data item is it?

◆ External_Delta_P_index

template<class ELEMENT >
unsigned SpineGravityTractionElement< ELEMENT >::External_Delta_P_index
protected

The index in the external data at which the Delta_p data is stored

◆ External_node

template<class ELEMENT >
Vector< unsigned > SpineGravityTractionElement< ELEMENT >::External_node
protected

Vector to keep track of the external data associated with each bulk node

◆ External_u_local_eqn

template<class ELEMENT >
DenseMatrix< int > SpineGravityTractionElement< ELEMENT >::External_u_local_eqn
protected

Array to hold the local eqn number information for the external data (other nodes in the bulk element)

◆ G_pt

template<class ELEMENT >
Vector< double > * SpineGravityTractionElement< ELEMENT >::G_pt
private

Pointer to global gravity Vector.

◆ invca_data_pt

template<class ELEMENT >
Data* SpineGravityTractionElement< ELEMENT >::invca_data_pt
protected

Pointer to the Data item that stores the capillary number.

◆ ReInvFr_pt

template<class ELEMENT >
double * SpineGravityTractionElement< ELEMENT >::ReInvFr_pt
private

Pointer to the global Reynold number divided by the Froude number.

◆ U_local_eqn

template<class ELEMENT >
DenseMatrix< int > SpineGravityTractionElement< ELEMENT >::U_local_eqn
protected

Array to hold local eqn number information for veloc: U_local_eqn(jnod,i) = local equation number or < 0 if pinned

◆ Viscosity_Ratio_pt

template<class ELEMENT >
double * SpineGravityTractionElement< ELEMENT >::Viscosity_Ratio_pt
private

Pointer to the viscosity ratio (relative to the viscosity used in the definition of the Reynolds number)

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