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)

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)

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	GetArtificialViscosity (const Array< OneD, Array< OneD, NekDouble > > &physfield, Array< OneD, NekDouble > &mu)
	Calculate the artificial viscosity. More...

void	SetElmtHP (const Array< OneD, NekDouble > &hOverP)
	Set h/p scaling. 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...

Array< OneD, NekDouble >	m_hOverP
	h/p scaling More...

Detailed Description

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

Definition at line 48 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.

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

     : ArtificialDiffusion(pSession, pFields, spacedim)
 {
     m_session->LoadParameter ("SensorOffset",  m_offset, 1);
 }

◆ ~NonSmoothShockCapture()

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

inlineprivatevirtual

Definition at line 81 of file NonSmoothShockCapture.h.

References m_offset.

81 {};

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 55 of file NonSmoothShockCapture.h.

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

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

◆ v_GetArtificialViscosity()

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

protectedvirtual

Implements Nektar::ArtificialDiffusion.

Definition at line 56 of file NonSmoothShockCapture.cpp.

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

 {
     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_hOverP[e];
         Vmath::Smul(nElmtPoints, LambdaElmt, tmp = mu + physOffset, 1,
             tmp = mu + physOffset, 1);
     }
 }

Friends And Related Function Documentation

◆ MemoryManager< NonSmoothShockCapture >

friend class MemoryManager< NonSmoothShockCapture >

friend

Definition at line 52 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 67 of file NonSmoothShockCapture.h.

◆ m_offset

int Nektar::NonSmoothShockCapture::m_offset

private

Parameters.

Definition at line 81 of file NonSmoothShockCapture.h.

Referenced by NonSmoothShockCapture(), v_GetArtificialViscosity(), and ~NonSmoothShockCapture().

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Private Member Functions

Private Attributes

Friends

Additional Inherited Members