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

#include <linearised_axisymmetric_fluid_interface_elements.h>

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

Public Member Functions

 PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement (FiniteElement *const &bulk_element_pt, const int &face_index)
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 Return the jacobian. More...
 
void fill_in_generic_residual_contribution_interface (Vector< double > &residuals, DenseMatrix< double > &jacobian, unsigned flag)
 
void output (std::ostream &outfile)
 Overload the output function. More...
 
void output (std::ostream &outfile, const unsigned &n_plot)
 Output the element. More...
 
void output (FILE *file_pt)
 Overload the C-style output function. More...
 
void output (FILE *file_pt, const unsigned &n_plot)
 C-style Output function: x,y,[z],u,v,[w],p in tecplot format. More...
 
void output_interface_position (std::ostream &outfile, const unsigned &nplot)
 Output just the interface position (base + perturbation) More...
 
void output_perturbation_to_interface (std::ostream &outfile, const unsigned &nplot)
 Output the perturbation to the interface position. More...
 
double interpolated_H (const Vector< double > &s, const unsigned &i) const
 
- Public Member Functions inherited from oomph::Hijacked< PerturbedSpineElement< FaceGeometry< ELEMENT > > >
 Hijacked ()
 Constructor, call the constructors of the base elements. More...
 
 Hijacked (FiniteElement *const &element_pt, const int &face_index)
 Constructor used for hijacking face elements. More...
 
 Hijacked (FiniteElement *const &element_pt, const int &face_index, const unsigned &id=0)
 
Data * hijack_internal_value (const unsigned &n, const unsigned &i, const bool &return_data=true)
 
Data * hijack_external_value (const unsigned &n, const unsigned &i, const bool &return_data=true)
 
Data * hijack_nodal_value (const unsigned &n, const unsigned &i, const bool &return_data=true)
 
Data * hijack_nodal_position_value (const unsigned &n, const unsigned &i, const bool &return_data=true)
 
Data * hijack_nodal_spine_value (const unsigned &n, const unsigned &i, const bool &return_data=true)
 
void assign_local_eqn_numbers (const bool &store_local_dof_pt)
 
void get_residuals (Vector< double > &residuals)
 
void get_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
- Public Member Functions inherited from oomph::PerturbedSpineElement< ELEMENT >
 PerturbedSpineElement ()
 Constructor, call the constructor of the base element. More...
 
 PerturbedSpineElement (FiniteElement *const &element_pt, const int &face_index)
 Constructor used for spine face elements. More...
 
 ~PerturbedSpineElement ()
 Destructor, clean up the storage allocated to the local equation numbers. More...
 
int spine_local_eqn (const unsigned &n, const unsigned &i)
 
virtual int get_local_eqn_number_corresponding_to_geometric_dofs (const unsigned &n, const unsigned &i)
 
- Public Member Functions inherited from oomph::ElementWithSpecificMovingNodes< ELEMENT, PerturbedSpineNode >
void describe_local_dofs (std::ostream &out, const std::string &current_string) const
 
void describe_local_dofs (std::ostream &out, std::string &curr_str)
 Unique final overrider for describe_dofs. More...
 
 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...
 
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 complete_setup_of_dependencies ()
 
void assign_all_generic_local_eqn_numbers (const bool &store_local_dof_pt)
 
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 node_update ()
 
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_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 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::PerturbedSpineFiniteElement
 PerturbedSpineFiniteElement ()
 Empty constructor. More...
 
virtual ~PerturbedSpineFiniteElement ()
 Emtpty virtual destructor. More...
 
- Public Member Functions inherited from oomph::HijackedElementBase
 HijackedElementBase ()
 
virtual ~HijackedElementBase ()
 Destructor, destroy the storage for the equation numbers. More...
 
void unhijack_all_data ()
 
const double & residual_multiplier () const
 Return the value of the residual multiplier. More...
 
double *& residual_multiplier_pt ()
 Return the pointer to the residual multiplier. More...
 
- Public Member Functions inherited from oomph::LinearisedAxisymmetricFluidInterfaceElement
 LinearisedAxisymmetricFluidInterfaceElement ()
 Constructor, set the default values of the booleans and pointers (null) More...
 
virtual double sigma (const Vector< double > &s_local)
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Calculate the residuals by calling the generic residual contribution. More...
 
const double & ca () const
 Return the value of the Capillary number. More...
 
double *& ca_pt ()
 Return a pointer to the Capillary number. More...
 
const double & st () const
 Return the value of the Strouhal number. More...
 
double *& st_pt ()
 Return a pointer to the Strouhal number. More...
 
const int & azimuthal_mode_number () const
 Azimuthal mode number k in e^ik(theta) decomposition. More...
 
int *& azimuthal_mode_number_pt ()
 Pointer to azimuthal mode number k in e^ik(theta) decomposition. More...
 
double u (const unsigned &n, const unsigned &i)
 
double interpolated_u (const Vector< double > &s, const unsigned &i)
 Calculate the i-th velocity component at the local coordinate s. 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...
 

Private Member Functions

int kinematic_local_eqn (const unsigned &n, const unsigned &i)
 
void hijack_kinematic_conditions (const Vector< unsigned > &bulk_node_number)
 
double dH_dt (const unsigned &n, const unsigned &i) const
 

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
 
- 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)
 
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 Member Functions inherited from oomph::HijackedElementBase
void hijack_global_eqn (long *const &global_eqn_pt)
 
void unhijack_global_eqn (long *const &global_eqn_pt)
 
- Protected Member Functions inherited from oomph::LinearisedAxisymmetricFluidInterfaceElement
virtual void add_additional_residual_contributions (Vector< double > &residuals, DenseMatrix< double > &jacobian, const unsigned &flag, const Shape &psif, const DShape &dpsifds, const Vector< double > &interpolated_n, const double &r, const double &W, const double &J)
 
double dXhat_dt (const unsigned &n, const unsigned &i) const
 
- Protected Member Functions inherited from oomph::FaceElement
void add_additional_values (const Vector< unsigned > &nadditional_values, const unsigned &id)
 
- 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::HijackedElementBase
std::set< long * > * Hijacked_global_eqn_number_pt
 
Vector< int > * Hijacked_local_eqn_number_pt
 
double * Residual_multiplier_pt
 
- Protected Attributes inherited from oomph::LinearisedAxisymmetricFluidInterfaceElement
Vector< unsigned > U_index_interface
 Index at which the i-th velocity component is stored. More...
 
Vector< unsigned > Xhat_index_interface
 
- 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
 
- Static Protected Attributes inherited from oomph::FiniteElement
static const unsigned Default_Initial_Nvalue = 0
 Default value for the number of values at a node. More...
 
static const double Node_location_tolerance = 1.0e-14
 
static const unsigned N2deriv [] = {0, 1, 3, 6}
 
- Static Protected Attributes inherited from oomph::GeneralisedElement
static DenseMatrix< double > Dummy_matrix
 
static std::deque< double * > Dof_pt_deque
 
- Static Protected Attributes inherited from oomph::HijackedElementBase
static double Default_residual_multiplier = 0.0
 

Detailed Description

template<class ELEMENT>
class oomph::PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement< ELEMENT >

Linearised axisymmetric interface elements that are used with a spine mesh, i.e. the mesh deformation is handled by Kistler & Scriven's "method of spines". These elements are FaceElements of bulk Fluid elements and the particular type of fluid element is passed as a template parameter to the element. It shouldn't matter whether the passed element is the underlying (fixed) element or the templated SpineElement<Element>. In the case of steady problems, an additional volume constaint must be imposed to select a particular solution from an otherwise inifinte number. This constraint is associated with an external pressure degree of freedom, which must be passed to the element as external data. If the element is a free surface, Free_surface = true, then the external pressure is the pressure in the inviscid external fluid; otherwise, the pressure degree of freedom must be passed from an element in the adjoining fluid.

Constructor & Destructor Documentation

◆ PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement()

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

Constructor, the arguments are a pointer to the "bulk" element the local coordinate that is fixed on the face, and whether it is the upper or lower limit of that coordinate.

323  :
324  Hijacked<PerturbedSpineElement<FaceGeometry<ELEMENT> > >(),
325  LinearisedAxisymmetricFluidInterfaceElement()
326  {
327  // Attach the geometrical information to the element
328  bulk_element_pt->build_face_element(face_index,this);
329 
330  // Find the index at which the velocity unknowns are stored
331  // from the bulk element
332  ELEMENT* cast_element_pt = dynamic_cast<ELEMENT*>(bulk_element_pt);
333 
334  // We must have six velocity components
335  this->U_index_interface.resize(6);
336  for(unsigned i=0;i<6;i++)
337  {
338  this->U_index_interface[i]=cast_element_pt->u_index_lin_axi_nst(i);
339  }
340 
341  // We must have four components of the perturbations to the
342  // nodal positions
343  this->Xhat_index_interface.resize(4);
344  for(unsigned i=0;i<4;i++)
345  {
346  this->Xhat_index_interface[i] =
347  cast_element_pt->xhat_index_lin_axi_nst(i);
348  }
349  }
i
int i
Definition: BiCGSTAB_step_by_step.cpp:9
oomph::FaceElement::face_index
int & face_index()
Definition: elements.h:4626
oomph::FaceElement::bulk_element_pt
FiniteElement *& bulk_element_pt()
Pointer to higher-dimensional "bulk" element.
Definition: elements.h:4735
oomph::FiniteElement::build_face_element
virtual void build_face_element(const int &face_index, FaceElement *face_element_pt)
Definition: elements.cc:5132
oomph::LinearisedAxisymmetricFluidInterfaceElement::Xhat_index_interface
Vector< unsigned > Xhat_index_interface
Definition: linearised_axisymmetric_fluid_interface_elements.h:75
oomph::LinearisedAxisymmetricFluidInterfaceElement::U_index_interface
Vector< unsigned > U_index_interface
Index at which the i-th velocity component is stored.
Definition: linearised_axisymmetric_fluid_interface_elements.h:71
oomph::LinearisedAxisymmetricFluidInterfaceElement::LinearisedAxisymmetricFluidInterfaceElement
LinearisedAxisymmetricFluidInterfaceElement()
Constructor, set the default values of the booleans and pointers (null)
Definition: linearised_axisymmetric_fluid_interface_elements.h:150

References oomph::FiniteElement::build_face_element(), oomph::FaceElement::bulk_element_pt(), oomph::FaceElement::face_index(), i, oomph::LinearisedAxisymmetricFluidInterfaceElement::U_index_interface, and oomph::LinearisedAxisymmetricFluidInterfaceElement::Xhat_index_interface.

Member Function Documentation

◆ dH_dt()

template<class ELEMENT >
double oomph::PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement< ELEMENT >::dH_dt ( const unsigned &  n,
const unsigned &  i 
) const
inlineprivate

i-th component of dH/dt at local node n. Uses suitably interpolated value for hanging nodes.

289  {
290  // Get the data's timestepper
291  TimeStepper* time_stepper_pt = this->node_pt(n)->time_stepper_pt();
292 
293  // Upcast from general node to PerturbedSpineNode
294  PerturbedSpineNode* perturbed_spine_node_pt =
295  dynamic_cast<PerturbedSpineNode*>(this->node_pt(n));
296 
297  // Initialise dH/dt
298  double dHdt = 0.0;
299 
300  // Loop over the timesteps, if there is a non Steady timestepper
301  if(!time_stepper_pt->is_steady())
302  {
303  // Number of timsteps (past & present)
304  const unsigned n_time = time_stepper_pt->ntstorage();
305 
306  // Add the contributions to the time derivative
307  for(unsigned t=0;t<n_time;t++)
308  {
309  dHdt += time_stepper_pt->weight(1,t)*
310  perturbed_spine_node_pt->perturbed_spine_pt()->height(t,i);
311  }
312  }
313 
314  return dHdt;
315  }
oomph::Data::time_stepper_pt
TimeStepper *& time_stepper_pt()
Return the pointer to the timestepper.
Definition: nodes.h:238
oomph::FiniteElement::node_pt
Node *& node_pt(const unsigned &n)
Return a pointer to the local node n.
Definition: elements.h:2175
oomph::GeomObject::time_stepper_pt
TimeStepper *& time_stepper_pt()
Definition: geom_objects.h:192
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::is_steady
bool is_steady() const
Definition: timesteppers.h:389
plotPSD.t
t
Definition: plotPSD.py:36

References oomph::PerturbedSpine::height(), i, oomph::TimeStepper::is_steady(), oomph::FiniteElement::node_pt(), oomph::TimeStepper::ntstorage(), oomph::PerturbedSpineNode::perturbed_spine_pt(), plotPSD::t, oomph::GeomObject::time_stepper_pt(), oomph::Data::time_stepper_pt(), and oomph::TimeStepper::weight().

◆ fill_in_contribution_to_jacobian()

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

Return the jacobian.

Reimplemented from oomph::GeneralisedElement.

354  {
355  // Call the generic residuals routine with the flag set to 1
356  fill_in_generic_residual_contribution_interface(residuals,jacobian,1);
357  // Call the generic routine to handle the spine variables
358  // PerturbedSpineElement<FaceGeometry<ELEMENT> >::
359  this->fill_in_jacobian_from_geometric_data(jacobian);
360  }
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::PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement::fill_in_generic_residual_contribution_interface
void fill_in_generic_residual_contribution_interface(Vector< double > &residuals, DenseMatrix< double > &jacobian, unsigned flag)
Definition: linearised_axisymmetric_fluid_interface_elements.cc:216

References oomph::PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement< ELEMENT >::fill_in_generic_residual_contribution_interface(), and oomph::ElementWithMovingNodes::fill_in_jacobian_from_geometric_data().

◆ fill_in_generic_residual_contribution_interface()

template<class ELEMENT >
void oomph::PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement< ELEMENT >::fill_in_generic_residual_contribution_interface ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian,
unsigned  flag 
)
virtual

This function is sort-of-hacky: its purpose is to overload the fill_in_generic_residual_contribution_interface(..) function in the LinearisedAxisymmetricFluidInterfaceElement base class, since I have made the maths too general and assumed not only spines but spines which always point vertically. When this is done properly this function should be removed so that the function in the base class is used, and that function should be implemented with the proper general maths which does not assume spines or anything.

Calculate the residuals for the linearised axisymmetric interface element. WARNING: THIS FUNCTION ASSUMES THAT THE PERTURBATION TO THE INTERFACE IS IN A VERTICAL DIRECTION ONLY

Reimplemented from oomph::LinearisedAxisymmetricFluidInterfaceElement.

220  {
221  // Determine number of nodes in the element
222  const unsigned n_node = this->nnode();
223 
224  // Set up memory for the shape functions and their derivative w.r.t. the
225  // local coordinate s
226  Shape psif(n_node);
227  DShape dpsifds(n_node,1);
228 
229  // Set up memory for the test functions and their derivative w.r.t. the
230  // local coordinate s
231  Shape testf(n_node);
232  DShape dtestfds(n_node,1);
233 
234  // Storage for the local coordinate
235  Vector<double> s(1);
236 
237  // Storage for the local coordinate in the parent element which
238  // corresponds to the local coordinate in this element
239  Vector<double> s_parent(2);
240 
241  // Get a pointer to the parent element
242  LinearisedAxisymmetricNavierStokesEquations* bulk_el_pt =
243  dynamic_cast<LinearisedAxisymmetricNavierStokesEquations*>
244  (bulk_element_pt());
245 
246  // Determine the number of nodes in the parent element
247  const unsigned n_node_parent = bulk_el_pt->nnode();
248 
249  // Set up memory for the parent test functions and their derivatives
250  Shape testf_parent(n_node_parent);
251  DShape dtestfdx_parent(n_node_parent,2);
252 
253  // Determine number of integration points
254  const unsigned n_intpt = this->integral_pt()->nweight();
255 
256  // Get physical variables from the element
257  const double Ca = ca();
258  const double St = st();
259  const int k = azimuthal_mode_number();
260 
261  // Integers used to store the local equation and unknown numbers
262  int local_eqn = 0, local_unknown = 0;
263 
264  // Loop over the integration points
265  for(unsigned ipt=0;ipt<n_intpt;ipt++)
266  {
267  // Get the local coordinate at the integration point
268  s[0] = integral_pt()->knot(ipt,0);
269 
270  // Get the corresponding local coordinate in the parent (bulk) element
271  this->get_local_coordinate_in_bulk(s,s_parent);
272 
273  // Get the integral weight
274  const double w = this->integral_pt()->weight(ipt);
275 
276  // Call the derivatives of the shape function
277  this->dshape_local_at_knot(ipt,psif,dpsifds);
278 
279  // Set test functions equal to shape functions
280  testf = psif;
281  dtestfds = dpsifds;
282 
283  // Find the test functions and derivatives of the parent
284  (void)bulk_el_pt->dshape_eulerian(s_parent,testf_parent,dtestfdx_parent);
285 
286  // Define and zero the tangent Vectors
287  double interpolated_t1[2] = {0.0,0.0};
288  Vector<double> interpolated_x(2,0.0);
289  double interpolated_dx_dt[2] = {0.0,0.0};
290 
291  // Provide storage for velocity unknowns
292  // (only need four of them in these equations)
293  double interpolated_UC = 0.0;
294  double interpolated_US = 0.0;
295  double interpolated_WC = 0.0;
296  double interpolated_WS = 0.0;
297 
298  // Provide storage for perturbed spine "heights"
299  double interpolated_HC = 0.0;
300  double interpolated_HS = 0.0;
301 
302  // Provide storage for derivatives of perturbed spine "heights" w.r.t.
303  // the local coordinate
304  double interpolated_dHCds = 0.0;
305  double interpolated_dHSds = 0.0;
306 
307  // Provide storage for derivatives of perturbed spine "heights" w.r.t.
308  // time
309  double interpolated_dHCdt = 0.0;
310  double interpolated_dHSdt = 0.0;
311 
312  // Loop over the shape functions
313  for(unsigned l=0;l<n_node;l++)
314  {
315  // Upcast from general node to PerturbedSpineNode
316  PerturbedSpineNode* perturbed_spine_node_pt =
317  dynamic_cast<PerturbedSpineNode*>(this->node_pt(l));
318 
319  // Cache the shape function and its derivative
320  const double psif_ = psif(l);
321  const double dpsifds_ = dpsifds(l,0);
322 
323  // Calculate interpolated velocity components
324  interpolated_UC += u(l,0)*psif_;
325  interpolated_US += u(l,1)*psif_;
326  interpolated_WC += u(l,2)*psif_;
327  interpolated_WS += u(l,3)*psif_;
328 
329  // Calculate interpolated perturbed spine "heights"
330  interpolated_HC +=
331  perturbed_spine_node_pt->perturbed_spine_pt()->height(0)*psif_;
332  interpolated_HS +=
333  perturbed_spine_node_pt->perturbed_spine_pt()->height(1)*psif_;
334 
335  // Calculate derivatives of interpolated perturbed spine "heights"
336  interpolated_dHCds +=
337  perturbed_spine_node_pt->perturbed_spine_pt()->height(0)*dpsifds_;
338  interpolated_dHSds +=
339  perturbed_spine_node_pt->perturbed_spine_pt()->height(1)*dpsifds_;
340  interpolated_dHCdt += this->dH_dt(l,0)*psif_;
341  interpolated_dHSdt += this->dH_dt(l,1)*psif_;
342 
343  // Loop over directional components
344  for(unsigned i=0;i<2;i++)
345  {
346  interpolated_x[i] += this->nodal_position(l,i)*psif_;
347  interpolated_dx_dt[i] += this->dnodal_position_dt(l,i)*psif_;
348 
349  // Calculate the tangent vector
350  interpolated_t1[i] += this->nodal_position(l,i)*dpsifds_;
351  }
352  }
353 
354 
355  // Get velocities from base flow problem
356  // -------------------------------------
357 
358  // Allocate storage for the velocity components of the base state
359  // solution (initialise to zero)
360  Vector<double> base_flow_u(3,0.0);
361 
362  // Get the user-defined base state solution velocity components
363  bulk_el_pt->get_base_flow_u(
364  bulk_el_pt->node_pt(0)->time_stepper_pt()->time(),
365  ipt,
366  interpolated_x,
367  base_flow_u);
368 
369  // Cache base flow velocities
370  const double interpolated_ur = base_flow_u[0];
371  const double interpolated_utheta = base_flow_u[2];
372 
373  // The first positional coordinate is the radial coordinate
374  const double r = interpolated_x[0];
375 
376  // Calculate the length of the tangent Vector
377  const double tlength = interpolated_t1[0]*interpolated_t1[0] +
378  interpolated_t1[1]*interpolated_t1[1];
379 
380  // Set the Jacobian of the line element
381  const double J = sqrt(tlength);
382 
383  // Normalise the tangent Vector
384  interpolated_t1[0] /= J; interpolated_t1[1] /= J;
385 
386  // Now calculate the normal Vector
387  Vector<double> interpolated_n(2);
388  outer_unit_normal(ipt,interpolated_n);
389 
390  // Also get the (possibly variable) surface tension
391  const double Sigma = this->sigma(s);
392 
393  // Loop over the test functions
394  for(unsigned l=0;l<n_node;l++)
395  {
396  // Cache test function and its derivative
397  const double testf_ = testf(l);
398  const double dtestfds_ = dtestfds(l,0);
399 
400  // =================================================
401  // START OF DYNAMIC BOUNDARY CONDITION CONTRIBUTIONS
402  // =================================================
403 
404  // -------------------------------------------------------------
405  // Contribution to first (radial) momentum equation: cosine part
406  // -------------------------------------------------------------
407 
408  // Get local equation number of first velocity value at this node
409  local_eqn = this->nodal_local_eqn(l,this->U_index_interface[0]);
410 
411  // If it's not a boundary condition
412  if(local_eqn >= 0)
413  {
414  residuals[local_eqn] -=
415  (Sigma/Ca)*interpolated_t1[1]*interpolated_dHCds*testf_*w;
416  residuals[local_eqn] +=
417  (Sigma/Ca)*r*interpolated_t1[1]*interpolated_t1[0]
418  *interpolated_dHCds*dtestfds_*w/J;
419  residuals[local_eqn] +=
420  (Sigma/Ca)*k*k*J*interpolated_t1[1]*interpolated_t1[0]
421  *interpolated_HC*testf_*w/r;
422 
423  // Calculate the Jacobian
424  // ----------------------
425  if(flag)
426  {
427  // Loop over the shape functions again
428  for(unsigned l2=0;l2<n_node;l2++)
429  {
430  // Cache the shape function and its derivative
431  const double psif_ = psif(l2);
432  const double dpsifds_ = dpsifds(l2,0);
433 
434  // Perturbed spine "height" (cosine part) H_k^C
435  local_unknown = kinematic_local_eqn(l2,0);
436  if(local_unknown >= 0)
437  {
438  jacobian(local_eqn,local_unknown) -=
439  (Sigma/Ca)*interpolated_t1[1]*dpsifds_*testf_*w;
440  jacobian(local_eqn,local_unknown) +=
441  (Sigma/Ca)*r*interpolated_t1[1]*interpolated_t1[0]
442  *dpsifds_*dtestfds_*w/J;
443  jacobian(local_eqn,local_unknown) +=
444  (Sigma/Ca)*k*k*J*interpolated_t1[1]*interpolated_t1[0]
445  *psif_*testf_*w/r;
446  }
447 
448  // Perturbed spine "height" (sine part) H_K^S
449  // has no contribution
450 
451  } // End of loop over shape functions
452  } // End of Jacobian calculation
453  } // End of if not boundary condition statement
454 
455  // ------------------------------------------------------------
456  // Contribution to second (radial) momentum equation: sine part
457  // ------------------------------------------------------------
458 
459  // Get local equation number of second velocity value at this node
460  local_eqn = this->nodal_local_eqn(l,this->U_index_interface[1]);
461 
462  // If it's not a boundary condition
463  if(local_eqn >= 0)
464  {
465  residuals[local_eqn] -=
466  (Sigma/Ca)*interpolated_t1[1]*interpolated_dHSds*testf_*w;
467  residuals[local_eqn] +=
468  (Sigma/Ca)*r*interpolated_t1[1]*interpolated_t1[0]
469  *interpolated_dHSds*dtestfds_*w/J;
470  residuals[local_eqn] += (Sigma/Ca)*k*k*J*interpolated_t1[1]
471  *interpolated_t1[0]*interpolated_HS*testf_*w/r;
472 
473  // Calculate the Jacobian
474  // ----------------------
475  if(flag)
476  {
477  // Loop over the shape functions again
478  for(unsigned l2=0;l2<n_node;l2++)
479  {
480  // Cache the shape function and its derivative
481  const double psif_ = psif(l2);
482  const double dpsifds_ = dpsifds(l2,0);
483 
484  // Perturbed spine "height" (cosine part) H_k^C
485  // has no contribution
486 
487  // Perturbed spine "height" (sine part) H_K^S
488  local_unknown = kinematic_local_eqn(l2,1);
489  if(local_unknown >= 0)
490  {
491  jacobian(local_eqn,local_unknown) -=
492  (Sigma/Ca)*interpolated_t1[1]*dpsifds_*testf_*w;
493  jacobian(local_eqn,local_unknown) +=
494  (Sigma/Ca)*r*interpolated_t1[1]*interpolated_t1[0]
495  *dpsifds_*dtestfds_*w/J;
496  jacobian(local_eqn,local_unknown) +=
497  (Sigma/Ca)*k*k*J*interpolated_t1[1]*interpolated_t1[0]
498  *psif_*testf_*w/r;
499  }
500  } // End of loop over shape functions
501  } // End of Jacobian calculation
502  } // End of if not boundary condition statement
503 
504  // ------------------------------------------------------------
505  // Contribution to third (axial) momentum equation: cosine part
506  // ------------------------------------------------------------
507 
508  // Get local equation number of third velocity value at this node
509  local_eqn = this->nodal_local_eqn(l,this->U_index_interface[2]);
510 
511  // If it's not a boundary condition
512  if(local_eqn >= 0)
513  {
514  residuals[local_eqn] -=
515  (Sigma/Ca)*r*interpolated_t1[0]*interpolated_t1[0]
516  *interpolated_dHCds*dtestfds_*w/J;
517  residuals[local_eqn] -= (Sigma/Ca)*k*k*J*interpolated_t1[0]
518  *interpolated_t1[0]*interpolated_HC*testf_*w/r;
519 
520  // Calculate the Jacobian
521  // ----------------------
522  if(flag)
523  {
524  // Loop over the shape functions again
525  for(unsigned l2=0;l2<n_node;l2++)
526  {
527  // Cache the shape function and its derivative
528  const double psif_ = psif(l2);
529  const double dpsifds_ = dpsifds(l2,0);
530 
531  // Perturbed spine "height" (cosine part) H_k^C
532  local_unknown = kinematic_local_eqn(l2,0);
533  if(local_unknown >= 0)
534  {
535  jacobian(local_eqn,local_unknown) -=
536  (Sigma/Ca)*r*interpolated_t1[0]*interpolated_t1[0]
537  *dpsifds_*dtestfds_*w/J;
538  jacobian(local_eqn,local_unknown) -=
539  (Sigma/Ca)*k*k*J*interpolated_t1[0]*interpolated_t1[0]
540  *psif_*testf_*w/r;
541  }
542 
543  // Perturbed spine "height" (sine part) H_K^S
544  // has no contribution
545 
546  } // End of loop over shape functions
547  } // End of Jacobian calculation
548  } // End of if not boundary condition statement
549 
550  // -----------------------------------------------------------
551  // Contribution to fourth (axial) momentum equation: sine part
552  // -----------------------------------------------------------
553 
554  // Get local equation number of fourth velocity value at this node
555  local_eqn = this->nodal_local_eqn(l,this->U_index_interface[3]);
556 
557  // If it's not a boundary condition
558  if(local_eqn >= 0)
559  {
560  residuals[local_eqn] -=
561  (Sigma/Ca)*r*interpolated_t1[0]*interpolated_t1[0]
562  *interpolated_dHSds*dtestfds_*w/J;
563  residuals[local_eqn] -=
564  (Sigma/Ca)*k*k*J*interpolated_t1[0]*interpolated_t1[0]
565  *interpolated_HS*testf_*w/r;
566 
567  // Calculate the Jacobian
568  // ----------------------
569  if(flag)
570  {
571  // Loop over the shape functions again
572  for(unsigned l2=0;l2<n_node;l2++)
573  {
574  // Cache the shape function and its derivative
575  const double psif_ = psif(l2);
576  const double dpsifds_ = dpsifds(l2,0);
577 
578  // Perturbed spine "height" (cosine part) H_k^C
579  // has no contribution
580 
581  // Perturbed spine "height" (sine part) H_K^S
582  local_unknown = kinematic_local_eqn(l2,1);
583  if(local_unknown >= 0)
584  {
585  jacobian(local_eqn,local_unknown) -=
586  (Sigma/Ca)*r*interpolated_t1[0]*interpolated_t1[0]*dpsifds_
587  *dtestfds_*w/J;
588  jacobian(local_eqn,local_unknown) -=
589  (Sigma/Ca)*k*k*J*interpolated_t1[0]*interpolated_t1[0]*psif_
590  *testf_*w/r;
591  }
592  } // End of loop over shape functions
593  } // End of Jacobian calculation
594  } // End of if not boundary condition statement
595 
596  // ----------------------------------------------------------------
597  // Contribution to fifth (azimuthal) momentum equation: cosine part
598  // ----------------------------------------------------------------
599 
600  // Get local equation number of fifth velocity value at this node
601  local_eqn = this->nodal_local_eqn(l,this->U_index_interface[4]);
602 
603  // If it's not a boundary condition
604  if(local_eqn >= 0)
605  {
606  residuals[local_eqn] -=
607  (Sigma/Ca)*2.0*k*J*interpolated_t1[1]*interpolated_t1[0]
608  *interpolated_HS*testf_*w/(r*r);
609  residuals[local_eqn] +=
610  (Sigma/Ca)*k*interpolated_t1[1]*interpolated_HS*dtestfds_*w/r;
611  residuals[local_eqn] +=
612  (Sigma/Ca)*k*interpolated_t1[1]*interpolated_dHSds*testf_*w/r;
613 
614  // Calculate the Jacobian
615  // ----------------------
616  if(flag)
617  {
618  // Loop over the shape functions again
619  for(unsigned l2=0;l2<n_node;l2++)
620  {
621  // Cache the shape function and its derivative
622  const double psif_ = psif(l2);
623  const double dpsifds_ = dpsifds(l2,0);
624 
625  // Perturbed spine "height" (cosine part) H_k^C
626  // has no contribution
627 
628  // Perturbed spine "height" (sine part) H_K^S
629  local_unknown = kinematic_local_eqn(l2,1);
630  if(local_unknown >= 0)
631  {
632  jacobian(local_eqn,local_unknown) -=
633  (Sigma/Ca)*2.0*k*J*interpolated_t1[1]*interpolated_t1[0]
634  *psif_*testf_*w/(r*r);
635  jacobian(local_eqn,local_unknown) +=
636  (Sigma/Ca)*k*interpolated_t1[1]*psif_*dtestfds_*w/r;
637  jacobian(local_eqn,local_unknown) +=
638  (Sigma/Ca)*k*interpolated_t1[1]*dpsifds_*testf_*w/r;
639  }
640  } // End of loop over shape functions
641  } // End of Jacobian calculation
642  } // End of if not boundary condition statement
643 
644  // --------------------------------------------------------------
645  // Contribution to sixth (azimuthal) momentum equation: sine part
646  // --------------------------------------------------------------
647 
648  // Get local equation number of sixth velocity value at this node
649  local_eqn = this->nodal_local_eqn(l,this->U_index_interface[5]);
650 
651  // If it's not a boundary condition
652  if(local_eqn >= 0)
653  {
654  residuals[local_eqn] +=
655  (Sigma/Ca)*2.0*k*J*interpolated_t1[1]*interpolated_t1[0]
656  *interpolated_HC*testf_*w/(r*r);
657  residuals[local_eqn] -=
658  (Sigma/Ca)*k*interpolated_t1[1]*interpolated_HC*dtestfds_*w/r;
659  residuals[local_eqn] -=
660  (Sigma/Ca)*k*interpolated_t1[1]*interpolated_dHCds*testf_*w/r;
661 
662  // Calculate the Jacobian
663  // ----------------------
664  if(flag)
665  {
666  // Loop over the shape functions again
667  for(unsigned l2=0;l2<n_node;l2++)
668  {
669  // Cache the shape function and its derivative
670  const double psif_ = psif(l2);
671  const double dpsifds_ = dpsifds(l2,0);
672 
673  // Perturbed spine "height" (cosine part) H_k^C
674  local_unknown = kinematic_local_eqn(l2,0);
675  if(local_unknown >= 0)
676  {
677  jacobian(local_eqn,local_unknown) +=
678  (Sigma/Ca)*2.0*k*J*interpolated_t1[1]*interpolated_t1[0]
679  *psif_*testf_*w/(r*r);
680  jacobian(local_eqn,local_unknown) -=
681  (Sigma/Ca)*k*interpolated_t1[1]*psif_*dtestfds_*w/r;
682  jacobian(local_eqn,local_unknown) -=
683  (Sigma/Ca)*k*interpolated_t1[1]*dpsifds_*testf_*w/r;
684  }
685 
686  // Perturbed spine "height" (sine part) H_K^S
687  // has no contribution
688 
689  } // End of loop over shape functions
690  } // End of Jacobian calculation
691  } // End of if not boundary condition statement
692 
693  // ===================================================
694  // START OF KINEMATIC BOUNDARY CONDITION CONTRIBUTIONS
695  // ===================================================
696 
697  // Using the same shape functions for the spines, so can stay inside
698  // the loop over test functions
699 
700  // ------------------------------------------------
701  // First kinematic boundary condition (cosine part)
702  // ------------------------------------------------
703 
704  // Get local equation number of first perturbed spine "height"
705  // (this is the cosine part H_k^C)
706  local_eqn = kinematic_local_eqn(l,0);
707 
708  // If the spine is not a boundary condition
709  if(local_eqn >= 0)
710  {
711  residuals[local_eqn] -=
712  r*interpolated_ur*interpolated_dHCds*testf_*w;
713 
714  residuals[local_eqn] -=
715  k*interpolated_utheta*interpolated_t1[0]*interpolated_HS*testf_*w*J;
716 
717  residuals[local_eqn] -=
718  r*interpolated_t1[1]*interpolated_UC*testf_*w*J;
719 
720  residuals[local_eqn] +=
721  r*interpolated_t1[0]*interpolated_WC*testf_*w*J;
722 
723  residuals[local_eqn] -=
724  St*r*interpolated_t1[0]*interpolated_dHCdt*testf_*w*J;
725 
726  // Calculate the Jacobian
727  // ----------------------
728  if(flag)
729  {
730  // Loop over the shape functions again
731  for(unsigned l2=0;l2<n_node;l2++)
732  {
733  // Cache the shape function and its derivative
734  const double psif_ = psif(l2);
735  const double dpsifds_ = dpsifds(l2,0);
736 
737  // Radial velocity component (cosine part) U_k^C
738  local_unknown = this->nodal_local_eqn(l2,U_index_interface[0]);
739  if(local_unknown >= 0)
740  {
741  jacobian(local_eqn,local_unknown) -=
742  r*interpolated_t1[1]*psif_*testf_*w*J;
743  }
744 
745  // Axial velocity component (cosine part) W_k^C
746  local_unknown = nodal_local_eqn(l2,U_index_interface[2]);
747  if(local_unknown >= 0)
748  {
749  jacobian(local_eqn,local_unknown) +=
750  r*interpolated_t1[0]*psif_*testf_*w*J;
751  }
752 
753  // Perturbed spine "height" (cosine part) H_k^C
754  local_unknown = kinematic_local_eqn(l2,0);
755  if(local_unknown >= 0)
756  {
757  jacobian(local_eqn,local_unknown) -=
758  r*interpolated_ur*dpsifds_*testf_*w;
759 
760  jacobian(local_eqn,local_unknown) -=
761  St*r*interpolated_t1[0]
762  *node_pt(l2)->time_stepper_pt()->weight(1,0)*psif_*testf_*w*J;
763  }
764 
765  // Perturbed spine "height" (sine part) H_k^S
766  local_unknown = kinematic_local_eqn(l2,1);
767  if(local_unknown >= 0)
768  {
769  jacobian(local_eqn,local_unknown) -=
770  k*interpolated_utheta*interpolated_t1[0]*psif_*testf_*w*J;
771  }
772 
773  } // End of loop over shape functions
774  } // End of Jacobian contribution
775  } // End of if not boundary condition statement
776 
777  // -----------------------------------------------
778  // Second kinematic boundary condition (sine part)
779  // -----------------------------------------------
780 
781  // Get local equation number of second perturbed spine "height"
782  // (this is the cosine part H_k^S)
783  local_eqn = kinematic_local_eqn(l,1);
784 
785  // If the spine is not a boundary condition
786  if(local_eqn >= 0)
787  {
788  residuals[local_eqn] -=
789  r*interpolated_ur*interpolated_dHSds*testf_*w;
790 
791  residuals[local_eqn] +=
792  k*interpolated_utheta*interpolated_t1[0]*interpolated_HC*testf_*w*J;
793 
794  residuals[local_eqn] -=
795  r*interpolated_t1[1]*interpolated_US*testf_*w*J;
796 
797  residuals[local_eqn] +=
798  r*interpolated_t1[0]*interpolated_WS*testf_*w*J;
799 
800  residuals[local_eqn] -=
801  St*r*interpolated_t1[0]*interpolated_dHSdt*testf_*w*J;
802 
803  // Add in the jacobian
804  if(flag)
805  {
806  // Loop over velocity shape functions
807  for(unsigned l2=0;l2<n_node;l2++)
808  {
809  // Cache the shape function and its derivative
810  const double psif_ = psif(l2);
811  const double dpsifds_ = dpsifds(l2,0);
812 
813  // Radial velocity component (sine part) U_k^S
814  local_unknown = this->nodal_local_eqn(l2,U_index_interface[1]);
815  if(local_unknown >= 0)
816  {
817  jacobian(local_eqn,local_unknown) -=
818  r*interpolated_t1[1]*psif_*testf_*w*J;
819  }
820 
821  // Axial velocity component (sine part) W_k^S
822  local_unknown = nodal_local_eqn(l2,U_index_interface[3]);
823  if(local_unknown >= 0)
824  {
825  jacobian(local_eqn,local_unknown) +=
826  r*interpolated_t1[0]*psif_*testf_*w*J;
827  }
828 
829  // Perturbed spine "height" (cosine part) H_k^C
830  local_unknown = kinematic_local_eqn(l2,0);
831  if(local_unknown >= 0)
832  {
833  jacobian(local_eqn,local_unknown) +=
834  k*interpolated_utheta*interpolated_t1[0]*psif_*testf_*w*J;
835  }
836 
837  // Perturbed spine "height" (sine part) H_k^S
838  local_unknown = kinematic_local_eqn(l2,1);
839  if(local_unknown >= 0)
840  {
841  jacobian(local_eqn,local_unknown) -=
842  r*interpolated_ur*dpsifds_*testf_*w;
843 
844  jacobian(local_eqn,local_unknown) -=
845  St*r*interpolated_t1[0]
846  *node_pt(l2)->time_stepper_pt()->weight(1,0)*psif_*testf_*w*J;
847  }
848 
849  } // End of loop over shape functions
850  } // End of Jacobian contribution
851  } // End of if not boundary condition statement
852 
853  } // End of loop over test functions
854 
855  } // End of loop over integration points
856 
857  } // End of fill_in_generic_residual_contribution_interface
sqrt
AnnoyingScalar sqrt(const AnnoyingScalar &x)
Definition: AnnoyingScalar.h:134
J
JacobiRotation< float > J
Definition: Jacobi_makeJacobi.cpp:3
w
RowVector3d w
Definition: Matrix_resize_int.cpp:3
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::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::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::dnodal_position_dt
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.
Definition: elements.h:2333
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::LinearisedAxisymmetricFluidInterfaceElement::sigma
virtual double sigma(const Vector< double > &s_local)
Definition: linearised_axisymmetric_fluid_interface_elements.h:165
oomph::LinearisedAxisymmetricFluidInterfaceElement::ca
const double & ca() const
Return the value of the Capillary number.
Definition: linearised_axisymmetric_fluid_interface_elements.h:176
oomph::LinearisedAxisymmetricFluidInterfaceElement::st
const double & st() const
Return the value of the Strouhal number.
Definition: linearised_axisymmetric_fluid_interface_elements.h:182
oomph::LinearisedAxisymmetricFluidInterfaceElement::azimuthal_mode_number
const int & azimuthal_mode_number() const
Azimuthal mode number k in e^ik(theta) decomposition.
Definition: linearised_axisymmetric_fluid_interface_elements.h:188
oomph::LinearisedAxisymmetricFluidInterfaceElement::u
double u(const unsigned &n, const unsigned &i)
Definition: linearised_axisymmetric_fluid_interface_elements.h:197
oomph::PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement::dH_dt
double dH_dt(const unsigned &n, const unsigned &i) const
Definition: linearised_axisymmetric_fluid_interface_elements.h:288
oomph::PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement::kinematic_local_eqn
int kinematic_local_eqn(const unsigned &n, const unsigned &i)
Definition: linearised_axisymmetric_fluid_interface_elements.h:267
s
RealScalar s
Definition: level1_cplx_impl.h:130
k
char char char int int * k
Definition: level2_impl.h:374
GlobalParameters::Sigma
double Sigma
Non-dim Stefan Boltzmann constant.
Definition: stefan_boltzmann.cc:115
Global_Physical_Variables::St
double St
Strouhal number.
Definition: elastic_two_layer_interface_axisym.cc:66
Global_Physical_Variables::Ca
double Ca
Capillary number.
Definition: fibre.cc:61
UniformPSDSelfTest.r
r
Definition: UniformPSDSelfTest.py:20

References Global_Physical_Variables::Ca, oomph::FiniteElement::dshape_eulerian(), oomph::LinearisedAxisymmetricNavierStokesEquations::get_base_flow_u(), oomph::PerturbedSpine::height(), i, J, k, oomph::FiniteElement::nnode(), oomph::FiniteElement::node_pt(), oomph::PerturbedSpineNode::perturbed_spine_pt(), UniformPSDSelfTest::r, s, calibrate::sigma, GlobalParameters::Sigma, sqrt(), Global_Physical_Variables::St, oomph::TimeStepper::time(), oomph::Data::time_stepper_pt(), and w.

Referenced by oomph::PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement< ELEMENT >::fill_in_contribution_to_jacobian().

◆ hijack_kinematic_conditions()

template<class ELEMENT >
void oomph::PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement< ELEMENT >::hijack_kinematic_conditions ( const Vector< unsigned > &  bulk_node_number)
inlineprivatevirtual

Hijacking the kinematic condition corresponds to hijacking the spine heights.

Implements oomph::LinearisedAxisymmetricFluidInterfaceElement.

276  {
277  // Loop over all the passed nodes
278  for(Vector<unsigned>::const_iterator it=bulk_node_number.begin();
279  it!=bulk_node_number.end();++it)
280  {
281  // Hijack the spine heights and delete the returned data objects
282  delete this->hijack_nodal_spine_value(*it,0);
283  }
284  }
oomph::FaceElement::bulk_node_number
unsigned & bulk_node_number(const unsigned &n)
Definition: elements.h:4825
oomph::Hijacked< PerturbedSpineElement< FaceGeometry< ELEMENT > > >::hijack_nodal_spine_value
Data * hijack_nodal_spine_value(const unsigned &n, const unsigned &i, const bool &return_data=true)
Definition: hijacked_elements.h:279

References oomph::FaceElement::bulk_node_number(), and oomph::Hijacked< PerturbedSpineElement< FaceGeometry< ELEMENT > > >::hijack_nodal_spine_value().

◆ interpolated_H()

template<class ELEMENT >
double oomph::PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement< ELEMENT >::interpolated_H ( const Vector< double > &  s,
const unsigned &  i 
) const
inline

Return the i-th component of the FE interpolated perturbed surface height (i=0 is cosine part, i=1 is sine part) at local coordinate s

402  {
403  // Determine number of nodes in the element
404  const unsigned n_node = nnode();
405 
406  // Provide storage for local shape functions
407  Shape psi(n_node);
408 
409  // Find values of shape functions
410  shape(s,psi);
411 
412  // Initialise value of H
413  double interpolated_H = 0.0;
414 
415  // Loop over the shape functions and sum
416  for(unsigned l=0;l<n_node;l++)
417  {
418  // Upcast from general node to PerturbedSpineNode
419  PerturbedSpineNode* perturbed_spine_node_pt =
420  dynamic_cast<PerturbedSpineNode*>(this->node_pt(l));
421 
422  // Calculate interpolated i-th component of perturbed spine "heights"
423  interpolated_H +=
424  perturbed_spine_node_pt->perturbed_spine_pt()->height(i)*psi[l];
425  }
426 
427  return(interpolated_H);
428  }
oomph::FiniteElement::shape
virtual void shape(const Vector< double > &s, Shape &psi) const =0
oomph::PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement::interpolated_H
double interpolated_H(const Vector< double > &s, const unsigned &i) const
Definition: linearised_axisymmetric_fluid_interface_elements.h:400

References oomph::PerturbedSpine::height(), i, oomph::FiniteElement::nnode(), oomph::FiniteElement::node_pt(), oomph::PerturbedSpineNode::perturbed_spine_pt(), s, and oomph::FiniteElement::shape().

◆ kinematic_local_eqn()

template<class ELEMENT >
int oomph::PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement< ELEMENT >::kinematic_local_eqn ( const unsigned &  n,
const unsigned &  i 
)
inlineprivatevirtual

In these elements, the kinematic condition is the equation used to determine the "order epsilon" contributions to the free surface "height". We have two sets of unknowns, HC and HS, and we therefore have two kinematic conditions. Overload the function accordingly.

Implements oomph::LinearisedAxisymmetricFluidInterfaceElement.

268  {
269  return this->spine_local_eqn(n,i);
270  }
oomph::PerturbedSpineElement::spine_local_eqn
int spine_local_eqn(const unsigned &n, const unsigned &i)
Definition: perturbed_spines.h:438

References i, and oomph::PerturbedSpineElement< ELEMENT >::spine_local_eqn().

◆ output() [1/4]

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

Overload the C-style output function.

Reimplemented from oomph::LinearisedAxisymmetricFluidInterfaceElement.

383 {FiniteElement::output(file_pt);}
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::PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement< ELEMENT >::output ( FILE *  file_pt,
const unsigned &  n_plot 
)
inlinevirtual

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

Reimplemented from oomph::LinearisedAxisymmetricFluidInterfaceElement.

387  {LinearisedAxisymmetricFluidInterfaceElement::output(file_pt,n_plot);}
oomph::LinearisedAxisymmetricFluidInterfaceElement::output
void output(std::ostream &outfile)
Overload the output functions.
Definition: linearised_axisymmetric_fluid_interface_elements.h:224

References oomph::LinearisedAxisymmetricFluidInterfaceElement::output().

◆ output() [3/4]

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

Overload the output function.

Reimplemented from oomph::LinearisedAxisymmetricFluidInterfaceElement.

376 {FiniteElement::output(outfile);}

References oomph::FiniteElement::output().

◆ output() [4/4]

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

Output the element.

Reimplemented from oomph::LinearisedAxisymmetricFluidInterfaceElement.

380  {LinearisedAxisymmetricFluidInterfaceElement::output(outfile,n_plot);}

References oomph::LinearisedAxisymmetricFluidInterfaceElement::output().

◆ output_interface_position()

template<class ELEMENT >
void oomph::PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement< ELEMENT >::output_interface_position ( std::ostream &  outfile,
const unsigned &  nplot 
)

Output just the interface position (base + perturbation)

Output the combined interface position (base position + perturbation) in tecplot format, using nplot points in each coordinate direction. WARNING: THIS FUNCTION ASSUMES THAT THE PERTURBATION TO THE INTERFACE IS IN A VERTICAL DIRECTION ONLY

172  {
173  // Provide storage for vector of local coordinates
174  Vector<double> s(1);
175 
176  // Tecplot header info
177  outfile << tecplot_zone_string(nplot);
178 
179  // Determine number of plot points
180  const unsigned n_plot_points = nplot_points(nplot);
181 
182  // Loop over plot points
183  for(unsigned iplot=0;iplot<n_plot_points;iplot++)
184  {
185  // Get local coordinates of plot point
186  get_s_plot(iplot,nplot,s);
187 
188  // Output the global r coordinate
189  outfile << interpolated_x(s,0) << " ";
190 
191  // Output the actual interface position (unperturbed+perturbed)
192  // We must ensure that we increase the precision for this!
193  outfile.precision(14);
194  outfile << (interpolated_x(s,1) + interpolated_H(s,0)) << " "
195  << (interpolated_x(s,1) + interpolated_H(s,1)) << " ";
196  outfile.precision(6);
197 
198  outfile << std::endl;
199  }
200  outfile << std::endl;
201 
202  // Write tecplot footer (e.g. FE connectivity lists)
203  write_tecplot_zone_footer(outfile,nplot);
204  }
oomph::FiniteElement::tecplot_zone_string
virtual std::string tecplot_zone_string(const unsigned &nplot) const
Definition: elements.h:3161
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
oomph::FiniteElement::write_tecplot_zone_footer
virtual void write_tecplot_zone_footer(std::ostream &outfile, const unsigned &nplot) const
Definition: elements.h:3174

References s.

Referenced by PerturbedStateProblem< BASE_ELEMENT, PERTURBED_ELEMENT >::doc_solution().

◆ output_perturbation_to_interface()

template<class ELEMENT >
void oomph::PerturbedSpineLinearisedAxisymmetricFluidInterfaceElement< ELEMENT >::output_perturbation_to_interface ( std::ostream &  outfile,
const unsigned &  nplot 
)

Output the perturbation to the interface position.

//////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////// Output the perturbation to the base interface position (computed by these interface elements) in tecplot format, using nplot points in each coordinate direction. WARNING: THIS FUNCTION ASSUMES THAT THE PERTURBATION TO THE INTERFACE IS IN A VERTICAL DIRECTION ONLY

130  {
131  // Provide storage for vector of local coordinates
132  Vector<double> s(1);
133 
134  // Tecplot header info
135  outfile << tecplot_zone_string(nplot);
136 
137  // Determine number of plot points
138  const unsigned n_plot_points = nplot_points(nplot);
139 
140  // Loop over plot points
141  for(unsigned iplot=0;iplot<n_plot_points;iplot++)
142  {
143  // Get local coordinates of plot point
144  get_s_plot(iplot,nplot,s);
145 
146  // Output the global r coordinate
147  outfile << interpolated_x(s,0) << " ";
148 
149  // Output perturbation to interface
150  outfile << interpolated_H(s,0) << " " << interpolated_H(s,1) << " ";
151 
152  outfile << std::endl;
153  }
154  outfile << std::endl;
155 
156  // Write tecplot footer (e.g. FE connectivity lists)
157  write_tecplot_zone_footer(outfile,nplot);
158  }

References s.

Referenced by PerturbedStateProblem< BASE_ELEMENT, PERTURBED_ELEMENT >::doc_solution().

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