ATHESIS
presented by
Guilhem DIF-PRADALIER
in candidacy for the degree of
Doctor of Philosophy of the Universite´ de Provence (Aix-Marseille I), France Recommended for acceptance by the Institut de Recherche sur la Fusion par confinement Magnetique´
Speciality: Plasma physics
First-principle description of collisional gyrokinetic turbulence in tokamak plasmas
Board of examiners for public PhD viva voce on October, 30th 2008:
Peter BEYER Professor Thesis advisor Jean JACQUINOT High Commissioner scientific counsellor Chair Yanick SARAZIN Physicist (HDR) CEA advisor Eric SONNENDRUCKER¨ Professor Invited member Laurent VILLARD Professor Referee Howard WILSON Professor Referee ar ne rien savoir, ce n’est rien, ne rien vouloir savoir Cnon plus, mais ne rien pouvoir savoir, savoir ne rien pouvoir savoir, voila` par ou` passe la paix dans l’ameˆ du chercheur incurieux.
Samuel Beckett, Molloy, 1951
out jadis. Jamais rien d’autre. D’essaye.´ De rate.´ TN’importe. Essayer encore. Rater encore. Rater mieux.
Samuel Beckett, Cap au pire, 1982 Contents
Acknowledgements / Remerciements i
Forewords iii
Notice: directions for use vii
1 Some introductory words on nuclear fusion, plasma physics and around 1 1.1 The basic physics of nuclear fusion ...... 2 1.1.1 What is a fusion plasma ? ...... 2 1.1.2 Basic principle: mass–energy equivalence ...... 4 1.2 When applied to human enterprises ...... 5
2 The eternal problem of representation: which model for the plasma ? 11 2.1 Browsing the models ...... 11 2.2 The kinetic model ...... 15 2.3 The fluid model ...... 18 2.3.1 Collective behaviour ...... 20 2.3.2 A parametric closure ...... 25 2.3.3 A closure based on entropy production ...... 29
3 Calculating the plasma behaviour as an initial value problem 35 3.1 A reduction of the kinetic description: the gyrokinetic model ...... 36 3.1.1 A special set of gyrokinetic equations to study plasma turbulence . . . . 36 3.1.2 Numerics and the second life of first-principle modelling ...... 45 3.2 Initial state v.s. equilibrium: the central role of the electric field ...... 49 3.2.1 Non-canonical initial distribution functions generate poloidal sheared flows 52 3.2.2 Canonical equilibrium in GYSELA ...... 54 3.2.3 Early self-organising dynamics due to initial conditions ...... 56 3.2.4 Saturation processes for the poloidal electric field ...... 59 3.2.5 Nonlinear evolution: impact on transport in a steady state of turbulence . 61 3.2.6 Discussion and conclusion ...... 64
4 Adding binary interactions to collective behaviour 67 4.1 Towards a kinetic “neoclassical” theory of collisions ...... 68 4.1.1 General requirements for a kinetic theory of collisions ...... 68 4.1.2 Some essential models of collisions ...... 72 4.2 Incorporating neoclassical theory into gyrokinetic modeling ...... 77 4.2.1 Main neoclassical results ...... 77 4.2.2 A minimal collision operator containing conventional neoclassical theory 81 4.2.3 Discussion: which physics can we address ? ...... 82 4.2.4 Recovering the neoclassical results ...... 84 4.3 Gyrokinetic turbulence calculated with neoclassical effects ...... 91 4.3.1 Collisional enhancement of turbulent transport –a cartography in (R/LT , ν?) space ...... 92 4.3.2 Evidence of turbulent generation of poloidal momentum ...... 96
Conclusion 99
A Important features about Fourier–Laplace transforms 103
B Establishing the dispersion function 105
C The analytic continuation issue 107
D Establishing the Plemelj formulae 111
E Landau damping: the historical derivation 113
F Closures in real space 115
G Establishing the inverse Fourier Transform of sgn(k) 117
H The numerical impact of non locality 119
I Marginal stability of the drift kinetic instability 121
J Recovering the particle transverse fluid current 123
K Short-time dynamics governed by the curvature and grad-B drifts 125
L Analysis of the departure between the analytical prediction and simulations 127
M Onset of GAMs with qualitative agreement on the frequency 129
N Parallel flows, a robust saturation mechanism for the charge polarisation 131
O Numerical implementation of the collision operator in GYSELA 133
Bibliography 135
Index 140 Acknowledgements / Remerciements
UESOIENTICICHALHEUREUSEMENTREMERCIES´ ... Q . . . tout d’abord les membres du jury pour la disponibilite´ –elle rime ici avec amabilite,´ cel´ erit´ e´ et travail durant des conges–´ dont ils ont fait preuve : je pense a` Peter Beyer mon directeur de these,` jury le jour prec´ edant,´ jury encore quelques jours apres` ; Jean Jacquinot pour avoir accepte´ de presider´ la soutenance, malgre´ un emploi du temps charge´ et pour son appreciation´ declar´ ee´ du manuscrit ; Eric Sonnendrucker¨ pour etreˆ venu de loin ; mes rapporteurs pour leur relecture atten- tive, sous forte contrainte temporelle de ce memeˆ manuscrit : Laurent Villard pour son amabilite´ et un aller, puis retour immediat´ vers Lausanne ou` l’enseignement attendait ; Howard Wilson pour sa cravate et le rapprochement franco-anglais auquel il a consenti par sa venue appreci´ ee´ ; . . . l’Universite´ de Provence, en la personne d’Anne-Marie Bajetto pour un “carton-plein” de derogations´ accordees´ ; . . . le DRFC–IRFM, notamment en les personnes d’Alain Becoulet´ et Tuong Hoang, pour m’avoir accueilli ces trois annees´ de these,` quatre mois de stage, alors ignorant de plasma, pour leur sou- tien, leurs encouragements et leur cordialite,´ en entretien ou plus fortuite ; . . . honneur aux dames et a` Virginie (Grandgirard), si Precieuse´ pour Platine et Jade –pour ses nuits parfois, weeks-end souvent avec effraction chez des amis pour glanerˆ un WIFI, surtout vers La Fin : “jobs non chaines´ : Virginie re-soumet !”, pour ses gateaux,ˆ son sourire...... ET pour Glaciere` TM ; Chantal (Passeron) pour une grande aide informatique, malgre´ un oeil inquiet a` toutes mes, entrees´ ; Gloria (Falchetto) pour ses visites souriantes et son tiramisu –superbe ; Ma- rina (Becoulet),´ Emmanuelle (Tzitrone), Remi´ (Dumont), Maurizio (Ottaviani) pour avoir egay´ e´ le cafe´ ; Guido (Huysmans) pour sa patience ”ehhh, puis-je t’emprunter les noeuds 16 a` 32 de Norma ?” et ses indications vestimentaires pour le jury. Enfin, des remerciements particuliers a` Phlippe (Ghendrih) pour sa cravate le Jour J, son oreille attentive et ses tres` appreci´ ees´ discussions “mainstream”, et surtout autres ; George–Lazenby–Xavier (Garbet) pour la lueur d’inter´ etˆ gratifi- ante du 33e transparent ce memeˆ Jour J, parce qu’il a souvent des idees,´ qu’elles sont souvent tres` bonnes et qu’il les dispense avec classe ; Yanick enfin (Sarazin) qui comprend mon attachement au canard et aux ballons pas ronds, parce qu’il court vraiment tres` bien, memeˆ en montagne, memeˆ a` Varenna, mais surtout parce que j’ai beaucoup appris et qu’il a beaucoup donne.´ . . Pardon, pour le costard aussi, encore ce memeˆ Jour J ; . . . les thesard(e)s´ / post-doc cotoyˆ es´ ces trois annees.´ Encore honneur a` la dame : Chris- tine (Nguyen), l’“inspiratrice de mes jour(nee)s”´ , toute “patrie–science–et–gloire”, elle com- prendra ; Roland (Belessa), un autre “patrie–science–et–gloire”,, plus axe´ sur le “parfois brutal, toujours loyal” –j’espere` que tu es bien la` ou` tu es maintenant ; Stan’ (Pamela)–le Big Ben de 16h15, moi qui suis souvent / parfois / si peu en retard– pour son sourire, sa fraˆıcheur et pour ses cours de phonetique´ : Evaˆaˆaˆaˆ Pamelaˆaˆaˆaˆaˆ ; Emiel, Guillaume, Guigui, Manu, Paolo, Guido “der´ e-d´ er´ e”,´ Fred´ eric,´ Ale’, Alessandro, Sophie, Ozg¨ ur¨ ; des mentions speciales´ : a` mon ancien et regrette´ collegue` de bureau Nico’ (Dubuit) –non Christine, tu n’es pas visee–´ maˆıtre es-Linux,` docteur es-plasma` et l’air patibulaire. . . “mais presque” –quelqu’un qui fait aimer le Nord ; aux patouillades (Tamain) de l’adjudant, tres` appreci´ ees,´ trop breves,` memeˆ s’il ne fallait vraiment pas que l’on parle politique ; a` un grand footballeur (Nardon) ; a` un bon ami, longtemps un voisin, un–mec–bien–(Gerbaud) qui vit maintenant entoure´ d’eau, en milieu hostile (pour les papilles du moins), froid et blafardˆ –a` tres` bientotˆ au soleil de Californie ; et enfin a` mon coach-surfer pref´ er´ e,´ l’ami Steve (Jaeger), pour son dehanch´ e´ inimitable, ses rouflaquettes et ton intelligence –encore quelqu’un qui ferait aimer le Nord, et qui aime le Sud ; . . . Aix-en-Provence maintenant, charmante cite.´ Charmante par ses habitant(e)s et surtout ceux-la` : Nico, Hel´ ene,` Roy, Ben´ e,´ JC, Malika et Herve´ (Fontan) pour sa ”grande gu..le”, ses qualites´ de DJ et de tolerance´ , un autre Herve´ (Roche) pour sa sensibilite,´ son Mac, ses vins et surtout pour son amitie,´ Melissou,´ , Axel (Torre) “alors-la` mec je t’arreteˆ tout de suite” –prends soin de Manu, Seeeergiooooooo (Chiarenza) pour son tact dans la discussion en public, pour les ’xxx–Watts’, les ’xxx-Points’, pour ses anniversaires, parce qu’il a un coeur gros comme c¸a –c’est normal, c’est un sportif ! mais surtout pour Elsa ; Manu (Rouge),´ pour les “bocs-bocs”, pour tout ce que je ne peux pas ecrire´ ici, et bien bien bien plus encore –on se reverra l’ami, et enfin Elsa pour l’“avenir...a` venir” ? . . . mon lecteur –qui ne lira pas les 100 pages et quelques suivantes– et qui est dej´ a` fatigue,´ je le sens. Mais je ne peux oublier ma famille, son soutien, mes absences (prolongees),´ sa tolerance,´ l’Exemple vivant de ma mere,` de mon pere,` de ma soeur en devenir et de mon frere.` . . parce qu’il est mon Grand frere.` . . et que je compatis aux moments semblables de redaction´ que tu vis : courage.
ESREGRETS... D . . . pour la journee´ de RTT de Jean Jacquinot ce Jeudi 30 Octobre, pour n’avoir pas pu me servir de Glaciere` TM pour emmener outre Atlantique des produits civilises´ (foie gras, magret, aile, manchon, gesier,´ cou farc¸i –en un mot: canard, petite´ , aligot et vin) et enfin parce que trois annees,´ entoure´ de si bonnes gens, c’est bien court. Forewords
c An investment in knowledge always pays the best interest. d Benjamin Franklin (Richard Saunders), Poor Richard, 1732.
O PROFITABLY PRODUCE electric energy from thermonuclear fusion reactions, one must confine i.e. sustain far for thermodynamic equilibrium and over long periods of time a Trarefied reactive medium, house of an unrelenting turbulent motion. At thermonuclear tem- peratures, the ambient thermal energy is by orders of magnitude larger than the ionisation energy. The elemental charged constituents of this medium may be considered as separate; the medium itself is termed a plasma. Its confinement is allowed by strong magnetic fields, exerting a contain- ment pressure over it. These fields generate an intricate immaterial magnetic topology of which the plasma particles are aware. As charged constituents, they indeed flow along the magnetic field lines, which, while winding up in space, generate three-dimensional containment surfaces. In a first approach, should these magnetic surfaces remain closed and nested, particles would not es- cape the configuration otherwise than by colliding with one another. Also, collisions alone would allow for a heat conduction across these surfaces, leading to a net energetic loss for the plasma. A key concept to quantify the achievable quality of this confinement is the so-called energy confinement time τE, ratio of the energetic content of the plasma over the power losses. In a tokamak, these losses essentially occur by conduction across the confinement surfaces. The con- 2 finement time may therefore be understood as a diffusion time a /χT , where a is a typical size of the plasma and χT its thermal diffusivity. Experimental measurements in the past thirty years of actual heat fluxes in magnetically confined devices have led to the following situation: the measured diffusivities for the ions are at least one order of magnitude higher as compared to the collisional expectations, and two orders of magnitude higher for the electrons. Tradition found a name for the situation: these transport phenomena are anomalous. The blame for this detrimental situation is most certainly to be put on the small-scale turbu- lence which originates from the various small-scale plasma instabilities. These instabilities feed on the density, temperature, velocity or magnetic field gradients transverse to the confinement sur- faces. These gradients cannot but exist and stem from the very idea of magnetically confining over long periods of time a far for equilibrium reactive medium. A key question to be answered for present devices and predicted for future reactors is therefore the decay or sustenance of plasma density and energy content. Extrapolation techniques may be based on semi-empirical scaling laws; the idea being (i) to squeeze into a limited set of dimensionless parameters a hole range of the supposedly relevant physical phenomena. Then (ii) write critical quantities, such as the energy confinement time as a function of these parameters. This function is chosen such that it minimises the dispersion of experimental data gathered from existing machines. For instance, according to present interna- −2.83 −0.50 −0.10 0.97 −2.56 −0.55 −2.72 tional databases, τE scales as follows: τE ∝ ρ? β ν? M q κ . The predicted confinement time for Iter should be around 6.0s. This law is represented in Fig.1. Here, ρ? is the normalised gyroradius and provides information on the required size of the device. In- formation on the achievable pressure is provided by the β parameter, ratio of the plasma kinetic iv Forewords
Figure 1: Thermal energy confinement time for ELMy data in the ITER H–mode database.
pressure over the magnetic pressure –thus smaller than unity. The collisionality of the medium is contained in the ν? parameter, a normalised frequency of collisions. Also, M provides informa- tion on the ion mass and thus on the dilution of the plasma, q on the ‘winding-up’ of the magnetic field lines, on the compactness of the device and κ on the elongation of the plasma. A trend takes shape (iii) which allows one to foresee behaviours in ranges of parameters which actually do not exist –see the red Iter dot in Fig.1. Such laws are interesting as trends, are interesting for engineers and for design purposes. Genuine complexity (and beauty) of the underlying physics is not addressed though. Ab initio theories, while describing turbulence (ρ?), magnetohydrodynam- ics (β) or binary interactions (ν?) from first principles, may provide a comprehensive knowledge of the physics at play. They alone would validate the scaling law approach in the relevant foreseen ranges of parameters and e.g. quantitatively seek for enhanced confinement regimes. Above all, they do seek comprehension. The problem then arises of how to describe the plasma, this medium composed of a formidable number of interacting charged particles, across which large-scale electric and magnetic fields flow and characterised by unrelenting self-organising motions. Given a distribution of charges and currents, the dynamics of the electromagnetic field is well described by the Maxwell theory. On the other hand, as the fields fluctuate, they back react on the density of charges and currents within the plasma which then adjusts to the perturbation and further contributes to modify the fields. A self-consistent description of the dynamics of the plasma thus involves to conjointly calculate the Maxwell equations and the response of the plasma. Essentially, while only adopting distinct angles, two frameworks are widely considered so as to calculate this latter response. They rest on the idea that the macroscopic properties of the plasma may be inferred from the microscopic description of its elementary entities and proceed from one another by further projections. In that spirit, the considerable number of degrees of freedom of the system –associated with the very large number of particles– contracts into the statistical concept of distribution function, which in turn may be projected on a few thermodynamic macroscopic quantities such as temperature, density, pressure, heat flux or entropy. The temporal dynamics of the distribution function is described in the framework of the kinetic theory whereas the fluid theory focuses on the temporal evolution of these latter macroscopic quantities. These may all be deduced from projections of the distribution function in velocity space. Practically, a fluid model consists in an open set of equations for these thermodynamic quantities to which is adjoined an Forewords v auxiliary equation –the closure relation, which may be freely chosen such that it closes the latter set of equations and makes it self-consistent. The velocity information contained in a six-dimentional kinetic phase space is smoothed out in fluid models, though it may partly re-enter the model through an adequate ad hoc choice of the closure relation. As such, fluid models have special trouble self-consistently describing, e.g. weakly collisional systems, wave–particle collective resonant effects, local trapping phenomena, or finite orbit width effects. Nonetheless, it is worth trying to refine them by advanced closure choices or multi-fluid descriptions –though doing so one may easily reach a somehow inordinate level of complexity– since they are numerically less demanding and at present time they alone are able to simulate the plasma up to its confinement time. As a hot rarefied medium, a tokamak plasma is crucially affected by wave–particle resonances, trapping phenomena, orbit effects and is intrinsically low collisional. As such, any fully self- consisent description of its response is therefore kinetic, and includes collisions. Though they lead to usually smaller transport than the one generated by the turbulence, one may not conclude on collisions being negligible. The discovery of high-confinement regimes –which seem universal and reproducible features of fusion devices and as such are the foreseen operating scenari for Iter– has led to a renewed interest in collisional theories. These regimes are indeed characterised by a local strong reduction of the turbulent activity; in that case, the transport processes may get very close to being collisional. Binary interactions are also important features for a kinetic theory since they alone may provide a correct description of equilibrium flows. Such large-scale sheared flows are of special relevance for the saturation mechanisms of the turbulence, potentially leading to high-confinement regimes. At last, there seems to exist a non trivial interaction between collisions and turbulence in tokamak-relevant operating regimes. Through these processes, collisions might indirectly have a strong influence on the actual level of transport; elucidating this interplay is a subject of great current focus for first-principle modelling.
∴
This dissertation starts in chapter 1 with a comprehensive introduction to nuclear fusion, its basic physics, goals and means. It especially defines the concept of a fusion plasma and some of its essential physical properties. The following chapter 2 discusses some fundamental con- cepts of statistical physics. It introduces the kinetic and the fluid frameworks, compares them and highlights their respective strengths and limitations. The end of the chapter is dedicated to the fluid theory. It presents two new sets of closure relations for fluid equations which retain impor- tant pieces of physics, relevant in the weakly collisional tokamak regimes: collective resonances which lead to Landau damping and entropy production. Nonetheless, since the evolution of the turbulence is intrinsically nonlinear and deeply influenced by velocity space effects, a kinetic col- lisional description is most relevant. First focussing on the kinetic aspect, chapter 3 introduces the so-called gyrokinetic framework along with the numerical solver –the GYSELA code– which will be used troughout this dissertation. Very generically, code solving is an initial value problem. The impact on turbulent nonlinear evolution of out of equilibrium initial conditions is discussed while studying transient flows, self-organising dynamics and memory effects due to initial conditions. This dissertation introduces an operational definition, now of routine use in the GYSELA code, for the initial state and concludes on the special importance of the accurate calculation of the ra- dial electric field. The GYSELA framework is further extended in chapter 4 to describe Coulomb collisions. The implementation of a collision operator acting on the full distribution function is presented. Its successful confrontation to collisional theory (neoclassical theory) is also shown. GYSELA is now part of the few gyrokinetic codes which can self-consistently address the inter- play between turbulence and collisions. While investigating this interaction, new results on the generation of poloidal momentum are presented.
Directions for use:
Throughout the writing of this dissertation, my constant concern was to try to put into perspective practical aspects of the defence of a thesis with wider problems in physics that epitomise the elusive quality physicists would call beauty. Not pretending to epistemology, the openings of the different chapters start with general discussions; technicality is momentarily swept aside and bridges with more fundamental concepts are tried to be unveiled. As the reader though goes deeper into each chapter, he might notice the intricate links which exist between the ‘ethereal’ ideas and the more practical aspects of a day by day work. If this happens to be, my attempt was successful. The reader on the jump may at first reading skip the paragraphs in smaller fonts and concentrate on the ends of chapters. This manuscript was also written so as to be a working document. As such, I made a wide use of cross-references, appendices and index entries. At last, the reader is not required to be familiar with plasma physics, statistical mechanics or any of the epistemological disputes briefly mentioned throughout this manuscript to be able to read it; I have endeavoured to make this dissertation as self-contained as I could.
Chapter 1
Some introductory words on nuclear fusion, plasma physics and around
c Just a few months before, the Korean war had broken out, and for the first time I saw direct confrontation with the communists. It was too disturbing. The cold war looked as if it were about to get hot. I knew then I had to reverse my earlier position. If I didn’t work on the [H] bomb, somebody else would – and I had thought if I were around Los Alamos I might still be a force for disarmament. So I agreed to join in developing the H-bomb. It seemed quite logical. But sometimes I wish I were more consistent an idealist. d H. Bethe, reported by S. S. Schweber, In the Shadow of the Bomb: Bethe, Oppenheimer, and the Moral Responsibility of the Scientist, 2000.
c OW are we constituted as subjects of our own knowledge ? How are we constituted as subjects who exercise or submit to power relations ? How are we constituted as moral Hsubjects of our own actions ?” This is Foucault’s translation of Kant’s famous threefold question: “What can I know ? What should I will ? and, What may I reasonably hope for ?” The central question is the moral one, which makes Foucault add that Les Lumieres` or the Kantian Aufklarung¨ have to be considered not “as a theory, a doctrine, nor even as a permanent body of knowledge which is accumulating; it has to be conceived as an attitude, an ethos, a philosophical life in which the critique of what we are is at one and the same time the historical analysis of the limits that are imposed on us and an experiment with the possibility of going beyond them.” As a legacy from Les Lumieres` or from Aufklarung¨ , knowledge is now highly determinant of political relations and social life; as a legacy from recent history, modern physics –and especially physics related to massive energy production– is a cornerstone of such knowledge. With quantum theory, necessity intrinsically arose to criticise the knowledge being built, the mere presence of the physicist, or of its measuring apparatus, fundamentally influencing the result of the measurement. With nuclear physics came mighty power, and politics. This dissertation does not forget that it is part of a half-century long fostered research effort to understand, master and mass produce a reliable source of energy for peaceful purposes. As such, this research is prominently political and the researcher fundamentally in the world or, following Merlot-Ponty’s word, fundamentally “at the world”. Which is not, to a young mind, one of its less beauties. The present chapter is written for the non specialist and gives some comprehensive basic definitions and general ideas for a better comprehension of the whole dissertation. It especially provides the reader with the basic physics of nuclear fusion, as it is thought to happen in Nature 2 Some introductory words on nuclear fusion, plasma physics and around and by which means it is investigated on Earth.
1.1 The basic physics of nuclear fusion
S a basic physics concept, nuclear fusion essentially requires to beat the Coulomb repulsive Aforce between atomic particles so that the attractive, strong nuclear force can bind these into a new nucleus. Risking a definition:
Nuclear fusion (or simply fusion) denotes the elementary process, accompanied by the release or absorption of energy, by which multiple atomic particles join together to form a heavier nucleus. If the energy to initiate the reaction comes from accelerating one of these particles, the process is called beam–target fusion; if both nuclei are accelerated, it is beam– beam fusion. If the nuclei are part of a plasma sustained out of thermal equilibrium, one speaks of thermonuclear fusion, which represents today’s most promising approach for a controlled source of energy production on Earth.
Throughout this dissertation, fusion is only considered in this latter way and investigated through the lens of plasma physics. This field bubbles with a tremendous variety of phenomena (turbu- lence, intermittency, bifurcation, chaos) due to the intrinsic nature of the reactive medium –the plasma. Complexity further increases since these phenomena are also highly dependent on e.g. (i) the atomic composition of the plasma –deuterium, tritium, alpha particles, other– (ii) its geome- try, (iii) boundary conditions, (iv) external sources, etc. Inversely to self-sustained fission-like chain reactions, the essential physical problem with fusion consists in finding the means to sustain through time a far from equilibrium unstable turbulent plasma, i.e. tame its more violent insta- bilities, understand / control the transport processes which underlay its dynamics, and ultimately predict its statistical behaviour –essentially: predict its transport processes. Before proceeding, let us discuss the central concept of fusion plasmas.
1.1.1 What is a fusion plasma ? The fundamental forces in Nature tend to organise matter into structures, which are perennial as long as their binding energies are large as compared to ambient thermal energy. If this is not the case, these structures would ultimately decompose into their elementary charged constituents. By no means though, are these constituents free; they keep being fundamentally affected by each others’ electromagnetic fields, and at a lesser level by binary interactions between each others. Collective motions of great strength and complexity are intrinsic to such an ensemble, predomi- nantly characterised by the excitation of a very broad variety of collective dynamical modes. Such a medium develops an intrinsic magnetism, which is a mere consequence of the un- relenting motion of its interacting charges. In turn, magnetism creates anisotropy, meaning that the properties of this ensemble in the direction parallel to the magnetic field are different from those in the transverse direction. This anisotropy contributes to the multi-scale dynamics of such a medium, i.e. involves a broad range of spatially and temporally correlated different scales: (i) from the fast electron density plasma oscillations at ωp or the fast Larmor gyromotion ωc along a magnetic field line, to the intermediate bounce frequencies ωb of particles trapped in local inho- mogeneities of the magnetic field, collisional frequencies ν, and down to the inverse of the energy −1 confinement time τE which characterises the rate at which a system loses energy to its environ- ment. Concerning (ii) the spatial scales, they range from the nuclear level of the strong force λSF or the tunneling length λtunn. under the Coulomb potential barrier which allow fusion to occur up to the system’s own ones and to the large mean free path λmfp characteristic of weakly collisional media. An important intermediate scale λD for neutral plasmas is named after its first discoverer Debye and defines the typical length beyond which the medium is almost electrically neutral and 1.1 The basic physics of nuclear fusion 3
−1 ωp ωc,e ωc,i ωb,e ωturb ωb,i νii νei τE νee
5 1011 109 107 106 105 105 10 1 1 0.2 (s−1)
λSF λtunn. λD ρc,e δb,e ρc,i δb,i a R λmfp
10−15 10−12 10−4 10−4 5 10−4 10−3 5 10−3 1 3 103 (m)
the former collective effects dominant due to large-scale long-range electromagnetic interactions. This introduces the important property of quasi-neutrality, which is the macroscopic manifesta- tion of the existence of an electric sheath; at and above this Debye length, electric fields are low –explaining why λD is also often referred to as a ‘screening’ length– and the medium is electri- cally neutral, yet polarisable. Inversely, at shorter lengths1, back-pulling strong electric fields and binary collisional interactions are the rule. A further physical interpretation is deeply connected to collision phenomena and shall be discussed in much detail in the last chapter 4, p.72.
Such an ensemble is termed a magnetised thermonuclear plasma, or simply a plasma in the following of this dissertation. Though such a definition can be much refined, it has the practical benefit to provide the reader with a list of the essential characteristics I shall refer to in the remainder of this manuscript. A more general definition for a plasma is quite loose and difficult since it increasingly applies to media which increasingly deviate from the concept of gas (see Table 1.1), finally ending up formalising a distinct state of matter, of which a thermonuclear plasma can be seen as a fulfilment.
Table 1.1: Intrinsic complexity of the thermonuclear plasma state as compared to gas.
Property Gas Plasma Electrical Very low Usually very high (higher than copper at ther- Conductivity monuclear temperatures ∼ T 3/2) Independently One; all gas particles behave Two or three: electrons, ions, and neutrals; dis- acting species in a similar way, influenced by tinguished by their charge, they can often behave gravity and binary collisions independently, allowing for new instabilities Velocity Maxwellian (far from the Often non-Maxwellian; collisional interactions distribution boundaries) since collisions are weak, external forcing can drive the plasma are dominant; few hot tails of far from local equilibrium, with a significant pop- fast particles ulation of fast particles Interactions Binary; three-body collisions Binary: within a sphere of Debye; collective out- are extremely rare side with long-range interactions due to electric and magnetic forces.
This somewhat highly versatile state of matter seems to be widely prevalent throughout the universe which has the double merit to be extremely flattering to plasma physics, and quite im- possible to disprove (or verify). To our present knowledge, matter is always in the plasma state
1The phase space volume encompassed within this so-called Debye sphere still contains a sufficiently large number of particles to allow for a statistical approach. This result gives a further idea of the broad variety of scales involved. 4 Some introductory words on nuclear fusion, plasma physics and around whenever fusion is achieved. Following through the initial goal of this section to provide the reader with a comprehensive physical background, let us now discuss the basic physics of nuclear fusion.
1.1.2 Basic principle: mass–energy equivalence Quantum chromodynamics describe the way the strong interaction (also called strong force or colour force) holds quarks and gluons together to form nucleons (protons and neutrons). Inversely, the nuclear force2 (also called nucleon–nucleon interaction or residual strong force) refers to the force between two or more nucleons and is responsible for binding of protons and neutrons into atomic nuclei.
Nuclear binding energy is derived from the nuclear force and is the energy required to disassemble a nucleus into free unbound neutrons and protons, so that the relative distances of the particles from each other are infinite (or in practise far enough to postulate a scale separation) so that the nuclear force can no longer cause the particles to interact.
Given the famous mass–energy equivalence: E = mc2, E being the energy, m the rest mass and c the celerity of light, and since a bound system is at a lower energy level than its unbound con- stituents, its mass is less than the total mass of its unbound constituents. A common misconception of this relation would lead to the disturbing idea that mass could ‘disappear’ in an isolated system, which is false. Binding energy mass transforms into heat, and as such, the system needs to be thermodynamically ‘opened’ –cooled in other words– for the so-called ‘deficit’ mass to appear. This mass difference after binding may range from a small fraction of the total mass for systems with low binding energies –essentially the lightweight elements of the periodic table– to an easily measurable fraction for systems with high binding energies. Iron and nickel nuclei have the largest binding energies per nucleon of all nuclei. These results are summarised in Fig. 1.1. To illustrate this idea, let us recall that the atomic mass unit –1.000000 u– is defined to be 1/12 of the mass of a 12C atom. The atomic mass of an hydrogen atom is 1.007825 u; each nucleon in 12C has thus lost, on average, about 0.8% of its mass in the form of binding energy. From the latter discussion readily stems the following rule:
Nuclear transformations (fusion and fission) which tend to produce high binding energy nuclei –middle weight nuclei close to 56F e– are generally exoergic reactions.
Important for human energy production means, this principle also plays a central role in astro- physics, through the so-called process of nucleosynthesis, first intuited in the late 200s by Atkinson and Houtermans: Nucleosynthesis is the process of creating new atomic nuclei from preexisting nucleons. The very first process of primordial nucleosynthesis –nucleogenesis– is believed to have occurred during the quark-gluon plasma of the Big Bang. All subsequent nucleosynthesis of nuclei heavier than lithium or beryllium –including all carbon, oxygen, . . . – is believed to occur primarily in stars, either by nuclear fusion or fission. Without considering advanced capture processes of nucleons –e.g. S-process, R-process for neu- tron capture; Rp-process for rapid proton capture– or photodisintegration of existing nuclei, nucle- osynthesis could have hardly produced the vast observed quantities of nuclei heavier than nickel or iron, even in massive stars. The basis for these theories was coined by Hoyle in 1946. Be- fore him, Bethe (1936–1938) and von Weizsacker¨ (1938) successively pictured nuclear fusion
2It was thought before the advent of QCD that the nuclear force was a fundamental force directly acting on the nucleons. The verification of the quark model led to consider this phenomenon as a residual side-effect of the truly fundamental strong interaction acting directly on quarks and gluons inside the nucleons. 1.2 When applied to human enterprises 5
Figure 1.1: Average binding energy per nucleon (in MeV) against the number of nucleons in nucleus for abundant isotopes. A few important ones for the purposes of nuclear fusion and nuclear fission are marked (see Fig. 1.3), as well as iron 56F e, which sits at the top of this graph and cannot yield energy from either fusion or fission.
processes which could account for the heat, luminosity and gravitational balance of most main- sequence stars, including the sun. These latter processes, known as chains or cycles, are only able to create elements up to iron and nickel. Known as the ‘proton–proton’ chain (Bethe) or the CNO cycle (Bethe, von Weizsacker),¨ they are often refered to as hydrogen burning (see Fig. 1.2). They are complemented with a successive set of reactions involving heavier elements, namely helium burning, carbon burning, neon burning, oxygen burning and silicon burning.
1.2 When applied to human enterprises
c Quand maˆıtrisera-t-on la fusion nucleaire,´ en laboratoire ? En 1979, apres` dej´ a` 30 ans de recherche, l’ignition n’avait et´ e´ obtenue que dans la bombe a` hydrogene;` toutefois j’ecrivis´ dans Pour la Sci- ence que la fusion pourrait ˆetre obtenue en laboratoire en moins de dix ans et que la construction de centrales electriques´ a` fusion suiv- rait de peu. [. . . ] Vingt ans ont passe,´ les physiciens courent toujours apres` le Graal de la fusion, et l’ignition, pensent-ils, est encore pour dans dix ans. [. . . ] Je crois que la fusion thermonucleaire´ [l’ignition] sera declench´ ee´ en laboratoire dans moins de dix ans. d Gerold Yonas, Pour la Science, no252, 1998.
Artificial fusion on Earth cannot mimic celestial thermonuclear reactions since the critical mass for gravitational-induced fusion is approximately 1/10th of the sun’s mass, roughly explain- ing why lighter celestial bodies are not stars. A central problem in this human attempt rests on the question of confinement, since gravitation is irrelevant in the laboratory. 6 Some introductory words on nuclear fusion, plasma physics and around
Figure 1.2: The Bethe ‘proton–proton’ chain (a) and the Bethe–von Weizsacker¨ CNO cycle (b) (sometimes also referred to as the ‘carbon cycle’). The former is believed to be the main reaction chain in most mean-sequence stars lighter or comparable to the sun; the latter is thought to be dominant in stars more massive than the sun.
Confinement stands for the means by which a thermonuclear plasma is constrained into a chosen volume of phase space to allow, at times, the Coulomb repulsive force to be overcome and fusion processes to occur. It encompasses both a control of the loss channels of particles (geometric control) and a heat and energy control, the latter being absolutely crucial. NB: Most of the elementary fusion processes occur through quantum tunneling.
Building upon the nuclear transmutation experiments of Rutherford done a few years earlier, fu- sion of deuterium nuclei was first observed by Oliphant and Harteck in 1932. Research in fusion for military purposes began in the early 19400s, as part of the Manhattan Project, but was not suc- cessful until 1952. Research for civilian purposes remained scarce and classified until 1958 and the ‘Atoms for peace’ conference in Geneva. And with it a new lease on life for plasma physics. Through fifty years, many ideas on confinement strategies have arisen, former conventional ideas now being almost abandoned and inversely. The now ‘non conventional’ approaches in civilian world are essentially ‘cold’ or ‘locally cold’–based principles such as palladium-induced “cold” fusion, muon-catalised fusion, accelerator-based light-ion fusion, sonoluminescence-based “bubble” fusion, antimatter-initialised fusion and, to my knowledge, the latest-born: pyroelec- tric fusion. Most of the latter become irrelevant with respect to the criterion of a possible mass production of electric energy. More convenional methods involve thermonuclear plasmas and are essentially twofold: (i) a military-coloured (at least in the EU and the USA) inertial confinement strategy in which the driver is a laser, an ion or electron beam, a ‘Fusor’ or a ‘Z-pinch’ or (ii) a magnetically confined method which also has its own variations, e.g. the reversed-field pinch, field-reversed configuration, levitated dipole, Spheromak, Riggatron, Stellarator and Tokamak.
Tokamak is the Russian acronym for тороидальная камера с магнитными катуш- ками –Toroidalnaya kamera s magnitnymi katushkami, litterally “toric magnetic chamber”. The concept has been developed in the 500s in USSR, among other magnetic configurations. It is nowadays considered, along with the Stellarator, as the most promising approach for energy production. Its physics represents a cornerstone for a wide class of fundamental problems related to turbulence, chaos and magnetism. 1.2 When applied to human enterprises 7
Figure 1.3: Subplot (a)– Fusion reaction rate (average of cross-section times speed) as a function of temperature. The average is over Maxwellian ion distributions with the appropriate temperature (data from the NRL Plasma Formulary, 2006). Subplot (b)– The believed fusion process between deuterium and tritium, Eq.(1.1a); lithium is here, as in a reactor, so as to generate in situ the tritium.
Of course, with each of the latterly mentioned means of confinement, dedicated problems arise, but the basic nuclear reactions remain. Amongst the great variety of conceivable fusion reactions, preference will go3 to those with the highest reaction rates τ, i.e. reactions having the highest probability of fusion processes per unit of time and volume: τ = nAnBhσvi, hσvi being the probability of a fusion reaction for two reactant nuclei of respective densities nA and nB. Figure 1.3–(a) shows the fusion reaction rate as a function of temperature for three of the most common reactions. The maximal values of the respective reaction rates constrain for each reaction the optimal operating temperature. The deuterium (D) – tritium (T) reaction Eq.(1.1a):
[5He]? D + T −→ 4He(3.5 MeV ) + n(14.1 MeV ) (1.1a) n +6Li −→ 4He(2.1 MeV ) + T (2.7 MeV ) (1.1b) n +7Li −→ 4He + T + n − 2.5 MeV (1.1c) thus readily has both (i) a higher probability of fusion than other reactions commonly considered and (ii) at a lower temperature, which makes it the best elementary choice for practical use. Helium has an extremely low mass per nucleon and therefore is energetically favoured as a fusion product. As a noble gas, it is chemically almost inert and harmless in the case of a confinement accident. Tritium has a short radioactive half-life of 12.32 years, with a subsequent very low abundance, which requires its breeding from lithium4 to sustain the deuterium–tritium latter fuel cycle, Eqs.(1.1b) and (1.1c). In the former, the reactant neutron is supplied by Eq.(1.1a). which also produces most of the useful energy. The reaction with 6Li is exoergic, also providing a small energy gain. These processes are illustrated in Fig. 1.3–(b). The relative abundance of isotopes
3The straightforward practical requirements for fusion reactions consist in (i) being exoergic, (ii) involving low Z nuclei, (iii) being a two body reaction –three body collisions are too improbable (iv) conserving protons and neutrons since cross sections for the weak interaction are too small and (v) having at least two products to conserve momentum and energy without relying on the electromagnetic force. 4The supply of lithium is more limited than that of deuterium, but still large enough to supply the world’s energy demand at an extrapolated increased rate for thousands of years. 8 Some introductory words on nuclear fusion, plasma physics and around
6Li/7Li ≈ 0.08 also requires to consider the third reaction Eq.(1.1c), which is unfavourable by at least two means: it is endoergic and does not consume the neutron –with radio-protection problems associated, even though it also breeds tritium. The neutron and radioactivity problems open the way to advanced prospective fusion reactions. A criterion for a reaction being ‘advanced’ –as compared to the D–T reaction, is that it increasingly relaxes the two following constraints, regardless of the intrinsic difficulty of the reaction:
Neutronicity is the fraction of fusion energy released as neutrons and gives a precious indica- tion of the magnitude of problems associated with neutrons: radiation damage, biological shielding, remote handling and safety;
Electrical efficiency is an increasing function of the proportion of charged products since di- rect conversion to electricity by charge exchange is much more efficient (> 98%) –though difficult– than conventional electricity production via thermal cycles (< 50%).
In that spirit, the four most widely considered advanced fuel cycles read:
D + D −→ T (1.01 MeV ) + p+(3.02 MeV ) (50%) (1.2) −→ 3He(0.82 MeV ) + n(2.45 MeV ) (50%) D +3He −→ 4He(3.6 MeV ) + p+(14.7 MeV ) (1.3) p+ +6Li −→ 4He(1.7 MeV ) +3He(2.3 MeV ) (1.4) p+ +11B −→ 3 4He + 8.7 MeV (1.5)
Other cycles involving 3He would require its supply either from other nuclear reactions or from extraterrestrial sources, such as the surface of the moon or the atmosphere of the gas giant planets and are not considered here. Quantitative comparison between Eq.(1.1a) and Eqs.(1.2), (1.3) and (1.5) is displayed in Table 1.2.
Table 1.2: Increasingly ‘advanced’ fusion reactions with the criteria of decreasing neutronicity and increasing electrical efficiency.
Temperature Reaction rate Fusion Charged Reaction Neutronicity (keV) ≡ 1/Z energy products energy D + T 13.6 1 17.6 3.5 0.80 D + D 15 1 12.5 4.2 0.66 D + 3He 58 1/2 18.3 18.3 ∼ 0.05 p+ + 11B 123 1/5 8.7 8.7 ∼ 0.001
A man-induced fusion reaction should trivially (i) produce more energy than it is fuelled with and (ii) be self-sustained. To quantify these ideas, a useful concept is the fusion energy gain factor, usually expressed with the symbol Q and ratio of fusion power produced in a nuclear fusion reactor to the power required to maintain the plasma in steady state.
Ignition refers to the conditions under which a plasma can be maintained by fusion re- actions without external energy input. Therefore, the goal of ignition corresponds to infinite Q and is not a necessary condition for a practical reactor.
The historical work which provides a criterion to reach ignition for a power producing thermonu- clear reactor was derived by Lawson in 1955, published in 1957 [77] and specifies a minimum 1.2 When applied to human enterprises 9 value on the so-called ‘triple product’ of density n, plasma temperature T and energy confinement time τE: 2 12kB T 21 −3 nT τE ≥ ≈ 10 keV sm (1.6) EFP hσvi where kB is the Boltzmann constant, EFP the energy of the charged fusion products and hσvi the probability of a fusion reaction for two reactant nuclei (cf. page 7). This so-called Lawson criterion can be physically interpretated as the ratio of the energetic contents of the plasma to its losses (or sources in a steady state). In practice, reaching ignition is not an absolute constraint for a viable power plant. Finite values for the fusion energy gain factor close to Q = 20 are often thought to be enough for a reactor. In Iter, whose objective is to demonstrate the practical feasability of a fusion-based production of energy on an industrial level, the goal is to reach Q = 10 and then demonstrate that a significant fraction of the heating power can come from fusion heating power.
Breakeven refers to the condition Q = 1 and means that 20% (for D–T fusion) of the heating power comes from fusion reactions. Above Q = 5 the fusion heating power is greater than the external heating power.
Chapter 2
The eternal problem of representation: which model for the plasma ?
c If you shut your door to all errors truth will be shut out. d Rabindranath Tagore, Stray Birds, 1916.
HYSICS essentially remains a back and forth confrontation of reality-induced man-forged models to man-probed ‘intrinsic’ reality. The discussion of whether this intrinsic, far-off P reality –the “truth” of Tagore– could be unveiled or does even exist is beyond the scope of the present dissertation. Such also applies to the question whether this reality is “written in a math- ematical language” and best described by Numbers. As a working hypothesis, let us postulate that indeed an ‘intrinsic’ reality does exist which can be accurately mathematised using differential equations, should our fate be that it shall always remain an unreachable ideal. A more immedi- ate problem for the physicist –some would say more relevant– is to build upon the transduction through perception or experience everyone does operate from the natural world and generate a representation. Even though this transduction process is very likely man-dependent, the represen- tation can become generic and a model emerge. The present chapter aims at providing the reader with the basics of some essentials models –the so-called kinetic and fluid models– which have proven most fecund over the years.
2.1 Browsing the models
c It is fair to say that plasma turbulence lacks the elements of simplicity, clarity and universality which have attracted many re- searchers to the study of high Reynolds number fluid turbulence. In contrast to the famous example, plasma turbulence is a problem in the dynamics of a multi-scale and complex system. d P. H. Diamond, S. I. Itoh, K. Itoh, Series: Self-Organization and Structure Formation in Turbulent Plasmas, to be published.
Arguably, a model for the plasma should at least satisfy the following four primal conditions, regardless of any of its later technical details: (i) have steady states which are special solutions of the time-independant evolution equation. This allows one to describe e.g. Langmuir waves, to define a thermodynamical equilib- rium characterised by temperatures, energies, . . . 12 The eternal problem of representation: which model for the plasma ?
(ii) be linear in time and space, i.e. satisfy a superposition principle to allow for undu- latory phenomena (interference, diffraction, resonances. . . ) (iii) be first order in time so as to satisfy causality. In the kinetic framework –which is of primal interest to us throughout the remainder of this dissertation– a quantitative explanation of this idea can be found in section 2.2, page 17. It is deeply linked with the Born interpretation of quantum mechanics.
So as not to renounce trying to understand Nature (even as an unreacheable ideal) rather than being satisfied with the increasing power physics can provide over it is highly dependent on the applicability of a causality principle in physics. Following Grete Hermann, it could be enounced as: Nothing is produced in Nature which is not due to anterior processes – i.e. which does not follow these processes with necessity. This statement applies to all physical characteristics which could be acknowledged.
Such a principle does not overlap with the concept of physical determinism since it would hastily assimilate the possibility to know that causal chains do underlay all perceived phe- nomena and the possibility to formulate predictions on them from primal hypotheses. Deter- minism is not contained in the concept of cause and some authors wrongly foretold the death of causality because of the emergence of chance. Going even further, Cournot did show that strict determinism does not even identify with absolute prediction capabilities as those of a Laplace demiurge for instance. Being able to predict a succession of events is not contained in the idea of causality and is only a means by which causal chains postulated by a theory can be justified or discarded in the light of experiments. The causality principle does only emphasise the necessary suc- cession of events and not the human possibility to evaluate them beforehand.
(iv) satisfy a ‘Correspondence’-like principle, i.e. encompass the Boltzmann kinetic theory of gases in the low temperature and high collisionality limits or the Navier–Stokes fluid mechanics in the limit of a vanishing electromagnetic field and highly collisional medium.
The fundamental idea of projection underlying the whole following discussion is very simple, has been known for long and is very likely better enunciated elsewhere. The consequences and re- sulting concepts though are non trivial. May the expert reader forgive their unsatisfactory wording; the goal of the following discussion is to provide the reader with some very basic ideas to better apprehend the profound unity behind the apparent great amount of available distinct models. The following ten pages could easily be skipped at first reading.
Three traditional approaches While only adopting distinct angles, three great possible frameworks are often introduced. They rest on the idea that the macroscopic properties of the plasma can be inferred from the microscopic description of its elementary entities. Building upon this idea, the successive models which will be briefly discussed are increasingly simplified (this simplification can go well beyond the fluid de- scription as enunciated in the following pages). They are successively obtained from one another by further projections, i.e. (i) a further decrease of the number of free parameters necessary to describe central concepts, (ii) a further cristallisation around a few central concepts which in turn (iii) often get increasingly abstract. Following this idea of projection, a disputable goal of physics is perhaps to explore the difficult (and somehow contradictory) focal point of (i) being techni- cally increasingly simple and (ii) conceptually increasingly abstract so as to put in one, or a few dense concepts an increasing proximity with the experiment. It is in that spirit that by successive 2.1 Browsing the models 13 projections the immensely large number of degrees of freedom of the system –associated with its very large number of particles– contracts into the statistical concept of distribution function which in turn can be projected on a few thermodynamic macroscopic quantities (temperature, density, pressure, volume, entropy. . . ). New fundamental concepts appear, deeply connected to the definition taken for the system: is it isolated, closed or open? is its evolution irreversible? is coarse graining a relevant concept (only a few macroscopic variables from the detailed microscopic state of the system are relevant for its physical description)?
An open system is a system with input, which changes its behavior in response to conditions outside its boundaries. Inversely, an isolated system is a system without input, whose behaviour is entirely explainable from within its boundaries. Stricto sensu, only the universe is rigorously isolated; isolation is often postulated –in plasma physics, core plasmas which this dissertation is about are always assumed to have no interactions with e.g. the coil system which generates the magnetic field. This assumption is somehow crude since work must be exchanged between the plasma and the coils, should it only be to sustain the former out of equilibrium. Inbetween, a closed system can exchange heat and work, but not matter with its environment. Very generally, systems may be variously closed to matter or energy and/or to work, to information and/or to organisation. It is sometimes stated that systems closed to energy are autark, systems closed to information are independent and systems closed to organisation are autonomous. At last, as emphasised by Krippendorff, whether or not a system has outputs does not enter the distinction between open and closed systems. Systems without output are non-knowable by an external observer (e.g. black holes) and systems without inputs are not controllable.
These questions are fundamental since thermodynamics and statistical mechanics occupy centre stage in the edifice of modern physics. Since their inception in the 19th century, they have seen spectacular developments, extending their range of application far beyond what their founding fathers had envisaged. However, at the same time it has also been well-known ever since that both theories face serious conceptual difficulties. Many of these became most palpable when the relation between the two theories was studied. Thermodynamics was (and is) believed to reduce to statistical mechanics, but Loschmidt’s reversibility objection and Zermelo’s recurrence objection soon showed that the two theories do not enter into any simple reductive relationship. This idea will not be discussed any further though it very interestingly emphasises how crucial the problem of representation is. The first model which shall now be described is at the same time the more intuitive and the more untractable. Each particle is individually described and its evolution intrinsically coupled to all other particles. The complete dynamical description of a such system therefore requires to solve the N coupled equations for the N particles the system is composed of and thus embrace at any given time the 6N degrees of freedom of the associated phase space.
As shown e.g. by Poincare´ [94], the description of geometrical space requires a three- dimensional basis q = (e1, e2, e3), defining what is usually referred to as real space. Con- jointly, velocity space v = (v1, v2, v3) –or momentum space p in the general case– can be defined, which allows one to introduce the six-dimensional (q, p) phase space. First introduced by Gibbs in 1901, it is a dimensionaly minimal space where all the possible trajectories of a dynamical system are represented.
As a many body problem and though making use of the most up-to-date computational techniques, it remains untractable for a system involving over a billion particles –which is far from the require- ments of the physical world since, typically, N ∼ 1023 in Iter. 14 The eternal problem of representation: which model for the plasma ?
The other approaches are statistical. Let us first introduce the notion of a system, i.e. a given ensemble of plasma particles. Within such an ensemble, given the great complexity of the inter- actions exerted on e.g one of its sub-systems by all the other parts of it, the statistical approach postulates that during a sufficiently long time, this sub-system will meet all its possible states a great number of times. More precisely, let dqdp be the infinitesimal volume of phase space occu- pied during the interval of time δt by this sub-system. The statistical hypothesis says that during a sufficiently long interval of time T the complex and entangled nature of the phase trajectory of this sub-system will make it loopback this elementary volume of phase space dqdp a great number of times. When T goes to infinitiy, δt/T is finite and defines the elementary probability dP to find at an arbitrary time the sub-system in the elementary volume dqdp: δt lim = dP = f(q, p) dqdp (2.1) T →∞ T where the fundamental quantity f is a function of all phase space coordinates and represents the density of probability in phase space. It is the physical analogue of the mathematical concept of a measure. As such it is positive definite and satisfies the normalisation condition: R f dqdp = 1. More commonly, f is referred to as the distribution function of the system and does not depend on the initial state of any of its sub-systems when T is sufficiently large. Indeed, during this interval of time, any of the states our sub-system has been through can be called the initial state since statistically equivalent to all the others. This very important property allows one to find the statistical distribution function for sub-systems without solving the initial mechanical problem with initial values given for the whole system. The fundamental problem of statistical physics therefore reduces for any physical system to the determination of its distribution function. At this stage, the reader could be concerned the temporal dynamics has vanished within this statistical description. It has not in a most interesting way. To plainly explain this idea, let us shortly go back to classical mechanics. In this theory, the knowledge of its initial conditions entirely determine the later dynamics of a system. For systems with a great number of degrees of freedom, these initial conditions are incompletely known –and even though they were, their mere number would hardly allow for an analytical solution of the mechanical problem. Statistical mechanics allows one to go beyond these latter imperfections since its mere concern for a given system is to determine its behaviour on average and does not depend on an exhaustive calculation of all its possible phase space trajectories. By definition of the distribution function, any average quantity O¯ of O(q, p) is simply obtained as the ensemble average over all possible states: Z O¯ = O(q, p) f(q, p) dqdp (2.2)
This formulation allows one to avoid following in time the temporal evolution of the physical quantity O(q, p). Statistical physics proceeds as follows:
(i) time is frozen; at t = t0, a snapshot of the system is taken (ii) an ensemble average is performed over all possible states which allows one to define at t0 the macroscopic averaged quantity O¯ of the physical quantity O(q, p) (iii) time is advanced infinitesimally further to t0 + δt and then frozen again, etc. In this approach, the ensemble average is the fundamental concept to define a macroscopic dynam- ical function and is not to be related to any other concept considered more fundamental. To answer a possible source of puzzlement for the reader, let us emphasise that this is the most fundamental approach of statistical physics since the ergodic assumption is sidestepped and time re-enters the problem rather as a parameter than as a true dynamical variable. For more details on this special interpretation of statistical mechanics, see e.g. Refs.[2], vol.2, p.386 or [76], p.14. From a practical point of view though, it is often very difficult –if not impossible, see e.g. discussion p.50 and subsequent definition– to define average quantities or equilibrium states via 2.2 The kinetic model 15 ensemble averaging procedures (e.g. numerically, this would require to perform way too many costly almost identical simulations, see following chapters). So practically, one would favour time averages. To rigourously be able to do this, since one does want to follow in time the variations of the physical quantity O(q, p), it is possible to construct a time-dependant quantity O(t) such that: 1 Z T O¯ = lim O(t)dt (2.3) T →∞ T 0 since from Eq.(2.1), the statistical average is equivalent to the temporal average. Output quantities of e.g. numerical simulations are such constructed O(t) functions, which explains why such a temporal approach of statistical mechanics is still of some special relevance. To finish with gen- eral principles, let us also emphasise that inversely to univocal classical mechanics, the statistical approach has an inbuilt randomness which is not due to the very nature of the physical objects which are considered but only to the fact that the macroscopic quantities are obtained from a lim- ited set of data as compared to the complete mechanical problem. Nonetheless, in practice, given the√ typical sizes of the considered systems, this randomness is totally negligible (it decreases as 1/ N , N being the number of particles) and the predictions of statistical physics univocal. According to the latter discussion, a first-principle description of the plasma dynamics is achieved when writing down: (i) an evolution equation for the distribution function coupled to (ii) self-consistent evolution equations for the electromagnetic fields E and B The laws of electrodynamics coined by Maxwell –step (ii)– were the essential discovery of the XIXth century and since then remain a well-defined problem. On the other hand, a clear descrip- tion of step (i) i.e. an unambiguous characterisation of the response of the plasma is of more recent history and to some extent is still undergoing. Throughout the two following sections, I shall briefly describe with decreasing complexity the essential theories coined to characterise the response of the plasma.
2.2 The kinetic model
An interesting starting point is the Liouville equation which is valid in the absence of dissipation. In such a framework, it is easy to show using hamiltonian dynamics that (i) the number of systems dN in an infinitesimal region of phase space is constant through time and (ii) using the integral invariants of Poincare´ (e.g. see [49], p.394) the latter infinitesimal volume of phase space dV itself dN is constant. Therefore, the density of systems D = dV is also constant and its total derivative vanishes: dD = 0 ⇔ ∂ − [H, ] D = 0 (2.4) dt t which constitutes the Liouville theorem. The formal solution of this equation is:
−H(t−t0) D(t) = e D(t0) (2.5) where t0 is an arbitrary time (initial time). The hamiltonian H is the propagator and satisfies the Hamilton equations:
∂tq ≡ q˙ = ∂pH (2.6)
∂tp ≡ p˙ = −∂qH (2.7)
The exponential form of the solution readily stems from the causality principle, as will be empha- sised page 17. A hand-waving argument why evolution equations must be first order in time is the following: starting from t0 = 0 < t1 < t2, one must satisfy: H(t2 − t1) + H(t1) = H(t2), 16 The eternal problem of representation: which model for the plasma ? stating that the final state is statically independant of the chosen evolution path. This idea is deeply connected to the Maupertuis’ principle of least action. From this property, there exists a function ℵ(H) such that: D(t2) = ℵ [H(t2)] D(0) and D(t2) = ℵ [H(t2 − t1)] ℵ [H(t1)] D(0). Therefore: ℵ is the exponential function as in Eq.(2.5) and the Liouville equation is first order in time. An equivalent and perhaps more geometric way to understand this theorem rests upon volume transformations of a region D(t = t0) → D(t) of phase space. Let us introduce x = (q, p) t t 2 and g the one parameter group of transformation such that: g (x) =t→0 x + v(x) t + O(t ). The function v(x) represents a velocity in phase space and the volume D(t) transforms such t R that: D(t) = g D(0). The elementary volume V(t) of phase space reads: V(t) = D(t) Idτ = R t R t D(0) Iog Jgt dx = D(0) 1. det ∂g /∂x dx, where I represents the unit tensor and Jgt the Jaco- t 2 bian of the transformation, which takes the simple form: det ∂g /∂x =t→0 1 + t ∇ · v + O(t ). The rate of change of an elementary volume of phase space then reads: Z dV = ∇ · v dx (2.8) dt t=t0 D(t0) which is the so-called Reynolds lemma. Since v represents the velocity of the elementary domain
(q, p), the divergence ∇ · v = ∂pi p˙i + ∂qi q˙i vanishes and so does the right-hand side of the latter equation. In the light of Eqs.(2.4) and (2.8) comes the following interpretations of the Louville theorem. In the absence of dissipation:
the volume of phase space occupied by an ensemble of trajectories is con- served through time. the distribution function of an ensemble of particles is constant along its or equivalently1: phase trajectories.
Liouville’s Theorem shows that, for closed systems, the local density of microstates following a particle path through phase space is constant as viewed by an observer moving with the ensemble –i.e. the total or convective time derivative vanishes. Thus, if the microstates are uniformly dis- tributed in phase space initially, they will remain so at all times. Liouville’s theorem ensures that the notion of time average makes sense, but ergodicity does not follow from Liouville’s theorem. This theorem does not hold in the presence of dissipation and, strictly speaking, is only valid fo closed systems or on short periods of time as compared to the typical time of energy exchange with the outside. Therefore, this approach must be extended to more general systems. Different approaches have been considered to go beyond the Liouville theorem –the so-called Klimontovitch or Lenard–Balescu approaches, or the BBGKY hierarchy which starts from the Liouville theorem. I shall not describe these ideas here, a complete derivation can be found e.g. in [91]. Essentially, these approaches all converge to what is known as a kinetic equation:
N df X = C(f) ⇔ ∂ f + (q˙ ∂ + p˙ ∂ ) f = C(f) dt t i qi i pi i=1
⇔ ∂tf − [H, f] = C(f) (2.9) where H is the hamiltonian Eqs.2.4–(2.7) and C is the collision operator. Though at this stage H and C are still unspecified, Eq.(2.9) is the generic form for a kinetic equation2 Such an equa- tion is sometimes referred to as a master equation –or the Botzmann equation in physics– i.e.
1This conservation property of f along the phase space trajectories is an important property which I shall come back to in section 3.1.2, page 46, since it has a practical numerical use. 2It is possible to go beyond the already very general hypotheses presiding over the derivation of this equation and to write even more general non-Markovian evolution equations; this aspect is far beyond the scope of the present dissertation and will not be discussed here. 2.2 The kinetic model 17 a phenomenological set of first-order differential equations describing the time evolution of the probability of a system to occupy each one of a discrete set of states. In probability theory, this identifies the evolution as a continuous-time Markov process. Over the past century, this equation has given rise to a considerable amount of works, herewith giving rise to an almost as large amount of limit models, often named after their discoverers and according to which version of the collision operator it contains (see section 4.1.2, p.72 for further details on the collision operator). Some of which are of great significance and have provided at a given time new physical insight. (i) A key insight applied by Boltzmann was to determine the collision term resulting solely from two-body collisions between particles which are assumed to be uncorrelated prior to the col- lision. This assumption was referred to by Boltzmann as the Stosszahl Ansatz, and is also known as the molecular chaos assumption; (ii) The Fokker–Planck equation was first coined to statistically describe the Brownian motion of a particle in a fluid, focusing on the time evolution of the probability density function of the position of this particle. It has since been generalised to other observables, especially the distribution function, and in its latter form is the fundamental equation for a collisional plasma. According to the form the collision operator it contains, Eq.(2.9) is either referred to a the Bolt- mann equation or the Fokker–Planck equation. Since his 1931 work on Markov processes, this same equation is also known as the Kolmogorov forward equation. (iii) In hot plasmas, it is of wide use not to consider collisions since the mean free path of a particle is usually very large and its collision rate very low. This simplified model was first studied by Vlasov in the 500s. For obvious reasons, the so-called Vlasov equation df/dt = 0 is also sometimes referred to as the collisionless Boltzmann equation. (iv) This Vlasov equation, though simplified, is significantly difficult to handle and still beyond today’s most up-to-date computational possibilities. In strongly magnetised plasmas, as in the case of fusion, it is possible to further simplify the Vlasov equation into the gyrokinetic equation. Since the past ten years essentially, this latter approach has triggered a considerable amount of interest in the plasma community. It also is the central focus of this dissertation. A more precise description of this model is especially given in section 3.1. (v) As detailed in the last two chapters of this dissertation, Coulomb binary collisions can have a significant impact on the response of the plasma and therefore must be modelled3. The collisional gyrokinetic equation, when coupled to the Maxwell equations, provides the state-of- the-art first-principle description of the evolution of a magnetised plasma. All the mentioned equations do satisfy the latter four primal conditions: they all (i) have steady states f0: [H, f0] + C(f0) = 0 or simply [H, f0] = 0. Also, since (ii) the Poisson bracket and the collision operator are bilinear forms, the superposition principle is straightforwardly met. All equations are (iii) first order in time, thus satisfying a causality principle. This idea needs some more explanation; I show below the weak form: non(∂t) ⇒ non(% > 0) ⇒ non physical.
An equation of the form: ∂tf − [H, f] = C(f) has a positive definite probability density, which is not the case for equations with higher order time derivatives.
Main line of the proof: In action-angle coordinates (α, J), equation: ∂tf −[H, f] = C, averaged over all angles and integrated by part can be recast into a continuity equation: ∂t% + ∇ · j = D with D a dissipative term, % = hfiα and ∇ · j = ∂JhH∂αfiα. Since f is positive definite, % is positive definite. Inversely, higher order time derivatives lead to probability densities written as time derivatives of hfiα, which can be negative. Nota bene: this has some deep connections with the Born interpretation of quantum mechanics. Let Ψ? be the ? ? complex conjugate of Ψ, the quantum state vector. Then ∂tΨ = G ⇒ Ψ (∂tΨ) = Ψ G and ∂tΨ = G ⇒ ? ? ? ? ? ? ? ∂tΨ = G . Thus: (∂tΨ ) Ψ = G Ψ. Summing the two contributions: ∂t (Ψ Ψ) = Ψ G + G Ψ. With very reasonable assumptions on G, the latter equation can be recast into a continuity equation: ∂t% + ∇ · j = 0 with ? ? a positive definite probability density % = Ψ Ψ and a probability current j = ~/m = (Ψ ∇Ψ). Inversely, a 2nd ? ? second-order equation in time ∂ttΨ = K(t) leads to a probability density % ∼ (Ψ ∂tΨ + Ψ∂tΨ ) which can take negative values. A more complete discussion can be found in the works of Born.
? 3In particular, an open question is especially whether lim (collisionality → 0) ≡ (collisionality = 0)? In particular, there always exists an ‘effective collisionality’ in numerical models, see discussions pp.26, 39 and 74. 18 The eternal problem of representation: which model for the plasma ?
At last, the problem (iv) of the ‘correspondence’ principle is non trivial and most interesting since it provides a greater insight on the true richness embedded within the kinetic description of a magnetised plasma as compared to a fluid. I shall now detail this idea.
2.3 The fluid model
This latter kinetic approach is still both analytically and numerically sorely tractable. The third possible description proceeds from a further simplification of the above model. The so-called fluid approach projects the kinetic equation Eq.(2.9) on the infinite polynomial velocity basis 2 k 4 {1, v, v , . . . , v ,... } and focuses on the moments of the kinetic equation . The moment Mk of order k is an integral over velocity space of v⊗kf, where ⊗k denotes a tensorial product of order i ≤ k. The first moments are the density (k = 0), the flow velocity (k = 1), the pressure (k = 2) and the heat flux (k = 3): Z +∞ n = d3v f (2.10) −∞ 1 Z +∞ u = d3v vf (2.11) n −∞ Z +∞ ¯ ¯ 3 P¯ ≡ p I + Π¯ = m d v (v − u) ⊗ (v − u) f (2.12) −∞ m Z +∞ q = d3v |v − u|2 (v − u) f (2.13) 2 −∞ where p = T r(P¯ )/3 is the isotropic scalar pressure and Π¯ the viscous stress tensor which embeds the off-diagonal contribution of the pressure tensor and its anisotropy. The result of this projection is an infinite set of equations, formally written: Z 3 ⊗k h i ∀k ∈ N, d v v Eq.(2.9) = g(Mk−1, Mk, Mk+1) (2.14) in which every kth–order fluid equation couples the order k and k − 1 moments to the following one k + 1, leading to the infinite fluid hierarchy. Should this infinite set of equations be solved, the approach would be equivalent to the kinetic model. As a greater simplicity is seeked, the fundamental fluid problem lies on how to truncate this infinite hierarchy, i.e. close the set of Eqns.(2.14) at a finite order kc < ∞, yet retaining the fairest possible amount of physics. This is referred to as the ‘closure problem’ and is detailed in the following sections.
It is not by chance that the fluid terminology has here been chosen. This model is indeed very closely related to what people usually refer to as fluids, i.e. systems which satisfy some form of Navier–Stokes equation. The first odd moment equation of Eq.(2.9), describing the momentum flux: ¯ m∂t (nu) + m∇ · (nu ⊗ u) + ∇ · P¯ − en (E + u × B) = F (2.15) displays the analogue of the usual convective derivative of a fluid particle and the latter equation can be recast in the usual Navier–Stokes way (with no particle sources):
mn (∂tu + u · ∇) u = −∇p + en (E + u × B) + F + R (2.16)
4Other projective basis are possible, especially those based on Hermite or Laguerre polynomials. The convergence is faster in the sense that for a given accuracy as compared to the kinetic model, less moments are needed. However, the physical interpretation of these moments is less clear than with the velocity basis. This idea will not be disussed in further detail here. The interested reader could e.g. refer to Ref.[48], p.104 for Hermite polynomials and to Ref.[106] for the Laguerre projection basis. 2.3 The fluid model 19 where F is the collisional friction force: Z F = m d3v v C (2.17) and R contains the viscous stresses. Though it is not the usual means by which the Navier– Stokes equation is usually derived, Eq.(2.16) shows that the ‘correspondance’ principle is indeed satisfied. The kinetic statistical approach does embed the usual fluid conservation principles of mass, momentum and energy. The derivation of the equations of fluid mechanics is usually done via the Reynolds transport theorem which states that the changes of some intensive property L defined over a control volume D must be equal to what is exchanged through the boundaries of the volume ∂D plus what is created (or consumed) by sources and sinks S inside the control volume:
d Z Z Z dτ L = − dσ L u · n + dτ S (2.18) dt D ∂D D
From this relation, readily stems the following equation on L:
∂tL + ∇ · (Lu) = S (2.19) which is the generic form of a continuity equation. Taking L = n, the density, gives the first even moment equation of the kinetic equation Eq.(2.9), describing the evolution of density:
∂tn + ∇ · (nu) = S (2.20)
In fluid mechanics, fluxes are often written in an isotropic way: Df ∆n. Here arises an essen- tial difference between a hydrodynamic fluid and a magnetised plasma since the magnetic field p creates an intrinsic anisotropy, which favors a transverse transport, often written: S = D ∆⊥n (transverse and parallel refer to the direction of the magnetic field lines). The most elemental form of the Navier–Stokes equation is obtained when Eq.(2.19) is applied to momentum L = ρu, with ρ = nm. Behind this idea are the hypotheses that: (i) it is possible to define a so-called fluid particle –composed of all particles inside the control volume– which, at any given time, is (ii) at thermodynamical equilibrium. This elementary volume is highly collisional –the mean free path is very small (∼ 10−9 m) as compared to the size of this pseudo-‘particle’– and therefore local thermodynamical equilibrium is always postulated. Inversely, in a plasma, the Debye length is often comparable to the electron Larmor radius (see page 3). It is therefore delicate to define a plasma particle in the same way than in a fluid, though the closely related notion of quasi-particle is at the root of the Particle In Cell (PIC) approach –see p.46. It is also rather inaccurate to postulate that a plasma is always close to local thermodynamical equilibrium (for a dedicated discussion on how to define a kinetic equilibrium, see next chapter). To briefly summarise this discussion, as compared to fluid mechanics:
• the ‘plasma mechanics’ lack two essential elements of simplicity, namely:
– the high collisionality of fluids which allows one to define local thermodynamical equilibrium inside a ‘fluid particle’; – the isotropy of fluids –with the notable exception of highly stratified geophysical fluids– which may be linked to their low or vanishing electro-magnetic field;
• . . . while containing the Navier–Stokes nonlinearities and therefore hydrodynamic turbu- lence and Kolmogorov cascade theory. This is the classic paradigm to approach the ‘prob- lem of turbulence’ as enounced by Kolmogorov in 1941, namely the question of the spectral 20 The eternal problem of representation: which model for the plasma ?
energy transfer and distribution between scales5.
The theory of hydrodynamic turbulence puts most emphasis on the question of spectral transfer and excitation; plasma kinetic turbulence theory complements this approach while also featuring the central questions of (i) the back-reaction of small-scale turbulence on all scales of the dis- tribution function and (ii) the response of the plasma to a seed perturbation. The fluid models of plasma physics derived from this kinetic description describe plasmas in terms of smoothed quantities in real space. All the velocity information contained in a phase space kinetic descrip- tion is smoothed out, though it can partly re-enter the model through an adequate ad hoc choice of the closure relation. Especially, fluid models have special trouble self-consistently describing systems which are (i) weakly collisional (ii) far from thermodynamic equilibrium; problems arise when also trying to (iii) resolve wave–particle resonance effects (e.g. collective behaviour such as Landau damping), (iv) local trapping phenomena (crucial for a weakly collisional description (i), see chapter 4), (v) finite orbit width effects or (vi) capturing velocity space structures like beams or double layers. Nonetheless, fluid models remain of wide use since: (i) they are numerically less demanding –until now they alone would allow to reach the plasma confinement time, given the available numerical resources, (ii) they allow for an interesting ‘physical intuition’ and as such give considerable insight on plasma physics, (iii) they interestingly model magnetohydrody- namic instabilities (the plasma is treated as a single fluid governed by a combination of Maxwell and Navier–Stokes equations) which are crucial for plasma stability and at last (iv) they can be much refined by advanced closure choices or multi-fluid descriptions, though possibly reaching a somehow inordinate level of complexity.
Before moving to the true heart of this dissertation based on a kinetic description of the plasma, I shall now detail some advanced closure ideas which allow to incorporate the Landau damping kinetic effect within a fluid model. On this specific question, the interested reader could also refer to e.g. Refs.[28, 55]. The following section 2.3.1 aims to provide the reader with the physical characteristics of the Landau damping mechanism. A new fluid closure embedding this latter mechanism is proposed in specified limits in section 2.3.2 and associated appendices while section 2.3.3 tackes a similar problem from the standing point of the entropy: it is intented to derive a fluid closure which intrinsically satisfies an H–theorem. For further details on this theorem, see discussion p.73 and Eq.(4.7).
2.3.1 Collective behaviour In this section, let us consider the simple Landau model very commonly studied so as to emphasise some basic and fundamental properties of a plasma kinetic description. It deals with the collision- less dynamics of non relativistic electrons, embedded in a vanishing magnetic field. The problem is further assumed to be electrostatic. Such a supposition is usual in tokamak plasmas, where most of the turbulent transport is governed by the electric potential. In the end, the ions are supposed at rest, with a constant density n0. The electron distribution function f and the electric potential φ satisfy the following normalised set of equations:
∂tf + v ∂xf + ∂xφ ∂vf = 0 (2.21a) Z +∞ 2 ∂xφ = fdv − 1 (2.21b) −∞
5Navier–Stokes turbulence is one of the last great open problems of classical physics; it is also of great interest in a purely mathematical sense. Somewhat surprisingly, given the wide range of practical uses for the Navier–Stokes equation, mathematicians have not yet proven that in three dimensions solutions always exist (existence), or that if they do exist they do not contain any infinities, singularities or discontinuities (smoothness). These are called the Navier– Stokes existence and smoothness problems. The Clay Mathematics Institute has called this one of the seven most important open problems in mathematics, and offered a $1, 000, 000 prize for a solution or a counter-example. 2.3 The fluid model 21
q 0kB T where distances are normalised to the Debye screening length λ = 2 and time is nor- D 4πnee −1 malised to ωp = λ/vT . The problem is two-dimensional in phase space (x, v). It is nonlinear, due to the last term in the Vlasov equation. In the remainder of this section, we shall restrict our- selves to the linear characteristics of Eqs.(2.21a) and (2.21b). For historical reasons, this system has deserved much attention since Landau first showed that it captures some essential features of the plasma kinetic properties, namely the resonant interactions between waves and particles [75]. He has solved Eqs.(2.21a) and (2.21b) for small-amplitude perturbations near an homogeneous 6 Maxwellian equilibrium , characterised by a vanishing electric field φeq = 0. Doing so, the Vlasov equation constrains the equilibrium distribution function feq to depend on velocity v only: f(x, v, t) = feq(v) + f˜(x, v, t) φ(x, t) = φ˜(x, t) where hf˜i = 0 by definition; the brackets referring to time and space averages. Moreover, in the limit of small perturbations, namely for f˜ feq, the nonlinear term can be restricted to ∂xφ ∂vfeq, the nonlinearity ∂xφ ∂vf˜ being of the second order. In this case Eqs.(2.21a) and (2.21b) reduce to their linearised expressions:
∂tf˜+ v ∂xf˜+ ∂xφ ∂vfeq = 0 (2.22a) Z +∞ 2 ˜ ∂xφ = fdv (2.22b) −∞
Consistently with the φeq = 0 hypothesis, the electron and ion equilibrium densities are equal.
The Fourier–Laplace issue: an initial value problem Physically, the idea is to look at the time and space response of a small-perturbed plasma, given an initial homogeneous equilibrium. Mathematically, since Eqs.(2.22a) and (2.22b) are linear and feq only depends on velocity, plane waves of the form exp{i(kx − ωt)} are eigenvectors of the linear problem. Therefore, in this linear framework, the following Fourier expansion seems natural:
P P ˆ ı(kx−ωt) f(x, v, t) = feq(v) + ∗ ∗ f (v) e k∈R ω∈R k,ω P P ı(kx−ωt) (2.23) φ(x, t) = ∗ ∗ φˆ e k∈R ω∈R k,ω and in the linear framework, the intricate form for ∂xφ ∂vf: