Nektar::LEEUpwindSolver Class Reference

#include <LEEUpwindSolver.h>

Inheritance diagram for Nektar::LEEUpwindSolver:

Static Public Member Functions
static RiemannSolverSharedPtr	create (const LibUtilities::SessionReaderSharedPtr &pSession)

Static Public Attributes
static std::string	solverName

Protected Member Functions
	LEEUpwindSolver (const LibUtilities::SessionReaderSharedPtr &pSession)
virtual void	v_PointSolve (NekDouble pL, NekDouble rhoL, NekDouble rhouL, NekDouble rhovL, NekDouble rhowL, NekDouble pR, NekDouble rhoR, NekDouble rhouR, NekDouble rhovR, NekDouble rhowR, 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 &rhouF, NekDouble &rhovF, NekDouble &rhowF)
	Lax-Friedrichs Riemann solver. More...
Protected Member Functions inherited from Nektar::LEESolver
	LEESolver (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)
Protected Member Functions inherited from Nektar::AcousticSolver
	AcousticSolver (const LibUtilities::SessionReaderSharedPtr &pSession)
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 LEEUpwindSolver.h.

Constructor & Destructor Documentation

LEEUpwindSolver()

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

protected

Definition at line 52 of file LEEUpwindSolver.cpp.

    : LEESolver(pSession)
{
}

Member Function Documentation

create()

static RiemannSolverSharedPtr Nektar::LEEUpwindSolver::create ( const LibUtilities::SessionReaderSharedPtr &  pSession )

inlinestatic

Definition at line 50 of file LEEUpwindSolver.h.

    {
        return RiemannSolverSharedPtr(new LEEUpwindSolver(pSession));
    }

v_PointSolve()

void Nektar::LEEUpwindSolver::v_PointSolve ( NekDouble  pL, NekDouble  rhoL, NekDouble  rhouL, NekDouble  rhovL, NekDouble  rhowL, NekDouble  pR, NekDouble  rhoR, NekDouble  rhouR, NekDouble  rhovR, NekDouble  rhowR, 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 &  rhouF, NekDouble &  rhovF, NekDouble &  rhowF  )

protectedvirtual

Lax-Friedrichs Riemann solver.

Parameters

pL	Perturbation pressure left state
rhoL	Perturbation density left state
pR	Perturbation pressure right state
rhoR	Perturbation density right state
rhouL	x perturbation velocity component left state
rhouR	x perturbation velocity component right state
rhovL	y perturbation velocity component left state
rhovR	y perturbation velocity component right state
rhowL	z perturbation velocity component left state
rhowR	z perturbation velocity component right state
c0sqL	Base pressure left state
c0sqR	Base pressure right state
rho0L	Base density left state
rho0R	Base density right state
u0L	Base x velocity component left state
u0R	Base x velocity component right state
v0L	Base y velocity component left state
v0R	Base y velocity component right state
w0L	Base z velocity component left state
w0R	Base z velocity component right state
pF	Computed Riemann flux for perturbation pressure
rhoF	Computed Riemann flux for perturbation density
rhouF	Computed Riemann flux for x perturbation velocity component
rhovF	Computed Riemann flux for y perturbation velocity component
rhowF	Computed Riemann flux for z perturbation velocity component

Implements Nektar::AcousticSolver.

Definition at line 87 of file LEEUpwindSolver.cpp.

{
    boost::ignore_unused(rho0L, v0L, w0L, rho0R, v0R, w0R);

    // Speed of sound
    NekDouble c0L = sqrt(c0sqL);
    NekDouble c0R = sqrt(c0sqR);
    NekDouble c0M = (c0L + c0R) / 2.0;

    NekDouble u0M = (u0L + u0R) / 2.0;

    pF    = 0.0;
    rhoF  = 0.0;
    rhouF = 0.0;
    rhovF = 0.0;
    rhowF = 0.0;

    // lambda_1,2,3
    if (u0M > 0)
    {
        rhoF  = rhoF + u0L * (c0sqL * rhoL - pL) / c0sqL;
        rhovF = rhovF + rhovL * u0L;
        rhowF = rhowF + rhowL * u0L;
    }
    else
    {
        rhoF  = rhoF + u0R * (c0sqR * rhoR - pR) / c0sqR;
        rhovF = rhovF + rhovR * u0R;
        rhowF = rhowF + rhowR * u0R;
    }

    // lambda_4
    if (u0M - c0M > 0)
    {
        pF   = pF + 0.5 * (c0L - u0L) * (rhouL * c0L - pL);
        rhoF = rhoF + 0.5 * (c0L - u0L) * (rhouL * c0sqL - c0L * pL) /
                          pow(c0sqL, 3.0 / 2.0);
        rhouF = rhouF + 0.5 * (c0L - u0L) * (-rhouL * c0L + pL) / c0L;
    }
    else
    {
        pF   = pF + 0.5 * (c0R - u0R) * (rhouR * c0R - pR);
        rhoF = rhoF + 0.5 * (c0R - u0R) * (rhouR * c0sqR - c0R * pR) /
                          pow(c0sqR, 3.0 / 2.0);
        rhouF = rhouF + 0.5 * (c0R - u0R) * (-rhouR * c0R + pR) / c0R;
    }

    // lambda_5
    if (u0M + c0M > 0)
    {
        pF   = pF + 0.5 * (c0L + u0L) * (rhouL * c0L + pL);
        rhoF = rhoF + 0.5 * (c0L + u0L) * (rhouL * c0sqL + c0L * pL) /
                          pow(c0sqL, 3.0 / 2.0);
        rhouF = rhouF + 0.5 * (c0L + u0L) * (rhouL * c0L + pL) / c0L;
    }
    else
    {
        pF   = pF + 0.5 * (c0R + u0R) * (rhouR * c0R + pR);
        rhoF = rhoF + 0.5 * (c0R + u0R) * (rhouR * c0sqR + c0R * pR) /
                          pow(c0sqR, 3.0 / 2.0);
        rhouF = rhouF + 0.5 * (c0R + u0R) * (rhouR * c0R + pR) / c0R;
    }
}

Member Data Documentation

solverName

std::string Nektar::LEEUpwindSolver::solverName

static

Initial value:

=
    SolverUtils::GetRiemannSolverFactory().RegisterCreatorFunction(
        "LEEUpwind", LEEUpwindSolver::create, "Upwind Solver for LEE")

Definition at line 56 of file LEEUpwindSolver.h.