[adambarfi]

Originally Posted by slash89

Hi all,
i got the same problem. I am using the buoyantSimpleRadiationFoam. The problems is always at the second time step, when solving the G file. Any suggestions to fix this problem?

Thank you,

Best regards

Greetings slash89,

please attach your log file, in this way we can help you much more easily

[slash89]

Here is the log file. I am sorry but i could not upload the file and i don't know why!

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : buoyantSimpleRadiationFoam -parallel
Date : Dec 05 2014
Time : 18:58:35
Host : "millegradi-nb"
PID : 4112
Case : /home/bolzo/TERMIGNONI/run/prove_solver/prova_5
nProcs : 4
Slaves :
3
(
"millegradi-nb.4113"
"millegradi-nb.4114"
"millegradi-nb.4115"
)

Pstream initialized with:
floatTransfer : 0
nProcsSimpleSum : 0
commsType : nonBlocking
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

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

Create mesh for time = 0


Reading g
Reading thermophysical properties

Selecting thermodynamics package hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
sigmak 1;
sigmaEps 1.3;
Prt 1;
}

Calculating field g.h

Reading field p_rgh

Selecting radiationModel P1
Selecting absorptionEmissionModel constantAbsorptionEmission
Selecting scatterModel constantScatter

SIMPLE: convergence criteria
field p_rgh tolerance 0.01
field U tolerance 0.001
field h tolerance 0.001
field G tolerance 0.001
field "(k|epsilon|omega)" tolerance 0.001


Starting time loop

Time = 1

DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.00766817, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.00749132, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 0.00791072, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 0.00766938, No Iterations 1
DICPCG: Solving for G, Initial residual = 1, Final residual = 0.096472, No Iterations 335
DICPCG: Solving for p_rgh, Initial residual = 0.999948, Final residual = 0.00829887, No Iterations 450
time step continuity errors : sum local = 0.145697, global = -0.0018876, cumulative = -0.0018876
rho max/min : 79.6338 1.1739
DILUPBiCG: Solving for epsilon, Initial residual = 0.12008, Final residual = 0.00390459, No Iterations 1
bounding epsilon, min: -5.64769 max: 1881.23 average: 13.0219
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 0.0865008, No Iterations 1
ExecutionTime = 23.58 s ClockTime = 23 s

Time = 2

DILUPBiCG: Solving for Ux, Initial residual = 0.998573, Final residual = 0.0137559, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.999799, Final residual = 0.0141418, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.999881, Final residual = 0.0139621, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 0.013323, No Iterations 1

Thak you!

[adambarfi]

Are you sure about the BCs you exerted on your geometry? and I wanted the log file not the fist two iteration log

it seems that your cumulative time step error will grow up and the range of rho variation I think is very much!!!

time step continuity errors : sum local = 0.145697, global = -0.0018876, cumulative = -0.0018876
rho max/min : 79.6338 1.1739

if you're sure about the BCs, then you should check your schemes.

[slash89]

Sorry but i don't understand which file you need. Do you need the G file?

Best regards

[adambarfi]

just run your case and get the log file:

solverName >log

and the errors appeared in your terminal.
such as: http://www.cfd-online.com/Forums/ope...tml#post375086

[slash89]

--> FOAM FATAL ERROR: 
[2] Maximum number of iterations exceeded
[2] 
[2]     From function specieThermo<Thermo>::T(scalar f, scalar T0, scalar (specieThermo<Thermo>::*F)(const scalar) const, scalar (specieThermo<Thermo>::*dFdT)(const scalar) const) const
[2]     in file /home/opencfd/OpenFOAM/OpenFOAM-2.1.1/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 69.
[2] 
FOAM parallel run aborting
[2] 
[2] #0  Foam::error::printStack(Foam::Ostream&)[0] 
[0] 
[0] --> FOAM FATAL ERROR: 
[0] Maximum number of iterations exceeded
[0] 
[0]     From function specieThermo<Thermo>::T(scalar f, scalar T0, scalar (specieThermo<Thermo>::*F)(const scalar) const, scalar (specieThermo<Thermo>::*dFdT)(const scalar) const) const
[0]     in file /home/opencfd/OpenFOAM/OpenFOAM-2.1.1/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 69.
[0] 
FOAM parallel run aborting
[0] 
[0] #0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[2] #1  Foam::error::abort() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #1  Foam::error::abort() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[2] #2  Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::T(double, double, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #2  Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::T(double, double, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
[2] #3  Foam::hPsiThermo<Foam::pureMixture<Foam::constTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::calculate() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
[0] #3  Foam::hPsiThermo<Foam::pureMixture<Foam::constTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::calculate() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
[0] #4   in "/opt/openfoam2Foam::hPsiThermo<Foam::pureMixture<Foam::constTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::correct()11/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
[2] #4  Foam::hPsiThermo<Foam::pureMixture<Foam::constTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::correct() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
[0] #5   in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
[2] #5  

[0]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantSimpleRadiationFoam"
[0] #6  __libc_start_main[2]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantSimpleRadiationFoam"
[2] #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #7   in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #7  

[0]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantSimpleRadiationFoam"
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD 
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
[2]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantSimpleRadiationFoam"
--------------------------------------------------------------------------
mpirun has exited due to process rank 0 with PID 2812 on
node millegradi-nb exiting without calling "finalize". This may
have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).
--------------------------------------------------------------------------
[millegradi-nb:02811] 1 more process has sent help message help-mpi-api.txt / mpi-abort
[millegradi-nb:02811] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

Hope this is what we need!

[adambarfi]

The problem is probably related to the pressure field. I guess you have negative pressure, so the thermophysical library crashes.

check the following:

1. check your grids, checkmesh
2. use 1st order scheme for temporal discretization
3. modify the tolerance tol_ from 1.0e-4 to something higher e.g. 1.0e-3, but this is not a good advice. thanks to dmoroian
4. at last, modify the maxIter_ from 100 to something larger. you can find the implementation process in http://www.cfd-online.com/Forums/ope...tml#post179437

hope these help you

[slash89]

It still does not work. At this point the problem is in the BCs. Can I upload them here?

[adambarfi]

Originally Posted by slash89

It still does not work. At this point the problem is in the BCs. Can I upload them here?

it's the best thing you can do!

[slash89]

0.tar.gz


Thak you!

[slash89]

nobody can help??

[masoudsh]

hi

i have this problem
can anyone help me?
if i change the mesh ,solve?
or if i remove the energy equation

best regards

masoud

--> FOAM FATAL ERROR:
Maximum number of iterations exceeded

From function specieThermo<Thermo>::T(scalar f, scalar T0, scalar (specieThermo<Thermo>::*F)(const scalar) const, scalar (specieThermo<Thermo>::*dFdT)(const scalar) const) const
in file /home/opencfd/OpenFOAM/OpenFOAM-2.0.1/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 67.

[derekm]

Need a lot more info

[masoudsh]

hi Derek

what do you want to know?
did you see this problem later?

best regards
masoud

[masoudsh]

Hi

I got the sam problem,Can anyone solve it?
I have this problem in my project,I do anything such as mesh,Bc , ... but doesn't work.
if I find anything tell here ,please help me if the problem solve

Best Regards
Masoud

[Mojtaba.a]

Originally Posted by masoudsh

Hi

I got the sam problem,Can anyone solve it?
I have this problem in my project,I do anything such as mesh,Bc , ... but doesn't work.
if I find anything tell here ,please help me if the problem solve

Best Regards
Masoud

Hi Masoud,
Take a look at this:
http://openfoam.ir/questions/questio...AE%D8%B7%D8%A7

Best.