CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Spalart-Allmaras and "Non-Conservative Diffusion"

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

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

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


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

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