Class KEpsilonModel

java.lang.Object
org.episteme.natural.physics.classical.matter.fluids.KEpsilonModel
public class KEpsilonModel extends Object
k-epsilon turbulence model for RANS simulations. Standard k-ε model equations: Dk/Dt = P_k - ε + ∇·((ν + ν_t/σ_k)∇k) Dε/Dt = C_ε1 * ε/k * P_k - C_ε2 * ε²/k + ∇·((ν + ν_t/σ_ε)∇ε) where ν_t = C_μ * k²/ε (turbulent viscosity) *

Reference:
Zeigler, B. P., Praehofer, H., invalid input: '&' Kim, T. G. (2000). Theory of Modeling and Simulation. Academic Press.

Since:
1.0
Author:
Silvere Martin-Michiellot, Gemini AI (Google DeepMind)

  • Field Summary

    Fields
    Modifier and Type
    Field
    Description
    static final double
    C_EPS1
     
    static final double
    C_EPS2
     
    static final double
    C_MU
     
    static final double
    SIGMA_EPS
     
    static final double
    SIGMA_K
     

  • Constructor Summary

    Constructors
    Constructor
    Description
    KEpsilonModel(int width, int height, double molecularViscosity)
     

  • Method Summary

    Modifier and Type
    Method
    Description
    static double[]
    estimateInletConditions(double velocity, double turbulenceIntensity, double lengthScale)
    Estimates inlet k and ε from turbulence intensity and length scale. k = 1.5 * (U * I)² ε = C_μ^(3/4) * k^(3/2) / L
    double[][]
    getEffectiveViscosity()
    Returns effective viscosity (molecular + turbulent).
    double[][]
    getEpsilon()
     
    double[][]
    getK()
     
    double[][]
    getTurbulentViscosity()
     
    void
    setInletConditions(double inletK, double inletEpsilon)
    Sets inlet turbulence conditions.
    void
    setVelocityGradients(double[][] dudx, double[][] dudy, double[][] dvdx, double[][] dvdy)
    Updates velocity gradients (call before each step).
    void
    step(double dt)
    Performs one time step of k-ε model.

    Methods inherited from class Object

    clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

  • Field Details

  • Constructor Details

    • KEpsilonModel

      public KEpsilonModel(int width, int height, double molecularViscosity)

  • Method Details

    • setInletConditions

      public void setInletConditions(double inletK, double inletEpsilon)
      Sets inlet turbulence conditions.
      Parameters:
      inletK - Inlet turbulent kinetic energy
      inletEpsilon - Inlet dissipation rate

    • estimateInletConditions

      public static double[] estimateInletConditions(double velocity, double turbulenceIntensity, double lengthScale)
      Estimates inlet k and ε from turbulence intensity and length scale. k = 1.5 * (U * I)² ε = C_μ^(3/4) * k^(3/2) / L

    • setVelocityGradients

      public void setVelocityGradients(double[][] dudx, double[][] dudy, double[][] dvdx, double[][] dvdy)
      Updates velocity gradients (call before each step).

    • step

      public void step(double dt)
      Performs one time step of k-ε model.
      Parameters:
      dt - Time step

    • getEffectiveViscosity

      public double[][] getEffectiveViscosity()
      Returns effective viscosity (molecular + turbulent).

    • getK

      public double[][] getK()

    • getEpsilon

      public double[][] getEpsilon()

    • getTurbulentViscosity

      public double[][] getTurbulentViscosity()