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

#include <AUSM3Solver.h>

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

Static Public Member Functions

static RiemannSolverSharedPtr create ()
 

Static Public Attributes

static std::string solverName
 

Protected Member Functions

 AUSM3Solver ()
 
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)
 AUSM3 Riemann solver. More...
 
double M1Function (int A, double M)
 
double M2Function (int A, double M)
 
double M4Function (int A, double beta, double M)
 
double P5Function (int A, double alpha, double M)
 
- 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::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...
 

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::CompressibleSolver
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...
 

Detailed Description

Definition at line 43 of file AUSM3Solver.h.

Constructor & Destructor Documentation

Nektar::AUSM3Solver::AUSM3Solver ( )
protected

Definition at line 46 of file AUSM3Solver.cpp.

Referenced by create().

47  {
48 
49  }

Member Function Documentation

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

Definition at line 46 of file AUSM3Solver.h.

References AUSM3Solver().

47  {
49  new AUSM3Solver());
50  }
boost::shared_ptr< RiemannSolver > RiemannSolverSharedPtr
A shared pointer to an EquationSystem object.
double Nektar::AUSM3Solver::M1Function ( int  A,
double  M 
)
protected

Definition at line 142 of file AUSM3Solver.cpp.

Referenced by M4Function(), and P5Function().

143  {
144  double out;
145 
146  if (A == 0)
147  {
148  out = 0.5 * (M + fabs(M));
149  }
150  else
151  {
152  out = 0.5 * (M - fabs(M));
153  }
154 
155  return out;
156  }
double Nektar::AUSM3Solver::M2Function ( int  A,
double  M 
)
protected

Definition at line 158 of file AUSM3Solver.cpp.

Referenced by M4Function(), and P5Function().

159  {
160  double out;
161 
162  if (A == 0)
163  {
164  out = 0.25 * (M + 1.0) * (M + 1.0);
165  }
166  else
167  {
168  out = -0.25 * (M - 1.0) * (M - 1.0);
169  }
170 
171  return out;
172  }
double Nektar::AUSM3Solver::M4Function ( int  A,
double  beta,
double  M 
)
protected

Definition at line 174 of file AUSM3Solver.cpp.

References M1Function(), and M2Function().

Referenced by v_PointSolve().

175  {
176  double out;
177 
178  if (fabs(M) >= 1.0)
179  {
180  out = M1Function(A, M);
181  }
182  else
183  {
184  out = M2Function(A, M);
185 
186  if (A == 0)
187  {
188  out *= 1.0 - 16.0 * beta * M2Function(1, M);
189  }
190  else
191  {
192  out *= 1.0 + 16.0 * beta * M2Function(0, M);
193  }
194  }
195 
196  return out;
197  }
double M1Function(int A, double M)
double M2Function(int A, double M)
double Nektar::AUSM3Solver::P5Function ( int  A,
double  alpha,
double  M 
)
protected

Definition at line 199 of file AUSM3Solver.cpp.

References M1Function(), and M2Function().

Referenced by v_PointSolve().

200  {
201  double out;
202 
203  if (fabs(M) >= 1.0)
204  {
205  out = (1.0 / M) * M1Function(A, M);
206  }
207  else
208  {
209  out = M2Function(A, M);
210 
211  if (A == 0)
212  {
213  out *= (2.0 - M) - 16.0 * alpha * M * M2Function(1, M);
214  }
215  else
216  {
217  out *= (-2.0 - M) + 16.0 * alpha * M * M2Function(0, M);
218  }
219  }
220 
221  return out;
222  }
double M1Function(int A, double M)
double M2Function(int A, double M)
void Nektar::AUSM3Solver::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

AUSM3 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 AUSM3Solver.cpp.

References M4Function(), Nektar::SolverUtils::RiemannSolver::m_params, and P5Function().

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  NekDouble cL = sqrt(gamma * pL / rhoL);
91  NekDouble cR = sqrt(gamma * pR / rhoR);
92 
93  // Average speeds of sound
94  NekDouble cA = 0.5 * (cL + cR);
95 
96  // Local Mach numbers
97  NekDouble ML = uL / cA;
98  NekDouble MR = uR / cA;
99 
100  // Parameters for specify the upwinding
101  // Note: if fa = 1 then AUSM3 = AUSM3
102  NekDouble Mco = 0.01;
103  NekDouble Mtilde = 0.5 * (ML * ML + MR * MR);
104  NekDouble Mo = std::min(1.0, std::max(Mtilde, Mco*Mco));
105  NekDouble fa = Mo * (2.0 - Mo);
106  NekDouble beta = 0.125;
107  NekDouble alpha = 0.1875;
108  NekDouble sigma = 1.0;
109  NekDouble Kp = 0.25;
110  NekDouble Ku = 0.75;
111  NekDouble rhoA = 0.5 * (rhoL + rhoR);
112  NekDouble Mp = -(Kp / fa) * ((pR - pL) / (rhoA * cA * cA)) *
113  std::max(1.0 - sigma * Mtilde, 0.0);
114 
115  NekDouble Mbar = M4Function(0, beta, ML) +
116  M4Function(1, beta, MR) + Mp;
117 
118  NekDouble pu = -2.0 * Ku * rhoA * cA * cA * (MR - ML) *
119  P5Function(0, alpha, ML) * P5Function(1, alpha, MR);
120 
121  NekDouble pbar = pL * P5Function(0, alpha, ML) +
122  pR * P5Function(1, alpha, MR) + pu;
123 
124  if (Mbar >= 0.0)
125  {
126  rhof = cA * Mbar * rhoL;
127  rhouf = cA * Mbar * rhoL * uL + pbar;
128  rhovf = cA * Mbar * rhoL * vL;
129  rhowf = cA * Mbar * rhoL * wL;
130  Ef = cA * Mbar * (EL + pL);
131  }
132  else
133  {
134  rhof = cA * Mbar * rhoR;
135  rhouf = cA * Mbar * rhoR * uR + pbar;
136  rhovf = cA * Mbar * rhoR * vR;
137  rhowf = cA * Mbar * rhoR * wR;
138  Ef = cA * Mbar * (ER + pR);
139  }
140  }
double P5Function(int A, double alpha, double M)
double NekDouble
double M4Function(int A, double beta, double M)
std::map< std::string, RSParamFuncType > m_params
Map of parameter function types.

Member Data Documentation

std::string Nektar::AUSM3Solver::solverName
static
Initial value:

Definition at line 52 of file AUSM3Solver.h.