[npatricia]

Greetings, I am kindly seeking help on my challenge. I am using the scalarCodedSource fvOption to add a solar source with expression given as; q_solar=I(1-ρ) A_pond With surface heat transfer due to mainly radiation and evaporation, where; The thermal radiant heat transfer is given by; q_thermal=h_r A_pond (T_sky-T_pond ) h_r=4εσ((T_pond+T_sky)/2)^3 T_sky=T_a (0.55+0.61〖P_a〗^0.5 )^0.25 P_a=φexp(18.403-(3885/(T_a-43.15))) Heat transfer due to evaporation modeled as; q_evap=h_fg A_pond m_w Where; m_w=h_d (w_air-w_surf ) , h_d=h_c/(C_p 〖Le〗^(2/3) ) ,and Le=α/D_AB I am stuck on how to write this in the program. Not so good in C++. Kindly seeking your help. Thank you. PATRICIA

[npatricia]

Greetings,

I have tried to include my equations in the codedfvOption, but keep getting errors.
Some one kindly help me figure this out please.

scalarCodedSourceCoeffs
{
selectionMode all;
fields (T);

codeInclude
#{
// Need this for Stefan-Boltzmann constant
#include "constants.H"
#};

codeCorrect
#{
#};

codeAddSup
#{
using constanthysicoChemical::sigma;

if (not isActive())
return;

// Looking up model parameters
scalar I = coeffs().lookupOrDefault<scalar>("I", 781);
scalar rhow = coeffs().lookupOrDefault<scalar>("rho", 1057);
scalar Cp = coeffs().lookupOrDefault<scalar>("Cp", 3.785);
scalar Rt = coeffs().lookupOrDefault<scalar>("Rt", 0.03);
scalar ew = coeffs().lookupOrDefault<scalar>("ew", 0.96);
scalar Ta = coeffs().lookupOrDefault<scalar>("Ta", 30);


// Getting source vector from equation matrix
scalarField& src = eqn.source();
const volScalarField& T = eqn.psi();

// Iterating through patch
word top_patch_name = coeffs().lookupOrDefault<word>("patch", "top");
const fvPatch& pp = mesh().boundary()[top_patch_name];
forAll(pp, i) {
label cell_i = pp.faceCells()[i];
scalar Apond = pp.magSf()[i];
scalar Tpond = T[cell_i];

// Heating
src[cell_i] -= I*(1-Rt)*Apond/rhow/Cp;

// Cooling
scalar Pa = 0.9*exp(18.403-(3885/(Ta-43.15)));
scalar Tsky = Ta*pow(0.55+0.61*sqrt(Pa)), 1.0/4);
scalar hr = sigma.value()*4*ew*cbrt((Tpond+Tsky)/2);
scalar qr = hr*(Tsky-Tpond);
src[cell_i] -= qr*Apond/rhow/Cp;
}

the error message is as follows;
Time = 0.1

DILUPBiCGStab: Solving for Ux, Initial residual = 1.50379e-12, Final residual = 1.50379e-12, No Iterations 0
DILUPBiCGStab: Solving for Uy, Initial residual = 2.30853e-12, Final residual = 2.30853e-12, No Iterations 0
DILUPBiCGStab: Solving for Uz, Initial residual = 1.86985e-12, Final residual = 1.86985e-12, No Iterations 0
Using dynamicCode for fvOption:: radiativeHeatExchange at line 29 in "/home/patricia/OpenFoam/OpenFOAM-5.0/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/solarpond/constant/fvOptions.radiation.scalarCodedSourceCoeffs"
Creating new library in "dynamicCode/radiativeHeatExchange/platforms/linux64GccDPInt32Opt/lib/libradiativeHeatExchange_4176b5e457e204505411c7a33 07a577e60b97f24.so"
Invoking "wmake -s libso /home/patricia/OpenFoam/OpenFOAM-5.0/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/solarpond/dynamicCode/radiativeHeatExchange"
wmake libso /home/patricia/OpenFoam/OpenFOAM-5.0/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/solarpond/dynamicCode/radiativeHeatExchange
ln: ./lnInclude
wmkdep: codedFvOptionTemplate.C
Ctoo: codedFvOptionTemplate.C
/home/patricia/OpenFoam/OpenFOAM-5.0/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/solarpond/constant/fvOptions.radiation.scalarCodedSourceCoeffs: In member function ‘virtual void Foam::fv::radiativeHeatExchangeFvOptionscalarSourc e::addSup(Foam::fvMatrix<double>&, Foam::label)’:
/home/patricia/OpenFoam/OpenFOAM-5.0/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/solarpond/constant/fvOptions.radiation.scalarCodedSourceCoeffs:75:51: error: no matching function for call to ‘pow(double)’
In file included from /opt/openfoam5/src/OpenFOAM/lnInclude/fieldTypes.H:33:0,
from /opt/openfoam5/src/finiteVolume/lnInclude/fvMatricesFwd.H:32,
from /opt/openfoam5/src/finiteVolume/lnInclude/fvOption.H:47,
from /opt/openfoam5/src/fvOptions/lnInclude/cellSetOption.H:54,
from codedFvOptionTemplate.H:148,
from codedFvOptionTemplate.C:26:
/opt/openfoam5/src/OpenFOAM/lnInclude/label.H:73:7: note: candidate: Foam::label Foampow(Foam::label, Foam::label)
label pow(label a, label b);
^
/opt/openfoam5/src/OpenFOAM/lnInclude/label.H:73:7: note: candidate expects 2 arguments, 1 provided
In file included from /opt/openfoam5/src/OpenFOAM/lnInclude/floatScalar.H:38:0,
from /opt/openfoam5/src/OpenFOAM/lnInclude/scalar.H:39,
from /opt/openfoam5/src/OpenFOAM/lnInclude/fieldTypes.H:34,
from /opt/openfoam5/src/finiteVolume/lnInclude/fvMatricesFwd.H:32,
from /opt/openfoam5/src/finiteVolume/lnInclude/fvOption.H:47,
from /opt/openfoam5/src/fvOptions/lnInclude/cellSetOption.H:54,
from codedFvOptionTemplate.H:148,
from codedFvOptionTemplate.C:26:
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate: double Foampow(double, double)
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:78:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(double, double, double)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate expects 2 arguments, 1 provided
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:78:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(double, double, double)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate: double Foampow(double, float)
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:79:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(double, double, float)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate expects 2 arguments, 1 provided
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:79:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(double, double, float)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate: double Foampow(float, double)
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:80:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(double, float, double)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate expects 2 arguments, 1 provided
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:80:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(double, float, double)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate: double Foampow(float, float)
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:81:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(float, float, float)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate expects 2 arguments, 1 provided
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:81:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(float, float, float)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate: double Foampow(double, int)
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:82:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(double, double, int)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate expects 2 arguments, 1 provided
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:82:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(double, double, int)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate: double Foampow(int, double)
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:83:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(double, int, double)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate expects 2 arguments, 1 provided
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:83:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(double, int, double)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate: double Foampow(double, long int)
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:84:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(double, double, long)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate expects 2 arguments, 1 provided
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:84:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(double, double, long)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate: double Foampow(long int, double)
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:85:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(double, long, double)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate expects 2 arguments, 1 provided
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:85:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(double, long, double)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate: double Foampow(float, int)
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:86:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(float, float, int)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate expects 2 arguments, 1 provided
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:86:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(float, float, int)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate: double Foampow(int, float)
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:87:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(float, int, float)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate expects 2 arguments, 1 provided
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:87:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(float, int, float)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate: double Foampow(float, long int)
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:88:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(float, float, long)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate expects 2 arguments, 1 provided
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:88:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(float, float, long)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate: double Foampow(long int, float)
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:89:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(float, long, float)
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate expects 2 arguments, 1 provided
inline double pow(const type1 s, const type2 e) \
^
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:89:1: note: in expansion of macro ‘MAXMINPOW’
MAXMINPOW(float, long, float)
^
In file included from /opt/openfoam5/src/OpenFOAM/lnInclude/IOstream.H:48:0,
from /opt/openfoam5/src/OpenFOAM/lnInclude/Ostream.H:39,
from /opt/openfoam5/src/OpenFOAM/lnInclude/OSstream.H:39,
from /opt/openfoam5/src/OpenFOAM/lnInclude/messageStream.H:216,
from /opt/openfoam5/src/OpenFOAM/lnInclude/error.H:51,
from /opt/openfoam5/src/OpenFOAM/lnInclude/VectorSpaceI.H:26,
from /opt/openfoam5/src/OpenFOAM/lnInclude/VectorSpace.H:230,
from /opt/openfoam5/src/OpenFOAM/lnInclude/Vector.H:44,
from /opt/openfoam5/src/OpenFOAM/lnInclude/vector.H:39,
from /opt/openfoam5/src/OpenFOAM/lnInclude/fieldTypes.H:35,
from /opt/openfoam5/src/finiteVolume/lnInclude/fvMatricesFwd.H:32,
from /opt/openfoam5/src/finiteVolume/lnInclude/fvOption.H:47,
from /opt/openfoam5/src/fvOptions/lnInclude/cellSetOption.H:54,
from codedFvOptionTemplate.H:148,
from codedFvOptionTemplate.C:26:
/opt/openfoam5/src/OpenFOAM/lnInclude/uLabel.H:72:8: note: candidate: Foam::uLabel Foampow(Foam::uLabel, Foam::uLabel)
uLabel pow(uLabel a, uLabel b);
^
/opt/openfoam5/src/OpenFOAM/lnInclude/uLabel.H:72:8: note: candidate expects 2 arguments, 1 provided
In file included from /opt/openfoam5/src/OpenFOAM/lnInclude/VectorSpace.H:230:0,
from /opt/openfoam5/src/OpenFOAM/lnInclude/Vector.H:44,
from /opt/openfoam5/src/OpenFOAM/lnInclude/vector.H:39,
from /opt/openfoam5/src/OpenFOAM/lnInclude/fieldTypes.H:35,
from /opt/openfoam5/src/finiteVolume/lnInclude/fvMatricesFwd.H:32,
from /opt/openfoam5/src/finiteVolume/lnInclude/fvOption.H:47,
from /opt/openfoam5/src/fvOptions/lnInclude/cellSetOption.H:54,
from codedFvOptionTemplate.H:148,
from codedFvOptionTemplate.C:26:
/opt/openfoam5/src/OpenFOAM/lnInclude/VectorSpaceI.H:373:43: note: candidate: template<class Form, class Cmpt, unsigned char Ncmpts> typename FoampowProduct<Form, 0u>::type Foampow(const Foam::VectorSpace<Form, Cmpt, Ncmpts>&, typename FoampowProduct<Form, 0u>::type)
inline typename powProduct<Form, 0>::type pow





Thank you
PATRICIA

[tomf]

Hi Patricia,

Next time please use the Code tags from the Advanced editing menu.

These lines (and many others like them) caught my attention:

/opt/openfoam5/src/OpenFOAM/lnInclude/label.H:73:7: note: candidate expects 2 arguments, 1 provided
/opt/openfoam5/src/OpenFOAM/lnInclude/doubleFloat.H:72:15: note: candidate: double Foampow(double, double)
inline double pow(const type1 s, const type2 e)

It basically states that the pow function you used had the wrong syntax. Closer inspection yields that you have a incorrectly located bracket.

scalar Tsky = Ta*pow(0.55+0.61*sqrt(Pa)), 1.0/4);

I guess it will work if you change this to:

scalar Tsky = Ta*pow(0.55+0.61*sqrt(Pa), 1.0/4);

Hope this helps.
Regards,
Tom

[npatricia]

Thank you so much sir @Tom.
My apologies for the long lines, was not familiar with the code tags.
Like you advised, the error was with the brackets, and was able to correct it. However am getting a new error as follows;
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::PBiCGStab::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5 ? at ??:?
#6 ? at ??:?
#7 ? at ??:?
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9 ? at ??:?
Floating point exception (core dumped)

I don't get what this means and how to resolve it.
Kindly adivise.

Thank you
PATRICIA

[npatricia]

Part of the log file with the error is as follows
Starting time loop

Time = 0.1

DILUPBiCGStab: Solving for Ux, Initial residual = 1.50379e-12, Final residual = 1.50379e-12, No Iterations 0
DILUPBiCGStab: Solving for Uy, Initial residual = 2.30853e-12, Final residual = 2.30853e-12, No Iterations 0
DILUPBiCGStab: Solving for Uz, Initial residual = 1.86985e-12, Final residual = 1.86985e-12, No Iterations 0
Using dynamicCode for fvOption:: radiativeHeatExchange at line 29 in "/home/patricia/OpenFoam/OpenFOAM-5.0/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/solarpond/constant/fvOptions.radiation.scalarCodedSourceCoeffs"
Creating new library in "dynamicCode/radiativeHeatExchange/platforms/linux64GccDPInt32Opt/lib/libradiativeHeatExchange_5808c7f10c6ee6371f8d633ea d4b61675b5788db.so"
Invoking "wmake -s libso /home/patricia/OpenFoam/OpenFOAM-5.0/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/solarpond/dynamicCode/radiativeHeatExchange"
wmake libso /home/patricia/OpenFoam/OpenFOAM-5.0/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/solarpond/dynamicCode/radiativeHeatExchange
ln: ./lnInclude
wmkdep: codedFvOptionTemplate.C
Ctoo: codedFvOptionTemplate.C
ld: /home/patricia/OpenFoam/OpenFOAM-5.0/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/solarpond/dynamicCode/radiativeHeatExchange/../platforms/linux64GccDPInt32Opt/lib/libradiativeHeatExchange_5808c7f10c6ee6371f8d633ea d4b61675b5788db.so
Selecting finite volume options model type radiativeHeatExchange
Source: radiativeHeatExchange
- selecting all cells
- selected 4000 cell(s) with volume 4
DILUPBiCGStab: Solving for T, Initial residual = 1, Final residual = 0.0414164, No Iterations 1
DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 0.00369374, No Iterations 7
time step continuity errors : sum local = 1.47713e+17, global = -5.975, cumulative = -5.975
ExecutionTime = 0.05 s ClockTime = 12 s

Time = 0.2

DILUPBiCGStab: Solving for Ux, Initial residual = 0.142672, Final residual = 0.00310035, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 0.050394, Final residual = 0.000837148, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 0.142672, Final residual = 0.00310035, No Iterations 1
DILUPBiCGStab: Solving for T, Initial residual = 0.205486, Final residual = 0.0142236, No Iterations 1
DICPCG: Solving for p_rgh, Initial residual = 0.999951, Final residual = 6.55022, No Iterations 1001
time step continuity errors : sum local = 1.38877e+18, global = -80.6, cumulative = -86.575
ExecutionTime = 0.14 s ClockTime = 12 s

Time = 0.3

DILUPBiCGStab: Solving for Ux, Initial residual = 0.24404, Final residual = 7.16863e-05, No Iterations 13
DILUPBiCGStab: Solving for Uy, Initial residual = 0.197246, Final residual = 4.50391e-06, No Iterations 13
DILUPBiCGStab: Solving for Uz, Initial residual = 0.24404, Final residual = 2.1367e-05, No Iterations 13
DILUPBiCGStab: Solving for T, Initial residual = 0.297251, Final residual = 0.015471, No Iterations 1
DICPCG: Solving for p_rgh, Initial residual = 0.959089, Final residual = 28022, No Iterations 1001
time step continuity errors : sum local = 1.69267e+42, global = 1.62361e+26, cumulative = 1.62361e+26
ExecutionTime = 0.23 s ClockTime = 12 s

Time = 0.4

DILUPBiCGStab: Solving for Ux, Initial residual = 0.252988, Final residual = 0.000193031, No Iterations 14
DILUPBiCGStab: Solving for Uy, Initial residual = 0.47049, Final residual = 0.00642101, No Iterations 10
DILUPBiCGStab: Solving for Uz, Initial residual = 0.252988, Final residual = 0.000260956, No Iterations 13
DILUPBiCGStab: Solving for T, Initial residual = 0.323635, Final residual = 0.0114384, No Iterations 1
DICPCG: Solving for p_rgh, Initial residual = 0.943869, Final residual = 2.77328e+16, No Iterations 1001
time step continuity errors : sum local = 5.35021e+77, global = -2.30677e+61, cumulative = -2.30677e+61
ExecutionTime = 0.33 s ClockTime = 12 s

Time = 0.5

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::PBiCGStab::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5 ? at ??:?
#6 ? at ??:?
#7 ? at ??:?
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9 ? at ??:?
Floating point exception (core dumped)

[tomf]

This looks like some kind of set-up error (1001 pressure iterations: no convergence). Please have a look at this thread.

Regards,
Tom

[npatricia]

Kindly seeking advice on how to include mass transfer in my pond model with OpenFoam. I wish to have a model where density is a function of both temperature and concentration. in so doing, i have to add a concentration equation, and as well include it in the boussinesq approximation.

Kindly help me please.

Thank you
PATRICIA