|
[Sponsors] |
Help with understanding production term in RHS! |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
January 17, 2023, 13:47 |
Help with understanding production term in RHS!
|
#1 |
New Member
S Ram Kumar
Join Date: Jan 2019
Location: India
Posts: 12
Rep Power: 7 |
Hello there!
I am trying to understand the differences between the versions of the 2 eqn Wilcox k-omega model, over the years. I am using the NASA turbulence models site for this. My doubt here is rather trivial, related to the Production term of the TKE. More specifically, I am trying to understand the mathematical breakdown for how the P_{inc} - Production term for incompressible flow is arrived at, in the below notes page at point 4. https://turbmodels.larc.nasa.gov/noteonrunning.html It would be of great help if someone could explain the step-wise breakdown for the P_{inc} clearly. I tried but in vain ! Request to see attached image of my write up. Thanks in advance! |
|
January 17, 2023, 14:30 |
|
#2 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,764
Rep Power: 71 |
Quote:
Not sure about your doubt, are you asking about the role of the trace of the velocity gradient tensor, that is Div v ? Generally, for incompressible flow we use the symmetric gradient with zero trace. The isotropic term is included in a modified pressure. |
||
January 17, 2023, 18:31 |
|
#3 |
Senior Member
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,665
Rep Power: 65 |
There are a few ways to arrive at it but what you need is to derive the transport equation for k. The transport equation for k is a sub-equation that lives inside the total energy equation.
Probably a more straightforward way is to take the dot product of time-accurate velocity with the momentum equation (i.e. take navier-stokes and multiply it by velocity) and then take the Reynolds average. You can then isolate the turbulent kinetic energy from the resulting expression after you subtract the time-averaged bulk energy (the average momentum equation multiplied by the average velocity). A Reynolds stress term will appear and after applying your eddy viscosity model you will arrive at the result you have printed. It is theoretically very simple to extend the validity of the production term into a compressible formulation (the formula is identical), but you need an eddy viscosity model that holds under these compressible actions. |
|
January 17, 2023, 22:33 |
|
#4 |
New Member
S Ram Kumar
Join Date: Jan 2019
Location: India
Posts: 12
Rep Power: 7 |
Sorry if I confused you, that's not my doubt.
What i wish to understand is how Production term P=(tau_(ij))*(do u_i)/(do x_j) becomes 2*mu_t*S_(ij)*S_(ij) for incompressible flow. As a first step of this process, i understand clearly that (tau_(ij))=2*mu_t*(S_(ij)) for incompressible flow. Where S_(ij)=0.5((do u_i)/(do x_j)+(do u_j)/(do x_i)) But I'm not able to deduce how they arrived at the final expression for P_inc in step 4. Request to see the attached screenshot for your reference. https://drive.google.com/file/d/10Pf...w?usp=drivesdk |
|
January 17, 2023, 23:13 |
|
#5 |
Senior Member
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,665
Rep Power: 65 |
||
January 18, 2023, 04:36 |
|
#6 |
Senior Member
|
You seem to miss the very fundamental fact that a double contraction between a tensor and a symmetric tensor equals the double contraction between the symmetric part of the first tensor and the originally symmetric tensor.
|
|
January 19, 2023, 13:42 |
Got it!
|
#7 | |
New Member
S Ram Kumar
Join Date: Jan 2019
Location: India
Posts: 12
Rep Power: 7 |
Quote:
Then each term becomes individually representative in that way, as shown. |
||
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Q: How to access turbulent production term G from TKE in scalarCodedSource | Kbshariff | OpenFOAM Running, Solving & CFD | 4 | January 11, 2022 14:57 |
How to calculate Production and Dissipation term | ht2017 | OpenFOAM Programming & Development | 0 | October 7, 2021 23:14 |
Replacing production term in Spalart-Allmaras via UDF | radix | Fluent UDF and Scheme Programming | 0 | February 10, 2017 11:18 |
understanding the term: fvm::Sp(fvc::div(phi), h) | dominik_christ | OpenFOAM Programming & Development | 3 | May 7, 2011 08:14 |
Numerical method for the k-epsilon production term | carlslack | Main CFD Forum | 0 | January 23, 2003 16:13 |