IAEA TM on Theory of Plasmas Instabilities, Trieste 2-4 / 03 / 2005Y. Sarazin 1 Kinetic & Fluid Kinetic & Fluid descriptions descriptions of interchange of interchange turbulence turbulence Y. Sarazin, E. Fleurence, X. Garbet, Ph. Ghendrih, V. Grandgirard, M. Ottaviani Association Euratom-CEA, CEA/DSM/DRFC Cadarache, France P. Bertrand LPMIA-Université Henri Poincaré, Vandœuvre-lès-Nancy, France
16
Embed
Kinetic & Fluid descriptions of interchange turbulence
Kinetic & Fluid descriptions of interchange turbulence. Y. Sarazin, E. Fleurence, X. Garbet, Ph. Ghendrih, V. Grandgirard, M. Ottaviani Association Euratom-CEA, CEA/DSM/DRFC Cadarache, France P. Bertrand LPMIA-Université Henri Poincaré, Vandœuvre-l è s-Nancy, France. c . fluid. kinetic. - PowerPoint PPT Presentation
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
IAEA TM on Theory of Plasmas Instabilities, Trieste 2-4 / 03 / 2005 Y. Sarazin 1
IAEA TM on Theory of Plasmas Instabilities, Trieste 2-4 / 03 / 2005 Y. Sarazin 5
2 neglected
Fluid Fluid descriptiondescription
2 first moments of Vlasov evolution of density & pressure
Closure:
IAEA TM on Theory of Plasmas Instabilities, Trieste 2-4 / 03 / 2005 Y. Sarazin 6
Different linear stability Different linear stability diagrammesdiagrammes
Threshold instability: *Tc (1 k2) d
kinetic
(1) (1+k2) d fluid
Vanishing relative discrepancy for large density gradients *n
at *n
0
IAEA TM on Theory of Plasmas Instabilities, Trieste 2-4 / 03 / 2005 Y. Sarazin 7
2 distributions at energies
Fluid Fluid F at 2 F at 2 energiesenergies
Constraint: same numerics to treat fluid & kinetic descriptions
Ensures dissipation at small scalesstability of
Equivalent to 1 closure
in the limit 1
cf [V. Grandgirard, 2004 & this conference]
IAEA TM on Theory of Plasmas Instabilities, Trieste 2-4 / 03 / 2005 Y. Sarazin 8
Adjustable linear Adjustable linear propertiesproperties
3 degrees of freedom in fluid: T0, and D
Linear fluid properties can mimic kinetic ones:Linear fluid properties can mimic kinetic ones:
1. Linear threshold:
*Tc d(1 k2)
F(T0, , D)
adequate choice
2. Unstable spectrum width(or maximum growth rate)
at *n
0
IAEA TM on Theory of Plasmas Instabilities, Trieste 2-4 / 03 / 2005 Y. Sarazin 9
(analogous to "Dimits graph" with the same code)
Non linear discrepancy: QNon linear discrepancy: Qflfl Q Qkinkin
Heat turbulent transport larger in fluid than kinetic
by orders of magnitude
Suggest non linear threshold (Dimits upshift ?)
Transition not understood: ZF unchanged
(amplitude & dynamics)
IAEA TM on Theory of Plasmas Instabilities, Trieste 2-4 / 03 / 2005 Y. Sarazin 10
(analogous to "Dimits graph" with the same code)
Zonal Flows DO NOT Zonal Flows DO NOT
explainexplain
the whole differencethe whole difference Larger turbulent flux when ZF artificially suppressed
Difference still present between kinetic & fluid (orders of magnitude)
Note similar T profiles
for similar fluxes
IAEA TM on Theory of Plasmas Instabilities, Trieste 2-4 / 03 / 2005 Y. Sarazin 11
Quantifying the departure from FQuantifying the departure from FMaxwell
Projection on the basis of Laguerre polynomials Lp
(standard approach for neoclassical transport)
with
CorrespondanceCorrespondance
kth Fluid moment Polynomes L1 … Lk
IAEA TM on Theory of Plasmas Instabilities, Trieste 2-4 / 03 / 2005 Y. Sarazin 12
2 fluid moments are not 2 fluid moments are not enoughenough
Ortho-normal basis Lp
Slow convergence towards 0
Suggest any fluid description
of the problem
should account for
high order moments Mk (k>2)
May explain why fluid & kinetic results are still different w/o ZFMay explain why fluid & kinetic results are still different w/o ZF
IAEA TM on Theory of Plasmas Instabilities, Trieste 2-4 / 03 / 2005 Y. Sarazin 13
Alternative closure: entropy production Alternative closure: entropy production ratesrates
SQL governed by QL transport
Closure fulfils 2nd principle
.
Weights WQL :
Kinetic: infinity of resonances
Fluid: only 2
Main ideas:
Fluid closure: Q ( P + Peq ) Toperator: r(ky) + i i(ky)
IAEA TM on Theory of Plasmas Instabilities, Trieste 2-4 / 03 / 2005 Y. Sarazin 14
Similar linear behaviour Similar linear behaviour w/o ad-hoc dissipation
Same threshold as in kinetic: *Tcfl kyvd(1+k
2) *Tc
kin
Stability of small scales: implies i / ky 0
Similar linear spectra what about non linear behaviour ?
IAEA TM on Theory of Plasmas Instabilities, Trieste 2-4 / 03 / 2005 Y. Sarazin 15
ConclusionConclusionss
Looking for adequate fluid closures: what degree of convergence kinetic-fluid is requested? (, spectrum, dynamics, …)
Same numerical tool applied to 2D interchange model
1st closure: weak departure from FM (Q P T)– Linear properties can be made comparable (D required)– Fluid transport Kinetic transport– Non linear upshift not captured by ZF only (Non linear upshift not captured by ZF only ( Dimits) Dimits)– Possible explanation: large number of fluid moments required
22ndnd closure: Q closure: Q (P (P P Peq eq ) T) T– Target: balance entropy production rates Target: balance entropy production rates – Linear properties are similar
IAEA TM on Theory of Plasmas Instabilities, Trieste 2-4 / 03 / 2005 Y. Sarazin 16