Nektar::SolverUtils::ALERotate Struct Reference

#include <ALEHelper.h>

Inheritance diagram for Nektar::SolverUtils::ALERotate:

Public Member Functions
	ALERotate (SpatialDomains::ZoneBaseShPtr zone)
void	v_UpdateGridVel (const NekDouble time, Array< OneD, MultiRegions::ExpListSharedPtr > &fields, Array< OneD, Array< OneD, NekDouble > > &gridVelocity) final
Public Member Functions inherited from Nektar::SolverUtils::ALEBase
virtual	~ALEBase ()=default
void	UpdateGridVel (const NekDouble time, Array< OneD, MultiRegions::ExpListSharedPtr > &fields, Array< OneD, Array< OneD, NekDouble > > &gridVelocity)

Private Attributes
SpatialDomains::ZoneRotateShPtr	m_zone

Detailed Description

Definition at line 144 of file ALEHelper.h.

Constructor & Destructor Documentation

ALERotate()

Nektar::SolverUtils::ALERotate::ALERotate

(

SpatialDomains::ZoneBaseShPtr

zone

)

Definition at line 339 of file ALEHelper.cpp.

    : m_zone(std::static_pointer_cast<SpatialDomains::ZoneRotate>(zone))
{
}

Member Function Documentation

v_UpdateGridVel()

void Nektar::SolverUtils::ALERotate::v_UpdateGridVel ( const NekDouble  time, Array< OneD, MultiRegions::ExpListSharedPtr > &  fields, Array< OneD, Array< OneD, NekDouble > > &  gridVelocity  )

finalvirtual

Implements Nektar::SolverUtils::ALEBase.

Definition at line 344 of file ALEHelper.cpp.

{
    boost::ignore_unused(time, fields, gridVelocity);
 
    auto angVel = m_zone->GetAngularVel(time);
    auto axis   = m_zone->GetAxis();
    auto origin = m_zone->GetOrigin();
 
    auto exp = fields[0]->GetExp();
 
    auto elements = m_zone->GetElements();
    for (auto &el : elements)
    {
        int indx       = fields[0]->GetElmtToExpId(el->GetGlobalID());
        int offset     = fields[0]->GetPhys_Offset(indx);
        auto expansion = (*exp)[indx];
 
        int nq = expansion->GetTotPoints();
 
        Array<OneD, NekDouble> xc(nq, 0.0), yc(nq, 0.0), zc(nq, 0.0);
        Array<OneD, NekDouble> norm(3, 0.0);
        expansion->GetCoords(xc, yc, zc);
        for (int i = 0; i < nq; ++i)
        {
            // Vector from origin to point
            NekDouble xpointMinOrigin = xc[i] - origin(0);
            NekDouble ypointMinOrigin = yc[i] - origin(1);
            NekDouble zpointMinOrigin = zc[i] - origin(2);
 
            // Vector orthogonal to plane formed by axis and point
            // We negate angVel here as by convention a positive angular
            // velocity is counter-clockwise
            norm[0] = (ypointMinOrigin * axis[2] - zpointMinOrigin * axis[1]) *
                      (-angVel);
            norm[1] = (zpointMinOrigin * axis[0] - xpointMinOrigin * axis[2]) *
                      (-angVel);
            norm[2] = (xpointMinOrigin * axis[1] - ypointMinOrigin * axis[0]) *
                      (-angVel);
            for (int j = 0; j < gridVelocity.size(); ++j)
            {
                gridVelocity[j][offset + i] = norm[j];
            }
        }
    }
}

References Nektar::MovementTests::axis, m_zone, and Nektar::MovementTests::origin.

Member Data Documentation

m_zone

SpatialDomains::ZoneRotateShPtr Nektar::SolverUtils::ALERotate::m_zone

private

Definition at line 154 of file ALEHelper.h.

Referenced by v_UpdateGridVel().