CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

fvOptions

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 17, 2018, 10:39
Default fvOptions
  #1
Member
 
PATRICIA NAKANWAGI
Join Date: May 2017
Posts: 47
Rep Power: 5
npatricia is on a distinguished road
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 is offline   Reply With Quote

Old   April 19, 2018, 17:12
Default Error message
  #2
Member
 
PATRICIA NAKANWAGI
Join Date: May 2017
Posts: 47
Rep Power: 5
npatricia is on a distinguished road
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
npatricia is offline   Reply With Quote

Old   April 20, 2018, 11:56
Default
  #3
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Breda, Netherlands
Posts: 556
Rep Power: 22
tomf will become famous soon enough
Send a message via MSN to tomf Send a message via Skype™ to 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:

Code:
/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.

Code:
scalar Tsky = Ta*pow(0.55+0.61*sqrt(Pa)), 1.0/4);
I guess it will work if you change this to:

Code:
scalar Tsky = Ta*pow(0.55+0.61*sqrt(Pa), 1.0/4);
Hope this helps.
Regards,
Tom
tomf is offline   Reply With Quote

Old   April 20, 2018, 17:51
Default
  #4
Member
 
PATRICIA NAKANWAGI
Join Date: May 2017
Posts: 47
Rep Power: 5
npatricia is on a distinguished road
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 is offline   Reply With Quote

Old   April 20, 2018, 18:01
Default
  #5
Member
 
PATRICIA NAKANWAGI
Join Date: May 2017
Posts: 47
Rep Power: 5
npatricia is on a distinguished road
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)
npatricia is offline   Reply With Quote

Old   April 23, 2018, 04:30
Default
  #6
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Breda, Netherlands
Posts: 556
Rep Power: 22
tomf will become famous soon enough
Send a message via MSN to tomf Send a message via Skype™ to tomf
This looks like some kind of set-up error (1001 pressure iterations: no convergence). Please have a look at this thread.

Regards,
Tom
tomf is offline   Reply With Quote

Old   May 23, 2018, 05:21
Default Inquiry
  #7
Member
 
PATRICIA NAKANWAGI
Join Date: May 2017
Posts: 47
Rep Power: 5
npatricia is on a distinguished road
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
npatricia is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
fvMatrix, fvOptions and SuSp: automatic implicit/explicit source-term treatment Zeppo OpenFOAM Programming & Development 6 May 5, 2020 12:30
New output variable for source term in fvoptions - without changing the solver vincent.clary OpenFOAM Programming & Development 2 June 26, 2018 05:21
Can I use fvOptions to couple a solid region and a fluid region? titanchao OpenFOAM Running, Solving & CFD 1 March 7, 2018 07:46
How to set fvOptions yurifrey OpenFOAM Pre-Processing 5 February 22, 2016 18:14
[swak4Foam] Setting BC for a passive scalar (groovy vs fvOptions) Tobi OpenFOAM Community Contributions 0 May 23, 2013 14:53


All times are GMT -4. The time now is 01:56.