Nektar++
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Static Public Member Functions | Static Public Attributes | Protected Member Functions | Private Attributes | Friends | List of all members
Nektar::SkewSymmetricAdvection Class Reference

#include <SkewSymmetricAdvection.h>

Inheritance diagram for Nektar::SkewSymmetricAdvection:
Inheritance graph
[legend]
Collaboration diagram for Nektar::SkewSymmetricAdvection:
Collaboration graph
[legend]

Static Public Member Functions

static
SolverUtils::AdvectionSharedPtr 		create (std::string)
 Creates an instance of this class.

Static Public Attributes

static std::string className
 Name of class.
static std::string className2

Protected Member Functions

 SkewSymmetricAdvection ()
virtual ~SkewSymmetricAdvection ()
virtual void v_InitObject (LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields)
 Initialises the advection object.
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)
 Advects a vector field.
- Protected Member Functions inherited from Nektar::SolverUtils::Advection
virtual SOLVER_UTILS_EXPORT void v_SetBaseFlow (const Array< OneD, Array< OneD, NekDouble > > &inarray)
 Overrides the base flow used during linearised advection.

Private Attributes

MultiRegions::CoeffState m_CoeffState
bool m_homogen_dealiasing
bool m_SingleMode
bool m_HalfMode

Friends

class MemoryManager< SkewSymmetricAdvection >

Additional Inherited Members

- Public Member Functions inherited from Nektar::SolverUtils::Advection
SOLVER_UTILS_EXPORT void InitObject (LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields)
 Interface function to initialise the advection object.
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)
 Interface function to advect the vector field.
template<typename FuncPointerT , typename ObjectPointerT >
void SetFluxVector (FuncPointerT func, ObjectPointerT obj)
 Set the flux vector callback function.
void SetRiemannSolver (RiemannSolverSharedPtr riemann)
 Set a Riemann solver object for this advection object.
void SetFluxVector (AdvectionFluxVecCB fluxVector)
 Set the flux vector callback function.
void SetBaseFlow (const Array< OneD, Array< OneD, NekDouble > > &inarray)
 Set the base flow used for linearised advection objects.
- Protected Attributes inherited from Nektar::SolverUtils::Advection
AdvectionFluxVecCB m_fluxVector
 Callback function to the flux vector (set when advection is in conservative form).
RiemannSolverSharedPtr m_riemann
 Riemann solver for DG-type schemes.
int m_spaceDim
 Storage for space dimension. Used for homogeneous extension.

Detailed Description

Definition at line 45 of file SkewSymmetricAdvection.h.

Constructor & Destructor Documentation

Nektar::SkewSymmetricAdvection::SkewSymmetricAdvection ( )
protected

Definition at line 48 of file SkewSymmetricAdvection.cpp.

:
Advection()
{
}
Nektar::SkewSymmetricAdvection::~SkewSymmetricAdvection ( )
protectedvirtual

Definition at line 58 of file SkewSymmetricAdvection.cpp.

{
}

Member Function Documentation

static SolverUtils::AdvectionSharedPtr Nektar::SkewSymmetricAdvection::create ( std::string  )
inlinestatic

Creates an instance of this class.

Definition at line 52 of file SkewSymmetricAdvection.h.

{
return MemoryManager<SkewSymmetricAdvection>::AllocateSharedPtr();
}
void Nektar::SkewSymmetricAdvection::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 
)
protectedvirtual

Advects a vector field.

Implements Nektar::SolverUtils::Advection.

Definition at line 82 of file SkewSymmetricAdvection.cpp.

References ASSERTL0, ASSERTL1, Nektar::MultiRegions::DirCartesianMap, m_CoeffState, m_homogen_dealiasing, Vmath::Neg(), Vmath::Smul(), Vmath::Vadd(), Vmath::Vmul(), and Vmath::Vvtvp().

{
int nqtot = fields[0]->GetTotPoints();
ASSERTL1(nConvectiveFields == inarray.num_elements(),"Number of convective fields and Inarray are not compatible");
Array<OneD, NekDouble > Deriv = Array<OneD, NekDouble> (nqtot*nConvectiveFields);
for(int n = 0; n < nConvectiveFields; ++n)
{
// use dimension of Velocity vector to dictate dimension of operation
int ndim = advVel.num_elements();
// ToDo: here we should add a check that V has right dimension
int nPointsTot = fields[0]->GetNpoints();
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:
fields[0]->PhysDeriv(inarray[n],gradV0);
Vmath::Vmul(nPointsTot,gradV0,1,advVel[0],1,outarray[n],1);
Vmath::Vmul(nPointsTot,inarray[n],1,advVel[0],1,gradV0,1);
fields[0]->PhysDeriv(gradV0,tmp);
Vmath::Vadd(nPointsTot,tmp,1,outarray[n],1,outarray[n],1);
Vmath::Smul(nPointsTot,0.5,outarray[n],1,outarray[n],1);
break;
case 2:
gradV1 = Array<OneD, NekDouble> (nPointsTot);
fields[0]->PhysDeriv(inarray[n],gradV0,gradV1);
Vmath::Vmul (nPointsTot,gradV0,1,advVel[0],1,outarray[n],1);
Vmath::Vvtvp(nPointsTot,gradV1,1,advVel[1],1,outarray[n],1,outarray[n],1);
Vmath::Vmul(nPointsTot,inarray[n],1,advVel[0],1,gradV0,1);
Vmath::Vmul(nPointsTot,inarray[n],1,advVel[1],1,gradV1,1);
fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[0],gradV0,tmp);
Vmath::Vadd(nPointsTot,tmp,1,outarray[n],1,outarray[n],1);
fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[1],gradV1,tmp);
Vmath::Vadd(nPointsTot,tmp,1,outarray[n],1,outarray[n],1);
Vmath::Smul(nPointsTot,0.5,outarray[n],1,outarray[n],1);
break;
case 3:
gradV1 = Array<OneD, NekDouble> (nPointsTot);
gradV2 = Array<OneD, NekDouble> (nPointsTot);
fields[0]->PhysDeriv(inarray[n],gradV0,gradV1,gradV2);
//outarray[n] = 1/2(u*du/dx + v*du/dy + w*du/dz + duu/dx + duv/dy + duw/dz)
if(m_homogen_dealiasing == true && fields[0]->GetWaveSpace() == false)
{
fields[0]->DealiasedProd(advVel[0],gradV0,gradV0,m_CoeffState);
fields[0]->DealiasedProd(advVel[1],gradV1,gradV1,m_CoeffState);
fields[0]->DealiasedProd(advVel[2],gradV2,gradV2,m_CoeffState);
Vmath::Vadd(nPointsTot,gradV0,1,gradV1,1,outarray[n],1);
Vmath::Vadd(nPointsTot,gradV2,1,outarray[n],1,outarray[n],1);
fields[0]->DealiasedProd(inarray[n],advVel[0],gradV0,m_CoeffState);
fields[0]->DealiasedProd(inarray[n],advVel[1],gradV1,m_CoeffState);
fields[0]->DealiasedProd(inarray[n],advVel[2],gradV2,m_CoeffState);
fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[0],gradV0,tmp);
Vmath::Vadd(nPointsTot,tmp,1,outarray[n],1,outarray[n],1);
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,0.5,outarray[n],1,outarray[n],1);
}
else if(fields[0]->GetWaveSpace() == true && m_homogen_dealiasing == false)
{
Up = Array<OneD, NekDouble> (nPointsTot);
//vector reused to avoid even more memory requirements
//names may be misleading
fields[0]->HomogeneousBwdTrans(gradV0,tmp);
Vmath::Vmul(nPointsTot,tmp,1,advVel[0],1,outarray[n],1); // + u*du/dx
fields[0]->HomogeneousBwdTrans(gradV1,tmp);
Vmath::Vvtvp(nPointsTot,tmp,1,advVel[1],1,outarray[n],1,outarray[n],1);// + v*du/dy
fields[0]->HomogeneousBwdTrans(gradV2,tmp);
Vmath::Vvtvp(nPointsTot,tmp,1,advVel[2],1,outarray[n],1,outarray[n],1);// + w*du/dz
fields[0]->HomogeneousBwdTrans(inarray[n],Up);
Vmath::Vmul(nPointsTot,Up,1,advVel[0],1,gradV0,1);
Vmath::Vmul(nPointsTot,Up,1,advVel[1],1,gradV1,1);
Vmath::Vmul(nPointsTot,Up,1,advVel[2],1,gradV2,1);
fields[0]->SetWaveSpace(false);
fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[0],gradV0,tmp);//duu/dx
Vmath::Vadd(nPointsTot,tmp,1,outarray[n],1,outarray[n],1);
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,0.5,outarray[n],1,tmp,1);
fields[0]->HomogeneousFwdTrans(tmp,outarray[n]);
}
else if(fields[0]->GetWaveSpace() == false && m_homogen_dealiasing == false)
{
Vmath::Vmul(nPointsTot,gradV0,1,advVel[0],1,outarray[n],1);
Vmath::Vvtvp(nPointsTot,gradV1,1,advVel[1],1,outarray[n],1,outarray[n],1);
Vmath::Vvtvp(nPointsTot,gradV2,1,advVel[2],1,outarray[n],1,outarray[n],1);
Vmath::Vmul(nPointsTot,inarray[n],1,advVel[0],1,gradV0,1);
Vmath::Vmul(nPointsTot,inarray[n],1,advVel[1],1,gradV1,1);
Vmath::Vmul(nPointsTot,inarray[n],1,advVel[2],1,gradV2,1);
fields[0]->PhysDeriv(MultiRegions::DirCartesianMap[0],gradV0,tmp);
Vmath::Vadd(nPointsTot,tmp,1,outarray[n],1,outarray[n],1);
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,0.5,outarray[n],1,outarray[n],1);
}
else
{
ASSERTL0(false, "Dealiasing is not allowed in combination "
"with the Skew-Symmetric advection form for "
"efficiency reasons.");
}
break;
default:
ASSERTL0(false,"dimension unknown");
}
Vmath::Neg(nqtot,outarray[n],1);
}
}
void Nektar::SkewSymmetricAdvection::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
pSessionSession information.
pFieldsExpansion lists for scalar fields.

Reimplemented from Nektar::SolverUtils::Advection.

Definition at line 66 of file SkewSymmetricAdvection.cpp.

References Nektar::MultiRegions::eLocal, m_CoeffState, m_HalfMode, m_homogen_dealiasing, and m_SingleMode.

{
Advection::v_InitObject(pSession, pFields);
m_CoeffState = MultiRegions::eLocal;
m_homogen_dealiasing = pSession->DefinesSolverInfo("dealiasing");
pSession->MatchSolverInfo("ModeType","SingleMode",m_SingleMode,false);
pSession->MatchSolverInfo("ModeType","HalfMode",m_HalfMode,false);
}

Friends And Related Function Documentation

friend class MemoryManager< SkewSymmetricAdvection >
friend

Definition at line 49 of file SkewSymmetricAdvection.h.

Member Data Documentation

string Nektar::SkewSymmetricAdvection::className
static
Initial value:
SolverUtils::GetAdvectionFactory().RegisterCreatorFunction(
"SkewSymmetric",
SkewSymmetricAdvection::create)

Name of class.

Definition at line 57 of file SkewSymmetricAdvection.h.

std::string Nektar::SkewSymmetricAdvection::className2
static

Definition at line 58 of file SkewSymmetricAdvection.h.

MultiRegions::CoeffState Nektar::SkewSymmetricAdvection::m_CoeffState
private

Definition at line 79 of file SkewSymmetricAdvection.h.

Referenced by v_Advect(), and v_InitObject().

bool Nektar::SkewSymmetricAdvection::m_HalfMode
private

Definition at line 82 of file SkewSymmetricAdvection.h.

Referenced by v_InitObject().

bool Nektar::SkewSymmetricAdvection::m_homogen_dealiasing
private

Definition at line 80 of file SkewSymmetricAdvection.h.

Referenced by v_Advect(), and v_InitObject().

bool Nektar::SkewSymmetricAdvection::m_SingleMode
private

Definition at line 81 of file SkewSymmetricAdvection.h.

Referenced by v_InitObject().

Generated on Wed Mar 4 2015 16:33:34 for Nektar++ by   doxygen 1.8.1.2