#include <Diffusion.h>

Inheritance diagram for Nektar::SolverUtils::Diffusion:

Public Member Functions
virtual SOLVER_UTILS_EXPORT	~Diffusion ()

SOLVER_UTILS_EXPORT void	InitObject (LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields)

SOLVER_UTILS_EXPORT void	Diffuse (const std::size_t nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const Array< OneD, Array< OneD, NekDouble > > &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble > > &pBwd=NullNekDoubleArrayOfArray)

SOLVER_UTILS_EXPORT void	DiffuseCoeffs (const std::size_t nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const Array< OneD, Array< OneD, NekDouble > > &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble > > &pBwd=NullNekDoubleArrayOfArray)
	Similar with Diffusion::Diffuse(): calculate diffusion flux The difference is in the outarray: it is the coefficients of basis for DiffuseCoeffs() it is the physics on quadrature points for Diffuse() More...

SOLVER_UTILS_EXPORT void	Diffuse (const std::size_t nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, NekDouble time, const Array< OneD, Array< OneD, NekDouble > > &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble > > &pBwd=NullNekDoubleArrayOfArray)

SOLVER_UTILS_EXPORT void	DiffuseCoeffs (const std::size_t nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const Array< OneD, Array< OneD, NekDouble > > &pFwd, const Array< OneD, Array< OneD, NekDouble > > &pBwd, TensorOfArray3D< NekDouble > &qfield, Array< OneD, int > &nonZeroIndex)

SOLVER_UTILS_EXPORT void	DiffuseCoeffs (const std::size_t nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, NekDouble time, const Array< OneD, Array< OneD, NekDouble > > &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble > > &pBwd=NullNekDoubleArrayOfArray)

SOLVER_UTILS_EXPORT void	DiffuseCalcDerivative (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &inarray, TensorOfArray3D< NekDouble > &qfields, const Array< OneD, Array< OneD, NekDouble >> &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble >> &pBwd=NullNekDoubleArrayOfArray)

SOLVER_UTILS_EXPORT void	DiffuseVolumeFlux (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &inarray, TensorOfArray3D< NekDouble > &qfields, TensorOfArray3D< NekDouble > &VolumeFlux, Array< OneD, int > &nonZeroIndex=NullInt1DArray)
	Diffusion Volume FLux. More...

SOLVER_UTILS_EXPORT void	DiffuseTraceFlux (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &inarray, TensorOfArray3D< NekDouble > &qfields, TensorOfArray3D< NekDouble > &VolumeFlux, Array< OneD, Array< OneD, NekDouble > > &TraceFlux, const Array< OneD, Array< OneD, NekDouble >> &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble >> &pBwd=NullNekDoubleArrayOfArray, Array< OneD, int > &nonZeroIndex=NullInt1DArray)
	Diffusion term Trace Flux. More...

SOLVER_UTILS_EXPORT void	DiffuseTraceSymmFlux (const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &inarray, const TensorOfArray3D< NekDouble > &qfield, const TensorOfArray3D< NekDouble > &VolumeFlux, TensorOfArray3D< NekDouble > &SymmFlux, const Array< OneD, Array< OneD, NekDouble >> &pFwd, const Array< OneD, Array< OneD, NekDouble >> &pBwd, Array< OneD, int > &nonZeroIndex, Array< OneD, Array< OneD, NekDouble > > &solution_Aver=NullNekDoubleArrayOfArray, Array< OneD, Array< OneD, NekDouble > > &solution_jump=NullNekDoubleArrayOfArray)

SOLVER_UTILS_EXPORT void	AddDiffusionSymmFluxToCoeff (const std::size_t nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, TensorOfArray3D< NekDouble > &qfield, TensorOfArray3D< NekDouble > &VolumeFlux, Array< OneD, Array< OneD, NekDouble > > &outarray, const Array< OneD, Array< OneD, NekDouble > > &pFwd, const Array< OneD, Array< OneD, NekDouble > > &pBwd)
	Add symmetric flux to field in coeff space. More...

SOLVER_UTILS_EXPORT void	AddDiffusionSymmFluxToPhys (const std::size_t nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, TensorOfArray3D< NekDouble > &qfield, TensorOfArray3D< NekDouble > &VolumeFlux, Array< OneD, Array< OneD, NekDouble > > &outarray, const Array< OneD, Array< OneD, NekDouble > > &pFwd, const Array< OneD, Array< OneD, NekDouble > > &pBwd)
	Add symmetric flux to field in coeff physical space. More...

SOLVER_UTILS_EXPORT void	AddSymmFluxIntegralToOffDiag (const int nConvectiveFields, const int nDim, const int nPts, const int nTracePts, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, const int > &nonZeroIndex, TensorOfArray3D< NekDouble > &Fwdflux, TensorOfArray3D< NekDouble > &Bwdflux, Array< OneD, Array< OneD, NekDouble > > &outarray)

SOLVER_UTILS_EXPORT void	FluxVec (TensorOfArray3D< NekDouble > &fluxvector)

template<typename FuncPointerT , typename ObjectPointerT >
void	SetFluxVector (FuncPointerT func, ObjectPointerT obj)

SOLVER_UTILS_EXPORT void	SetFluxVector (DiffusionFluxVecCB fluxVector)

template<typename FuncPointerT , typename ObjectPointerT >
void	SetFluxVectorNS (FuncPointerT func, ObjectPointerT obj)

void	SetFluxVectorNS (DiffusionFluxVecCBNS fluxVector)

template<typename FuncPointerT , typename ObjectPointerT >
void	SetFluxPenaltyNS (FuncPointerT func, ObjectPointerT obj)

void	SetFluxPenaltyNS (DiffusionFluxPenaltyNS flux)

template<typename FuncPointerT , typename ObjectPointerT >
void	SetDiffusionFluxCons (FuncPointerT func, ObjectPointerT obj)

void	SetDiffusionFluxCons (DiffusionFluxCons flux)

template<typename FuncPointerT , typename ObjectPointerT >
void	SetDiffusionFluxConsTrace (FuncPointerT func, ObjectPointerT obj)

void	SetDiffusionFluxConsTrace (DiffusionFluxCons flux)

template<typename FuncPointerT , typename ObjectPointerT >
void	SetArtificialDiffusionVector (FuncPointerT func, ObjectPointerT obj)

template<typename FuncPointerT , typename ObjectPointerT >
void	SetSpecialBndTreat (FuncPointerT func, ObjectPointerT obj)

template<typename FuncPointerT , typename ObjectPointerT >
void	SetCalcViscosity (FuncPointerT func, ObjectPointerT obj)

template<typename FuncPointerT , typename ObjectPointerT >
void	SetDiffusionSymmFluxCons (FuncPointerT func, ObjectPointerT obj)

void	SetHomoDerivs (Array< OneD, Array< OneD, NekDouble > > &deriv)

virtual TensorOfArray3D< NekDouble > &	GetFluxTensor ()

SOLVER_UTILS_EXPORT void	GetAVmu (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, NekDouble > &muvar, Array< OneD, NekDouble > &MuVarTrace)
	Get the mu of artifical viscosity(AV) More...

SOLVER_UTILS_EXPORT void	ConsVarAveJump (const std::size_t nConvectiveFields, const size_t npnts, const Array< OneD, const Array< OneD, NekDouble > > &vFwd, const Array< OneD, const Array< OneD, NekDouble > > &vBwd, Array< OneD, Array< OneD, NekDouble > > &aver, Array< OneD, Array< OneD, NekDouble > > &jump)
	Get the average and jump value of conservative variables on trace. More...

SOLVER_UTILS_EXPORT const Array< OneD, const Array< OneD, NekDouble > > &	GetTraceNormal ()
	Get trace normal. More...

Public Attributes
Array< OneD, NekDouble >	m_divVel
	Params for Ducros sensor. More...

Array< OneD, NekDouble >	m_divVelSquare

Array< OneD, NekDouble >	m_curlVelSquare

Protected Member Functions
virtual SOLVER_UTILS_EXPORT void	v_InitObject (LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields)

virtual SOLVER_UTILS_EXPORT void	v_Diffuse (const std::size_t nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const Array< OneD, Array< OneD, NekDouble > > &pFwd, const Array< OneD, Array< OneD, NekDouble > > &pBwd)

virtual SOLVER_UTILS_EXPORT void	v_DiffuseCoeffs (const std::size_t nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const Array< OneD, Array< OneD, NekDouble > > &pFwd, const Array< OneD, Array< OneD, NekDouble > > &pBwd)

virtual SOLVER_UTILS_EXPORT void	v_DiffuseCoeffs (const std::size_t nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const Array< OneD, Array< OneD, NekDouble > > &vFwd, const Array< OneD, Array< OneD, NekDouble > > &vBwd, TensorOfArray3D< NekDouble > &qfield, Array< OneD, int > &nonZeroIndex)

virtual SOLVER_UTILS_EXPORT void	v_ConsVarAveJump (const std::size_t nConvectiveFields, const size_t npnts, const Array< OneD, const Array< OneD, NekDouble > > &vFwd, const Array< OneD, const Array< OneD, NekDouble > > &vBwd, Array< OneD, Array< OneD, NekDouble > > &aver, Array< OneD, Array< OneD, NekDouble > > &jump)

virtual SOLVER_UTILS_EXPORT void	v_DiffuseCalcDerivative (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, TensorOfArray3D< NekDouble > &qfields, const Array< OneD, Array< OneD, NekDouble > > &pFwd, const Array< OneD, Array< OneD, NekDouble > > &pBwd)
	Diffusion Flux, calculate the physical derivatives. More...

virtual SOLVER_UTILS_EXPORT void	v_DiffuseVolumeFlux (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &inarray, TensorOfArray3D< NekDouble > &qfields, TensorOfArray3D< NekDouble > &VolumeFlux, Array< OneD, int > &nonZeroIndex)
	Diffusion Volume Flux. More...

virtual SOLVER_UTILS_EXPORT void	v_DiffuseTraceFlux (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &inarray, TensorOfArray3D< NekDouble > &Qfields, TensorOfArray3D< NekDouble > &VolumeFlux, Array< OneD, Array< OneD, NekDouble > > &TraceFlux, const Array< OneD, Array< OneD, NekDouble >> &pFwd, const Array< OneD, Array< OneD, NekDouble >> &pBwd, Array< OneD, int > &nonZeroIndex)
	Diffusion term Trace Flux. More...

virtual SOLVER_UTILS_EXPORT void	v_DiffuseTraceSymmFlux (const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &inarray, const TensorOfArray3D< NekDouble > &qfield, const TensorOfArray3D< NekDouble > &VolumeFlux, TensorOfArray3D< NekDouble > &SymmFlux, const Array< OneD, Array< OneD, NekDouble >> &pFwd, const Array< OneD, Array< OneD, NekDouble >> &pBwd, Array< OneD, int > &nonZeroIndex, Array< OneD, Array< OneD, NekDouble > > &solution_Aver, Array< OneD, Array< OneD, NekDouble > > &solution_jump)

virtual SOLVER_UTILS_EXPORT void	v_AddDiffusionSymmFluxToCoeff (const std::size_t nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, TensorOfArray3D< NekDouble > &qfield, TensorOfArray3D< NekDouble > &VolumeFlux, Array< OneD, Array< OneD, NekDouble > > &outarray, const Array< OneD, Array< OneD, NekDouble > > &pFwd, const Array< OneD, Array< OneD, NekDouble > > &pBwd)

virtual SOLVER_UTILS_EXPORT void	v_AddDiffusionSymmFluxToPhys (const std::size_t nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, TensorOfArray3D< NekDouble > &qfield, TensorOfArray3D< NekDouble > &VolumeFlux, Array< OneD, Array< OneD, NekDouble > > &outarray, const Array< OneD, Array< OneD, NekDouble > > &pFwd, const Array< OneD, Array< OneD, NekDouble > > &pBwd)

virtual void	v_SetHomoDerivs (Array< OneD, Array< OneD, NekDouble > > &deriv)

virtual TensorOfArray3D< NekDouble > &	v_GetFluxTensor ()

virtual SOLVER_UTILS_EXPORT const Array< OneD, const Array< OneD, NekDouble > > &	v_GetTraceNormal ()

virtual SOLVER_UTILS_EXPORT void	v_GetPrimVar (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &primVar)
	Compute primary derivatives. More...

SOLVER_UTILS_EXPORT void	GetDivCurl (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const TensorOfArray3D< NekDouble > &pVarDer)
	Compute divergence and curl squared. More...

Detailed Description

Definition at line 155 of file Diffusion.h.

Constructor & Destructor Documentation

◆ ~Diffusion()

Diffusion::~Diffusion ( )

inlinevirtual

Definition at line 163 of file Diffusion.h.

164 {};

Member Function Documentation

◆ AddDiffusionSymmFluxToCoeff()

SOLVER_UTILS_EXPORT void Nektar::SolverUtils::Diffusion::AddDiffusionSymmFluxToCoeff	(	const std::size_t	nConvectiveFields,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		TensorOfArray3D< NekDouble > &	qfield,
		TensorOfArray3D< NekDouble > &	VolumeFlux,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		const Array< OneD, Array< OneD, NekDouble > > &	pFwd,
		const Array< OneD, Array< OneD, NekDouble > > &	pBwd
	)

inline

Add symmetric flux to field in coeff space.

Definition at line 285 of file Diffusion.h.

             {
                 v_AddDiffusionSymmFluxToCoeff(nConvectiveFields, fields,
                     inarray, qfield, VolumeFlux, outarray, pFwd, pBwd);
             }

References v_AddDiffusionSymmFluxToCoeff().

Referenced by Nektar::SolverUtils::DiffusionIP::v_DiffuseCoeffs().

◆ AddDiffusionSymmFluxToPhys()

SOLVER_UTILS_EXPORT void Nektar::SolverUtils::Diffusion::AddDiffusionSymmFluxToPhys	(	const std::size_t	nConvectiveFields,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		TensorOfArray3D< NekDouble > &	qfield,
		TensorOfArray3D< NekDouble > &	VolumeFlux,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		const Array< OneD, Array< OneD, NekDouble > > &	pFwd,
		const Array< OneD, Array< OneD, NekDouble > > &	pBwd
	)

inline

Add symmetric flux to field in coeff physical space.

Definition at line 300 of file Diffusion.h.

             {
                 v_AddDiffusionSymmFluxToPhys(nConvectiveFields, fields, inarray,
                     qfield, VolumeFlux, outarray, pFwd, pBwd);
             }

References v_AddDiffusionSymmFluxToPhys().

◆ AddSymmFluxIntegralToOffDiag()

void Diffusion::AddSymmFluxIntegralToOffDiag	(	const int	nConvectiveFields,
		const int	nDim,
		const int	nPts,
		const int	nTracePts,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, const int > &	nonZeroIndex,
		TensorOfArray3D< NekDouble > &	Fwdflux,
		TensorOfArray3D< NekDouble > &	Bwdflux,
		Array< OneD, Array< OneD, NekDouble > > &	outarray
	)

Definition at line 137 of file Diffusion.cpp.

         {
             boost::ignore_unused(nTracePts);
             
             int nCoeffs =   outarray[nConvectiveFields-1].size();
             Array<OneD, NekDouble> tmpCoeff(nCoeffs,0.0);
             Array<OneD, Array<OneD, NekDouble> > tmpfield(nDim);
  
             for(int i = 0;i<nDim;i++)
             {
                 tmpfield[i]    =   Array<OneD, NekDouble>(nPts,0.0);
             }
             int nv = 0;
  
             for(int j=0;j<nonZeroIndex.size();j++)
             {
                 nv  =   nonZeroIndex[j];
                 MultiRegions::ExpListSharedPtr tracelist = fields[nv]->GetTrace();
                 for(int nd=0;nd<nDim;nd++)
                 {
                     Vmath::Zero(nPts,tmpfield[nd],1);
  
                     tracelist->MultiplyByQuadratureMetric(Fwdflux[nd][nv],
                                                           Fwdflux[nd][nv]);
                     tracelist->MultiplyByQuadratureMetric(Bwdflux[nd][nv],
                                                           Bwdflux[nd][nv]);
  
                     fields[nv]->AddTraceQuadPhysToOffDiag(Fwdflux[nd][nv],
                                                 Bwdflux[nd][nv],tmpfield[nd]);
                     fields[nv]->DivideByQuadratureMetric(tmpfield[nd],tmpfield[nd]);
                 }
                 fields[nv]->IProductWRTDerivBase(tmpfield,tmpCoeff);
                 Vmath::Vadd(nCoeffs,tmpCoeff,1,outarray[nv],1,outarray[nv],1);
             }
         }

References Vmath::Vadd(), and Vmath::Zero().

◆ ConsVarAveJump()

SOLVER_UTILS_EXPORT void Nektar::SolverUtils::Diffusion::ConsVarAveJump	(	const std::size_t	nConvectiveFields,
		const size_t	npnts,
		const Array< OneD, const Array< OneD, NekDouble > > &	vFwd,
		const Array< OneD, const Array< OneD, NekDouble > > &	vBwd,
		Array< OneD, Array< OneD, NekDouble > > &	aver,
		Array< OneD, Array< OneD, NekDouble > > &	jump
	)

inline

Get the average and jump value of conservative variables on trace.

Definition at line 446 of file Diffusion.h.

             {
                 v_ConsVarAveJump(nConvectiveFields, npnts, vFwd, vBwd, aver,
                     jump);
             }

References v_ConsVarAveJump().

Referenced by Nektar::SolverUtils::DiffusionIP::CalcTraceNumFlux().

◆ Diffuse() [1/2]

void Diffusion::Diffuse	(	const std::size_t	nConvectiveFields,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		const Array< OneD, Array< OneD, NekDouble > > &	pFwd = `NullNekDoubleArrayOfArray`,
		const Array< OneD, Array< OneD, NekDouble > > &	pBwd = `NullNekDoubleArrayOfArray`
	)

Definition at line 60 of file Diffusion.cpp.

         {
             v_Diffuse(nConvectiveFields, fields, inarray, outarray, pFwd, pBwd);
         }

References v_Diffuse().

◆ Diffuse() [2/2]

void Diffusion::Diffuse	(	const std::size_t	nConvectiveFields,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		NekDouble	time,
		const Array< OneD, Array< OneD, NekDouble > > &	pFwd = `NullNekDoubleArrayOfArray`,
		const Array< OneD, Array< OneD, NekDouble > > &	pBwd = `NullNekDoubleArrayOfArray`
	)

Definition at line 89 of file Diffusion.cpp.

         {
             m_time  =    time;
             v_Diffuse(nConvectiveFields, fields, inarray, outarray, pFwd, pBwd);
         }

References m_time, and v_Diffuse().

◆ DiffuseCalcDerivative()

void Diffusion::DiffuseCalcDerivative	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble >> &	inarray,
		TensorOfArray3D< NekDouble > &	qfields,
		const Array< OneD, Array< OneD, NekDouble >> &	pFwd = `NullNekDoubleArrayOfArray`,
		const Array< OneD, Array< OneD, NekDouble >> &	pBwd = `NullNekDoubleArrayOfArray`
	)

Definition at line 224 of file Diffusion.cpp.

         {
             v_DiffuseCalcDerivative(fields, inarray, qfields, pFwd, pBwd);
         }

References v_DiffuseCalcDerivative().

Referenced by Nektar::SolverUtils::DiffusionLDG::v_DiffuseCoeffs(), and Nektar::DiffusionLDGNS::v_DiffuseCoeffs().

◆ DiffuseCoeffs() [1/3]

SOLVER_UTILS_EXPORT void Nektar::SolverUtils::Diffusion::DiffuseCoeffs	(	const std::size_t	nConvectiveFields,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		const Array< OneD, Array< OneD, NekDouble > > &	pFwd,
		const Array< OneD, Array< OneD, NekDouble > > &	pBwd,
		TensorOfArray3D< NekDouble > &	qfield,
		Array< OneD, int > &	nonZeroIndex
	)

inline

Definition at line 199 of file Diffusion.h.

             {
                 v_DiffuseCoeffs(nConvectiveFields, fields, inarray, outarray,
                                 pFwd, pBwd,qfield,nonZeroIndex);
             }

References v_DiffuseCoeffs().

◆ DiffuseCoeffs() [2/3]

void Diffusion::DiffuseCoeffs	(	const std::size_t	nConvectiveFields,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		const Array< OneD, Array< OneD, NekDouble > > &	pFwd = `NullNekDoubleArrayOfArray`,
		const Array< OneD, Array< OneD, NekDouble > > &	pBwd = `NullNekDoubleArrayOfArray`
	)

Similar with Diffusion::Diffuse(): calculate diffusion flux The difference is in the outarray: it is the coefficients of basis for DiffuseCoeffs() it is the physics on quadrature points for Diffuse()

Definition at line 77 of file Diffusion.cpp.

         {
             v_DiffuseCoeffs(nConvectiveFields, fields, inarray,
                             outarray, pFwd, pBwd);
         }

References v_DiffuseCoeffs().

◆ DiffuseCoeffs() [3/3]

void Diffusion::DiffuseCoeffs	(	const std::size_t	nConvectiveFields,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		NekDouble	time,
		const Array< OneD, Array< OneD, NekDouble > > &	pFwd = `NullNekDoubleArrayOfArray`,
		const Array< OneD, Array< OneD, NekDouble > > &	pBwd = `NullNekDoubleArrayOfArray`
	)

Definition at line 102 of file Diffusion.cpp.

         {
             m_time  =    time;
             v_DiffuseCoeffs(nConvectiveFields, fields, inarray, outarray,
                             pFwd, pBwd);
         }

References m_time, and v_DiffuseCoeffs().

◆ DiffuseTraceFlux()

SOLVER_UTILS_EXPORT void Nektar::SolverUtils::Diffusion::DiffuseTraceFlux	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble >> &	inarray,
		TensorOfArray3D< NekDouble > &	qfields,
		TensorOfArray3D< NekDouble > &	VolumeFlux,
		Array< OneD, Array< OneD, NekDouble > > &	TraceFlux,
		const Array< OneD, Array< OneD, NekDouble >> &	pFwd = `NullNekDoubleArrayOfArray`,
		const Array< OneD, Array< OneD, NekDouble >> &	pBwd = `NullNekDoubleArrayOfArray`,
		Array< OneD, int > &	nonZeroIndex = `NullInt1DArray`
	)

inline

Diffusion term Trace Flux.

Definition at line 249 of file Diffusion.h.

             {
                 v_DiffuseTraceFlux(fields, inarray, qfields,
                     VolumeFlux, TraceFlux, pFwd, pBwd, nonZeroIndex);
             }

References v_DiffuseTraceFlux().

Referenced by Nektar::SolverUtils::DiffusionLDG::v_DiffuseCoeffs(), Nektar::DiffusionLDGNS::v_DiffuseCoeffs(), and Nektar::SolverUtils::DiffusionIP::v_DiffuseCoeffs().

◆ DiffuseTraceSymmFlux()

SOLVER_UTILS_EXPORT void Nektar::SolverUtils::Diffusion::DiffuseTraceSymmFlux	(	const int	nConvectiveFields,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble >> &	inarray,
		const TensorOfArray3D< NekDouble > &	qfield,
		const TensorOfArray3D< NekDouble > &	VolumeFlux,
		TensorOfArray3D< NekDouble > &	SymmFlux,
		const Array< OneD, Array< OneD, NekDouble >> &	pFwd,
		const Array< OneD, Array< OneD, NekDouble >> &	pBwd,
		Array< OneD, int > &	nonZeroIndex,
		Array< OneD, Array< OneD, NekDouble > > &	solution_Aver = `NullNekDoubleArrayOfArray`,
		Array< OneD, Array< OneD, NekDouble > > &	solution_jump = `NullNekDoubleArrayOfArray`
	)

inline

Definition at line 266 of file Diffusion.h.

             {
                 v_DiffuseTraceSymmFlux(nConvectiveFields,fields,inarray,qfield,
                                        VolumeFlux,SymmFlux,pFwd,pBwd,
                                        nonZeroIndex,solution_Aver,solution_jump);
             }

References v_DiffuseTraceSymmFlux().

◆ DiffuseVolumeFlux()

SOLVER_UTILS_EXPORT void Nektar::SolverUtils::Diffusion::DiffuseVolumeFlux	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble >> &	inarray,
		TensorOfArray3D< NekDouble > &	qfields,
		TensorOfArray3D< NekDouble > &	VolumeFlux,
		Array< OneD, int > &	nonZeroIndex = `NullInt1DArray`
	)

inline

Diffusion Volume FLux.

Definition at line 236 of file Diffusion.h.

             {
                 v_DiffuseVolumeFlux(fields, inarray, qfields,
                                     VolumeFlux, nonZeroIndex);
             }

References v_DiffuseVolumeFlux().

Referenced by Nektar::SolverUtils::DiffusionLDG::v_DiffuseCoeffs(), Nektar::DiffusionLDGNS::v_DiffuseCoeffs(), and Nektar::SolverUtils::DiffusionIP::v_DiffuseCoeffs().

◆ FluxVec()

SOLVER_UTILS_EXPORT void Nektar::SolverUtils::Diffusion::FluxVec ( TensorOfArray3D< NekDouble > & fluxvector )

◆ GetAVmu()

void Diffusion::GetAVmu	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		Array< OneD, NekDouble > &	muvar,
		Array< OneD, NekDouble > &	MuVarTrace
	)

Get the mu of artifical viscosity(AV)

Definition at line 115 of file Diffusion.cpp.

         {
             size_t nTracePts = fields[0]->GetTrace()->GetTotPoints();
  
             Array<OneD, NekDouble> Fwd{nTracePts, 0.0};
             Array<OneD, NekDouble> Bwd{nTracePts, 0.0};
  
             m_ArtificialDiffusionVector(inarray, muvar);
  
             // BwdMuvar is left to be 0.0 according to DiffusionLDG.cpp
             fields[0]->GetFwdBwdTracePhys(muvar, Fwd, Bwd, false);
  
             for (int k = 0; k < nTracePts; ++k)
             {
                 MuVarTrace[k] = 0.5 * (Fwd[k] + Bwd[k]) ;
             }
         }

References m_ArtificialDiffusionVector.

Referenced by Nektar::SolverUtils::DiffusionIP::v_DiffuseVolumeFlux().

◆ GetDivCurl()

void Diffusion::GetDivCurl	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const TensorOfArray3D< NekDouble > &	pVarDer
	)

protected

Compute divergence and curl squared.

Compute and store scalars needed for shock sensor, i.e. divergence and curl squared.

Parameters

dimensions,points,primary variables derivative

Definition at line 319 of file Diffusion.cpp.

         {
             int nDim = fields[0]->GetCoordim(0);
             int nPts = fields[0]->GetTotPoints();
             // Compute and store scalars needed for shock sensor
             Array<OneD, NekDouble>              tmp3(nPts, 0.0);
             Array<OneD, NekDouble>              tmp4(nPts, 0.0);
  
             // Zero the variables for the current iteration step
             Vmath::Zero(nPts, m_divVel,        1);
             Vmath::Zero(nPts, m_divVelSquare,  1);
             Vmath::Zero(nPts, m_curlVelSquare, 1);
             
             // div vel
             for (int j = 0; j < nDim; ++j)
             {
                 Vmath::Vadd(nPts, m_divVel, 1, pVarDer[j][j], 1, m_divVel, 1);
             }
             // (div vel)**2
             Vmath::Vmul(nPts, m_divVel, 1, m_divVel, 1, m_divVelSquare, 1);
  
             // curl vel
             if ( nDim > 2 )
             {
                 // curl[0] 3/2
                 Vmath::Vadd(nPts, tmp4, 1, pVarDer[2][1], 1, tmp4, 1);
                 // curl[0]-2/3
                 Vmath::Vcopy(nPts, pVarDer[1][2], 1, tmp3, 1);
                 Vmath::Neg(nPts, tmp3, 1);
                 Vmath::Vadd(nPts, tmp4, 1, tmp3, 1, tmp4, 1);
                 // square curl[0]
                 Vmath::Vmul(nPts, tmp4, 1, tmp4, 1, tmp3, 1);
                 Vmath::Vadd(nPts, m_curlVelSquare, 1, tmp3, 1, m_curlVelSquare, 1);
  
                 tmp4 = Array<OneD, NekDouble>(nPts, 0.0);
                 // curl[1] 3/1
                 Vmath::Vadd(nPts, tmp4, 1, pVarDer[2][0], 1, tmp4, 1);
                 // curl[1]-1/3
                 Vmath::Vcopy(nPts, pVarDer[0][2], 1, tmp3, 1);
                 Vmath::Neg(nPts, tmp3, 1);
                 Vmath::Vadd(nPts, tmp4, 1, tmp3, 1, tmp4, 1);
                 // square curl[1]
                 Vmath::Vmul(nPts, tmp4, 1, tmp4, 1, tmp3, 1);
                 Vmath::Vadd(nPts, m_curlVelSquare, 1, tmp3, 1, m_curlVelSquare, 1);
  
                 tmp4 = Array<OneD, NekDouble>(nPts, 0.0);
                 // curl[2] 1/2
                 Vmath::Vadd(nPts, tmp4, 1, pVarDer[0][1], 1, tmp4, 1);
                 // curl[2]-2/1
                 Vmath::Vcopy(nPts, pVarDer[1][0], 1, tmp3, 1);
                 Vmath::Neg(nPts, tmp3, 1);
                 Vmath::Vadd(nPts, tmp4, 1, tmp3, 1, tmp4, 1);
                 // square curl[2]
                 Vmath::Vmul(nPts, tmp4, 1, tmp4, 1, tmp3, 1);
                 Vmath::Vadd(nPts, m_curlVelSquare, 1, tmp3, 1, m_curlVelSquare, 1);
  
             }
             else if ( nDim > 1 )
             {
                 // curl[2] 1/2
                 Vmath::Vadd(nPts, tmp4, 1, pVarDer[0][1], 1, tmp4, 1);
                 // curl[2]-2/1
                 Vmath::Vcopy(nPts, pVarDer[1][0], 1, tmp3, 1);
                 Vmath::Neg(nPts, tmp3, 1);
                 Vmath::Vadd(nPts, tmp4, 1, tmp3, 1, tmp4, 1);
                 // square curl[2]
                 Vmath::Vmul(nPts, tmp4, 1, tmp4, 1, tmp3, 1);
                 Vmath::Vadd(nPts, m_curlVelSquare, 1, tmp3, 1, m_curlVelSquare, 1);
             }
             else
             {
                 Vmath::Fill(nPts, 0.0, m_curlVelSquare, 1);
             }
         }

References Vmath::Fill(), m_curlVelSquare, m_divVel, m_divVelSquare, Vmath::Neg(), Vmath::Vadd(), Vmath::Vcopy(), Vmath::Vmul(), and Vmath::Zero().

◆ GetFluxTensor()

virtual TensorOfArray3D<NekDouble>& Nektar::SolverUtils::Diffusion::GetFluxTensor ( )

inlinevirtual

Definition at line 434 of file Diffusion.h.

             {
                 return v_GetFluxTensor();
             }

References v_GetFluxTensor().

◆ GetTraceNormal()

SOLVER_UTILS_EXPORT const Array<OneD, const Array<OneD, NekDouble> >& Nektar::SolverUtils::Diffusion::GetTraceNormal ( )

inline

Get trace normal.

Definition at line 460 of file Diffusion.h.

             {
                 return v_GetTraceNormal();
             }

References v_GetTraceNormal().

◆ InitObject()

void Diffusion::InitObject	(	LibUtilities::SessionReaderSharedPtr	pSession,
		Array< OneD, MultiRegions::ExpListSharedPtr >	pFields
	)

Definition at line 47 of file Diffusion.cpp.

         {
             v_InitObject(pSession, pFields);
  
             // Div curl storage
             int nPts = pFields[0]->GetTotPoints();
             m_divVel        = Array<OneD, NekDouble>(nPts, 0.0);
             m_divVelSquare  = Array<OneD, NekDouble>(nPts, 0.0);
             m_curlVelSquare = Array<OneD, NekDouble>(nPts, 0.0);
         }

References m_curlVelSquare, m_divVel, m_divVelSquare, and v_InitObject().

◆ SetArtificialDiffusionVector()

template<typename FuncPointerT , typename ObjectPointerT >

void Nektar::SolverUtils::Diffusion::SetArtificialDiffusionVector	(	FuncPointerT	func,
		ObjectPointerT	obj
	)

inline

Definition at line 399 of file Diffusion.h.

             {
                 m_ArtificialDiffusionVector = std::bind(
                     func, obj, std::placeholders::_1, std::placeholders::_2);
             }

References m_ArtificialDiffusionVector.

◆ SetCalcViscosity()

template<typename FuncPointerT , typename ObjectPointerT >

void Nektar::SolverUtils::Diffusion::SetCalcViscosity	(	FuncPointerT	func,
		ObjectPointerT	obj
	)

inline

Definition at line 414 of file Diffusion.h.

             {
                 m_CalcViscosity = std::bind(
                     func, obj, std::placeholders::_1, std::placeholders::_2);
             }

References m_CalcViscosity.

◆ SetDiffusionFluxCons() [1/2]

void Nektar::SolverUtils::Diffusion::SetDiffusionFluxCons ( DiffusionFluxCons flux )

inline

Definition at line 378 of file Diffusion.h.

             {
                 m_FunctorDiffusionfluxCons =   flux;
             }

References m_FunctorDiffusionfluxCons.

◆ SetDiffusionFluxCons() [2/2]

template<typename FuncPointerT , typename ObjectPointerT >

void Nektar::SolverUtils::Diffusion::SetDiffusionFluxCons	(	FuncPointerT	func,
		ObjectPointerT	obj
	)

inline

Definition at line 369 of file Diffusion.h.

             {
                 m_FunctorDiffusionfluxCons = std::bind(
                     func, obj, std::placeholders::_1, std::placeholders::_2,
                                std::placeholders::_3, std::placeholders::_4,
                                std::placeholders::_5, std::placeholders::_6,
                                std::placeholders::_7);
             }

References m_FunctorDiffusionfluxCons.

◆ SetDiffusionFluxConsTrace() [1/2]

void Nektar::SolverUtils::Diffusion::SetDiffusionFluxConsTrace ( DiffusionFluxCons flux )

inline

Definition at line 393 of file Diffusion.h.

             {
                 m_FunctorDiffusionfluxConsTrace =   flux;
             }

References m_FunctorDiffusionfluxConsTrace.

◆ SetDiffusionFluxConsTrace() [2/2]

template<typename FuncPointerT , typename ObjectPointerT >

void Nektar::SolverUtils::Diffusion::SetDiffusionFluxConsTrace	(	FuncPointerT	func,
		ObjectPointerT	obj
	)

inline

Definition at line 384 of file Diffusion.h.

             {
                 m_FunctorDiffusionfluxConsTrace = std::bind(
                     func, obj, std::placeholders::_1, std::placeholders::_2,
                                std::placeholders::_3, std::placeholders::_4,
                                std::placeholders::_5, std::placeholders::_6,
                                std::placeholders::_7);
             }

References m_FunctorDiffusionfluxConsTrace.

◆ SetDiffusionSymmFluxCons()

template<typename FuncPointerT , typename ObjectPointerT >

void Nektar::SolverUtils::Diffusion::SetDiffusionSymmFluxCons	(	FuncPointerT	func,
		ObjectPointerT	obj
	)

inline

Definition at line 421 of file Diffusion.h.

             {
                 m_FunctorSymmetricfluxCons = std::bind(
                     func, obj, std::placeholders::_1, std::placeholders::_2,
                                std::placeholders::_3, std::placeholders::_4,
                                std::placeholders::_5, std::placeholders::_6);
             }

References m_FunctorSymmetricfluxCons.

◆ SetFluxPenaltyNS() [1/2]

void Nektar::SolverUtils::Diffusion::SetFluxPenaltyNS ( DiffusionFluxPenaltyNS flux )

inline

Definition at line 363 of file Diffusion.h.

             {
                 m_fluxPenaltyNS = flux;
             }

References m_fluxPenaltyNS.

◆ SetFluxPenaltyNS() [2/2]

template<typename FuncPointerT , typename ObjectPointerT >

void Nektar::SolverUtils::Diffusion::SetFluxPenaltyNS	(	FuncPointerT	func,
		ObjectPointerT	obj
	)

inline

Definition at line 356 of file Diffusion.h.

             {
                 m_fluxPenaltyNS = std::bind(
                     func, obj, std::placeholders::_1, std::placeholders::_2,
                     std::placeholders::_3);
             }

References m_fluxPenaltyNS.

◆ SetFluxVector() [1/2]

SOLVER_UTILS_EXPORT void Nektar::SolverUtils::Diffusion::SetFluxVector ( DiffusionFluxVecCB fluxVector )

inline

Definition at line 336 of file Diffusion.h.

             {
                 m_fluxVector = fluxVector;
             }

References m_fluxVector.

◆ SetFluxVector() [2/2]

template<typename FuncPointerT , typename ObjectPointerT >

void Nektar::SolverUtils::Diffusion::SetFluxVector	(	FuncPointerT	func,
		ObjectPointerT	obj
	)

inline

Definition at line 329 of file Diffusion.h.

             {
                 m_fluxVector = std::bind(
                     func, obj, std::placeholders::_1, std::placeholders::_2,
                     std::placeholders::_3);
             }

References m_fluxVector.

◆ SetFluxVectorNS() [1/2]

void Nektar::SolverUtils::Diffusion::SetFluxVectorNS ( DiffusionFluxVecCBNS fluxVector )

inline

Definition at line 350 of file Diffusion.h.

             {
                 m_fluxVectorNS = fluxVector;
             }

References m_fluxVectorNS.

◆ SetFluxVectorNS() [2/2]

template<typename FuncPointerT , typename ObjectPointerT >

void Nektar::SolverUtils::Diffusion::SetFluxVectorNS	(	FuncPointerT	func,
		ObjectPointerT	obj
	)

inline

Definition at line 343 of file Diffusion.h.

             {
                 m_fluxVectorNS = std::bind(
                     func, obj, std::placeholders::_1, std::placeholders::_2,
                     std::placeholders::_3);
             }

References m_fluxVectorNS.

◆ SetHomoDerivs()

void Nektar::SolverUtils::Diffusion::SetHomoDerivs ( Array< OneD, Array< OneD, NekDouble > > & deriv )

inline

Definition at line 429 of file Diffusion.h.

             {
                 v_SetHomoDerivs(deriv);
             }

References v_SetHomoDerivs().

◆ SetSpecialBndTreat()

template<typename FuncPointerT , typename ObjectPointerT >

void Nektar::SolverUtils::Diffusion::SetSpecialBndTreat	(	FuncPointerT	func,
		ObjectPointerT	obj
	)

inline

Definition at line 407 of file Diffusion.h.

             {
                 m_SpecialBndTreat = std::bind(
                     func, obj, std::placeholders::_1);
             }

References m_SpecialBndTreat.

◆ v_AddDiffusionSymmFluxToCoeff()

virtual SOLVER_UTILS_EXPORT void Nektar::SolverUtils::Diffusion::v_AddDiffusionSymmFluxToCoeff	(	const std::size_t	nConvectiveFields,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		TensorOfArray3D< NekDouble > &	qfield,
		TensorOfArray3D< NekDouble > &	VolumeFlux,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		const Array< OneD, Array< OneD, NekDouble > > &	pFwd,
		const Array< OneD, Array< OneD, NekDouble > > &	pBwd
	)

inlineprotectedvirtual

Definition at line 568 of file Diffusion.h.

             {
                 boost::ignore_unused(nConvectiveFields, fields, inarray, qfield,
                                         VolumeFlux, outarray, pFwd, pBwd);
  
             }

Referenced by AddDiffusionSymmFluxToCoeff().

◆ v_AddDiffusionSymmFluxToPhys()

virtual SOLVER_UTILS_EXPORT void Nektar::SolverUtils::Diffusion::v_AddDiffusionSymmFluxToPhys	(	const std::size_t	nConvectiveFields,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		TensorOfArray3D< NekDouble > &	qfield,
		TensorOfArray3D< NekDouble > &	VolumeFlux,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		const Array< OneD, Array< OneD, NekDouble > > &	pFwd,
		const Array< OneD, Array< OneD, NekDouble > > &	pBwd
	)

inlineprotectedvirtual

Definition at line 582 of file Diffusion.h.

             {
                 boost::ignore_unused(nConvectiveFields, fields, inarray, qfield,
                                         VolumeFlux, outarray, pFwd, pBwd);
  
             }

Referenced by AddDiffusionSymmFluxToPhys().

◆ v_ConsVarAveJump()

void Diffusion::v_ConsVarAveJump	(	const std::size_t	nConvectiveFields,
		const size_t	npnts,
		const Array< OneD, const Array< OneD, NekDouble > > &	vFwd,
		const Array< OneD, const Array< OneD, NekDouble > > &	vBwd,
		Array< OneD, Array< OneD, NekDouble > > &	aver,
		Array< OneD, Array< OneD, NekDouble > > &	jump
	)

protectedvirtual

Definition at line 242 of file Diffusion.cpp.

         {
             boost::ignore_unused(nConvectiveFields, npnts, vFwd, vBwd,
                                  aver, jump);
             NEKERROR(ErrorUtil::efatal, "v_ConsVarAveJump not defined");
         }

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by ConsVarAveJump().

◆ v_Diffuse()

void Diffusion::v_Diffuse	(	const std::size_t	nConvectiveFields,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		const Array< OneD, Array< OneD, NekDouble > > &	pFwd,
		const Array< OneD, Array< OneD, NekDouble > > &	pBwd
	)

protectedvirtual

Reimplemented in Nektar::SolverUtils::DiffusionLFR, Nektar::SolverUtils::Diffusion3DHomogeneous1D, Nektar::DiffusionLDGNS, Nektar::SolverUtils::DiffusionLFRNS, and Nektar::SolverUtils::DiffusionLDG.

Definition at line 182 of file Diffusion.cpp.

         {
             boost::ignore_unused(nConvectiveFields, fields, inarray, outarray,
                                     pFwd, pBwd);
             NEKERROR(ErrorUtil::efatal, "v_Diffuse not defined");
         }

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by Diffuse().

◆ v_DiffuseCalcDerivative()

void Diffusion::v_DiffuseCalcDerivative	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		TensorOfArray3D< NekDouble > &	qfields,
		const Array< OneD, Array< OneD, NekDouble > > &	pFwd,
		const Array< OneD, Array< OneD, NekDouble > > &	pBwd
	)

protectedvirtual

Diffusion Flux, calculate the physical derivatives.

Reimplemented in Nektar::DiffusionLDGNS, and Nektar::SolverUtils::DiffusionLDG.

Definition at line 255 of file Diffusion.cpp.

         {
             boost::ignore_unused(fields, inarray, qfields,
                                  pFwd, pBwd);
             NEKERROR(ErrorUtil::efatal, "Not defined for function DiffuseVolumeFLux.");
         }

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by DiffuseCalcDerivative().

◆ v_DiffuseCoeffs() [1/2]

void Diffusion::v_DiffuseCoeffs	(	const std::size_t	nConvectiveFields,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		const Array< OneD, Array< OneD, NekDouble > > &	pFwd,
		const Array< OneD, Array< OneD, NekDouble > > &	pBwd
	)

protectedvirtual

Reimplemented in Nektar::DiffusionLDGNS, and Nektar::SolverUtils::DiffusionLDG.

Definition at line 195 of file Diffusion.cpp.

         {
             boost::ignore_unused(nConvectiveFields, fields, inarray,
                                     outarray, pFwd, pBwd);
             NEKERROR(ErrorUtil::efatal, "v_DiffuseCoeffs not defined");
         }

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by DiffuseCoeffs().

◆ v_DiffuseCoeffs() [2/2]

void Diffusion::v_DiffuseCoeffs	(	const std::size_t	nConvectiveFields,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		const Array< OneD, Array< OneD, NekDouble > > &	vFwd,
		const Array< OneD, Array< OneD, NekDouble > > &	vBwd,
		TensorOfArray3D< NekDouble > &	qfield,
		Array< OneD, int > &	nonZeroIndex
	)

protectedvirtual

Definition at line 208 of file Diffusion.cpp.

         {
             boost::ignore_unused(nConvectiveFields, fields, inarray,
                                  outarray, vFwd, vBwd,qfield,nonZeroIndex);
             NEKERROR(ErrorUtil::efatal,"v_DiffuseCoeffs not defined");
         }

References Nektar::ErrorUtil::efatal, and NEKERROR.

◆ v_DiffuseTraceFlux()

void Diffusion::v_DiffuseTraceFlux	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble >> &	inarray,
		TensorOfArray3D< NekDouble > &	Qfields,
		TensorOfArray3D< NekDouble > &	VolumeFlux,
		Array< OneD, Array< OneD, NekDouble > > &	TraceFlux,
		const Array< OneD, Array< OneD, NekDouble >> &	pFwd,
		const Array< OneD, Array< OneD, NekDouble >> &	pBwd,
		Array< OneD, int > &	nonZeroIndex
	)

protectedvirtual

Diffusion term Trace Flux.

Reimplemented in Nektar::DiffusionLDGNS, and Nektar::SolverUtils::DiffusionLDG.

Definition at line 279 of file Diffusion.cpp.

         {
             boost::ignore_unused(fields, inarray, qfields, VolumeFlux,
                                  TraceFlux, pFwd, pBwd, nonZeroIndex);
             NEKERROR(ErrorUtil::efatal, "Not defined function DiffuseTraceFLux.");
         }

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by DiffuseTraceFlux().

◆ v_DiffuseTraceSymmFlux()

virtual SOLVER_UTILS_EXPORT void Nektar::SolverUtils::Diffusion::v_DiffuseTraceSymmFlux	(	const int	nConvectiveFields,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble >> &	inarray,
		const TensorOfArray3D< NekDouble > &	qfield,
		const TensorOfArray3D< NekDouble > &	VolumeFlux,
		TensorOfArray3D< NekDouble > &	SymmFlux,
		const Array< OneD, Array< OneD, NekDouble >> &	pFwd,
		const Array< OneD, Array< OneD, NekDouble >> &	pBwd,
		Array< OneD, int > &	nonZeroIndex,
		Array< OneD, Array< OneD, NekDouble > > &	solution_Aver,
		Array< OneD, Array< OneD, NekDouble > > &	solution_jump
	)

inlineprotectedvirtual

Definition at line 549 of file Diffusion.h.

             {
                 boost::ignore_unused(nConvectiveFields, fields, inarray, qfield,
                                      VolumeFlux, SymmFlux, pFwd, pBwd,
                                      nonZeroIndex,solution_Aver,solution_jump);
                 ASSERTL0(false, "Not defined function v_DiffuseTraceSymmFlux.");
             }

References ASSERTL0.

Referenced by DiffuseTraceSymmFlux().

◆ v_DiffuseVolumeFlux()

void Diffusion::v_DiffuseVolumeFlux	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble >> &	inarray,
		TensorOfArray3D< NekDouble > &	qfields,
		TensorOfArray3D< NekDouble > &	VolumeFlux,
		Array< OneD, int > &	nonZeroIndex
	)

protectedvirtual

Diffusion Volume Flux.

Reimplemented in Nektar::DiffusionLDGNS, Nektar::SolverUtils::DiffusionLDG, and Nektar::SolverUtils::DiffusionIP.

Definition at line 267 of file Diffusion.cpp.

         {
             boost::ignore_unused(fields, inarray, qfields, VolumeFlux,
                                  nonZeroIndex);
             NEKERROR(ErrorUtil::efatal, "Not defined for function DiffuseVolumeFLux.");
         }

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by DiffuseVolumeFlux().

◆ v_GetFluxTensor()

virtual TensorOfArray3D<NekDouble>& Nektar::SolverUtils::Diffusion::v_GetFluxTensor ( )

inlineprotectedvirtual

Reimplemented in Nektar::DiffusionLDGNS, and Nektar::SolverUtils::DiffusionLFRNS.

Definition at line 603 of file Diffusion.h.

             {
                 static TensorOfArray3D<NekDouble> tmp;
                 return tmp;
             }

Referenced by GetFluxTensor().

◆ v_GetPrimVar()

void Diffusion::v_GetPrimVar	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble >> &	inarray,
		Array< OneD, Array< OneD, NekDouble >> &	primVar
	)

protectedvirtual

Compute primary derivatives.

Compute primary variables.

Reimplemented in Nektar::DiffusionLDGNS.

Definition at line 298 of file Diffusion.cpp.

         {
  
             int nDim = fields[0]->GetCoordim(0);
             int nPts = fields[0]->GetTotPoints();
             for(int i = 0; i < nDim; ++i)
             {
                 primVar[i] = Array<OneD, NekDouble>(nPts, 0.0);
                 Vmath::Vdiv(nPts, inarray[i+1], 1, inarray[0], 1, primVar[i], 1);
             }
         }

References Vmath::Vdiv().

◆ v_GetTraceNormal()

const Array< OneD, const Array< OneD, NekDouble > > & Diffusion::v_GetTraceNormal ( )

protectedvirtual

Reimplemented in Nektar::SolverUtils::DiffusionIP.

Definition at line 236 of file Diffusion.cpp.

         {
             NEKERROR(ErrorUtil::efatal,"v_GetTraceNormal not defined");
             return NullNekDoubleArrayOfArray;
         }

References Nektar::ErrorUtil::efatal, NEKERROR, and Nektar::NullNekDoubleArrayOfArray.

Referenced by GetTraceNormal().

◆ v_InitObject()

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

inlineprotectedvirtual

Reimplemented in Nektar::DiffusionLDGNS, Nektar::SolverUtils::DiffusionLFRNS, Nektar::SolverUtils::DiffusionLFR, Nektar::SolverUtils::DiffusionLDG, Nektar::SolverUtils::DiffusionIP, and Nektar::SolverUtils::Diffusion3DHomogeneous1D.

Definition at line 478 of file Diffusion.h.

             {
                 boost::ignore_unused(pSession, pFields);
             };

Referenced by InitObject().

◆ v_SetHomoDerivs()

virtual void Nektar::SolverUtils::Diffusion::v_SetHomoDerivs ( Array< OneD, Array< OneD, NekDouble > > & deriv )

inlineprotectedvirtual

Reimplemented in Nektar::DiffusionLDGNS, and Nektar::SolverUtils::DiffusionLFRNS.

Definition at line 597 of file Diffusion.h.

             {
                 boost::ignore_unused(deriv);
             }

Referenced by SetHomoDerivs().

Member Data Documentation

◆ m_ArtificialDiffusionVector

DiffusionArtificialDiffusion Nektar::SolverUtils::Diffusion::m_ArtificialDiffusionVector

protected

Definition at line 469 of file Diffusion.h.

Referenced by GetAVmu(), SetArtificialDiffusionVector(), Nektar::SolverUtils::DiffusionIP::v_DiffuseVolumeFlux(), and Nektar::SolverUtils::DiffusionIP::v_InitObject().

◆ m_CalcViscosity

CalcViscosity Nektar::SolverUtils::Diffusion::m_CalcViscosity

protected

Definition at line 473 of file Diffusion.h.

Referenced by SetCalcViscosity().

◆ m_curlVelSquare

Array<OneD, NekDouble> Nektar::SolverUtils::Diffusion::m_curlVelSquare

Definition at line 161 of file Diffusion.h.

Referenced by GetDivCurl(), and InitObject().

◆ m_divVel

Array<OneD, NekDouble> Nektar::SolverUtils::Diffusion::m_divVel

Params for Ducros sensor.

Definition at line 159 of file Diffusion.h.

Referenced by GetDivCurl(), and InitObject().

◆ m_divVelSquare

Array<OneD, NekDouble> Nektar::SolverUtils::Diffusion::m_divVelSquare

Definition at line 160 of file Diffusion.h.

Referenced by GetDivCurl(), and InitObject().

◆ m_fluxPenaltyNS

DiffusionFluxPenaltyNS Nektar::SolverUtils::Diffusion::m_fluxPenaltyNS

protected

Definition at line 468 of file Diffusion.h.

Referenced by Nektar::DiffusionLDGNS::NumericalFluxO2(), and SetFluxPenaltyNS().

◆ m_fluxVector

DiffusionFluxVecCB Nektar::SolverUtils::Diffusion::m_fluxVector

protected

Definition at line 466 of file Diffusion.h.

Referenced by SetFluxVector(), and Nektar::SolverUtils::DiffusionLDG::v_DiffuseVolumeFlux().

◆ m_fluxVectorNS

DiffusionFluxVecCBNS Nektar::SolverUtils::Diffusion::m_fluxVectorNS

protected

Definition at line 467 of file Diffusion.h.

Referenced by SetFluxVectorNS(), Nektar::SolverUtils::DiffusionLFRNS::v_Diffuse(), and Nektar::DiffusionLDGNS::v_DiffuseVolumeFlux().

◆ m_FunctorDiffusionfluxCons

DiffusionFluxCons Nektar::SolverUtils::Diffusion::m_FunctorDiffusionfluxCons

protected

Definition at line 470 of file Diffusion.h.

Referenced by SetDiffusionFluxCons(), and Nektar::SolverUtils::DiffusionIP::v_DiffuseVolumeFlux().

◆ m_FunctorDiffusionfluxConsTrace

DiffusionFluxCons Nektar::SolverUtils::Diffusion::m_FunctorDiffusionfluxConsTrace

protected

Definition at line 471 of file Diffusion.h.

Referenced by Nektar::SolverUtils::DiffusionIP::CalcTraceNumFlux(), and SetDiffusionFluxConsTrace().

◆ m_FunctorSymmetricfluxCons

DiffusionSymmFluxCons Nektar::SolverUtils::Diffusion::m_FunctorSymmetricfluxCons

protected

Definition at line 474 of file Diffusion.h.

Referenced by Nektar::SolverUtils::DiffusionIP::CalcTraceSymFlux(), and SetDiffusionSymmFluxCons().

◆ m_SpecialBndTreat

SpecialBndTreat Nektar::SolverUtils::Diffusion::m_SpecialBndTreat

protected

Definition at line 472 of file Diffusion.h.

Referenced by SetSpecialBndTreat(), and Nektar::SolverUtils::DiffusionIP::v_ConsVarAveJump().

◆ m_time

NekDouble Nektar::SolverUtils::Diffusion::m_time =0.0

protected

Definition at line 476 of file Diffusion.h.

Referenced by Diffuse(), and DiffuseCoeffs().

Protected Attributes
DiffusionFluxVecCB	m_fluxVector

DiffusionFluxVecCBNS	m_fluxVectorNS

DiffusionFluxPenaltyNS	m_fluxPenaltyNS

DiffusionArtificialDiffusion	m_ArtificialDiffusionVector

DiffusionFluxCons	m_FunctorDiffusionfluxCons

DiffusionFluxCons	m_FunctorDiffusionfluxConsTrace

SpecialBndTreat	m_SpecialBndTreat

CalcViscosity	m_CalcViscosity

DiffusionSymmFluxCons	m_FunctorSymmetricfluxCons

NekDouble	m_time =0.0

Public Member Functions

Public Attributes

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ~Diffusion()

Member Function Documentation

◆ AddDiffusionSymmFluxToCoeff()

◆ AddDiffusionSymmFluxToPhys()

◆ AddSymmFluxIntegralToOffDiag()

◆ ConsVarAveJump()

◆ Diffuse() [1/2]

◆ Diffuse() [2/2]

◆ DiffuseCalcDerivative()

◆ DiffuseCoeffs() [1/3]

◆ DiffuseCoeffs() [2/3]

◆ DiffuseCoeffs() [3/3]

◆ DiffuseTraceFlux()

◆ DiffuseTraceSymmFlux()

◆ DiffuseVolumeFlux()

◆ FluxVec()

◆ GetAVmu()

◆ GetDivCurl()

◆ GetFluxTensor()

◆ GetTraceNormal()

◆ InitObject()

◆ SetArtificialDiffusionVector()

◆ SetCalcViscosity()

◆ SetDiffusionFluxCons() [1/2]

◆ SetDiffusionFluxCons() [2/2]

◆ SetDiffusionFluxConsTrace() [1/2]

◆ SetDiffusionFluxConsTrace() [2/2]

◆ SetDiffusionSymmFluxCons()

◆ SetFluxPenaltyNS() [1/2]

◆ SetFluxPenaltyNS() [2/2]

◆ SetFluxVector() [1/2]

◆ SetFluxVector() [2/2]

◆ SetFluxVectorNS() [1/2]

◆ SetFluxVectorNS() [2/2]

◆ SetHomoDerivs()

◆ SetSpecialBndTreat()

◆ v_AddDiffusionSymmFluxToCoeff()

◆ v_AddDiffusionSymmFluxToPhys()

◆ v_ConsVarAveJump()

◆ v_Diffuse()

◆ v_DiffuseCalcDerivative()

◆ v_DiffuseCoeffs() [1/2]

◆ v_DiffuseCoeffs() [2/2]

◆ v_DiffuseTraceFlux()

◆ v_DiffuseTraceSymmFlux()

◆ v_DiffuseVolumeFlux()

◆ v_GetFluxTensor()

◆ v_GetPrimVar()

◆ v_GetTraceNormal()

◆ v_InitObject()

◆ v_SetHomoDerivs()

Member Data Documentation

◆ m_ArtificialDiffusionVector

◆ m_CalcViscosity

◆ m_curlVelSquare

◆ m_divVel

◆ m_divVelSquare

◆ m_fluxPenaltyNS

◆ m_fluxVector

◆ m_fluxVectorNS

◆ m_FunctorDiffusionfluxCons

◆ m_FunctorDiffusionfluxConsTrace

◆ m_FunctorSymmetricfluxCons

◆ m_SpecialBndTreat

◆ m_time