CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   MRFSimpleFoam generates stack overflow (http://www.cfd-online.com/Forums/openfoam-solving/97905-mrfsimplefoam-generates-stack-overflow.html)

danny123 February 28, 2012 07:19

MRFSimpleFoam generates stack overflow
 
Hi,

I have problem on MRFSimpleFoam. The calculation starts just fine. After some iterations, it generates the following error:


Time = 0.435

smoothSolver: Solving for Ux, Initial residual = 0.000747915, Final residual = 8.79793e-06, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.000770294, Final residual = 9.06574e-06, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.00062951, Final residual = 7.33717e-06, No Iterations 4
GAMG: Solving for p, Initial residual = 0.00230964, Final residual = 9.97597e-05, No Iterations 3
time step continuity errors : sum local = 2.99323e-09, global = -9.6741e-22, cumulative = 8.70791e-20
ExecutionTime = 797.74 s ClockTime = 798 s

Time = 0.4375

#0 Foam::error::printStack(Foam::Ostream&) in "/home/dw/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigSegv::sigHandler(int) in "/home/dw/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::tmp<Foam::Field<Foam::outerProduct<Foam::Vec tor<double>, Foam::Vector<double> >::type> > Foam::operator*<Foam::Vector<double>, Foam::Vector<double> >(Foam::UList<Foam::Vector<double> > const&, Foam::tmp<Foam::Field<Foam::Vector<double> > > const&) in "/home/dw/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4 Foam::fv::gaussGrad<Foam::Vector<double> >::correctBoundaryConditions(Foam::GeometricField< Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>&) in "/home/dw/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 Foam::fv::gaussGrad<Foam::Vector<double> >::calcGrad(Foam::GeometricField<Foam::Vector<doub le>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const in "/home/dw/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6 Foam::fv::gradScheme<Foam::Vector<double> >::grad(Foam::GeometricField<Foam::Vector<double >, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const in "/home/dw/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#7 Foam::tmp<Foam::GeometricField<Foam::outerProduct< Foam::Vector<double>, Foam::Vector<double> >::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<Foam::Vector<double> >(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/home/dw/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#8 Foam::tmp<Foam::GeometricField<Foam::outerProduct< Foam::Vector<double>, Foam::Vector<double> >::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<Foam::Vector<double> >(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/home/dw/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#9 Foam::incompressible::RASModels::laminar::divDevRe ff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/home/dw/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#10
in "/home/dw/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/MRFSimpleFoam"
#11 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12
in "/home/dw/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/MRFSimpleFoam"
Speicherzugriffsfehler

It seems to me that error appears when OpenFoam is about to write the 0.45 time folder. The error seems reproducible. The total time is set to 2 s. My system is Ubuntu 11.4 on Intel I7 980 (6 dual core).

Is there anybody who can help?

Thanks,

Daniel

danny123 March 7, 2012 11:43

get around
 
It seems that nobody is really interested in this. Meanwhile I found a workaround. I ran the case in parallel using my 12 CPU that I have. This worked fine. Recomposing the case brought up a minor error in the 0.5 directory U field. One caracter was wrong (removed it).

So, it seems that running parallel is more stable.


All times are GMT -4. The time now is 04:43.