sutherland.hpp

#ifndef ABLATELIBRARY_SUTHERLAND_HPP
#define ABLATELIBRARY_SUTHERLAND_HPP
 
#include <eos/eos.hpp>
#include <memory>
#include "transportModel.hpp"
 
namespace ablate::eos::transport {
 
class Sutherland : public TransportModel {
   private:
    const std::shared_ptr<eos::EOS> eos;
 
    const std::vector<TransportProperty> enabledProperties;
 
    // constant values
    inline const static PetscReal pr = 0.707;
    inline const static PetscReal muo = 1.716e-5;
    inline const static PetscReal to = 273.e+0;
    inline const static PetscReal so = 111.e+0;
    inline const static PetscReal sc = 0.707;
 
    static PetscErrorCode SutherlandConductivityFunction(const PetscReal conserved[], PetscReal* property, void* ctx);
    static PetscErrorCode SutherlandConductivityTemperatureFunction(const PetscReal conserved[], PetscReal temperature, PetscReal* property, void* ctx);
    static PetscErrorCode SutherlandViscosityFunction(const PetscReal conserved[], PetscReal* property, void* ctx);
    static PetscErrorCode SutherlandViscosityTemperatureFunction(const PetscReal conserved[], PetscReal temperature, PetscReal* property, void* ctx);
    static PetscErrorCode SutherlandDiffusivityFunction(const PetscReal conserved[], PetscReal* property, void* ctx);
    static PetscErrorCode SutherlandDiffusivityTemperatureFunction(const PetscReal conserved[], PetscReal temperature, PetscReal* property, void* ctx);
    static PetscErrorCode SutherlandZeroFunction(const PetscReal conserved[], PetscReal* property, void* ctx);
    static PetscErrorCode SutherlandZeroTemperatureFunction(const PetscReal conserved[], PetscReal temperature, PetscReal* property, void* ctx);
   public:
    explicit Sutherland(std::shared_ptr<eos::EOS> eos, const std::vector<TransportProperty>& enabledProperties = {});
 
    [[nodiscard]] ThermodynamicFunction GetTransportFunction(TransportProperty property, const std::vector<domain::Field>& fields) const override;
 
    [[nodiscard]] ThermodynamicTemperatureFunction GetTransportTemperatureFunction(TransportProperty property, const std::vector<domain::Field>& fields) const override;
};
}  // namespace ablate::eos::transport
 
#endif  // ABLATELIBRARY_SUTHERLAND_HPP