Rectangular domain with circular whole DRAIG: This looks like a redefinition of the RectangleWithHoleAndAnnularRegionDomain in src/meshes but it creates 8 macro-elements instead of 4 macro-elements and creates an annular region around the cylinder. It's probably a good idea to rename this class to avoid ambiguity and a name clash... More...

#include <rectangle_with_moving_cylinder_mesh.h>

Inheritance diagram for oomph::RectangleWithHoleAndAnnularRegionDomain:

Public Member Functions
	RectangleWithHoleAndAnnularRegionDomain (GeomObject *cylinder_pt, const double &annular_region_radius, const double &length)
	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.

	~RectangleWithHoleAndAnnularRegionDomain ()
	Destructor: Empty; macro elements are deleted in base class destructor.

void	project_point_on_cylinder_to_annular_boundary (const unsigned &time, const Vector< double > &xi, Vector< double > &r)
	Helper function that, given the Lagrangian coordinate, xi, (associated with a point on the cylinder), returns the corresponding point on the outer boundary of the annular region (where the inner boundary is prescribed by the boundary of the cylinder)

void	project_point_on_cylinder_to_annular_boundary (const double &time, const Vector< double > &xi, Vector< double > &r)
	Helper function that, given the Lagrangian coordinate, xi, (associated with a point on the cylinder), returns the corresponding point on the outer boundary of the annular region (where the inner boundary is prescribed by the boundary of the cylinder)

void	linear_interpolate (const Vector< double > &left, const Vector< double > &right, const double &s, Vector< double > &f)
	Helper function to interpolate linearly between the "right" and "left" points; $s \in [-1,1]$ .

void	macro_element_boundary (const double &time, const unsigned &m, const unsigned &direction, const Vector< double > &s, Vector< double > &f)
	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 specified discrete time level (time=0: present; time>0: previous)

void	macro_element_boundary (const unsigned &time, const unsigned &m, const unsigned &direction, const Vector< double > &s, Vector< double > &f)
	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 specified discrete time level (time=0: present; time>0: previous)

Private Attributes
Vector< double >	Lower_left
	Lower left corner of rectangle.

Vector< double >	Lower_right
	Lower right corner of rectangle.

Vector< double >	Lower_mid_left
	Where the "radial" line from circle meets lower boundary on left.

Vector< double >	Lower_mid_right
	Where the "radial" line from circle meets lower boundary on right.

Vector< double >	Upper_left
	Upper left corner of rectangle.

Vector< double >	Upper_right
	Upper right corner of rectangle.

Vector< double >	Upper_mid_left
	Where the "radial" line from circle meets upper boundary on left.

Vector< double >	Upper_mid_right
	Where the "radial" line from circle meets upper boundary on right.

GeomObject *	Cylinder_pt
	Pointer to geometric object that represents the central cylinder.

double	Annular_region_radius
	The radius of the outer boundary of the annular region whose inner boundary is described by Cylinder_pt and outer boundary is a circle with radius lying between half the length of the bounding box and the radius cylinder.

Detailed Description

Rectangular domain with circular whole DRAIG: This looks like a redefinition of the RectangleWithHoleAndAnnularRegionDomain in src/meshes but it creates 8 macro-elements instead of 4 macro-elements and creates an annular region around the cylinder. It's probably a good idea to rename this class to avoid ambiguity and a name clash...

Definition at line 53 of file rectangle_with_moving_cylinder_mesh.h.

Constructor & Destructor Documentation

◆ RectangleWithHoleAndAnnularRegionDomain()

oomph::RectangleWithHoleAndAnnularRegionDomain::RectangleWithHoleAndAnnularRegionDomain	(	GeomObject *	cylinder_pt,
		const double &	annular_region_radius,
		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.

Definition at line 61 of file rectangle_with_moving_cylinder_mesh.h.

References Lower_left, Lower_mid_left, Lower_mid_right, Lower_right, Upper_left, Upper_mid_left, Upper_mid_right, and Upper_right.

◆ ~RectangleWithHoleAndAnnularRegionDomain()

oomph::RectangleWithHoleAndAnnularRegionDomain::~RectangleWithHoleAndAnnularRegionDomain ( )

inline

Destructor: Empty; macro elements are deleted in base class destructor.

Definition at line 116 of file rectangle_with_moving_cylinder_mesh.h.

Member Function Documentation

◆ linear_interpolate()

void oomph::RectangleWithHoleAndAnnularRegionDomain::linear_interpolate	(	const Vector< double > &	left,
		const 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]$ .

Definition at line 138 of file rectangle_with_moving_cylinder_mesh.h.

Referenced by macro_element_boundary(), and macro_element_boundary().

◆ macro_element_boundary() [1/2]

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

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 specified discrete time level (time=0: present; time>0: previous)

Definition at line 131 of file rectangle_with_moving_cylinder_mesh.template.cc.

References Cylinder_pt, linear_interpolate(), Lower_mid_left, Lower_mid_right, project_point_on_cylinder_to_annular_boundary(), Upper_mid_left, and Upper_mid_right.

◆ macro_element_boundary() [2/2]

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

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 specified discrete time level (time=0: present; time>0: previous)

Definition at line 521 of file rectangle_with_moving_cylinder_mesh.template.cc.

References Cylinder_pt, linear_interpolate(), Lower_mid_left, Lower_mid_right, project_point_on_cylinder_to_annular_boundary(), Upper_mid_left, and Upper_mid_right.

◆ project_point_on_cylinder_to_annular_boundary() [1/2]

void oomph::RectangleWithHoleAndAnnularRegionDomain::project_point_on_cylinder_to_annular_boundary	(	const double &	time,
		const Vector< double > &	xi,
		Vector< double > &	r
	)

Helper function that, given the Lagrangian coordinate, xi, (associated with a point on the cylinder), returns the corresponding point on the outer boundary of the annular region (where the inner boundary is prescribed by the boundary of the cylinder)

Definition at line 87 of file rectangle_with_moving_cylinder_mesh.template.cc.

References Annular_region_radius, and Cylinder_pt.

◆ project_point_on_cylinder_to_annular_boundary() [2/2]

void oomph::RectangleWithHoleAndAnnularRegionDomain::project_point_on_cylinder_to_annular_boundary	(	const unsigned &	time,
		const Vector< double > &	xi,
		Vector< double > &	r
	)

Helper function that, given the Lagrangian coordinate, xi, (associated with a point on the cylinder), returns the corresponding point on the outer boundary of the annular region (where the inner boundary is prescribed by the boundary of the cylinder)

Rectangular domain with circular whole.

Definition at line 42 of file rectangle_with_moving_cylinder_mesh.template.cc.

References Annular_region_radius, and Cylinder_pt.

Referenced by macro_element_boundary(), and macro_element_boundary().

Member Data Documentation

◆ Annular_region_radius

double oomph::RectangleWithHoleAndAnnularRegionDomain::Annular_region_radius

private

The radius of the outer boundary of the annular region whose inner boundary is described by Cylinder_pt and outer boundary is a circle with radius lying between half the length of the bounding box and the radius cylinder.

Definition at line 204 of file rectangle_with_moving_cylinder_mesh.h.

Referenced by project_point_on_cylinder_to_annular_boundary(), and project_point_on_cylinder_to_annular_boundary().