oomph::RectangleWithHoleDomain Class Reference

Rectangular domain with circular whole. More...

#include <rectangle_with_hole_domain.h>

Inheritance diagram for oomph::RectangleWithHoleDomain:

Public Member Functions
	RectangleWithHoleDomain (GeomObject *cylinder_pt, const double &length)
	~RectangleWithHoleDomain ()
	Destructor: Empty; cleanup done in base class. More...
void	linear_interpolate (Vector< double > left, Vector< double > right, const double &s, Vector< double > &f)
void	macro_element_boundary (const unsigned &time, const unsigned &m, const unsigned &direction, const Vector< double > &s, Vector< double > &f)
Public Member Functions inherited from oomph::Domain
	Domain ()
	Constructor. More...
	Domain (const Domain &)=delete
	Broken copy constructor. More...
void	operator= (const Domain &)=delete
	Broken assignment operator. More...
virtual	~Domain ()
MacroElement *	macro_element_pt (const unsigned &i)
	Access to i-th macro element. More...
unsigned	nmacro_element ()
	Number of macro elements in domain. More...
void	output (const std::string &filename, const unsigned &nplot)
	Output macro elements. More...
void	output (std::ostream &outfile, const unsigned &nplot)
	Output macro elements. More...
virtual void	macro_element_boundary (const double &t, const unsigned &i_macro, const unsigned &i_direct, const Vector< double > &s, Vector< double > &f)
void	macro_element_boundary (const unsigned &i_macro, const unsigned &i_direct, const Vector< double > &s, Vector< double > &f)
void	output_macro_element_boundaries (const std::string &filename, const unsigned &nplot)
	Output all macro element boundaries as tecplot zones. More...
void	output_macro_element_boundaries (std::ostream &outfile, const unsigned &nplot)
	Output all macro element boundaries as tecplot zones. More...
virtual void	dmacro_element_boundary (const unsigned &t, const unsigned &i_macro, const unsigned &i_direct, const Vector< double > &s, Vector< double > &f)
virtual void	dmacro_element_boundary (const double &t, const unsigned &i_macro, const unsigned &i_direct, const Vector< double > &s, Vector< double > &f)
void	dmacro_element_boundary (const unsigned &i_macro, const unsigned &i_direct, const Vector< double > &s, Vector< double > &f)
virtual void	d2macro_element_boundary (const unsigned &t, const unsigned &i_macro, const unsigned &i_direct, const Vector< double > &s, Vector< double > &f)
virtual void	d2macro_element_boundary (const double &t, const unsigned &i_macro, const unsigned &i_direct, const Vector< double > &s, Vector< double > &f)
void	d2macro_element_boundary (const unsigned &i_macro, const unsigned &i_direct, const Vector< double > &s, Vector< double > &f)

Private Attributes
Vector< double >	Lower_left
	Lower left corner of rectangle. More...
Vector< double >	Lower_right
	Lower right corner of rectangle. More...
Vector< double >	Lower_mid_left
	Where the "radial" line from circle meets lower boundary on left. More...
Vector< double >	Lower_mid_right
	Where the "radial" line from circle meets lower boundary on right. More...
Vector< double >	Upper_left
	Upper left corner of rectangle. More...
Vector< double >	Upper_right
	Upper right corner of rectangle. More...
Vector< double >	Upper_mid_left
	Where the "radial" line from circle meets upper boundary on left. More...
Vector< double >	Upper_mid_right
	Where the "radial" line from circle meets upper boundary on right. More...
GeomObject *	Cylinder_pt
	Pointer to geometric object that represents the central cylinder. More...

Additional Inherited Members
Protected Attributes inherited from oomph::Domain
Vector< MacroElement * >	Macro_element_pt
	Vector of pointers to macro elements. More...

Detailed Description

Rectangular domain with circular whole.

Constructor & Destructor Documentation

RectangleWithHoleDomain()

oomph::RectangleWithHoleDomain::RectangleWithHoleDomain ( GeomObject *  cylinder_pt, const double &  length  )

inline

Constructor. Pass pointer to geometric object that represents the cylinder, the length of the (square) domain. The GeomObject must be parametrised such that \(\zeta \in [0,2\pi]\) sweeps around the circumference in anticlockwise direction.

      : Cylinder_pt(cylinder_pt)
    {
      // Vertices of rectangle
      Lower_left.resize(2);
      Lower_left[0] = -0.5 * length;
      Lower_left[1] = -0.5 * length;
 
      Upper_left.resize(2);
      Upper_left[0] = -0.5 * length;
      Upper_left[1] = 0.5 * length;
 
      Lower_right.resize(2);
      Lower_right[0] = 0.5 * length;
      Lower_right[1] = -0.5 * length;
 
      Upper_right.resize(2);
      Upper_right[0] = 0.5 * length;
      Upper_right[1] = 0.5 * length;
 
 
      // Coordinates of points where the "radial" lines from central
      // cylinder meet the upper and lower boundaries
      Lower_mid_left.resize(2);
      Lower_mid_left[0] = -0.5 * length;
      Lower_mid_left[1] = -0.5 * length;
 
      Upper_mid_left.resize(2);
      Upper_mid_left[0] = -0.5 * length;
      Upper_mid_left[1] = 0.5 * length;
 
      Lower_mid_right.resize(2);
      Lower_mid_right[0] = 0.5 * length;
      Lower_mid_right[1] = -0.5 * length;
 
      Upper_mid_right.resize(2);
      Upper_mid_right[0] = 0.5 * length;
      Upper_mid_right[1] = 0.5 * length;
 
 
      // There are four macro elements
      Macro_element_pt.resize(4);
 
      // Build the 2D macro elements
      for (unsigned i = 0; i < 4; i++)
      {
        Macro_element_pt[i] = new QMacroElement<2>(this, i);
      }
    }

References i, Lower_left, Lower_mid_left, Lower_mid_right, Lower_right, oomph::Domain::Macro_element_pt, Upper_left, Upper_mid_left, Upper_mid_right, and Upper_right.

~RectangleWithHoleDomain()

oomph::RectangleWithHoleDomain::~RectangleWithHoleDomain ( )

inline

Destructor: Empty; cleanup done in base class.

{}

Member Function Documentation

linear_interpolate()

void oomph::RectangleWithHoleDomain::linear_interpolate ( Vector< double >  left, Vector< double >  right, const double &  s, Vector< double > &  f  )

inline

Helper function to interpolate linearly between the "right" and "left" points; \( s \in [-1,1] \)

    {
      for (unsigned i = 0; i < 2; i++)
      {
        f[i] = left[i] + (right[i] - left[i]) * 0.5 * (s + 1.0);
      }
    }

References f(), i, and s.

Referenced by macro_element_boundary().

macro_element_boundary()

void oomph::RectangleWithHoleDomain::macro_element_boundary ( const unsigned &  time, const unsigned &  m, const unsigned &  direction, const Vector< double > &  s, Vector< double > &  f  )

inlinevirtual

Parametrisation of macro element boundaries: f(s) is the position vector to macro-element m's boundary in the specified direction [N/S/E/W] at the specfied discrete time level (time=0: present; time>0: previous)

Implements oomph::Domain.

    {
#ifdef WARN_ABOUT_SUBTLY_CHANGED_OOMPH_INTERFACES
      // Warn about time argument being moved to the front
      OomphLibWarning(
        "Order of function arguments has changed between versions 0.8 and 0.85",
        "RectangleWithHoleDomain::macro_element_boundary(...)",
        OOMPH_EXCEPTION_LOCATION);
#endif
 
      // Lagrangian coordinate along surface of cylinder
      Vector<double> xi(1);
 
      // Point on circle
      Vector<double> point_on_circle(2);
 
      using namespace QuadTreeNames;
 
      // Switch on the macro element
      switch (m)
      {
          // Macro element 0, is is immediately left of the cylinder
        case 0:
 
          switch (direction)
          {
            case N:
              xi[0] = 3.0 * atan(1.0);
              Cylinder_pt->position(time, xi, point_on_circle);
              linear_interpolate(Upper_mid_left, point_on_circle, s[0], f);
              break;
 
            case S:
              xi[0] = -3.0 * atan(1.0);
              Cylinder_pt->position(time, xi, point_on_circle);
              linear_interpolate(Lower_mid_left, point_on_circle, s[0], f);
              break;
 
            case W:
              linear_interpolate(Lower_mid_left, Upper_mid_left, s[0], f);
              break;
 
            case E:
              xi[0] = 5.0 * atan(1.0) - 2.0 * atan(1.0) * 0.5 * (1.0 + s[0]);
              Cylinder_pt->position(time, xi, f);
              break;
 
            default:
 
              std::ostringstream error_stream;
              error_stream << "Direction is incorrect:  " << direction
                           << std::endl;
              throw OomphLibError(error_stream.str(),
                                  OOMPH_CURRENT_FUNCTION,
                                  OOMPH_EXCEPTION_LOCATION);
          }
 
          break;
 
        // Macro element 1, is immediately above the cylinder
        case 1:
 
          switch (direction)
          {
            case N:
              linear_interpolate(Upper_mid_left, Upper_mid_right, s[0], f);
              break;
 
            case S:
              xi[0] = 3.0 * atan(1.0) - 2.0 * atan(1.0) * 0.5 * (1.0 + s[0]);
              Cylinder_pt->position(time, xi, f);
              break;
 
            case W:
              xi[0] = 3.0 * atan(1.0);
              Cylinder_pt->position(time, xi, point_on_circle);
              linear_interpolate(point_on_circle, Upper_mid_left, s[0], f);
              break;
 
            case E:
              xi[0] = 1.0 * atan(1.0);
              Cylinder_pt->position(time, xi, point_on_circle);
              linear_interpolate(point_on_circle, Upper_mid_right, s[0], f);
              break;
 
            default:
 
              std::ostringstream error_stream;
              error_stream << "Direction is incorrect:  " << direction
                           << std::endl;
              throw OomphLibError(error_stream.str(),
                                  OOMPH_CURRENT_FUNCTION,
                                  OOMPH_EXCEPTION_LOCATION);
          }
 
          break;
 
          // Macro element 2, is immediately right of the cylinder
        case 2:
 
          switch (direction)
          {
            case N:
              xi[0] = 1.0 * atan(1.0);
              Cylinder_pt->position(time, xi, point_on_circle);
              linear_interpolate(point_on_circle, Upper_mid_right, s[0], f);
              break;
 
            case S:
              xi[0] = -1.0 * atan(1.0);
              Cylinder_pt->position(time, xi, point_on_circle);
              linear_interpolate(point_on_circle, Lower_mid_right, s[0], f);
              break;
 
            case W:
              xi[0] = -atan(1.0) + 2.0 * atan(1.0) * 0.5 * (1.0 + s[0]);
              Cylinder_pt->position(time, xi, f);
              break;
 
            case E:
              linear_interpolate(Lower_mid_right, Upper_mid_right, s[0], f);
              break;
 
            default:
 
              std::ostringstream error_stream;
              error_stream << "Direction is incorrect:  " << direction
                           << std::endl;
              throw OomphLibError(error_stream.str(),
                                  OOMPH_CURRENT_FUNCTION,
                                  OOMPH_EXCEPTION_LOCATION);
          }
 
          break;
 
          // Macro element 3, is immediately below cylinder
        case 3:
 
          switch (direction)
          {
            case N:
              xi[0] = -3.0 * atan(1.0) + 2.0 * atan(1.0) * 0.5 * (1.0 + s[0]);
              Cylinder_pt->position(time, xi, f);
              break;
 
            case S:
              linear_interpolate(Lower_mid_left, Lower_mid_right, s[0], f);
              break;
 
            case W:
              xi[0] = -3.0 * atan(1.0);
              Cylinder_pt->position(time, xi, point_on_circle);
              linear_interpolate(Lower_mid_left, point_on_circle, s[0], f);
              break;
 
            case E:
              xi[0] = -1.0 * atan(1.0);
              Cylinder_pt->position(time, xi, point_on_circle);
              linear_interpolate(Lower_mid_right, point_on_circle, s[0], f);
              break;
 
            default:
 
              std::ostringstream error_stream;
              error_stream << "Direction is incorrect:  " << direction
                           << std::endl;
              throw OomphLibError(error_stream.str(),
                                  OOMPH_CURRENT_FUNCTION,
                                  OOMPH_EXCEPTION_LOCATION);
          }
 
          break;
 
        default:
 
          std::ostringstream error_stream;
          error_stream << "Wrong macro element number" << m << std::endl;
          throw OomphLibError(error_stream.str(),
                              OOMPH_CURRENT_FUNCTION,
                              OOMPH_EXCEPTION_LOCATION);
      }
    }

References Eigen::bfloat16_impl::atan(), Cylinder_pt, Global_Physical_Variables::E, f(), linear_interpolate(), Lower_mid_left, Lower_mid_right, m, N, OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, oomph::GeomObject::position(), s, oomph::QuadTreeNames::S, Upper_mid_left, Upper_mid_right, and oomph::QuadTreeNames::W.

Member Data Documentation

Cylinder_pt

GeomObject * RectangleWithHoleDomain::Cylinder_pt

private

Pointer to geometric object that represents the central cylinder.

Referenced by macro_element_boundary().

Lower_left

Vector< double > RectangleWithHoleDomain::Lower_left

private

Lower left corner of rectangle.

Referenced by RectangleWithHoleDomain().

Lower_mid_left

Vector< double > RectangleWithHoleDomain::Lower_mid_left

private

Where the "radial" line from circle meets lower boundary on left.

Referenced by macro_element_boundary(), and RectangleWithHoleDomain().

Lower_mid_right

Vector< double > RectangleWithHoleDomain::Lower_mid_right

private

Where the "radial" line from circle meets lower boundary on right.

Referenced by macro_element_boundary(), and RectangleWithHoleDomain().

Lower_right

Vector< double > RectangleWithHoleDomain::Lower_right

private

Lower right corner of rectangle.

Referenced by RectangleWithHoleDomain().

Upper_left

Vector< double > RectangleWithHoleDomain::Upper_left

private

Upper left corner of rectangle.

Referenced by RectangleWithHoleDomain().

Upper_mid_left

Vector< double > RectangleWithHoleDomain::Upper_mid_left

private

Where the "radial" line from circle meets upper boundary on left.

Referenced by macro_element_boundary(), and RectangleWithHoleDomain().

Upper_mid_right

Vector< double > RectangleWithHoleDomain::Upper_mid_right

private

Where the "radial" line from circle meets upper boundary on right.

Referenced by macro_element_boundary(), and RectangleWithHoleDomain().

Upper_right

Vector< double > RectangleWithHoleDomain::Upper_right

private

Upper right corner of rectangle.

Referenced by RectangleWithHoleDomain().

---

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

• rectangle_with_hole_domain.h
• adaptive_hopf_with_separate_meshes.cc
• flow_past_cylinder.cc