Nektar::SolverUtils::ALEHelper Class Reference

#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.
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)
SOLVER_UTILS_EXPORT void	ResetMatricesNormal (Array< OneD, Array< OneD, NekDouble > > &traceNormals)
SOLVER_UTILS_EXPORT void	UpdateNormalsFlag ()
const Array< OneD, const Array< OneD, NekDouble > > &	GetGridVelocity ()
bool &	GetUpdateNormalsFlag ()
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)
SOLVER_UTILS_EXPORT void	ExtraFldOutputGrid (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_meshDistorted = false
bool	m_implicitALESolver = false
bool	m_updateNormals = false
NekDouble	m_prevStageTime = 0.0
int	m_spaceDim

Detailed Description

Definition at line 86 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 135 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 152 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().

ExtraFldOutputGrid()

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

Definition at line 497 of file ALEHelper.cpp.

{
    int nCoeffs = m_fieldsALE[0]->GetNcoeffs();
    int nPhys   = m_fieldsALE[0]->GetTotPoints();
    // Adds extra output variables for grid velocity
    std::string gridXName[3] = {"gridX", "gridY", "gridZ"};
    Array<OneD, Array<OneD, NekDouble>> Xc(m_spaceDim);
    for (int i = 0; i < m_spaceDim; ++i)
    {
        Xc[i] = Array<OneD, NekDouble>(nPhys, 0.0);
    }
 
    switch (m_spaceDim)
    {
        case 1:
        {
            m_fieldsALE[0]->GetCoords(Xc[0]);
            break;
        }
        case 2:
        {
            m_fieldsALE[0]->GetCoords(Xc[0], Xc[1]);
            break;
        }
        case 3:
        {
            m_fieldsALE[0]->GetCoords(Xc[0], Xc[1], Xc[2]);
            break;
        }
    }
 
    for (int i = 0; i < m_spaceDim; ++i)
    {
        Array<OneD, NekDouble> coeffs(nCoeffs);
        m_fieldsALE[0]->FwdTransLocalElmt(Xc[i], coeffs);
        fieldcoeffs.emplace_back(coeffs);
        variables.emplace_back(gridXName[i]);
    }
}

References m_fieldsALE, and m_spaceDim.

Referenced by Nektar::CompressibleFlowSystem::v_ExtraFldOutput(), and Nektar::NavierStokesCFE::v_ExtraFldOutput().

ExtraFldOutputGridVelocity()

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

Definition at line 480 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(), Nektar::CompressibleFlowSystem::v_ExtraFldOutput(), and Nektar::NavierStokesCFE::v_ExtraFldOutput().

GetGridVelocity()

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

inline

Definition at line 119 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 232 of file ALEHelper.cpp.

{
    return m_gridVelocityTrace;
}

References m_gridVelocityTrace.

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

GetUpdateNormalsFlag()

bool & Nektar::SolverUtils::ALEHelper::GetUpdateNormalsFlag ( )

inline

Definition at line 124 of file ALEHelper.h.

    {
        return m_updateNormals;
    }

References m_updateNormals.

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);
 
    // Updates trace grid velocity
    for (int i = 0; i < m_gridVelocityTrace.size(); ++i)
    {
        m_fieldsALE[0]->ExtractTracePhys(m_gridVelocity[i],
                                         m_gridVelocityTrace[i], true);
    }
}

References Nektar::SpatialDomains::eFixed, Nektar::SpatialDomains::eNone, Nektar::SpatialDomains::eRotate, Nektar::SpatialDomains::eTranslate, m_ALEs, m_ALESolver, m_fieldsALE, m_gridVelocity, m_gridVelocityTrace, 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 172 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");
 
    ResetMatricesNormal(traceNormals);
    UpdateNormalsFlag();
    // 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], true);
    }
 
    // 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(), m_fieldsALE, m_gridVelocity, m_gridVelocityTrace, m_prevStageTime, ResetMatricesNormal(), Nektar::LibUtilities::Timer::Start(), Nektar::LibUtilities::Timer::Stop(), UpdateNormalsFlag(), and v_UpdateGridVelocity().

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

ResetMatricesNormal()

void Nektar::SolverUtils::ALEHelper::ResetMatricesNormal

(

Array< OneD, Array< OneD, NekDouble > > &

traceNormals

)

Definition at line 209 of file ALEHelper.cpp.

{
    // Now update matrices and normal for each movement zone
    for (auto &ALE : m_ALEs)
    {
        ALE->ResetMatricesNormal(traceNormals, m_fieldsALE);
    }
}

References m_ALEs, and m_fieldsALE.

Referenced by MoveMesh().

UpdateNormalsFlag()

void Nektar::SolverUtils::ALEHelper::UpdateNormalsFlag

(

)

Definition at line 219 of file ALEHelper.cpp.

{
    // check all update nromal definitions and return true if any are true
    for (auto &ALE : m_ALEs)
    {
        if (ALE->UpdateNormalsFlag())
        {
            m_updateNormals = ALE->UpdateNormalsFlag();
            break;
        }
    }
}

References m_ALEs, and m_updateNormals.

Referenced by MoveMesh().

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.

{
}

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

{
    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, 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 97 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 143 of file ALEHelper.h.

Referenced by InitObject(), ResetMatricesNormal(), UpdateNormalsFlag(), and v_UpdateGridVelocity().

m_ALESolver

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

protected

Definition at line 146 of file ALEHelper.h.

Referenced by Nektar::UnsteadyAdvection::DoOdeProjection(), Nektar::CompressibleFlowSystem::DoOdeProjection(), Nektar::CompressibleFlowSystem::GetFluxVector(), Nektar::CompressibleFlowSystem::InitAdvection(), InitObject(), Nektar::UnsteadyAdvection::v_ExtraFldOutput(), Nektar::UnsteadyAdvectionDiffusion::v_ExtraFldOutput(), Nektar::CompressibleFlowSystem::v_ExtraFldOutput(), Nektar::NavierStokesCFE::v_ExtraFldOutput(), and v_UpdateGridVelocity().

m_fieldsALE

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

protected

Definition at line 140 of file ALEHelper.h.

Referenced by ALEDoElmtInvMass(), ALEDoElmtInvMassBwdTrans(), ExtraFldOutputGrid(), ExtraFldOutputGridVelocity(), InitObject(), MoveMesh(), ResetMatricesNormal(), 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 141 of file ALEHelper.h.

Referenced by ExtraFldOutputGridVelocity(), Nektar::UnsteadyAdvection::GetFluxVector(), Nektar::CompressibleFlowSystem::GetFluxVector(), GetGridVelocity(), Nektar::UnsteadyAdvection::GetNormalVel(), InitObject(), 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 142 of file ALEHelper.h.

Referenced by GetGridVelocityTrace(), InitObject(), 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 150 of file ALEHelper.h.

Referenced by Nektar::CFSImplicit::v_ALEInitObject().

m_meshDistorted

bool Nektar::SolverUtils::ALEHelper::m_meshDistorted = false

protected

Definition at line 148 of file ALEHelper.h.

Referenced by Nektar::CompressibleFlowSystem::DoAdvection(), Nektar::CFSImplicit::DoImplicitSolve(), Nektar::UnsteadyAdvection::DoOdeRhs(), Nektar::UnsteadyAdvectionDiffusion::DoOdeRhs(), Nektar::CompressibleFlowSystem::DoOdeRhs(), Nektar::CompressibleFlowSystem::SetBoundaryConditions(), Nektar::CompressibleFlowSystem::v_ALEInitObject(), Nektar::CFSImplicit::v_ALEInitObject(), Nektar::NavierStokesCFE::v_DoDiffusion(), and Nektar::SolverUtils::UnsteadySystem::v_DoSolve().

m_prevStageTime

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

protected

Definition at line 154 of file ALEHelper.h.

Referenced by MoveMesh().

m_spaceDim

int Nektar::SolverUtils::ALEHelper::m_spaceDim

protected

Definition at line 155 of file ALEHelper.h.

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

m_updateNormals

bool Nektar::SolverUtils::ALEHelper::m_updateNormals = false

protected

Definition at line 152 of file ALEHelper.h.

Referenced by GetUpdateNormalsFlag(), and UpdateNormalsFlag().