Public Member Functions | Public Attributes | Protected Member Functions | Protected Attributes | List of all members
oomph::AxisymmetricNavierStokesTractionElement< ELEMENT > Class Template Reference

#include <axisym_fluid_traction_elements.h>

+ Inheritance diagram for oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >:

Public Member Functions

 AxisymmetricNavierStokesTractionElement (FiniteElement *const &element_pt, const int &face_index)
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Return the residuals. More...
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 Fill in contribution from Jacobian. More...
 
double zeta_nodal (const unsigned &n, const unsigned &k, const unsigned &i) const
 
void output (std::ostream &outfile)
 Output function. More...
 
unsigned nscalar_paraview () const
 
void scalar_value_paraview (std::ofstream &file_out, const unsigned &k, const unsigned &nplot) const
 
std::string scalar_name_paraview (const unsigned &i) const
 
void output (std::ostream &outfile, const unsigned &n_plot)
 Output function. More...
 
void output (FILE *file_pt)
 C_style output function. More...
 
void output (FILE *file_pt, const unsigned &n_plot)
 C-style output function. More...
 
void traction (const double &time, const Vector< double > &s, Vector< double > &traction)
 
- 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...
 
- 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_local_dofs (std::ostream &out, const std::string &current_string) const
 
virtual void describe_nodal_local_dofs (std::ostream &out, const std::string &current_string) const
 
virtual void assign_all_generic_local_eqn_numbers (const bool &store_local_dof_pt)
 
Node *& node_pt (const unsigned &n)
 Return a pointer to the local node n. More...
 
Node *const & node_pt (const unsigned &n) const
 Return a pointer to the local node n (const version) More...
 
unsigned nnode () const
 Return the number of nodes. More...
 
virtual unsigned nnode_1d () const
 
double raw_nodal_position (const unsigned &n, const unsigned &i) const
 
double raw_nodal_position (const unsigned &t, const unsigned &n, const unsigned &i) const
 
double raw_dnodal_position_dt (const unsigned &n, const unsigned &i) const
 
double raw_dnodal_position_dt (const unsigned &n, const unsigned &j, const unsigned &i) const
 
double raw_nodal_position_gen (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double raw_nodal_position_gen (const unsigned &t, const unsigned &n, const unsigned &k, const unsigned &i) const
 
double raw_dnodal_position_gen_dt (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double raw_dnodal_position_gen_dt (const unsigned &j, const unsigned &n, const unsigned &k, const unsigned &i) const
 
double nodal_position (const unsigned &n, const unsigned &i) const
 
double nodal_position (const unsigned &t, const unsigned &n, const unsigned &i) const
 
double dnodal_position_dt (const unsigned &n, const unsigned &i) const
 Return the i-th component of nodal velocity: dx/dt at local node n. More...
 
double dnodal_position_dt (const unsigned &n, const unsigned &j, const unsigned &i) const
 
double nodal_position_gen (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double nodal_position_gen (const unsigned &t, const unsigned &n, const unsigned &k, const unsigned &i) const
 
double dnodal_position_gen_dt (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double dnodal_position_gen_dt (const unsigned &j, const unsigned &n, const unsigned &k, const unsigned &i) const
 
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
 
virtual Node * construct_node (const unsigned &n)
 
virtual Node * construct_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 
virtual Node * construct_boundary_node (const unsigned &n)
 
virtual Node * construct_boundary_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 
int get_node_number (Node *const &node_pt) const
 
virtual Node * get_node_at_local_coordinate (const Vector< double > &s) const
 
double raw_nodal_value (const unsigned &n, const unsigned &i) const
 
double raw_nodal_value (const unsigned &t, const unsigned &n, const unsigned &i) const
 
double nodal_value (const unsigned &n, const unsigned &i) const
 
double nodal_value (const unsigned &t, const unsigned &n, const unsigned &i) const
 
unsigned dim () const
 
virtual ElementGeometry::ElementGeometry element_geometry () const
 Return the geometry type of the element (either Q or T usually). More...
 
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 node_update ()
 
virtual void identify_field_data_for_interactions (std::set< std::pair< Data *, unsigned >> &paired_field_data)
 
virtual void identify_geometric_data (std::set< Data * > &geometric_data_pt)
 
virtual double s_min () const
 Min value of local coordinate. More...
 
virtual double s_max () const
 Max. value of local coordinate. More...
 
double size () const
 
virtual double compute_physical_size () const
 
virtual void point_output_data (const Vector< double > &s, Vector< double > &data)
 
void point_output (std::ostream &outfile, const Vector< double > &s)
 
virtual unsigned nplot_points_paraview (const unsigned &nplot) const
 
virtual unsigned nsub_elements_paraview (const unsigned &nplot) const
 
void output_paraview (std::ofstream &file_out, const unsigned &nplot) const
 
virtual void write_paraview_output_offset_information (std::ofstream &file_out, const unsigned &nplot, unsigned &counter) const
 
virtual void write_paraview_type (std::ofstream &file_out, const unsigned &nplot) const
 
virtual void write_paraview_offsets (std::ofstream &file_out, const unsigned &nplot, unsigned &offset_sum) const
 
virtual void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt) const
 
virtual void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt) const
 
virtual void output (const unsigned &t, std::ostream &outfile, const unsigned &n_plot) const
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 Output an exact solution over the element. More...
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt)
 Output a time-dependent exact solution over the element. More...
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, const SolutionFunctorBase &exact_soln) const
 Output a time-dependent exact solution over the element. More...
 
virtual void get_s_plot (const unsigned &i, const unsigned &nplot, Vector< double > &s, const bool &shifted_to_interior=false) const
 
virtual std::string tecplot_zone_string (const unsigned &nplot) const
 
virtual void write_tecplot_zone_footer (std::ostream &outfile, const unsigned &nplot) const
 
virtual void write_tecplot_zone_footer (FILE *file_pt, const unsigned &nplot) const
 
virtual unsigned nplot_points (const unsigned &nplot) const
 
virtual void compute_error (FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 Calculate the norm of the error and that of the exact solution. More...
 
virtual void compute_error (FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 Calculate the norm of the error and that of the exact solution. More...
 
virtual void compute_error (FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_abs_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error)
 
void integrate_fct (FiniteElement::SteadyExactSolutionFctPt integrand_fct_pt, Vector< double > &integral)
 Integrate Vector-valued function over element. More...
 
void integrate_fct (FiniteElement::UnsteadyExactSolutionFctPt integrand_fct_pt, const double &time, Vector< double > &integral)
 Integrate Vector-valued time-dep function over element. More...
 
virtual void build_face_element (const int &face_index, FaceElement *face_element_pt)
 
virtual unsigned self_test ()
 
virtual unsigned get_bulk_node_number (const int &face_index, const unsigned &i) const
 
virtual int face_outer_unit_normal_sign (const int &face_index) const
 Get the sign of the outer unit normal on the face given by face_index. More...
 
virtual unsigned nnode_on_face () const
 
void face_node_number_error_check (const unsigned &i) const
 Range check for face node numbers. More...
 
virtual CoordinateMappingFctPt face_to_bulk_coordinate_fct_pt (const int &face_index) const
 
virtual BulkCoordinateDerivativesFctPt bulk_coordinate_derivatives_fct_pt (const int &face_index) const
 
- Public Member Functions inherited from oomph::GeneralisedElement
 GeneralisedElement ()
 Constructor: Initialise all pointers and all values to zero. More...
 
virtual ~GeneralisedElement ()
 Virtual destructor to clean up any memory allocated by the object. More...
 
 GeneralisedElement (const GeneralisedElement &)=delete
 Broken copy constructor. More...
 
void operator= (const GeneralisedElement &)=delete
 Broken assignment operator. More...
 
Data *& internal_data_pt (const unsigned &i)
 Return a pointer to i-th internal data object. More...
 
Data *const & internal_data_pt (const unsigned &i) const
 Return a pointer to i-th internal data object (const version) More...
 
Data *& external_data_pt (const unsigned &i)
 Return a pointer to i-th external data object. More...
 
Data *const & external_data_pt (const unsigned &i) const
 Return a pointer to i-th external data object (const version) More...
 
unsigned long eqn_number (const unsigned &ieqn_local) const
 
int local_eqn_number (const unsigned long &ieqn_global) const
 
unsigned add_external_data (Data *const &data_pt, const bool &fd=true)
 
bool external_data_fd (const unsigned &i) const
 
void exclude_external_data_fd (const unsigned &i)
 
void include_external_data_fd (const unsigned &i)
 
void flush_external_data ()
 Flush all external data. More...
 
void flush_external_data (Data *const &data_pt)
 Flush the object addressed by data_pt from the external data array. More...
 
unsigned ninternal_data () const
 Return the number of internal data objects. More...
 
unsigned nexternal_data () const
 Return the number of external data objects. More...
 
unsigned ndof () const
 Return the number of equations/dofs in the element. More...
 
void dof_vector (const unsigned &t, Vector< double > &dof)
 Return the vector of dof values at time level t. More...
 
void dof_pt_vector (Vector< double * > &dof_pt)
 Return the vector of pointers to dof values. More...
 
void set_internal_data_time_stepper (const unsigned &i, TimeStepper *const &time_stepper_pt, const bool &preserve_existing_data)
 
void assign_internal_eqn_numbers (unsigned long &global_number, Vector< double * > &Dof_pt)
 
void describe_dofs (std::ostream &out, const std::string &current_string) const
 
void add_internal_value_pt_to_map (std::map< unsigned, double * > &map_of_value_pt)
 
virtual void assign_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual void complete_setup_of_dependencies ()
 
virtual void get_residuals (Vector< double > &residuals)
 
virtual void get_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
virtual void get_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &mass_matrix)
 
virtual void get_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
virtual void get_dresiduals_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam)
 
virtual void get_djacobian_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
 
virtual void get_djacobian_and_dmass_matrix_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam, DenseMatrix< double > &dmass_matrix_dparam)
 
virtual void get_hessian_vector_products (Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
 
virtual void get_inner_products (Vector< std::pair< unsigned, unsigned >> const &history_index, Vector< double > &inner_product)
 
virtual void get_inner_product_vectors (Vector< unsigned > const &history_index, Vector< Vector< double >> &inner_product_vector)
 
virtual void compute_norm (Vector< double > &norm)
 
virtual void compute_norm (double &norm)
 
virtual 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 Attributes

void(*&)(const double &time, const Vector< double > &x, const Vector< double > &n, Vector< double > &traction) traction_fct_pt ()
 Reference to the traction function pointer. More...
 

Protected Member Functions

virtual void get_traction (const double &time, const unsigned &intpt, const Vector< double > &x, const Vector< double > &n, Vector< double > &traction) const
 
void fill_in_contribution_to_residuals_axisymmetric_nst_traction (Vector< double > &residuals)
 Helper function that actually calculates the residuals. More...
 
- Protected Member Functions inherited from oomph::FaceElement
void add_additional_values (const Vector< unsigned > &nadditional_values, const unsigned &id)
 
- 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)
 
virtual void assign_additional_local_eqn_numbers ()
 
int internal_local_eqn (const unsigned &i, const unsigned &j) const
 
int external_local_eqn (const unsigned &i, const unsigned &j)
 
void fill_in_jacobian_from_internal_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_internal_by_fd (DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_external_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_external_by_fd (DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
virtual void update_before_internal_fd ()
 
virtual void reset_after_internal_fd ()
 
virtual void update_in_internal_fd (const unsigned &i)
 
virtual void reset_in_internal_fd (const unsigned &i)
 
virtual void update_before_external_fd ()
 
virtual void reset_after_external_fd ()
 
virtual void update_in_external_fd (const unsigned &i)
 
virtual void reset_in_external_fd (const unsigned &i)
 
virtual void fill_in_contribution_to_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &mass_matrix)
 
virtual void fill_in_contribution_to_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 Attributes

Vector< unsigned > U_index_axisymmetric_nst_traction
 Index at which the i-th velocity component is stored. More...
 
void(* Traction_fct_pt )(const double &time, const Vector< double > &x, const Vector< double > &n, Vector< double > &result)
 
- 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
 
- 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
 

Additional Inherited Members

- Public Types inherited from oomph::FiniteElement
typedef void(* SteadyExactSolutionFctPt) (const Vector< double > &, Vector< double > &)
 
typedef void(* UnsteadyExactSolutionFctPt) (const double &, const Vector< double > &, Vector< double > &)
 
- 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 oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >

A class for elements that allow the imposition of an applied traction in the axisym Navier Stokes eqns. The geometrical information can be read from the FaceGeometry<ELEMENT> class and and thus, we can be generic enough without the need to have a separate equations class.

Constructor & Destructor Documentation

◆ AxisymmetricNavierStokesTractionElement()

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

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

118  : FaceGeometry<ELEMENT>(), FaceElement()
119  {
120  // Attach the geometrical information to the element. N.B. This function
121  // also assigns nbulk_value from the required_nvalue of the bulk element
122  element_pt->build_face_element(face_index, this);
123 
124  // Find the dimension of the problem
125  unsigned n_dim = element_pt->nodal_dimension();
126 
127  // Find the index at which the velocity unknowns are stored
128  ELEMENT* cast_element_pt = dynamic_cast<ELEMENT*>(element_pt);
129  this->U_index_axisymmetric_nst_traction.resize(n_dim + 1);
130  for (unsigned i = 0; i < n_dim + 1; i++)
131  {
132  this->U_index_axisymmetric_nst_traction[i] =
133  cast_element_pt->u_index_axi_nst(i);
134  }
135 
136  // Zero traction
137  Traction_fct_pt =
138  &AxisymmetricNavierStokesTractionElementHelper::Zero_traction_fct;
139  }
i
int i
Definition: BiCGSTAB_step_by_step.cpp:9
oomph::AxisymmetricNavierStokesTractionElement::U_index_axisymmetric_nst_traction
Vector< unsigned > U_index_axisymmetric_nst_traction
Index at which the i-th velocity component is stored.
Definition: axisym_fluid_traction_elements.h:77
oomph::AxisymmetricNavierStokesTractionElement::Traction_fct_pt
void(* Traction_fct_pt)(const double &time, const Vector< double > &x, const Vector< double > &n, Vector< double > &result)
Definition: axisym_fluid_traction_elements.h:84
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::AxisymmetricNavierStokesTractionElementHelper::Zero_traction_fct
void Zero_traction_fct(const double &time, const Vector< double > &x, const Vector< double > &N, Vector< double > &load)
Default load function (zero traction)
Definition: axisym_fluid_traction_elements.cc:47

References oomph::FiniteElement::build_face_element(), oomph::FaceElement::face_index(), i, oomph::FiniteElement::nodal_dimension(), oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::Traction_fct_pt, and oomph::AxisymmetricNavierStokesTractionElementHelper::Zero_traction_fct().

Member Function Documentation

◆ fill_in_contribution_to_jacobian()

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

Fill in contribution from Jacobian.

Reimplemented from oomph::GeneralisedElement.

162  {
163  // Call the residuals
164  fill_in_contribution_to_residuals_axisymmetric_nst_traction(residuals);
165  }
oomph::AxisymmetricNavierStokesTractionElement::fill_in_contribution_to_residuals_axisymmetric_nst_traction
void fill_in_contribution_to_residuals_axisymmetric_nst_traction(Vector< double > &residuals)
Helper function that actually calculates the residuals.
Definition: axisym_fluid_traction_elements.h:447

References oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::fill_in_contribution_to_residuals_axisymmetric_nst_traction().

◆ fill_in_contribution_to_residuals()

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

Return the residuals.

Reimplemented from oomph::GeneralisedElement.

154  {
155  fill_in_contribution_to_residuals_axisymmetric_nst_traction(residuals);
156  }

References oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::fill_in_contribution_to_residuals_axisymmetric_nst_traction().

◆ fill_in_contribution_to_residuals_axisymmetric_nst_traction()

template<class ELEMENT >
void oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::fill_in_contribution_to_residuals_axisymmetric_nst_traction ( Vector< double > &  residuals)
protected

Helper function that actually calculates the residuals.

Return the residuals for the AxisymmetricNavierStokesTractionElement equations

449  {
450  // Find out how many nodes there are
451  unsigned n_node = nnode();
452 
453  // Get continuous time from timestepper of first node
454  double time = node_pt(0)->time_stepper_pt()->time_pt()->time();
455 
456 #ifdef PARANOID
457  // Find out how many positional dofs there are
458  unsigned n_position_type = this->nnodal_position_type();
459  if (n_position_type != 1)
460  {
461  throw OomphLibError("AxisymmetricNavierStokes is not yet implemented for "
462  "more than one position type",
463  OOMPH_CURRENT_FUNCTION,
464  OOMPH_EXCEPTION_LOCATION);
465  }
466 #endif
467 
468  // Find out the dimension of the node
469  unsigned n_dim = this->nodal_dimension();
470 
471  // Cache the nodal indices at which the velocity components are stored
472  unsigned u_nodal_index[n_dim + 1];
473  for (unsigned i = 0; i < n_dim + 1; i++)
474  {
475  u_nodal_index[i] = this->U_index_axisymmetric_nst_traction[i];
476  }
477 
478  // Integer to hold the local equation number
479  int local_eqn = 0;
480 
481  // Set up memory for the shape functions
482  // Note that in this case, the number of lagrangian coordinates is always
483  // equal to the dimension of the nodes
484  Shape psi(n_node);
485  DShape dpsids(n_node, n_dim - 1);
486 
487  // Set the value of n_intpt
488  unsigned n_intpt = integral_pt()->nweight();
489 
490  // Loop over the integration points
491  for (unsigned ipt = 0; ipt < n_intpt; ipt++)
492  {
493  // Get the integral weight
494  double w = integral_pt()->weight(ipt);
495 
496  // Only need to call the local derivatives
497  dshape_local_at_knot(ipt, psi, dpsids);
498 
499  // Calculate the Eulerian and Lagrangian coordinates
500  Vector<double> interpolated_x(n_dim, 0.0);
501 
502  // Also calculate the surface Vectors (derivatives wrt local coordinates)
503  DenseMatrix<double> interpolated_A(n_dim - 1, n_dim, 0.0);
504 
505  // Calculate positions and derivatives
506  for (unsigned l = 0; l < n_node; l++)
507  {
508  // Loop over directions
509  for (unsigned i = 0; i < n_dim; i++)
510  {
511  // Calculate the Eulerian coords
512  const double x_local = nodal_position(l, i);
513  interpolated_x[i] += x_local * psi(l);
514 
515  // Loop over LOCAL derivative directions, to calculate the tangent(s)
516  for (unsigned j = 0; j < n_dim - 1; j++)
517  {
518  interpolated_A(j, i) += x_local * dpsids(l, j);
519  }
520  }
521  }
522 
523  // Now find the local metric tensor from the tangent Vectors
524  DenseMatrix<double> A(n_dim - 1);
525  for (unsigned i = 0; i < n_dim - 1; i++)
526  {
527  for (unsigned j = 0; j < n_dim - 1; j++)
528  {
529  // Initialise surface metric tensor to zero
530  A(i, j) = 0.0;
531 
532  // Take the dot product
533  for (unsigned k = 0; k < n_dim; k++)
534  {
535  A(i, j) += interpolated_A(i, k) * interpolated_A(j, k);
536  }
537  }
538  }
539 
540  // Get the outer unit normal
541  Vector<double> interpolated_normal(n_dim);
542  outer_unit_normal(ipt, interpolated_normal);
543 
544  // Find the determinant of the metric tensor
545  double Adet = 0.0;
546  switch (n_dim)
547  {
548  case 2:
549  Adet = A(0, 0);
550  break;
551  case 3:
552  Adet = A(0, 0) * A(1, 1) - A(0, 1) * A(1, 0);
553  break;
554  default:
555  throw OomphLibError(
556  "Wrong dimension in AxisymmetricNavierStokesTractionElement",
557  "AxisymmetricNavierStokesTractionElement::fill_in_contribution_to_"
558  "residuals()",
559  OOMPH_EXCEPTION_LOCATION);
560  }
561 
562  // Premultiply the weights and the square-root of the determinant of
563  // the metric tensor
564  double W = w * sqrt(Adet);
565 
566  // Now calculate the load
567  Vector<double> traction(n_dim + 1);
568  get_traction(time, ipt, interpolated_x, interpolated_normal, traction);
569 
570  // Loop over the test functions, nodes of the element
571  for (unsigned l = 0; l < n_node; l++)
572  {
573  // Loop over the velocity components
574  for (unsigned i = 0; i < n_dim + 1; i++)
575  {
576  // Equation number
577  local_eqn = this->nodal_local_eqn(l, u_nodal_index[i]);
578  /*IF it's not a boundary condition*/
579  if (local_eqn >= 0)
580  {
581  // Add the loading terms to the residuals
582  residuals[local_eqn] -=
583  traction[i] * psi(l) * interpolated_x[0] * W;
584  }
585  }
586  } // End of loop over shape functions
587  } // End of loop over integration points
588  }
sqrt
AnnoyingScalar sqrt(const AnnoyingScalar &x)
Definition: AnnoyingScalar.h:134
w
RowVector3d w
Definition: Matrix_resize_int.cpp:3
A
Matrix< SCALARA, Dynamic, Dynamic, opt_A > A
Definition: bench_gemm.cpp:47
oomph::AxisymmetricNavierStokesTractionElement::traction
void traction(const double &time, const Vector< double > &s, Vector< double > &traction)
Definition: axisym_fluid_traction_elements.h:420
oomph::AxisymmetricNavierStokesTractionElement::get_traction
virtual void get_traction(const double &time, const unsigned &intpt, const Vector< double > &x, const Vector< double > &n, Vector< double > &traction) const
Definition: axisym_fluid_traction_elements.h:95
oomph::Data::time_stepper_pt
TimeStepper *& time_stepper_pt()
Return the pointer to the timestepper.
Definition: nodes.h:238
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::interpolated_x
double interpolated_x(const Vector< double > &s, const unsigned &i) const
Definition: elements.h:4528
oomph::FiniteElement::node_pt
Node *& node_pt(const unsigned &n)
Return a pointer to the local node n.
Definition: elements.h:2175
oomph::FiniteElement::nnodal_position_type
unsigned nnodal_position_type() const
Definition: elements.h:2463
oomph::FiniteElement::dshape_local_at_knot
virtual void dshape_local_at_knot(const unsigned &ipt, Shape &psi, DShape &dpsids) const
Definition: elements.cc:3239
oomph::FiniteElement::nodal_local_eqn
int nodal_local_eqn(const unsigned &n, const unsigned &i) const
Definition: elements.h:1432
oomph::FiniteElement::nnode
unsigned nnode() const
Return the number of nodes.
Definition: elements.h:2210
oomph::FiniteElement::integral_pt
Integral *const & integral_pt() const
Return the pointer to the integration scheme (const version)
Definition: elements.h:1963
oomph::FiniteElement::nodal_position
double nodal_position(const unsigned &n, const unsigned &i) const
Definition: elements.h:2317
oomph::FiniteElement::nodal_dimension
unsigned nodal_dimension() const
Return the required Eulerian dimension of the nodes in this element.
Definition: elements.h:2484
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::TimeStepper::time_pt
Time *const & time_pt() const
Access function for the pointer to time (const version)
Definition: timesteppers.h:572
oomph::Time::time
double & time()
Return the current value of the continuous time.
Definition: timesteppers.h:123
k
char char char int int * k
Definition: level2_impl.h:374
oomph::QuadTreeNames::W
@ W
Definition: quadtree.h:63
OOMPH_EXCEPTION_LOCATION
#define OOMPH_EXCEPTION_LOCATION
Definition: oomph_definitions.h:61
OOMPH_CURRENT_FUNCTION
#define OOMPH_CURRENT_FUNCTION
Definition: oomph_definitions.h:86
j
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2

References i, j, k, OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, sqrt(), w, and oomph::QuadTreeNames::W.

Referenced by oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::fill_in_contribution_to_jacobian(), and oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::fill_in_contribution_to_residuals().

◆ get_traction()

template<class ELEMENT >
virtual void oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::get_traction ( const double &  time,
const unsigned &  intpt,
const Vector< double > &  x,
const Vector< double > &  n,
Vector< double > &  traction 
) const
inlineprotectedvirtual

Get the traction vector: Pass number of integration point (dummy), Eulerian coordinate and normal vector and return the load vector (not all of the input arguments will be required for all specific load functions but the list should cover all cases). This function is virtual so it can be overloaded for FSI.

100  {
101  Traction_fct_pt(time, x, n, traction);
102  }
n
const unsigned n
Definition: CG3DPackingUnitTest.cpp:11
plotDoE.x
list x
Definition: plotDoE.py:28

References n, oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::traction(), oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::Traction_fct_pt, and plotDoE::x.

Referenced by oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::output(), and oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::scalar_value_paraview().

◆ nscalar_paraview()

template<class ELEMENT >
unsigned oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::nscalar_paraview ( ) const
inlinevirtual

Number of scalars/fields output by this element. Reimplements broken virtual function in base class.

Reimplemented from oomph::FiniteElement.

189  {
190  // Number of dimensions
191  unsigned n_dim = this->nodal_dimension();
192 
193  return 2 * (n_dim + 1);
194  }

References oomph::FiniteElement::nodal_dimension().

◆ output() [1/4]

template<class ELEMENT >
void oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::output ( FILE *  file_pt)
inlinevirtual

C_style output function.

Reimplemented from oomph::FiniteElement.

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

References oomph::FiniteElement::output().

◆ output() [2/4]

template<class ELEMENT >
void oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::output ( FILE *  file_pt,
const unsigned &  n_plot 
)
inlinevirtual

C-style output function.

Reimplemented from oomph::FiniteElement.

397  {
398  FiniteElement::output(file_pt, n_plot);
399  }

References oomph::FiniteElement::output().

◆ output() [3/4]

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

Output function.

Reimplemented from oomph::FiniteElement.

181  {
182  unsigned nplot = 5;
183  output(outfile, nplot);
184  }
oomph::AxisymmetricNavierStokesTractionElement::output
void output(std::ostream &outfile)
Output function.
Definition: axisym_fluid_traction_elements.h:180

◆ output() [4/4]

template<class ELEMENT >
void oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::output ( std::ostream &  outfile,
const unsigned &  n_plot 
)
inlinevirtual

Output function.

Reimplemented from oomph::FiniteElement.

311  {
312  // Number of dimensions
313  unsigned n_dim = this->nodal_dimension();
314 
315  // Find out how many nodes there are
316  const unsigned n_node = nnode();
317 
318  // Get continuous time from timestepper of first node
319  double time = node_pt(0)->time_stepper_pt()->time_pt()->time();
320 
321  // Set up memory for the shape functions
322  Shape psi(n_node);
323 
324  // Local and global coordinates
325  Vector<double> s(n_dim - 1);
326  Vector<double> interpolated_x(n_dim);
327 
328  // Tecplot header info
329  outfile << this->tecplot_zone_string(n_plot);
330 
331  // Loop over plot points
332  unsigned num_plot_points = this->nplot_points(n_plot);
333  for (unsigned iplot = 0; iplot < num_plot_points; iplot++)
334  {
335  // Get local coordinates of plot point
336  this->get_s_plot(iplot, n_plot, s);
337 
338  // Outer unit normal
339  Vector<double> unit_normal(n_dim);
340  outer_unit_normal(s, unit_normal);
341 
342  // Find the shape functions
343  shape(s, psi);
344 
345  // Initialise to zero
346  for (unsigned i = 0; i < n_dim; i++)
347  {
348  interpolated_x[i] = 0.0;
349  }
350 
351  // Calculate stuff
352  for (unsigned l = 0; l < n_node; l++)
353  {
354  // Loop over directions
355  for (unsigned i = 0; i < n_dim; i++)
356  {
357  interpolated_x[i] += this->nodal_position(l, i) * psi[l];
358  }
359  }
360 
361  // Get the imposed traction
362  Vector<double> traction(3);
363 
364  // Dummy integration point
365  unsigned ipt = 0;
366  get_traction(time, ipt, interpolated_x, unit_normal, traction);
367 
368  // Output the x,y,..
369  for (unsigned i = 0; i < n_dim; i++)
370  {
371  outfile << interpolated_x[i] << " ";
372  }
373 
374  // Output the traction components
375  for (unsigned i = 0; i < n_dim + 1; i++)
376  {
377  outfile << traction[i] << " ";
378  }
379 
380  // Output normal
381  for (unsigned i = 0; i < n_dim; i++)
382  {
383  outfile << unit_normal[i] << " ";
384  }
385  outfile << std::endl;
386  }
387  }
oomph::FiniteElement::tecplot_zone_string
virtual std::string tecplot_zone_string(const unsigned &nplot) const
Definition: elements.h:3161
oomph::FiniteElement::shape
virtual void shape(const Vector< double > &s, Shape &psi) const =0
oomph::FiniteElement::get_s_plot
virtual void get_s_plot(const unsigned &i, const unsigned &nplot, Vector< double > &s, const bool &shifted_to_interior=false) const
Definition: elements.h:3148
oomph::FiniteElement::nplot_points
virtual unsigned nplot_points(const unsigned &nplot) const
Definition: elements.h:3186
s
RealScalar s
Definition: level1_cplx_impl.h:130

References oomph::FiniteElement::get_s_plot(), oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::get_traction(), i, oomph::FaceElement::interpolated_x(), oomph::FiniteElement::nnode(), oomph::FiniteElement::nodal_dimension(), oomph::FiniteElement::nodal_position(), oomph::FiniteElement::node_pt(), oomph::FiniteElement::nplot_points(), oomph::FaceElement::outer_unit_normal(), s, oomph::FiniteElement::shape(), oomph::FiniteElement::tecplot_zone_string(), oomph::Time::time(), oomph::TimeStepper::time_pt(), oomph::Data::time_stepper_pt(), and oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::traction().

◆ scalar_name_paraview()

template<class ELEMENT >
std::string oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::scalar_name_paraview ( const unsigned &  i) const
inlinevirtual

Name of the i-th scalar field. Default implementation returns V1 for the first one, V2 for the second etc. Can (should!) be overloaded with more meaningful names in specific elements.

Reimplemented from oomph::FiniteElement.

283  {
284  // Number of dimensions
285  unsigned n_dim = this->nodal_dimension();
286 
287  // Traction components
288  if (i < n_dim + 1)
289  {
290  return "Traction component " + StringConversion::to_string(i);
291  }
292  // Normals
293  else if (i < 2 * n_dim + 1 && i >= n_dim + 1)
294  {
295  return "Normal " + StringConversion::to_string(i % (n_dim + 1));
296  }
297  // Never get here
298  else
299  {
300  std::stringstream error_stream;
301  error_stream
302  << "Axisymmetric Fluid Traction Navier-Stokes Elements only store "
303  << 2 * (n_dim + 1) << " fields " << std::endl;
304  throw OomphLibError(
305  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
306  }
307  }
oomph::StringConversion::to_string
std::string to_string(T object, unsigned float_precision=8)
Definition: oomph_utilities.h:189

References i, oomph::FiniteElement::nodal_dimension(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, and oomph::StringConversion::to_string().

◆ scalar_value_paraview()

template<class ELEMENT >
void oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::scalar_value_paraview ( std::ofstream &  file_out,
const unsigned &  k,
const unsigned &  nplot 
) const
inlinevirtual

Write values of the k-th scalar field at the plot points. Needs to be implemented for each new specific element type.

Reimplemented from oomph::FiniteElement.

201  {
202  // Number of dimensions
203  unsigned n_dim = this->nodal_dimension();
204 
205  // Find out how many nodes there are
206  const unsigned n_node = nnode();
207 
208  // Get continuous time from timestepper of first node
209  double time = node_pt(0)->time_stepper_pt()->time_pt()->time();
210 
211  // Set up memory for the shape functions
212  Shape psi(n_node);
213 
214  // Local and global coordinates
215  Vector<double> s(n_dim - 1);
216  Vector<double> interpolated_x(n_dim);
217 
218  // Loop over plot points
219  unsigned num_plot_points = this->nplot_points_paraview(nplot);
220  for (unsigned iplot = 0; iplot < num_plot_points; iplot++)
221  {
222  // Get local coordinates of plot point
223  this->get_s_plot(iplot, nplot, s);
224 
225  // Outer unit normal
226  Vector<double> unit_normal(n_dim);
227  outer_unit_normal(s, unit_normal);
228 
229  // Find the shape functions
230  shape(s, psi);
231 
232  // Initialise to zero
233  for (unsigned i = 0; i < n_dim; i++)
234  {
235  interpolated_x[i] = 0.0;
236  }
237 
238  // Calculate stuff
239  for (unsigned l = 0; l < n_node; l++)
240  {
241  // Loop over directions
242  for (unsigned i = 0; i < n_dim; i++)
243  {
244  interpolated_x[i] += this->nodal_position(l, i) * psi[l];
245  }
246  }
247 
248  // Get the imposed traction
249  Vector<double> traction(3);
250 
251  // Dummy integration point
252  unsigned ipt = 0;
253  get_traction(time, ipt, interpolated_x, unit_normal, traction);
254 
255  // Traction components
256  if (k < n_dim + 1)
257  {
258  file_out << traction[k] << std::endl;
259  }
260  // Advection Diffusion
261  else if (k < 2 * n_dim + 1 && k >= n_dim + 1)
262  {
263  file_out << unit_normal[k] << std::endl;
264  }
265  // Never get here
266  else
267  {
268  std::stringstream error_stream;
269  error_stream
270  << "Axisymmetric Fluid Traction Navier-Stokes Elements only store "
271  << 2 * (n_dim + 1) << " fields " << std::endl;
272  throw OomphLibError(error_stream.str(),
273  OOMPH_CURRENT_FUNCTION,
274  OOMPH_EXCEPTION_LOCATION);
275  }
276  }
277  }
oomph::FiniteElement::nplot_points_paraview
virtual unsigned nplot_points_paraview(const unsigned &nplot) const
Definition: elements.h:2862

References oomph::FiniteElement::get_s_plot(), oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::get_traction(), i, oomph::FaceElement::interpolated_x(), k, oomph::FiniteElement::nnode(), oomph::FiniteElement::nodal_dimension(), oomph::FiniteElement::nodal_position(), oomph::FiniteElement::node_pt(), oomph::FiniteElement::nplot_points_paraview(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, oomph::FaceElement::outer_unit_normal(), s, oomph::FiniteElement::shape(), oomph::Time::time(), oomph::TimeStepper::time_pt(), oomph::Data::time_stepper_pt(), and oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::traction().

◆ traction()

template<class ELEMENT >
void oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::traction ( const double &  time,
const Vector< double > &  s,
Vector< double > &  traction 
)

Compute traction vector at specified local coordinate Should only be used for post-processing; ignores dependence on integration point!

//////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////// Compute traction vector at specified local coordinate Should only be used for post-processing; ignores dependence on integration point!

422  {
423  unsigned n_dim = this->nodal_dimension();
424 
425  // Position vector
426  Vector<double> x(n_dim);
427  interpolated_x(s, x);
428 
429  // Outer unit normal (only in r and z direction!)
430  Vector<double> unit_normal(n_dim);
431  outer_unit_normal(s, unit_normal);
432 
433  // Dummy
434  unsigned ipt = 0;
435 
436  // Traction vector
437  get_traction(time, ipt, x, unit_normal, traction);
438  }

References s, and plotDoE::x.

Referenced by oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::get_traction(), oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::output(), and oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::scalar_value_paraview().

◆ zeta_nodal()

template<class ELEMENT >
double oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::zeta_nodal ( const unsigned &  n,
const unsigned &  k,
const unsigned &  i 
) const
inlinevirtual

Specify the value of nodal zeta from the face geometry The "global" intrinsic coordinate of the element when viewed as part of a geometric object should be given by the FaceElement representation, by default (needed to break indeterminacy if bulk element is SolidElement)

Reimplemented from oomph::FiniteElement.

175  {
176  return FaceElement::zeta_nodal(n, k, i);
177  }
oomph::FaceElement::zeta_nodal
double zeta_nodal(const unsigned &n, const unsigned &k, const unsigned &i) const
Definition: elements.h:4497

References i, k, n, and oomph::FaceElement::zeta_nodal().

Member Data Documentation

◆ traction_fct_pt

template<class ELEMENT >
void(*&)(const double& time, const Vector<double>& x, const Vector<double>& n, Vector<double>& traction) oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::traction_fct_pt()
inline

Reference to the traction function pointer.

147  {
148  return Traction_fct_pt;
149  }

Referenced by PressureWaveFSIProblem< FLUID_ELEMENT, SOLID_ELEMENT >::create_fluid_traction_elements().

◆ Traction_fct_pt

template<class ELEMENT >
void(* oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::Traction_fct_pt) (const double &time, const Vector< double > &x, const Vector< double > &n, Vector< double > &result)
protected

Pointer to an imposed traction function. Arguments: Eulerian coordinate; outer unit normal; applied traction. (Not all of the input arguments will be required for all specific load functions but the list should cover all cases)

Referenced by oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::AxisymmetricNavierStokesTractionElement(), and oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::get_traction().

◆ U_index_axisymmetric_nst_traction

template<class ELEMENT >
Vector<unsigned> oomph::AxisymmetricNavierStokesTractionElement< ELEMENT >::U_index_axisymmetric_nst_traction
protected

Index at which the i-th velocity component is stored.

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