[wavefunction]

Hello to everyone reading this.

I am having trouble adding a source term in MULES. I want to modify interMixingFoam to accomodate phase change. So I want to add the source terms -mDot/rho1 to alpha1Eqn and mDot/rho2 to alpha2Eqn.

mDot = D*rhogp*mag(fvc::grad(Y))*mag(fvc::grad(alpha1))/(scalar(1) - Y)

where D - diffusion coefficient, Y - mass fraction, rhogp - gas phase density.

Without adding the source terms to MULES, my case using my custom solver is diverging immensely. On running gdb, I found that the problem lies at line 442 of MULESTemplates.C.

Please help me.

Program received signal SIGFPE, Arithmetic exception.
0x00000000004900b1 in Foam::MULES::limiter<double, Foam::geometricOneField, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>, Foam::zeroField> (
allLambda=..., rDeltaT=@0x7fffffff4cb0: 400000000, rho=..., psi=...,
phiBD=..., phiCorr=..., Sp=..., Su=..., psiMax=1, psiMin=0)
at /home/saicharanb56/OpenFOAM/OpenFOAM-5.x/src/finiteVolume/lnInclude/MULESTemplates.C:442
442 max(min

[Pavithra]

Hello,

I am facing the same problem. Did you have any success in solving this issue ?

Kindly please help me.

Thank You.