[qibintj]

hello, Foamers

I am trying to modify the scalarTransportFoam solver TO solve these following equation.



# include "readSIMPLEControls.H"
volScalarField nu1=DO2*pow(DT*(1-S),Tor);
volScalarField nu2=Kp*S*S*S/Nu*(-1*Sigma)*0.5*pow((DT/Kp),0.5)*(-4.24*S+3.789*S*S);
volScalarField jc=(1-S)*Tcu*(T/CO2)*exp(-alpha*F*Eta/(R*T1));
for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
{
solve
(
//fvm::ddt(T)
fvm::div(nu1,T)
+(1/(4*F))*jc
);
solve
(
//fvm::ddt(S)
fvm::laplacian(nu2,S)
-(1/(2*F))*jc*M
);
}

}


there is always problem.

[COLOR="rgb(139, 0, 0)"]
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::diagonalPreconditioner::diagonalPrecondition er(Foam::lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam::lduMatrix:reconditioner::addsymMatrixConst ructorToTable<Foam::diagonalPreconditioner>::New(F oam::lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#5 Foam::lduMatrix:reconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#6 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam171/lib/linuxGccDPOpt/libfiniteVolume.so"
#8
in "/opt/openfoam171/applications/bin/linuxGccDPOpt/myFoam"
#9
in "/opt/openfoam171/applications/bin/linuxGccDPOpt/myFoam"
#10 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#11
in "/opt/openfoam171/applications/bin/linuxGccDPOpt/myFoam"
Floating point exception
[/COLOR]

please help me to solve this problem. very appreciated.