From 1e6d8d73a0b7fa9da4143f612b8174e468f65b04 Mon Sep 17 00:00:00 2001
From: Kutalmis Bercin <kutalmis.bercin@esi-group.com>
Date: Mon, 9 Dec 2024 15:27:13 +0000
Subject: [PATCH] BUG: OwenRyley: correct the modelling constant (fixes #3278)
---
.../filmSeparationModels/OwenRyleyModel/OwenRyleyModel.C | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/regionFaModels/liquidFilm/subModels/kinematic/injectionModel/filmSeparation/filmSeparationModels/OwenRyleyModel/OwenRyleyModel.C b/src/regionFaModels/liquidFilm/subModels/kinematic/injectionModel/filmSeparation/filmSeparationModels/OwenRyleyModel/OwenRyleyModel.C
index aed0e6a7e48..3db260f1b82 100644
--- a/src/regionFaModels/liquidFilm/subModels/kinematic/injectionModel/filmSeparation/filmSeparationModels/OwenRyleyModel/OwenRyleyModel.C
+++ b/src/regionFaModels/liquidFilm/subModels/kinematic/injectionModel/filmSeparation/filmSeparationModels/OwenRyleyModel/OwenRyleyModel.C
@@ -195,7 +195,7 @@ tmp<scalarField> OwenRyleyModel::netForce() const
const scalar R2 = R1 + h[i];
// Inertial force (OR:Eq. 11; '2' is an exponent in the Eq.)
- const scalar Fi = -72.0/60.0*h[i]*rho[i]*magSqrU[i]*invR1[i];
+ const scalar Fi = -4.0/3.0*h[i]*rho[i]*magSqrU[i]*invR1[i];
// Body force (OR:Eq. 11)
const scalar Fb =
--
GitLab