how to set up MULES for void fraction with source term??
 

Hi Foamers,
I want to solve this equation:

fvm::ddt(T)
+ fvm::div(phi, alpha)
==
S

where S is a source term. I saw the alphaEqn.H file from interPhaseChangeFoam
and I think that I have to use this correct function from CMULES.H file

template<class RhoType, class SpType, class SuType>
76 void correct
77 (
78 const RhoType& rho,
79 volScalarField& psi,
80 const surfaceScalarField& phi,
81 surfaceScalarField& phiCorr,
82 const SpType& Sp,
83 const SuType& Su,
84 const scalar psiMax,
85 const scalar psiMin
86 );

Could you tell me what should be inserted as Sp and Su in this function in my case?

Hi Foamers,
I want to solve this equation:

fvm::ddt(alpha)
+ fvm::div(phi, alpha)
==
S

where S is a source term. I saw the alphaEqn.H file from interPhaseChangeFoam
and I think that I have to use this correct function from CMULES.H file

template<class RhoType, class SpType, class SuType>
76 void correct
77 (
78 const RhoType& rho,
79 volScalarField& psi,
80 const surfaceScalarField& phi,
81 surfaceScalarField& phiCorr,
82 const SpType& Sp,
83 const SuType& Su,
84 const scalar psiMax,
85 const scalar psiMin
86 );

Could you tell me what should be inserted as Sp and Su in this function in my case?