Non Smooth artificial diffusion for shock capture for compressible flow problems. More...

#include <NonSmoothShockCapture.h>

Inheritance diagram for Nektar::NonSmoothShockCapture:

Static Public Member Functions
static ArtificialDiffusionSharedPtr	create (const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const int spacedim)
	Creates an instance of this class. More...

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

Protected Member Functions
virtual void	v_GetArtificialViscosity (const Array< OneD, Array< OneD, NekDouble >> &physfield, Array< OneD, NekDouble > &mu) override

Protected Member Functions inherited from Nektar::ArtificialDiffusion
	ArtificialDiffusion (const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const int spacedim)
	Constructor. More...

virtual void	v_DoArtificialDiffusion (const Array< OneD, const Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)

virtual void	v_DoArtificialDiffusionCoeff (const Array< OneD, const Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)

void	GetFluxVector (const Array< OneD, Array< OneD, NekDouble >> &inarray, const Array< OneD, Array< OneD, Array< OneD, NekDouble >>> &qfield, Array< OneD, Array< OneD, Array< OneD, NekDouble >>> &viscousTensor)
	Return the flux vector for the artificial viscosity operator. More...

Private Member Functions
	NonSmoothShockCapture (const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const int spacedim)

virtual	~NonSmoothShockCapture (void)

Private Attributes
int	m_offset
	Parameters. More...

Friends
class	MemoryManager< NonSmoothShockCapture >

Additional Inherited Members
Public Member Functions inherited from Nektar::ArtificialDiffusion
virtual	~ArtificialDiffusion ()

void	DoArtificialDiffusion (const Array< OneD, const Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
	Apply the artificial diffusion. More...

void	DoArtificialDiffusionCoeff (const Array< OneD, const Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
	Apply the artificial diffusion the outarray is in coeff space. More...

void	GetArtificialViscosity (const Array< OneD, Array< OneD, NekDouble >> &physfield, Array< OneD, NekDouble > &mu)
	Calculate the artificial viscosity. More...

Protected Attributes inherited from Nektar::ArtificialDiffusion
LibUtilities::SessionReaderSharedPtr	m_session
	Session reader. More...

Array< OneD, MultiRegions::ExpListSharedPtr >	m_fields
	Array of fields. More...

VariableConverterSharedPtr	m_varConv
	Auxiliary object to convert variables. More...

SolverUtils::DiffusionSharedPtr	m_diffusion
	LDG Diffusion operator. More...

NekDouble	m_mu0
	Constant scaling. More...

Detailed Description

Non Smooth artificial diffusion for shock capture for compressible flow problems.

Definition at line 47 of file NonSmoothShockCapture.h.

Constructor & Destructor Documentation

◆ NonSmoothShockCapture()

Nektar::NonSmoothShockCapture::NonSmoothShockCapture	(	const LibUtilities::SessionReaderSharedPtr &	pSession,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		const int	spacedim
	)

private

Definition at line 47 of file NonSmoothShockCapture.cpp.

     : ArtificialDiffusion(pSession, pFields, spacedim)
 {
     m_session->LoadParameter("SensorOffset", m_offset, 1);
     // init h/p
     m_varConv->SetElmtMinHP(m_fields);
 }

References Nektar::ArtificialDiffusion::m_fields, m_offset, Nektar::ArtificialDiffusion::m_session, and Nektar::ArtificialDiffusion::m_varConv.

◆ ~NonSmoothShockCapture()

virtual Nektar::NonSmoothShockCapture::~NonSmoothShockCapture ( void )

inlineprivatevirtual

Definition at line 78 of file NonSmoothShockCapture.h.

78 {};

Member Function Documentation

◆ create()

static ArtificialDiffusionSharedPtr Nektar::NonSmoothShockCapture::create	(	const LibUtilities::SessionReaderSharedPtr &	pSession,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		const int	spacedim
	)

inlinestatic

Creates an instance of this class.

Definition at line 53 of file NonSmoothShockCapture.h.

     {
         ArtificialDiffusionSharedPtr p =
             MemoryManager<NonSmoothShockCapture>::AllocateSharedPtr(
                 pSession, pFields, spacedim);
         return p;
     }

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), and CellMLToNektar.cellml_metadata::p.

◆ v_GetArtificialViscosity()

void Nektar::NonSmoothShockCapture::v_GetArtificialViscosity	(	const Array< OneD, Array< OneD, NekDouble >> &	physfield,
		Array< OneD, NekDouble > &	mu
	)

overrideprotectedvirtual

Implements Nektar::ArtificialDiffusion.

Definition at line 58 of file NonSmoothShockCapture.cpp.

 {
     int nPts      = m_fields[0]->GetTotPoints();
     int nElements = m_fields[0]->GetExpSize();
  
     // Determine the maximum wavespeed
     Array<OneD, NekDouble> Lambda(nPts, 0.0);
     Array<OneD, NekDouble> soundspeed(nPts, 0.0);
     Array<OneD, NekDouble> absVelocity(nPts, 0.0);
     m_varConv->GetSoundSpeed(physfield, soundspeed);
     m_varConv->GetAbsoluteVelocity(physfield, absVelocity);
     Vmath::Vadd(nPts, absVelocity, 1, soundspeed, 1, Lambda, 1);
  
     // Compute sensor based on rho
     Array<OneD, NekDouble> Sensor(nPts, 0.0);
     m_varConv->GetSensor(m_fields[0], physfield, Sensor, mu, m_offset);
  
     // Scale AV
     Array<OneD, NekDouble> tmp;
     for (int e = 0; e < nElements; e++)
     {
         int physOffset  = m_fields[0]->GetPhys_Offset(e);
         int nElmtPoints = m_fields[0]->GetExp(e)->GetTotPoints();
  
         // Get max wavespeed per element
         NekDouble LambdaElmt = 0.0;
         LambdaElmt = Vmath::Vmax(nElmtPoints, tmp = Lambda + physOffset, 1);
  
         // Scale viscosity by the maximum wave speed and h/p
         LambdaElmt *= m_mu0 * m_varConv->GetElmtMinHP()[e];
         Vmath::Smul(nElmtPoints, LambdaElmt, tmp = mu + physOffset, 1,
                     tmp = mu + physOffset, 1);
     }
 }

References Nektar::ArtificialDiffusion::m_fields, Nektar::ArtificialDiffusion::m_mu0, m_offset, Nektar::ArtificialDiffusion::m_varConv, Vmath::Smul(), Vmath::Vadd(), and Vmath::Vmax().

Friends And Related Function Documentation

◆ MemoryManager< NonSmoothShockCapture >

friend class MemoryManager< NonSmoothShockCapture >

friend

Definition at line 1 of file NonSmoothShockCapture.h.

Member Data Documentation

◆ className

std::string Nektar::NonSmoothShockCapture::className

static

Initial value:

=
    GetArtificialDiffusionFactory().RegisterCreatorFunction(
        "NonSmooth", NonSmoothShockCapture::create,
        "NonSmooth artificial diffusion for shock capture.")

Name of the class.

Definition at line 65 of file NonSmoothShockCapture.h.

◆ m_offset

int Nektar::NonSmoothShockCapture::m_offset

private

Parameters.

Definition at line 81 of file NonSmoothShockCapture.h.

Referenced by NonSmoothShockCapture(), and v_GetArtificialViscosity().

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Private Member Functions

Private Attributes

Friends

Additional Inherited Members