From 1b262422621426bf04e76e1225980e68e57c94d5 Mon Sep 17 00:00:00 2001
From: schottenhamml <helen.schottenhamml@fau.de>
Date: Thu, 6 Feb 2025 09:43:03 +0100
Subject: [PATCH] Fix logarithm in Musker's law.
---
src/lbmpy/boundaries/wall_function_models.py | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/lbmpy/boundaries/wall_function_models.py b/src/lbmpy/boundaries/wall_function_models.py
index 4f627787..e1962673 100644
--- a/src/lbmpy/boundaries/wall_function_models.py
+++ b/src/lbmpy/boundaries/wall_function_models.py
@@ -157,7 +157,7 @@ class MuskerLaw(ImplicitWallFunctionModel):
def law(u_p, y_p):
arctan = sp.Float(5.424) * sp.atan(sp.Float(0.119760479041916168) * y_p - sp.Float(0.488023952095808383))
logarithm = (sp.Float(0.434) * sp.log((y_p + sp.Float(10.6)) ** sp.Float(9.6)
- / (y_p ** 2 - sp.Float(8.15) * y_p + sp.Float(86)) ** 2, 10))
+ / (y_p ** 2 - sp.Float(8.15) * y_p + sp.Float(86)) ** 2))
return (arctan + logarithm - sp.Float(3.50727901936264842)) - u_p
u_plus = velocity_symbol / self.u_tau[0]
--
GitLab