thermophoreticDiffusion.hpp

#ifndef ABLATELIBRARY_THERMOPHORETICDIFFUSION_HPP
#define ABLATELIBRARY_THERMOPHORETICDIFFUSION_HPP
 
#include "eos/transport/transportModel.hpp"
#include "finiteVolume/fluxCalculator/fluxCalculator.hpp"
#include "flowProcess.hpp"
 
namespace ablate::finiteVolume::processes {
 
class ThermophoreticDiffusion : public FlowProcess {
   private:
    const std::shared_ptr<eos::transport::TransportModel> transportModel;
 
    // store the thermodynamicTemperatureFunction to compute viscosity
    eos::ThermodynamicTemperatureFunction viscosityTemperatureFunction;
 
   public:
    explicit ThermophoreticDiffusion(std::shared_ptr<eos::transport::TransportModel> transportModel);
 
    void Setup(ablate::finiteVolume::FiniteVolumeSolver& flow) override;
 
   private:
    static PetscErrorCode ThermophoreticDiffusionEnergyFlux(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 ThermophoreticDiffusionVariableFlux(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);
};
 
}  // namespace ablate::finiteVolume::processes
#endif  // ABLATELIBRARY_SPECIESTRANSPORT_HPP