[PhD Quals] Plasma Moments Equations, vector calc

Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
2 replies · 3K views
erin85
Messages
8
Reaction score
0

Homework Statement



Derive the moments equations for a plasma. I specifically need help developing the 2nd moment equation (the Energy Balance). I can't get one of the terms...

Homework Equations



The Boltzmann Equation:
[tex]\frac{\partial f_{\sigma}}{\partial t} + \vec{v}\cdot \nabla f_{\sigma}+\frac{e_{\sigma}}{m_{\sigma}}(\vec{E}+\vec{v}\times\vec{B})\cdot \nabla_{v}f_{\sigma} = C_{\sigma} + S_{\sigma}[/tex]

How to take the moment: Multiply by zn(v), then integrate over velocity (d3v), where z2=[tex]\frac{1}{2}m(\vec{v}\cdot\vec{v})[/tex].

Other misc. definitions:

[tex]n_{\sigma}=\int f_{\sigma}(\vec{v})d^3v[/tex]
[tex]v_{\sigma}=\int f_{\sigma}(\vec{v})\vec{v}d^3v/n_{\sigma}[/tex]
and in general,
[tex](A)_{\sigma}=\int f_{\sigma}(\vec{v})A(\vec{v})d^3v[/tex]

Also, [tex]\vec{Q}_{\sigma} = (1/2)n_{\sigma}m_{\sigma}[(\vec{v}\cdot\vec{v})\vec{v}]_{\sigma} =[/tex]energy flow,
[tex]TrM_{\sigma}= n_{\sigma}m_{\sigma}(\vec{v}\cdot\vec{v})_{\sigma}[/tex] = scalar trace of momentum stress tensor.

The 2nd moment equation is:

[tex]\frac{1}{2}\frac{\partial}{\partial t}(TrM_{\sigma})+\nabla \cdot \vec{Q}_{\sigma} = n_{\sigma}e{\sigma}\vec{v}_{\sigma}\cdot\vec{E}+R^2_{\sigma}+S^2_{\sigma}[/tex].

Don't worry about the C's, S's, and R's... those definitions are just straightforward, but annoying to type.

One last thing: [tex]\vec{A}\cdot\nabla f = \nabla \cdot f\vec{A}-f\nabla\cdot\vec{A}.[/tex]

The Attempt at a Solution



I have no trouble getting the first two terms... but the one with the electric field in it is tricky. I get every term going to zero...

Just taking the part of the equation with E (not v cross B yet), and a little simplfication:

[tex]\int \frac{e_{\sigma}}{2}(\vec{v}\cdot\vec{v})(\vec{E}\cdot\nabla_v f_{\sigma})d^3v[/tex]

Grouping the v's and the E together, and using the vector identity above,
=[tex]\frac{e_{\sigma}}{2}\int d^3v [\nabla_v\cdot f_{\sigma}((\vec{v}\cdot\vec{v})\vec{E}) - f_{\sigma}\nabla_v\cdot((\vec{v}\cdot\vec{v})\vec{E})][/tex]

Breaking up the second term, just using the chain rule... I get a term that's del dot a dot product, which should go to zero, and a del dot E, and E is not a function of v, so should also go to zero... Maybe breaking this up via chain rule is not correct.

The first term I also feel goes to zero... the integral and the del basically cancel, so I evaluate what's left at +_ infinity.

Any help please? I am having a lot of trouble with this... and it took a freaking long time to type, so it would be really excellent to get an answer.
 
Physics news on Phys.org
erin85 said:

Homework Statement



Derive the moments equations for a plasma. I specifically need help developing the 2nd moment equation (the Energy Balance). I can't get one of the terms...

Homework Equations



The Boltzmann Equation:
[tex]\frac{\partial f_{\sigma}}{\partial t} + \vec{v}\cdot \nabla f_{\sigma}+\frac{e_{\sigma}}{m_{\sigma}}(\vec{E}+\vec{v}\times\vec{B})\cdot \nabla_{v}f_{\sigma} = C_{\sigma} + S_{\sigma}[/tex]

How to take the moment: Multiply by zn(v), then integrate over velocity (d3v), where z2=[tex]\frac{1}{2}m(\vec{v}\cdot\vec{v})[/tex].

Other misc. definitions:

[tex]n_{\sigma}=\int f_{\sigma}(\vec{v})d^3v[/tex]
[tex]v_{\sigma}=\int f_{\sigma}(\vec{v})\vec{v}d^3v/n_{\sigma}[/tex]
and in general,
[tex](A)_{\sigma}=\int f_{\sigma}(\vec{v})A(\vec{v})d^3v[/tex]

Also, [tex]\vec{Q}_{\sigma} = (1/2)n_{\sigma}m_{\sigma}[(\vec{v}\cdot\vec{v})\vec{v}]_{\sigma} =[/tex]energy flow,
[tex]TrM_{\sigma}= n_{\sigma}m_{\sigma}(\vec{v}\cdot\vec{v})_{\sigma}[/tex] = scalar trace of momentum stress tensor.

The 2nd moment equation is:

[tex]\frac{1}{2}\frac{\partial}{\partial t}(TrM_{\sigma})+\nabla \cdot \vec{Q}_{\sigma} = n_{\sigma}e{\sigma}\vec{v}_{\sigma}\cdot\vec{E}+R^2_{\sigma}+S^2_{\sigma}[/tex].

Don't worry about the C's, S's, and R's... those definitions are just straightforward, but annoying to type.

One last thing: [tex]\vec{A}\cdot\nabla f = \nabla \cdot f\vec{A}-f\nabla\cdot\vec{A}.[/tex]

The Attempt at a Solution



I have no trouble getting the first two terms... but the one with the electric field in it is tricky. I get every term going to zero...

Just taking the part of the equation with E (not v cross B yet), and a little simplfication:

[tex]\int \frac{e_{\sigma}}{2}(\vec{v}\cdot\vec{v})(\vec{E}\cdot\nabla_v f_{\sigma})d^3v[/tex]

Grouping the v's and the E together, and using the vector identity above,
=[tex]\frac{e_{\sigma}}{2}\int d^3v [\nabla_v\cdot f_{\sigma}((\vec{v}\cdot\vec{v})\vec{E}) - f_{\sigma}\nabla_v\cdot((\vec{v}\cdot\vec{v})\vec{E})][/tex]

Breaking up the second term, just using the chain rule... I get a term that's del dot a dot product, which should go to zero, and a del dot E, and E is not a function of v, so should also go to zero... Maybe breaking this up via chain rule is not correct.

The first term I also feel goes to zero... the integral and the del basically cancel, so I evaluate what's left at +_ infinity.

Any help please? I am having a lot of trouble with this... and it took a freaking long time to type, so it would be really excellent to get an answer.

Hmm... let's see...

[tex] \int d^3 v \frac{m}{2}v^2\frac{e}{m}\vec E \cdot \frac{\partial f}{\partial \vec v}[/tex]

[tex] =<br /> \frac{e}{2}\vec E\cdot\int d^3v v^2 \frac{\partial f}{\partial \vec v}[/tex]

[tex] =<br /> -\frac{e}{2}\vec E \cdot<br /> \int d^3 v f \frac{\partial}{\partial \vec v}v^2[/tex]

[tex] =-{e}\vec E\cdot \int d^3 v f\vec v[/tex]

Also, just wondering... what school are these quals for?
 
Hi, thanks for the response. Sorry, I didn't notice... physicsforum usually e-mails me when I get a response... What rule is that for the transition from the 2nd to 3rd line? Not quite sure I'm getting that step. Although it is almost 3 in the morning, and I've been studying for a while...

Also, this quals is for Georgia Tech, the nuclear engineering department.