From c6220c162e22c15f77dcf740dc7fdf81948af36c Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Jan=20G=C3=A4rtner?= <>
Date: Fri, 9 May 2025 13:14:53 +0100
Subject: [PATCH] BUG: RosinRammler: correct the order of scale and shape
parameters (fixes #3342)
Co-authored-by:Kutalmis Bercin <kutalmis.bercin@esi-group.com>
---
.../distributionModels/RosinRammler/RosinRammler.C | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/src/lagrangian/distributionModels/RosinRammler/RosinRammler.C b/src/lagrangian/distributionModels/RosinRammler/RosinRammler.C
index 8cd242a961f..13a41a138d7 100644
--- a/src/lagrangian/distributionModels/RosinRammler/RosinRammler.C
+++ b/src/lagrangian/distributionModels/RosinRammler/RosinRammler.C
@@ -100,13 +100,13 @@ Foam::scalar Foam::distributionModels::RosinRammler::sample() const
Foam::scalar Foam::distributionModels::RosinRammler::meanValue() const
{
// (C:Eq. 5)
- const scalar a = scalar(1)/lambda_ + scalar(1);
- const scalar qMax = pow(maxValue_/n_, lambda_);
- const scalar qMin = pow(minValue_/n_, lambda_);
+ const scalar a = scalar(1)/n_ + scalar(1);
+ const scalar qMax = pow(maxValue_/lambda_, n_);
+ const scalar qMin = pow(minValue_/lambda_, n_);
const scalar gMax = Math::incGamma_P(a, qMax);
const scalar gMin = Math::incGamma_P(a, qMin);
- return n_/(exp(-qMin) - exp(-qMax))*(gMax - gMin);
+ return lambda_/(exp(-qMin) - exp(-qMax))*(gMax - gMin);
}
--
GitLab