#include <AlternateSkewAdvection.h>

Inheritance diagram for Nektar::AlternateSkewAdvection:

Static Public Member Functions
static SolverUtils::AdvectionSharedPtr	create (std::string)
	Creates an instance of this class. More...

Static Public Attributes
static std::string	className
	Name of class. More...

static std::string	className2

Protected Member Functions
	AlternateSkewAdvection ()

virtual	~AlternateSkewAdvection ()

virtual void	v_InitObject (LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields)
	Initialises the advection object. More...

virtual 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)
	Advects a vector field. More...

Protected Member Functions inherited from Nektar::SolverUtils::Advection
virtual SOLVER_UTILS_EXPORT void	v_AdvectVolumeFlux (const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &advVel, const Array< OneD, Array< OneD, NekDouble >> &inarray, TensorOfArray3D< NekDouble > &pVolumeFlux, const NekDouble &time)
	Advects Volume Flux. More...

virtual SOLVER_UTILS_EXPORT void	v_AdvectTraceFlux (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 >> &pTraceFlux, const NekDouble &time, const Array< OneD, Array< OneD, NekDouble >> &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble >> &pBwd=NullNekDoubleArrayOfArray)
	Advects Trace Flux. More...

virtual SOLVER_UTILS_EXPORT void	v_AdvectCoeffs (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)

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

virtual SOLVER_UTILS_EXPORT void	v_AddVolumJacToMat (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const int &nConvectiveFields, const TensorOfArray5D< NekDouble > &ElmtJacArray, Array< OneD, Array< OneD, SNekBlkMatSharedPtr >> &gmtxarray)

Private Attributes
int	m_advectioncalls

bool	m_SingleMode

bool	m_HalfMode

Friends
class	MemoryManager< AlternateSkewAdvection >

Additional Inherited Members
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...

SOLVER_UTILS_EXPORT void	AdvectVolumeFlux (const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble >> &pAdvVel, const Array< OneD, Array< OneD, NekDouble >> &pInarray, TensorOfArray3D< NekDouble > &pVolumeFlux, const NekDouble &pTime)
	Interface function to advect the Volume field. More...

SOLVER_UTILS_EXPORT void	AdvectTraceFlux (const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble >> &pAdvVel, const Array< OneD, Array< OneD, NekDouble >> &pInarray, Array< OneD, Array< OneD, NekDouble >> &pTraceFlux, const NekDouble &pTime, const Array< OneD, Array< OneD, NekDouble >> &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble >> &pBwd=NullNekDoubleArrayOfArray)
	Interface function to advect the Trace field. More...

SOLVER_UTILS_EXPORT void	AdvectCoeffs (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)
	Similar with Advection::Advect(): calculate the advection flux The difference is in the outarray: it is the coefficients of basis for AdvectCoeffs() it is the physics on quadrature points for Advect() 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)

SOLVER_UTILS_EXPORT void	AddVolumJacToMat (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const int &nConvectiveFields, const TensorOfArray5D< NekDouble > &ElmtJacArray, Array< OneD, Array< OneD, SNekBlkMatSharedPtr >> &gmtxarray)

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)

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 43 of file AlternateSkewAdvection.h.

Constructor & Destructor Documentation

◆ AlternateSkewAdvection()

Nektar::AlternateSkewAdvection::AlternateSkewAdvection ( )

protected

Constructor. Creates ...

Parameters

Definition at line 53 of file AlternateSkewAdvection.cpp.

                                                : Advection()
 {
 }

◆ ~AlternateSkewAdvection()

Nektar::AlternateSkewAdvection::~AlternateSkewAdvection ( )

protectedvirtual

Definition at line 57 of file AlternateSkewAdvection.cpp.

58 {

59 }

Member Function Documentation

◆ create()

static SolverUtils::AdvectionSharedPtr Nektar::AlternateSkewAdvection::create ( std::string )

inlinestatic

Creates an instance of this class.

Definition at line 50 of file AlternateSkewAdvection.h.

     {
         return MemoryManager<AlternateSkewAdvection>::AllocateSharedPtr();
     }

References Nektar::MemoryManager< DataType >::AllocateSharedPtr().

◆ v_Advect()

void Nektar::AlternateSkewAdvection::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`
	)

protectedvirtual

Advects a vector field.

Implements Nektar::SolverUtils::Advection.

Definition at line 69 of file AlternateSkewAdvection.cpp.

 {
     // use dimension of Velocity vector to dictate dimension of operation
     int ndim       = advVel.size();
     int nPointsTot = fields[0]->GetNpoints();
     Array<OneD, Array<OneD, NekDouble>> velocity(ndim);
     for (int i = 0; i < ndim; ++i)
     {
         if (fields[i]->GetWaveSpace() && !m_SingleMode && !m_HalfMode)
         {
             velocity[i] = Array<OneD, NekDouble>(nPointsTot, 0.0);
             fields[i]->HomogeneousBwdTrans(advVel[i], velocity[i]);
         }
         else
         {
             velocity[i] = advVel[i];
         }
     }
     for (int n = 0; n < nConvectiveFields; ++n)
     {
         // ToDo: here we should add a check that V has right dimension
         Array<OneD, NekDouble> gradV0, gradV1, gradV2, tmp, Up;
  
         gradV0 = Array<OneD, NekDouble>(nPointsTot);
         tmp    = Array<OneD, NekDouble>(nPointsTot);
  
         // Evaluate V\cdot Grad(u)
         switch (ndim)
         {
             case 1:
                 if (m_advectioncalls % 2 == 0)
                 {
                     fields[0]->PhysDeriv(inarray[n], gradV0);
                     Vmath::Vmul(nPointsTot, gradV0, 1, velocity[0], 1,
                                 outarray[n], 1);
                 }
                 else
                 {
                     Vmath::Vmul(nPointsTot, inarray[n], 1, velocity[0], 1,
                                 gradV0, 1);
                     fields[0]->PhysDeriv(gradV0, outarray[n]);
                 }
                 Vmath::Smul(nPointsTot, 0.5, outarray[n], 1, outarray[n],
                             1); // must be mult by 0.5????
                 break;
             case 2:
                 gradV1 = Array<OneD, NekDouble>(nPointsTot);
                 if (m_advectioncalls % 2 == 0)
                 {
                     fields[0]->PhysDeriv(inarray[n], gradV0, gradV1);
                     Vmath::Vmul(nPointsTot, gradV0, 1, velocity[0], 1,
                                 outarray[n], 1);
                     Vmath::Vvtvp(nPointsTot, gradV1, 1, velocity[1], 1,
                                  outarray[n], 1, outarray[n], 1);
                 }
                 else
                 {
                     Vmath::Vmul(nPointsTot, inarray[n], 1, velocity[0], 1,
                                 gradV0, 1);
                     Vmath::Vmul(nPointsTot, inarray[n], 1, velocity[1], 1,
                                 gradV1, 1);
                     fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[0],
                                          gradV0, outarray[n]);
                     fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[1],
                                          gradV1, tmp);
                     Vmath::Vadd(nPointsTot, tmp, 1, outarray[n], 1, outarray[n],
                                 1);
                 }
                 Vmath::Smul(nPointsTot, 1.0, outarray[n], 1, outarray[n],
                             1); // must be mult by 0.5????
                 break;
             case 3:
                 gradV1 = Array<OneD, NekDouble>(nPointsTot);
                 gradV2 = Array<OneD, NekDouble>(nPointsTot);
  
                 // outarray[n] = 1/2(u*du/dx + v*du/dy + w*du/dz + duu/dx +
                 // duv/dy + duw/dz)
  
                 if (fields[0]->GetWaveSpace() == true)
                 {
                     if (m_advectioncalls % 2 == 0)
                     {
                         // vector reused to avoid even more memory requirements
                         // names may be misleading
                         fields[0]->PhysDeriv(inarray[n], gradV0, gradV1,
                                              gradV2);
                         fields[0]->HomogeneousBwdTrans(gradV0, tmp);
                         Vmath::Vmul(nPointsTot, tmp, 1, velocity[0], 1,
                                     outarray[n], 1); // + u*du/dx
                         fields[0]->HomogeneousBwdTrans(gradV1, tmp);
                         Vmath::Vvtvp(nPointsTot, tmp, 1, velocity[1], 1,
                                      outarray[n], 1, outarray[n],
                                      1); // + v*du/dy
                         fields[0]->HomogeneousBwdTrans(gradV2, tmp);
                         Vmath::Vvtvp(nPointsTot, tmp, 1, velocity[2], 1,
                                      outarray[n], 1, outarray[n],
                                      1); // + w*du/dz
                     }
                     else
                     {
                         Up = Array<OneD, NekDouble>(nPointsTot);
                         fields[0]->HomogeneousBwdTrans(inarray[n], Up);
                         Vmath::Vmul(nPointsTot, Up, 1, velocity[0], 1, gradV0,
                                     1);
                         Vmath::Vmul(nPointsTot, Up, 1, velocity[1], 1, gradV1,
                                     1);
                         Vmath::Vmul(nPointsTot, Up, 1, velocity[2], 1, gradV2,
                                     1);
  
                         fields[0]->SetWaveSpace(false);
                         fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[0],
                                              gradV0, outarray[n]); // duu/dx
                         fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[1],
                                              gradV1, tmp); // duv/dy
                         Vmath::Vadd(nPointsTot, tmp, 1, outarray[n], 1,
                                     outarray[n], 1);
                         fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[2],
                                              gradV2, tmp); // duw/dz
                         Vmath::Vadd(nPointsTot, tmp, 1, outarray[n], 1,
                                     outarray[n], 1);
                         fields[0]->SetWaveSpace(true);
                     }
  
                     Vmath::Smul(nPointsTot, 1.0, outarray[n], 1, tmp,
                                 1); // must be mult by 0.5????
                     fields[0]->HomogeneousFwdTrans(tmp, outarray[n]);
                 }
                 else
                 {
                     if (m_advectioncalls % 2 == 0)
                     {
                         fields[0]->PhysDeriv(inarray[n], gradV0, gradV1,
                                              gradV2);
                         Vmath::Vmul(nPointsTot, gradV0, 1, velocity[0], 1,
                                     outarray[n], 1);
                         Vmath::Vvtvp(nPointsTot, gradV1, 1, velocity[1], 1,
                                      outarray[n], 1, outarray[n], 1);
                         Vmath::Vvtvp(nPointsTot, gradV2, 1, velocity[2], 1,
                                      outarray[n], 1, outarray[n], 1);
                     }
                     else
                     {
                         Vmath::Vmul(nPointsTot, inarray[n], 1, velocity[0], 1,
                                     gradV0, 1);
                         Vmath::Vmul(nPointsTot, inarray[n], 1, velocity[1], 1,
                                     gradV1, 1);
                         Vmath::Vmul(nPointsTot, inarray[n], 1, velocity[2], 1,
                                     gradV2, 1);
                         fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[0],
                                              gradV0, outarray[n]);
                         fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[1],
                                              gradV1, tmp);
                         Vmath::Vadd(nPointsTot, tmp, 1, outarray[n], 1,
                                     outarray[n], 1);
                         fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[2],
                                              gradV2, tmp);
                         Vmath::Vadd(nPointsTot, tmp, 1, outarray[n], 1,
                                     outarray[n], 1);
                     }
                     Vmath::Smul(nPointsTot, 1.0, outarray[n], 1, outarray[n],
                                 1); // must be mult by 0.5????
                 }
                 break;
             default:
                 ASSERTL0(false, "dimension unknown");
         }
     }
 }

References ASSERTL0, Nektar::MultiRegions::DirCartesianMap, m_advectioncalls, m_HalfMode, m_SingleMode, Vmath::Smul(), Vmath::Vadd(), Vmath::Vmul(), and Vmath::Vvtvp().

◆ v_InitObject()

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

protectedvirtual

Initialises the advection object.

This function should be overridden in derived classes to initialise the specific advection data members. However, this base class function should be called as the first statement of the overridden function to ensure the base class is correctly initialised in order.

Parameters

pSession	Session information.
pFields	Expansion lists for scalar fields.

Reimplemented from Nektar::SolverUtils::Advection.

Definition at line 61 of file AlternateSkewAdvection.cpp.

 {
     pSession->MatchSolverInfo("ModeType", "SingleMode", m_SingleMode, false);
     pSession->MatchSolverInfo("ModeType", "HalfMode", m_HalfMode, false);
 }

References m_HalfMode, and m_SingleMode.

Friends And Related Function Documentation

◆ MemoryManager< AlternateSkewAdvection >

friend class MemoryManager< AlternateSkewAdvection >

friend

Definition at line 1 of file AlternateSkewAdvection.h.

Member Data Documentation

◆ className

string Nektar::AlternateSkewAdvection::className

static

Initial value:

=
    SolverUtils::GetAdvectionFactory().RegisterCreatorFunction(
        "AlternateSkew", AlternateSkewAdvection::create,
        "Alternating Skew Symmetric")

Name of class.

Definition at line 55 of file AlternateSkewAdvection.h.

◆ className2

std::string Nektar::AlternateSkewAdvection::className2

static

Definition at line 56 of file AlternateSkewAdvection.h.

◆ m_advectioncalls

int Nektar::AlternateSkewAdvection::m_advectioncalls

private

Definition at line 79 of file AlternateSkewAdvection.h.

Referenced by v_Advect().

◆ m_HalfMode

bool Nektar::AlternateSkewAdvection::m_HalfMode

private

Definition at line 81 of file AlternateSkewAdvection.h.

Referenced by v_Advect(), and v_InitObject().

◆ m_SingleMode

bool Nektar::AlternateSkewAdvection::m_SingleMode

private

Definition at line 80 of file AlternateSkewAdvection.h.

Referenced by v_Advect(), and v_InitObject().

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Private Attributes

Friends

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ AlternateSkewAdvection()

◆ ~AlternateSkewAdvection()

Member Function Documentation

◆ create()

◆ v_Advect()

◆ v_InitObject()

Friends And Related Function Documentation

◆ MemoryManager< AlternateSkewAdvection >

Member Data Documentation

◆ className

◆ className2

◆ m_advectioncalls

◆ m_HalfMode

◆ m_SingleMode