surface_tension.hpp

Go to the documentation of this file.

 
#ifndef PECLET_FLOW_VOF_SURFACE_TENSION_HPP
#define PECLET_FLOW_VOF_SURFACE_TENSION_HPP
 
#include <Kokkos_Core.hpp>
#include <Kokkos_MathematicalConstants.hpp>
#include <Kokkos_NumericTraits.hpp>
 
#include "vof/curvature.hpp"  // CurvatureBranch
 
namespace peclet::flow::vof {
 
KOKKOS_INLINE_FUNCTION bool csfKappaDefined(double branch) {
  const int b = static_cast<int>(branch + 0.5);
  return b >= static_cast<int>(kCurvHf) && b <= static_cast<int>(kCurvPvReduced);
}
 
KOKKOS_INLINE_FUNCTION bool csfFaceCurvature(double k0, double b0, double k1, double b1,
                                             double& kf) {
  const bool d0 = csfKappaDefined(b0), d1 = csfKappaDefined(b1);
  if (d0 && d1) {
    kf = 0.5 * (k0 + k1);
    return true;
  }
  if (d0) {
    kf = k0;
    return true;
  }
  if (d1) {
    kf = k1;
    return true;
  }
  kf = 0.0;
  return false;
}
 
KOKKOS_INLINE_FUNCTION double csfFaceForce(double sigma, double kf, double dC, double hGrad) {
  return sigma * kf * dC / hGrad;
}
 
KOKKOS_INLINE_FUNCTION double capillaryDt(double rhoSum, double h, double sigma) {
  if (!(sigma > 0.0))
    return Kokkos::Experimental::infinity_v<double>;
  return Kokkos::sqrt(rhoSum * h * h * h / (4.0 * Kokkos::numbers::pi_v<double> * sigma));
}
 
}  // namespace peclet::flow::vof
 
#endif  // PECLET_FLOW_VOF_SURFACE_TENSION_HPP