SEARCH · Engineering Papers

Results for “flow instability”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Crossed beam energy transfer in the presence of laser speckle ponderomotive self-focusing and nonlinear sound waves

Crossed beam energy transfer, CBET, is investigated by taking into account the speckle structure of smoothed laser beams that overlap in a plasma with an inhomogeneous flow profile. Using the two-dimensional simulation code Harmony, it is shown how ponderomotive self-focusing of laser speckles in crossing beams can significantly affect the transfer of energy from one beam to the other. The role of plasma flow in speckle self-focusing is investigated and revisited, in particular its consequences in terms of redirection and increasing angular spread of the laser beams due to beam bending and plasma-induced smoothing, respectively. In close-to-sonic flow, the onset of self-focusing in the beam speckle structure occurs at considerably lower beam intensities than expected for the case without flow. CBET and speckle self-focusing can, hence, occur together when two crossed beams with equal frequency resonantly exchange energy via their ponderomotively driven density perturbations flowing with sound speed. From the simulations, it is found that consequences of ponderomotive self-focusing can be expected above an average intensity threshold scaling as $I_L~ 2 × 10^{14} W cm^{-2}(λ_0/1 μm)^{-2}(T_e/ keV)$, with an impact on the spatial and temporal coherence of the transmitted light. The density perturbations due to the ponderomotive force of the crossing beams can locally be enhanced in self-focusing speckles, partly leading to shock-like structures. These structures eventually increase the effect of plasma-induced smoothing and are at the origin of the stronger angular spread.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Evolution and feedback of AGN jets of different cosmic ray composition

Jet feedback from active galactic nuclei (AGN) is one of the most promising mechanisms for suppressing cooling flows in cool-core clusters. However, the composition of AGN jets and bubbles remains uncertain; they could be thermally dominated, or dominated by cosmic ray proton (CRp), cosmic ray electron (CRe), or magnetic energy. In this work, we investigate the evolution and feedback effects of CRp and CRe dominated jets by conducting 3D magnetohydrodynamic simulations of AGN jet-inflated bubbles in the intracluster medium using the FLASH code. We present the evolution of their energies, dynamics, and heating, and model their expected cavity-power versus radio-luminosity relation (P cav –L R ). We find that bubbles inflated by CRe dominated jets follow a very similar dynamical evolution to CRp dominated bubbles even though CRe within bubbles suffer significantly stronger synchrotron and inverse-Compton cooling. This is because, as CRe lose their energy, the jet-inflated bubbles quickly become thermally dominated within ∼30 Myr. Their total energy stops decreasing with CR energy and evolves similarly to CRp dominated bubbles. The ability of CRe and CRp dominated bubbles to heat the intracluster medium is also comparable; the cold gas formed via local thermal instabilities is well suppressed in both cases. The CRp and CRe bubbles follow different evolutionary trajectories on the P cav –L R plane, but the values are broadly consistent with observed ranges for FR-I sources. We also discuss observational techniques that have potential for constraining the composition of AGN jets and bubbles.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Fractional Laplacian spectral approach to turbulence in a dusty plasma monolayer

This work presents an analytical investigation of anomalous diffusion and turbulence in a dusty plasma monolayer, where energy transport across scales leads to the spontaneous formation of spatially disordered patterns. Many-body simulations of 10,000-particle dusty plasma monolayers are used to demonstrate how the global dynamics depend on the statistical properties of the dust assembly for realistic laboratory conditions. We find that disorder due to variations in the dust size distribution and charge-driven nonlocal interactions resulting in anomalous dust diffusion are key factors for the onset of instabilities. The resulting dynamics exhibit features of inertial turbulence over slightly more than half a decade of scales proportional or smaller than the Debye shielding length. These processes are examined analytically using a recently developed Fractional Laplacian Spectral (FLS) technique, which identifies the active energy channels as a function of scale, disorder concentration, and features of the nonlocal-interactions. The predictions from the theoretical (spectral) analysis demonstrate agreement with the results from the many-body (kinetic) simulations, thus providing a powerful tool for the study of active turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Energetic Particles - microturbulence interaction thrust (EP SIWG white paper)

This white paper outlines a strategic approach to Energetic Particle (EP) research, aimed at applying the results to the future burning devices such as ITER. Our goal is to deepen the understanding of the interplay between plasma background microturbulence and diverse species of energetic ions in burning plasma devices. This strategy is poised to advance our knowledge in this crucial area of fusion science, paving the way for significant breakthroughs in plasma physics and fusion energy research. We consider the confinement of energetic ions such as auxiliary heating beams, minority ions from ICRH, and fusion product alpha particles. Recent studies have underscored the potential of EP/microturbulence interaction to enhance fusion plasma performance by stabilizing microturbulence. However, it also poses risks, such as exacerbating Alfvénic eigenmode (AE) instabilities, which could compromise the sustainability of plasma discharge by causing EP losses. The interplay between EPs and background microturbulence is evident in the phenomenon of effective pitch angle scattering, a crucial aspect of the quasilinear (QL) theory. This theory is integral to developing numerically-efficient yet comprehensive and self-consistent approaches, recently employed to investigate the relaxation of energetic particle populations in the holistic modeling of fusion-grade plasmas. The primary objectives of this white paper thrust encompass a dual focus. First is that we need to delve into the mechanisms through which microturbulence engenders effective pitch angle scattering. This investigation will entail the examination of the presence of known or self-consistently predicted spectra of modes accountable for microturbulence. In these explorations, electrostatic microturbulence serves as a logical initial stride towards achieving the trust's objectives. This endeavor is anticipated to yield formulations expressing the parametric dependencies of the effective pitch angle scattering frequency on variables such as thermal plasma electron and/or ion temperatures, as well as their respective thermal conductivities. In our second objective, we aim to delve into the intricate formation of zonal flow (ZF) structures amidst the complex interplay of microturbulence and Alfvenic eigenmodes (AEs). This endeavor poses greater challenges as we seek to unravel the macroscopic manifestations influenced by microturbulence, termed as zonal structures (ZS), stemming from microturbulence-induced ZF. Drawing from this understanding, we anticipate employing a QL approach to yield comprehensive insights into the distribution function of energetic particles (EP) within phase space. This method entails resolving the dominant multidimensional phase space diffusion processes while effectively averaging over the rapid ballistic responses. Nevertheless, substantial strides remain imperative to realize a comprehensive whole-device modeling framework. This entails meticulous verification and validation exercises against experimental observations, as well as rigorous benchmarking against theoretical frameworks and numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Two self-similar Reynolds-stress transport models with anisotropic eddy viscosity

Two Reynolds-averaged Navier-Stokes models with full Reynolds-stress transport (RST) and tensor eddy viscosity are presented. These new models represent RST extensions of the $k−2L−a−\mathscr{C}$ and $k−ϕ−L−a−\mathscr{C}$ models by Morgan. Self-similarity analysis is used to derive constraints on model coefficients required to reproduce expected growth parameters for a variety of canonical flows, including Rayleigh-Taylor (RT) and Kelvin-Helmholtz (KH) mixing layers. Both models are then applied in one-dimensional simulation of RT and KH mixing layers, and the expected self-similar growth rates and anisotropy are obtained. Next, models are applied in two-dimensional simulation of the so-called “tilted rocket rig” inclined RT experiment and in simulation of a shock-accelerated localized patch of turbulence. Here it is found that RST is required to capture the qualitative growth of the shock-accelerated patch, and an anisotropic eddy viscosity provides substantial improvement over a Boussinesq treatment for the tilted rocket rig problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Local step-flow dynamics in thin film growth with desorption

Desorption of deposited species plays a role in determining the evolution of surface morphology during crystal growth when the desorption time constant is short compared with the time to diffuse to a defect site, step edge, or kink. However, experiments to directly test the predictions of these effects are lacking. Novel techniques such as in situ coherent x-ray scattering can provide significant new information. Herein we present x-ray photon correlation spectroscopy (XPCS) measurements during diindenoperylene (DIP) vapor deposition on thermally oxidized silicon surfaces. DIP forms a nearly complete two-dimensional first layer over the range of temperatures studied (40–120 ⁢°C), followed by mounded growth during subsequent deposition. Local step flow within mounds was observed, and we find that there was a terrace-length-dependent behavior of the step edge dynamics. This led to unstable growth with rapid roughening (β>0.5) and deviation from a symmetric error-function-like height profile. At high temperatures, the grooves between the mounds tend to close up leading to nearly flat polycrystalline films. In conclusion, numerical analysis based on a (1+1)-dimensional model suggests that terrace-length dependent desorption of deposited ad-molecules is an essential cause of the step dynamics, and it influences the morphology evolution.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Fluctuation dynamo in a weakly collisional plasma

The turbulent amplification of cosmic magnetic fields depends upon the material properties of the host plasma. In many hot, dilute astrophysical systems, such as the intracluster medium (ICM) of galaxy clusters, the rarity of particle–particle collisions allows departures from local thermodynamic equilibrium. These departures – pressure anisotropies – exert anisotropic viscous stresses on the plasma motions that inhibit their ability to stretch magnetic-field lines. We present an extensive numerical study of the fluctuation dynamo in a weakly collisional plasma using magnetohydrodynamic (MHD) equations endowed with a field-parallel viscous (Braginskii) stress. When the stress is limited to values consistent with a pressure anisotropy regulated by firehose and mirror instabilities, the Braginskii-MHD dynamo largely resembles its MHD counterpart, particularly when the magnetic field is dynamically weak. If instead the parallel viscous stress is left unabated – a situation relevant to recent kinetic simulations of the fluctuation dynamo and, we argue, to the early stages of the dynamo in a magnetized ICM – the dynamo changes its character, amplifying the magnetic field while exhibiting many characteristics reminiscent of the saturated state of the large-Prandtl-number ( ${Pm}\gtrsim {1}$ ) MHD dynamo. We construct an analytic model for the Braginskii-MHD dynamo in this regime, which successfully matches simulated dynamo growth rates and magnetic-energy spectra. A prediction of this model, confirmed by our numerical simulations, is that a Braginskii-MHD plasma without pressure-anisotropy limiters will not support a dynamo if the ratio of perpendicular and parallel viscosities is too small. This ratio reflects the relative allowed rates of field-line stretching and mixing, the latter of which promotes resistive dissipation of the magnetic field. Finally, in all cases that do exhibit a viable dynamo, the generated magnetic field is organized into folds that persist into the saturated state and bias the chaotic flow to acquire a scale-dependent spectral anisotropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Exceptional fracture toughness of CrCoNi-based medium- and high-entropy alloys at 20 kelvin

Medium- and high-entropy alloys based on the CrCoNi-system have been shown to display outstanding strength, tensile ductility and fracture toughness (damage-tolerance properties), especially at cryogenic temperatures. Here we examine the J Ic and (back-calculated) K JIc fracture toughness values of the face-centered cubic, equiatomic CrCoNi and CrMnFeCoNi alloys at 20 K. At flow stress values of ~1.5 GPa, crack-initiation K JIc toughnesses were found to be exceptionally high, respectively 235 and 415 MPa√m for CrMnFeCoNi and CrCoNi, with the latter displaying a crack-growth toughness K ss exceeding 540 MPa√m after 2.25 mm of stable cracking, which to our knowledge is the highest such value ever reported. Characterization of the crack-tip regions in CrCoNi by scanning electron and transmission electron microscopy reveal deformation structures at 20 K that are quite distinct from those at higher temperatures and involve heterogeneous nucleation, but restricted growth, of stacking faults and fine nanotwins, together with transformation to the hexagonal closed-packed phase. The coherent interfaces of these features can promote both the arrest and transmission of dislocations to generate respectively strength and ductility which strongly contributes to sustained strain hardening. Indeed, we believe that these nominally single-phase, concentrated solid-solution alloys develop their fracture resistance through a progressive synergy of deformation mechanisms, including dislocation glide, stacking-fault formation, nano-twinning and eventually in situ phase transformation, all of which serve to extend continuous strain hardening which simultaneously elevates strength and ductility (by delaying plastic instability), leading to truly exceptional resistance to fracture.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Magnetic flux emergence in binary neutron star remnants

Using high-resolution AthenaK simulations of a twisted toroidal flux tube, we study the flux emergence of magnetic structures in the shear layer of a hot massive neutron star typical of a binary neutron star remnant. High-resolution simulations demonstrate that magnetic buoyant instabilities allow for emergence only for extremely large magnetic fields significantly exceeding 10 17 G, and more typical fields around 10 16 G are instead dominated by hydrodynamic effects. Because merger remnants tend to be stable against hydrodynamic convection, our work places strong limitations on the mechanisms by which massive binary neutron star remnants can produce the magnetically-driven outflows needed to power jets.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Power-transfer and fixed-point analysis of sawtooth simulations of a current-carrying stellarator

Power-transfer and fixed-point analysis of previous NIMROD simulations (Roberds et al. , Phys. Plasmas , vol. 23, issue 9, 2016, 092513) improved the understanding of the effect of 3D (non-axisymmetric equilibrium) magnetic fields on sawtooth oscillations in the Compact Toroidal Hybrid (CTH) experiment. Computing the locations of order-1 fixed points, their Greene's residues, and local values for the rotational transform results in a description of CTH sawteeth consistent with Kadomtsev's model. A power-transfer analysis quantifies the distribution of energy among toroidal Fourier modes and their nonlinear interactions. The Lorentz power transfer drives sawtooth growth, and it is unambiguously interpreted as the flow of energy from toroidal mode $n'$ to mode $n$ , catalysed by $\boldsymbol {B}_{n-n'}$ . It has been reported previously that the CTH sawtooth frequency increases with the 3D field strength. This is attributed to an increased growth rate of the internal kink that drives sawtooth oscillations. Here, 3D fields remove energy from the kink, eliminating the possibility that these fields are an additional energy source that drives growth. Instead, 3D fields catalyse energy transfer from large-to-small scales, where magnetic reconnection is stronger. It is proposed that this energy transfer increases the reconnection rate at small scales, which is consistent with the increased growth rate observed at higher 3D field strengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Growing neoclassical tearing modes seeded via transient-induced-multimode interactions

Nonlinear extended magnetohydrodynamics simulations demonstrating seeding of neoclassical tearing modes (NTMs) via magnetohydrodynamics-transient-induced multimode interactions are presented. Simulations of NTMs are enabled by two recent NIMROD code developments: the implementation of heuristic neoclassical stresses and the application of transient magnetic perturbations (MPs) at the boundary. NTMs are driven unstable by the inherently pressure driven kinetic bootstrap current, which arises due to collisional viscosity between passing and trapped electrons. These simulations use heuristic closures that model the neoclassical electron and ion stresses. NTM growth requires a seed island, which is generated by a transiently applied MP in simulations. The capability is demonstrated using kinetic-based reconstructions with flow of a DIII-D ITER Baseline Scenario discharge (La Haye et al., in Proceedings IAEA FEC, 2020). The applied MP seeds a 2/1 NTM that grows in two phases: a slow growth phase followed by a faster robust growth phase. Additionally, an evolving sequence of higher order core modes are excited at first. Power transfer analysis shows that nonlinear interactions between the core modes and the 2/1 helps drive the initial slow growth. Once the induced 2/1 magnetic island reaches a critical width, the NTM transitions to faster robust growth, which is well described by the nonlinear modified Rutherford equation. Furthermore, this work highlights the role of nonlinear mode coupling in seeding NTMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Nonlinear coupling of electromagnetic and electron acoustic waves in multi-species degenerate astrophysical plasma

Nonlinear wave-coupling is studied in a multispecies degenerate astrophysical plasma consisting of two electron species (at different temperatures): a highly degenerate main component plus a smaller classical relativistic flow immersed in a static neutralizing ion background. It is shown that the high frequency electromagnetic waves through their strong nonlinear interactions with the electron-acoustic waves [sustained by a multielectron component (degenerate) plasma surrounding a compact astrophysical object] can scatter to lower frequencies so that the radiation observed faraway will be spectrally shifted downward. It is also shown that, under definite conditions, the electromagnetic waves could settle into stationary solitonic states. It is expected that the effects of such structures may persist as detectable signatures in forms of modulated micropulses in the radiation observed far away from the accreting compact object. Both these effects will advance our abilities to interpret the radiation coming out of the compact objects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Optical Imaging of Laser-Driven Fast Electron Weibel-like Filamentation in Overcritical Density Plasma

We report on the measurement of filamented transport of laser-generated fast electron beams in near-critical density plasma. A relativistic intensity long-wave-infrared laser irradiated a hydrodynamically shaped helium gas flow at an electron density n e ≃ 10 25 m − 3 , generating a large flux of fast electrons that propagated beyond the critical surface. The beam-to-background electron density ratio was sufficiently high to drive growth of Weibel-like filamentation, which was measured by optical probing to extend up to 800 μ m with radii ∼ 10 μ m . Particle-in-cell simulations reproduce the main features of the filamentation generation, suggesting that collisionless processes are dominant in these interactions. Expansion of the filaments after formation infers a fast electron heated plasma temperature ∼ 400 eV in the overcritical density plasma. Published by the American Physical Society 2025

43 PARTICLE ACCELERATORS↗

DOE OSTI2020

MHD stability and disruptions in the SPARC tokamak

SPARC is being designed to operate with a normalized beta of $β_N$ = 1.0, a normalized density of $n_G$ = 0.37 and a safety factor of $q_{95}$ ≈ 3.4, providing a comfortable margin to their respective disruption limits. Further, a low beta poloidal $β_p$ = 0.19 at the safety factor $\textit{q}$ = 2 surface reduces the drive for neoclassical tearing modes, which together with a frozen-in classically stable current profile might allow access to a robustly tearing-free operating space. Although the inherent stability is expected to reduce the frequency of disruptions, the disruption loading is comparable to and in some cases higher than that of ITER. The machine is being designed to withstand the predicted unmitigated axisymmetric halo current forces up to 50 MN and similarly large loads from eddy currents forced to flow poloidally in the vacuum vessel. Runaway electron (RE) simulations using GO+CODE show high flattop-to-RE current conversions in the absence of seed losses, although NIMROD modelling predicts losses of ~80 %; self-consistent modelling is ongoing. A passive RE mitigation coil designed to drive stochastic RE losses is being considered and COMSOL modelling predicts peak normalized fields at the plasma of order 10 –2 that rises linearly with a change in the plasma current. Massive material injection is planned to reduce the disruption loading. A data-driven approach to predict an oncoming disruption and trigger mitigation is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Two Mechanisms Limiting the Emitted Electron Current from a Cathode to an Anode

It is known that the current of emitted electrons flowing through a plasma can saturate upon formation of a potential well adjacent to the cathode (the “space charge effect”). Here, we demonstrate another saturation mechanism that will often set a more restrictive limit on the global current. When “backflow saturation” occurs, the cathode sheath weakens to allow emitted electrons that already entered the plasma to backflow to the cathode. This effect could not be captured by studies modelling the cathode sheath by itself because its origin is coupled to processes in the interior plasma and anode sheath. By modeling a full plasma diode, we show that depending on conditions the global current can be limited in four ways; by backflow alone, by space charge alone, by both mechanisms in a stable cooperative form, or by both in a competing oscillatory form. Published by the American Physical Society 2025

Laboratory plasma↗

DOE OSTI2022

Optimum Injector Parameters for Thermoacoustic Stability in a Multi-Nozzle Can Combustion System

Abstract High-frequency transverse instabilities are an important concern in can combustor configurations. In these configurations which are typically operated with multiple injectors around a central injector, each injector is subjected to different parts of the acoustic mode shape and thus respond differently for the same instability mode. Recent work by the author has modeled the response of premixed flames to excitation by natural high-frequency transverse modes in a can combustor both in the center and outer nozzles. The stability of these acoustically non-compact flames was assessed using the Rayleigh criterion (Rayleigh Integral denoted as RI) and not the overall unsteady heat release as is the case for compact flames. Several key control parameters were studied, namely — flame angle, swirling strength, nozzle location. For non-axisymmetric modes such as the commonly occurring 1-T mode, both radial and azimuthal offsets of the nozzle location affected stability. The framework was applied to an optimization study to identify the optimal combination of parameters that minimizes RI for the different nozzles in the multi-nozzle system. In this study, a N-around-1 configuration was studied, and the results indicated that the different nozzles needed to be operated at different flame angles and swirl numbers to result in an overall minimum RI. However, the specific response of the different injectors was not considered. The helical mode distribution at each injector varies as we azimuthally go around the combustor’s injector distribution and thus the most amplified mode and the resulting flame response would be different. To minimize RI, it is important to determine the injector configurations that result in a hydrodynamic profile that minimizes the individual RI for each nozzle. The resulting relationship between the injector’s flow and local hydrodynamics can then be used in a hydrodynamics study of an individual injector so that the most optimal injector is chosen depending on its location in the combustor dump plane.

Acharya, Vishal↗

DOE OSTI2023

Measuring scale-dependent shape anisotropy by coarse-graining: Application to inhomogeneous Rayleigh-Taylor turbulence

Here, we generalize the “filtering spectrum” to probe scales along different directions by spatial coarse-graining. This multidimensional filtering spectrum quantifies the spectral content of flows that are not necessarily homogeneous. From multidimensional spectral information, we propose a simple metric for shape anisotropy at various scales. The method is applied to simulations of two-dimensional (2D) and 3D Rayleigh-Taylor (RT) turbulence, which is inhomogeneous and anisotropic. We show that 3D RT has clear shape anisotropy at large scales with approximately 4:3 vertical to horizontal aspect ratio but tends toward isotropy at small scales as expected. In sharp contrast, we find that RT in 2D simulations, which are still the main modeling framework for many applications, is isotropic at large scales and its shape anisotropy increases at smaller scales where structures tend to be horizontally elongated. While this may be surprising, it is consistent with recent results that large-scale isotropy in 2D RT is due to the generation of a large-scale overturning circulation via an upscale cascade, while small-scale anisotropy is due to the stable stratification resultant from such overturning and the inefficient mixing in 2D.

42 ENGINEERING↗

DOE OSTI2022

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗