zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

lagrangian: SizeDistribution: Improve logging and prevent crash when cloud is empty

Browse Source
This commit is contained in:
Will Bainbridge
2025-04-15 10:11:28 +01:00
parent 86e126a7bc
commit 188910325f
1 changed files with 21 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

25
src/lagrangian/parcel/submodels/CloudFunctionObjects/SizeDistribution/SizeDistribution.C
View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org \\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2023 OpenFOAM Foundation \\ / A nd | Copyright (C) 2023-2025 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -32,10 +32,24 @@ License
template<class CloudType> template<class CloudType>
void Foam::SizeDistribution<CloudType>::write() void Foam::SizeDistribution<CloudType>::write()
{ {
const scalar d0 = this->owner().Dmin(), d1 = this->owner().Dmax(); // Check that there are some parcels
const label nParcels =
returnReduce(this->owner().nParcels(), sumOp<label>());
if (nParcels == 0)
{
Info<< type() << ": Not writing the distribution as the cloud "
<< "is empty" << nl << endl;
return;
}
// Quit if there is no distribution of diameters // Check that there is a non-zero range of diameters
if (d1 == d0) return; const scalar d0 = this->owner().Dmin(), d1 = this->owner().Dmax();
if (d1 == d0)
{
Info<< type() << ": Not writing the distribution as the cloud "
<< "has uniform particle diameters" << nl << endl;
return;
}
// The x-axis is linearly spaced between the limiting diameters // The x-axis is linearly spaced between the limiting diameters
scalarField ds(nPoints_); scalarField ds(nPoints_);
@ -77,6 +91,9 @@ void Foam::SizeDistribution<CloudType>::write()
parcelPDF.last() *= 2; parcelPDF.last() *= 2;
parcelPDF /= sum(parcelPDF)*(d1 - d0)/(nPoints_ - 1); parcelPDF /= sum(parcelPDF)*(d1 - d0)/(nPoints_ - 1);
Info<< type() << ": Writing the distribution to "
<< this->writeTimeDir() << nl << endl;
// Write // Write
if (Pstream::master()) if (Pstream::master())
{ {