CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM

turbulence in compressible cases

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   April 3, 2020, 05:56
Default turbulence in compressible cases
New Member
Edmund Teller
Join Date: Feb 2020
Posts: 12
Rep Power: 3
edmund32 is on a distinguished road
Hi all,
I've been having some trouble with understanding the turbulence models in the compressible case, I have a cell that is essentially a cuboid with a central pillar cut out of its central (long) axis, so a slice through it would look like:

 side view:
           ____  ____
           |   | |   |
           |   | |   |
           |   | |   |
           |   | |   |
           |___| |___|
top view:
                             |  ____  |
                             |  |  |  |
                             |  |__|  |
Excuse the ASCII art. So the central pillar patches are defined as walls with zeroGradient boundary conditions for the pressure and noSlip conditions for velocity, while the outer top and side patches are open boundaries, with zeroGradient boundary conditions for both the pressure and velocity.

So, the solver I've been running is essentially mhdFoam combined with rhoPisofoam, as in the paper by Carlos Xisto (, so there are magnetic fields applied to generate the velocities.
The issue I've been having is that the velocities always diverge: all of the turbulence properties I set in the /0/ folder (k,nut,omega for RAN omegaSST and k,nut for LES) seem to be calculated, so their values immediately starts to decrease from the internalField values I set on the first iteration and turbulence seems to have less of an effect. For the walls I have a kqRWallFunction, nutkWallFunction, omegaWallFunction set for k, nut, and omega respectively, and all have zeroGradient set for the outer patches, with internal fields set to 1e5, 9e-3, 1.11e7 values for the internalField of k,nut,omega respectively.

My question is how do I increase the effect of the turbulence on the velocities to damp these instabilities? I don't really understand where to set constants to increase the effect of turbulence since every variable (k,nut,omega) seems to be calculated, should I be setting the outer patches to fixedValue? Should I just increase the viscosity mu in thermophysicalProperties?
Thank you.

Last edited by edmund32; April 3, 2020 at 06:08. Reason: ASCII art got messed up
edmund32 is offline   Reply With Quote

Old   April 8, 2020, 04:44
New Member
Edmund Teller
Join Date: Feb 2020
Posts: 12
Rep Power: 3
edmund32 is on a distinguished road
bumping for response
edmund32 is offline   Reply With Quote


compressible komegasst, compressible les, mhdfoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
New Phase Compressible Turbulence Model for multiphase system vishal3 OpenFOAM Programming & Development 0 November 24, 2015 01:17
Compressible Isotropic Turbulence Lokanayaki Main CFD Forum 0 December 23, 2012 10:21
Compressible Turbulence Model Implementations shaun OpenFOAM Running, Solving & CFD 4 October 21, 2009 12:55
the homogeneous isotrotic compressible turbulence. loong Main CFD Forum 0 March 30, 2006 01:18
LES turbulence for compressible flow jana FLUENT 0 June 17, 2005 14:14

All times are GMT -4. The time now is 10:47.