Public Member Functions | Protected Member Functions | Private Attributes | Static Private Attributes | List of all members
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT > Class Template Reference
+ Inheritance diagram for oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >:

Public Member Functions

 SpineLineMarangoniSurfactantFluidInterfaceElement (FiniteElement *const &element_pt, const int &face_index)
 
double bi ()
 Return the Biot number. More...
 
double ma ()
 Return the Marangoni number. More...
 
double beta ()
 Return the Elasticity number. More...
 
double peclet_s ()
 Return the surface peclect number. More...
 
double peclet_strouhal_s ()
 Return the surface peclect strouhal number. More...
 
double *& ma_pt ()
 Access function for pointer to the Marangoni number. More...
 
double *& bi_pt ()
 Access function for pointer to the Biot number. More...
 
double *& beta_pt ()
 Access function for pointer to the Elasticity number. More...
 
double *& peclet_s_pt ()
 Access function for pointer to the surface Peclet number. More...
 
double *& peclet_strouhal_s_pt ()
 Access function for pointer to the surface Peclet x Strouhal number. More...
 
void output (std::ostream &outfile, const unsigned &n_plot)
 
 SpineLineMarangoniSurfactantFluidInterfaceElement (FiniteElement *const &element_pt, const int &face_index)
 
double bi ()
 Return the Biot number. More...
 
double ma ()
 Return the Marangoni number. More...
 
double beta ()
 Return the Elasticity number. More...
 
double peclet_s ()
 Return the surface peclect number. More...
 
double peclet_strouhal_s ()
 Return the surface peclect strouhal number. More...
 
double d ()
 Return the diffusion ratio. More...
 
double k ()
 Return the reaction ratio. More...
 
double *& ma_pt ()
 Access function for pointer to the Marangoni number. More...
 
double *& bi_pt ()
 Access function for pointer to the Biot number. More...
 
double *& beta_pt ()
 Access function for pointer to the Elasticity number. More...
 
double *& peclet_s_pt ()
 Access function for pointer to the surface Peclet number. More...
 
double *& peclet_strouhal_s_pt ()
 Access function for pointer to the surface Peclet x Strouhal number. More...
 
double *& d_pt ()
 Access function for pointer to the diffusion ratios. More...
 
double *& k_pt ()
 Access function for pointer to the reaction ratios. More...
 
void output (std::ostream &outfile, const unsigned &n_plot)
 
- Public Member Functions inherited from oomph::SpineLineFluidInterfaceElement< ELEMENT >
 SpineLineFluidInterfaceElement (FiniteElement *const &element_pt, const int &face_index)
 
- Public Member Functions inherited from oomph::SpineUpdateFluidInterfaceElement< FluidInterfaceElement, LineDerivatives, ELEMENT >
 SpineUpdateFluidInterfaceElement (FiniteElement *const &element_pt, const int &face_index, const unsigned &id=0)
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 Calculate the contribution to the residuals and the jacobian. More...
 
void add_additional_residual_contributions_interface (Vector< double > &residuals, DenseMatrix< double > &jacobian, const unsigned &flag, const Shape &psif, const DShape &dpsifds, const DShape &dpsifdS, const DShape &dpsifdS_div, const Vector< double > &s, const Vector< double > &interpolated_x, const Vector< double > &interpolated_n, const double &W, const double &J)
 Helper function to calculate the additional contributions These are those filled in by the particular equations. More...
 
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. More...
 
virtual FluidInterfaceBoundingElement * make_bounding_element (const int &face_index)
 
- Public Member Functions inherited from oomph::Hijacked< SpineElement< 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::SpineElement< FaceGeometry< ELEMENT > >
 SpineElement ()
 Constructor, call the constructor of the base element. More...
 
 SpineElement (FiniteElement *const &element_pt, const int &face_index)
 Constructor used for spine face elements. More...
 
 ~SpineElement ()
 Destructor, clean up the storage allocated to the local equation numbers. More...
 
int spine_local_eqn (const unsigned &n)
 
- Public Member Functions inherited from oomph::ElementWithSpecificMovingNodes< ELEMENT, NODE_TYPE >
void describe_local_dofs (std::ostream &out, const std::string &current_string) const
 
 ElementWithSpecificMovingNodes ()
 Constructor, call the constructor of the base element. More...
 
 ElementWithSpecificMovingNodes (FiniteElement *const &element_pt, const int &face_index)
 Constructor used for face elements. More...
 
 ~ElementWithSpecificMovingNodes ()
 Empty Destructor,. More...
 
void describe_local_dofs (std::ostream &out, std::string &curr_str)
 Unique final overrider for describe_dofs. More...
 
Node * construct_node (const unsigned &n)
 
Node * construct_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 Overloaded node allocation for unsteady problems. More...
 
Node * construct_boundary_node (const unsigned &n)
 Overload the node assignment routine to assign boundary nodes. More...
 
Node * construct_boundary_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 Overloaded boundary node allocation for unsteady problems. More...
 
void assign_all_generic_local_eqn_numbers (const bool &store_local_dof_pt)
 
void get_jacobian_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::SpineFiniteElement
 SpineFiniteElement ()
 Empty constructor. More...
 
virtual ~SpineFiniteElement ()
 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::FluidInterfaceElement
 FluidInterfaceElement ()
 Constructor, set the default values of the booleans and pointers (null) More...
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Calculate the residuals by calling the generic residual contribution. More...
 
const double & ca () const
 The value of the Capillary number. More...
 
double *& ca_pt ()
 Pointer to the Capillary number. More...
 
const double & st () const
 The value of the Strouhal number. More...
 
double *& st_pt ()
 The pointer to the Strouhal number. More...
 
double u (const unsigned &j, const unsigned &i)
 Return the i-th velocity component at local node j. More...
 
double interpolated_u (const Vector< double > &s, const unsigned &i)
 Calculate the i-th velocity component at the local coordinate s. More...
 
double pext () const
 Return the value of the external pressure. More...
 
void set_external_pressure_data (Data *external_pressure_data_pt)
 
void set_external_pressure_data (Data *external_pressure_data_pt, const unsigned &index_of_external_pressure_value)
 
- 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::LineDerivatives
 LineDerivatives ()
 

Protected Member Functions

double interpolated_T (const Vector< double > &s)
 Get the temperature. More...
 
double interpolated_C (const Vector< double > &s)
 Get the surfactant concentration. More...
 
double dcdt_surface (const unsigned &l) const
 The time derivative of the surface concentration. More...
 
double sigma (const Vector< double > &s)
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 Calculate the contribution to the jacobian. More...
 
void add_additional_residual_contributions_interface (Vector< double > &residuals, DenseMatrix< double > &jacobian, const unsigned &flag, const Shape &psif, const DShape &dpsifds, const DShape &dpsifdS, const DShape &dpsifds_div, const Vector< double > &s, const Vector< double > &interpolated_x, const Vector< double > &interpolated_n, const double &W, const double &J)
 
void fill_in_contribution_to_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
double interpolated_T (const Vector< double > &s)
 Get the temperature. More...
 
double interpolated_C (const Vector< double > &s)
 Get the surfactant concentration. More...
 
double dcdt_surface (const unsigned &l) const
 The time derivative of the surface concentration. More...
 
double sigma (const Vector< double > &s)
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 Calculate the contribution to the jacobian. More...
 
void add_additional_residual_contributions_interface (Vector< double > &residuals, DenseMatrix< double > &jacobian, const unsigned &flag, const Shape &psif, const DShape &dpsifds, const DShape &dpsifdS, const DShape &dpsifds_div, const Vector< double > &s, const Vector< double > &interpolated_x, const Vector< double > &interpolated_n, const double &W, const double &J)
 
void fill_in_contribution_to_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
- Protected Member Functions inherited from oomph::SpineUpdateFluidInterfaceElement< FluidInterfaceElement, LineDerivatives, ELEMENT >
double compute_surface_derivatives (const Shape &psi, const DShape &dpsids, const DenseMatrix< double > &interpolated_t, const Vector< double > &interpolated_x, DShape &surface_gradient, DShape &surface_divergence)
 
- 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)
 
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_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::FluidInterfaceElement
int pext_local_eqn ()
 
virtual void fill_in_generic_residual_contribution_interface (Vector< double > &residuals, DenseMatrix< double > &jacobian, unsigned flag)
 
- 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::LineDerivatives
double compute_surface_derivatives (const Shape &psi, const DShape &dpsids, const DenseMatrix< double > &interpolated_t, const Vector< double > &interpolated_x, DShape &surface_gradient, DShape &surface_divergence)
 Fill in the specific surface derivative calculations. More...
 

Private Attributes

double * Bi_pt
 Pointer to a Biot number. More...
 
double * Ma_pt
 Pointer to a Marangoni number. More...
 
double * Beta_pt
 Pointer to an Elasticity number. More...
 
double * Peclet_S_pt
 Pointer to Surface Peclet number. More...
 
double * Peclet_Strouhal_S_pt
 Pointer to the surface Peclect Strouhal number. More...
 
unsigned T_index
 Index at which the temperature is stored at the nodes. More...
 
unsigned C_index
 
double * D_pt
 Pointer to the diffusion ratios. More...
 
double * K_pt
 Pointer to the reaction ratios. More...
 
unsigned C_bulk_index
 Index at which the bulk concentration is stored at the nodes. More...
 

Static Private Attributes

static double Default_Physical_Constant_Value = 1.0
 Default value of the physical constants. More...
 

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 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::FluidInterfaceElement
Vector< unsigned > U_index_interface
 Nodal index at which the i-th velocity component is stored. More...
 
int External_data_number_of_external_pressure
 
Data * Pext_data_pt
 Pointer to the Data item that stores the external pressure. More...
 
unsigned Index_of_external_pressure_value
 Which of the values in Pext_data_pt stores the external pressure. More...
 
- 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::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >

Spine-based Marangoni surface tension elements that add a linear dependence on temperature and concentration of a surface chemical to the surface tension, which decreases with increasing temperature. In addition, this element adds a flux contribution to the advection-diffusion equation to represent heat loss at the free surface. This introduces the Biot number.

Constructor & Destructor Documentation

◆ SpineLineMarangoniSurfactantFluidInterfaceElement() [1/2]

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

Constructor that passes the bulk element and face index down to the underlying

490  :
491  SpineLineFluidInterfaceElement<ELEMENT>
492  (element_pt,face_index)
493  {
494  //Initialise the values
495  Bi_pt = &Default_Physical_Constant_Value;
496  Ma_pt = &Default_Physical_Constant_Value;
497  Beta_pt = &Default_Physical_Constant_Value;
498  Peclet_S_pt = &Default_Physical_Constant_Value;
499  Peclet_Strouhal_S_pt = &Default_Physical_Constant_Value;
500 
501  //Cast the bulk element
502  ELEMENT* cast_element_pt = dynamic_cast<ELEMENT*>(element_pt);
503  //Now find the index at which the temperature is stored from the
504  //advection-diffusion part of the bulk element
505  T_index = cast_element_pt->u_index_adv_diff();
506 
507  //Add the additional surfactant terms to these surface elements
508 
509  //Read out the number of nodes on the face
510  //For some reason I need to specify the this pointer here(!)
511  unsigned n_node_face = this->nnode();
512  //Set the additional data values in the face
513  //There is one additional values at each node --- the lagrange multiplier
514  Vector<unsigned> additional_data_values(n_node_face);
515  for(unsigned i=0;i<n_node_face;i++) additional_data_values[i] = 1;
516  //Resize the data arrays accordingly
517  this->resize_nodes(additional_data_values);
518 
519  //The C_index is the new final value
520  //HACK HERE
521  C_index = this->node_pt(0)->nvalue()-1;
522  }
i
int i
Definition: BiCGSTAB_step_by_step.cpp:9
oomph::Data::nvalue
unsigned nvalue() const
Return number of values stored in data object (incl pinned ones).
Definition: nodes.h:483
oomph::FaceElement::face_index
int & face_index()
Definition: elements.h:4626
oomph::FaceElement::resize_nodes
void resize_nodes(Vector< unsigned > &nadditional_data_values)
Definition: elements.h:4882
oomph::FiniteElement::node_pt
Node *& node_pt(const unsigned &n)
Return a pointer to the local node n.
Definition: elements.h:2175
oomph::FiniteElement::nnode
unsigned nnode() const
Return the number of nodes.
Definition: elements.h:2210
oomph::SpineLineFluidInterfaceElement
Definition: specific_node_update_interface_elements.h:556
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::Peclet_S_pt
double * Peclet_S_pt
Pointer to Surface Peclet number.
Definition: surfactant.cc:80
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::Peclet_Strouhal_S_pt
double * Peclet_Strouhal_S_pt
Pointer to the surface Peclect Strouhal number.
Definition: surfactant.cc:83
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::Ma_pt
double * Ma_pt
Pointer to a Marangoni number.
Definition: surfactant.cc:74
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::T_index
unsigned T_index
Index at which the temperature is stored at the nodes.
Definition: surfactant.cc:86
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::C_index
unsigned C_index
Definition: surfactant.cc:90
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::Bi_pt
double * Bi_pt
Pointer to a Biot number.
Definition: surfactant.cc:71
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::Default_Physical_Constant_Value
static double Default_Physical_Constant_Value
Default value of the physical constants.
Definition: surfactant.cc:93
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::Beta_pt
double * Beta_pt
Pointer to an Elasticity number.
Definition: surfactant.cc:77

References oomph::MultiDomainBoussinesqHelper::Default_Physical_Constant_Value, and i.

◆ SpineLineMarangoniSurfactantFluidInterfaceElement() [2/2]

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

Constructor that passes the bulk element and face index down to the underlying

580  :
581  SpineLineFluidInterfaceElement<ELEMENT>
582  (element_pt,face_index)
583  {
584  //Initialise the values
585  Bi_pt = &Default_Physical_Constant_Value;
586  Ma_pt = &Default_Physical_Constant_Value;
587  Beta_pt = &Default_Physical_Constant_Value;
588  Peclet_S_pt = &Default_Physical_Constant_Value;
589  Peclet_Strouhal_S_pt = &Default_Physical_Constant_Value;
590  K_pt = &Default_Physical_Constant_Value;
591  D_pt = &Default_Physical_Constant_Value;
592 
593  //Cast the bulk element
594  ELEMENT* cast_element_pt = dynamic_cast<ELEMENT*>(element_pt);
595  //Now find the index at which the temperature is stored from the
596  //advection-diffusion part of the bulk element
597  T_index = cast_element_pt->c_index_adv_diff_react(0);
598 
599  //Find the index at which the bulk surfactant concentration is stored
600  C_bulk_index = cast_element_pt->c_index_adv_diff_react(1);
601 
602  //Add the additional surfactant terms to these surface elements
603 
604  //Read out the number of nodes on the face
605  //For some reason I need to specify the this pointer here(!)
606  unsigned n_node_face = this->nnode();
607  //Set the additional data values in the face
608  //There is one additional values at each node --- the lagrange multiplier
609  Vector<unsigned> additional_data_values(n_node_face);
610  for(unsigned i=0;i<n_node_face;i++) additional_data_values[i] = 1;
611  //Resize the data arrays accordingly
612  this->resize_nodes(additional_data_values);
613 
614  //The C_index is the new final value
615  //HACK HERE
616  C_index = this->node_pt(0)->nvalue()-1;
617  }
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::C_bulk_index
unsigned C_bulk_index
Index at which the bulk concentration is stored at the nodes.
Definition: soluble_surfactant.cc:95
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::K_pt
double * K_pt
Pointer to the reaction ratios.
Definition: soluble_surfactant.cc:89
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::D_pt
double * D_pt
Pointer to the diffusion ratios.
Definition: soluble_surfactant.cc:86

References oomph::MultiDomainBoussinesqHelper::Default_Physical_Constant_Value, and i.

Member Function Documentation

◆ add_additional_residual_contributions_interface() [1/2]

template<class ELEMENT >
void oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::add_additional_residual_contributions_interface ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian,
const unsigned &  flag,
const Shape &  psif,
const DShape &  dpsifds,
const DShape &  dpsifdS,
const DShape &  dpsifds_div,
const Vector< double > &  s,
const Vector< double > &  interpolated_x,
const Vector< double > &  interpolated_n,
const double &  W,
const double &  J 
)
inlineprotectedvirtual

Overload the Helper function to calculate the residuals and jacobian entries. This particular function ensures that the additional entries are calculated inside the integration loop

Reimplemented from oomph::FluidInterfaceElement.

321  {
322  //Find the index at which the temperature is stored
323  unsigned t_index = this->T_index;
324 
325  //Find out how many nodes there are
326  unsigned n_node = this->nnode();
327 
328  //Read out the Bi number
329  const double Bi = this->bi();
330 
331  //Now calculate the temperature at this point
332  //Assuming the same shape functions are used (which they are)
333  double T = 0.0;
334  for(unsigned l=0;l<n_node;l++)
335  {
336  T += this->nodal_value(l,t_index)*psif(l);
337  }
338 
339  //Storage for the local equation numbers and unknowns
340  int local_eqn = 0, local_unknown = 0;
341 
342  //Now we add the flux term to the appropriate residuals
343  for(unsigned l=0;l<n_node;l++)
344  {
345  //Read out the apprporiate local equation
346  local_eqn = this->nodal_local_eqn(l,t_index);
347 
348  //If not a boundary condition
349  if(local_eqn >= 0)
350  {
351  residuals[local_eqn] -= (1.0 + Bi*(T-1.0))*psif(l)*W*J;
352 
353  //We also need to worry about the jacobian terms
354  if(flag)
355  {
356  //Loop over the nodes again
357  for(unsigned l2=0;l2<n_node;l2++)
358  {
359  //Get the unknown
360  local_unknown = this->nodal_local_eqn(l2,t_index);
361  //If not a boundary condition
362  if(local_unknown >= 0)
363  {
364  jacobian(local_eqn,local_unknown) -= Bi*psif(l2)*psif(l)*W*J;
365  }
366  }
367  }
368  }
369  } //End of loop over the nodes
370 
371 
372  //Surface advection-diffusion equation
373 
374  //Find the index at which the temperature is stored
375  unsigned c_index = this->C_index;
376  Vector<unsigned> u_index = this->U_index_interface;
377 
378  //Read out the surface peclect number
379  const double Pe_s = this->peclet_s();
380  const double PeSt_s = this->peclet_strouhal_s();
381 
382  //Now calculate the concentration at this point
383  //Assuming the same shape functions are used (which they are)
384  double interpolated_C = 0.0;
385  double interpolated_dCds = 0.0;
386  double dCdt = 0.0;
387  //The tangent vector
388  const unsigned ndim = this->node_pt(0)->ndim();
389  Vector<double> interpolated_tangent(ndim,0.0);
390  Vector<double> interpolated_u(ndim,0.0);
391  Vector<double> interpolated_duds(ndim,0.0);
392  Vector<double> mesh_velocity(ndim,0.0);
393 
394 
395  if(ndim != u_index.size())
396  {
397  throw OomphLibError("Dimension Incompatibility",
398  OOMPH_CURRENT_FUNCTION,
399  OOMPH_EXCEPTION_LOCATION);
400  }
401 
402  for(unsigned l=0;l<n_node;l++)
403  {
404  const double psi = psif(l);
405  const double dpsi = dpsifds(l,0);
406  interpolated_C += this->nodal_value(l,c_index)*psi;
407  interpolated_dCds += this->nodal_value(l,c_index)*dpsi;
408  dCdt += dcdt_surface(l)*psi;
409  for(unsigned i=0;i<ndim;i++)
410  {
411  interpolated_tangent[i] += this->nodal_position(l,i)*dpsi;
412  interpolated_u[i] += this->nodal_value(l,u_index[i])*psi;
413  interpolated_duds[i] += this->nodal_value(l,u_index[i])*dpsi;
414  mesh_velocity[i] += this->dnodal_position_dt(l,i)*psi;
415  }
416  }
417 
418  double u_tangent = 0.0, t_length = 0.0;
419  for(unsigned i=0;i<ndim;i++)
420  {
421  u_tangent += (interpolated_u[i] - mesh_velocity[i])*interpolated_tangent[i];
422  t_length += interpolated_tangent[i]*interpolated_tangent[i];
423  }
424 
425  //Now we add the flux term to the appropriate residuals
426  for(unsigned l=0;l<n_node;l++)
427  {
428  //Read out the apprporiate local equation
429  local_eqn = this->nodal_local_eqn(l,c_index);
430 
431  //If not a boundary condition
432  if(local_eqn >= 0)
433  {
434  //Time derivative term
435  residuals[local_eqn] -= PeSt_s*dCdt*psif(l)*W*J;
436 
437  //Diffusion term
438  residuals[local_eqn] -= interpolated_dCds*dpsifds(l,0)*W/J;
439 
440  //Advection term in new formulation
441  residuals[local_eqn] -=
442  Pe_s*(u_tangent*interpolated_dCds +
443  interpolated_C*
444  (interpolated_tangent[0]*interpolated_duds[0]
445  + interpolated_tangent[1]*interpolated_duds[1]))*psif(l)*W/J;
446 
447  //We also need to worry about the jacobian terms
448  if(flag)
449  {
450  //Loop over the nodes again
451  for(unsigned l2=0;l2<n_node;l2++)
452  {
453  //Loop over the velocity components
454  for(unsigned i2=0;i2<ndim;i2++)
455  {
456  //Get the unknown
457  local_unknown = this->nodal_local_eqn(l2,u_index[i2]);
458  //If not a boundary condition
459  if(local_unknown >= 0)
460  {
461  jacobian(local_eqn,local_unknown)
462  -= Pe_s*interpolated_tangent[i2]*
463  (psif(l2)*interpolated_dCds +
464  interpolated_C*dpsifds(l2,0))*psif(l)*W/J;
465  }
466  }
467  }
468  }
469  }
470  } //End of loop over the nodes
471 
472  }
J
JacobiRotation< float > J
Definition: Jacobi_makeJacobi.cpp:3
oomph::FiniteElement::nodal_value
double nodal_value(const unsigned &n, const unsigned &i) const
Definition: elements.h:2593
oomph::FiniteElement::nodal_local_eqn
int nodal_local_eqn(const unsigned &n, const unsigned &i) const
Definition: elements.h:1432
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::FluidInterfaceElement::U_index_interface
Vector< unsigned > U_index_interface
Nodal index at which the i-th velocity component is stored.
Definition: interface_elements.h:338
oomph::FluidInterfaceElement::interpolated_u
double interpolated_u(const Vector< double > &s, const unsigned &i)
Calculate the i-th velocity component at the local coordinate s.
Definition: interface_elements.cc:442
oomph::GeomObject::ndim
unsigned ndim() const
Access function to # of Eulerian coordinates.
Definition: geom_objects.h:177
oomph::Node::ndim
unsigned ndim() const
Return (Eulerian) spatial dimension of the node.
Definition: nodes.h:1054
oomph::OomphLibError
Definition: oomph_definitions.h:222
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::peclet_strouhal_s
double peclet_strouhal_s()
Return the surface peclect strouhal number.
Definition: surfactant.cc:537
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::interpolated_C
double interpolated_C(const Vector< double > &s)
Get the surfactant concentration.
Definition: surfactant.cc:126
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::dcdt_surface
double dcdt_surface(const unsigned &l) const
The time derivative of the surface concentration.
Definition: surfactant.cc:154
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::peclet_s
double peclet_s()
Return the surface peclect number.
Definition: surfactant.cc:534
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::bi
double bi()
Return the Biot number.
Definition: surfactant.cc:525
Global_Physical_Variables::Bi
double Bi
Biot number.
Definition: thermal_fibre.cc:85
oomph::Global_Physical_Variables::Pe_s
double Pe_s
Surface Peclet number.
Definition: refineable_two_layer_soluble_surfactant.cc:133
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

References Global_Physical_Variables::Bi, i, J, OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, oomph::Global_Physical_Variables::Pe_s, and oomph::QuadTreeNames::W.

◆ add_additional_residual_contributions_interface() [2/2]

template<class ELEMENT >
void oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::add_additional_residual_contributions_interface ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian,
const unsigned &  flag,
const Shape &  psif,
const DShape &  dpsifds,
const DShape &  dpsifdS,
const DShape &  dpsifds_div,
const Vector< double > &  s,
const Vector< double > &  interpolated_x,
const Vector< double > &  interpolated_n,
const double &  W,
const double &  J 
)
inlineprotectedvirtual

Overload the Helper function to calculate the residuals and jacobian entries. This particular function ensures that the additional entries are calculated inside the integration loop

Reimplemented from oomph::FluidInterfaceElement.

330  {
331  //Find the index at which the temperature is stored
332  unsigned t_index = this->T_index;
333 
334  //Find out how many nodes there are
335  unsigned n_node = this->nnode();
336 
337  //Read out the Bi number
338  const double Bi = this->bi();
339 
340  //Now calculate the temperature at this point
341  //Assuming the same shape functions are used (which they are)
342  double T = 0.0;
343  for(unsigned l=0;l<n_node;l++)
344  {
345  T += this->nodal_value(l,t_index)*psif(l);
346  }
347 
348  //Storage for the local equation numbers and unknowns
349  int local_eqn = 0, local_unknown = 0;
350 
351  //Now we add the flux term to the appropriate residuals
352  for(unsigned l=0;l<n_node;l++)
353  {
354  //Read out the apprporiate local equation
355  local_eqn = this->nodal_local_eqn(l,t_index);
356 
357  //If not a boundary condition
358  if(local_eqn >= 0)
359  {
360  residuals[local_eqn] -= (1.0 + Bi*(T-1.0))*psif(l)*W*J;
361 
362  //We also need to worry about the jacobian terms
363  if(flag)
364  {
365  //Loop over the nodes again
366  for(unsigned l2=0;l2<n_node;l2++)
367  {
368  //Get the unknown
369  local_unknown = this->nodal_local_eqn(l2,t_index);
370  //If not a boundary condition
371  if(local_unknown >= 0)
372  {
373  jacobian(local_eqn,local_unknown) -= Bi*psif(l2)*psif(l)*W*J;
374  }
375  }
376  }
377  }
378  } //End of loop over the nodes
379 
380 
381  //Bulk flux condition
382  //Find the index at which the bulk is stored
383  unsigned c_bulk_index = this->C_bulk_index;
384 
385  //Find the index at which the surface surfactant
386  unsigned c_index = this->C_index;
387 
388 
389  //Now calculate the bulk concentration at this point
390  //Assuming the same shape functions are used (which they are)
391  double C_bulk = 0.0;
392  double C = 0.0;
393  for(unsigned l=0;l<n_node;l++)
394  {
395  C_bulk += this->nodal_value(l,c_bulk_index)*psif(l);
396  C += this->nodal_value(l,c_index)*psif(l);
397  }
398 
399  //Get the reaction ratio
400  const double K = this->k();
401 
402  //The transport between the two layers is given by the flux
403  double flux = K*C_bulk - C;
404 
405  //Now we add the flux term to the appropriate residuals
406  for(unsigned l=0;l<n_node;l++)
407  {
408  //Read out the apprporiate local equation
409  local_eqn = this->nodal_local_eqn(l,c_bulk_index);
410 
411  //If not a boundary condition
412  if(local_eqn >= 0)
413  {
414  residuals[local_eqn] -= flux*psif(l)*W*J;
415 
416  //We also need to worry about the jacobian terms
417  if(flag)
418  {
419  //Loop over the nodes again
420  for(unsigned l2=0;l2<n_node;l2++)
421  {
422  //Get the unknown
423  local_unknown = this->nodal_local_eqn(l2,c_bulk_index);
424  //If not a boundary condition
425  if(local_unknown >= 0)
426  {
427  jacobian(local_eqn,local_unknown) -= K*psif(l2)*psif(l)*W*J;
428  }
429 
430  local_unknown = this->nodal_local_eqn(l2,c_index);
431  //If not a boundary condition
432  if(local_unknown >= 0)
433  {
434  jacobian(local_eqn,local_unknown) += psif(l2)*psif(l)*W*J;
435  }
436  }
437  }
438  }
439  } //End of loop over the nodes
440 
441 
442 
443  //Surface advection-diffusion equation
444  //Get the diffusion ratio
445  const double D = d();
446 
447  //Find the index at which the temperature is stored
448  //unsigned c_index = this->C_index;
449  Vector<unsigned> u_index = this->U_index_interface;
450 
451  //Read out the surface peclect number
452  const double Pe_s = this->peclet_s();
453  const double PeSt_s = this->peclet_strouhal_s();
454 
455  //Now calculate the concentration at this point
456  //Assuming the same shape functions are used (which they are)
457  double interpolated_C = 0.0;
458  double interpolated_dCds = 0.0;
459  double dCdt = 0.0;
460  //The tangent vector
461  const unsigned ndim = this->node_pt(0)->ndim();
462  Vector<double> interpolated_tangent(ndim,0.0);
463  Vector<double> interpolated_u(ndim,0.0);
464  Vector<double> interpolated_duds(ndim,0.0);
465  Vector<double> mesh_velocity(ndim,0.0);
466 
467  if(ndim != u_index.size())
468  {
469  throw OomphLibError("Dimension Incompatibility",
470  OOMPH_CURRENT_FUNCTION,
471  OOMPH_EXCEPTION_LOCATION);
472  }
473 
474  for(unsigned l=0;l<n_node;l++)
475  {
476  const double psi = psif(l);
477  const double dpsi = dpsifds(l,0);
478  interpolated_C += this->nodal_value(l,c_index)*psi;
479  interpolated_dCds += this->nodal_value(l,c_index)*dpsi;
480  dCdt += dcdt_surface(l)*psi;
481  for(unsigned i=0;i<ndim;i++)
482  {
483  interpolated_tangent[i] += this->nodal_position(l,i)*dpsi;
484  interpolated_u[i] += this->nodal_value(l,u_index[i])*psi;
485  interpolated_duds[i] += this->nodal_value(l,u_index[i])*dpsi;
486  mesh_velocity[i] += this->dnodal_position_dt(l,i)*psif(l);
487  }
488  }
489 
490  double u_tangent = 0.0, t_length = 0.0;
491  for(unsigned i=0;i<ndim;i++)
492  {
493  u_tangent +=
494  (interpolated_u[i] - mesh_velocity[i])*interpolated_tangent[i];
495  t_length += interpolated_tangent[i]*interpolated_tangent[i];
496  }
497 
498  //Now we add the flux term to the appropriate residuals
499  for(unsigned l=0;l<n_node;l++)
500  {
501  //Read out the apprporiate local equation
502  local_eqn = this->nodal_local_eqn(l,c_index);
503 
504  //If not a boundary condition
505  if(local_eqn >= 0)
506  {
507  //Time derivative term
508  residuals[local_eqn] -= PeSt_s*dCdt*psif(l)*W*J;
509 
510  //Diffusion term
511  residuals[local_eqn] -= interpolated_dCds*dpsifds(l,0)*W/J;
512 
513  //Advection term in new formulation
514  residuals[local_eqn] -=
515  Pe_s*(u_tangent*interpolated_dCds +
516  interpolated_C*
517  (interpolated_tangent[0]*interpolated_duds[0]
518  + interpolated_tangent[1]*interpolated_duds[1]))*psif(l)*W/J;
519 
520  //Now add the flux term
521  residuals[local_eqn] += D*flux*psif(l)*W*J;
522  //We also need to worry about the jacobian terms
523  if(flag)
524  {
525  //Loop over the nodes again
526  for(unsigned l2=0;l2<n_node;l2++)
527  {
528  //Loop over the velocity components
529  for(unsigned i2=0;i2<ndim;i2++)
530  {
531  //Get the unknown
532  local_unknown = this->nodal_local_eqn(l2,u_index[i2]);
533  //If not a boundary condition
534  if(local_unknown >= 0)
535  {
536  jacobian(local_eqn,local_unknown)
537  -= Pe_s*interpolated_tangent[i2]*
538  (psif(l2)*interpolated_dCds +
539  interpolated_C*dpsifds(l2,0))*psif(l)*W/J;
540  }
541  }
542 
543  //Get the unknown
544  local_unknown = this->nodal_local_eqn(l2,c_bulk_index);
545  //If not a boundary condition
546  if(local_unknown >= 0)
547  {
548  jacobian(local_eqn,local_unknown) += D*K*psif(l2)*psif(l)*W*J;
549  }
550 
551  local_unknown = this->nodal_local_eqn(l2,c_index);
552  //If not a boundary condition
553  if(local_unknown >= 0)
554  {
555  jacobian(local_eqn,local_unknown) -= D*psif(l2)*psif(l)*W*J;
556  }
557  }
558  }
559  }
560  } //End of loop over the nodes
561 
562  }
C
Matrix< Scalar, Dynamic, Dynamic > C
Definition: bench_gemm.cpp:49
oomph::Matrix
Definition: matrices.h:74
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::d
double d()
Return the diffusion ratio.
Definition: soluble_surfactant.cc:635
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::k
double k()
Return the reaction ratio.
Definition: soluble_surfactant.cc:638
Global_Physical_Variables::D
Vector< double > D(2, 1.0)
Diffusivity (identically one from our non-dimensionalisation)
PlanarWave::K
double K
Wave number.
Definition: sphere_scattering.cc:115
ProblemParameters::flux
void flux(const double &time, const Vector< double > &x, double &flux)
Get flux applied along boundary x=0.
Definition: pretend_melt.cc:59

References Global_Physical_Variables::Bi, Global_Physical_Variables::D, ProblemParameters::flux(), i, J, k, PlanarWave::K, OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, oomph::Global_Physical_Variables::Pe_s, and oomph::QuadTreeNames::W.

◆ beta() [1/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::beta ( )
inline

Return the Elasticity number.

531 {return *Beta_pt;}

◆ beta() [2/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::beta ( )
inline

Return the Elasticity number.

626 {return *Beta_pt;}

◆ beta_pt() [1/2]

template<class ELEMENT >
double* & oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::beta_pt ( )
inline

Access function for pointer to the Elasticity number.

546 {return Beta_pt;}

◆ beta_pt() [2/2]

template<class ELEMENT >
double* & oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::beta_pt ( )
inline

Access function for pointer to the Elasticity number.

647 {return Beta_pt;}

◆ bi() [1/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::bi ( )
inline

Return the Biot number.

525 {return *Bi_pt;}

◆ bi() [2/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::bi ( )
inline

Return the Biot number.

620 {return *Bi_pt;}

◆ bi_pt() [1/2]

template<class ELEMENT >
double* & oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::bi_pt ( )
inline

Access function for pointer to the Biot number.

543 {return Bi_pt;}

◆ bi_pt() [2/2]

template<class ELEMENT >
double* & oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::bi_pt ( )
inline

Access function for pointer to the Biot number.

644 {return Bi_pt;}

◆ d()

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::d ( )
inline

Return the diffusion ratio.

635 {return *D_pt;}

◆ d_pt()

template<class ELEMENT >
double* & oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::d_pt ( )
inline

Access function for pointer to the diffusion ratios.

656 {return D_pt;}

◆ dcdt_surface() [1/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::dcdt_surface ( const unsigned &  l) const
inlineprotected

The time derivative of the surface concentration.

155  {
156  // Get the data's timestepper
157  TimeStepper* time_stepper_pt= this->node_pt(l)->time_stepper_pt();
158 
159  //Initialise dudt
160  double dcdt=0.0;
161  //Loop over the timesteps, if there is a non Steady timestepper
162  if (time_stepper_pt->type()!="Steady")
163  {
164  //Find the index at which the variable is stored
165  const unsigned c_index = C_index;
166 
167  // Number of timsteps (past & present)
168  const unsigned n_time = time_stepper_pt->ntstorage();
169 
170  for(unsigned t=0;t<n_time;t++)
171  {
172  dcdt += time_stepper_pt->weight(1,t)*this->nodal_value(t,l,c_index);
173  }
174  }
175  return dcdt;
176  }
oomph::Data::time_stepper_pt
TimeStepper *& time_stepper_pt()
Return the pointer to the timestepper.
Definition: nodes.h:238
oomph::GeomObject::time_stepper_pt
TimeStepper *& time_stepper_pt()
Definition: geom_objects.h:192
oomph::TimeStepper
Definition: timesteppers.h:231
oomph::TimeStepper::ntstorage
unsigned ntstorage() const
Definition: timesteppers.h:601
oomph::TimeStepper::weight
virtual double weight(const unsigned &i, const unsigned &j) const
Access function for j-th weight for the i-th derivative.
Definition: timesteppers.h:594
oomph::TimeStepper::type
std::string type() const
Definition: timesteppers.h:490
plotPSD.t
t
Definition: plotPSD.py:36

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

◆ dcdt_surface() [2/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::dcdt_surface ( const unsigned &  l) const
inlineprotected

The time derivative of the surface concentration.

164  {
165  // Get the data's timestepper
166  TimeStepper* time_stepper_pt= this->node_pt(l)->time_stepper_pt();
167 
168  //Initialise dudt
169  double dcdt=0.0;
170  //Loop over the timesteps, if there is a non Steady timestepper
171  if (time_stepper_pt->type()!="Steady")
172  {
173  //Find the index at which the variable is stored
174  const unsigned c_index = C_index;
175 
176  // Number of timsteps (past & present)
177  const unsigned n_time = time_stepper_pt->ntstorage();
178 
179  for(unsigned t=0;t<n_time;t++)
180  {
181  dcdt += time_stepper_pt->weight(1,t)*this->nodal_value(t,l,c_index);
182  }
183  }
184  return dcdt;
185  }

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

◆ fill_in_contribution_to_jacobian() [1/2]

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

Calculate the contribution to the jacobian.

Fill in the contribution to the residuals

Reimplemented from oomph::GeneralisedElement.

211  {
212  //Call the generic routine with the flag set to 1
213  this->fill_in_generic_residual_contribution_interface(residuals,jacobian,1);
214  {
215  //Use finite differences to handle temperature variations
216  const unsigned n_node = this->nnode();
217  //Find the number of dofs in the element
218  const unsigned n_dof = this->ndof();
219  //Create newres vector
220  Vector<double> newres(n_dof);
221 
222  //Integer storage for local unknown
223  int local_unknown=0;
224 
225  //Use the default finite difference step
226  const double fd_step = this->Default_fd_jacobian_step;
227 
228  //Loop over the nodes
229  for(unsigned n=0;n<n_node;n++)
230  {
231  //Get the number of values stored at the node
232  unsigned t_index = this->T_index;
233 
234  //Get the local equation number
235  local_unknown = this->nodal_local_eqn(n,t_index);
236  //If it's not pinned
237  if(local_unknown >= 0)
238  {
239  //Store a pointer to the nodal data value
240  double *value_pt = this->node_pt(n)->value_pt(t_index);
241 
242  //Save the old value of the Nodal data
243  double old_var = *value_pt;
244 
245  //Increment the value of the Nodal data
246  *value_pt += fd_step;
247 
248  //Calculate the new residuals
249  this->get_residuals(newres);
250 
251  //Do finite differences
252  for(unsigned m=0;m<n_dof;m++)
253  {
254  double sum = (newres[m] - residuals[m])/fd_step;
255  //Stick the entry into the Jacobian matrix
256  jacobian(m,local_unknown) = sum;
257  }
258 
259  //Reset the Nodal data
260  *value_pt = old_var;
261  }
262  }
263 
264  //Use finite differences to handle the concentration variations
265  //Loop over the nodes again
266  for(unsigned n=0;n<n_node;n++)
267  {
268  //Get the number of values stored at the node
269  unsigned c_index = this->C_index;
270 
271  //Get the local equation number
272  local_unknown = this->nodal_local_eqn(n,c_index);
273  //If it's not pinned
274  if(local_unknown >= 0)
275  {
276  //Store a pointer to the nodal data value
277  double *value_pt = this->node_pt(n)->value_pt(c_index);
278 
279  //Save the old value of the Nodal data
280  double old_var = *value_pt;
281 
282  //Increment the value of the Nodal data
283  *value_pt += fd_step;
284 
285  //Calculate the new residuals
286  this->get_residuals(newres);
287 
288  //Do finite differences
289  for(unsigned m=0;m<n_dof;m++)
290  {
291  double sum = (newres[m] - residuals[m])/fd_step;
292  //Stick the entry into the Jacobian matrix
293  jacobian(m,local_unknown) = sum;
294  }
295 
296  //Reset the Nodal data
297  *value_pt = old_var;
298  }
299  }
300  }
301 
302  //Call the generic routine to handle the spine variables
303  SpineElement<FaceGeometry<ELEMENT> >::
304  fill_in_jacobian_from_geometric_data(jacobian);
305  }
n
const unsigned n
Definition: CG3DPackingUnitTest.cpp:11
oomph::Data::value_pt
double * value_pt(const unsigned &i) const
Definition: nodes.h:324
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::FluidInterfaceElement::fill_in_generic_residual_contribution_interface
virtual void fill_in_generic_residual_contribution_interface(Vector< double > &residuals, DenseMatrix< double > &jacobian, unsigned flag)
Definition: interface_elements.cc:471
oomph::GeneralisedElement::Default_fd_jacobian_step
static double Default_fd_jacobian_step
Definition: elements.h:1198
oomph::GeneralisedElement::ndof
unsigned ndof() const
Return the number of equations/dofs in the element.
Definition: elements.h:835
oomph::Hijacked< SpineElement< FaceGeometry< ELEMENT > > >::get_residuals
void get_residuals(Vector< double > &residuals)
Definition: hijacked_elements.h:377
oomph::SpineElement
Definition: spines.h:477
m
int * m
Definition: level2_cplx_impl.h:294

References GlobalFct::get_residuals(), m, and n.

◆ fill_in_contribution_to_jacobian() [2/2]

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

Calculate the contribution to the jacobian.

Fill in the contribution to the residuals

Reimplemented from oomph::GeneralisedElement.

220  {
221  //Call the generic routine with the flag set to 1
222  this->fill_in_generic_residual_contribution_interface(residuals,jacobian,1);
223  {
224  //Use finite differences to handle temperature variations
225  const unsigned n_node = this->nnode();
226  //Find the number of dofs in the element
227  const unsigned n_dof = this->ndof();
228  //Create newres vector
229  Vector<double> newres(n_dof);
230 
231  //Integer storage for local unknown
232  int local_unknown=0;
233 
234  //Use the default finite difference step
235  const double fd_step = this->Default_fd_jacobian_step;
236 
237  //Loop over the nodes
238  for(unsigned n=0;n<n_node;n++)
239  {
240  //Get the number of values stored at the node
241  unsigned t_index = this->T_index;
242 
243  //Get the local equation number
244  local_unknown = this->nodal_local_eqn(n,t_index);
245  //If it's not pinned
246  if(local_unknown >= 0)
247  {
248  //Store a pointer to the nodal data value
249  double *value_pt = this->node_pt(n)->value_pt(t_index);
250 
251  //Save the old value of the Nodal data
252  double old_var = *value_pt;
253 
254  //Increment the value of the Nodal data
255  *value_pt += fd_step;
256 
257  //Calculate the new residuals
258  this->get_residuals(newres);
259 
260  //Do finite differences
261  for(unsigned m=0;m<n_dof;m++)
262  {
263  double sum = (newres[m] - residuals[m])/fd_step;
264  //Stick the entry into the Jacobian matrix
265  jacobian(m,local_unknown) = sum;
266  }
267 
268  //Reset the Nodal data
269  *value_pt = old_var;
270  }
271  }
272 
273  //Use finite differences to handle the concentration variations
274  //Loop over the nodes again
275  for(unsigned n=0;n<n_node;n++)
276  {
277  //Get the number of values stored at the node
278  unsigned c_index = this->C_index;
279 
280  //Get the local equation number
281  local_unknown = this->nodal_local_eqn(n,c_index);
282  //If it's not pinned
283  if(local_unknown >= 0)
284  {
285  //Store a pointer to the nodal data value
286  double *value_pt = this->node_pt(n)->value_pt(c_index);
287 
288  //Save the old value of the Nodal data
289  double old_var = *value_pt;
290 
291  //Increment the value of the Nodal data
292  *value_pt += fd_step;
293 
294  //Calculate the new residuals
295  this->get_residuals(newres);
296 
297  //Do finite differences
298  for(unsigned m=0;m<n_dof;m++)
299  {
300  double sum = (newres[m] - residuals[m])/fd_step;
301  //Stick the entry into the Jacobian matrix
302  jacobian(m,local_unknown) = sum;
303  }
304 
305  //Reset the Nodal data
306  *value_pt = old_var;
307  }
308  }
309  }
310 
311  //Call the generic routine to handle the spine variables
312  SpineElement<FaceGeometry<ELEMENT> >::
313  fill_in_jacobian_from_geometric_data(jacobian);
314  }

References GlobalFct::get_residuals(), m, and n.

◆ fill_in_contribution_to_jacobian_and_mass_matrix() [1/2]

template<class ELEMENT >
void oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::fill_in_contribution_to_jacobian_and_mass_matrix ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian,
DenseMatrix< double > &  mass_matrix 
)
inlineprotectedvirtual

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

Reimplemented from oomph::GeneralisedElement.

480  {
481  //Add the contribution to the jacobian
482  this->fill_in_contribution_to_jacobian(residuals,jacobian);
483  //No mass matrix terms, but should probably do kinematic bit here
484  }
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::fill_in_contribution_to_jacobian
void fill_in_contribution_to_jacobian(Vector< double > &residuals, DenseMatrix< double > &jacobian)
Calculate the contribution to the jacobian.
Definition: surfactant.cc:209

References oomph::fill_in_contribution_to_jacobian().

◆ fill_in_contribution_to_jacobian_and_mass_matrix() [2/2]

template<class ELEMENT >
void oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::fill_in_contribution_to_jacobian_and_mass_matrix ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian,
DenseMatrix< double > &  mass_matrix 
)
inlineprotectedvirtual

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

Reimplemented from oomph::GeneralisedElement.

570  {
571  //Add the contribution to the jacobian
572  this->fill_in_contribution_to_jacobian(residuals,jacobian);
573  //No mass matrix terms, but should probably do kinematic bit here
574  }

References oomph::fill_in_contribution_to_jacobian().

◆ interpolated_C() [1/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::interpolated_C ( const Vector< double > &  s)
inlineprotected

Get the surfactant concentration.

127  {
128  //Find number of nodes
129  unsigned n_node = this->nnode();
130 
131  //Get the nodal index at which the unknown is stored
132  const unsigned c_index = C_index;
133 
134  //Local shape function
135  Shape psi(n_node);
136 
137  //Find values of shape function
138  this->shape(s,psi);
139 
140  //Initialise value of C
141  double C = 0.0;
142 
143  //Loop over the local nodes and sum
144  for(unsigned l=0;l<n_node;l++)
145  {
146  C += this->nodal_value(l,c_index)*psi(l);
147  }
148 
149  return(C);
150  }
oomph::FiniteElement::shape
virtual void shape(const Vector< double > &s, Shape &psi) const =0
oomph::Shape
Definition: shape.h:76

References oomph::OneDimLagrange::shape().

◆ interpolated_C() [2/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::interpolated_C ( const Vector< double > &  s)
inlineprotected

Get the surfactant concentration.

136  {
137  //Find number of nodes
138  unsigned n_node = this->nnode();
139 
140  //Get the nodal index at which the unknown is stored
141  const unsigned c_index = C_index;
142 
143  //Local shape function
144  Shape psi(n_node);
145 
146  //Find values of shape function
147  this->shape(s,psi);
148 
149  //Initialise value of C
150  double C = 0.0;
151 
152  //Loop over the local nodes and sum
153  for(unsigned l=0;l<n_node;l++)
154  {
155  C += this->nodal_value(l,c_index)*psi(l);
156  }
157 
158  return(C);
159  }

References oomph::OneDimLagrange::shape().

◆ interpolated_T() [1/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::interpolated_T ( const Vector< double > &  s)
inlineprotected

Get the temperature.

99  {
100  //Find number of nodes
101  unsigned n_node = this->nnode();
102 
103  //Get the nodal index at which the unknown is stored
104  const unsigned t_index = T_index;
105 
106  //Local shape function
107  Shape psi(n_node);
108 
109  //Find values of shape function
110  this->shape(s,psi);
111 
112  //Initialise value of t
113  double T = 0.0;
114 
115  //Loop over the local nodes and sum
116  for(unsigned l=0;l<n_node;l++)
117  {
118  T += this->nodal_value(l,t_index)*psi(l);
119  }
120 
121  return(T);
122  }

References oomph::OneDimLagrange::shape().

◆ interpolated_T() [2/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::interpolated_T ( const Vector< double > &  s)
inlineprotected

Get the temperature.

108  {
109  //Find number of nodes
110  unsigned n_node = this->nnode();
111 
112  //Get the nodal index at which the unknown is stored
113  const unsigned t_index = T_index;
114 
115  //Local shape function
116  Shape psi(n_node);
117 
118  //Find values of shape function
119  this->shape(s,psi);
120 
121  //Initialise value of t
122  double T = 0.0;
123 
124  //Loop over the local nodes and sum
125  for(unsigned l=0;l<n_node;l++)
126  {
127  T += this->nodal_value(l,t_index)*psi(l);
128  }
129 
130  return(T);
131  }

References oomph::OneDimLagrange::shape().

◆ k()

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::k ( )
inline

Return the reaction ratio.

638 {return *K_pt;}

◆ k_pt()

template<class ELEMENT >
double* & oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::k_pt ( )
inline

Access function for pointer to the reaction ratios.

659 {return K_pt;}

◆ ma() [1/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::ma ( )
inline

Return the Marangoni number.

528 {return *Ma_pt;}

◆ ma() [2/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::ma ( )
inline

Return the Marangoni number.

623 {return *Ma_pt;}

◆ ma_pt() [1/2]

template<class ELEMENT >
double* & oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::ma_pt ( )
inline

Access function for pointer to the Marangoni number.

540 {return Ma_pt;}

◆ ma_pt() [2/2]

template<class ELEMENT >
double* & oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::ma_pt ( )
inline

Access function for pointer to the Marangoni number.

641 {return Ma_pt;}

◆ output() [1/2]

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

Output the element data — pass (some measure of) the number of plot points per element

Reimplemented from oomph::FiniteElement.

555 {
556  //Set output Vector
557  Vector<double> s(1);
558 
559  //Tecplot header info
560  outfile << "ZONE I=" << n_plot << std::endl;
561 
562  //Loop over plot points
563  for(unsigned l=0;l<n_plot;l++)
564  {
565  s[0] = -1.0 + l*2.0/(n_plot-1);
566 
567  //Output the x,y,u,v
568  for(unsigned i=0;i<2;i++) outfile << this->interpolated_x(s,i) << " ";
569  for(unsigned i=0;i<2;i++) outfile << this->interpolated_u(s,i) << " ";
570  //Output a dummy pressure
571  outfile << 0.0 << " ";
572  //Output the temperature
573  outfile << interpolated_T(s) << " "
574  << interpolated_C(s) << std::endl;
575  }
576  outfile << std::endl;
577 }
oomph::FaceElement::interpolated_x
double interpolated_x(const Vector< double > &s, const unsigned &i) const
Definition: elements.h:4528
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::interpolated_T
double interpolated_T(const Vector< double > &s)
Get the temperature.
Definition: surfactant.cc:98
s
RealScalar s
Definition: level1_cplx_impl.h:130

References i, and s.

◆ output() [2/2]

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

Output the element data — pass (some measure of) the number of plot points per element

Reimplemented from oomph::FiniteElement.

663 {
664  //Set output Vector
665  Vector<double> s(1);
666 
667  //Tecplot header info
668  outfile << "ZONE I=" << n_plot << std::endl;
669 
670  //Loop over plot points
671  for(unsigned l=0;l<n_plot;l++)
672  {
673  s[0] = -1.0 + l*2.0/(n_plot-1);
674 
675  //Output the x,y,u,v
676  for(unsigned i=0;i<2;i++) outfile << this->interpolated_x(s,i) << " ";
677  for(unsigned i=0;i<2;i++) outfile << this->interpolated_u(s,i) << " ";
678  //Output a dummy pressure
679  outfile << 0.0 << " ";
680  //Output the temperature
681  outfile << interpolated_T(s) << " "
682  << interpolated_C(s) << std::endl;
683  }
684  outfile << std::endl;
685 }

References i, and s.

◆ peclet_s() [1/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::peclet_s ( )
inline

Return the surface peclect number.

534 {return *Peclet_S_pt;}

◆ peclet_s() [2/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::peclet_s ( )
inline

Return the surface peclect number.

629 {return *Peclet_S_pt;}

◆ peclet_s_pt() [1/2]

template<class ELEMENT >
double* & oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::peclet_s_pt ( )
inline

Access function for pointer to the surface Peclet number.

549 {return Peclet_S_pt;}

◆ peclet_s_pt() [2/2]

template<class ELEMENT >
double* & oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::peclet_s_pt ( )
inline

Access function for pointer to the surface Peclet number.

650 {return Peclet_S_pt;}

◆ peclet_strouhal_s() [1/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::peclet_strouhal_s ( )
inline

Return the surface peclect strouhal number.

537 {return *Peclet_Strouhal_S_pt;}

◆ peclet_strouhal_s() [2/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::peclet_strouhal_s ( )
inline

Return the surface peclect strouhal number.

632 {return *Peclet_Strouhal_S_pt;}

◆ peclet_strouhal_s_pt() [1/2]

template<class ELEMENT >
double* & oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::peclet_strouhal_s_pt ( )
inline

Access function for pointer to the surface Peclet x Strouhal number.

552 {return Peclet_Strouhal_S_pt;}

◆ peclet_strouhal_s_pt() [2/2]

template<class ELEMENT >
double* & oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::peclet_strouhal_s_pt ( )
inline

Access function for pointer to the surface Peclet x Strouhal number.

653 {return Peclet_Strouhal_S_pt;}

◆ sigma() [1/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::sigma ( const Vector< double > &  s)
inlineprotectedvirtual

The surface tension function is linear in the temperature with constant of proportionality equal to the Marangoni number.

Reimplemented from oomph::FluidInterfaceElement.

182  {
183  //Find the number of shape functions
184  const unsigned n_node = this->nnode();
185  //Now get the shape fuctions at the local coordinate
186  Shape psi(n_node);
187  this->shape(s,psi);
188 
189  //Now interpolate the temperature and surfactant concentration
190  double T = 0.0, C=0.0;
191  for(unsigned l=0;l<n_node;l++)
192  {
193  T += this->nodal_value(l,T_index)*psi(l);
194  C += this->nodal_value(l,C_index)*psi(l);
195  }
196 
197  //Get the Marangoni, Capillary and Elasticity numbers
198  double Ma = this->ma();
199  double Ca = this->ca();
200  double Beta = this->beta();
201  //Return the variable surface tension
202  //The additional multiplication by Ca will cancel with the 1/Ca
203  //in the underlying equations
204  return (1.0 - Ca*Ma*T - Ca*Beta*(C-1.0));
205  }
oomph::FluidInterfaceElement::ca
const double & ca() const
The value of the Capillary number.
Definition: interface_elements.h:473
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::beta
double beta()
Return the Elasticity number.
Definition: surfactant.cc:531
oomph::SpineLineMarangoniSurfactantFluidInterfaceElement::ma
double ma()
Return the Marangoni number.
Definition: surfactant.cc:528
Global_Physical_Variables::Beta
double Beta
Surface Elasticity number.
Definition: surfactant.cc:622
Global_Physical_Variables::Ca
double Ca
Capillary number.
Definition: fibre.cc:61
oomph::Global_Physical_Variables::Ma
double Ma
Surfactant Parameters.
Definition: refineable_two_layer_soluble_surfactant.cc:127

References beta, Global_Physical_Variables::Beta, Global_Physical_Variables::Ca, oomph::Global_Physical_Variables::Ma, and oomph::OneDimLagrange::shape().

◆ sigma() [2/2]

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::sigma ( const Vector< double > &  s)
inlineprotectedvirtual

The surface tension function is linear in the temperature with constant of proportionality equal to the Marangoni number.

Reimplemented from oomph::FluidInterfaceElement.

191  {
192  //Find the number of shape functions
193  const unsigned n_node = this->nnode();
194  //Now get the shape fuctions at the local coordinate
195  Shape psi(n_node);
196  this->shape(s,psi);
197 
198  //Now interpolate the temperature and surfactant concentration
199  double T = 0.0, C=0.0;
200  for(unsigned l=0;l<n_node;l++)
201  {
202  T += this->nodal_value(l,T_index)*psi(l);
203  C += this->nodal_value(l,C_index)*psi(l);
204  }
205 
206  //Get the Marangoni, Capillary and Elasticity numbers
207  double Ma = this->ma();
208  double Ca = this->ca();
209  double Beta = this->beta();
210  //Return the variable surface tension
211  //The additional multiplication by Ca will cancel with the 1/Ca
212  //in the underlying equations
213  return (1.0 - Ca*Ma*T - Ca*Beta*(C-1.0));
214  }
Global::Beta
double Beta
Definition: ff_step.cc:156

References beta, Global::Beta, Global_Physical_Variables::Ca, oomph::Global_Physical_Variables::Ma, and oomph::OneDimLagrange::shape().

Member Data Documentation

◆ Beta_pt

template<class ELEMENT >
double * oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::Beta_pt
private

Pointer to an Elasticity number.

◆ Bi_pt

template<class ELEMENT >
double * oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::Bi_pt
private

Pointer to a Biot number.

◆ C_bulk_index

template<class ELEMENT >
unsigned oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::C_bulk_index
private

Index at which the bulk concentration is stored at the nodes.

◆ C_index

template<class ELEMENT >
unsigned oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::C_index
private

Index at which the surfactant concentration is stored at the nodes

◆ D_pt

template<class ELEMENT >
double* oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::D_pt
private

Pointer to the diffusion ratios.

◆ Default_Physical_Constant_Value

template<class ELEMENT >
double oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::Default_Physical_Constant_Value = 1.0
staticprivate

Default value of the physical constants.

◆ K_pt

template<class ELEMENT >
double* oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::K_pt
private

Pointer to the reaction ratios.

◆ Ma_pt

template<class ELEMENT >
double * oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::Ma_pt
private

Pointer to a Marangoni number.

◆ Peclet_S_pt

template<class ELEMENT >
double * oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::Peclet_S_pt
private

Pointer to Surface Peclet number.

◆ Peclet_Strouhal_S_pt

template<class ELEMENT >
double * oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::Peclet_Strouhal_S_pt
private

Pointer to the surface Peclect Strouhal number.

◆ T_index

template<class ELEMENT >
unsigned oomph::SpineLineMarangoniSurfactantFluidInterfaceElement< ELEMENT >::T_index
private

Index at which the temperature is stored at the nodes.

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