#include <ALEHelper.h>

Inheritance diagram for Nektar::SolverUtils::ALEHelper:

Public Member Functions
virtual	~ALEHelper ()=default

virtual SOLVER_UTILS_EXPORT void	v_ALEInitObject (int spaceDim, Array< OneD, MultiRegions::ExpListSharedPtr > &fields)

SOLVER_UTILS_EXPORT void	InitObject (int spaceDim, Array< OneD, MultiRegions::ExpListSharedPtr > &fields)

virtual SOLVER_UTILS_EXPORT void	v_UpdateGridVelocity (const NekDouble &time)

virtual SOLVER_UTILS_EXPORT void	v_ALEPreMultiplyMass (Array< OneD, Array< OneD, NekDouble > > &fields)

SOLVER_UTILS_EXPORT void	ALEDoElmtInvMass (Array< OneD, Array< OneD, NekDouble > > &traceNormals, Array< OneD, Array< OneD, NekDouble > > &fields, NekDouble time)
	Update m_fields with u^n by multiplying by inverse mass matrix. That's then used in e.g. checkpoint output and L^2 error calculation. More...

SOLVER_UTILS_EXPORT void	ALEDoElmtInvMassBwdTrans (const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray)

SOLVER_UTILS_EXPORT void	MoveMesh (const NekDouble &time, Array< OneD, Array< OneD, NekDouble > > &traceNormals)

const Array< OneD, const Array< OneD, NekDouble > > &	GetGridVelocity ()

SOLVER_UTILS_EXPORT const Array< OneD, const Array< OneD, NekDouble > > &	GetGridVelocityTrace ()

SOLVER_UTILS_EXPORT void	ExtraFldOutputGridVelocity (std::vector< Array< OneD, NekDouble > > &fieldcoeffs, std::vector< std::string > &variables)

Protected Attributes
Array< OneD, MultiRegions::ExpListSharedPtr >	m_fieldsALE

Array< OneD, Array< OneD, NekDouble > >	m_gridVelocity

Array< OneD, Array< OneD, NekDouble > >	m_gridVelocityTrace

std::vector< ALEBaseShPtr >	m_ALEs

bool	m_ALESolver = false

bool	m_ImplicitALESolver = false

NekDouble	m_prevStageTime = 0.0

int	m_spaceDim

Detailed Description

Definition at line 49 of file ALEHelper.h.

Constructor & Destructor Documentation

◆ ~ALEHelper()

virtual Nektar::SolverUtils::ALEHelper::~ALEHelper ( )

virtualdefault

Member Function Documentation

◆ ALEDoElmtInvMass()

void Nektar::SolverUtils::ALEHelper::ALEDoElmtInvMass	(	Array< OneD, Array< OneD, NekDouble > > &	traceNormals,
		Array< OneD, Array< OneD, NekDouble > > &	fields,
		NekDouble	time
	)

Update m_fields with u^n by multiplying by inverse mass matrix. That's then used in e.g. checkpoint output and L^2 error calculation.

Definition at line 131 of file ALEHelper.cpp.

{
    // @TODO: Look at geometric factor and junk only what is needed
    // @TODO: Look at collections and see if they offer a speed up
    for (int i = 0; i < m_fieldsALE.size(); ++i)
    {
        m_fieldsALE[i]->MultiplyByElmtInvMass(
            fields[i],
            m_fieldsALE[i]->UpdateCoeffs()); // @TODO: Potentially matrix free?
        m_fieldsALE[i]->BwdTrans(m_fieldsALE[i]->GetCoeffs(),
                                 m_fieldsALE[i]->UpdatePhys());
    }
}

References m_fieldsALE.

Referenced by Nektar::SolverUtils::UnsteadySystem::v_DoSolve().

◆ ALEDoElmtInvMassBwdTrans()

void Nektar::SolverUtils::ALEHelper::ALEDoElmtInvMassBwdTrans	(	const Array< OneD, const Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	outarray
	)

Definition at line 148 of file ALEHelper.cpp.

{
    const int nc   = m_fieldsALE[0]->GetNcoeffs();
    int nVariables = inarray.size();
 
    // General idea is that we are time-integrating the quantity (Mu), so we
    // need to multiply input by inverse mass matrix to get coefficients u,
    // and then backwards transform to physical space so we can apply the DG
    // operator.
    Array<OneD, NekDouble> tmp(nc);
    for (int i = 0; i < nVariables; ++i)
    {
        outarray[i] = Array<OneD, NekDouble>(m_fieldsALE[0]->GetNpoints());
        m_fieldsALE[i]->MultiplyByElmtInvMass(inarray[i], tmp);
        m_fieldsALE[i]->BwdTrans(tmp, outarray[i]);
    }
}

References m_fieldsALE.

Referenced by Nektar::CFSImplicit::DoImplicitSolve(), Nektar::UnsteadyAdvection::DoOdeRhs(), Nektar::UnsteadyAdvectionDiffusion::DoOdeRhs(), Nektar::CompressibleFlowSystem::DoOdeRhs(), and Nektar::CompressibleFlowSystem::SetBoundaryConditions().

◆ ExtraFldOutputGridVelocity()

void Nektar::SolverUtils::ALEHelper::ExtraFldOutputGridVelocity	(	std::vector< Array< OneD, NekDouble > > &	fieldcoeffs,
		std::vector< std::string > &	variables
	)

Definition at line 389 of file ALEHelper.cpp.

{
    int nCoeffs = m_fieldsALE[0]->GetNcoeffs();
    // Adds extra output variables for grid velocity
    std::string gridVarName[3] = {"gridVx", "gridVy", "gridVz"};
    for (int i = 0; i < m_spaceDim; ++i)
    {
        Array<OneD, NekDouble> gridVel(nCoeffs, 0.0);
        m_fieldsALE[0]->FwdTransLocalElmt(m_gridVelocity[i], gridVel);
        fieldcoeffs.emplace_back(gridVel);
        variables.emplace_back(gridVarName[i]);
    }
}

References m_fieldsALE, m_gridVelocity, and m_spaceDim.

Referenced by Nektar::UnsteadyAdvection::v_ExtraFldOutput(), Nektar::UnsteadyAdvectionDiffusion::v_ExtraFldOutput(), and Nektar::CompressibleFlowSystem::v_ExtraFldOutput().

◆ GetGridVelocity()

const Array< OneD, const Array< OneD, NekDouble > > & Nektar::SolverUtils::ALEHelper::GetGridVelocity ( )

inline

Definition at line 76 of file ALEHelper.h.

    {
        return m_gridVelocity;
    }

References m_gridVelocity.

◆ GetGridVelocityTrace()

const Array< OneD, const Array< OneD, NekDouble > > & Nektar::SolverUtils::ALEHelper::GetGridVelocityTrace ( )

Definition at line 288 of file ALEHelper.cpp.

{
    return m_gridVelocityTrace;
}

References m_gridVelocityTrace.

Referenced by Nektar::CompressibleFlowSystem::InitAdvection().

◆ InitObject()

void Nektar::SolverUtils::ALEHelper::InitObject	(	int	spaceDim,
		Array< OneD, MultiRegions::ExpListSharedPtr > &	fields
	)

Definition at line 48 of file ALEHelper.cpp.

{
    // Create ALE objects for each interface zone
    if (fields[0]->GetGraph() !=
        nullptr) // homogeneous graphs are missing the graph data
    {
        for (auto &zone : fields[0]->GetGraph()->GetMovement()->GetZones())
        {
            switch (zone.second->GetMovementType())
            {
                case SpatialDomains::MovementType::eFixed:
                    m_ALEs.emplace_back(ALEFixedShPtr(
                        MemoryManager<ALEFixed>::AllocateSharedPtr(
                            zone.second)));
                    break;
                case SpatialDomains::MovementType::eTranslate:
                    m_ALEs.emplace_back(ALETranslateShPtr(
                        MemoryManager<ALETranslate>::AllocateSharedPtr(
                            zone.second)));
                    m_ALESolver = true;
                    break;
                case SpatialDomains::MovementType::eRotate:
                    m_ALEs.emplace_back(ALERotateShPtr(
                        MemoryManager<ALERotate>::AllocateSharedPtr(
                            zone.second)));
                    m_ALESolver = true;
                    break;
                case SpatialDomains::MovementType::eNone:
                    WARNINGL0(false,
                              "Zone cannot have movement type of 'None'.");
                default:
                    break;
            }
        }
    }
 
    // Update grid velocity
    v_UpdateGridVelocity(0);
}

References Nektar::SpatialDomains::eFixed, Nektar::SpatialDomains::eNone, Nektar::SpatialDomains::eRotate, Nektar::SpatialDomains::eTranslate, m_ALEs, m_ALESolver, v_UpdateGridVelocity(), and WARNINGL0.

Referenced by Nektar::UnsteadyAdvection::v_ALEInitObject(), Nektar::UnsteadyAdvectionDiffusion::v_ALEInitObject(), Nektar::CompressibleFlowSystem::v_ALEInitObject(), and Nektar::CFSImplicit::v_ALEInitObject().

◆ MoveMesh()

void Nektar::SolverUtils::ALEHelper::MoveMesh	(	const NekDouble &	time,
		Array< OneD, Array< OneD, NekDouble > > &	traceNormals
	)

Definition at line 168 of file ALEHelper.cpp.

{
    // Only move if timestepped
    if (time == m_prevStageTime)
    {
        return;
    }
 
    auto curvedEdges = m_fieldsALE[0]->GetGraph()->GetCurvedEdges();
    auto curvedFaces = m_fieldsALE[0]->GetGraph()->GetCurvedFaces();
 
    LibUtilities::Timer timer;
    timer.Start();
    m_fieldsALE[0]->GetGraph()->GetMovement()->PerformMovement(
        time); // @TODO: Moved out of loop!
    timer.Stop();
    timer.AccumulateRegion("Movement::PerformMovement");
 
    // The order of the resets below is v important to avoid errors
    for (auto &field : m_fieldsALE)
    {
        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.
    for (auto &zone : m_fieldsALE[0]->GetGraph()->GetMovement()->GetZones())
    {
        if (zone.second->GetMoved())
        {
            auto conEl = zone.second->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[m_fieldsALE[0]->GetShapeDimension() -
                                       1]) // This only takes the trace elements
            {
                m_fieldsALE[0]
                    ->GetTrace()
                    ->GetExpFromGeomId(i->GetGlobalID())
                    ->Reset();
            }
        }
    }
 
    for (auto &field : m_fieldsALE)
    {
        for (auto &zone : field->GetGraph()->GetMovement()->GetZones())
        {
            if (zone.second->GetMoved())
            {
                auto conEl = zone.second->GetConstituentElements();
                // Loop over zone elements expansions and rebuild geometric
                // factors
                for (const auto &i :
                     conEl[0]) // This only takes highest dimensioned elements
                {
                    field->GetExpFromGeomId(i->GetGlobalID())->Reset();
                }
            }
        }
    }
 
    for (auto &zone : m_fieldsALE[0]->GetGraph()->GetMovement()->GetZones())
    {
        if (zone.second->GetMoved())
        {
            auto conEl = zone.second->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 = m_fieldsALE[0]
                                 ->GetExpFromGeomId(i->GetGlobalID())
                                 ->GetNtraces();
                for (int j = 0; j < nfaces; ++j)
                {
                    m_fieldsALE[0]
                        ->GetExpFromGeomId(i->GetGlobalID())
                        ->ComputeTraceNormal(j);
                }
            }
        }
    }
 
    for (auto &field : m_fieldsALE)
    {
        // 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
    m_fieldsALE[0]->GetTrace()->GetNormals(traceNormals);
 
    // Recompute grid velocity.
    v_UpdateGridVelocity(time);
 
    // Updates trace grid velocity
    for (int i = 0; i < m_gridVelocityTrace.size(); ++i)
    {
        m_fieldsALE[0]->ExtractTracePhys(m_gridVelocity[i],
                                         m_gridVelocityTrace[i]);
    }
 
    // Set the flag to exchange coords in InterfaceMapDG to true
    m_fieldsALE[0]->GetGraph()->GetMovement()->GetCoordExchangeFlag() = true;
 
    m_prevStageTime = time;
}

References Nektar::LibUtilities::Timer::AccumulateRegion(), Nektar::Collections::eNoCollection, FilterPython_Function::field, m_fieldsALE, m_gridVelocity, m_gridVelocityTrace, m_prevStageTime, Nektar::LibUtilities::Timer::Start(), Nektar::LibUtilities::Timer::Stop(), and v_UpdateGridVelocity().

Referenced by Nektar::UnsteadyAdvection::DoOdeProjection(), and Nektar::CompressibleFlowSystem::DoOdeProjection().

◆ v_ALEInitObject()

void Nektar::SolverUtils::ALEHelper::v_ALEInitObject	(	int	spaceDim,
		Array< OneD, MultiRegions::ExpListSharedPtr > &	fields
	)

virtual

Reimplemented in Nektar::UnsteadyAdvection, Nektar::UnsteadyAdvectionDiffusion, Nektar::CompressibleFlowSystem, and Nektar::CFSImplicit.

Definition at line 42 of file ALEHelper.cpp.

45{

46}

Referenced by Nektar::SolverUtils::UnsteadySystem::v_InitObject().

◆ v_ALEPreMultiplyMass()

void Nektar::SolverUtils::ALEHelper::v_ALEPreMultiplyMass ( Array< OneD, Array< OneD, NekDouble > > & fields )

virtual

Definition at line 108 of file ALEHelper.cpp.

{
    if (m_ALESolver)
    {
        const int nm = m_fieldsALE[0]->GetNcoeffs();
        MultiRegions::GlobalMatrixKey mkey(StdRegions::eMass);
 
        // Premultiply each field by the mass matrix
        for (int i = 0; i < m_fieldsALE.size(); ++i)
        {
            fields[i] = Array<OneD, NekDouble>(nm);
            m_fieldsALE[i]->GeneralMatrixOp(mkey, m_fieldsALE[i]->GetCoeffs(),
                                            fields[i]);
        }
    }
}

References Nektar::StdRegions::eMass, m_ALESolver, and m_fieldsALE.

Referenced by Nektar::SolverUtils::UnsteadySystem::v_DoSolve().

◆ v_UpdateGridVelocity()

void Nektar::SolverUtils::ALEHelper::v_UpdateGridVelocity ( const NekDouble & time )

virtual

Definition at line 90 of file ALEHelper.cpp.

{
 
    // Reset grid velocity to 0
    for (int i = 0; i < m_spaceDim; ++i)
    {
        std::fill(m_gridVelocity[i].begin(), m_gridVelocity[i].end(), 0.0);
    }
    if (m_ALESolver)
    {
        // Now update for each movement zone, adding the grid velocities
        for (auto &ALE : m_ALEs)
        {
            ALE->UpdateGridVel(time, m_fieldsALE, m_gridVelocity);
        }
    }
}

References m_ALEs, m_ALESolver, m_fieldsALE, m_gridVelocity, and m_spaceDim.

Referenced by InitObject(), MoveMesh(), and Nektar::SolverUtils::UnsteadySystem::v_DoInitialise().

Member Data Documentation

◆ m_ALEs

std::vector<ALEBaseShPtr> Nektar::SolverUtils::ALEHelper::m_ALEs

protected

Definition at line 91 of file ALEHelper.h.

Referenced by InitObject(), and v_UpdateGridVelocity().

◆ m_ALESolver

bool Nektar::SolverUtils::ALEHelper::m_ALESolver = false

protected

Definition at line 92 of file ALEHelper.h.

Referenced by Nektar::CompressibleFlowSystem::DoAdvection(), Nektar::CFSImplicit::DoImplicitSolve(), Nektar::UnsteadyAdvection::DoOdeProjection(), Nektar::CompressibleFlowSystem::DoOdeProjection(), Nektar::UnsteadyAdvection::DoOdeRhs(), Nektar::UnsteadyAdvectionDiffusion::DoOdeRhs(), Nektar::CompressibleFlowSystem::DoOdeRhs(), Nektar::CompressibleFlowSystem::GetFluxVector(), Nektar::CompressibleFlowSystem::InitAdvection(), InitObject(), Nektar::CompressibleFlowSystem::SetBoundaryConditions(), v_ALEPreMultiplyMass(), Nektar::NavierStokesCFE::v_DoDiffusion(), Nektar::SolverUtils::UnsteadySystem::v_DoSolve(), Nektar::UnsteadyAdvection::v_ExtraFldOutput(), Nektar::UnsteadyAdvectionDiffusion::v_ExtraFldOutput(), Nektar::CompressibleFlowSystem::v_ExtraFldOutput(), and v_UpdateGridVelocity().

◆ m_fieldsALE

Array<OneD, MultiRegions::ExpListSharedPtr> Nektar::SolverUtils::ALEHelper::m_fieldsALE

protected

Definition at line 88 of file ALEHelper.h.

Referenced by ALEDoElmtInvMass(), ALEDoElmtInvMassBwdTrans(), ExtraFldOutputGridVelocity(), MoveMesh(), Nektar::UnsteadyAdvection::v_ALEInitObject(), Nektar::UnsteadyAdvectionDiffusion::v_ALEInitObject(), Nektar::CompressibleFlowSystem::v_ALEInitObject(), Nektar::CFSImplicit::v_ALEInitObject(), v_ALEPreMultiplyMass(), and v_UpdateGridVelocity().

◆ m_gridVelocity

Array<OneD, Array<OneD, NekDouble> > Nektar::SolverUtils::ALEHelper::m_gridVelocity

protected

Definition at line 89 of file ALEHelper.h.

Referenced by ExtraFldOutputGridVelocity(), Nektar::UnsteadyAdvection::GetFluxVector(), Nektar::CompressibleFlowSystem::GetFluxVector(), GetGridVelocity(), Nektar::UnsteadyAdvection::GetNormalVel(), MoveMesh(), Nektar::UnsteadyAdvection::v_ALEInitObject(), Nektar::UnsteadyAdvectionDiffusion::v_ALEInitObject(), Nektar::CompressibleFlowSystem::v_ALEInitObject(), Nektar::CFSImplicit::v_ALEInitObject(), Nektar::CompressibleFlowSystem::v_GetMaxStdVelocity(), and v_UpdateGridVelocity().

◆ m_gridVelocityTrace

Array<OneD, Array<OneD, NekDouble> > Nektar::SolverUtils::ALEHelper::m_gridVelocityTrace

protected

Definition at line 90 of file ALEHelper.h.

Referenced by GetGridVelocityTrace(), Nektar::NavierStokesCFE::InitObject_Explicit(), MoveMesh(), Nektar::UnsteadyAdvection::v_ALEInitObject(), Nektar::UnsteadyAdvectionDiffusion::v_ALEInitObject(), Nektar::CompressibleFlowSystem::v_ALEInitObject(), Nektar::CFSImplicit::v_ALEInitObject(), and Nektar::CompressibleFlowSystem::v_InitObject().

◆ m_ImplicitALESolver

bool Nektar::SolverUtils::ALEHelper::m_ImplicitALESolver = false

protected

Definition at line 93 of file ALEHelper.h.

Referenced by Nektar::CompressibleFlowSystem::SetBoundaryConditions(), and Nektar::CFSImplicit::v_ALEInitObject().

◆ m_prevStageTime

NekDouble Nektar::SolverUtils::ALEHelper::m_prevStageTime = 0.0

protected

Definition at line 94 of file ALEHelper.h.

Referenced by MoveMesh().

◆ m_spaceDim

int Nektar::SolverUtils::ALEHelper::m_spaceDim

protected

Definition at line 95 of file ALEHelper.h.

Referenced by ExtraFldOutputGridVelocity(), Nektar::UnsteadyAdvection::v_ALEInitObject(), Nektar::UnsteadyAdvectionDiffusion::v_ALEInitObject(), Nektar::CompressibleFlowSystem::v_ALEInitObject(), Nektar::CFSImplicit::v_ALEInitObject(), and v_UpdateGridVelocity().

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ~ALEHelper()

Member Function Documentation

◆ ALEDoElmtInvMass()

◆ ALEDoElmtInvMassBwdTrans()

◆ ExtraFldOutputGridVelocity()

◆ GetGridVelocity()

◆ GetGridVelocityTrace()

◆ InitObject()

◆ MoveMesh()

◆ v_ALEInitObject()

◆ v_ALEPreMultiplyMass()

◆ v_UpdateGridVelocity()

Member Data Documentation

◆ m_ALEs

◆ m_ALESolver

◆ m_fieldsALE

◆ m_gridVelocity

◆ m_gridVelocityTrace

◆ m_ImplicitALESolver

◆ m_prevStageTime

◆ m_spaceDim