[amolganjare@gmail.com]

I am trying to use the utility "reconstructPointDistribution" of OpenQBMM version 9 with OpenFOAM-v2206.
I have used a tutorial case of Tyalor Dispersion solved with pbeTransportFoam. When I use generateMoment utility it works fine.

I want to use reconstructPointDistribution utility on result or say on the initial "0" time, I have kept the pointDistributionDict in system directory.



When I run the command reconstructPointDistribution I get following error.



[amol.ganjare@mercury TaylorDispersion]$ reconstructPointDistribution
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2206 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : _76d719d1e6-20220624 OPENFOAM=2206 version=v2206
Arch : "LSB;label=32;scalar=64"
Exec : reconstructPointDistribution
Date : Dec 23 2022
Time : 16:23:14
Host : mercury
PID : 32539
I/O : uncollated
Case : /nas1/amol.ganjare/OF_2206/OpenFOAM-v2206/modules/OpenQBMM/tutorials/pbeTransportFoam/TaylorDispersion
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

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

--> FOAM Warning :
From static Foam::IOstreamOption::compressionType Foam::IOstreamOption::compressionEnum(const Foam::word&, Foam::IOstreamOption::compressionType)
in file db/IOstreams/IOstreams/IOstreamOption.C at line 115
Unknown compression specifier 'uncompressed', using compression off
Create mesh for time = 42

Reading pointDistributionDict

Reading quadratureProperties.populationBalance
--> FOAM Warning :
From virtual void Foam:robes::findElements(const Foam::fvMesh&)
in file probes/probes.C at line 230
Did not find location (0 0.5 0.5) in any cell. Skipping location.
Reading moment.0.populationBalanceMean
Reading moment.1.populationBalanceMean
Reading moment.2.populationBalanceMean
Reading moment.3.populationBalanceMean
Reading moment.4.populationBalanceMean
Selecting extendedMomentInversion: lognormal
Selecting univariateMomentInversion: Gauss


--> FOAM FATAL ERROR: (openfoam-2206)
The zero-order moment is negative.
Moment set: 5{-1e+300}

From void Foam::univariateMomentSet::checkRealizability(bool )
in file univariateMomentSet/univariateMomentSet.C at line 208.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::simpleExit(int, bool) at ??:?
#2 Foam::error::exiting(int, bool) at ??:?
#3 Foam::univariateMomentSet::checkRealizability(bool ) at ??:?
#4 ? at ??:?
#5 __libc_start_main in /lib64/libc.so.6
#6 ? at ??:?
Aborted (core dumped)



I get similar error on other case and data as well.
Error is same for all the cases with different initial moments.
--> FOAM FATAL ERROR: (openfoam-2206)
The zero-order moment is negative.
Moment set: 5{-1e+300}

Is there anything that I am missing out?