|
[Sponsors] |
September 7, 2017, 20:13 |
Dam Break sigFpe Error
|
#1 |
Member
Andrew O. Winter
Join Date: Aug 2015
Location: Seattle, WA, USA
Posts: 78
Rep Power: 11 |
Hello CFD-Online Community,
I'm trying to update to OpenFOAM-5.0 from OpenFOAM-2.2.1 and want to re-validate my case setup using data from a dam break experiment performed in a wave flume with a bridge deck to measure the forces and pressures applied to the structure. I'm using the interFoam solver and running this case in parallel using 8 8-core machines in the University of Washington's HPC network, Hyak, which uses Red Hat Linux (kernel version 2.6.32-696.3.2.el6.x86_64). OpenFOAM-5.0 was compiled using gcc 5.1.0 with openmpi 1.10.1. The initial state of this setup was reproduced successfully in OpenFOAM as shown below. Initial state (domain is truncated in the image) As far as I can tell, I've updated all of my case files adequately; however, it appears that I'm having issues with the solution blowing up when the wave is just about to impact the bridge deck as shown below. Blow up state (domain is truncated in the image) Does anyone here have any clues as to why this may be occurring? I'm wondering if the solution schemes changing from one OpenFOAM version to another is hurting me somehow. I used the interFoam dambreak tutorial setup as an example for how to modify my old case to work with OpenFOAM-5.0, but am not noticing any red flags in the scripts. I've attached a .zip of the case files prior to being run as well as directly posting the log file, error output, fvSchemes, and fvSolution, which seemed like the most relevant things to post in full here. Thanks in advance! Andrew Winter Log file (I omitted the beginning where it seemed things were going fine and the pressure force in the x-direction was less than 1 N) Code:
PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 2.87216e-08, Final residual = 1.64545e-14, No Iterations 1 Phase-1 volume fraction = 0.305209 Min(alpha.water) = -6.89942e-06 Max(alpha.water) = 1.00001 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.305209 Min(alpha.water) = -6.89936e-06 Max(alpha.water) = 1.00001 DICPCG: Solving for p_rgh, Initial residual = 0.0183283, Final residual = 0.000705997, No Iterations 3 time step continuity errors : sum local = 9.64653e-12, global = 4.09114e-15, cumulative = -8.99631e-09 DICPCG: Solving for p_rgh, Initial residual = 0.0138203, Final residual = 0.000603006, No Iterations 6 time step continuity errors : sum local = 8.44461e-12, global = 3.59057e-15, cumulative = -8.9963e-09 DICPCG: Solving for p_rgh, Initial residual = 0.0362388, Final residual = 9.56855e-08, No Iterations 552 time step continuity errors : sum local = 1.40785e-15, global = 5.722e-17, cumulative = -8.9963e-09 smoothSolver: Solving for epsilon, Initial residual = 4.61478e-07, Final residual = 4.14502e-09, No Iterations 1 smoothSolver: Solving for k, Initial residual = 2.64957e-07, Final residual = 4.28379e-08, No Iterations 1 ExecutionTime = 11958.4 s ClockTime = 11981 s forces bridgeLoads write: sum of forces: pressure : (232.409 4.84786 0.850411) viscous : (3.50913 1.77603 -0.00446115) porous : (0 0 0) sum of moments: pressure : (1.66139 -113.588 1.14428) viscous : (0.884317 -1.75135 0.135472) porous : (0 0 0) Courant Number mean: 4.89028e-06 max: 25.0716 Interface Courant Number mean: 2.32166e-07 max: 0.000914802 deltaT = 3.44935e-09 Time = 3.18936 PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 5.72776e-10, Final residual = 5.72776e-10, No Iterations 0 Phase-1 volume fraction = 0.305209 Min(alpha.water) = -6.89936e-06 Max(alpha.water) = 1.00001 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.305209 Min(alpha.water) = -6.89936e-06 Max(alpha.water) = 1.00001 DICPCG: Solving for p_rgh, Initial residual = 0.983134, Final residual = 0.0430521, No Iterations 7 time step continuity errors : sum local = 1.35661e-11, global = -8.2689e-18, cumulative = -8.9963e-09 DICPCG: Solving for p_rgh, Initial residual = 0.415889, Final residual = 0.0161375, No Iterations 6 time step continuity errors : sum local = 1.76416e-11, global = -6.73244e-18, cumulative = -8.9963e-09 DICPCG: Solving for p_rgh, Initial residual = 0.520079, Final residual = 9.88526e-08, No Iterations 615 time step continuity errors : sum local = 2.77983e-16, global = 2.135e-17, cumulative = -8.9963e-09 smoothSolver: Solving for epsilon, Initial residual = 0.923994, Final residual = 4.41476e-07, No Iterations 4 smoothSolver: Solving for k, Initial residual = 0.710582, Final residual = 5.52735e-07, No Iterations 22 ExecutionTime = 11983.3 s ClockTime = 12006 s forces bridgeLoads write: sum of forces: pressure : (-2.09321e+06 -3.7909e+06 -29224.4) viscous : (-140514 -90938.7 665.627) porous : (0 0 0) sum of moments: pressure : (-1.89788e+06 1.0424e+06 -306745) viscous : (-45362.4 70040.7 -6928.93) porous : (0 0 0) Courant Number mean: 1.02963e-07 max: 71.2109 Interface Courant Number mean: 4.63004e-09 max: 1.82454e-05 deltaT = 2.42193e-11 Time = 3.18936 PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 4.02161e-12, Final residual = 4.02161e-12, No Iterations 0 Phase-1 volume fraction = 0.305209 Min(alpha.water) = -6.89936e-06 Max(alpha.water) = 1.00001 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.305209 Min(alpha.water) = -6.89936e-06 Max(alpha.water) = 1.00001 DICPCG: Solving for p_rgh, Initial residual = 0.99995, Final residual = 0.0420752, No Iterations 7 time step continuity errors : sum local = 1.2412e-11, global = 6.14918e-21, cumulative = -8.9963e-09 DICPCG: Solving for p_rgh, Initial residual = 1.17775e-05, Final residual = 5.12902e-07, No Iterations 4 time step continuity errors : sum local = 0.000107511, global = 2.1704e-18, cumulative = -8.9963e-09 DICPCG: Solving for p_rgh, Initial residual = 1.10584e-05, Final residual = 9.06598e-08, No Iterations 18 time step continuity errors : sum local = 1.37982e-05, global = 7.69188e-21, cumulative = -8.9963e-09 smoothSolver: Solving for epsilon, Initial residual = 0.00121858, Final residual = 8.73656e-14, No Iterations 1 smoothSolver: Solving for k, Initial residual = 0.998659, Final residual = 1.13328e-06, No Iterations 1000 ExecutionTime = 12017.4 s ClockTime = 12040 s forces bridgeLoads write: sum of forces: pressure : (1.68331e+22 8.35418e+19 1.82866e+15) viscous : (-6.66721e+06 -7.07181e+07 1.06096e+08) porous : (0 0 0) sum of moments: pressure : (4.14398e+19 -8.41021e+21 7.99315e+19) viscous : (-3.57783e+07 -5.15335e+06 -5.19118e+06) porous : (0 0 0) Courant Number mean: 0.000182021 max: 463651 Interface Courant Number mean: 3.32243e-09 max: 0.000700647 deltaT = 2.6118e-17 Time = 3.18936 PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 1.73072e-16, Final residual = 1.73072e-16, No Iterations 0 Phase-1 volume fraction = 0.305209 Min(alpha.water) = -6.89936e-06 Max(alpha.water) = 1.00001 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.305209 Min(alpha.water) = -6.89936e-06 Max(alpha.water) = 1.00001 DICPCG: Solving for p_rgh, Initial residual = 0.946778, Final residual = 0.0397426, No Iterations 3 time step continuity errors : sum local = 1.86432e-10, global = -2.60064e-25, cumulative = -8.9963e-09 DICPCG: Solving for p_rgh, Initial residual = 0.000111561, Final residual = 5.08022e-06, No Iterations 3 time step continuity errors : sum local = 1.04018e-07, global = 2.12946e-22, cumulative = -8.9963e-09 DICPCG: Solving for p_rgh, Initial residual = 5.61518e-05, Final residual = 9.39591e-08, No Iterations 64 time step continuity errors : sum local = 1.60715e-09, global = -7.45538e-14, cumulative = -8.99638e-09 smoothSolver: Solving for epsilon, Initial residual = 0.998169, Final residual = 3.42172e-07, No Iterations 11 smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.939896, No Iterations 1000 ExecutionTime = 12055.2 s ClockTime = 12078 s forces bridgeLoads write: sum of forces: pressure : (1.49918e+27 -1.21963e+30 -7.08408e+28) viscous : (-7.45654e+18 -5.13495e+19 1.62365e+19) porous : (0 0 0) sum of moments: pressure : (-6.09217e+29 4.99264e+27 -9.87025e+28) viscous : (-2.57415e+19 2.4802e+18 -3.71744e+18) porous : (0 0 0) Courant Number mean: 4.32858e-06 max: 12460.3 Interface Courant Number mean: 2.47922e-11 max: 4.67598e-06 deltaT = 1.04805e-21 Time = 3.18936 PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 2.68028e-16, Final residual = 2.68028e-16, No Iterations 0 Phase-1 volume fraction = 0.305209 Min(alpha.water) = -6.89936e-06 Max(alpha.water) = 1.00001 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.305209 Min(alpha.water) = -6.89936e-06 Max(alpha.water) = 1.00001 DICPCG: Solving for p_rgh, Initial residual = 0.858654, Final residual = 0.0348155, No Iterations 5 time step continuity errors : sum local = 6.31934e-12, global = -3.59149e-18, cumulative = -8.99638e-09 DICPCG: Solving for p_rgh, Initial residual = 1.4193e-11, Final residual = 1.4193e-11, No Iterations 0 time step continuity errors : sum local = 7.46006e+12, global = -0.000702905, cumulative = -0.000702914 DICPCG: Solving for p_rgh, Initial residual = 1.4193e-11, Final residual = 1.4193e-11, No Iterations 0 time step continuity errors : sum local = 7.46009e+12, global = 0.00030909, cumulative = -0.000393824 smoothSolver: Solving for epsilon, Initial residual = 0.888711, Final residual = 2.12041e-08, No Iterations 3 bounding epsilon, min: -1.55108e+38 max: 2.0259e+66 average: 1.00553e+59 smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.827953, No Iterations 1000 bounding k, min: -2.25591e+06 max: 2.18236e+57 average: 1.70514e+50 ExecutionTime = 12089.6 s ClockTime = 12112 s forces bridgeLoads write: sum of forces: pressure : (-2.48095e+64 -1.97533e+62 -1.13171e+61) viscous : (6.46368e+43 1.02607e+43 5.6584e+41) porous : (0 0 0) sum of moments: pressure : (-9.86727e+61 1.23903e+64 -4.16812e+62) viscous : (5.12643e+42 -3.23491e+43 9.09553e+41) porous : (0 0 0) Courant Number mean: 3.73024e+12 max: 6.1465e+22 Interface Courant Number mean: 6.32826e-16 max: 8.40748e-11 deltaT = 8.52559e-45 Time = 3.18936 PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 2.16154e-17, Final residual = 2.16154e-17, No Iterations 0 Phase-1 volume fraction = 0.305209 Min(alpha.water) = -6.89936e-06 Max(alpha.water) = 1.00001 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.305209 Min(alpha.water) = -6.89936e-06 Max(alpha.water) = 1.00001 DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 0.0487684, No Iterations 12 time step continuity errors : sum local = 1.55181e-12, global = -3.48026e-26, cumulative = -0.000393824 DICPCG: Solving for p_rgh, Initial residual = 1.97172e-23, Final residual = 1.97172e-23, No Iterations 0 time step continuity errors : sum local = 1.2961e+10, global = 3.59592e-07, cumulative = -0.000393465 DICPCG: Solving for p_rgh, Initial residual = 3.19389e-23, Final residual = 3.19389e-23, No Iterations 0 time step continuity errors : sum local = 2.09949e+10, global = 9.81777e-07, cumulative = -0.000392483 -------------------------------------------------------------------------- mpirun noticed that process rank 56 with PID 0 on node n0885 exited on signal 8 (Floating point exception). -------------------------------------------------------------------------- Code:
[56] #0 Foam::error::printStack(Foam::Ostream&) in "/gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" [56] #1 Foam::sigFpe::sigHandler(int) in "/gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" [56] #2 ? in "/lib64/libc.so.6" [56] #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) in "/gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" [56] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const in "/gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" [56] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" [56] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so" [56] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/bin/interFoam" [56] #8 Foam::fvMatrix<double>::solve() in "/gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/bin/interFoam" [56] #9 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) in "/gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so" [56] #10 Foam::RASModels::kEpsilon<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() in "/gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so" [56] #11 ? in "/gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/bin/interFoam" [56] #12 __libc_start_main in "/lib64/libc.so.6" [56] #13 ? in "/gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/bin/interFoam" [n0885:12114] *** Process received signal *** [n0885:12114] Signal: Floating point exception (8) [n0885:12114] Signal code: (-6) [n0885:12114] Failing at address: 0x573de00002f52 [n0885:12114] [ 0] /lib64/libc.so.6(+0x32510)[0x2b4b72740510] [n0885:12114] [ 1] /lib64/libc.so.6(gsignal+0x35)[0x2b4b72740495] [n0885:12114] [ 2] /lib64/libc.so.6(+0x32510)[0x2b4b72740510] [n0885:12114] [ 3] /gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam22symGaussSeidelSmoother6smoothERKNS_4wordERNS_5FieldIdEERKNS_9lduMatrixERKS5_RKNS_10FieldFieldIS4_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEEhi+0x485)[0x2b4b715e7bf5] [n0885:12114] [ 4] /gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam22symGaussSeidelSmoother6smoothERNS_5FieldIdEERKS2_hi+0x28)[0x2b4b715e80d8] [n0885:12114] [ 5] /gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam12smoothSolver5solveERNS_5FieldIdEERKS2_h+0x4a3)[0x2b4b715ddd43] [n0885:12114] [ 6] /gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x14c)[0x2b4b6ed4e90c] [n0885:12114] [ 7] interFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x110)[0x4ceff0] [n0885:12114] [ 8] interFoam(_ZN4Foam8fvMatrixIdE5solveEv+0x8f)[0x4cf24f] [n0885:12114] [ 9] /gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam5solveIdEENS_17SolverPerformanceIT_EERKNS_3tmpINS_8fvMatrixIS2_EEEE+0x1f)[0x2b4b6df4f39f] [n0885:12114] [10] /gscratch/motley/shared/OpenFOAM500/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASModels8kEpsilonINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE7correctEv+0x8d4)[0x2b4b6dfa9664] [n0885:12114] [11] interFoam[0x4c1019] [n0885:12114] [12] /lib64/libc.so.6(__libc_start_main+0xfd)[0x2b4b7272cd1d] [n0885:12114] [13] interFoam[0x43f9c9] [n0885:12114] *** End of error message *** [n0881][[52882,1],24][/sw/contrib/openmpi/1.10.1/build_gcc_5.1.0/../openmpi-1.10.1/ompi/mca/btl/tcp/btl_tcp_frag.c:237:mca_btl_tcp_frag_recv] mca_btl_tcp_frag_recv: readv failed: Connection reset by peer (104) [n0884][[52882,1],48][/sw/contrib/openmpi/1.10.1/build_gcc_5.1.0/../openmpi-1.10.1/ompi/mca/btl/tcp/btl_tcp_frag.c:237:mca_btl_tcp_frag_recv] mca_btl_tcp_frag_recv: readv failed: Connection reset by peer (104) [n0882][[52882,1],35][/sw/contrib/openmpi/1.10.1/build_gcc_5.1.0/../openmpi-1.10.1/ompi/mca/btl/tcp/btl_tcp_frag.c:237:mca_btl_tcp_frag_recv] mca_btl_tcp_frag_recv: readv failed: Connection reset by peer (104) [n0884][[52882,1],51][/sw/contrib/openmpi/1.10.1/build_gcc_5.1.0/../openmpi-1.10.1/ompi/mca/btl/tcp/btl_tcp_frag.c:237:mca_btl_tcp_frag_recv] mca_btl_tcp_frag_recv: readv failed: Connection reset by peer (104) fvSchemes Code:
ddtSchemes { default Euler; } gradSchemes { default Gauss linear; } divSchemes { div(rhoPhi,U) Gauss linearUpwind grad(U); div(phi,alpha) Gauss vanLeer; div(phirb,alpha) Gauss linear; div(phi,k) Gauss upwind; div(phi,epsilon) Gauss upwind; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } Code:
solvers { "alpha.water.*" { nAlphaCorr 2; nAlphaSubCycles 1; cAlpha 1; MULESCorr yes; nLimiterIter 3; solver smoothSolver; smoother symGaussSeidel; tolerance 1e-8; relTol 0; } "pcorr.*" { solver PCG; preconditioner DIC; tolerance 1e-5; relTol 0; } p_rgh { solver PCG; preconditioner DIC; tolerance 1e-07; relTol 0.05; } p_rghFinal { $p_rgh; relTol 0; } "(U|k|epsilon).*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-06; relTol 0; minIter 1; } } PIMPLE { momentumPredictor no; nOuterCorrectors 1; nCorrectors 3; nNonOrthogonalCorrectors 0; } relaxationFactors { equations { ".*" 1; } } |
|
May 24, 2018, 04:35 |
|
#2 |
New Member
chunping Ren
Join Date: May 2018
Posts: 1
Rep Power: 0 |
The relaxation factors (in fvSolution) are too high. Try to use less than 0.5
|
|
May 24, 2018, 15:32 |
|
#3 |
Member
Andrew O. Winter
Join Date: Aug 2015
Location: Seattle, WA, USA
Posts: 78
Rep Power: 11 |
Thanks for the tip! I completely overlooked the fact that the relaxation factor was set to 1. I'll try it out when I have some time early next week and if all goes well I'll post again to confirm that the issue is resolved.
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
DPM udf error | haghshenasfard | FLUENT | 0 | April 13, 2016 07:35 |
[swak4Foam] GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh | gschaider | OpenFOAM Community Contributions | 300 | October 29, 2014 19:00 |
OpenFOAM without MPI | kokizzu | OpenFOAM Installation | 4 | May 26, 2014 10:17 |
Problem with compile the setParabolicInlet | ivanyao | OpenFOAM Running, Solving & CFD | 6 | September 5, 2008 21:50 |
DecomposePar links against liblamso0 with OpenMPI | jens_klostermann | OpenFOAM Bugs | 11 | June 28, 2007 18:51 |