collocated_varrho_point.hpp

Go to the documentation of this file.

 
#ifndef PECLET_FLOW_COLLOCATED_VARRHO_POINT_HPP
#define PECLET_FLOW_COLLOCATED_VARRHO_POINT_HPP
 
#include <Kokkos_Core.hpp>
 
namespace peclet::flow::varrho {
 
// Momentum-weighted face velocity: (rho_L u_L + rho_R u_R) / (rho_L + rho_R).
KOKKOS_INLINE_FUNCTION double faceMassWeighted(double uL, double uR, double rL, double rR) {
  return (rL * uL + rR * uR) / (rL + rR);
}
 
// The cell value reconstructed from its two faces along one axis: ½ (o_lo a_lo + o_hi a_hi).
KOKKOS_INLINE_FUNCTION double cellFromFaces(double aLo, double aHi, double oLo, double oHi) {
  return 0.5 * (oLo * aLo + oHi * aHi);
}
 
// The weight sum of cellFromFaces for the same two openness values: ½ (o_lo + o_hi).
KOKKOS_INLINE_FUNCTION double weightSum(double oLo, double oHi) {
  return 0.5 * (oLo + oHi);
}
 
}  // namespace peclet::flow::varrho
 
#endif  // PECLET_FLOW_COLLOCATED_VARRHO_POINT_HPP