wmpearsonIVfunctions.h

Go to the documentation of this file.

#ifndef WMPEARSONIVFUNCTIONS_H
#define WMPEARSONIVFUNCTIONS_H
 
namespace WmPearsonIVFunctions
{
enum
{
    ixx = 4,
    ixy = 5,
    ixz = 6,
    iyy = 7,
    iyz = 8,
    izz = 9
};
 
template<typename real> inline void PrimitiveToConserved(const real* qnc, real* q)
{
    // Convert non-conserved (density, velocity, pressure, heat flux) to conserved
    // variables (density, momentum, energy, conserved heat flux) Density
    q[0] = qnc[0];
 
    // Velocity
    q[1] = qnc[1] * qnc[0];
    q[2] = qnc[2] * qnc[0];
    q[3] = qnc[3] * qnc[0];
 
    // Energy
    q[ixx] = qnc[ixx] + qnc[1] * q[1];
    q[ixy] = qnc[ixy] + qnc[1] * q[2];
    q[ixz] = qnc[ixz] + qnc[1] * q[3];
    q[iyy] = qnc[iyy] + qnc[2] * q[2];
    q[iyz] = qnc[iyz] + qnc[2] * q[3];
    q[izz] = qnc[izz] + qnc[3] * q[3];
 
    // (Conserved) Heat flux
    const real ejj = q[ixx] + q[iyy] + q[izz];
    q[10] = 2. * (qnc[10] + qnc[ixx] * qnc[1] + qnc[ixy] * qnc[2] + qnc[ixz] * qnc[3]) +
            ejj * qnc[1];
    q[11] = 2. * (qnc[11] + qnc[ixy] * qnc[1] + qnc[iyy] * qnc[2] + qnc[iyz] * qnc[3]) +
            ejj * qnc[2];
    q[12] = 2. * (qnc[12] + qnc[ixz] * qnc[1] + qnc[iyz] * qnc[2] + qnc[izz] * qnc[3]) +
            ejj * qnc[3];
}
 
template<typename real>
inline void InitializeIsotropic(const real density,
                                const real velocity[3],
                                const real pressure,
                                real q[13])
{
    // Initializes 13 conserved moments from density, velocity and pressure
    real qnc[13] = {0};
    // Density
    qnc[0] = density;
 
    // Velocity
    qnc[1] = velocity[0];
    qnc[2] = velocity[1];
    qnc[3] = velocity[2];
 
    // Pressure
    qnc[ixx] = pressure;
    qnc[ixy] = 0.0;
    qnc[ixz] = 0.0;
    qnc[iyy] = pressure;
    qnc[iyz] = 0.0;
    qnc[izz] = pressure;
 
    // Heat flux
    qnc[10] = 0.0;
    qnc[11] = 0.0;
    qnc[12] = 0.0;
 
    // Convert to conserved quantities
    PrimitiveToConserved(qnc, q);
}
 
} // namespace WmPearsonIVFunctions
 
#endif // WMPEARSONIVFUNCTIONS_H