-
Notifications
You must be signed in to change notification settings - Fork 190
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Implementing k-epsilon #2278
Comments
This would be a major new feature for Oceananigans! The implementation would probably resemble what we've done for "CATKE", which has a prognostic TKE variable: I would vote to implement Burchard and Umlauf's "generic length scale" equation. In this formulation there's a TKE variable (like CATKE), and also a "length scale" variable, which can either be epsilon, omega, or Mellor-Yamada's "q^2 l". So we get all of those models with one implementation: https://www.ingentaconnect.com/content/jmr/jmr/2003/00000061/00000002/art00004 The implementation might also want to take into account / support a few other formulations like I think @simone-silvestri could be interested in this kind of thing. It's a major undertaking! I think we should refactor |
By the way, the |
Hi @iurt! Very nice that you are thinking to implement a new parametrization! To start I would try to implement the additional prognostic equations, decoupled from the main solution. Once that is done and validated we can couple them to the momentum equations through the definition of a new "eddy viscosity" which uses values from the newly defined tracers |
Hi,
I just would like to give my humble contribution opening this issue.
I might not be able to do that because I am just a newbie in Julia, but it would be very nice to have k-epsilon as an option for turbulence closure.
Questions:
Where to start?
Any idea on how to implement that?
The text was updated successfully, but these errors were encountered: