You should average all times. There is only one average.
Technically the answer is neither. What you are looking for is the population quantity in statistics. What you can obtain is a sample estimate of the population quantity. The real time average is from time going from negative infinity to positive infinity. |
Quote:
The time averaging must be theoretically extended up to T->Infinite. Therefore, in practice you have to compute from 0 to T, but T must be enough large. |
Quote:
This leads me to another question. Since k is a time varying parameter and involves the mean of the squared u' (UPrime2mean as OpenFoam calls it), does it mean that this mean is not computed as poulation average but instead just the average up until the specific time t that I am computing k for? |
Quote:
|
Quote:
So is the (u'u')_bar term a time varying parameter? |
Quote:
[f_fil(x,t)]_bar = F(x) |
UMean and UPrime2Mean are averages from 0 to T and both of these are sample estimates for population statistics. There are infinite number of possible such estimators. You only really care about the best estimator that you get in the very last time dir because it is averaged over the longest time duration. Getting the population value is not possible unless you have the entire population ensemble, which means running your simulation until T=>infinity.
Hence, it's nearly impossible to directly calculate the time-varying u' and time-varying k at runtime because you don't know (and can never know) the population statistics. You must calculate it until T=>infinity, save U at all the times you care about, and then post-process it a posteriori. If you waned to do this manually you would run your simulation with fieldAverage option set. Go to the very last time dir and find UMean. Copy (and overwrite) this latest UMean into all previous time dirs. Then do something like: Code:
foamCalc addSubtract U subtract -field UMean -resultName UPrime And then calculate k. |
Quote:
However, I am not sure how k can be time varying when it only depends on the 3 directional components of UPrime2Mean which are not time varying. What am I missing? |
Quote:
|
Quote:
Apologies, I did not really get what you meant by "[f_fil(x,t)]_bar = F(x)" but thanks for helping me recall that in LES it is a filtering operation instead of averaging. I completely missed that. Anyhow, according to Sagaut's book on LES for incompressible flow, the resolved turbulence kinetic energy is simply half of the squared of filtered velocity (page 51). Since the velocity data I have is for the given grid, the filtered velocity in my case is simply this velocity data. Does this sound like a valid statement? |
Quote:
Your LES database provided for sure the filtered velocity. From that you can deduce the k (unsteady) but also the approximation for the RANS (steady) velocity and the fluctuations. My previous notation means you can apply the time averaging to the filtered variable and the result is a steady field. |
The instantaneous k depends on the three direction components of u', not uPrime2Mean. UPrime2Mean is the root mean square of u', it is a statistic like UMean. k inherits the temporal characteristics of u'.
You can get the time-averaged k from UPrime2Mean, but not the instantaneous k. And UPrime2Mean is ,,, etc. |
Quote:
Quote:
Thank you FMDenaro and LuckyTran for your kind replied. I have eventually figured it out with your help and managed to calculate the TKE of my dataset. |
All times are GMT -4. The time now is 23:35. |