SEARCH · Engineering Papers
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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
This presentation will be presented at FY 2020 SDRD Annual (Webex) Meeting scheduled for September 23–24, 2020 High-Fidelity Dynamic Neutron Imaging for Subcritical Experiments and Other Applications LO-005-19, Year 2
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Electron and X-ray beams originating from compact laser-wakefield accelerators have very small source sizes that are typically on the micrometre scale. Therefore, the beam divergences are relatively high, which makes it difficult to preserve their high quality during transport to applications. To improve on this, tremendous efforts have been invested in controlling the divergence of the electron beams, but no mechanism for generating collimated X-ray beams has yet been demonstrated experimentally. Here we propose and realize a scheme where electron bunches undergoing focusing in a dense, passive plasma lens can emit X-ray pulses with divergences approaching the incoherent limit. Furthermore, compared with conventional betatron emission, the divergence of this so-called plasma lens radiation is reduced by more than an order of magnitude in solid angle, while maintaining a similar number of emitted photons per electron. This X-ray source offers the possibility of producing brilliant and collimated few-femtosecond X-ray pulses for ultra-fast science, in particular for studies based on X-ray diffraction and absorption spectroscopy.
High energy density physics is the field of physics dedicated to the study of matter and plasmas in extreme conditions of temperature, densities and pressures. It encompasses multiple disciplines such as material science, planetary science, laboratory and astrophysical plasma science. For the latter, high energy density states can be accompanied by extreme radiation environments and super-strong magnetic fields. The creation of high energy density states in the laboratory consists in concentrating/depositing large amounts of energy in a reduced mass, typically solid material sample or dense plasma, over a time shorter than the typical timescales of heat conduction and hydrodynamic expansion. Laser-generated, high current–density ion beams constitute an important tool for the creation of high energy density states in the laboratory. Focusing plasma devices, such as cone-targets are necessary in order to focus and direct these intense beams towards the heating sample or dense plasma, while protecting the proton generation foil from the harsh environments typical of an integrated high-power laser experiment. A full understanding of the ion beam dynamics in focusing devices is therefore necessary in order to properly design and interpret the numerous experiments in the field. In this work, we report a detailed investigation of large-scale, kilojoule-class laser-generated ion beam dynamics in focusing devices and we demonstrate that high-brilliance ion beams compress magnetic fields to amplitudes exceeding tens of kilo-Tesla, which in turn play a dominant role in the focusing process, resulting either in a worsening or enhancement of focusing capabilities depending on the target geometry.
Accurate phase diagrams of multicomponent plasmas are required for the modeling of dense stellar plasmas, such as those found in the cores of white dwarf stars and the crusts of neutron stars. Those phase diagrams have been computed using a variety of standard techniques, which suffer from physical and computational limitations. Here we present an efficient and accurate method that overcomes the drawbacks of previously used approaches. In particular, finite-size effects are avoided as each phase is calculated separately; the plasma electrons and volume changes are explicitly taken into account; and arbitrary analytic fits to simulation data as well as particle insertions are avoided. Furthermore, no simulations at “uninteresting” state conditions, i.e., away from the phase coexistence curves, are required, which improves the efficiency of the technique. The method consists of an adaptation of the so-called Gibbs-Duhem integration approach to electron-ion plasmas, where the coexistence curve is determined by direct numerical integration of its underlying Clapeyron equation. The thermodynamics properties of the coexisting phases are evaluated separately using Monte Carlo simulations in the isobaric semigrand canonical ensemble (NPT Δ μ ). We describe this Monte Carlo-based Clapeyron integration method, including its basic physical and numerical principles, our extension to electron-ion plasmas, and our numerical implementation. We illustrate its applicability and benefits with the calculation of the melting curve of dense carbon-oxygen plasmas under conditions relevant for the cores of white dwarf stars and provide analytic fits to implement this new melting curve in white dwarf models. While this work focuses on the liquid-solid phase boundary of dense two-component plasmas, a wider range of physical systems and phase boundaries are within the scope of the Clapeyron integration method, which had until now only been applied to simple model systems of neutral particles.
Kr L-shell spectroscopy modeling results are discussed in this paper, focusing on the n = 4 to n = 2 line transitions of Be- and Li-like Kr ions. Collisional radiative atomic kinetic and Stark-broadened spectral line shape calculations show electron temperature Te and density ne sensitivity in the spectrum. The combination of the Te dependence due to the relative intensity of Be-like to Li-like line emissions in the range from 1.5 to 3 keV and the ne sensitivity from the Stark broadening effect on the line shapes in the range from 5×1023 to 2×1024/ cc results in a spectrum that can be employed to diagnose Te and ne. Two different collisional radiative atomic kinetic models i.e., Prismspect [J. J. MacFarlane, et al., Int. Fusion Sci. Appl. Conf. Proc. 457 (2003)] and ABAKO [Florido, et al., PRE, 80, 056402 (2009)] produce similar results in level populations and spectra. In x-ray spectroscopy of implosion cores, this Kr L-shell spectrum may prove useful in an intermediate Te range in which Ar is too ionized for its K-shell to be of diagnostic value and Kr is not ionized enough for its K-shell emission to be useful.
Strongly localized, propagating plasma density structures that are capable of crossing magnetic field lines are known as “blobs.” Here we demonstrate a novel mechanism for the formation and propagation of an ion gyroradius-scale blob-cavity structure at the interface between a super-Alfvénic laser-produced plasma (LPP) and an ambient magnetized plasma. The LPP self-focuses along the edge of the diamagnetic cavity which results in a dense, jetlike structure as compared to ballistic motion. This collimated flow couples momentum to the ambient plasma through a collisionless process known as Larmor coupling. The Larmor electric fields locally displace the ambient ions forming a blob above the LPP flow. In the region between a gyrating blob and collimated LPP flow, a secondary cavity of expelled magnetic field forms. In conclusion, these findings are supported by particle-in-cell simulations that replicate the blob formation mechanism and provide insight to similar processes in space, astrophysical, and laboratory settings characterized by ion kinetic scales.
Ablation pressure is a key parameter controlling laser drive performance in direct-drive inertial confinement fusion implosions. In this work, we develop a theoretical framework to clarify the role of thermal conduction in determining how ablation pressure scales with laser intensity, wavelength, and target size. We focus on a steady-state, spherically expanding plasma produced when an over-dense pellet is irradiated by a high-intensity laser. Scaling laws for ablation pressure and the mass ablation rate are derived with respect to laser intensity, wavelength, and target size. As a result, the theoretical predictions are found to be in good agreement with one-dimensional hydrodynamics simulations in spherical geometry.
Current and future applications of intense proton sources abound, including radiography, cancer therapy, warm dense matter generation, and inertial confinement fusion. With increasingly efficient acceleration and focusing mechanisms, proton current densities may soon approach and exceed 10 10 A/cm 2 , e.g., via intense laser drivers. Simulations have previously shown that in this current density regime, beam-induced field generation plays a significant role in beam transport through dense plasmas. Here, we present a theoretical model for the generation of resistive magnetic fields by intense proton beam transport through solid density plasmas. The theoretical evolution of the magnetic field profile is calculated using an analytic model for aluminum resistivity, heat capacity, and stopping power, applicable from cold matter to hot plasma. The effects of various beam and material parameters on the field are investigated and explained for both monoenergetic and Maxwellian proton beams. For a proton beam with Maxwellian temperature 5 MeV and total energy 10 J, the model calculates resistive magnetic fields up to 150 T in aluminum. In conclusion, the calculated field profiles from several beam cases are compared with 2D hybrid particle-in-cell simulations, with good agreement found in magnitude and time scale.
The fast development of quantum technologies over the last decades has offered a glimpse to a future where the quantum properties of multi particle systems might be more fully understood. In particular, quantum computing might prove crucial to explore many aspects of high energy physics unaccessible to classical methods. In this talk, we will describe how one can use digital quantum computers to study the evolution of QCD jets in quark gluons plasmas. We construct a quantum circuit to study single particle evolution in a dense QCD medium. Focusing on the jet quenching parameter $\hat{q}$, we present some early numerical results for a small quantum circuit. Future extensions of this strategy are also addressed.
In this paper we discuss design considerations and beam dynamics challenges associated with laser-driven plasma-based accelerators as applied to multi-TeV-scale linear colliders. Plasma accelerators provide ultra-high gradients and ultra-short bunches, offering the potential for compact linacs and reduced power requirements. We show that stable, efficient acceleration with beam quality preservation is possible in the nonlinear bubble regime of laser-plasma accelerators using beam shaping. Ion motion, naturally occuring for dense beams (i.e., low emittance and high energy) severely damps transverse beam instabilities. Coulomb scattering by the background ions is considered and it is shown that the strong focusing in the plasma strongly suppresses scattering-induced emittance growth. Finally, Betatron radiation emission from the transverse motion of the beam in the plasma will result in beam power loss and energy spread growth; however for sub-100 nm emittances, the beam power loss and energy spread growth will be sub-percent for multi-TeV-class plasma linacs.
Molecular dynamics (MD) simulations are a powerful tool for the calculation of transport properties in mixtures. Not only are MD simulations capable of treating multicomponent systems, they are also applicable over a wide range of temperatures and densities. In plasma physics, this is particularly important for applications such as inertial confinement fusion. While many studies have focused on the effect of plasma coupling on transport properties, here we focus on the effects of mixing. We compute the thermodynamic factor, a measure of ideal/non-ideal mixing, for three binary ionic mixtures. Here, we consider mixtures of hydrogen and carbon, hydrogen and argon, and argon and carbon, each at 500 randomly generated state points in the warm dense matter and plasma regimes. The calculated thermodynamic factors indicate different mixing behavior across phase space, which can significantly affect the corresponding mutual diffusion coefficients. As MD simulations are still computationally expensive, we apply modern data science tools to predict the thermodynamic factor over a large phase space. Further, we propose a more accurate approximation to the mutual diffusion coefficient than the commonly applied Darken relation.
We study the modifications to the fragmentation pattern of partons into charged particles in the presence of a hot and dense quark gluon plasma. To this end, we analyze the perturbative renormalization group equations of the track functions, which describe the energy fraction carried by charged hadrons. Focusing on pure Yang-Mills theory, we compute the lowest-order moments of the medium-modified track functions, which are found to be sensitive to the reduced phase space for emissions in the medium and to energy loss. We use the extracted moments to calculate the energy energy correlator (EEC) on tracks in the collinear limit. The EEC on medium-evolved tracks does not differ qualitatively from the EEC on vacuum tracks despite being sensitive to the color decoherence transition and suppressing the distribution due to quenching, as seen in other jet observables. Published by the American Physical Society 2024
This documents is the final technical report for our grant entitled "Path Integral Monte Carlo Simulations of Iron Plasmas" that focused on developing path integral Monte Carlo (PIMC) computer simulations. This techniques will be developed to study plasmas composed of heavier elements including iron and other third row elements. Equations of state (EOS) and transport properties will be derived in the regime of warm dense matter (WDM) and dense plasmas where existing first-principles methods cannot be applied. While standard density functional theory (DFT) has been used to accurately predict the structure of many solids and liquids up to temperatures on the order of 100,000 K, this method is not applicable at much higher temperature because the number of partially occupied electronic orbitals reaches intractably large numbers or the use of finite-temperature free energy functionals in orbital-free DFT introduces an uncontrolled approximation. Here we focus on PIMC methods that become more and more efficient with increasing temperatures and still include all electronic correlation effects. In this approach, electronic excitations increase the efficiency rather than reduce it. While it had commonly been assumed this method could only be applied to elements without core electrons, we showed that PIMC with free-particle nodes works well for first-row elements (PRL 108 (2012) 115502). Most recently, we extended the applicability range of all-electron PIMC to second-row elements by adopting localized nodal surfaces (PRL 115 (2015) 176403). To simulate third-row elements efficiently under WDM conditions, we propose a new method to remove core electrons by introducing pseudo-nodes. We explain our approach step by step and present preliminary results. We focus our method development on getting PIMC simulations of iron to work because of its fundamental importance for WDM and astrophysics. Then we move on to krypton and copper-doped beryllium, a ICF ablator material. We plan to continue working on key second-row material such as Na, Mg, MgO, Al, silica, and silicon-doped plastic ablators. Our collaborators at LLNL, will use our PIMC EOS data both as comparisons to existing semi-empirical, EOS-generating schemes, and as input for continuum radiation hydrodynamics simulations. We will establish an efficient pipeline from PIMC to macroscopic continuum studies of materials response. An emphasis will be placed on benchmarking such methods for plasmas of heavy elements at the very high temperatures (~100 eV) and low densities that are generated when Hohlraum radiation heats the ablator material in indirect drive laser experiments. Results from changes to the EOS will be of immeasurable importance to the designers at the National Ignition Facility (NIF) and at other facilities. Starting with our EOS of Cu-doped Be, our second collaborator at LLE, will perform real-time simulations of laser fusion experiments at the Omega laser and at the NIF to determine how sensitive the compression path depends on the ablator EOS. Since our collaborator also has experience in performing orbital-free DFT calculations, we propose to compare predictions from this method with PIMC results. In joint publications, we plan to analyze the accuracy of different free-energy functionals in order to understand why existing orbital-free DFT calculations do not predict compression peaks along the shock Hugoniot curve that we see with PIMC. The peaks are caused by the ionization of various electron shells. Their accurate characterization is important to compare with experimental results. We will break new ground by developing PIMC techniques that can simulate iron and all other third row elements in the plasma and WDM regimes. We introduce the concept of pseudo-nodes for the efficient treatment of the core-electrons. The EOS and transport properties will be derived and published online in the form of a new WDM database. Our PIMC EOS calculations will benchmark and possibly replace semi-analytical EOS tables like QEOS or SESAME, which will impact the hydrocode simulation community and will affect the design of NIF targets. Our PIMC results will help to improve the accuracy of orbital-free DFT simulations.
We report the first fully kinetic, quantum study of kinetic electrostatic electron nonlinear (KEEN) waves, showing that quantum diffraction systematically erodes the classical trapping mechanism, narrows harmonic locking to the fundamental, and hastens post-drive decay. Electrons are evolved with a second-order Strang-split 1D1V Wigner–Poisson solver that couples conservative semi-Lagrangian WENO advection to an analytic Fourier space update for the non-local Wigner term, while ions remain classical. We focus on collisionless dynamics in a weakly coupled regime, providing a controlled baseline before collisional extensions. Short, frequency-tuned ponderomotive pulses drive KEEN formation in a uniform Maxwellian plasma; as the dimensionless quantum parameter H rises from the classical limit to values relevant to warm-dense matter, doped semiconductors, and 2D electron systems, the drive threshold increases, higher harmonics are damped, trapped electron vortices diffuse, and the subplasma electrostatic energy relaxes to a lower stationary level, as confirmed by continuous wavelet analysis. These microscopic changes carry macroscopic weight. Ignition-scale capsules now compress matter to regimes where the electron de Broglie wavelength rivals the Debye length, making classical kinetic descriptions insufficient. By extending KEEN physics into this quantum domain, our results offer a potential diagnostic of non-equilibrium electron dynamics for next-generation inertial-confinement designs and high-energy-density platforms, indicating that predictive fusion modeling may benefit from the integration of kinetic fidelity with quantum effects.
High repetition rate laser systems enable new strategies for diagnosing plasma behavior with large datasets. Here, we define an ensemble technique that relies on randomized targeting of x-ray tracer micro-stripes. On each shot, a high-intensity laser pulse is focused on a solid target with Ti tracer stripes embedded in an Al foil, randomly targeting a micro-stripe, a portion of a stripe, or a gap between stripes. High-resolution, time-integrated x-ray spectrometers capture line emission from the portion of the micro-stripe that is heated to sufficiently high electron temperatures. Accumulation of many such cases is used to construct ensemble distributions of x-ray line intensities that encompass all relative offsets of the laser focus to the micro-stripe centers. Synthetic intensity distributions are likewise generated using collisional-radiative modeling. Bayesian fitting of modeled to measured intensity distributions establishes the most likely radial temperature profiles, enabling comparison to hydrodynamic models and calling into question the cylindrical symmetry of these micro-stripe-embedded systems. Ensemble techniques have significant potential for high-energy-density plasma diagnostics, especially with the advent of high repetition rate experiments.