Encapsulates the artificial diffusion used in shock capture. More...

#include <ArtificialDiffusion.h>

Inheritance diagram for Nektar::ArtificialDiffusion:

Collaboration diagram for Nektar::ArtificialDiffusion:

Public Member Functions
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...

Protected Member Functions
	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_GetArtificialViscosity (const Array< OneD, Array< OneD, NekDouble > > &physfield, Array< OneD, NekDouble > &mu)=0

Protected Attributes
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_FacL
	Parameters. More...

NekDouble	m_FacH

NekDouble	m_hFactor

NekDouble	m_C1

NekDouble	m_C2

NekDouble	m_mu0

NekDouble	m_Skappa

NekDouble	m_Kappa

Detailed Description

Encapsulates the artificial diffusion used in shock capture.

Definition at line 69 of file ArtificialDiffusion.h.

Constructor & Destructor Documentation

virtual Nektar::ArtificialDiffusion::~ArtificialDiffusion ( )

inlinevirtual

Definition at line 72 of file ArtificialDiffusion.h.

72 {}

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

protected

Constructor.

Definition at line 52 of file ArtificialDiffusion.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::LibUtilities::NekFactory< tKey, tBase, >::CreateInstance(), GetArtificialViscosity(), Nektar::SolverUtils::GetDiffusionFactory(), m_C1, m_C2, m_diffusion, m_FacH, m_FacL, m_fields, m_hFactor, m_Kappa, m_mu0, m_session, m_Skappa, and m_varConv.

         : m_session(pSession),
         m_fields(pFields)
 {
     m_session->LoadParameter ("FL",            m_FacL,          0.0);
     m_session->LoadParameter ("FH",            m_FacH,          0.0);
     m_session->LoadParameter ("hFactor",       m_hFactor,       1.0);
     m_session->LoadParameter ("C1",            m_C1,            3.0);
     m_session->LoadParameter ("C2",            m_C2,            5.0);
     m_session->LoadParameter ("mu0",           m_mu0,           1.0);
     m_session->LoadParameter ("Skappa",        m_Skappa,        -2.048);
     m_session->LoadParameter ("Kappa",         m_Kappa,         0.0);
 
     // Create auxiliary object to convert variables
     m_varConv = MemoryManager<VariableConverter>::AllocateSharedPtr(
                 m_session, spacedim);
 
     m_diffusion = SolverUtils::GetDiffusionFactory()
                                     .CreateInstance("LDG", "LDG");
     m_diffusion->SetArtificialDiffusionVector(
                         &ArtificialDiffusion::GetArtificialViscosity, this);
     m_diffusion->InitObject (m_session, m_fields);
 }

Member Function Documentation

void Nektar::ArtificialDiffusion::DoArtificialDiffusion	(	const Array< OneD, const Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	outarray
	)

Apply the artificial diffusion.

Definition at line 82 of file ArtificialDiffusion.cpp.

References v_DoArtificialDiffusion().

 {
     v_DoArtificialDiffusion(inarray, outarray);
 }

void Nektar::ArtificialDiffusion::GetArtificialViscosity	(	const Array< OneD, Array< OneD, NekDouble > > &	physfield,
		Array< OneD, NekDouble > &	mu
	)

Calculate the artificial viscosity.

Definition at line 121 of file ArtificialDiffusion.cpp.

References v_GetArtificialViscosity().

Referenced by ArtificialDiffusion().

 {
     v_GetArtificialViscosity(physfield, mu);
 }

void Nektar::ArtificialDiffusion::v_DoArtificialDiffusion	(	const Array< OneD, const Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	outarray
	)

protectedvirtual

Reimplemented in Nektar::SmoothShockCapture.

Definition at line 92 of file ArtificialDiffusion.cpp.

References m_diffusion, m_fields, and Vmath::Vadd().

Referenced by DoArtificialDiffusion().

 {
     int i;
     int nvariables = inarray.num_elements();
     int npoints    = m_fields[0]->GetNpoints();
 
     Array<OneD, Array<OneD, NekDouble> > outarrayDiff(nvariables);
 
     for (i = 0; i < nvariables; ++i)
     {
         outarrayDiff[i] = Array<OneD, NekDouble>(npoints, 0.0);
     }
 
     m_diffusion->Diffuse(nvariables, m_fields, inarray, outarrayDiff);
 
     for (i = 0; i < nvariables; ++i)
     {
         Vmath::Vadd(npoints,
                     outarray[i], 1,
                     outarrayDiff[i], 1,
                     outarray[i], 1);
     }
 }

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