Nektar++
Static Public Member Functions | Static Public Attributes | Protected Member Functions | Protected Attributes | List of all members
Nektar::SolverUtils::DiffusionLDG Class Reference

#include <DiffusionLDG.h>

Inheritance diagram for Nektar::SolverUtils::DiffusionLDG:
[legend]

Static Public Member Functions

static DiffusionSharedPtr create (std::string diffType)
 

Static Public Attributes

static std::string type
 

Protected Member Functions

 DiffusionLDG ()
 
virtual void v_InitObject (LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields) override
 
virtual void v_Diffuse (const std::size_t nConvective, 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) override
 
virtual void v_DiffuseCoeffs (const std::size_t nConvective, 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) override
 
virtual 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) override
 Diffusion Flux, calculate the physical derivatives. More...
 
virtual 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) override
 Diffusion Volume Flux. More...
 
virtual 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) override
 Diffusion term Trace Flux. More...
 
void NumFluxforScalar (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &ufield, TensorOfArray3D< NekDouble > &uflux, const Array< OneD, Array< OneD, NekDouble >> &pFwd, const Array< OneD, Array< OneD, NekDouble >> &pBwd)
 
void ApplyScalarBCs (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const std::size_t var, const Array< OneD, const NekDouble > &ufield, const Array< OneD, const NekDouble > &Fwd, const Array< OneD, const NekDouble > &Bwd, Array< OneD, NekDouble > &penaltyflux)
 
void NumFluxforVector (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &ufield, TensorOfArray3D< NekDouble > &qfield, Array< OneD, Array< OneD, NekDouble >> &qflux)
 Build the numerical flux for the 2nd order derivatives todo: add variable coeff and h dependence to penalty term. More...
 
void ApplyVectorBCs (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const std::size_t var, const std::size_t dir, const Array< OneD, const NekDouble > &qfield, const Array< OneD, const NekDouble > &qFwd, const Array< OneD, const NekDouble > &qBwd, Array< OneD, NekDouble > &penaltyflux)
 
- Protected Member Functions inherited from Nektar::SolverUtils::Diffusion
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_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

std::string m_shockCaptureType
 
NekDouble m_C11
 Coefficient of penalty term. More...
 
Array< OneD, Array< OneD, NekDouble > > m_traceNormals
 
LibUtilities::SessionReaderSharedPtr m_session
 
- Protected Attributes inherited from Nektar::SolverUtils::Diffusion
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
 

Additional Inherited Members

- Public Member Functions inherited from Nektar::SolverUtils::Diffusion
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 inherited from Nektar::SolverUtils::Diffusion
Array< OneD, NekDouble > m_divVel
 Params for Ducros sensor. More...
 
Array< OneD, NekDouble > m_divVelSquare
 
Array< OneD, NekDouble > m_curlVelSquare
 

Detailed Description

Definition at line 46 of file DiffusionLDG.h.

Constructor & Destructor Documentation

◆ DiffusionLDG()

Nektar::SolverUtils::DiffusionLDG::DiffusionLDG ( )
protected

Definition at line 49 of file DiffusionLDG.cpp.

50 {
51 }

Referenced by create().

Member Function Documentation

◆ ApplyScalarBCs()

void Nektar::SolverUtils::DiffusionLDG::ApplyScalarBCs ( const Array< OneD, MultiRegions::ExpListSharedPtr > &  fields,
const std::size_t  var,
const Array< OneD, const NekDouble > &  ufield,
const Array< OneD, const NekDouble > &  Fwd,
const Array< OneD, const NekDouble > &  Bwd,
Array< OneD, NekDouble > &  penaltyflux 
)
protected

Definition at line 275 of file DiffusionLDG.cpp.

281 {
282  boost::ignore_unused(ufield, Bwd);
283  // Number of boundary regions
284  std::size_t nBndRegions = fields[var]->GetBndCondExpansions().size();
285  std::size_t cnt = 0;
286  for (std::size_t i = 0; i < nBndRegions; ++i)
287  {
288  if (fields[var]->GetBndConditions()[i]->GetBoundaryConditionType() ==
289  SpatialDomains::ePeriodic)
290  {
291  continue;
292  }
293 
294  // Number of boundary expansion related to that region
295  std::size_t nBndEdges =
296  fields[var]->GetBndCondExpansions()[i]->GetExpSize();
297 
298  // Weakly impose boundary conditions by modifying flux values
299  for (std::size_t e = 0; e < nBndEdges; ++e)
300  {
301  std::size_t nBndEdgePts = fields[var]
302  ->GetBndCondExpansions()[i]
303  ->GetExp(e)
304  ->GetTotPoints();
305 
306  std::size_t id1 =
307  fields[var]->GetBndCondExpansions()[i]->GetPhys_Offset(e);
308 
309  std::size_t id2 = fields[0]->GetTrace()->GetPhys_Offset(
310  fields[0]->GetTraceMap()->GetBndCondIDToGlobalTraceID(cnt++));
311 
312  // AV boundary conditions
313  if (boost::iequals(
314  fields[var]->GetBndConditions()[i]->GetUserDefined(),
315  "Wall") ||
316  boost::iequals(
317  fields[var]->GetBndConditions()[i]->GetUserDefined(),
318  "Symmetry") ||
319  boost::iequals(
320  fields[var]->GetBndConditions()[i]->GetUserDefined(),
321  "WallViscous") ||
322  boost::iequals(
323  fields[var]->GetBndConditions()[i]->GetUserDefined(),
324  "WallAdiabatic"))
325  {
326  Vmath::Vcopy(nBndEdgePts, &Fwd[id2], 1, &penaltyflux[id2], 1);
327  }
328  // For Dirichlet boundary condition: uflux = g_D
329  else if (fields[var]
330  ->GetBndConditions()[i]
331  ->GetBoundaryConditionType() ==
332  SpatialDomains::eDirichlet)
333  {
334  Vmath::Vcopy(
335  nBndEdgePts,
336  &(fields[var]->GetBndCondExpansions()[i]->GetPhys())[id1],
337  1, &penaltyflux[id2], 1);
338  }
339  // For Neumann boundary condition: uflux = u+
340  else if ((fields[var]->GetBndConditions()[i])
341  ->GetBoundaryConditionType() ==
342  SpatialDomains::eNeumann)
343  {
344  Vmath::Vcopy(nBndEdgePts, &Fwd[id2], 1, &penaltyflux[id2], 1);
345  }
346  }
347  }
348 }
Vmath::Vcopy
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1255

References Nektar::SpatialDomains::eDirichlet, Nektar::SpatialDomains::eNeumann, Nektar::SpatialDomains::ePeriodic, and Vmath::Vcopy().

Referenced by NumFluxforScalar().

◆ ApplyVectorBCs()

void Nektar::SolverUtils::DiffusionLDG::ApplyVectorBCs ( const Array< OneD, MultiRegions::ExpListSharedPtr > &  fields,
const std::size_t  var,
const std::size_t  dir,
const Array< OneD, const NekDouble > &  qfield,
const Array< OneD, const NekDouble > &  qFwd,
const Array< OneD, const NekDouble > &  qBwd,
Array< OneD, NekDouble > &  penaltyflux 
)
protected

Diffusion: Imposing weak boundary condition for q with flux uflux = g_D on Dirichlet boundary condition uflux = u_Fwd on Neumann boundary condition

Definition at line 415 of file DiffusionLDG.cpp.

422 {
423  boost::ignore_unused(qfield, qBwd);
424 
425  std::size_t nBndRegions = fields[var]->GetBndCondExpansions().size();
426  std::size_t cnt = 0;
427 
428  for (std::size_t i = 0; i < nBndRegions; ++i)
429  {
430  if (fields[var]->GetBndConditions()[i]->GetBoundaryConditionType() ==
431  SpatialDomains::ePeriodic)
432  {
433  continue;
434  }
435  std::size_t nBndEdges =
436  fields[var]->GetBndCondExpansions()[i]->GetExpSize();
437 
438  // Weakly impose boundary conditions by modifying flux values
439  for (std::size_t e = 0; e < nBndEdges; ++e)
440  {
441  std::size_t nBndEdgePts = fields[var]
442  ->GetBndCondExpansions()[i]
443  ->GetExp(e)
444  ->GetTotPoints();
445 
446  std::size_t id1 =
447  fields[var]->GetBndCondExpansions()[i]->GetPhys_Offset(e);
448 
449  std::size_t id2 = fields[0]->GetTrace()->GetPhys_Offset(
450  fields[0]->GetTraceMap()->GetBndCondIDToGlobalTraceID(cnt++));
451 
452  // AV boundary conditions
453  if (boost::iequals(
454  fields[var]->GetBndConditions()[i]->GetUserDefined(),
455  "Wall") ||
456  boost::iequals(
457  fields[var]->GetBndConditions()[i]->GetUserDefined(),
458  "Symmetry") ||
459  boost::iequals(
460  fields[var]->GetBndConditions()[i]->GetUserDefined(),
461  "WallViscous") ||
462  boost::iequals(
463  fields[var]->GetBndConditions()[i]->GetUserDefined(),
464  "WallAdiabatic"))
465  {
466  Vmath::Zero(nBndEdgePts, &penaltyflux[id2], 1);
467  }
468  // For Dirichlet boundary condition:
469  // qflux = q+ - C_11 (u+ - g_D) (nx, ny)
470  else if (fields[var]
471  ->GetBndConditions()[i]
472  ->GetBoundaryConditionType() ==
473  SpatialDomains::eDirichlet)
474  {
475  Vmath::Vmul(nBndEdgePts, &m_traceNormals[dir][id2], 1,
476  &qFwd[id2], 1, &penaltyflux[id2], 1);
477  }
478  // For Neumann boundary condition: qflux = g_N
479  else if ((fields[var]->GetBndConditions()[i])
480  ->GetBoundaryConditionType() ==
481  SpatialDomains::eNeumann)
482  {
483  Vmath::Vmul(
484  nBndEdgePts, &m_traceNormals[dir][id2], 1,
485  &(fields[var]->GetBndCondExpansions()[i]->GetPhys())[id1],
486  1, &penaltyflux[id2], 1);
487  }
488  }
489  }
490 }
Nektar::SolverUtils::DiffusionLDG::m_traceNormals
Array< OneD, Array< OneD, NekDouble > > m_traceNormals
Definition: DiffusionLDG.h:65
Vmath::Vmul
void Vmul(int n, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
Multiply vector z = x*y.
Definition: Vmath.cpp:209
Vmath::Zero
void Zero(int n, T *x, const int incx)
Zero vector.
Definition: Vmath.cpp:492

References Nektar::SpatialDomains::eDirichlet, Nektar::SpatialDomains::eNeumann, Nektar::SpatialDomains::ePeriodic, m_traceNormals, Vmath::Vmul(), and Vmath::Zero().

Referenced by NumFluxforVector().

◆ create()

static DiffusionSharedPtr Nektar::SolverUtils::DiffusionLDG::create ( std::string  diffType)
inlinestatic

Definition at line 49 of file DiffusionLDG.h.

50  {
51  boost::ignore_unused(diffType);
52  return DiffusionSharedPtr(new DiffusionLDG());
53  }
Nektar::SolverUtils::DiffusionSharedPtr
std::shared_ptr< SolverUtils::Diffusion > DiffusionSharedPtr
A shared pointer to an EquationSystem object.
Definition: Diffusion.h:550

References DiffusionLDG().

◆ NumFluxforScalar()

void Nektar::SolverUtils::DiffusionLDG::NumFluxforScalar ( const Array< OneD, MultiRegions::ExpListSharedPtr > &  fields,
const Array< OneD, Array< OneD, NekDouble >> &  ufield,
TensorOfArray3D< NekDouble > &  uflux,
const Array< OneD, Array< OneD, NekDouble >> &  pFwd,
const Array< OneD, Array< OneD, NekDouble >> &  pBwd 
)
protected

Definition at line 226 of file DiffusionLDG.cpp.

232 {
233  std::size_t nTracePts = fields[0]->GetTrace()->GetTotPoints();
234  std::size_t nvariables = fields.size();
235  std::size_t nDim = fields[0]->GetCoordim(0);
236 
237  Array<OneD, NekDouble> Fwd{nTracePts};
238  Array<OneD, NekDouble> Bwd{nTracePts};
239  Array<OneD, NekDouble> fluxtemp{nTracePts, 0.0};
240 
241  // Get the sign of (v \cdot n), v = an arbitrary vector
242  // Evaluate upwind flux:
243  // uflux = \hat{u} \phi \cdot u = u^{(+,-)} n
244  for (std::size_t i = 0; i < nvariables; ++i)
245  {
246  // Compute Fwd and Bwd value of ufield of i direction
247  if (pFwd == NullNekDoubleArrayOfArray ||
248  pBwd == NullNekDoubleArrayOfArray)
249  {
250  fields[i]->GetFwdBwdTracePhys(ufield[i], Fwd, Bwd);
251  }
252  else
253  {
254  Fwd = pFwd[i];
255  Bwd = pBwd[i];
256  }
257 
258  // Upwind
259  Vmath::Vcopy(nTracePts, Fwd, 1, fluxtemp, 1);
260 
261  // Imposing weak boundary condition with flux
262  if (fields[0]->GetBndCondExpansions().size())
263  {
264  ApplyScalarBCs(fields, i, ufield[i], Fwd, Bwd, fluxtemp);
265  }
266 
267  for (std::size_t j = 0; j < nDim; ++j)
268  {
269  Vmath::Vmul(nTracePts, m_traceNormals[j], 1, fluxtemp, 1,
270  uflux[j][i], 1);
271  }
272  }
273 }
Nektar::SolverUtils::DiffusionLDG::ApplyScalarBCs
void ApplyScalarBCs(const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const std::size_t var, const Array< OneD, const NekDouble > &ufield, const Array< OneD, const NekDouble > &Fwd, const Array< OneD, const NekDouble > &Bwd, Array< OneD, NekDouble > &penaltyflux)
Definition: DiffusionLDG.cpp:275
Nektar::NullNekDoubleArrayOfArray
static Array< OneD, Array< OneD, NekDouble > > NullNekDoubleArrayOfArray
Definition: SharedArray.hpp:901

References ApplyScalarBCs(), m_traceNormals, Nektar::NullNekDoubleArrayOfArray, Vmath::Vcopy(), and Vmath::Vmul().

Referenced by v_DiffuseCalcDerivative().

◆ NumFluxforVector()

void Nektar::SolverUtils::DiffusionLDG::NumFluxforVector ( const Array< OneD, MultiRegions::ExpListSharedPtr > &  fields,
const Array< OneD, Array< OneD, NekDouble >> &  ufield,
TensorOfArray3D< NekDouble > &  qfield,
Array< OneD, Array< OneD, NekDouble >> &  qflux 
)
protected

Build the numerical flux for the 2nd order derivatives todo: add variable coeff and h dependence to penalty term.

Definition at line 354 of file DiffusionLDG.cpp.

359 {
360  std::size_t nTracePts = fields[0]->GetTrace()->GetTotPoints();
361  std::size_t nvariables = fields.size();
362  std::size_t nDim = qfield.size();
363 
364  Array<OneD, NekDouble> Fwd{nTracePts};
365  Array<OneD, NekDouble> Bwd{nTracePts};
366  Array<OneD, NekDouble> qFwd{nTracePts};
367  Array<OneD, NekDouble> qBwd{nTracePts};
368  Array<OneD, NekDouble> qfluxtemp{nTracePts, 0.0};
369  Array<OneD, NekDouble> uterm{nTracePts};
370 
371  // Evaulate upwind flux:
372  // qflux = \hat{q} \cdot u = q \cdot n - C_(11)*(u^+ - u^-)
373  for (std::size_t i = 0; i < nvariables; ++i)
374  {
375  // Generate Stability term = - C11 ( u- - u+ )
376  fields[i]->GetFwdBwdTracePhys(ufield[i], Fwd, Bwd);
377  Vmath::Vsub(nTracePts, Fwd, 1, Bwd, 1, uterm, 1);
378  Vmath::Smul(nTracePts, -m_C11, uterm, 1, uterm, 1);
379 
380  qflux[i] = Array<OneD, NekDouble>{nTracePts, 0.0};
381  for (std::size_t j = 0; j < nDim; ++j)
382  {
383  // Compute Fwd and Bwd value of ufield of jth direction
384  fields[i]->GetFwdBwdTracePhys(qfield[j][i], qFwd, qBwd);
385 
386  // Downwind
387  Vmath::Vcopy(nTracePts, qBwd, 1, qfluxtemp, 1);
388 
389  Vmath::Vmul(nTracePts, m_traceNormals[j], 1, qfluxtemp, 1,
390  qfluxtemp, 1);
391 
392  // Flux = {Fwd, Bwd} * (nx, ny, nz) + uterm * (nx, ny)
393  Vmath::Vadd(nTracePts, uterm, 1, qfluxtemp, 1, qfluxtemp, 1);
394 
395  // Imposing weak boundary condition with flux
396  if (fields[0]->GetBndCondExpansions().size())
397  {
398  ApplyVectorBCs(fields, i, j, qfield[j][i], qFwd, qBwd,
399  qfluxtemp);
400  }
401 
402  // q_hat \cdot n = (q_xi \cdot n_xi) or (q_eta \cdot n_eta)
403  // n_xi = n_x * tan_xi_x + n_y * tan_xi_y + n_z * tan_xi_z
404  // n_xi = n_x * tan_eta_x + n_y * tan_eta_y + n_z*tan_eta_z
405  Vmath::Vadd(nTracePts, qfluxtemp, 1, qflux[i], 1, qflux[i], 1);
406  }
407  }
408 }
Nektar::SolverUtils::DiffusionLDG::ApplyVectorBCs
void ApplyVectorBCs(const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const std::size_t var, const std::size_t dir, const Array< OneD, const NekDouble > &qfield, const Array< OneD, const NekDouble > &qFwd, const Array< OneD, const NekDouble > &qBwd, Array< OneD, NekDouble > &penaltyflux)
Definition: DiffusionLDG.cpp:415
Nektar::SolverUtils::DiffusionLDG::m_C11
NekDouble m_C11
Coefficient of penalty term.
Definition: DiffusionLDG.h:63
Vmath::Vadd
void Vadd(int n, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
Add vector z = x+y.
Definition: Vmath.cpp:359
Vmath::Smul
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*x.
Definition: Vmath.cpp:248
Vmath::Vsub
void Vsub(int n, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
Subtract vector z = x-y.
Definition: Vmath.cpp:419

References ApplyVectorBCs(), m_C11, m_traceNormals, Vmath::Smul(), Vmath::Vadd(), Vmath::Vcopy(), Vmath::Vmul(), and Vmath::Vsub().

Referenced by v_DiffuseTraceFlux().

◆ v_Diffuse()

void Nektar::SolverUtils::DiffusionLDG::v_Diffuse ( const std::size_t  nConvective,
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 
)
overrideprotectedvirtual

Reimplemented from Nektar::SolverUtils::Diffusion.

Definition at line 76 of file DiffusionLDG.cpp.

83 {
84  std::size_t nCoeffs = fields[0]->GetNcoeffs();
85 
86  Array<OneD, Array<OneD, NekDouble>> tmp{nConvectiveFields};
87  for (std::size_t i = 0; i < nConvectiveFields; ++i)
88  {
89  tmp[i] = Array<OneD, NekDouble>{nCoeffs, 0.0};
90  }
91 
92  DiffusionLDG::v_DiffuseCoeffs(nConvectiveFields, fields, inarray, tmp, pFwd,
93  pBwd);
94 
95  for (std::size_t i = 0; i < nConvectiveFields; ++i)
96  {
97  fields[i]->BwdTrans(tmp[i], outarray[i]);
98  }
99 }
Nektar::SolverUtils::DiffusionLDG::v_DiffuseCoeffs
virtual void v_DiffuseCoeffs(const std::size_t nConvective, 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) override
Definition: DiffusionLDG.cpp:101

References v_DiffuseCoeffs().

◆ v_DiffuseCalcDerivative()

void Nektar::SolverUtils::DiffusionLDG::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 
)
overrideprotectedvirtual

Diffusion Flux, calculate the physical derivatives.

Reimplemented from Nektar::SolverUtils::Diffusion.

Definition at line 165 of file DiffusionLDG.cpp.

171 {
172  std::size_t nConvectiveFields = fields.size();
173  std::size_t nDim = fields[0]->GetCoordim(0);
174  std::size_t nCoeffs = fields[0]->GetNcoeffs();
175  std::size_t nTracePts = fields[0]->GetTrace()->GetTotPoints();
176 
177  Array<OneD, NekDouble> tmp{nCoeffs};
178  TensorOfArray3D<NekDouble> flux{nDim};
179  for (std::size_t j = 0; j < nDim; ++j)
180  {
181  flux[j] = Array<OneD, Array<OneD, NekDouble>>{nConvectiveFields};
182  for (std::size_t i = 0; i < nConvectiveFields; ++i)
183  {
184  flux[j][i] = Array<OneD, NekDouble>{nTracePts, 0.0};
185  }
186  }
187 
188  NumFluxforScalar(fields, inarray, flux, pFwd, pBwd);
189 
190  for (std::size_t j = 0; j < nDim; ++j)
191  {
192  for (std::size_t i = 0; i < nConvectiveFields; ++i)
193  {
194  fields[i]->IProductWRTDerivBase(j, inarray[i], tmp);
195  Vmath::Neg(nCoeffs, tmp, 1);
196  fields[i]->AddTraceIntegral(flux[j][i], tmp);
197  fields[i]->SetPhysState(false);
198  fields[i]->MultiplyByElmtInvMass(tmp, tmp);
199  fields[i]->BwdTrans(tmp, qfield[j][i]);
200  }
201  }
202 }
Nektar::SolverUtils::DiffusionLDG::NumFluxforScalar
void NumFluxforScalar(const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &ufield, TensorOfArray3D< NekDouble > &uflux, const Array< OneD, Array< OneD, NekDouble >> &pFwd, const Array< OneD, Array< OneD, NekDouble >> &pBwd)
Definition: DiffusionLDG.cpp:226
Vmath::Neg
void Neg(int n, T *x, const int incx)
Negate x = -x.
Definition: Vmath.cpp:518

References Vmath::Neg(), and NumFluxforScalar().

◆ v_DiffuseCoeffs()

void Nektar::SolverUtils::DiffusionLDG::v_DiffuseCoeffs ( const std::size_t  nConvective,
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 
)
overrideprotectedvirtual

Reimplemented from Nektar::SolverUtils::Diffusion.

Definition at line 101 of file DiffusionLDG.cpp.

108 {
109  std::size_t nDim = fields[0]->GetCoordim(0);
110  std::size_t nPts = fields[0]->GetTotPoints();
111  std::size_t nCoeffs = fields[0]->GetNcoeffs();
112  std::size_t nTracePts = fields[0]->GetTrace()->GetTotPoints();
113 
114  Array<OneD, NekDouble> tmp{nCoeffs};
115 
116  TensorOfArray3D<NekDouble> qfield{nDim};
117  for (std::size_t j = 0; j < nDim; ++j)
118  {
119  qfield[j] = Array<OneD, Array<OneD, NekDouble>>{nConvectiveFields};
120  for (std::size_t i = 0; i < nConvectiveFields; ++i)
121  {
122  qfield[j][i] = Array<OneD, NekDouble>{nPts, 0.0};
123  }
124  }
125 
126  Array<OneD, Array<OneD, NekDouble>> traceflux{nConvectiveFields};
127  for (std::size_t i = 0; i < nConvectiveFields; ++i)
128  {
129  traceflux[i] = Array<OneD, NekDouble>{nTracePts, 0.0};
130  }
131 
132  DiffuseCalcDerivative(fields, inarray, qfield, pFwd, pBwd);
133 
134  // Initialize viscous tensor
135  Array<OneD, Array<OneD, Array<OneD, NekDouble>>> viscTensor{nDim};
136  for (std::size_t j = 0; j < nDim; ++j)
137  {
138  viscTensor[j] = Array<OneD, Array<OneD, NekDouble>>{nConvectiveFields};
139  for (std::size_t i = 0; i < nConvectiveFields; ++i)
140  {
141  viscTensor[j][i] = Array<OneD, NekDouble>{nPts, 0.0};
142  }
143  }
144  DiffuseVolumeFlux(fields, inarray, qfield, viscTensor);
145  DiffuseTraceFlux(fields, inarray, qfield, viscTensor, traceflux, pFwd,
146  pBwd);
147 
148  Array<OneD, Array<OneD, NekDouble>> qdbase{nDim};
149 
150  for (std::size_t i = 0; i < nConvectiveFields; ++i)
151  {
152  for (std::size_t j = 0; j < nDim; ++j)
153  {
154  qdbase[j] = viscTensor[j][i];
155  }
156  fields[i]->IProductWRTDerivBase(qdbase, tmp);
157 
158  Vmath::Neg(nCoeffs, tmp, 1);
159  fields[i]->AddTraceIntegral(traceflux[i], tmp);
160  fields[i]->SetPhysState(false);
161  fields[i]->MultiplyByElmtInvMass(tmp, outarray[i]);
162  }
163 }
Nektar::SolverUtils::Diffusion::DiffuseVolumeFlux
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.
Definition: Diffusion.h:206
Nektar::SolverUtils::Diffusion::DiffuseCalcDerivative
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)
Definition: Diffusion.cpp:194
Nektar::SolverUtils::Diffusion::DiffuseTraceFlux
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.
Definition: Diffusion.h:217

References Nektar::SolverUtils::Diffusion::DiffuseCalcDerivative(), Nektar::SolverUtils::Diffusion::DiffuseTraceFlux(), Nektar::SolverUtils::Diffusion::DiffuseVolumeFlux(), and Vmath::Neg().

Referenced by v_Diffuse().

◆ v_DiffuseTraceFlux()

void Nektar::SolverUtils::DiffusionLDG::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 
)
overrideprotectedvirtual

Diffusion term Trace Flux.

Reimplemented from Nektar::SolverUtils::Diffusion.

Definition at line 213 of file DiffusionLDG.cpp.

221 {
222  boost::ignore_unused(qfield, pFwd, pBwd, nonZeroIndex);
223  NumFluxforVector(fields, inarray, viscTensor, TraceFlux);
224 }
Nektar::SolverUtils::DiffusionLDG::NumFluxforVector
void NumFluxforVector(const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &ufield, TensorOfArray3D< NekDouble > &qfield, Array< OneD, Array< OneD, NekDouble >> &qflux)
Build the numerical flux for the 2nd order derivatives todo: add variable coeff and h dependence to p...
Definition: DiffusionLDG.cpp:354

References NumFluxforVector().

◆ v_DiffuseVolumeFlux()

void Nektar::SolverUtils::DiffusionLDG::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 
)
overrideprotectedvirtual

Diffusion Volume Flux.

Reimplemented from Nektar::SolverUtils::Diffusion.

Definition at line 204 of file DiffusionLDG.cpp.

209 {
210  boost::ignore_unused(fields, nonZeroIndex);
211  m_fluxVector(inarray, qfield, viscTensor);
212 }
Nektar::SolverUtils::Diffusion::m_fluxVector
DiffusionFluxVecCB m_fluxVector
Definition: Diffusion.h:402

References Nektar::SolverUtils::Diffusion::m_fluxVector.

◆ v_InitObject()

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

Reimplemented from Nektar::SolverUtils::Diffusion.

Definition at line 53 of file DiffusionLDG.cpp.

56 {
57  m_session = pSession;
58 
59  m_session->LoadSolverInfo("ShockCaptureType", m_shockCaptureType, "Off");
60 
61  // Set up penalty term for LDG
62  m_session->LoadParameter("LDGc11", m_C11, 1.0);
63 
64  // Setting up the normals
65  std::size_t nDim = pFields[0]->GetCoordim(0);
66  std::size_t nTracePts = pFields[0]->GetTrace()->GetTotPoints();
67 
68  m_traceNormals = Array<OneD, Array<OneD, NekDouble>>{nDim};
69  for (std::size_t i = 0; i < nDim; ++i)
70  {
71  m_traceNormals[i] = Array<OneD, NekDouble>{nTracePts};
72  }
73  pFields[0]->GetTrace()->GetNormals(m_traceNormals);
74 }
Nektar::SolverUtils::DiffusionLDG::m_session
LibUtilities::SessionReaderSharedPtr m_session
Definition: DiffusionLDG.h:66

References m_C11, m_session, m_shockCaptureType, and m_traceNormals.

Member Data Documentation

◆ m_C11

NekDouble Nektar::SolverUtils::DiffusionLDG::m_C11
protected

Coefficient of penalty term.

Definition at line 63 of file DiffusionLDG.h.

Referenced by NumFluxforVector(), and v_InitObject().

◆ m_session

LibUtilities::SessionReaderSharedPtr Nektar::SolverUtils::DiffusionLDG::m_session
protected

Definition at line 66 of file DiffusionLDG.h.

Referenced by v_InitObject().

◆ m_shockCaptureType

std::string Nektar::SolverUtils::DiffusionLDG::m_shockCaptureType
protected

Definition at line 60 of file DiffusionLDG.h.

Referenced by v_InitObject().

◆ m_traceNormals

Array<OneD, Array<OneD, NekDouble> > Nektar::SolverUtils::DiffusionLDG::m_traceNormals
protected

Definition at line 65 of file DiffusionLDG.h.

Referenced by ApplyVectorBCs(), NumFluxforScalar(), NumFluxforVector(), and v_InitObject().

◆ type

std::string Nektar::SolverUtils::DiffusionLDG::type
static
Initial value:
= GetDiffusionFactory().RegisterCreatorFunction(
"LDG", DiffusionLDG::create, "Local Discontinuous Galkerin")
Nektar::LibUtilities::NekFactory::RegisterCreatorFunction
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
Definition: NekFactory.hpp:198
Nektar::SolverUtils::DiffusionLDG::create
static DiffusionSharedPtr create(std::string diffType)
Definition: DiffusionLDG.h:49
Nektar::SolverUtils::GetDiffusionFactory
DiffusionFactory & GetDiffusionFactory()
Definition: Diffusion.cpp:41

Definition at line 55 of file DiffusionLDG.h.

Generated on Sat Mar 25 2023 12:36:00 for Nektar++ by doxygen 1.9.1