[siamak 60]

I made a new solver which compiles correctly and I created a configuration which works correctly with other solvers can anyone please help me on following error:


/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 6-e29811f5dff8
Exec : Mixedbonj
Date : Jun 22 2019
Time : 12:27:28
Host : "siamak-HP-ZBook-15-G3"
PID : 10086
I/O : uncollated
Case : /home/siamak/Desktop/validation/squareCavity
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


SIMPLE: Convergence criteria found
p_rgh: tolerance 1e-06
p: tolerance 1e-06
U: tolerance 1e-06
T: tolerance 1e-06
V: tolerance 1e-06
"(k|epsilon|omega)": tolerance 1e-05

Reading thermophysical properties

Reading field T

Reading field V

Reading field p_rgh

Reading field U

Reading field Teta

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Creating turbulence model

Selecting turbulence model type laminar
Selecting laminar stress model Stokes
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::dimensioned<double> Foam:perator/<double>(Foam::dimensioned<double> const&, Foam::dimensioned<double> const&) in "/home/siamak/OpenFOAM/siamak-6/platforms/linux64GccDPInt32Opt/bin/Mixedbonj"
#4 ? in "/home/siamak/OpenFOAM/siamak-6/platforms/linux64GccDPInt32Opt/bin/Mixedbonj"
#5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6 ? in "/home/siamak/OpenFOAM/siamak-6/platforms/linux64GccDPInt32Opt/bin/Mixedbonj"
Floating point exception (core dumped)


http://www.mediafire.com/file/y2s0ye...avity.zip/file http://www.mediafire.com/file/ctr57ii1opnk3lo/Mixedbonj.zip/file

first link is configuration file and the 2nd one is solver's.


thanks in advance

Siamak

[siamak 60]

I find that error comes from these lines and when I delete them no bugs occur:



// thermal capacity of nanofluid
dimensionedScalar Cpnf = ((Cv*rhoP+(1.0-Cv)*rhoRefbf)/rhonf) ;

//ratio of product of rhoCp of particle to rhoCp of nanofluid
dimensionedScalar rhoCp_ratio = (rhoP*Cpp/(rhoRefbf*Cpbf));



// Brownian diffusion coefficient of nanofluid
volScalarField DB
(
IOobject
(
"DB",
runTime.timeName(),
mesh
),
(KB*T/(3.0*3.14159265*miunf*dp))
);

// Thermal diffusion coefficient of nanofluid
volScalarField DT
(
IOobject
(
"DT",
runTime.timeName(),
mesh
),
((miunf*V*0.26*Knf)/(rhonf*rhok*(2.0*Knf+Kp)))
);



Also I find out problem is from denominator part ? all of denominator parts are dimensionedScalar elements !!!
can anyone guess something ?