#include <geom_obj_with_boundary.h>

Inheritance diagram for oomph::WarpedCircularDisk:

Public Member Functions
	WarpedCircularDisk (const double &epsilon, const unsigned &n, const double &z_offset=0.0)

	WarpedCircularDisk ()
	Empty default constructor. More...

	WarpedCircularDisk (const WarpedCircularDisk &dummy)=delete
	Broken copy constructor. More...

void	operator= (const WarpedCircularDisk &)=delete
	Broken assignment operator. More...

virtual	~WarpedCircularDisk ()
	Destructor. More...

double &	epsilon ()
	Access fct to amplitude of disk warping. More...

void	position (const Vector< double > &zeta, Vector< double > &r) const
	Position Vector at Lagrangian coordinate zeta. More...

void	position (const unsigned &t, const Vector< double > &zeta, Vector< double > &r) const

void	boundary_triad (const unsigned &b, const double &zeta_bound, Vector< double > &r, Vector< double > &tangent, Vector< double > &normal, Vector< double > &binormal)
	Boundary triad on boundary b at boundary coordinate zeta_bound. More...

Public Member Functions inherited from oomph::DiskLikeGeomObjectWithBoundaries
	DiskLikeGeomObjectWithBoundaries ()
	Constructor. More...

unsigned	nboundary () const
	How many boundaries do we have? More...

void	position_on_boundary (const unsigned &b, const double &zeta_bound, Vector< double > &r) const

void	zeta_on_boundary (const unsigned &b, const double &zeta_bound, Vector< double > &zeta) const

GeomObject *	boundary_parametrising_geom_object_pt (const unsigned &b) const

double	zeta_boundary_start (const unsigned &b) const

double	zeta_boundary_end (const unsigned &b) const

void	output_boundaries (const unsigned &nplot, std::ofstream &two_d_boundaries_file, std::ofstream &three_d_boundaries_file)

void	output_boundaries_and_triads (const unsigned &nplot, std::ofstream &two_d_boundaries_file, std::ofstream &three_d_boundaries_file, std::ofstream &boundaries_tangent_file, std::ofstream &boundaries_normal_file, std::ofstream &boundaries_binormal_file)

void	add_region_coordinates (const unsigned &r, Vector< double > &zeta_in_region)

std::map< unsigned, Vector< double > >	zeta_in_region () const
	Return map that stores zeta coordinates of points that identify regions. More...

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 unsigned	ngeom_data () const

virtual Data *	geom_data_pt (const unsigned &j)

virtual void	position (const double &t, const Vector< double > &zeta, Vector< double > &r) const

virtual void	dposition_dt (const Vector< double > &zeta, const unsigned &j, Vector< double > &drdt)

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

virtual void	locate_zeta (const Vector< double > &zeta, GeomObject *&sub_geom_object_pt, Vector< double > &s, const bool &use_coordinate_as_initial_guess=false)

virtual void	interpolated_zeta (const Vector< double > &s, Vector< double > &zeta) const

Private Member Functions
double	w (const double &r, const double &phi) const
	Vertical deflection. More...

double	dwdr (const double &r, const double &phi) const
	Deriv of vertical deflection w.r.t. radius. More...

double	dwdphi (const double &r, const double &phi) const
	Deriv of vertical deflection w.r.t. angle. More...

Private Attributes
double	Epsilon
	Amplitude of non-axisymmetric deformation. More...

unsigned	N
	Wavenumber of non-axisymmetric deformation. More...

double	Z_offset
	Vertical offset. More...

Additional Inherited Members
Protected Attributes inherited from oomph::DiskLikeGeomObjectWithBoundaries
Vector< double >	Zeta_boundary_start

Vector< double >	Zeta_boundary_end

Vector< GeomObject * >	Boundary_parametrising_geom_object_pt

std::map< unsigned, Vector< double > >	Zeta_in_region
	Map to store zeta coordinates of points that identify regions. 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

Detailed Description

////////////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////////// Warped disk in 3d: zeta[0]=x; zeta[1]=y (so it doesn't have coordinate singularities), with specification of two boundaries (b=0,1) that turn the whole thing into a circular disk.

Constructor & Destructor Documentation

◆ WarpedCircularDisk() [1/3]

oomph::WarpedCircularDisk::WarpedCircularDisk	(	const double &	epsilon,
		const unsigned &	n,
		const double &	z_offset = `0.0`
	)

inline

Constructor. Pass amplitude and azimuthal wavenumber of warping as arguments. Can specify vertical offset as final, optional argument.

       : Epsilon(epsilon), N(n), Z_offset(z_offset)
     {
       // How many boundaries do we have?
       unsigned nb = 2;
       Boundary_parametrising_geom_object_pt.resize(nb);
       Zeta_boundary_start.resize(nb, 0.0);
       Zeta_boundary_end.resize(nb, 0.0);
  
       // GeomObject<1,2> representing the first boundary
       Boundary_parametrising_geom_object_pt[0] = new Ellipse(1.0, 1.0);
       Zeta_boundary_start[0] = 0.0;
       Zeta_boundary_end[0] = MathematicalConstants::Pi;
  
       // GeomObject<1,2> representing the second boundary
       Boundary_parametrising_geom_object_pt[1] = new Ellipse(1.0, 1.0);
       Zeta_boundary_start[1] = MathematicalConstants::Pi;
       Zeta_boundary_end[1] = 2.0 * MathematicalConstants::Pi;
     }

References oomph::DiskLikeGeomObjectWithBoundaries::Boundary_parametrising_geom_object_pt, nb, oomph::MathematicalConstants::Pi, oomph::DiskLikeGeomObjectWithBoundaries::Zeta_boundary_end, and oomph::DiskLikeGeomObjectWithBoundaries::Zeta_boundary_start.

◆ WarpedCircularDisk() [2/3]

oomph::WarpedCircularDisk::WarpedCircularDisk ( )

inline

Empty default constructor.

     {
       throw OomphLibError("Don't call default constructor!",
                           OOMPH_CURRENT_FUNCTION,
                           OOMPH_EXCEPTION_LOCATION);
     }

References OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.

◆ WarpedCircularDisk() [3/3]

oomph::WarpedCircularDisk::WarpedCircularDisk ( const WarpedCircularDisk & dummy )

delete

Broken copy constructor.

◆ ~WarpedCircularDisk()

virtual oomph::WarpedCircularDisk::~WarpedCircularDisk ( )

inlinevirtual

Destructor.

     {
       unsigned n = nboundary();
       for (unsigned b = 0; b < n; b++)
       {
         delete Boundary_parametrising_geom_object_pt[b];
         Boundary_parametrising_geom_object_pt[b] = 0;
       }
     }

References b, oomph::DiskLikeGeomObjectWithBoundaries::Boundary_parametrising_geom_object_pt, n, and oomph::DiskLikeGeomObjectWithBoundaries::nboundary().

Member Function Documentation

◆ boundary_triad()

void oomph::WarpedCircularDisk::boundary_triad	(	const unsigned &	b,
		const double &	zeta_bound,
		Vector< double > &	r,
		Vector< double > &	tangent,
		Vector< double > &	normal,
		Vector< double > &	binormal
	)

inlinevirtual

Boundary triad on boundary b at boundary coordinate zeta_bound.

Reimplemented from oomph::DiskLikeGeomObjectWithBoundaries.

     {
       double phi = zeta_bound;
  
       // Position Vector
       r[0] = cos(phi);
       r[1] = sin(phi);
       r[2] = Z_offset + w(1.0, phi);
  
       Vector<double> dr_dr(3);
       dr_dr[0] = cos(phi);
       dr_dr[1] = sin(phi);
       dr_dr[2] = dwdr(1.0, phi);
       double inv_norm = 1.0 / sqrt(dr_dr[0] * dr_dr[0] + dr_dr[1] * dr_dr[1] +
                                    dr_dr[2] * dr_dr[2]);
  
       normal[0] = dr_dr[0] * inv_norm;
       normal[1] = dr_dr[1] * inv_norm;
       normal[2] = dr_dr[2] * inv_norm;
  
       Vector<double> dr_dphi(3);
       dr_dphi[0] = -sin(phi);
       dr_dphi[1] = cos(phi);
       dr_dphi[2] = dwdphi(1.0, phi);
  
       inv_norm = 1.0 / sqrt(dr_dphi[0] * dr_dphi[0] + dr_dphi[1] * dr_dphi[1] +
                             dr_dphi[2] * dr_dphi[2]);
  
       tangent[0] = dr_dphi[0] * inv_norm;
       tangent[1] = dr_dphi[1] * inv_norm;
       tangent[2] = dr_dphi[2] * inv_norm;
  
       binormal[0] = tangent[1] * normal[2] - tangent[2] * normal[1];
       binormal[1] = tangent[2] * normal[0] - tangent[0] * normal[2];
       binormal[2] = tangent[0] * normal[1] - tangent[1] * normal[0];
     }

References cos(), dwdphi(), dwdr(), WallFunction::normal(), UniformPSDSelfTest::r, sin(), sqrt(), w(), and Z_offset.

◆ dwdphi()

double oomph::WarpedCircularDisk::dwdphi	(	const double &	r,
		const double &	phi
	)		const

inlineprivate

Deriv of vertical deflection w.r.t. angle.

     {
       return -Epsilon * double(N) * sin(double(N) * phi) * pow(r, 2);
     }

References Epsilon, N, Eigen::bfloat16_impl::pow(), UniformPSDSelfTest::r, and sin().

Referenced by boundary_triad().

◆ dwdr()

double oomph::WarpedCircularDisk::dwdr	(	const double &	r,
		const double &	phi
	)		const

inlineprivate

Deriv of vertical deflection w.r.t. radius.

     {
       return Epsilon * cos(double(N) * phi) * 2.0 * r;
     }

References cos(), Epsilon, N, and UniformPSDSelfTest::r.

Referenced by boundary_triad().

◆ epsilon()

double& oomph::WarpedCircularDisk::epsilon ( )

inline

Access fct to amplitude of disk warping.

     {
       return Epsilon;
     }

References Epsilon.

◆ operator=()

void oomph::WarpedCircularDisk::operator= ( const WarpedCircularDisk & )

delete

Broken assignment operator.

◆ position() [1/2]

void oomph::WarpedCircularDisk::position	(	const unsigned &	t,
		const Vector< double > &	zeta,
		Vector< double > &	r
	)		const

inlinevirtual

Parametrised position on object: r(zeta). Evaluated at previous timestep. t=0: current time; t>0: previous timestep. Object is steady so calls time-independent version

Reimplemented from oomph::GeomObject.

     {
       position(zeta, r);
     }

References position(), UniformPSDSelfTest::r, and Eigen::zeta().

◆ position() [2/2]

void oomph::WarpedCircularDisk::position	(	const Vector< double > &	zeta,
		Vector< double > &	r
	)		const

inlinevirtual

Position Vector at Lagrangian coordinate zeta.

Implements oomph::GeomObject.

     {
       // Position Vector
       r[0] = zeta[0];
       r[1] = zeta[1];
       double radius = sqrt(r[0] * r[0] + r[1] * r[1]);
       double phi = atan2(r[1], r[0]);
       r[2] = Z_offset + w(radius, phi);
     }

References atan2(), UniformPSDSelfTest::r, UniformPSDSelfTest::radius, sqrt(), w(), Z_offset, and Eigen::zeta().

Referenced by main(), and position().

◆ w()

double oomph::WarpedCircularDisk::w	(	const double &	r,
		const double &	phi
	)		const

inlineprivate

Vertical deflection.

     {
       return Epsilon * cos(double(N) * phi) * pow(r, 2);
     }

References cos(), Epsilon, N, Eigen::bfloat16_impl::pow(), and UniformPSDSelfTest::r.

Referenced by boundary_triad(), and position().

Member Data Documentation

◆ Epsilon

double oomph::WarpedCircularDisk::Epsilon

private

Amplitude of non-axisymmetric deformation.

Referenced by dwdphi(), dwdr(), epsilon(), and w().

◆ N

unsigned oomph::WarpedCircularDisk::N

private

Wavenumber of non-axisymmetric deformation.

Referenced by dwdphi(), dwdr(), and w().

◆ Z_offset

double oomph::WarpedCircularDisk::Z_offset

private

Vertical offset.

Referenced by boundary_triad(), and position().

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

geom_obj_with_boundary.h

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ WarpedCircularDisk() [1/3]

◆ WarpedCircularDisk() [2/3]

◆ WarpedCircularDisk() [3/3]

◆ ~WarpedCircularDisk()

Member Function Documentation

◆ boundary_triad()

◆ dwdphi()

◆ dwdr()

◆ epsilon()

◆ operator=()

◆ position() [1/2]

◆ position() [2/2]

◆ w()

Member Data Documentation

◆ Epsilon

◆ N

◆ Z_offset