#include <AcousticSolver.h>

Inheritance diagram for Nektar::AcousticSolver:

Protected Member Functions
	AcousticSolver (const LibUtilities::SessionReaderSharedPtr &pSession)

virtual void	v_Solve (const int nDim, const Array< OneD, const Array< OneD, NekDouble >> &Fwd, const Array< OneD, const Array< OneD, NekDouble >> &Bwd, Array< OneD, Array< OneD, NekDouble >> &flux)

virtual void	v_PointSolve (NekDouble pL, NekDouble rhoL, NekDouble uL, NekDouble vL, NekDouble wL, NekDouble pR, NekDouble rhoR, NekDouble uR, NekDouble vR, NekDouble wR, NekDouble c0sqL, NekDouble rho0L, NekDouble u0L, NekDouble v0L, NekDouble w0L, NekDouble c0sqR, NekDouble rho0R, NekDouble u0R, NekDouble v0R, NekDouble w0R, NekDouble &pF, NekDouble &rhoF, NekDouble &uF, NekDouble &vF, NekDouble &wF)=0

void	GetRotBasefield (Array< OneD, Array< OneD, NekDouble >> &bfFwd, Array< OneD, Array< OneD, NekDouble >> &bfBwd)

Protected Member Functions inherited from Nektar::SolverUtils::RiemannSolver
SOLVER_UTILS_EXPORT	RiemannSolver (const LibUtilities::SessionReaderSharedPtr &pSession)

virtual SOLVER_UTILS_EXPORT	~RiemannSolver ()

virtual void	v_Solve (const int nDim, const Array< OneD, const Array< OneD, NekDouble > > &Fwd, const Array< OneD, const Array< OneD, NekDouble > > &Bwd, Array< OneD, Array< OneD, NekDouble > > &flux)=0

void	GenerateRotationMatrices (const Array< OneD, const Array< OneD, NekDouble > > &normals)
	Generate rotation matrices for 3D expansions. More...

void	FromToRotation (Array< OneD, const NekDouble > &from, Array< OneD, const NekDouble > &to, NekDouble *mat)
	A function for creating a rotation matrix that rotates a vector from into another vector to. More...

SOLVER_UTILS_EXPORT void	rotateToNormal (const Array< OneD, const Array< OneD, NekDouble > > &inarray, const Array< OneD, const Array< OneD, NekDouble > > &normals, const Array< OneD, const Array< OneD, NekDouble > > &vecLocs, Array< OneD, Array< OneD, NekDouble > > &outarray)
	Rotate a vector field to trace normal. More...

SOLVER_UTILS_EXPORT void	rotateFromNormal (const Array< OneD, const Array< OneD, NekDouble > > &inarray, const Array< OneD, const Array< OneD, NekDouble > > &normals, const Array< OneD, const Array< OneD, NekDouble > > &vecLocs, Array< OneD, Array< OneD, NekDouble > > &outarray)
	Rotate a vector field from trace normal. More...

SOLVER_UTILS_EXPORT bool	CheckScalars (std::string name)
	Determine whether a scalar has been defined in m_scalars. More...

SOLVER_UTILS_EXPORT bool	CheckVectors (std::string name)
	Determine whether a vector has been defined in m_vectors. More...

SOLVER_UTILS_EXPORT bool	CheckParams (std::string name)
	Determine whether a parameter has been defined in m_params. More...

SOLVER_UTILS_EXPORT bool	CheckAuxScal (std::string name)
	Determine whether a scalar has been defined in m_auxScal. More...

SOLVER_UTILS_EXPORT bool	CheckAuxVec (std::string name)
	Determine whether a vector has been defined in m_auxVec. More...

Additional Inherited Members
Public Member Functions inherited from Nektar::SolverUtils::RiemannSolver
SOLVER_UTILS_EXPORT void	Solve (const int nDim, const Array< OneD, const Array< OneD, NekDouble > > &Fwd, const Array< OneD, const Array< OneD, NekDouble > > &Bwd, Array< OneD, Array< OneD, NekDouble > > &flux)
	Perform the Riemann solve given the forwards and backwards spaces. More...

template<typename FuncPointerT , typename ObjectPointerT >
void	SetScalar (std::string name, FuncPointerT func, ObjectPointerT obj)

void	SetScalar (std::string name, RSScalarFuncType fp)

template<typename FuncPointerT , typename ObjectPointerT >
void	SetVector (std::string name, FuncPointerT func, ObjectPointerT obj)

void	SetVector (std::string name, RSVecFuncType fp)

template<typename FuncPointerT , typename ObjectPointerT >
void	SetParam (std::string name, FuncPointerT func, ObjectPointerT obj)

void	SetParam (std::string name, RSParamFuncType fp)

template<typename FuncPointerT , typename ObjectPointerT >
void	SetAuxScal (std::string name, FuncPointerT func, ObjectPointerT obj)

template<typename FuncPointerT , typename ObjectPointerT >
void	SetAuxVec (std::string name, FuncPointerT func, ObjectPointerT obj)

std::map< std::string, RSScalarFuncType > &	GetScalars ()

std::map< std::string, RSVecFuncType > &	GetVectors ()

std::map< std::string, RSParamFuncType > &	GetParams ()

Public Attributes inherited from Nektar::SolverUtils::RiemannSolver
int	m_spacedim

Protected Attributes inherited from Nektar::SolverUtils::RiemannSolver
bool	m_requiresRotation
	Indicates whether the Riemann solver requires a rotation to be applied to the velocity fields. More...

std::map< std::string, RSScalarFuncType >	m_scalars
	Map of scalar function types. More...

std::map< std::string, RSVecFuncType >	m_vectors
	Map of vector function types. More...

std::map< std::string, RSParamFuncType >	m_params
	Map of parameter function types. More...

std::map< std::string, RSScalarFuncType >	m_auxScal
	Map of auxiliary scalar function types. More...

std::map< std::string, RSVecFuncType >	m_auxVec
	Map of auxiliary vector function types. More...

Array< OneD, Array< OneD, NekDouble > >	m_rotMat
	Rotation matrices for each trace quadrature point. More...

Array< OneD, Array< OneD, Array< OneD, NekDouble > > >	m_rotStorage
	Rotation storage. More...

Detailed Description

Definition at line 47 of file AcousticSolver.h.

Constructor & Destructor Documentation

◆ AcousticSolver()

Nektar::AcousticSolver::AcousticSolver ( const LibUtilities::SessionReaderSharedPtr & pSession )

protected

Definition at line 46 of file RiemannSolvers/AcousticSolver.cpp.

References Nektar::SolverUtils::RiemannSolver::m_requiresRotation.

     : RiemannSolver(pSession)
 {
     m_requiresRotation = true;
 }

Member Function Documentation

◆ GetRotBasefield()

void Nektar::AcousticSolver::GetRotBasefield	(	Array< OneD, Array< OneD, NekDouble >> &	bfFwd,
		Array< OneD, Array< OneD, NekDouble >> &	bfBwd
	)

protected

Definition at line 117 of file RiemannSolvers/AcousticSolver.cpp.

References ASSERTL1, Nektar::SolverUtils::RiemannSolver::CheckVectors(), Nektar::SolverUtils::RiemannSolver::m_vectors, and Nektar::SolverUtils::RiemannSolver::rotateToNormal().

Referenced by v_Solve(), and Nektar::LEESolver::v_Solve().

 {
     ASSERTL1(CheckVectors("N"), "N not defined.");
     ASSERTL1(CheckVectors("basefieldFwdBwd"), "basefieldFwdBwd not defined.");
     const Array<OneD, const Array<OneD, NekDouble>> normals = m_vectors["N"]();
     const Array<OneD, const Array<OneD, NekDouble>> basefieldFwdBwd =
         m_vectors["basefieldFwdBwd"]();
 
     int nBF  = basefieldFwdBwd.num_elements() / 2;
     int nDim = normals.num_elements();
 
     Array<OneD, Array<OneD, NekDouble>> basefieldFwd(nBF);
     Array<OneD, Array<OneD, NekDouble>> basefieldBwd(nBF);
 
     for (int i = 0; i < nBF; i++)
     {
         int j           = nBF + i;
         basefieldFwd[i] = basefieldFwdBwd[i];
         basefieldBwd[i] = basefieldFwdBwd[j];
     }
 
     Array<OneD, Array<OneD, NekDouble>> baseVecLocs(1);
     baseVecLocs[0] = Array<OneD, NekDouble>(nDim);
     for (int i = 0; i < nDim; ++i)
     {
         baseVecLocs[0][i] = i + 2;
     }
     rotateToNormal(basefieldFwd, normals, baseVecLocs, bfFwd);
     rotateToNormal(basefieldBwd, normals, baseVecLocs, bfBwd);
 }

◆ v_PointSolve()

virtual void Nektar::AcousticSolver::v_PointSolve	(	NekDouble	pL,
		NekDouble	rhoL,
		NekDouble	uL,
		NekDouble	vL,
		NekDouble	wL,
		NekDouble	pR,
		NekDouble	rhoR,
		NekDouble	uR,
		NekDouble	vR,
		NekDouble	wR,
		NekDouble	c0sqL,
		NekDouble	rho0L,
		NekDouble	u0L,
		NekDouble	v0L,
		NekDouble	w0L,
		NekDouble	c0sqR,
		NekDouble	rho0R,
		NekDouble	u0R,
		NekDouble	v0R,
		NekDouble	w0R,
		NekDouble &	pF,
		NekDouble &	rhoF,
		NekDouble &	uF,
		NekDouble &	vF,
		NekDouble &	wF
	)

protectedpure virtual

Implemented in Nektar::APELaxFriedrichsSolver, Nektar::LEELaxFriedrichsSolver, Nektar::APEUpwindSolver, and Nektar::LEEUpwindSolver.

Referenced by v_Solve(), and Nektar::LEESolver::v_Solve().

◆ v_Solve()

void Nektar::AcousticSolver::v_Solve	(	const int	nDim,
		const Array< OneD, const Array< OneD, NekDouble >> &	Fwd,
		const Array< OneD, const Array< OneD, NekDouble >> &	Bwd,
		Array< OneD, Array< OneD, NekDouble >> &	flux
	)

protectedvirtual

Reimplemented in Nektar::LEESolver.

Definition at line 56 of file RiemannSolvers/AcousticSolver.cpp.

References GetRotBasefield(), and v_PointSolve().

 {
     int nTracePts = Fwd[0].num_elements();
 
     Array<OneD, Array<OneD, NekDouble>> bfFwd(nDim + 2);
     Array<OneD, Array<OneD, NekDouble>> bfBwd(nDim + 2);
     for (int i = 0; i < nDim + 2; i++)
     {
         bfFwd[i] = Array<OneD, NekDouble>(nTracePts);
         bfBwd[i] = Array<OneD, NekDouble>(nTracePts);
     }
 
     GetRotBasefield(bfFwd, bfBwd);
 
     int expDim = nDim;
     NekDouble vF, wF, rhoF;
 
     if (expDim == 1)
     {
         for (int i = 0; i < nTracePts; ++i)
         {
             v_PointSolve(
                   Fwd[0][i],         0.0,   Fwd[1][i], 0.0,  0.0,
                   Bwd[0][i],         0.0,   Bwd[1][i], 0.0,  0.0,
                 bfFwd[0][i], bfFwd[1][i], bfFwd[2][i], 0.0,  0.0,
                 bfBwd[0][i], bfBwd[1][i], bfBwd[2][i], 0.0,  0.0,
                  flux[0][i],        rhoF,  flux[1][i],  vF,   wF);
         }
     }
     else if (expDim == 2)
     {
         for (int i = 0; i < nTracePts; ++i)
         {
             v_PointSolve(
                   Fwd[0][i],         0.0,   Fwd[1][i],   Fwd[2][i],  0.0,
                   Bwd[0][i],         0.0,   Bwd[1][i],   Bwd[2][i],  0.0,
                 bfFwd[0][i], bfFwd[1][i], bfFwd[2][i], bfFwd[3][i],  0.0,
                 bfBwd[0][i], bfBwd[1][i], bfBwd[2][i], bfBwd[3][i],  0.0,
                  flux[0][i],        rhoF,  flux[1][i],  flux[2][i],   wF);
         }
     }
     else if (expDim == 3)
     {
         for (int i = 0; i < nTracePts; ++i)
         {
             v_PointSolve(
                   Fwd[0][i],         0.0,   Fwd[1][i],   Fwd[2][i],   Fwd[3][i],
                   Bwd[0][i],         0.0,   Bwd[1][i],   Bwd[2][i],   Bwd[3][i],
                 bfFwd[0][i], bfFwd[1][i], bfFwd[2][i], bfFwd[3][i], bfFwd[4][i],
                 bfBwd[0][i], bfBwd[1][i], bfBwd[2][i], bfBwd[3][i], bfBwd[4][i],
                  flux[0][i],        rhoF,  flux[1][i],  flux[2][i],  flux[3][i]);
         }
     }
 }

Protected Member Functions