Search NASA⌕ Search

Engineering topics

Bhattacharjee, A.

Publications and source records attributed to Bhattacharjee, A..

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An improved ten-moment closure for reconnection and instabilities

The integration of kinetic effects in fluid models is important for global simulations of Earth's magnetosphere. The use of the two-fluid ten-moment model, which includes the pressure tensor for both species, has had some success in simulating Ganymede and Mercury with a simple closure model. Here, we discuss a heat flux closure which accounts for some limitations of the earlier work while remaining computationally tractable. Comparisons with kinetic simulations for magnetic reconnection and lower-hybrid drift instabilities show good agreement with kinetic results and improvements on previous closure models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Necessary and sufficient conditions for quasisymmetry

A necessary and sufficient set of conditions for a quasisymmetric magnetic field in the form of constraint equations is derived from first principles. Without any assumption regarding the magnetohydrodynamic (MHD) equilibrium of the plasma, conditions for quasisymmetry are constructed starting from the single-particle Lagrangian to the leading order. The conditions presented in the paper are less restrictive than the set recently obtained by Burby et al. [“Some mathematics for quasi-symmetry,” arXiv:1912.06468 (2019)], and could facilitate ongoing efforts toward investigating the existence of global quasisymmetric MHD equilibria. Finally, it is also shown that quasisymmetry implies the existence of flux surfaces, regardless of whether the field corresponds to an MHD equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Kinetic simulations of piston-driven collisionless shock formation in magnetized laboratory plasmas

Laboratory laser experiments offer a novel approach to studying magnetized collisionless shocks, and a common method in recent experiments is to drive shocks using a laser-ablated piston plasma. However, current experimental capabilities are still limited to spatio-temporal scales on the order of shock formation, making it challenging to distinguish piston and shock dynamics. We present quasi-1D particle-in-cell simulations of piston-driven, magnetized collisionless shock formation using the code psc , which includes a model of laser-driven plasmas that can be well-matched to experimental conditions. The simulations cover a range of upstream and ablation parameters and yield several robust signatures of shock formation that can provide a reference for experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

An Extended MHD Study of the 16 October 2015 MMS Diffusion Region Crossing

The Magnetospheric Multiscale (MMS) mission has given us unprecedented access to high cadence particle and field data of magnetic reconnection at Earth's magnetopause. MMS first passed very near an X-line on 16 October 2015, the Burch event, and has since observed multiple X-line crossings. Subsequent 3-D particle-in-cell (PIC) modeling efforts of and comparison with the Burch event have revealed a host of novel physical insights concerning magnetic reconnection, turbulence-induced particle mixing, and secondary instabilities. In this study, we employ the Gkeyll simulation framework to study the Burch event with different classes of extended, multifluid magnetohydrodynamics (MHD), including models that incorporate important kinetic effects, such as the electron pressure tensor, with physics-based closure relations designed to capture linear Landau damping. Such fluid modeling approaches are able to capture different levels of kinetic physics in global simulations and are generally less costly than fully kinetic PIC. Herein we focus on the additional physics one can capture with increasing levels of fluid closure refinement via comparison with MMS data and existing PIC simulations. In particular, we find that the ten-moment model well captures the agyrotropic structure of the pressure tensor in the vicinity of the X-line and the magnitude of anisotropic electron heating observed in MMS and PIC simulations. However, the ten-moment model is found to have difficulty resolving the lower hybrid drift instability, which plays a fundamental role in heating and mixing electrons in the current layer.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Direct Observations of Particle Dynamics in Magnetized Collisionless Shock Precursors in Laser-Produced Plasmas

We present the first laboratory observations of time-resolved electron and ion velocity distributions in magnetized collisionless shock precursors. Thomson scattering of a probe laser beam was used to observe the interaction of a laser-driven, supersonic piston plasma expanding through an ambient plasma in an external magnetic field. From the Thomson-scattered spectra we measure time-resolved profiles of electron density, temperature, and ion flow speed, as well as spatially resolved magnetic fields from proton radiography. We observe direct evidence of the coupling between piston and ambient plasmas, including the acceleration of ambient ions driven by magnetic and pressure gradient electric fields, and deformation of the piston ion flow, key steps in the formation of magnetized collisionless shocks. Even before a shock has fully formed, we observe strong density compressions and electron heating associated with the pileup of piston ions. Furthermore, the results demonstrate that laboratory experiments can probe particle velocity distributions relevant to collisionless shocks, and can complement, and in some cases overcome, the limitations of similar measurements undertaken by spacecraft missions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Drift Instabilities in Thin Current Sheets Using a Two–Fluid Model With Pressure Tensor Effects

The integration of kinetic effects in fluid models is important for global simulations of the Earth's magnetosphere. We use a two–fluid 10–moment model, which includes the pressure tensor and has been used to study reconnection, to study the drift kink and lower hybrid drift instabilities. Using a nonlocal linear eigenmode analysis, we find that for the kink mode, the 10–moment model shows good agreement with kinetic calculations with the same closure model used in reconnection simulations, while the electromagnetic and electrostatic lower hybrid instabilities require modeling the effects of the ion resonance using a Landau fluid closure. As a result, comparisons with kinetic simulations and the implications of the results for global magnetospheric simulations are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Magnetohydrodynamical equilibria with current singularities and continuous rotational transform

We revisit the Hahm–Kulsrud–Taylor (HKT) problem, a classic prototype problem for studying resonant magnetic perturbations and 3D magnetohydrodynamical equilibria. We employ the boundary-layer techniques developed by Rosenbluth, Dagazian, and Rutherford (RDR) for the internal m = 1 kink instability, while addressing the subtle difference in the matching procedure for the HKT problem. Pedagogically, the essence of RDR's approach becomes more transparent in the reduced slab geometry of the HKT problem. We then compare the boundary-layer solution, which yields a current singularity at the resonant surface, to the numerical solution obtained using a flux-preserving Grad–Shafranov solver. The remarkable agreement between the solutions demonstrates the validity and universality of RDR's approach. Additionally, we show that RDR's approach consistently preserves the rotational transform, which hence stays continuous, contrary to a recent claim that RDR's solution contains a discontinuity in the rotational transform.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Nonlinear saturation of kinetic ballooning modes by zonal fields in toroidal plasmas

Kinetic ballooning modes (KBMs) are widely believed to play a critical role in disruptive dynamics as well as turbulent transport in magnetic fusion and space plasmas. While the nonlinear evolution of the ballooning modes has been proposed as a mechanism for “detonation” in various scenarios such as the edge localized modes in tokamaks, the role of the kinetic effects in such nonlinear dynamics remains largely unexplored. In this study, global gyrokinetic simulation results of KBM nonlinear behavior are presented. Instead of the finite-time singularity predicted by ideal magnetohydrodynamic theory, the kinetic instability is shown to develop into an intermediate nonlinear regime of exponential growth, followed by a nonlinear saturation regulated by spontaneously generated zonal fields. In the intermediate nonlinear regime, rapid growth of localized current sheets, which can induce magnetic reconnection, is observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

A Model of Solar Equilibrium: The Hydrodynamic Limit

Helioseismology has revealed the internal density and rotation profiles of the Sun. Yet, knowledge of its magnetic fields and meridional circulation is confined much closer to the surface, and latitudinal entropy gradients are below detectable limits. While numerical simulations can offer insight into the interior dynamics and help identify which ingredients are necessary to reproduce particular observations, some features of the Sun can be understood analytically from an equilibrium perspective. Examples of such features include: the 1D density profile arising from steady-state energy transport from the core to the surface, and the tilting of isorotation contours in the convection zone (CZ) due to baroclinic forcing. To help answer the question of which features can be explained by equilibrium, we propose analyzing stationary axisymmetric ideal magnetohydrodynamic flows in the solar regime. By prescribing an appropriate entropy at the surface, we recover a rotation profile that reasonably matches observations in the bulk of the CZ. Additionally, by including the effects of poloidal flow, we reproduce a feature that is reminiscent of the near surface shear layer. Furthermore, no tachocline-like feature is seen in hydrodynamic equilibrium, suggesting the importance of either dynamics or magnetic fields in its description.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Relativistic-electron-driven magnetic reconnection in the laboratory

Magnetic reconnection is a fundamental process occurring in many plasma systems. Magnetic field lines break and reconfigure in to a lower energy state, converting released magnetic field energy into plasma kinetic energy. Around some of the universe’s most energetic objects, such as gamma ray burst or active galactic nuclei, where the magnetic field energy exceeds the plasma rest mass energy, the most extreme magnetic reconnection in the relativistic regime is theorized. The pre- sented experiments and three-dimensional particle-in-cell modeling recreate in the laboratory the scaled plasma conditions necessary to access the relativistic electron regime and therefore approach conditions around these distant, inaccessible objects. High-power, ultrashort laser pulses focused to high-intensity (I > 2.5 × 1018 Wcm-2) on solid targets produces relativistic temperature electrons within the focal volume. The hot electrons are largely confined to the target surface and form a ra- dial surface current that generates a huge, expanding azimuthal magnetic field. Focusing two laser pulses in close proximity on the target surface leads to oppositely directed magnetic fields being driven together. The fast electron motion due to the magnetic reconnection is inferred using an ex- perimental x-ray imaging technique. The x-ray images enable the measurement of the reconnection layer dimensions and temporal duration. The reconnection rates implied from the aspect ratio of the reconnection layer, δ/L ≈ 0.3, was found to be consistent over a range of experimental pulse durations (40 fs–20 ps) and agreed with the modeling. Further experimental evidence for magnetic reconnection is the formation of a nonthermal electron population shown by the modeling to be accelerated in the reconnection layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Kinetic simulation of magnetic field generation and collisionless shock formation in expanding laboratory plasmas

Recent laboratory experiments with laser-produced plasmas have observed and studied a number of fundamental physical processes relevant to magnetized astrophysical plasmas, including magnetic reconnection, collisionless shocks, and magnetic field generation by Weibel instability, opening up new experimental platforms for laboratory astrophysics. We develop a fully kinetic simulation model for first-principles simulation of these systems including the dynamics of magnetic fields—magnetic field generation by the Biermann battery effect or Weibel instability; advection by the ion flow, Hall effect, and Nernst effect; and destruction of the field by dissipative mechanisms. Key dimensionless parameters describing the system are derived for scaling between kinetic simulation, recent experiments, and astrophysical plasmas. First, simulations are presented which model Biermann battery magnetic field generation in plasmas expanding from a thin target. Ablation of two neighboring plumes leads to the formation of a current sheet as the opposing Biermann-generated fields collide, modeling recent laser-driven magnetic reconnection experiments. Second, we simulate recent experiments on collisionless magnetized shock generation, by expanding a piston plasma into a pre-magnetized ambient plasma. For parameters considered, the Biermann effect generates additional magnetic fields in the curved shock front and thereby increases shock particle reflection. In conclusion, both cases show the importance of kinetic processes in the interaction of plasmas with magnetic fields and open opportunities to benchmark these important processes through comparison of theory and experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Regimes of magnetic reconnection in colliding laser-produced magnetized plasma bubbles

Here, we conduct a multiparametric study of driven magnetic reconnection relevant to recent experiments on colliding magnetized laser-produced plasmas using particle-in-cell simulations. Varying the background plasma density, plasma resistivity, and plasma bubble geometry, the 2D simulations demonstrate a rich variety of reconnection behaviors and show the coupling between magnetic reconnection and the global hydrodynamical evolution of the system. We consider both the collision between two radially expanding bubbles where reconnection is seeded by the pre-existing X-point and the collision between two flows in a quasi-1D geometry with initially anti-parallel fields where reconnection must be initiated by the tearing instability. At a baseline case of low-collisionality and low background density, the current sheet is strongly compressed to below scale of the ion-skin-depth scale, and rapid, multi-plasmoid reconnection results. Increasing the plasma resistivity, we observe a collisional slow-down of reconnection and stabilization of plasmoid instability for Lundquist numbers less than approximately S ~10 3 . Second, increasing the background plasma density modifies the compressibility of the plasma and can also slow down or even prevent reconnection, even in completely collisionless regimes, by preventing the current sheet from thinning down to the scale of the ion-skin depth. These results have implications for understanding recent and future experiments, and signatures for these processes for proton-radiography diagnostics of these experiments are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Biermann-Battery-Mediated Magnetic Reconnection in 3D Colliding Plasmas

Recent experiments have demonstrated magnetic reconnection between colliding plasma plumes, where the reconnecting magnetic fields were self-generated in the plasma by the Biermann-battery effect. Using fully kinetic 3D simulations, we show the full evolution of the magnetic fields and plasma in these experiments, including self-consistent magnetic field generation about the expanding plume. The collision of the two plasmas drives the formation of a current sheet, where reconnection occurs in a strongly time- and space-dependent manner, demonstrating a new 3D reconnection mechanism. Specifically, we observe a fast, vertically localized Biermann-mediated reconnection, an inherently 3D process where the temperature profile in the current sheet coupled with the out-of-plane ablation density profile conspires to break inflowing field lines, reconnecting the field downstream. Fast reconnection is sustained by both the Biermann effect and the traceless electron pressure tensor, where the development of plasmoids appears to modulate the contribution of the latter. In conclusion, we present a simple and general formulation to consider the relevance of Biermann-mediated reconnection in general astrophysical scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Using the maximum entropy distribution to describe electrons in reconnecting current sheets

Particle distributions in weakly collisional environments such as the magnetosphere have been observed to show deviations from the Maxwellian distribution. These can often be reproduced in kinetic simulations, but fluid models, which are used in global simulations of the magnetosphere, do not necessarily capture any of this. We apply the maximum entropy fluid closure of Levermore, which leads to well posed moment equations, to reconstruct particle distributions from a kinetic simulation in a reconnection region. Our results show that without information other than the moments, the model can reproduce the general structure of the distributions but not all of the finer details. In conclusion, the advantages of the closure over the traditional Grad closure are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Cross-verification of the global gyrokinetic codes GENE and XGC

A detailed cross-verification between two global gyrokinetic codes, the core continuum code GENE and the edge particle-in-cell code XGC, for the linear and nonlinear simulations of ion-temperature-gradient modes is carried out. With the recent developments in the edge gyrokinetics, it may be feasible someday to describe the whole tokamak plasma on turbulence timescales using a coupled gyrokinetic simulation model. Before pursuing this, the core code (GENE) and the edge code (XGC) must be carefully benchmarked with each other. The present verification provides a solid basis for future code coupling research. Also included in the benchmarking is the global particle-in-cell code ORB5, to raise the confidence in the quality of the obtained results. An excellent agreement between all three codes is obtained. Lastly, in order to facilitate a benchmark framework for other codes, we make a specific effort to provide all the relevant input parameters and precise details for each code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Algebraic motion of vertically displacing plasmas

In this paper, the vertical motion of a tokamak plasma is analytically modelled during its non-linear phase by a free-moving current-carrying rod inductively coupled to a set of fixed conducting wires or a cylindrical conducting shell. The solutions capture the leading term in a Taylor expansion of the Green's function for the interaction between the plasma column and the surrounding vacuum vessel. The plasma shape and profiles are assumed not to vary during the vertical drifting phase such that the plasma column behaves as a rigid body. In the limit of perfectly conducting structures, the plasma is prevented to come in contact with the wall due to steep effective potential barriers created by the induced Eddy currents. Resistivity in the wall allows the equilibrium point to drift towards the vessel on the slow timescale of flux penetration. The initial exponential motion of the plasma, understood as a resistive vertical instability, is succeeded by a non-linear “sinking” behaviour shown to be algebraic and decelerating. Finally, the acceleration of the plasma column often observed in experiments is thus concluded to originate from an early sharing of toroidal current between the core, the halo plasma, and the wall or from the thermal quench dynamics precipitating loss of plasma current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Continuum Reverberation Mapping of the Accretion Disks in Two Seyfert 1 Galaxies

We present optical continuum lags for two Seyfert 1 galaxies, MCG+08-11-011 and NGC 2617, using monitoring data from a reverberation mapping campaign carried out in 2014. Our light curves span the ugriz filters over four months, with median cadences of 1.0 and 0.6 days for MCG+08-11-011 and NGC 2617, respectively, combined with roughly daily X-ray and near-UV data from Swift for NGC 2617. We find lags consistent with geometrically thin accretion-disk models that predict a lag-wavelength relation of τ ∝ λ 4/3 . However, the observed lags are larger than predictions based on standard thin-disk theory by factors of 3.3 for MCG+08-11-011 and 2.3 for NGC 2617. Furthermore, these differences can be explained if the mass accretion rates are larger than inferred from the optical luminosity by a factor of 4.3 in MCG+08-11-011 and a factor of 1.3 in NGC 2617, although uncertainty in the SMBH masses determines the significance of this result. While the X-ray variability in NGC 2617 precedes the UV/optical variability, the long (2.6 day) lag is problematic for coronal reprocessing models.

79 ASTRONOMY AND ASTROPHYSICS↗