Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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.

At least 1,135 records · Page 63

Reversible Exsolution of Nanometric Fe 2 O 3 Particles in BaFe 2-x (PO 4 ) 2 (0 ≤ x ≤ 2/3):The Logic of Vacancy Ordering in Novel Metal-Depleted Two-Dimensional Lattices

We show here that the exsolution of Fe 2+ ions out of two-dimensional (2D) honeycomb layers of BaFe 2 (PO 4 ) 2 into iron-deficient BaFe 2–x (PO 4 ) 2 phases and nanometric α-Fe 2 O 3 (typically 50 nm diameter at the grain surface) is efficient and reversible until x = 2/3 in mild oxidizing/reducing conditions. It corresponds to the renewable conversion of 12 wt % of the initial mass into iron oxide. After analyzing single crystal X-ray diffraction data of intermediate members x = 2/7, x = 1/3, x = 1/2 and the ultimate Fe-depleted x = 2/3 term, we then observed a systematic full ordering between Fe ions and vacancies (V Fe ) that denote unprecedented easy in-plane metal diffusion driven by the Fe 2+ /Fe 3+ redox. Besides the discovery of a diversity of original depleted triangular ∞ {Fe 2/3+ 2–x O 6 } topologies, we propose a unified model correlating the x Fe-removal and the experimental Fe/V Fe ordering into periodic one-dimensional motifs paving the layers, gaining insights into predictive crystahemistry of complex low dimensional oxides. When we increased the x values it led to a progressive change of the materials from 2D ferromagnets (Fe2+) to 2D ferrimagnets (Fe2/3+) to antiferromagnets for x = 2/3 (Fe3+).

Alcover, Ignacio Blazquez↗

The effect of relative humidity and temperature on the response of stomatal conductance to vapor pressure deficit in tropical trees

Understanding how leaf gas exchange responds to changes in vapor pressure deficit (VPD) is key to predicting tropical forest resilience to climate change. Stomata regulate leaf water and CO2 diffusion, and respond to changes in temperature and relative humidity (RH), two drivers of VPD. At high temperatures, the cuticular pathway may also become significant and participate in the overall leaf conductance. Here, we measured gas exchange under light and dark conditions to investigate the stomatal and cuticular responses to temperature and RH on detached branches in five tropical tree species. Leaf conductance in the dark, when stomata are essentially closed, was not markedly impacted by temperature and RH, suggesting a minimal response of the cuticular pathway to these conditions. We compared six steady-state conductance models incorporating different effects of photosynthesis and evaporative demand on stomatal control. All models performed well (residual standard deviation, σ, < 0.025 mol m-2 s-1), but the best-fitting model (σ = 0.017 mol m-2 s-1) used a nonlinear relationship between photosynthesis and stomatal conductance and incorporated RH rather than vapor pressure difference. All models overestimated the steady-state conductance at high RH. Furthermore, leaf conductance immediately increased after a decrease in RH (wrong-way response), but not after an increase in leaf temperature. This suggests that the mechanisms underlying stomatal response to VPD need further investigation. The non-linear coupling between photosynthetic rate and stomatal conductance indicates a sharper physiological response than previously acknowledged.

Lamour, Julien↗

Microstructural and infiltration properties of woven preforms during chemical vapor infiltration

Interface-resolved direct numerical simulations (DNSs) of chemical vapor infiltration (CVI) have been performed over a range of furnace-operating conditions (Thiele moduli) and for practical woven preform geometries. A level-set method is used to resolve the geometry of the initial preform at tow scale. The interface between the vapor and solid phase is then evolved in time through the entire CVI densification cycle, fully resolving the time-varying topology between the two phases. In contrast to previous level-set methods for CVI simulation, the physical reaction and diffusion processes govern the level-set movement in the current approach. The surface deposition kinetics is described by the usual one-step model. In this paper, the DNS data are used to study the evolving porosity, surface-to-volume ratio, and flow infiltration properties (permeability and effective diffusivities). Additionally, comparisons are made to popularly-assumed structure functions and the standard, Kozeny–Carmen porous media model commonly employed in modeled CFD simulations of CVI. The virtual DNS experiments reveal a Thiele modulus and preform geometry (fabric layup) dependence which the existing microstructural and infiltration models are not able to describe throughout the entire densification process. The DNS-based, woven geometry-specific correlations can be applied directly to mean-field, furnace-scale CFD simulations.

36 MATERIALS SCIENCE↗

The Application of an Evolutionary Programming Process to a Simulation of the ETEX Large-Scale Airborne Dispersion Experiment

Abstract Airborne tracer simulations are typically performed using a dispersion model driven by a high-resolution meteorological model. Besides solving the dynamic equations of momentum, heat, and moisture on the resolved model grid, mesoscale models must account for subgrid-scale fluxes and other unresolved processes. These are estimated through parameterization schemes of eddy diffusion, convection, and surface interactions, and they make use of prescribed parameters set by the user. Such “free” model parameters are often poorly constrained, and a range of plausible values exists for each. Evolutionary programming (EP) is a process to improve the selection of the parameters. A population of simulations is first run with a different set of parameter values for each member, and the member judged most accurate is selected as the “parent” of a new “generation.” After a number of iterations, the simulations should approach a configuration that is best adapted to the atmospheric conditions. We apply the EP process to simulate the first release of the 1994 European Tracer Experiment (ETEX) project, which comprised two experiments in which a tracer was released in western France and sampled by an observing network. The EP process is used to improve a simulation of the RAMS mesoscale weather model, with weather data collected during ETEX being used to “score” the individual members according to how well each simulation matches the observations. The meteorological simulations from before and after application of the EP process are each used to force a dispersion model to create a simulation of the ETEX release, and substantial improvement is observed when these are validated against sampled tracer concentrations.

54 ENVIRONMENTAL SCIENCES↗

Connecting the dots: Semi-analytical and random walk numerical solutions of the diffusion–reaction equation with stochastic initial conditions

We study a system with bimolecular irreversible kinetic reaction A+B→∅ where the underlying transport of reactants is governed by diffusion, and the local reaction term is given by the law of mass action. We consider the case where the initial concentrations are given in terms of an average and a white noise perturbation. Our goal is to solve the diffusion–reaction equation which governs the system, and we tackle it with both analytical and numerical approaches. To obtain an analytical solution, we develop the equations of moments and solve them approximately. To obtain a numerical solution, we develop a grid-less Monte Carlo particle tracking approach, where diffusion is modeled by a random walk of the particles, and reaction is modeled by annihilation of particles. The probability of annihilation is derived analytically from the particles' co-location probability. We rigorously derive the relationship between the initial number of particles in the system and the amplitude of white noise represented by that number. This enables us to compare the particle simulations and the approximate analytical solution and offer an explanation of the late time discrepancies. - Graphical abstract:.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Phase Field Modeling of Pressure Induced Densification in Solid Electrolytes

Adoption of dense and homogeneous solid electrolytes can possibly mitigate the propagation of lithium dendrites and enable lithium metal anodes. Application of external pressure helps to minimize the sintering temperature in oxide ceramics and can potentially densify softer sulfide electrolytes even under room temperature conditions. Here, a previously developed phase field-based computational scheme for predicting the high-temperature sintering-induced densification of oxide ceramic solid electrolytes is extended in the present context to capture the influence of external pressure for densifying solid electrolytes. Two different bulk deformation mechanisms, namely, "reorganization" and "creep deformation," are dominant under external pressure, which is different from the surface and grain-boundary diffusion-induced densification of solid electrolytes that occurs during high temperature sintering. External pressure also increases the points of contact between the particles, which further enhances the propensity of diffusion-induced sintering process. Results obtained from simulations indicate that densification under external pressure is independent of the solid electrolyte particle morphology. Finally, a phase map is generated between applied pressure and temperature for achieving complete densification of oxide ceramics, which can possibly guide the synthesis of thin and dense solid electrolyte separators.

25 ENERGY STORAGE↗

Current challenges in the physics of white dwarf stars

White dwarfs are a class of stars with unique physical properties. They present many challenging problems whose solution requires advanced theories of dense matter, state-of-the-art experimental techniques, and extensive computing efforts. New ground- and space-based observatories will soon provide an increasingly detailed view of white dwarf stars and reveal new phenomena that will challenge our models. This review is an introduction to the nature of white dwarfs, the physical processes that determine their structure and evolution, and the physical conditions they span. Here we discuss a wide variety of currently unsolved or partially resolved problems in their constitutive physics that are broadly related to equations of state, transport processes and opacities.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling the molecular gas in NGC 6240

We present the first observations of H{sup 13}CN (1−0), H{sup 13}CO{sup +}(1−0), and SiO (2−1)in NGC 6240, obtained with the IRAM Plateau de Bure Interferometer. Combining a Markov Chain Monte Carlo code with Large Velocity Gradient (LVG) modeling, and with additional data from the literature, we simultaneously fit three gas phases and six molecular species to constrain the physical condition of the molecular gas, including mass−luminosity conversion factors. We find ∼10{sup 10}M{sub ⊙} of dense molecular gas in cold, dense clouds (T{sub k}∼10 K, n{sub H{sub 2}}∼10{sup 6} cm{sup −3}) with a volume filling factor <0.002, embedded in a shock heated molecular medium (T{sub k}∼2000 K, n{sub H{sub 2}}∼10{sup 3.6} cm{sup −3}), both surrounded by an extended diffuse phase (T{sub k}∼200 K, n{sub H{sub 2}}∼10{sup 2.5} cm{sup −3}). We derive a global α{sub CO}=1.5{sub 1.1}{sup 7.1} with gas masses log{sub 10}(M/[M{sub ⊙}])=10.1{sub 10.0}{sup 10.8}, dominated by the dense gas. We also find α{sub HCN}=32{sub 13}{sup 89}, which traces the cold, dense gas. The [{sup 12}C]/[{sup 13}C] ratio is only slightly elevated (98{sub 65}{sup 230}), contrary to the very high [CO]/[{sup 13}CO] ratio (300–500) reported in the literature. However, we find very high [HCN]/[H{sup 13}CN] and [HCO{sup +}]/[H{sup 13}CO{sup +}] abundance ratios (300{sub 200}{sup 500}) which we attribute to isotope fractionation in the cold, dense clouds.

79 ASTRONOMY AND ASTROPHYSICS↗

Theoretical and Experimental Study of Current from Non-Disintegrable Suspended Particles at a Rotating Disk Electrode

Understanding the current response at an electrode from suspended solid particles in an electrolyte is crucial for developing materials to be used in semi-solid electrodes for energy storage applications. Here, an analytical model is proposed to predict and understand the current response from non-disintegrable solid particles at a rotating disk electrode. The current is shown to be limited by a combination of ion diffusion within the solid particle and the mean residence time of the particle at the rotating disk electrode. This results in a relationship between current and angular frequency of I ∝ ω 3 / 4 , instead of the classical I ∝ ω 1 / 2 predicted by Levich theory. Specifically, the current response of Li 4 Ti 5 O 12 (LTO) microparticles suspended in a non-aqueous electrolyte of lithium hexafluorophosphate (LiPF 6 ) in ethylene carbonate: diethyl carbonate (EC:DEC) was determined experimentally and compared favorably with predictions from the proposed analytical model using fitting parameters consistent with the experimental conditions.

25 ENERGY STORAGE↗

A Simplified Approach Based on Cellular Automata for Describing Direct Reduced Iron Production in Different Reducing Conditions

A quick computation approach based on cellular automata is developed and implemented to describe the reduction of iron ore pellets by a mixture of reducing agents featured by different H 2 /CO ratios. The evolution of oxygen concentration inside the pellet is followed from the beginning to the end of contact between reducing agent and pellet. The variation of thermal state of pellets and gas mixture is computed based on their initial temperature, considering the heat involved and the convective heat exchange between pellet and gas mixture. The use of cellular automata and finite‐difference method to solve the diffusion equation point out the absence of any diffusion coefficient value, allowing to make the model fit the experimental trial, because the problem is that it is not ruled just by diffusion but also by the concentration variation of reducing agent inside the pellet due to porosity increasing during reduction. The updating of the reducing agents concentration implies a sharp decrease of oxygen concentration that the cellular automata model considers. The developed model is able to provide the in‐line control of reduction process and could be used to adjust the chemical concentration and temperature of injected reducing agents.

Metallurgy & Metallurgical Engineering↗

Advanced Technologies for Monitoring CO2 Saturation and Pore Pressure in Geologic Formations: Linking the Chemical and Physical Effects to Elastic and Transport Properties

Ultrasonic P- and S-wave velocities were measured over a range of confining pressures while injecting CO2 and brine into the samples. Pore fluid pressure was also varied and monitored together with porosity during injection. Effective medium models were developed to understand the mechanisms and impact of observed changes and to provide the means for implementation of the interpretation methodologies in the field. Ultrasonic P- and S-wave velocities in carbonate rocks show as much as 20-50% decrease after injection of the reactive CO2-brine mixture; the changes were caused by permanent changes to the rock elastic frame associated with dissolution of mineral. Velocity decreases were observed under both dry and fluid-saturated conditions, and the amount of change was correlated with the initial pore fabrics. Scanning Electron Microscope images of carbonate rock microstructures were taken before and after injection of CO2-rich water. The images reveal enlargement of the pores, dissolution of micrite (micron-scale calcite crystals), and pitting of grain surfaces caused by the fluid- solid chemical reactivity. The magnitude of the changes correlates with the rock microtexture – tight, high surface area samples showed the largest changes in permeability and smallest changes in porosity and elastic stiffness compared to those in rocks with looser texture and larger intergranular pore space. Changes to the pore space also occurred from flow of fine particles with the injected fluid. Carbonates with grain-coating materials, such as residual oil, experienced very little permanent change during injection. In the tight micrite/spar cement component, dissolution is controlled by diffusion: the mass transfer of products and reactants is thus slow and the fluid is expected to be close to thermodynamical equilibrium with the calcite, leading to very little dissolution, or even precipitation. In the microporous rounded micrite and macropores, dissolution is controlled by advection: because of an efficient mass transfer of reactants and products, the fluid remains acidic, far from thermodynamical equilibrium and the dissolution of calcite is important. These conclusions are consistent with the lab observations. Sandstones from the Tuscaloosa formation in Mississippi were also subjected to injection under representative in situ stress and pore pressure conditions. Again, both P- and S-wave velocities decreased with injection. Time-lapse SEM images indicated permanent changes induced in the sandstone microstructure by chamosite dissolution upon injection of CO2-rich brine. After injection, the sandstone showed an overall cleaner microstructure. Two main changes are involved: (a) clay dissolution between grains and at the grain contact and (b) rearrangement of grains due to compaction under pressure Theoretical and empirical models were developed to quantify the elastic changes associated with injection. Permanent changes to the rock frame resulted in seismic velocity-porosity trends that mimic natural diagenetic changes. Hence, when laboratory measurments are not available for a candidate site, these trends can be estimated from depth trends in well logs. New theoretical equations were developed to predict the changes in elastic moduli upon substitution of pore-filling material. These equations reduce to Gassmann’s equations for the case of constant frame properties, low seismic frequencies, and fluid changes in the pore space. The new models also predict the change dissolution or precipitation of mineral, which cannot be described with the conventional Gassmann theory.

Mavko, G.↗

A continuum dislocation-based model of wedge microindentation of single crystals

Recent Electron Backscatter Diffraction (EBSD) experiments have revealed the emergence of heterogeneous dislocation microstructures forming under a wedge indenter in fcc crystals, where micro-meter dislocation patterns challenge the predictions of traditional models of plasticity. In order to explain the formation of these features and develop a relationship between the force-displacement curve and the dislocation substructure, here we present a model of wedge indentation based on the continuum theory of dislocations. The model accounts for large deformation kinematics through the multiplicative split of the deformation gradient tensor, where the incompatible plastic component of deformation results from the flux of dislocations on different and interacting slips systems. Constitutive equations for dislocation fluxes are determined from a dissipative variational principle. As a result, each dislocation density satisfies an initial-boundary value problem with convective-diffusive character, which is coupled to the macroscopic stress and displacement fields governing the deformation process. Solution to the self-consistent continuum formulation is found using the finite element method. Computer simulations mimic the experimental conditions of wedge micro-indentation experiments into Ni single-crystals used by Kysar et al. (2010a). A comparison of overall dislocation density distribution and macroscopic mechanical response shows good overall agreement with the experimental results in terms of the detailed features of dislocation patterns and lattice rotations as well as the macroscopic force-displacement response.

36 MATERIALS SCIENCE↗

4D-STEM of Beam-Sensitive Materials

We report that a scanning electron nanobeam diffraction, or 4D-STEM (four-dimensional scanning transmission electron microscopy), is a flexible and powerful approach to elucidate structure from "soft" materials that are challenging to image in the transmission electron microscope because their structure is easily damaged by the electron beam. In a 4D-STEM experiment, a converged electron beam is scanned across the sample, and a pixelated camera records a diffraction pattern at each scan position. This four-dimensional data set can be mined for various analyses, producing maps of local crystal orientation, structural distortions, crystallinity, or different structural classes. Holding the sample at cryogenic temperatures minimizes diffusion of radicals and the resulting damage and disorder caused by the electron beam. The total fluence of incident electrons can easily be controlled during 4D-STEM experiments by careful use of the beam blanker, steering of the localized electron dose, and by minimizing the fluence in the convergent beam thus minimizing beam damage. This technique can be applied to both organic and inorganic materials that are known to be beam-sensitive; they can be highly crystalline, semicrystalline, mixed phase, or amorphous. One common example is the case for many organic materials that have a π-π stacking of polymer chains or rings on the order of 3.4-4.2 Å separation. If these chains or rings are aligned in some regions, they will produce distinct diffraction spots (as would other crystalline spacings in this range), though they may be weak or diffuse for disordered or weakly scattering materials. We can reconstruct the orientation of the π-π stacking, the degree of π-π stacking in the sample, and the domain size of the aligned regions. This Account summarizes illumination conditions and experimental parameters for 4D-STEM experiments with the goal of producing images of structural features for materials that are beam-sensitive. We will discuss experimental parameters including sample cooling, probe size and shape, fluence, and cameras. 4D-STEM has been applied to a variety of materials, not only as an advanced technique for model systems, but as a technique for the beginning microscopist to answer materials science questions. It is noteworthy that the experimental data acquisition does not require an aberration-corrected TEM but can be produced on a variety of instruments with the right attention to experimental parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Removing Numerical Pathologies in a Turbulence Parameterization Through Convergence Testing

Abstract Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e., clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, nonlinear artificial diffusion employed for improving numerical stability can introduce unintended small‐scale features into the solution of the model equations. Smoothing the limiters and using linear artificial diffusion reduces the noise and restores the expected first‐order convergence in CLUBB's solutions. These model reformulations enhance our confidence in the trustworthiness of solutions from CLUBB by eliminating the unphysical oscillations in high‐resolution simulations. The improvements in the results at coarser, near‐operational grid spacing and timestep are also seen in cumulus cloud and dry turbulence tests. In addition, convergence testing is proven to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

58 GEOSCIENCES↗

Multi-scale gyrokinetic simulations of an Alcator C-Mod, ELM-y H-mode plasma

High fidelity, multi-scale gyrokinetic simulations capable of capturing both ion (${k}_{\theta }{\rho }_{s}\sim { \mathcal O }(1.0)$) and electron-scale (${k}_{\theta }{\rho }_{e}\sim { \mathcal O }(1.0)$) turbulence were performed in the core of an Alcator C-Mod ELM-y H-mode discharge which exhibits reactor-relevant characteristics. These simulations, performed with all experimental inputs and realistic ion to electron mass ratio (${({m}_{i}/{m}_{e})}^{1/2}=60.0$) provide insight into the physics fidelity that may be needed for accurate simulation of the core of fusion reactor discharges. Three multi-scale simulations and series of separate ion and electron-scale simulations performed using the GYRO code are presented. As with earlier multi-scale results in L-mode conditions, both ion and multi-scale simulations results are compared with experimentally inferred ion and electron heat fluxes, as well as the measured values of electron incremental thermal diffusivities—indicative of the experimental electron temperature profile stiffness. Consistent with the L-mode results, cross-scale coupling is found to play an important role in the simulation of these H-mode conditions. Extremely stiff ion-scale transport is observed in these high-performance conditions which is shown to likely play and important role in the reproduction of measurements of perturbative transport. In conclusion, these results provide important insight into the role of multi-scale plasma turbulence in the core of reactor-relevant plasmas and establish important constraints on the fidelity of models needed for predictive simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

WaterTransport in PEM Fuel Cells: Advanced Modeling, Material Selection, Testing and Design Optimization

Water management in Proton Exchange Membrane, PEM, Fuel Cells is challenging because of the inherent conflicts between the requirements for efficient low and high power operation. Particularly at low powers, adequate water must be supplied to sufficiently humidify the membrane or protons will not move through it adequately and resistance losses will decrease the cell efficiency. At high power density operation, more water is produced at the cathode than is necessary for membrane hydration. This excess water must be removed effectively or it will accumulate in the Gas Diffusion Layers, GDLs, between the gas channels and catalysts, blocking diffusion paths for reactants to reach the catalysts and potentially flooding the electrode. As power density of the cells is increased, the challenges arising from water management are expected to become more difficult to overcome simply due to the increased rate of liquid water generation relative to fuel cell volume. Thus, effectively addressing water management based issues is a key challenge in successful application of PEMFC systems. In this project, CFDRC and our partners used a combination of experimental characterization, controlled experimental studies of important processes governing how water moves through the fuel cell materials, and detailed models and simulations to improve understanding of water management in operating hydrogen PEM fuel cells. The characterization studies provided key data that is used as inputs to all state-of-the-art models for commercially important GDL materials. Experimental studies and microscopic scale models of how water moves through the GDLs showed that the water follows preferential paths, not branching like a river, as it moves toward the surface of the material. Experimental studies and detailed models of water and airflow in fuel cells channels demonstrated that such models can be used as an effective design tool to reduce operating pressure drop in the channels and the associated costs and weight of blowers and pumps to force air and hydrogen gas through the fuel cell. Promising improvements to materials structure and surface treatments that can potentially aid in managing the distribution and removal of liquid water were developed; and improved steady-state and freeze-thaw performance was demonstrated for a fuel cell stack under the self-humidified operating conditions that are promising for stationary power generation with reduced operating costs.

30 DIRECT ENERGY CONVERSION↗

Direct Molecular Evidence of Proton Transfer and Mass Dynamics at the Electrode–Electrolyte Interface

Proton transfer has been widely regarded as a key step in many electrochemical and biological processes. However, direct molecular evidence has long been lacking. In this work, we chose the electrochemical oxidation of acetaminophen (APAP) as a model system and utilized in situ liquid time-of-flight secondary ion mass spectroscopy (ToF-SIMS) to molecularly examine proton solvation and transfer in this process. In addition, we successfully captured and identified the transient radical intermediate, providing solid molecular evidence to resolve an important debate in electron transfer-proton transfer oxidation mechanism of APAP. Moreover, the potential-dependent behaviors of both inert ions and electroactive species during the dynamic potential scanning were chemically monitored in real time and the mass diffusion mechanism regarding the electroactive and non-electroactive species was revealed under polarized conditions. The results are well consistent with our computer simulation. The observations in this work greatly improved the comprehensive understanding of proton transfer and mass dynamics occurring at the electrode-electrolyte interface involved in complex electrochemical processes.

Wang, Jun-Gang↗

Electron acceleration in laboratory-produced turbulent collisionless shocks

Astrophysical collisionless shocks are among the most powerful particle accelerators in the Universe. Generated by violent interactions of supersonic plasma flows with the interstellar medium, supernova remnant shocks are observed to amplify magnetic fields and accelerate electrons and protons to highly relativistic speeds. In the well-established model of diffusive shock acceleration, relativistic particles are accelerated by repeated shock crossings. However, this requires a separate mechanism that pre-accelerates particles to enable shock crossing. This is known as the ‘injection problem’, which is particularly relevant for electrons, and remains one of the most important puzzles in shock acceleration. In most astrophysical shocks, the details of the shock structure cannot be directly resolved, making it challenging to identify the injection mechanism. In this letter we report results from laser-driven plasma flow experiments, and related simulations, that probe the formation of turbulent collisionless shocks in conditions relevant to young supernova remnants. We show that electrons can be effectively accelerated in a first-order Fermi process by small-scale turbulence produced within the shock transition to relativistic non-thermal energies, helping overcome the injection problem. Our observations provide new insight into electron injection at shocks and open the way for controlled laboratory studies of the physics underlying cosmic accelerators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗