View Issue Details

IDProjectCategoryView StatusLast Update
0003675OpenFOAMBugpublic2021-05-18 14:36
Reporteroder Assigned Tohenry  
PrioritynormalSeverityminorReproducibilitysometimes
Status resolvedResolutionfixed 
PlatformLinux-5.5.16-100.fc30.x86_64OSFedoraOS Version30
Summary0003675: Undefined behavior in function Foam::ISstream::getLine
DescriptionThe undefined behavior occurs at the following line
https://github.com/OpenFOAM/OpenFOAM-8/blob/6a90c5e21f9f73f21de679dd775e62a783fbd5c3/src/OpenFOAM/db/IOstreams/Sstreams/ISstream.C#L778

Just a few lines above (L774), a line is read from a file into Foam::string s. Foam::string is a specialization of class std::string. However, if the line was empty, string s is empty. Calling s.back() on an empty string leads to undefined behavior, as stated here:
https://en.cppreference.com/w/cpp/string/basic_string/back

I believe this is also present in the OpenFoam-dev version, although there was a global continuation flag added there, which, in case it is true, leads to the same behavior.


Steps To ReproduceIn my instance there is no problem with a version compiled with optimization flags, however, when compiled with debugging flags (WM_COMPILE_OPTION=Debug in etc/bashrc) the program gets aborted. Thus, I marked the reproducibility to sometimes and severity as minor, although the undefined behavior is always there.

This part of code is used when reading the template file etc/codeTemplates/dynamicCode/codedFvOptionTemplate.C when creating the C++ source files with scalarCodedSource source term in system/fvOptions. The template file contains empty lines.
Additional InformationThe solution I would propose is to insert the following line just above the while loop in line L778:
if (s.empty()) return *this;

If the line is empty there is nothing to concatenate and no line to continue anyway.

If required I can provide a simplified test case with which I was debugging the code.
TagsNo tags attached.

Activities

henry

2021-05-17 18:29

manager   ~0012027

Can you provide a means to reproduce the problem?

oder

2021-05-18 10:21

reporter   ~0012032

I hope I did this correctly. Attached is the tarball with case files (issue-0003675-case.tar.gz) and a tarball with logs (issue-0003675-logs.tar.gz):

issue-0003675-logs/log.buoyantSimpleFoam.Debug ... Run with OF-8 from master branch, compiled with Debug flag. Gets aborted.
issue-0003675-logs/log.buoyantSimpleFoam.fix ... Run with OF-8 from master branch with proposed fix.
issue-0003675-logs/log.buoyantSimpleFoam.Opt ... Run with OF-8 from master branch, compiled with Opt flag. No problems.
issue-0003675-logs/log.buoyantSimpleFoam.system ... Run with our system version of OF-8. No problems.

I also attach the relevant output of the debugging session with gdb of OF-8 from master branch, compiled with debugging symbols (Debug flag).
issue-0003675-gdb.txt (5,606 bytes)   
$ gdb buoyantSimpleFoam

(gdb) break ISstream.C:778
No source file named ISstream.C.
Make breakpoint pending on future shared library load? (y or [n]) y
Breakpoint 1 (ISstream.C:778) pending.

(gdb) ignore 1 8
Will ignore next 8 crossings of breakpoint 1.

(gdb) run
(gdb) p s._M_dataplus._M_p
$1 = (std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> >::pointer) 0x9c7220 "    This file is part of OpenFOAM."

(gdb) p s.c_str()
$2 = 0x9c7220 "    This file is part of OpenFOAM."

(gdb) c
(gdb) p s._M_dataplus._M_p
$3 = (std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> >::pointer) 0x9c7220 ""

(gdb) si 14
0x00007ffff3231260 in std::__replacement_assert(char const*, int, char const*, char const*) () from /lib64/libstdc++.so.6

(gdb) disassemble
Dump of assembler code for function _ZSt20__replacement_assertPKciS0_S0_:
=> 0x00007ffff3231260 <+0>:     endbr64
   0x00007ffff3231264 <+4>:     push   %rax
   0x00007ffff3231265 <+5>:     pop    %rax
   0x00007ffff3231266 <+6>:     xor    %eax,%eax
   0x00007ffff3231268 <+8>:     mov    %rcx,%r8
   0x00007ffff323126b <+11>:    sub    $0x8,%rsp
   0x00007ffff323126f <+15>:    mov    %rdx,%rcx
   0x00007ffff3231272 <+18>:    mov    %esi,%edx
   0x00007ffff3231274 <+20>:    mov    %rdi,%rsi
   0x00007ffff3231277 <+23>:    lea    0xd3f52(%rip),%rdi        # 0x7ffff33051d0
   0x00007ffff323127e <+30>:    callq  0x7ffff3206250 <printf@plt>
   0x00007ffff3231283 <+35>:    callq  0x7ffff32066f0 <abort@plt>
End of assembler dump.

(gdb) ni 11
(gdb) disassemble
Dump of assembler code for function _ZSt20__replacement_assertPKciS0_S0_:
   0x00007ffff3231260 <+0>:     endbr64
   0x00007ffff3231264 <+4>:     push   %rax
   0x00007ffff3231265 <+5>:     pop    %rax
   0x00007ffff3231266 <+6>:     xor    %eax,%eax
   0x00007ffff3231268 <+8>:     mov    %rcx,%r8
   0x00007ffff323126b <+11>:    sub    $0x8,%rsp
   0x00007ffff323126f <+15>:    mov    %rdx,%rcx
   0x00007ffff3231272 <+18>:    mov    %esi,%edx
   0x00007ffff3231274 <+20>:    mov    %rdi,%rsi
   0x00007ffff3231277 <+23>:    lea    0xd3f52(%rip),%rdi        # 0x7ffff33051d0
   0x00007ffff323127e <+30>:    callq  0x7ffff3206250 <printf@plt>
=> 0x00007ffff3231283 <+35>:    callq  0x7ffff32066f0 <abort@plt>
End of assembler dump.

(gdb) bt
#0  0x00007ffff3231283 in std::__replacement_assert(char const*, int, char const*, char const*) () from /lib64/libstdc++.so.6
#1  0x00007ffff32b52a6 in std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> >::back() () from /lib64/libstdc++.so.6
#2  0x00007ffff39df778 in Foam::ISstream::getLine (this=0x7fffffff1680, s=...) at db/IOstreams/Sstreams/ISstream.C:778
#3  0x00007ffff3a1b886 in Foam::dynamicCode::copyAndFilter (is=..., os=..., mapping=...) at db/dynamicLibrary/dynamicCode/dynamicCode.C:129
#4  0x00007ffff3a1d439 in Foam::dynamicCode::copyOrCreateFiles (this=0x7fffffff1a40, verbose=true) at db/dynamicLibrary/dynamicCode/dynamicCode.C:454
#5  0x00007ffff3a1f8a6 in Foam::codedBase::createLibrary (this=0x9bc408, dynCode=..., context=...) at db/dynamicLibrary/codedBase/codedBase.C:190
#6  0x00007ffff3a20283 in Foam::codedBase::updateLibrary (this=0x9bc408) at db/dynamicLibrary/codedBase/codedBase.C:331
#7  0x00007ffff5e29fb4 in Foam::fv::CodedSource<double>::read (this=0x9bc220, dict=...) at sources/general/codedSource/CodedSourceIO.C:50
#8  0x00007ffff5e2965f in Foam::fv::CodedSource<double>::CodedSource (this=0x9bc220, name=..., modelType=..., dict=..., mesh=...) at sources/general/codedSource/CodedSource.C:132
#9  0x00007ffff5e28fa6 in Foam::fv::option::adddictionaryConstructorToTable<Foam::fv::CodedSource<double> >::New (name=..., modelType=..., dict=..., mesh=...) at /data2/nobackup/oder/src/OpenFOAM-8/src/finiteVolume/lnInclude/fvOption.H:106
#10 0x00007ffff7219f2f in Foam::fv::option::New (name=..., coeffs=..., mesh=...) at cfdTools/general/fvOptions/fvOption.C:97
#11 0x00007ffff721d33b in Foam::fv::optionList::reset (this=0x95e7a0, dict=...) at cfdTools/general/fvOptions/fvOptionList.C:132
#12 0x00007ffff721d11e in Foam::fv::optionList::optionList (this=0x95e7a0, mesh=..., dict=...) at cfdTools/general/fvOptions/fvOptionList.C:94
#13 0x00007ffff721e5f4 in Foam::fv::options::options (this=0x95e630, mesh=...) at cfdTools/general/fvOptions/fvOptions.C:97
#14 0x00007ffff721e795 in Foam::fv::options::New (mesh=...) at cfdTools/general/fvOptions/fvOptions.C:119
#15 0x0000000000478492 in main (argc=1, argv=0x7fffffffb3c8) at /data2/nobackup/oder/src/OpenFOAM-8/src/finiteVolume/lnInclude/createFvOptions.H:1

(gdb) up 4
(gdb) p is.pathname_._M_dataplus._M_p
$4 = (std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> >::pointer) 0x9c9f20 "/data2/nobackup/oder/src/OpenFOAM-8/etc/codeTemplates/dynamicCode/codedFvOptionTemplate.C"

(gdb) !head -n 10 /data2/nobackup/oder/src/OpenFOAM-8/etc/codeTemplates/dynamicCode/codedFvOptionTemplate.C | nl -b a
     1  /*---------------------------------------------------------------------------*\
     2    =========                 |
     3    \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
     4     \\    /   O peration     | Website:  https://openfoam.org
     5      \\  /    A nd           | Copyright (C) YEAR OpenFOAM Foundation
     6       \\/     M anipulation  |
     7  -------------------------------------------------------------------------------
     8  License
     9      This file is part of OpenFOAM.
    10

(gdb) ni

Program received signal SIGABRT, Aborted.
0x00007ffff2e7be35 in raise () from /lib64/libc.so.6
(gdb)

issue-0003675-gdb.txt (5,606 bytes)   

henry

2021-05-18 13:16

manager   ~0012033

I have been unable to reproduce the problem and so cannot test the fix in OpenFOAM-8:

commit 6d8538864cc3d80b285e2dbfa7c01097de9f0276 (HEAD -> master, origin/master, origin/HEAD)
Author: Henry Weller <http://cfd.direct>
Date: Tue May 18 13:10:50 2021 +0100

    ISstream: Added test for empty string before parsing continuation lines
    
    Resolves bug-report https://bugs.openfoam.org/view.php?id=3675

or OpenFOAM-dev:

commit 58afe5dccbc864afb41428a11b622d49d4c2c1f0 (HEAD -> master, origin/master, origin/HEAD)
Author: Henry Weller <http://cfd.direct>
Date: Tue May 18 13:12:03 2021 +0100

    ISstream: Added test for empty string before parsing continuation lines
    
    Resolves bug-report https://bugs.openfoam.org/view.php?id=3675

Could you test and report back?

oder

2021-05-18 14:24

reporter   ~0012034

I have tested it with OpenFOAM-8 and the fix works (log file attached).

As far as I have checked, in OpenFOAM-dev, getLine in this particular case is called with continuation == false, so s.back() is never called anyway. [1]

[1] https://github.com/OpenFOAM/OpenFOAM-dev/blob/40bc30c0f7f8394437737f5311b1d2960f3f8a06/src/OpenFOAM/db/dynamicLibrary/dynamicCode/dynamicCode.C#L130

Thank you!
log.buoyantSimpleFoam (13,035 bytes)   
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  8
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 8-6d8538864cc3
Exec   : buoyantSimpleFoam
Date   : May 18 2021
Time   : 15:04:35
Host   : "ealin27.private.vki.eu"
PID    : 50414
I/O    : uncollated
Case   : /data2/nobackup/oder/projects/Patricia/Kawamura/OF-bug
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
        U: tolerance 1e-05
        h: tolerance 1e-07

Reading thermophysical properties

Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleEnthalpy;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model LaunderSharmaKE
RAS
{
    RASModel        LaunderSharmaKE;
    turbulence      on;
    printCoeffs     on;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
}

Creating thermophysical transport model

Selecting thermophysical transport type RAS
Selecting RAS thermophysical transport model eddyDiffusivity

Reading g

Reading hRef
Calculating field g.h


Reading pRef
Reading field p_rgh

No MRF models present

Radiation model not active: radiationProperties not found
Selecting radiationModel none
Creating finite volume options from "system/fvOptions"

Selecting finite volume options model type meanVelocityForce
    Source: momentumSource
    - selecting all cells
    - selected 5000 cell(s) with volume 3.32168e-06
    Initial pressure gradient = 0

Selecting finite volume options model type scalarCodedSource
    Source: energySource
    - selecting all cells
    - selected 5000 cell(s) with volume 3.32168e-06
Using dynamicCode for fvOption:: sourceTime at line 42 in "/data2/nobackup/oder/projects/Patricia/Kawamura/OF-bug/system/fvOptions/energySource/scalarCodedSourceCoeffs"
Creating new library in "dynamicCode/sourceTime/platforms/linux64GccDPInt32Debug/lib/libsourceTime_934561a688b79242ec892890ff7f93af0123113a.so"
Invoking "wmake -s libso /data2/nobackup/oder/projects/Patricia/Kawamura/OF-bug/dynamicCode/sourceTime"
wmake libso /data2/nobackup/oder/projects/Patricia/Kawamura/OF-bug/dynamicCode/sourceTime
    ln: ./lnInclude
    wmkdep: codedFvOptionTemplate.C
    Ctoo: codedFvOptionTemplate.C
    ld: /data2/nobackup/oder/projects/Patricia/Kawamura/OF-bug/dynamicCode/sourceTime/../platforms/linux64GccDPInt32Debug/lib/libsourceTime_934561a688b79242ec892890ff7f93af0123113a.so

Starting time loop

Time = 0.01

DILUPBiCG:  Solving for Ux, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
Pressure gradient source: uncorrected Ubar = 0, pressure gradient = 6.60848e+07
Selecting finite volume options model type sourceTime
    Source: sourceTime
    - selecting all cells
    - selected 5000 cell(s) with volume 3.32168e-06
DILUPBiCG:  Solving for h, Initial residual = 0.999988, Final residual = 8.18579e-10, No Iterations 3
DICPCG:  Solving for p_rgh, Initial residual = 0.579293, Final residual = 6.44706e-09, No Iterations 69
time step continuity errors : sum local = 6.0289e-16, global = 4.70813e-21, cumulative = 4.70813e-21
Pressure gradient source: uncorrected Ubar = 0.00801752, pressure gradient = 5.28948e+07
DILUPBiCG:  Solving for epsilon, Initial residual = 0.143133, Final residual = 8.41709e-09, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 9.03117e-11, No Iterations 4
ExecutionTime = 0.28 s  ClockTime = 1 s

Time = 0.02

DILUPBiCG:  Solving for Ux, Initial residual = 0.672352, Final residual = 9.55384e-11, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.999932, Final residual = 1.29541e-09, No Iterations 3
Pressure gradient source: uncorrected Ubar = 483.414, pressure gradient = -5.62741e+06
DILUPBiCG:  Solving for h, Initial residual = 0.00217435, Final residual = 2.70714e-09, No Iterations 2
DICPCG:  Solving for p_rgh, Initial residual = 0.999801, Final residual = 6.17561e-09, No Iterations 120
time step continuity errors : sum local = 7.08284e-11, global = -9.84794e-16, cumulative = -9.8479e-16
Pressure gradient source: uncorrected Ubar = 439.854, pressure gradient = -353609
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0039612, Final residual = 8.81878e-09, No Iterations 3
DILUPBiCG:  Solving for k, Initial residual = 0.98768, Final residual = 3.32357e-10, No Iterations 5
ExecutionTime = 0.42 s  ClockTime = 2 s

Time = 0.03

DILUPBiCG:  Solving for Ux, Initial residual = 0.00429494, Final residual = 7.86621e-11, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.33209, Final residual = 1.52472e-09, No Iterations 3
Pressure gradient source: uncorrected Ubar = 0.00419748, pressure gradient = 80840.9
DILUPBiCG:  Solving for h, Initial residual = 3.15452e-05, Final residual = 2.31642e-09, No Iterations 2
DICPCG:  Solving for p_rgh, Initial residual = 0.875392, Final residual = 7.98824e-09, No Iterations 131
time step continuity errors : sum local = 8.53143e-13, global = 2.60271e-18, cumulative = -9.82187e-16
Pressure gradient source: uncorrected Ubar = 0.00874653, pressure gradient = 25900.2
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0206072, Final residual = 3.69845e-09, No Iterations 3
DILUPBiCG:  Solving for k, Initial residual = 0.0087236, Final residual = 7.39964e-11, No Iterations 3
ExecutionTime = 0.57 s  ClockTime = 2 s

Time = 0.04

DILUPBiCG:  Solving for Ux, Initial residual = 0.0323507, Final residual = 1.97488e-09, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.35906, Final residual = 5.54657e-11, No Iterations 3
Pressure gradient source: uncorrected Ubar = 0.0479449, pressure gradient = -5506.3
DILUPBiCG:  Solving for h, Initial residual = 0.00103375, Final residual = 2.45738e-10, No Iterations 2
DICPCG:  Solving for p_rgh, Initial residual = 0.601559, Final residual = 4.44023e-09, No Iterations 130
time step continuity errors : sum local = 1.91763e-12, global = -9.55238e-19, cumulative = -9.83142e-16
Pressure gradient source: uncorrected Ubar = 0.0473258, pressure gradient = -3005.45
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00438388, Final residual = 1.15083e-09, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.00363259, Final residual = 1.12763e-09, No Iterations 2
ExecutionTime = 0.73 s  ClockTime = 2 s

Time = 0.05

DILUPBiCG:  Solving for Ux, Initial residual = 0.0290605, Final residual = 4.37551e-10, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.130584, Final residual = 2.2614e-11, No Iterations 3
Pressure gradient source: uncorrected Ubar = 0.0586101, pressure gradient = -49381.7
DILUPBiCG:  Solving for h, Initial residual = 0.000664288, Final residual = 8.08869e-10, No Iterations 2
DICPCG:  Solving for p_rgh, Initial residual = 0.581652, Final residual = 3.66123e-09, No Iterations 129
time step continuity errors : sum local = 3.17337e-12, global = 2.36902e-18, cumulative = -9.80773e-16
Pressure gradient source: uncorrected Ubar = 0.0572639, pressure gradient = -45996.1
DILUPBiCG:  Solving for epsilon, Initial residual = 0.000789131, Final residual = 3.45826e-10, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.00144369, Final residual = 5.18166e-10, No Iterations 2
ExecutionTime = 0.87 s  ClockTime = 2 s

Time = 0.06

DILUPBiCG:  Solving for Ux, Initial residual = 0.0140218, Final residual = 4.48639e-10, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.0586515, Final residual = 5.57355e-09, No Iterations 2
Pressure gradient source: uncorrected Ubar = 0.0503411, pressure gradient = -63829.1
DILUPBiCG:  Solving for h, Initial residual = 0.000275275, Final residual = 2.64868e-10, No Iterations 2
DICPCG:  Solving for p_rgh, Initial residual = 0.405713, Final residual = 9.61165e-09, No Iterations 128
time step continuity errors : sum local = 3.91875e-12, global = -9.21501e-18, cumulative = -9.89988e-16
Pressure gradient source: uncorrected Ubar = 0.0496138, pressure gradient = -62553.8
DILUPBiCG:  Solving for epsilon, Initial residual = 0.000278952, Final residual = 1.80413e-10, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.000409126, Final residual = 1.22179e-10, No Iterations 2
ExecutionTime = 1.01 s  ClockTime = 2 s

Time = 0.07

DILUPBiCG:  Solving for Ux, Initial residual = 0.0080979, Final residual = 5.34282e-10, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.0256374, Final residual = 2.60029e-09, No Iterations 2
Pressure gradient source: uncorrected Ubar = 0.0333121, pressure gradient = -52904.3
DILUPBiCG:  Solving for h, Initial residual = 0.000188539, Final residual = 3.93459e-10, No Iterations 2
DICPCG:  Solving for p_rgh, Initial residual = 0.289926, Final residual = 4.70118e-09, No Iterations 125
time step continuity errors : sum local = 4.88752e-12, global = 6.92995e-18, cumulative = -9.83058e-16
Pressure gradient source: uncorrected Ubar = 0.0338645, pressure gradient = -53681.5
DILUPBiCG:  Solving for epsilon, Initial residual = 0.000152831, Final residual = 1.48652e-10, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.000173423, Final residual = 4.25529e-11, No Iterations 2
ExecutionTime = 1.16 s  ClockTime = 2 s

Time = 0.08

DILUPBiCG:  Solving for Ux, Initial residual = 0.00587445, Final residual = 5.7389e-10, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.0240492, Final residual = 2.71413e-09, No Iterations 2
Pressure gradient source: uncorrected Ubar = 0.0247687, pressure gradient = -37116.6
DILUPBiCG:  Solving for h, Initial residual = 0.000154987, Final residual = 5.59656e-10, No Iterations 2
DICPCG:  Solving for p_rgh, Initial residual = 0.359624, Final residual = 6.98975e-09, No Iterations 124
time step continuity errors : sum local = 5.65544e-12, global = -8.60527e-19, cumulative = -9.83919e-16
Pressure gradient source: uncorrected Ubar = 0.0261542, pressure gradient = -38606.8
DILUPBiCG:  Solving for epsilon, Initial residual = 8.66088e-05, Final residual = 1.31376e-10, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 7.56366e-05, Final residual = 1.63743e-11, No Iterations 2
ExecutionTime = 1.3 s  ClockTime = 2 s

Time = 0.09

DILUPBiCG:  Solving for Ux, Initial residual = 0.00639378, Final residual = 7.1454e-10, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.0248617, Final residual = 5.48245e-09, No Iterations 2
Pressure gradient source: uncorrected Ubar = 0.0242948, pressure gradient = -25390.3
DILUPBiCG:  Solving for h, Initial residual = 0.000151275, Final residual = 6.90081e-10, No Iterations 2
DICPCG:  Solving for p_rgh, Initial residual = 0.196092, Final residual = 3.0253e-09, No Iterations 122
time step continuity errors : sum local = 6.6295e-12, global = 1.87006e-19, cumulative = -9.83732e-16
Pressure gradient source: uncorrected Ubar = 0.0258247, pressure gradient = -26664
DILUPBiCG:  Solving for epsilon, Initial residual = 5.03314e-05, Final residual = 1.55761e-10, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 4.21344e-05, Final residual = 1.99701e-11, No Iterations 2
ExecutionTime = 1.44 s  ClockTime = 3 s

Time = 0.1

DILUPBiCG:  Solving for Ux, Initial residual = 0.00450563, Final residual = 5.9589e-10, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.0220232, Final residual = 5.62992e-12, No Iterations 3
Pressure gradient source: uncorrected Ubar = 0.0271459, pressure gradient = -17901.9
DILUPBiCG:  Solving for h, Initial residual = 0.000163334, Final residual = 7.56622e-10, No Iterations 2
DICPCG:  Solving for p_rgh, Initial residual = 0.160125, Final residual = 6.32183e-09, No Iterations 123
time step continuity errors : sum local = 8.39001e-12, global = 1.21477e-17, cumulative = -9.71584e-16
Pressure gradient source: uncorrected Ubar = 0.0284683, pressure gradient = -18791.6
DILUPBiCG:  Solving for epsilon, Initial residual = 2.81283e-05, Final residual = 8.87772e-11, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 3.35105e-05, Final residual = 4.72719e-11, No Iterations 2
ExecutionTime = 1.59 s  ClockTime = 3 s

End

log.buoyantSimpleFoam (13,035 bytes)   

Issue History

Date Modified Username Field Change
2021-05-17 18:24 oder New Issue
2021-05-17 18:29 henry Note Added: 0012027
2021-05-18 10:21 oder File Added: issue-0003675-logs.tar.gz
2021-05-18 10:21 oder File Added: issue-0003675-gdb.txt
2021-05-18 10:21 oder File Added: issue-0003675-case.tar.gz
2021-05-18 10:21 oder Note Added: 0012032
2021-05-18 13:16 henry Note Added: 0012033
2021-05-18 14:24 oder File Added: log.buoyantSimpleFoam
2021-05-18 14:24 oder Note Added: 0012034
2021-05-18 14:36 henry Assigned To => henry
2021-05-18 14:36 henry Status new => resolved
2021-05-18 14:36 henry Resolution open => fixed
2021-05-18 14:36 henry Fixed in Version => 8