CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > CFX

Heat flux and k-e wall boundary conditions

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 5, 2009, 06:47
Default Heat flux and k-e wall boundary conditions
  #1
New Member
 
Join Date: Mar 2009
Posts: 8
Rep Power: 8
Blob is on a distinguished road
I have questions regarding boundary conditions (BCs) when the k-e model is used.

1. No-slip wall BC: What are the boundary conditions for the k-e equations here? I read somewhere k is zero at walls. But if either one is zero, then singularity occurs due to the \frac{\epsilon}{k} and \frac{k^2}{\epsilon} terms. But obviously the solver doesn't blow up or anything. So what is happening here?

2. Heat flux BC: When I specify a fixed heat flux at a wall boundary, the laminar flux equation is q=-k \textunderscore thermal\frac{dT}{dx}. k_thermal is the conductivity. But since the k-e equation uses a turbulent thermal conductivity in the energy equation, what is k_thermal in this case? For analysis, I need to know \frac{dT}{dx}=-\frac{q}{k \textunderscore thermal} at the boundary which is the reason I need the k_thermal.

Thanks a lot in advance for your input and thoughts.
Blob is offline   Reply With Quote

Old   June 5, 2009, 08:08
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,807
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Hi,

This is all discussed in the documentation. It is also discussed in any decent CFD textbook.

But in summary:
1) When you integrate the turbulence model to the wall then k is zero and e is undefined. That is why you cannot integrate k-e turbulence models to the wall without some special treatment of the e equation. k-w based models are good here as w id defined at the wall so this family of models can be integrated to the wall.

So rather than integrating k-e to the wall it is usual to apply wall functions. This means that you model the shear of the boundary layer outside all the complicated near-wall stuff rather than directly resolving it in your simulation. Now k and e have finite values and k-e now has a wall boundary condition.

2) The turbulence model assumes the conductivity enhancement is proportional to the viscosity enhancement. The proportionality constant is the Turbulent Prandtl Number. From this assumption wall function boundary conditions can be derived for temperature at the wall. Also for a basic introduction try http://en.wikipedia.org/wiki/Turbulent_Prandtl_number

In any CFX simulation with heat transfer you should have the wall heat flux. I think the default turbulent Prandtl Number is 0.7 so if you don't have the turbulent enhancement of conductivity in the results you can work it out from the turbulent viscosity and the Turbulent Prandtl Number.

Glenn Horrocks
ghorrocks is offline   Reply With Quote

Old   June 5, 2009, 09:15
Default
  #3
New Member
 
Join Date: Mar 2009
Posts: 8
Rep Power: 8
Blob is on a distinguished road
Thank you for the pointers, Glenn. I have been going through the manuals trying to piece together the k-e model/wall function/turbulent energy equation. Your comments clarify the missing links.

I'm interested in the BCs because I'm trying to do a model reduction of the governing PDEs (momentum, energy, k-e) by Galerkin projecting them onto subspaces created with the CFD data using proper orthogonal/singular value decomposition. If all works, I'm left with a few ODEs that dynamically approximate the transient flow. Incorporating the BC effects into the ODEs has been tricky with the k-e model as my laminar counterpart works beautifully while the turbulent one doesn't at all. If you are familiar with this model reduction method, please let me know. I'll do a separate post on this as well.
Blob is offline   Reply With Quote

Old   June 6, 2009, 07:42
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,807
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Hi,

In that case I think you need a basic textbook in turbulence modelling. "Turbulence modelling for CFD" by Wilcox is my standard text and comes highly recommended (although he does have a very strong bias to the k-w turbulence model but we can forgive him for that because he developed most of it)

Glenn Horrocks
ghorrocks is offline   Reply With Quote

Old   June 8, 2009, 05:10
Default
  #5
New Member
 
Join Date: Mar 2009
Posts: 8
Rep Power: 8
Blob is on a distinguished road
Great! Thanks a lot!
Blob is offline   Reply With Quote

Reply

Tags
boundary condition, heat flux, k-e, turbulence

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Heat Flux Wall Boundary Confusion. Joee FLUENT 1 August 21, 2010 12:20
No results for solid domain Gary Holland CFX 10 March 13, 2009 04:30
Concentric tube heat exchanger (Air-Water) Young CFX 5 October 6, 2008 23:17
Convective Heat Transfer - Heat Exchanger Mark CFX 6 November 15, 2004 16:55
Multicomponent fluid Andrea CFX 2 October 11, 2004 05:12


All times are GMT -4. The time now is 03:00.