Nektar::SolverUtils::AdvectionWeakDG Class Reference

#include <AdvectionWeakDG.h>

Inheritance diagram for Nektar::SolverUtils::AdvectionWeakDG:

Public Member Functions
SOLVER_UTILS_EXPORT void	AdvectCoeffs (const int nConvective, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &advVel, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble &time, const Array< OneD, Array< OneD, NekDouble > > &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble > > &pBwd=NullNekDoubleArrayOfArray)
SOLVER_UTILS_EXPORT void	AdvectTraceFlux (const int nConvective, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &advVel, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &TraceFlux, const NekDouble &time, const Array< OneD, Array< OneD, NekDouble > > &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble > > &pBwd=NullNekDoubleArrayOfArray)
SOLVER_UTILS_EXPORT void	AdvectVolumeFlux (const Array< OneD, Array< OneD, NekDouble > > &inarray, TensorOfArray3D< NekDouble > &VolumeFlux)
Public Member Functions inherited from Nektar::SolverUtils::Advection
virtual SOLVER_UTILS_EXPORT	~Advection ()
SOLVER_UTILS_EXPORT void	InitObject (LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields)
	Interface function to initialise the advection object. More...
SOLVER_UTILS_EXPORT void	Advect (const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &advVel, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble &time, const Array< OneD, Array< OneD, NekDouble > > &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble > > &pBwd=NullNekDoubleArrayOfArray)
	Interface function to advect the vector field. More...
template<typename FuncPointerT , typename ObjectPointerT >
void	SetFluxVector (FuncPointerT func, ObjectPointerT obj)
	Set the flux vector callback function. More...
void	SetRiemannSolver (RiemannSolverSharedPtr riemann)
	Set a Riemann solver object for this advection object. More...
void	SetFluxVector (AdvectionFluxVecCB fluxVector)
	Set the flux vector callback function. More...
void	SetBaseFlow (const Array< OneD, Array< OneD, NekDouble > > &inarray, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields)
	Set the base flow used for linearised advection objects. More...
template<typename DataType , typename TypeNekBlkMatSharedPtr >
SOLVER_UTILS_EXPORT void	AddTraceJacToMat (const int nConvectiveFields, const int nSpaceDim, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, TypeNekBlkMatSharedPtr > &TracePntJacCons, Array< OneD, Array< OneD, TypeNekBlkMatSharedPtr > > &gmtxarray, const Array< OneD, TypeNekBlkMatSharedPtr > &TracePntJacGrad, const Array< OneD, Array< OneD, DataType > > &TracePntJacGradSign)
template<typename DataType , typename TypeNekBlkMatSharedPtr >
void	CalcJacobTraceInteg (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const int m, const int n, const Array< OneD, const TypeNekBlkMatSharedPtr > &PntJac, const Array< OneD, const Array< OneD, DataType > > &PntJacSign, Array< OneD, DNekMatSharedPtr > &TraceJacFwd, Array< OneD, DNekMatSharedPtr > &TraceJacBwd)
template<typename DataType , typename TypeNekBlkMatSharedPtr >
void	AddTraceJacToMat (const int nConvectiveFields, const int nSpaceDim, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, TypeNekBlkMatSharedPtr > &TracePntJacCons, Array< OneD, Array< OneD, TypeNekBlkMatSharedPtr > > &gmtxarray, const Array< OneD, TypeNekBlkMatSharedPtr > &TracePntJacGrad, const Array< OneD, Array< OneD, DataType > > &TracePntJacGradSign)

Static Public Member Functions
static AdvectionSharedPtr	create (std::string advType)

Static Public Attributes
static std::string	type

Protected Member Functions
	AdvectionWeakDG ()
void	v_InitObject (LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields) override
	Initialise AdvectionWeakDG objects and store them before starting the time-stepping. More...
void	v_Advect (const int nConvective, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &advVel, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble &time, const Array< OneD, Array< OneD, NekDouble > > &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble > > &pBwd=NullNekDoubleArrayOfArray) override
	Compute the advection term at each time-step using the Discontinuous Galerkin approach (DG). More...
Protected Member Functions inherited from Nektar::SolverUtils::Advection
virtual SOLVER_UTILS_EXPORT void	v_InitObject (LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields)
	Initialises the advection object. More...
virtual SOLVER_UTILS_EXPORT void	v_Advect (const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &advVel, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble &time, const Array< OneD, Array< OneD, NekDouble > > &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble > > &pBwd=NullNekDoubleArrayOfArray)=0
	Advects a vector field. More...
virtual SOLVER_UTILS_EXPORT void	v_SetBaseFlow (const Array< OneD, Array< OneD, NekDouble > > &inarray, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields)
	Overrides the base flow used during linearised advection. More...

Additional Inherited Members
Protected Attributes inherited from Nektar::SolverUtils::Advection
AdvectionFluxVecCB	m_fluxVector
	Callback function to the flux vector (set when advection is in conservative form). More...
RiemannSolverSharedPtr	m_riemann
	Riemann solver for DG-type schemes. More...
int	m_spaceDim
	Storage for space dimension. Used for homogeneous extension. More...

Detailed Description

Definition at line 42 of file AdvectionWeakDG.h.

Constructor & Destructor Documentation

AdvectionWeakDG()

Nektar::SolverUtils::AdvectionWeakDG::AdvectionWeakDG ( )

protected

Definition at line 47 of file AdvectionWeakDG.cpp.

{
}

Referenced by create().

Member Function Documentation

AdvectCoeffs()

void Nektar::SolverUtils::AdvectionWeakDG::AdvectCoeffs	(	const int	nConvective,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	advVel,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		const NekDouble &	time,
		const Array< OneD, Array< OneD, NekDouble > > &	pFwd = NullNekDoubleArrayOfArray,
		const Array< OneD, Array< OneD, NekDouble > > &	pBwd = NullNekDoubleArrayOfArray
	)

Definition at line 117 of file AdvectionWeakDG.cpp.

{
    size_t nPointsTot      = fields[0]->GetTotPoints();
    size_t nCoeffs         = fields[0]->GetNcoeffs();
    size_t nTracePointsTot = fields[0]->GetTrace()->GetTotPoints();
 
    Array<OneD, Array<OneD, Array<OneD, NekDouble>>> fluxvector(
        nConvectiveFields);
    // Allocate storage for flux vector F(u).
    for (int i = 0; i < nConvectiveFields; ++i)
    {
        fluxvector[i] = Array<OneD, Array<OneD, NekDouble>>{size_t(m_spaceDim)};
        for (int j = 0; j < m_spaceDim; ++j)
        {
            fluxvector[i][j] = Array<OneD, NekDouble>(nPointsTot, 0.0);
        }
    }
 
    LibUtilities::Timer timer;
    timer.Start();
    AdvectVolumeFlux(inarray, fluxvector);
    timer.Stop();
    timer.AccumulateRegion("AdvWeakDG:_fluxVector", 3);
 
    timer.Start();
    // Get the advection part (without numerical flux)
    for (int i = 0; i < nConvectiveFields; ++i)
    {
        Vmath::Fill(outarray[i].size(), 0.0, outarray[i], 1);
        fields[i]->IProductWRTDerivBase(fluxvector[i], outarray[i]);
    }
    timer.Stop();
    timer.AccumulateRegion("AdvWeakDG:_IProductWRTDerivBase", 10);
 
    Array<OneD, Array<OneD, NekDouble>> numflux{size_t(nConvectiveFields)};
 
    for (int i = 0; i < nConvectiveFields; ++i)
    {
        numflux[i] = Array<OneD, NekDouble>{nTracePointsTot, 0.0};
    }
 
    AdvectTraceFlux(nConvectiveFields, fields, advVel, inarray, numflux, time,
                    pFwd, pBwd);
 
    // Evaulate <\phi, \hat{F}\cdot n> - OutField[i]
    for (int i = 0; i < nConvectiveFields; ++i)
    {
        Vmath::Neg(nCoeffs, outarray[i], 1);
 
        timer.Start();
        fields[i]->AddTraceIntegral(numflux[i], outarray[i]);
        timer.Stop();
        timer.AccumulateRegion("AdvWeakDG:_AddTraceIntegral", 10);
    }
 
    if (!fields[0]->GetGraph()->GetMovement()->GetMoveFlag() ||
        fields[0]->GetGraph()->GetMovement()->GetImplicitALESolverFlag())
    {
        for (int i = 0; i < nConvectiveFields; ++i)
        {
            timer.Start();
            fields[i]->MultiplyByElmtInvMass(outarray[i], outarray[i]);
            timer.Stop();
            timer.AccumulateRegion("AdvWeakDG:_MultiplyByElmtInvMass", 10);
        }
    }
}

References Nektar::LibUtilities::Timer::AccumulateRegion(), AdvectTraceFlux(), AdvectVolumeFlux(), Vmath::Fill(), Nektar::SolverUtils::Advection::m_spaceDim, Vmath::Neg(), Nektar::LibUtilities::Timer::Start(), and Nektar::LibUtilities::Timer::Stop().

Referenced by v_Advect().

AdvectTraceFlux()

void Nektar::SolverUtils::AdvectionWeakDG::AdvectTraceFlux	(	const int	nConvective,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	advVel,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	TraceFlux,
		const NekDouble &	time,
		const Array< OneD, Array< OneD, NekDouble > > &	pFwd = NullNekDoubleArrayOfArray,
		const Array< OneD, Array< OneD, NekDouble > > &	pBwd = NullNekDoubleArrayOfArray
	)

Definition at line 195 of file AdvectionWeakDG.cpp.

{
    int nTracePointsTot = fields[0]->GetTrace()->GetTotPoints();
 
    ASSERTL1(m_riemann, "Riemann solver must be provided for AdvectionWeakDG.");
 
    // Store forwards/backwards space along trace space
    Array<OneD, Array<OneD, NekDouble>> Fwd(nConvectiveFields);
    Array<OneD, Array<OneD, NekDouble>> Bwd(nConvectiveFields);
 
    if (pFwd == NullNekDoubleArrayOfArray || pBwd == NullNekDoubleArrayOfArray)
    {
        for (int i = 0; i < nConvectiveFields; ++i)
        {
            Fwd[i] = Array<OneD, NekDouble>(nTracePointsTot, 0.0);
            Bwd[i] = Array<OneD, NekDouble>(nTracePointsTot, 0.0);
            fields[i]->GetFwdBwdTracePhys(inarray[i], Fwd[i], Bwd[i]);
        }
    }
    else
    {
        for (int i = 0; i < nConvectiveFields; ++i)
        {
            Fwd[i] = pFwd[i];
            Bwd[i] = pBwd[i];
        }
    }
 
    LibUtilities::Timer timer;
    timer.Start();
    m_riemann->Solve(m_spaceDim, Fwd, Bwd, TraceFlux);
    timer.Stop();
    timer.AccumulateRegion("AdvWeakDG:_Riemann", 10);
}

References Nektar::LibUtilities::Timer::AccumulateRegion(), ASSERTL1, Nektar::SolverUtils::Advection::m_riemann, Nektar::SolverUtils::Advection::m_spaceDim, Nektar::NullNekDoubleArrayOfArray, Nektar::LibUtilities::Timer::Start(), and Nektar::LibUtilities::Timer::Stop().

Referenced by AdvectCoeffs().

AdvectVolumeFlux()

SOLVER_UTILS_EXPORT void Nektar::SolverUtils::AdvectionWeakDG::AdvectVolumeFlux ( const Array< OneD, Array< OneD, NekDouble > > &  inarray, TensorOfArray3D< NekDouble > &  VolumeFlux  )

inline

Definition at line 74 of file AdvectionWeakDG.h.

    {
        m_fluxVector(inarray, VolumeFlux);
    }

References Nektar::SolverUtils::Advection::m_fluxVector.

Referenced by AdvectCoeffs().

create()

static AdvectionSharedPtr Nektar::SolverUtils::AdvectionWeakDG::create ( std::string  advType )

inlinestatic

Definition at line 45 of file AdvectionWeakDG.h.

    {
        return AdvectionSharedPtr(new AdvectionWeakDG());
    }

References AdvectionWeakDG().

v_Advect()

void Nektar::SolverUtils::AdvectionWeakDG::v_Advect ( const int  nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &  fields, const Array< OneD, Array< OneD, NekDouble > > &  advVel, const Array< OneD, Array< OneD, NekDouble > > &  inarray, Array< OneD, Array< OneD, NekDouble > > &  outarray, const NekDouble &  time, const Array< OneD, Array< OneD, NekDouble > > &  pFwd = NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble > > &  pBwd = NullNekDoubleArrayOfArray  )

overrideprotectedvirtual

Compute the advection term at each time-step using the Discontinuous Galerkin approach (DG).

Parameters

nConvectiveFields	Number of fields.
fields	Pointer to fields.
advVel	Advection velocities.
inarray	Solution at the previous time-step.
outarray	Advection term to be passed at the time integration class.

Implements Nektar::SolverUtils::Advection.

Definition at line 76 of file AdvectionWeakDG.cpp.

{
    LibUtilities::Timer timer1;
    timer1.Start();
 
    size_t nCoeffs = fields[0]->GetNcoeffs();
 
    Array<OneD, Array<OneD, NekDouble>> tmp{size_t(nConvectiveFields)};
    for (int i = 0; i < nConvectiveFields; ++i)
    {
        tmp[i] = Array<OneD, NekDouble>{nCoeffs, 0.0};
    }
    LibUtilities::Timer timer;
    timer.Start();
    AdvectCoeffs(nConvectiveFields, fields, advVel, inarray, tmp, time, pFwd,
                 pBwd);
    timer.Stop();
    timer.AccumulateRegion("AdvectCoeffs", 2);
 
    timer.Start();
    for (int i = 0; i < nConvectiveFields; ++i)
    {
        fields[i]->BwdTrans(tmp[i], outarray[i]);
    }
    timer.Stop();
    timer.AccumulateRegion("AdvWeakDG:_BwdTrans", 10);
    timer1.Stop();
    timer1.AccumulateRegion("AdvWeakDG:All", 10);
}

References Nektar::LibUtilities::Timer::AccumulateRegion(), AdvectCoeffs(), Nektar::LibUtilities::Timer::Start(), and Nektar::LibUtilities::Timer::Stop().

v_InitObject()

void Nektar::SolverUtils::AdvectionWeakDG::v_InitObject ( LibUtilities::SessionReaderSharedPtr  pSession, Array< OneD, MultiRegions::ExpListSharedPtr >  pFields  )

overrideprotectedvirtual

Initialise AdvectionWeakDG objects and store them before starting the time-stepping.

Parameters

pSession	Pointer to session reader.
pFields	Pointer to fields.

Reimplemented from Nektar::SolverUtils::Advection.

Definition at line 58 of file AdvectionWeakDG.cpp.

{
    Advection::v_InitObject(pSession, pFields);
}

References Nektar::SolverUtils::Advection::v_InitObject().

Member Data Documentation

type

std::string Nektar::SolverUtils::AdvectionWeakDG::type

static

Initial value:

=
    GetAdvectionFactory().RegisterCreatorFunction(
        "WeakDG", AdvectionWeakDG::create, "Weak DG")

Definition at line 50 of file AdvectionWeakDG.h.