#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)

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

Protected Attributes
DiffusionFluxVecCB	m_fluxVector

DiffusionFluxVecCBNS	m_fluxVectorNS

DiffusionFluxPenaltyNS	m_fluxPenaltyNS

DiffusionFluxCons	m_FunctorDiffusionfluxCons

DiffusionFluxCons	m_FunctorDiffusionfluxConsTrace

SpecialBndTreat	m_SpecialBndTreat

DiffusionSymmFluxCons	m_FunctorSymmetricfluxCons

NekDouble	m_time = 0.0

Detailed Description

Definition at line 127 of file Diffusion.h.

Constructor & Destructor Documentation

◆ ~Diffusion()

Diffusion::~Diffusion ( )

inlinevirtual

Definition at line 135 of file Diffusion.h.

135 {};

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 254 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 269 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 111 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 384 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 59 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 87 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 194 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 171 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 76 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 99 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 217 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 233 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 206 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().

◆ 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 281 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 378 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 396 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().

◆ SetDiffusionFluxCons() [1/2]

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

inline

Definition at line 339 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 331 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);
     }

References m_FunctorDiffusionfluxCons.

◆ SetDiffusionFluxConsTrace() [1/2]

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

inline

Definition at line 353 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 345 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);
     }

References m_FunctorDiffusionfluxConsTrace.

◆ SetDiffusionSymmFluxCons()

template<typename FuncPointerT , typename ObjectPointerT >

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

inline

Definition at line 365 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 325 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 318 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 299 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 293 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 312 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 305 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 373 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 359 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

Reimplemented in Nektar::SolverUtils::DiffusionIP.

Definition at line 496 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

Reimplemented in Nektar::SolverUtils::DiffusionIP.

Definition at line 509 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

Reimplemented in Nektar::SolverUtils::DiffusionIP.

Definition at line 210 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::DiffusionIP, Nektar::SolverUtils::Diffusion3DHomogeneous1D, Nektar::DiffusionLDGNS, Nektar::SolverUtils::DiffusionLFRNS, and Nektar::SolverUtils::DiffusionLDG.

Definition at line 153 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, Nektar::SolverUtils::DiffusionLDG, and Nektar::SolverUtils::DiffusionIP.

Definition at line 221 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::SolverUtils::DiffusionIP, Nektar::DiffusionLDGNS, and Nektar::SolverUtils::DiffusionLDG.

Definition at line 166 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

Reimplemented in Nektar::SolverUtils::DiffusionIP.

Definition at line 179 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, Nektar::SolverUtils::DiffusionLDG, and Nektar::SolverUtils::DiffusionIP.

Definition at line 242 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 477 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 232 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 528 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 260 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 204 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 412 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 523 of file Diffusion.h.

     {
         boost::ignore_unused(deriv);
     }

Referenced by SetHomoDerivs().

Member Data Documentation

◆ m_curlVelSquare

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

Definition at line 133 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 131 of file Diffusion.h.

Referenced by GetDivCurl(), and InitObject().

◆ m_divVelSquare

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

Definition at line 132 of file Diffusion.h.

Referenced by GetDivCurl(), and InitObject().

◆ m_fluxPenaltyNS

DiffusionFluxPenaltyNS Nektar::SolverUtils::Diffusion::m_fluxPenaltyNS

protected

Definition at line 404 of file Diffusion.h.

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

◆ m_fluxVector

DiffusionFluxVecCB Nektar::SolverUtils::Diffusion::m_fluxVector

protected

Definition at line 402 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 403 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 405 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 406 of file Diffusion.h.

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

◆ m_FunctorSymmetricfluxCons

DiffusionSymmFluxCons Nektar::SolverUtils::Diffusion::m_FunctorSymmetricfluxCons

protected

Definition at line 408 of file Diffusion.h.

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

◆ m_SpecialBndTreat

SpecialBndTreat Nektar::SolverUtils::Diffusion::m_SpecialBndTreat

protected

Definition at line 407 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 410 of file Diffusion.h.

Referenced by Diffuse(), and DiffuseCoeffs().

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()

◆ GetDivCurl()

◆ GetFluxTensor()

◆ GetTraceNormal()

◆ InitObject()

◆ 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_curlVelSquare

◆ m_divVel

◆ m_divVelSquare

◆ m_fluxPenaltyNS

◆ m_fluxVector

◆ m_fluxVectorNS

◆ m_FunctorDiffusionfluxCons

◆ m_FunctorDiffusionfluxConsTrace

◆ m_FunctorSymmetricfluxCons

◆ m_SpecialBndTreat

◆ m_time