I am trying to simulate decaying homogeneous isotropic turbulence. As initial condition I want a divergence-free vector-field, i.e, $\mathrm{div} U = 0$.

How do I initialize random velocity field in an uniform grid such that above condition is satisfied?

  • $\begingroup$ Welcome to SciComp.SE! Are you sure you mean $\mathrm{div} U =0$, i.e., that the sum of the derivatives vanish? This can obviously happen without the single derivatives being zero. $\endgroup$ May 13 '16 at 16:45
  • $\begingroup$ Thank you Christian Clason. yes I mean it right, I want to create a set of incompressible initial velocity field such that divU = 0 $\endgroup$
    – verito
    May 13 '16 at 18:37
  • $\begingroup$ OK, such a vector field is called divergence-free; I have taken the liberty to edit your question to make that explicit -- I hope you don't mind (otherwise you can revert the edit). $\endgroup$ May 13 '16 at 19:12
  • 2
    $\begingroup$ And googling "random divergence-free" gives a few hits: tandfonline.com/doi/abs/10.1080/15326349.2012.699756, martian-labs.com/martiantoolz/files/DFnoiseR.pdf, citeseerx.ist.psu.edu/viewdoc/…. Does any of that help? $\endgroup$ May 13 '16 at 19:13
  • $\begingroup$ It seems that there are two aspects to your question that both Bill Barth and Doug Lipinski address. The first aspect is the purely mathematical one, how to make a initial velocity field divergence free, which Bill points out can be achieved, basically, by using a vector identity. I believe this is similar (or equivalent) to Helmholtz vector decomposition. The second aspect is a physical one, will the divergence-free initial velocity field have characteristics representative of homogeneous isotropic turbulence? Doug points to references on how to address this second aspect. $\endgroup$
    – Charles
    May 18 '16 at 3:24

A purely random initial condition will not resemble realistic turbulence and will take quite some time to reach a realistic decay state. If you want good results it's more complicated than just ensuring zero divergence. Usually the velocity field is generated in spectral space to give the proper spectral energy distribution and is then transformed to physical space to initialize the simulation.

Section 4.1 of the following paper outlines one way to do this:
Huidan Yu, Sharath S. Girimaji, Li-Shi Luo, DNS and LES of decaying isotropic turbulence with and without frame rotation using lattice Boltzmann method, Journal of Computational Physics, Volume 209, Issue 2, 1 November 2005, Pages 599-616, ISSN 0021-9991, http://dx.doi.org/10.1016/j.jcp.2005.03.022.

This section begins:

We conduct the simulations in a 3D periodic cube with various resolutions $N^3$. The initial incompressible homogeneous isotropic velocity field $\vec{u}_0$ (${\nabla} \cdot \vec{u}_0 = 0$) is generated in spectral space $\kappa$ with the following energy spectrum in a prescribed range $\kappa_{\min} \le \kappa \le \kappa_{\max}$ and a random phase (cf. details in [47]):
\begin{align} \hat{E}(\kappa,0) = \left\{ \begin{array} 0.038\kappa^m \exp(-0.14\kappa^2) & \kappa\in[\kappa_{\min},\kappa_{\max}],\\ 0 & \kappa\notin[\kappa_{\min},\kappa_{\max}]. \end{array} \right.~~~~~~~~~~~~(16) \end{align}

Then the velocity field is transferred to physical space. In what follows, we use $m = 4$ or $2$ in Eq. (16) to investigate the effect of $m$ on the energy spectrum and other quantities.


[47] T. Miyauchi, T. Ishizu, Direct numerical simulation of homogeneous isotropic turbulence – decay of passive scalar fluctuation, in Preprint vol. JSME No. 914-2, 1991, pp. 166–168.

You might also want to look at:
Consistent initial conditions for the DNS of compressible turbulence Ristorcelli, J. R. and Blaisdell, G. A., Physics of Fluids, 9, 4-6 (1997), DOI:http://dx.doi.org/10.1063/1.869152
which could also be used for incompressible turbulence by using only the first step in their procedure which involves generating an appropriate divergence-free velocity field.

That paper chooses a solenoidal (zero divergence) velocity field with spectrum $E(\kappa) = A\kappa^4\exp(-2\kappa^2/\kappa_p^2)$ and random phase components. $A$ determines the velocity amplitude. $\kappa_p$ determines the peak in the spectral energy and is set equal to $12$ in that paper.

Note that these two methods are equivalent with appropriate choices of parameters.

As for the specifics of generating a zero-divergence velocity with a desired energy spectrum, that question has been asked here but has not been answered. Perhaps someone can go back and answer that old question now.

  • $\begingroup$ Thank you for your reply Doug. This is exactly what I want to do "decaying homogeneous isotropic turbulence". Although the paper states it clearly about the generation of velocity field. I am facing difficulty in how to transfer velocity field from spectral space to physical space. $\endgroup$
    – verito
    May 19 '16 at 8:50
  • $\begingroup$ I am aware of the fact that it is related to Fast Fourier Transforms but not really sure how it is done exactly. If you can briefly point out how it is done? This would be great help!! $\endgroup$
    – verito
    May 19 '16 at 8:54
  • $\begingroup$ @verito That question has been asked here but not answered. This related question may help you a bit but asks about the other direction (computing the energy spectrum from the velocity field) scicomp.stackexchange.com/questions/21360/… $\endgroup$ May 19 '16 at 13:17
  • $\begingroup$ The second paper I linked also has some details, including the fact that the Fourier components of the velocity field must be aligned with the wave vector $\vec{\kappa}$ for the field to be divergence-free. $\endgroup$ May 19 '16 at 13:18

Create two random scalar fields $f$ and $g$ and set the velocity to: $$ u=\nabla f \times \nabla g $$ which is guaranteed to be divergence free.


If you are in 2D and if you want more physical setups, I suggest you consider potential flows.

There are various ways to construct these potential flows which are always divergence free and which satisfy certain boundary conditions.

In the wikipedia article, you find the power laws, that describes flows around plates, edges, or in corners.

Also, there are other generating functions, that may give you, e.g., the flow around an idealized airfoil.

  • $\begingroup$ Thank you Jan. I am in 3D and looking for DNS simulation of DHIT as Doug Lipinski has pointed out. $\endgroup$
    – verito
    May 19 '16 at 8:59

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.