Skip to content
GitLab
Projects Groups Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Register
  • Sign in
  • dumux dumux
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 100
    • Issues 100
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 86
    • Merge requests 86
  • CI/CD
    • CI/CD
    • Pipelines
    • Jobs
    • Schedules
  • Deployments
    • Deployments
    • Environments
    • Releases
  • Packages and registries
    • Packages and registries
    • Container Registry
  • Monitor
    • Monitor
    • Incidents
  • Analytics
    • Analytics
    • Value stream
    • CI/CD
    • Repository
  • Wiki
    • Wiki
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
Collapse sidebar
  • dumux-repositoriesdumux-repositories
  • dumuxdumux
  • Merge requests
  • !2139

[rans] add density in all terms of the RANS models

  • Review changes

  • Download
  • Email patches
  • Plain diff
Merged Katharina Heck requested to merge feature/rans-include-density-changes into master May 26, 2020
  • Overview 14
  • Commits 9
  • Pipelines 0
  • Changes 29

This is to include variable density in the the 1-2eq turbulence models.


Explanation, using the K-omega model as an example:

Previously implemented (K):

\frac{\partial \left( k \right)}{\partial t}
+ \nabla \cdot \left( \mathbf{v} k \right)
- \nabla \cdot \left[ \left( \nu + \sigma_\textrm{k} \frac{k}{\omega}\right) \nabla k \right]
- (2 \nu_t \mathbf{S} \cdot \mathbf{S})
+ \beta_k k \omega = 0 

Introduced in this MR (K):

\frac{\partial \left( \rho k \right)}{\partial t}
+ \nabla \cdot \left( \rho k \mathbf{v} \right)
- \nabla \cdot \left[ \left( \mu +  \rho \sigma_\textrm{k} \frac{k}{\omega}\right) \nabla k \right]
- (2 \mu_t \mathbf{S} \cdot \mathbf{S})
+ \rho \beta_k k \omega = 0 

Previously implemented (\omega):

\frac{\partial \left( \omega \right)}{\partial t}
+\nabla \cdot \left( \mathbf{v} \omega \right)
-\nabla \cdot \left[ \left( \nu + \sigma_{\omega} \frac{k}{\omega}\right) \nabla \omega \right]
-\alpha \frac{\omega}{k} (2 \nu_t \mathbf{S} \cdot \mathbf{S})
+\beta_{\omega} \omega^2
-\frac{\sigma_d}{\omega} \nabla k \nabla \omega = 0

Introduced in this MR (\omega):

\frac{\partial \left( \rho \omega \right)}{\partial t}
+\nabla \cdot \left( \rho \mathbf{v} \omega \right)
-\nabla \cdot \left[ \left( \mu + \rho \sigma_{\omega} \frac{k}{\omega}\right) \nabla \omega \right]
-\alpha \frac{\omega}{k} (2 \mu_t \mathbf{S} \cdot \mathbf{S})
+\rho \beta_{\omega} \omega^2
-\rho \frac{\sigma_d}{\omega} \nabla k \nabla \omega = 0

  • add density in production limiter term The production limiter is also modified from:
20 \beta_K k \omega \text{  to  } 20 \rho \beta_K k \omega

  • test how much that changes results
  • This shouldn't effect any of the isothermal single component tests. All ni or nc tests should be affected.
  • All of the ni and nc tests are effected.
  • Update reference solutions.

Introduce this additionally for the:

  • kepsilon
  • low re kepsilon
  • one eq models

Fix failing kepsilon wall function test

  • fix kepsilon wall functions

relevant documentation: https://www.afs.enea.it/project/neptunius/docs/fluent/html/th/node99.htm#eq6.3.8bc


  • Check, if it is necessary to implement the additional term (MOVED TO #951 (closed))
- 2/3 \rho k \delta_{ij}

in the production term in the 2-eq models (see: Computational Fluid Dynamics, Versteeg& Malalasekera page 90)

\tau_{ij} = 2 \, \mu_t \, S_{ij} - \frac{2}{3}\rho k \delta_{ij}
Edited Oct 28, 2020 by Ned Coltman
Assignee
Assign to
Reviewers
Request review from
Time tracking
Source branch: feature/rans-include-density-changes