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

gaza · July 7, 2015, 15:23

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?

gaza · July 7, 2015, 15:28

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?

July 7, 2015, 15:23	how to set up MULES for void fraction with source term??	#1
gaza Senior Member Przemek Join Date: Jun 2011 Posts: 249 Rep Power: 16	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? __________________ best regards pblasiak

July 7, 2015, 15:28		#2
gaza Senior Member Przemek Join Date: Jun 2011 Posts: 249 Rep Power: 16	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? __________________ best regards pblasiak

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Source Term due to evaporation in energy transport equation	styleworker	OpenFOAM Programming & Development	3	September 7, 2022 04:09
[swak4Foam] funkyDoCalc with OF2.3 massflow	NiFl	OpenFOAM Community Contributions	14	November 25, 2020 04:30
[swak4Foam] swak4foam building problem	GGerber	OpenFOAM Community Contributions	54	April 24, 2015 17:02
[swak4Foam] build problem swak4Foam OF 2.2.0	mcathela	OpenFOAM Community Contributions	14	April 23, 2013 14:59
DecomposePar links against liblamso0 with OpenMPI	jens_klostermann	OpenFOAM Bugs	11	June 28, 2007 18:51