twoPhaseTransport.hpp

#ifndef ABLATELIBRARY_TWOPHASETRANSPORT_HPP
#define ABLATELIBRARY_TWOPHASETRANSPORT_HPP
 
#include "transportModel.hpp"
 
namespace ablate::eos::transport {
 
class TwoPhaseTransport : public TransportModel {
   private:
    const std::shared_ptr<eos::transport::TransportModel> alphaTransportModel;
    const std::shared_ptr<eos::transport::TransportModel> otherTransportModel;
 
    struct FunctionContext {
        // The function for the alpha property
        PetscErrorCode (*alphaFunction)(const PetscReal conserved[], PetscReal* property, void* ctx);
        void* alphaContext;
 
        // The function for the other property
        PetscErrorCode (*otherFunction)(const PetscReal conserved[], PetscReal* property, void* ctx);
        void* otherContext;
 
        // offset location for volume fraction in the solution field
        PetscInt vfOffset;
    };
 
    struct TemperatureFunctionContext {
        // The function for the alpha property
        PetscErrorCode (*alphaFunction)(const PetscReal conserved[], PetscReal T, PetscReal* property, void* ctx);
        void* alphaContext;
 
        // The function for the other property
        PetscErrorCode (*otherFunction)(const PetscReal conserved[], PetscReal T, PetscReal* property, void* ctx);
        void* otherContext;
 
        // offset location for volume fraction in the solution field
        PetscInt vfOffset;
    };
 
    mutable std::vector<ThermodynamicFunction> thermodynamicFunctionReference;
 
    mutable std::vector<ThermodynamicTemperatureFunction> thermodynamicTemperatureFunctionReference;
 
    static PetscErrorCode TwoPhaseFunction(const PetscReal conserved[], PetscReal* property, void* ctx);
    static PetscErrorCode TwoPhaseTemperatureFunction(const PetscReal conserved[], PetscReal temperature, PetscReal* property, void* ctx);
   public:
    explicit TwoPhaseTransport(std::shared_ptr<TransportModel> transportModelAlpha, std::shared_ptr<TransportModel> transportModelOther);
 
    [[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_TWOPHASETRANSPORT_HPP