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

Implementation details of SST k-omega model

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 1, 2008, 05:14
Default Implementation details of SST k-omega model
  #1
Niels
Guest
 
Posts: n/a
Hello,

I'm considering the SST k-omega model for my finite volume solver for the compressible Navier-Stokes eqs. Thinking about the implementation of this model, I came up with some details that I don't understand. I hope you can help:

(1) Boundary conditions:

-------------------------

I bought (yes, bought, I don't have access to journals) the AIAA paper 93-2906 by Menter which explains the model. It does not mention which boundary conditions to apply at various types of boundaries. I sought help in the CFD wiki article "Turbulence free-stream boundary conditions" but this seems to be more about how to specify initial conditions. In particular, I would like to know what to do at no-slip walls and inflows/outflows.

For the moment, my solver implements rather simple BC's. At no slip walls I specify zero flux of momentum and I intended to do the same for k and omega. At inflows and outflows I use a flux vector splitting in the spirit of Steger and Warming to convect ingoing and outgoing information according to the characteristics. I was also planning on using a similar approach for k and omega. But I lack the insight to assess if this is ok or just plain stupid.

(2) Computation of blending functions

--------------------------------------

Referring to the expressions here:

http://www.cfd-online.com/Wiki/SST_k-omega_model

there seems to me to be a potential for numerical problems whenever omega or y approach zero, because they occur as denominators (fx see F1). Would you normally in an implementation introduce some mechanism that detects values of omega close to zero and then switch to some better behaved expression (e.g. by replacing omega with a larger "safe" value) ?

For y (the wall distance) to actually become zero I would have to evaluate the blending function at wall points. Presumably, then, this means I should not do that. This makes sense in as much as the blending function should be 1 in the near wall region anyway, so there should be no need to evaluate it. It does make me wonder, however, if I really understand what is going on. After all, shouldn't the SST model be valid right down to the wall, and hence it should be possible to evaluate all its expressions at the wall ? Or is there some tacid assumption that whenever you need to evaluate at the wall you replace the blending function expression with a hard-coded 1 ?

Thanks, Martin
  Reply With Quote

Old   July 1, 2008, 09:22
Default Re: Implementation details of SST k-omega model
  #2
John
Guest
 
Posts: n/a
Just as a clue refer to Peric codes (consult his book too), there u find implementation of both of k-eps and k-omega, maybe helps.

peric code URL: ftp://ftp.springer.de/pub/technik/peric/peric.tar.gz

good luck
  Reply With Quote

Old   July 2, 2008, 02:32
Default Re: Implementation details of SST k-omega model
  #3
Niels
Guest
 
Posts: n/a
Thanks for your reply.

This code uses the original Wilcox model with wall functions though, so it's not directly applicable.

I'm still not sure what to do at walls and inlets/outlets. I spent a great deal of time examining what happens to omega as y approaches zero. It seems that omega approaches a value inversely proportional to the square of y. So now I'm thinking about just setting omega to this value at walls, and setting other turbulent quantities (k, and turbulent viscosity) to zero there. But again, I lack the insight to know if that makes sense-
  Reply With Quote

Old   July 3, 2008, 17:45
Default Re: Implementation details of SST k-omega model
  #4
couillu
Guest
 
Posts: n/a
Niels,

For Menter's papers (92 and 93), they are free : http://ntrs.nasa.gov/search.jsp

Regards
  Reply With Quote

Reply

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
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
Boundary Conditions for k omega SST dancfd OpenFOAM Pre-Processing 0 June 9, 2011 23:25
why SST model??? vijay CFX 3 May 22, 2008 03:05
user subroutine error CFDUSER CFX 2 December 9, 2006 07:31
user defined function cfduser CFX 0 April 29, 2006 10:58


All times are GMT -4. The time now is 02:26.