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
bool	v_UpdateNormalsFlag () final
void	v_ResetMatricesNormal (Array< OneD, Array< OneD, NekDouble > > &traceNormals, Array< OneD, MultiRegions::ExpListSharedPtr > &fields) 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)
void	ResetMatricesNormal (Array< OneD, Array< OneD, NekDouble > > &traceNormals, Array< OneD, MultiRegions::ExpListSharedPtr > &fields)
bool	UpdateNormalsFlag ()

Private Attributes
SpatialDomains::ZoneRotateShPtr	m_zone

Additional Inherited Members
Public Attributes inherited from Nektar::SolverUtils::ALEBase
bool	m_meshDistorted = false

Detailed Description

Definition at line 202 of file ALEHelper.h.

Constructor & Destructor Documentation

ALERotate()

Nektar::SolverUtils::ALERotate::ALERotate

(

SpatialDomains::ZoneBaseShPtr

zone

)

Definition at line 346 of file ALEHelper.cpp.

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

Member Function Documentation

v_ResetMatricesNormal()

void Nektar::SolverUtils::ALERotate::v_ResetMatricesNormal ( Array< OneD, Array< OneD, NekDouble > > &  traceNormals, Array< OneD, MultiRegions::ExpListSharedPtr > &  fields  )

finalvirtual

Implements Nektar::SolverUtils::ALEBase.

Definition at line 399 of file ALEHelper.cpp.

{
 
    auto curvedEdges = fields[0]->GetGraph()->GetCurvedEdges();
    auto curvedFaces = fields[0]->GetGraph()->GetCurvedFaces();
 
    // The order of the resets below is v important to avoid errors
    if (m_meshDistorted)
    {
        for (auto &field : fields)
        {
            field->ResetMatrices();
        }
    }
 
    // Loop over all elements and faces and edges and reset geometry
    // information. Only need to do this on the first field as the geometry
    // information is shared.
    if (m_zone->GetMoved() && m_zone->v_GetAngle() != 0)
    {
        auto conEl = m_zone->GetConstituentElements();
        for (const auto &i : conEl)
        {
            for (const auto &j : i)
            {
                j->ResetNonRecursive(curvedEdges, curvedFaces);
            }
        }
 
        // We need to rebuild geometric factors on the trace elements
        for (const auto &i : conEl[fields[0]->GetShapeDimension() -
                                   1]) // This only takes the trace elements
        {
            fields[0]->GetTrace()->GetExpFromGeomId(i->GetGlobalID())->Reset();
        }
    }
 
    if (m_zone->GetMoved() && m_zone->v_GetAngle() != 0)
    {
        auto conEl = m_zone->GetConstituentElements();
        // Loop over zone elements expansions and rebuild geometric
        // factors
        for (const auto &i :
             conEl[0]) // This only takes highest dimensioned elements
        {
            fields[0]->GetExpFromGeomId(i->GetGlobalID())->Reset();
        }
    }
 
    if (m_zone->GetMoved() && m_zone->v_GetAngle() != 0)
    {
        auto conEl = m_zone->GetConstituentElements();
        // Loop over zone elements expansions and rebuild geometric factors
        // and recalc trace normals
        for (const auto &i :
             conEl[0]) // This only takes highest dimensioned elements
        {
            int nfaces =
                fields[0]->GetExpFromGeomId(i->GetGlobalID())->GetNtraces();
            for (int j = 0; j < nfaces; ++j)
            {
                fields[0]
                    ->GetExpFromGeomId(i->GetGlobalID())
                    ->ComputeTraceNormal(j);
            }
        }
    }
 
    for (auto &field : fields)
    {
        // Reset collections (despite the default being eNoCollection it does
        // remember the last auto-tuned values), eNoImpType gives lots of output
        field->CreateCollections(Collections::eNoCollection);
    }
 
    // Reload new trace normals in to the solver cache
    fields[0]->GetTrace()->GetNormals(traceNormals);
}

References Nektar::Collections::eNoCollection, Nektar::SolverUtils::ALEBase::m_meshDistorted, and m_zone.

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 351 of file ALEHelper.cpp.

{
 
    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 m_zone.

v_UpdateNormalsFlag()

bool Nektar::SolverUtils::ALERotate::v_UpdateNormalsFlag ( )

inlinefinalvirtual

Implements Nektar::SolverUtils::ALEBase.

Definition at line 211 of file ALEHelper.h.

    {
        return true;
    }

Member Data Documentation

m_zone

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

private

Definition at line 221 of file ALEHelper.h.

Referenced by v_ResetMatricesNormal(), and v_UpdateGridVel().