This function evaluates q(k) (McComb, p61) by summing all wavevectors in a k-shell. Then we divide through by the volume of the box. More...
Go to the source code of this file.
Namespaces | |
Foam | |
Namespace for OpenFOAM. | |
Functions | |
void | writeEk (const volVectorField &U, const Kmesh &K) |
This function evaluates q(k) (McComb, p61) by summing all wavevectors in a k-shell. Then we divide through by the volume of the box.
Original source file writeEk.H
- to be accurate, by the volume of the ellipsoid enclosing the box (assume cells even in each direction). Finally, multiply the q(k) values by k^2 to give the full power spectrum E(k). Integrating this over the whole range gives the energy in turbulence.
Definition in file writeEk.H.