Skip to content
Snippets Groups Projects
Commit 1b262422 authored by Helen Schottenhamml's avatar Helen Schottenhamml
Browse files

Fix logarithm in Musker's law.

parent e7a8f3a3
Branches
No related tags found
1 merge request!185Fix Musker's law
Pipeline #73558 failed
Stage: pretest
Stage: test
Stage: docs
Hide whitespace changes
Inline Side-by-side
src/lbmpy/boundaries/wall_function_models.py
+ 1
− 1
View file @ 1b262422
...@@ -157,7 +157,7 @@ class MuskerLaw(ImplicitWallFunctionModel): ...@@ -157,7 +157,7 @@ class MuskerLaw(ImplicitWallFunctionModel):
def law(u_p, y_p): def law(u_p, y_p):
arctan = sp.Float(5.424) * sp.atan(sp.Float(0.119760479041916168) * y_p - sp.Float(0.488023952095808383)) 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) 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 return (arctan + logarithm - sp.Float(3.50727901936264842)) - u_p
u_plus = velocity_symbol / self.u_tau[0] u_plus = velocity_symbol / self.u_tau[0]
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment