Nektar++
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Static Public Member Functions | Static Public Attributes | Protected Member Functions | List of all members
Nektar::LaxFriedrichsSolver Class Reference

#include <LaxFriedrichsSolver.h>

Inheritance diagram for Nektar::LaxFriedrichsSolver:
Inheritance graph
[legend]
Collaboration diagram for Nektar::LaxFriedrichsSolver:
Collaboration graph
[legend]

Static Public Member Functions

static RiemannSolverSharedPtr create ()
 
static RiemannSolverSharedPtr create ()
 
static RiemannSolverSharedPtr create ()
 

Static Public Attributes

static std::string solverName
 

Protected Member Functions

 LaxFriedrichsSolver ()
 
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 p0, NekDouble rho0, NekDouble u0, NekDouble v0, NekDouble w0, NekDouble &pF, NekDouble &rhoF, NekDouble &uF, NekDouble &vF, NekDouble &wF)
 Lax-Friedrichs Riemann solver. More...
 
 LaxFriedrichsSolver ()
 
virtual void v_PointSolve (double rhoL, double rhouL, double rhovL, double rhowL, double EL, double rhoR, double rhouR, double rhovR, double rhowR, double ER, double &rhof, double &rhouf, double &rhovf, double &rhowf, double &Ef)
 Lax-Friedrichs Riemann solver. More...
 
 LaxFriedrichsSolver ()
 
virtual void v_PointSolve (double hL, double huL, double hvL, double hR, double huR, double hvR, double &hf, double &huf, double &hvf)
 Lax-Friedrichs Riemann solver. More...
 
- Protected Member Functions inherited from Nektar::NonlinearSWESolver
 NonlinearSWESolver ()
 
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_ArraySolve (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::SolverUtils::RiemannSolver
SOLVER_UTILS_EXPORT RiemannSolver ()
 
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...
 
- Protected Member Functions inherited from Nektar::CompressibleSolver
 CompressibleSolver ()
 
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_ArraySolve (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_PointSolveVisc (NekDouble rhoL, NekDouble rhouL, NekDouble rhovL, NekDouble rhowL, NekDouble EL, NekDouble EpsL, NekDouble rhoR, NekDouble rhouR, NekDouble rhovR, NekDouble rhowR, NekDouble ER, NekDouble EpsR, NekDouble &rhof, NekDouble &rhouf, NekDouble &rhovf, NekDouble &rhowf, NekDouble &Ef, NekDouble &Epsf)
 
- Protected Member Functions inherited from Nektar::APESolver
 APESolver ()
 
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)
 
Array< OneD, Array< OneD,
NekDouble > > 
GetRotBasefield ()
 

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::NonlinearSWESolver
bool m_pointSolve
 
- 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...
 
- Protected Attributes inherited from Nektar::CompressibleSolver
bool m_pointSolve
 

Detailed Description

Definition at line 48 of file APESolver/RiemannSolvers/LaxFriedrichsSolver.h.

Constructor & Destructor Documentation

Nektar::LaxFriedrichsSolver::LaxFriedrichsSolver ( )
protected

Definition at line 53 of file APESolver/RiemannSolvers/LaxFriedrichsSolver.cpp.

Referenced by create().

53  :
54  APESolver()
55 {
56 
57 }
Nektar::LaxFriedrichsSolver::LaxFriedrichsSolver ( )
protected
Nektar::LaxFriedrichsSolver::LaxFriedrichsSolver ( )
protected

Member Function Documentation

static RiemannSolverSharedPtr Nektar::LaxFriedrichsSolver::create ( )
inlinestatic

Definition at line 46 of file CompressibleFlowSolver/RiemannSolvers/LaxFriedrichsSolver.h.

References LaxFriedrichsSolver().

47  {
49  new LaxFriedrichsSolver());
50  }
boost::shared_ptr< RiemannSolver > RiemannSolverSharedPtr
A shared pointer to an EquationSystem object.
static RiemannSolverSharedPtr Nektar::LaxFriedrichsSolver::create ( )
inlinestatic

Definition at line 46 of file ShallowWaterSolver/RiemannSolvers/LaxFriedrichsSolver.h.

References LaxFriedrichsSolver().

47  {
49  new LaxFriedrichsSolver());
50  }
boost::shared_ptr< RiemannSolver > RiemannSolverSharedPtr
A shared pointer to an EquationSystem object.
static RiemannSolverSharedPtr Nektar::LaxFriedrichsSolver::create ( )
inlinestatic

Definition at line 51 of file APESolver/RiemannSolvers/LaxFriedrichsSolver.h.

52  {
54  }
boost::shared_ptr< RiemannSolver > RiemannSolverSharedPtr
A shared pointer to an EquationSystem object.
void Nektar::LaxFriedrichsSolver::v_PointSolve ( double  hL,
double  huL,
double  hvL,
double  hR,
double  huR,
double  hvR,
double &  hf,
double &  huf,
double &  hvf 
)
protectedvirtual

Lax-Friedrichs Riemann solver.

Parameters
hLWater depth left state.
hRWater depth right state.
huLx-momentum component left state.
huRx-momentum component right state.
hvLy-momentum component left state.
hvRy-momentum component right state.
hfComputed Riemann flux for density.
hufComputed Riemann flux for x-momentum component
hvfComputed Riemann flux for y-momentum component

Reimplemented from Nektar::NonlinearSWESolver.

Definition at line 64 of file ShallowWaterSolver/RiemannSolvers/LaxFriedrichsSolver.cpp.

References Nektar::SolverUtils::RiemannSolver::m_params, and sign.

68  {
69  static NekDouble g = m_params["gravity"]();
70 
71  // Left and right velocities
72  NekDouble uL = huL / hL;
73  NekDouble vL = hvL / hL;
74  NekDouble uR = huR / hR;
75  NekDouble vR = hvR / hR;
76 
77  // Left and right wave speeds
78  NekDouble cL = sqrt(g * hL);
79  NekDouble cR = sqrt(g * hR);
80 
81  // Square root of hL and hR.
82  NekDouble srL = sqrt(hL);
83  NekDouble srR = sqrt(hR);
84  NekDouble srLR = srL + srR;
85 
86  // Velocity Roe averages
87  NekDouble uRoe = (srL * uL + srR * uR) / srLR;
88  NekDouble cRoe = sqrt(0.5 * (cL * cL + cR * cR));
89 
90  // Minimum and maximum wave speeds
91  NekDouble S = std::max(uRoe+cRoe, std::max(uR+cR, -(uL-cL)));
92  NekDouble sign = 1.0;
93 
94  if(S == -(uL-cL))
95  {
96  sign = -1.0;
97  }
98 
99  // Lax-Friedrichs Riemann h flux
100  hf = 0.5 * ((huL + huR) - sign * S * (hR -hL));
101 
102  // Lax-Friedrichs Riemann hu flux
103  huf = 0.5 * ((hL * uL * uL + 0.5 * g * hL * hL +
104  hR * uR * uR + 0.5 * g * hR * hR) -
105  sign * S * (huR - huL));
106 
107  // Lax-Friedrichs Riemann hv flux
108  hvf = 0.5 * ((hL * uL * vL + hR * uR * vR) -
109  sign * S * (hvR - hvL));
110 
111  }
#define sign(a, b)
return the sign(b)*a
Definition: Polylib.cpp:22
double NekDouble
std::map< std::string, RSParamFuncType > m_params
Map of parameter function types.
void Nektar::LaxFriedrichsSolver::v_PointSolve ( double  rhoL,
double  rhouL,
double  rhovL,
double  rhowL,
double  EL,
double  rhoR,
double  rhouR,
double  rhovR,
double  rhowR,
double  ER,
double &  rhof,
double &  rhouf,
double &  rhovf,
double &  rhowf,
double &  Ef 
)
protectedvirtual

Lax-Friedrichs Riemann solver.

Parameters
rhoLDensity left state.
rhoRDensity right state.
rhouLx-momentum component left state.
rhouRx-momentum component right state.
rhovLy-momentum component left state.
rhovRy-momentum component right state.
rhowLz-momentum component left state.
rhowRz-momentum component right state.
ELEnergy left state.
EREnergy right state.
rhofComputed Riemann flux for density.
rhoufComputed Riemann flux for x-momentum component
rhovfComputed Riemann flux for y-momentum component
rhowfComputed Riemann flux for z-momentum component
EfComputed Riemann flux for energy.

Reimplemented from Nektar::CompressibleSolver.

Definition at line 70 of file CompressibleFlowSolver/RiemannSolvers/LaxFriedrichsSolver.cpp.

References Nektar::SolverUtils::RiemannSolver::m_params.

74  {
75  static NekDouble gamma = m_params["gamma"]();
76 
77  // Left and right velocities
78  NekDouble uL = rhouL / rhoL;
79  NekDouble vL = rhovL / rhoL;
80  NekDouble wL = rhowL / rhoL;
81  NekDouble uR = rhouR / rhoR;
82  NekDouble vR = rhovR / rhoR;
83  NekDouble wR = rhowR / rhoR;
84 
85  // Left and right pressures
86  NekDouble pL = (gamma - 1.0) *
87  (EL - 0.5 * (rhouL * uL + rhovL * vL + rhowL * wL));
88  NekDouble pR = (gamma - 1.0) *
89  (ER - 0.5 * (rhouR * uR + rhovR * vR + rhowR * wR));
90 
91  // Left and right enthalpy
92  NekDouble hL = (EL + pL) / rhoL;
93  NekDouble hR = (ER + pR) / rhoR;
94 
95  // Square root of rhoL and rhoR
96  NekDouble srL = sqrt(rhoL);
97  NekDouble srR = sqrt(rhoR);
98  NekDouble srLR = srL + srR;
99 
100  // Roe average state
101  NekDouble uRoe = (srL * uL + srR * uR) / srLR;
102  NekDouble vRoe = (srL * vL + srR * vR) / srLR;
103  NekDouble wRoe = (srL * wL + srR * wR) / srLR;
104  NekDouble hRoe = (srL * hL + srR * hR) / srLR;
105  NekDouble URoe = (uRoe * uRoe + vRoe * vRoe + wRoe * wRoe);
106  NekDouble cRoe = sqrt((gamma - 1.0)*(hRoe - 0.5 * URoe));
107 
108  // Maximum eigenvalue
109  URoe = fabs(uRoe) + cRoe;
110 
111  // Lax-Friedrichs flux formula
112  rhof = 0.5*(rhouL + rhouR - URoe*(rhoR - rhoL));
113  rhouf = 0.5*(pL + rhouL*uL + pR + rhouR*uR - URoe*(rhouR - rhouL));
114  rhovf = 0.5*(rhouL*vL + rhouR*vR - URoe*(rhovR - rhovL));
115  rhowf = 0.5*(rhouL*wL + rhouR*wR - URoe*(rhowR - rhowL));
116  Ef = 0.5*(uL*(EL + pL) + uR*(ER + pR) - URoe*(ER - EL));
117  }
double NekDouble
std::map< std::string, RSParamFuncType > m_params
Map of parameter function types.
void Nektar::LaxFriedrichsSolver::v_PointSolve ( NekDouble  pL,
NekDouble  rhoL,
NekDouble  uL,
NekDouble  vL,
NekDouble  wL,
NekDouble  pR,
NekDouble  rhoR,
NekDouble  uR,
NekDouble  vR,
NekDouble  wR,
NekDouble  p0,
NekDouble  rho0,
NekDouble  u0,
NekDouble  v0,
NekDouble  w0,
NekDouble pF,
NekDouble rhoF,
NekDouble uF,
NekDouble vF,
NekDouble wF 
)
protectedvirtual

Lax-Friedrichs Riemann solver.

Parameters
pLPerturbation pressure left state
rhoLPerturbation density left state
pRPerturbation pressure right state
rhoRPerturbation density right state
uLx perturbation velocity component left state
uRx perturbation velocity component right state
vLy perturbation velocity component left state
vRy perturbation velocity component right state
wLz perturbation velocity component left state
wRz perturbation velocity component right state
p0Base pressure
rho0Base density
u0Base x velocity component
v0Base y velocity component
w0Base z velocity component
pFComputed Riemann flux for perturbation pressure
rhoFComputed Riemann flux for perturbation density
uFComputed Riemann flux for x perturbation velocity component
vFComputed Riemann flux for y perturbation velocity component
wFComputed Riemann flux for z perturbation velocity component

Reimplemented from Nektar::APESolver.

Definition at line 83 of file APESolver/RiemannSolvers/LaxFriedrichsSolver.cpp.

References ASSERTL1, Nektar::SolverUtils::RiemannSolver::CheckParams(), and Nektar::SolverUtils::RiemannSolver::m_params.

88 {
89  ASSERTL1(CheckParams("Gamma"), "Gamma not defined.");
90  const NekDouble &gamma = m_params["Gamma"]();
91 
92  // Speed of sound
93  NekDouble c = sqrt(gamma * p0 / rho0);
94 
95  // max absolute eigenvalue of the jacobian of F_n1
96  NekDouble a_1_max = std::max(std::abs(u0 - c), std::abs(u0 + c));
97 
98  NekDouble pFL = rho0*uL + u0*pL/(c*c);
99  NekDouble uFL = pL/rho0 + u0*uL + v0*vL + w0*wL;
100  NekDouble vFL = 0;
101  NekDouble wFL = 0;
102 
103  NekDouble pFR = rho0*uR + u0*pR/(c*c);
104  NekDouble uFR = pR/rho0 + u0*uR + v0*vR + w0*wR;
105  NekDouble vFR = 0;
106  NekDouble wFR = 0;
107 
108  // assemble the face-normal fluxes
109  pF = 0.5*(pFL + pFR - a_1_max*(pR/(c*c) - pL/(c*c)));
110  uF = 0.5*(uFL + uFR - a_1_max*(uR - uL));
111  vF = 0.5*(vFL + vFR - a_1_max*(vR - vL));
112  wF = 0.5*(wFL + wFR - a_1_max*(wR - wL));
113 }
double NekDouble
std::map< std::string, RSParamFuncType > m_params
Map of parameter function types.
SOLVER_UTILS_EXPORT bool CheckParams(std::string name)
Determine whether a parameter has been defined in m_params.
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:218

Member Data Documentation

static std::string Nektar::LaxFriedrichsSolver::solverName
static
Initial value:
=
RegisterCreatorFunction("LaxFriedrichs", LaxFriedrichsSolver::create,
"Lax-Friedrichs Solver")

Definition at line 56 of file APESolver/RiemannSolvers/LaxFriedrichsSolver.h.