[Zi_Aun]

Hello,
Trying to derive the totalPressure boundary condition for a specific use, I found that the conditions for the selection between transonic and supersonic regimes for High-speed compressible flow are not consistent with the aerodynamic definition of these regimes. These conditions are associated with \gamma (specific heat ratios), where they should depend on the speed at least (Mach number).
See the script here below from: fields/fvPatchFields/derived/totalPressure/totalPressureFvPatchScalarField.C (OpenFOAM-5.x)

################################################## ##########################

        else
          {
              // High-speed compressible flow
   
              const fvPatchField<scalar>& psip =
                  patch().lookupPatchField<volScalarField, scalar>(psiName_);
   
              if (gamma_ > 1)
              {
                  scalar gM1ByG = (gamma_ - 1)/gamma_;
   
                  operator==
                  (
                      p0p
                     /pow
                      (
                          (1.0 + 0.5*psip*gM1ByG*(1.0 - pos0(phip))*magSqr(Up)),
                          1.0/gM1ByG
                      )
                  );
              }
              else
              {
                  operator==(p0p/(1.0 + 0.5*psip*(1.0 - pos0(phip))*magSqr(Up)));
              }
          }

################################################## #########################
This may induce errors, any comment or suggestion regarding the way to proceed with this?
Thanks.
Zi Aun.