oomph::AxisymmetricDerivatives Class Reference

#include <interface_elements.h>

Inheritance diagram for oomph::AxisymmetricDerivatives:

Public Member Functions
	AxisymmetricDerivatives ()

Protected Member Functions
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...

Detailed Description

Class that establishes the surface derivative functions for AxisymmetricInterfaceElements. These are defined in a separate class so that they can be used by other interface equation-type classes.

Constructor & Destructor Documentation

AxisymmetricDerivatives()

oomph::AxisymmetricDerivatives::AxisymmetricDerivatives ( )

inline

{}

Member Function Documentation

compute_surface_derivatives()

double oomph::AxisymmetricDerivatives::compute_surface_derivatives ( const Shape &  psi, const DShape &  dpsids, const DenseMatrix< double > &  interpolated_t, const Vector< double > &  interpolated_x, DShape &  dpsidS, DShape &  dpsidS_div  )

protected

Fill in the specific surface derivative calculations.

Specialise the surface derivatives for the axisymmetric interface case.

///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// /////////////////////////////////////////////////////////////////////

  {
    // Initially the same as the 2D case
    const unsigned n_shape = psi.nindex1();
    const unsigned n_dim = 2;
 
    // Calculate the only entry of the surface
    // metric tensor (square length of the tangent vector)
    double a11 = interpolated_t(0, 0) * interpolated_t(0, 0) +
                 interpolated_t(0, 1) * interpolated_t(0, 1);
 
    // Now set the derivative terms of the shape functions
    for (unsigned l = 0; l < n_shape; l++)
    {
      for (unsigned i = 0; i < n_dim; i++)
      {
        dpsidS(l, i) = dpsids(l, 0) * interpolated_t(0, i) / a11;
        // Set the standard components of the divergence
        dpsidS_div(l, i) = dpsidS(l, i);
      }
    }
 
    const double r = interpolated_x[0];
 
    // The surface divergence is different because we need
    // to take account of variation of the base vectors
    for (unsigned l = 0; l < n_shape; l++)
    {
      dpsidS_div(l, 0) += psi(l) / r;
    }
 
    // Return the jacobian, needs to be multiplied by r
    return r * sqrt(a11);
  }

References i, oomph::Shape::nindex1(), UniformPSDSelfTest::r, and sqrt().

---

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

• interface_elements.h
• interface_elements.cc