Quantcast
  • Register
PhysicsOverflow is a next-generation academic platform for physicists and astronomers, including a community peer review system and a postgraduate-level discussion forum analogous to MathOverflow.

Welcome to PhysicsOverflow! PhysicsOverflow is an open platform for community peer review and graduate-level Physics discussion.

Please help promote PhysicsOverflow ads elsewhere if you like it.

News

PO is now at the Physics Department of Bielefeld University!

New printer friendly PO pages!

Migration to Bielefeld University was successful!

Please vote for this year's PhysicsOverflow ads!

Please do help out in categorising submissions. Submit a paper to PhysicsOverflow!

... see more

Tools for paper authors

Submit paper
Claim Paper Authorship

Tools for SE users

Search User
Reclaim SE Account
Request Account Merger
Nativise imported posts
Claim post (deleted users)
Import SE post

Users whose questions have been imported from Physics Stack Exchange, Theoretical Physics Stack Exchange, or any other Stack Exchange site are kindly requested to reclaim their account and not to register as a new user.

Public \(\beta\) tools

Report a bug with a feature
Request a new functionality
404 page design
Send feedback

Attributions

(propose a free ad)

Site Statistics

205 submissions , 163 unreviewed
5,082 questions , 2,232 unanswered
5,353 answers , 22,789 comments
1,470 users with positive rep
820 active unimported users
More ...

  Why is there a transverse projection operator in the Langevin forcing used to drive turbulence?

+ 2 like - 0 dislike
3162 views

In equation (32) of this paper, the Fourir transformed correlator (spectral density) of the solenoidal and homogeneous in space and time "bare" forcing is given by\(\langle f_{\alpha}^0(\omega,k) f_{\beta}^0(\omega',k')\rangle = D_0(k)\delta(\omega+ \omega')\delta(k+k') P_{\alpha\beta}\)

where the transverse projection operator is given by

\(P_{\alpha\beta} = \delta_{\alpha\beta}-\frac{k_\alpha k_\beta}{k^2}\)

$D_0(k)$ is the spectrum of the bare forcing, and $\alpha$ and $\beta$ are indices for direction in space.

I understand the general form of this correlator, but what I do not understand is why the transverse projection operator $P_{\alpha\beta}$ is there. Does this have something to do with the fact that the forcing is assumed to be solenoidel? I am also not exactly sure about how a "solenoidal forcing" in this context is defined anyway...

asked Sep 10, 2014 in Theoretical Physics by Dilaton (6,240 points) [ revision history ]

1 Answer

+ 5 like - 0 dislike

This is explained nicely in arXiv:1012.0461 [cond-mat.stat-mech] (section 2.1) in the case where there are no thermal effects. I leave out the thermal effects in the following. Since you know your paper better than me you should be able to generalise this.

The projection comes from the incompressibility of the fluid. The Navier-Stokes equations read

$$ \partial_t v^\alpha + v^\beta \partial_\beta v^\alpha = \nu \partial_\beta \partial_\beta v^\alpha - \partial_\alpha P + f^\alpha \, ,  \qquad \partial_\beta v^\beta = 0 \, .$$

$\nu$ is the kinematic viscosity, $\partial_t = \partial/\partial t$, $\partial_\beta = \partial/\partial x_\beta$ and upper indices indicate the components of the corresponding vector. Repeated indices are summed over.

The pressure can be eliminated from these equations by taking the divergence of Navier-Stokes equations,

$$\partial_\alpha \left(\partial_t v^\alpha + v^\beta \partial_\beta v^\alpha = \nu \partial_\beta \partial_\beta v^\alpha - \partial_\alpha P + f^\alpha\right) \, .$$

Taking incompressibility into account, this leads to

$$\partial_\alpha \partial_\alpha P = -\partial_\alpha v^\beta \partial_\beta v^\alpha + \partial_\alpha f^\alpha \, , \text{ or } \, P = - \frac{1}{\nabla^2} \left(\partial_\gamma v^\beta \partial_\beta v^\gamma - \partial_\gamma f^\gamma \right)\, .$$

Then Navier-Stokes equations are written as

$$ \partial_t v^\alpha + \left(\delta_{\alpha,\gamma}-\partial_\alpha \frac{1}{\nabla^2} \partial_\gamma \right) v_\beta \partial_\beta v^\gamma = \nu \partial_\beta \partial_\beta v^\alpha +\left(\delta_{\alpha,\gamma}-\partial_\alpha \frac{1}{\nabla^2} \partial_\gamma \right) f^\gamma \, , $$

after inserting the second expression for the pressure. Finally, since the velocity field is incompressible,

$$ \left(\delta_{\alpha,\gamma}-\partial_\alpha \frac{1}{\nabla^2} \partial_\gamma \right) v^\gamma = v^\alpha \, ,$$

we can define new fields as

$$ \hat{v}^\alpha = \left(\delta_{\alpha,\gamma}-\partial_\alpha \frac{1}{\nabla^2} \partial_\gamma \right) v^\gamma \, , \text{ and } \hat{f}^\alpha = \left(\delta_{\alpha,\gamma}-\partial_\alpha \frac{1}{\nabla^2} \partial_\gamma \right) f^\gamma \, ,$$

in terms of which Navier-Stokes equations become

$$ \partial_t \hat{v}^\alpha + \hat{v}^\beta \partial_\beta \hat{v}^\alpha = \nu \partial_\beta \partial_\beta \hat{v}^\alpha + \hat{f}^\alpha \, ,  \qquad \partial_\beta v^\beta = 0 \, .$$

The pressure term has been completely eliminated. The price to pay is that the velocity field is implicitly constrained to be incompressible and that the force is replaced by $\hat{f}$.

If you compute the correlator of $\hat{f}$ in Fourier space you get,

$$ \langle \hat{f}^\alpha(p) \hat{f}^\beta(q) \rangle = \int_{x,y} \text{e}^{i px + q y} \left(\delta_{\alpha,\gamma_1}-\partial_\alpha \frac{1}{\nabla^2} \partial_{\gamma_1} \right)  \left(\delta_{\beta,\gamma_2}-\partial_\beta \frac{1}{\nabla^2} \partial_{\gamma_2} \right)  \langle f^{\gamma_1}(x) f^{\gamma_2}(y) \rangle$$

$$= \left(\delta_{\alpha,\gamma_1}- \frac{p^\alpha p^{\gamma_1}}{p^2} \right)  \left(\delta_{\beta,\gamma_2}-\frac{q^\beta q^{\gamma_2}}{q^2}  \right) \int_{x,y} \text{e}^{i px + q y}  \langle f^{\gamma_1}(x) f^{\gamma_2}(y) \rangle$$

Finally, inserting

$$ \int_{x,y} \text{e}^{i px + q y}  \langle f^{\gamma_1}(x) f^{\gamma_2}(y) \rangle = \delta_{\gamma_1,\gamma_2} \delta(p+q) D_0(p) \, ,$$

and contracting the indices leads to your expression for the correlator of $\hat{f}$.

answered Apr 8, 2015 by Steven Mathey (350 points) [ no revision ]

Your answer

Please use answers only to (at least partly) answer questions. To comment, discuss, or ask for clarification, leave a comment instead.
To mask links under text, please type your text, highlight it, and click the "link" button. You can then enter your link URL.
Please consult the FAQ for as to how to format your post.
This is the answer box; if you want to write a comment instead, please use the 'add comment' button.
Live preview (may slow down editor)   Preview
Your name to display (optional):
Privacy: Your email address will only be used for sending these notifications.
Anti-spam verification:
If you are a human please identify the position of the character covered by the symbol $\varnothing$ in the following word:
p$\hbar$ysicsOverflo$\varnothing$
Then drag the red bullet below over the corresponding character of our banner. When you drop it there, the bullet changes to green (on slow internet connections after a few seconds).
Please complete the anti-spam verification




user contributions licensed under cc by-sa 3.0 with attribution required

Your rights
...