les.hpp

#ifndef ABLATELIBRARY_LES_H
#define ABLATELIBRARY_LES_H
 
#include "eos/transport/transportModel.hpp"
#include "finiteVolume/fluxCalculator/fluxCalculator.hpp"
#include "flowProcess.hpp"
#include "navierStokesTransport.hpp"
 
namespace ablate::finiteVolume::processes {
 
class LES : public FlowProcess {
   private:
    const std::string tkeField;
 
    // constant values
    inline const static PetscReal c_k = 0.094;
    inline const static PetscReal c_e = 1.048;
    inline const static PetscReal c_p = 1004.0;
    inline const static PetscReal scT = 1.00;
    inline const static PetscReal prT = 1.00;
 
    // Store the numberComponents for each ev/species
    std::vector<PetscInt> numberComponents;
 
   public:
    explicit LES(std::string tkeField = {});
 
    void Setup(ablate::finiteVolume::FiniteVolumeSolver& flow) override;
 
   public:
    static PetscErrorCode LesMomentumFlux(PetscInt dim, const PetscFVFaceGeom* fg, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar field[], const PetscScalar grad[],
                                          const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar aux[], const PetscScalar gradAux[], PetscScalar flux[], void* ctx);
    static PetscErrorCode LesEnergyFlux(PetscInt dim, const PetscFVFaceGeom* fg, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar field[], const PetscScalar grad[],
                                        const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar aux[], const PetscScalar gradAux[], PetscScalar flux[], void* ctx);
 
    static PetscErrorCode LesTkeFlux(PetscInt dim, const PetscFVFaceGeom* fg, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar field[], const PetscScalar grad[],
                                     const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar aux[], const PetscScalar gradAux[], PetscScalar flux[], void* ctx);
 
    static PetscErrorCode LesEvFlux(PetscInt dim, const PetscFVFaceGeom* fg, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar field[], const PetscScalar grad[], const PetscInt aOff[],
                                    const PetscInt aOff_x[], const PetscScalar aux[], const PetscScalar gradAux[], PetscScalar flux[], void* ctx);
 
    static PetscErrorCode LesViscosity(PetscInt dim, const PetscFVFaceGeom* fg, const PetscScalar* densityField, const PetscReal turbulence, PetscReal& mut);
};
}  // namespace ablate::finiteVolume::processes
 
#endif  // ABLATELIBRARY_LES_H