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

Spalart-Allmaras and "Non-Conservative Diffusion"

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 15, 2013, 18:22
Default Spalart-Allmaras and "Non-Conservative Diffusion"
  #1
New Member
 
Peter Doings
Join Date: Sep 2012
Location: Portugal
Posts: 4
Rep Power: 13
PeterDoings is on a distinguished road
Hello.

I'm an undergrad student and I'm implementing the S-A model for incompressible flow in a FV solver;

I have access to STAR-CCM+'s documentation and there I found an entry explaining that they did not discretize explicitly one of parcels of the diffusion term but that they combined it with the conservative diffusion one. I tried to get to the formula on my own using the identity:

\nabla\cdot(\phi\nabla\psi)=\phi\nabla^{2}\psi + \nabla\phi\cdot\nabla\psi

And got a very similar result to STAR's:

\frac{1}{\sigma_\nu}[\int_S(\nu+\tilde{\nu})\nabla\nu \cdot dS + \int_V C_{b2} (\nabla \nu \cdot \nabla \nu)dV] = \frac{1+C_{b2}}{\sigma_\nu}  \int(\nu+\tilde{\nu})\nabla\tilde{\nu} \cdot dS -  \int_V\frac{C_{b2}}{\sigma_\nu}(\nu+\tilde{\nu})\nabla^2\tilde{\nu}dV

Except that in the documentation it appears as a \approx and the last term has, instead of the Laplacian, a \nabla\tilde{\nu}^2 which I assume is their mistake (is it?).
My questions are: a) is this common practice? and b) I'll treat the first term as a regular diffusion term with the diffusivity given by last iteration's \tilde{\nu}, but how can I best discretize the last one? Explicitly calculate the laplacian on a cell by cell basis, multiply by \nu+\tilde{\nu} and put the result in the right-hand side?

Any help apreciated.
PeterDoings is offline   Reply With Quote

Reply


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



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