Nektar::SpatialDomains::ZoneRotate Struct Reference

Rotating zone: Motion of every point around a given axis on an origin. More...

#include <Zones.h>

Inheritance diagram for Nektar::SpatialDomains::ZoneRotate:

Public Member Functions
	ZoneRotate (int id, const CompositeMap &domain, const int coordDim, const NekPoint< NekDouble > &origin, const DNekVec &axis, const LibUtilities::EquationSharedPtr &angularVelEqn)
virtual	~ZoneRotate ()=default
	Default destructor. More...
NekDouble	GetAngularVel (NekDouble &time) const
	Return the angular velocity of the zone at. More...
Public Member Functions inherited from Nektar::SpatialDomains::ZoneBase
	ZoneBase (MovementType type, int indx, CompositeMap domain, int coordDim)
	Constructor. More...
virtual	~ZoneBase ()=default
	Default destructor. More...
MovementType	GetMovementType () const
	Returns the type of movement. More...
CompositeMap	GetDomain () const
	Returns the domain the zone is on. More...
int &	GetId ()
	Returns the zone ID. More...
bool	Move (NekDouble time)
	Performs the movement of the zone at. More...
std::vector< GeometrySharedPtr > const &	GetElements () const
	Returns all highest dimension elements in the zone. More...
bool &	GetMoved ()
	Returns the flag which states if the zone has moved in this timestep. More...
void	ClearBoundingBoxes ()
	Clears all bounding boxes associated with the zones elements. More...

Protected Member Functions
virtual bool	v_Move (NekDouble time) final
	Virtual function for movement of the zone at. More...
virtual bool	v_Move (NekDouble time)
	Virtual function for movement of the zone at. More...

Protected Attributes
NekPoint< NekDouble >	m_origin
	Origin point rotation is performed around. More...
DNekVec	m_axis
	Axis rotation is performed around. More...
LibUtilities::EquationSharedPtr	m_angularVelEqn
	Equation defining angular velocity as a function of time. More...
DNekMat	m_W = DNekMat(3, 3, 0.0)
	W matrix Rodrigues' rotation formula, cross product of axis. More...
DNekMat	m_W2 = DNekMat(3, 3, 0.0)
	W^2 matrix Rodrigues' rotation formula, cross product of axis squared. More...
Protected Attributes inherited from Nektar::SpatialDomains::ZoneBase
MovementType	m_type = MovementType::eNone
	Type of zone movement. More...
int	m_id
	Zone ID. More...
CompositeMap	m_domain
	Zone domain. More...
std::vector< GeometrySharedPtr >	m_elements
	Vector of highest dimension zone elements. More...
bool	m_moved = true
	Moved flag. More...
int	m_coordDim
	Coordinate dimension. More...
std::vector< PointGeomSharedPtr >	m_verts
	Vector of all points in the zone. More...
std::vector< CurveSharedPtr >	m_curves
	Vector of all curves in the zone. More...
std::vector< PointGeom >	m_origVerts
	Vector of all points in the zone at initialisation. More...

Detailed Description

Rotating zone: Motion of every point around a given axis on an origin.

Definition at line 140 of file Zones.h.

Constructor & Destructor Documentation

ZoneRotate()

Nektar::SpatialDomains::ZoneRotate::ZoneRotate	(	int	id,
		const CompositeMap &	domain,
		const int	coordDim,
		const NekPoint< NekDouble > &	origin,
		const DNekVec &	axis,
		const LibUtilities::EquationSharedPtr &	angularVelEqn
	)

Constructor for rotating zones

Parameters

id	Zone ID
domain	Domain that the zone consists of
coordDim	Coordinate dimension
origin	Origin that the zone rotates about
axis	Axis that the zone rotates about
angularVelEqn	Equation for the angular velocity of rotation

Definition at line 141 of file Zones.cpp.

    : ZoneBase(MovementType::eRotate, id, domain, coordDim), m_origin(origin),
      m_axis(axis), m_angularVelEqn(angularVelEqn)
{
    // Construct rotation matrix
    m_W(0, 1) = -m_axis[2];
    m_W(0, 2) = m_axis[1];
    m_W(1, 0) = m_axis[2];
    m_W(1, 2) = -m_axis[0];
    m_W(2, 0) = -m_axis[1];
    m_W(2, 1) = m_axis[0];
 
    m_W2 = m_W * m_W;
}

References Nektar::SpatialDomains::eRotate, m_axis, m_W, and m_W2.

~ZoneRotate()

virtual Nektar::SpatialDomains::ZoneRotate::~ZoneRotate ( )

virtualdefault

Default destructor.

Member Function Documentation

GetAngularVel()

NekDouble Nektar::SpatialDomains::ZoneRotate::GetAngularVel

(

NekDouble &

time

)

const

Return the angular velocity of the zone at.

Parameters

time

Definition at line 180 of file Zones.cpp.

{
    return m_angularVelEqn->Evaluate(0, 0, 0, time);
}

References m_angularVelEqn.

Referenced by v_Move().

v_Move()

bool Nektar::SpatialDomains::ZoneRotate::v_Move ( NekDouble  time )

finalprotectedvirtual

Virtual function for movement of the zone at.

Parameters

time

Reimplemented from Nektar::SpatialDomains::ZoneBase.

Definition at line 186 of file Zones.cpp.

{
    // Currently only valid for constant angular velocity
    NekDouble angle = GetAngularVel(time) * time;
 
    // Identity matrix
    DNekMat rot(3, 3, 0.0);
    rot(0, 0) = 1.0;
    rot(1, 1) = 1.0;
    rot(2, 2) = 1.0;
 
    // Rodrigues' rotation formula in matrix form
    rot = rot + sin(angle) * m_W + (1 - cos(angle)) * m_W2;
 
    int cnt = 0;
    for (auto &vert : m_verts)
    {
        NekPoint<NekDouble> pnt = m_origVerts[cnt] - m_origin;
        DNekVec pntVec          = {pnt[0], pnt[1], pnt[2]};
 
        DNekVec newLoc = rot * pntVec;
 
        vert->UpdatePosition(newLoc(0) + m_origin[0], newLoc(1) + m_origin[1],
                             newLoc(2) + m_origin[2]);
        cnt++;
    }
 
    for (auto &curve : m_curves)
    {
        for (auto &vert : curve->m_points)
        {
            NekPoint<NekDouble> pnt = m_origVerts[cnt] - m_origin;
            DNekVec pntVec          = {pnt[0], pnt[1], pnt[2]};
 
            DNekVec newLoc = rot * pntVec;
 
            vert->UpdatePosition(newLoc(0) + m_origin[0],
                                 newLoc(1) + m_origin[1],
                                 newLoc(2) + m_origin[2]);
            cnt++;
        }
    }
 
    ClearBoundingBoxes();
 
    return true;
}

References Nektar::SpatialDomains::ZoneBase::ClearBoundingBoxes(), GetAngularVel(), Nektar::SpatialDomains::ZoneBase::m_curves, m_origin, Nektar::SpatialDomains::ZoneBase::m_origVerts, Nektar::SpatialDomains::ZoneBase::m_verts, m_W, and m_W2.

Member Data Documentation

m_angularVelEqn

LibUtilities::EquationSharedPtr Nektar::SpatialDomains::ZoneRotate::m_angularVelEqn

protected

Equation defining angular velocity as a function of time.

Definition at line 168 of file Zones.h.

Referenced by GetAngularVel().

m_axis

DNekVec Nektar::SpatialDomains::ZoneRotate::m_axis

protected

Axis rotation is performed around.

Definition at line 166 of file Zones.h.

Referenced by ZoneRotate().

m_origin

NekPoint<NekDouble> Nektar::SpatialDomains::ZoneRotate::m_origin

protected

Origin point rotation is performed around.

Definition at line 164 of file Zones.h.

Referenced by v_Move().

m_W

DNekMat Nektar::SpatialDomains::ZoneRotate::m_W = DNekMat(3, 3, 0.0)

protected

W matrix Rodrigues' rotation formula, cross product of axis.

Definition at line 170 of file Zones.h.

Referenced by v_Move(), and ZoneRotate().

m_W2

DNekMat Nektar::SpatialDomains::ZoneRotate::m_W2 = DNekMat(3, 3, 0.0)

protected

W^2 matrix Rodrigues' rotation formula, cross product of axis squared.

Definition at line 172 of file Zones.h.

Referenced by v_Move(), and ZoneRotate().