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.
Solitary zonal structures have recently been identified in gyrokinetic simulations of subcritical drift-wave (DW) turbulence with background shear flows.
Due to high fossil fuel imports, high electricity rates, an unreliable electricity grid, and 100% renewable electricity goals, demand response can play a crucial role for the Puerto Rico electricity grid. The Puerto Rico Demand Response Impact and Forecast Tool (PR-DRIFT) is a spreadsheet-based tool in which users can estimate the potential impacts of demand response, energy efficiency, and VRE and storage adoption in Puerto Rico from 2021 through 2040. The tool includes projections for solar, wind, and battery adoption based on released RFPs and energy targets (including Act-17 2019) to generate a projected net load profile for each hour through 2040. Based on user inputs and default assumptions, the tool also projects load profile impacts of demand response and energy efficiency, and is specifically focused on highlighting projected demand response technical potential. This demand response technical potential can be used by local utilities, regulators, program administrators, or researchers to help design demand response programs for larger impact.
Human-to-swine transmission of influenza A (H3N2) virus occurs repeatedly and plays a critical role in swine influenza A virus (IAV) evolution and diversity. Human seasonal H3 IAVs were introduced from human-to-swine in the 1990s in the United States and classified as 1990.1 and 1990.4 lineages; the 1990.4 lineage diversified into 1990.4.A–F clades. Additional introductions occurred in the 2010s, establishing the 2010.1 and 2010.2 lineages. Human zoonotic cases with swine IAV, known as variant viruses, have occurred from the 1990.4 and 2010.1 lineages, highlighting a public health concern. If a variant virus is antigenically drifted from current human seasonal vaccine (HuVac) strains, it may be chosen as a candidate virus vaccine (CVV) for pandemic preparedness purposes. We assessed the zoonotic risk of US swine H3N2 strains by performing phylogenetic analyses of recent swine H3 strains to identify the major contemporary circulating genetic clades. Representatives were tested in hemagglutination inhibition assays with ferret post-infection antisera raised against existing CVVs or HuVac viruses. The 1990.1, 1990.4.A, and 1990.4.B.2 clade viruses displayed significant loss in cross-reactivity to CVV and HuVac antisera, and interspecies transmission potential was subsequently investigated in a pig-to-ferret transmission study. Strains from the three lineages were transmitted from pigs to ferrets via respiratory droplets, but there were differential shedding profiles. These data suggest that existing CVVs may offer limited protection against swine H3N2 infection, and that contemporary 1990.4.A viruses represent a specific concern given their widespread circulation among swine in the United States and association with multiple zoonotic cases.
Partial shading of a solar module can induce a set of cells within the module to operate under reverse bias. Studies have shown that metal halide perovskite solar cells with a wide variety of compositions and contacts exhibit interesting behavior in reverse bias that includes both reversible performance loss and non-reversible degradation. In this paper, an advanced drift-diffusion approach incorporating an electrochemical term to explain the short-circuit, open circuit and fill factor losses that are experimentally measured after prolonged reverse bias is used. It is shown that holes can tunnel into the perovskite due to sharp band bending near the contact, accumulate within the bulk of the perovskite absorber, and trigger the oxidation of halides to form neutral halogens. The density of neutral halogens is much higher in reverse bias because there are hardly any electrons available to reduce the iodine. The resulting halogens act as bulk recombination centers. While the interstitial halogen density does decay when the cell is operated in forward bias, permanent degradation can occur if the iodine diffuses out of the perovskite layer. Finally, the ways in which changing parameters such as the mobile ion density or the series resistance at the contact can influence device performance and stability are discussed.
Recently the drift-reduced Landau fluid six-field turbulence model within the BOUT++ framework has been upgraded. In particular, this new model employs a new normalization, adds a volumetric flux-driven source option, the Landau fluid closure for parallel heat flux and a Laplacian inversion solver which is able to capture $n$ = 0 axisymmetric mode evolution in realistic tokamak configurations. As we report here, these improvements substantially extended model's capability to study a wider range of tokamak edge phenomena, and are essential to build a fully self-consistent edge turbulence model capable of both transient (e.g., ELM, disruption) and transport time-scale simulations.
We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.
Here we introduce fluctuating hydrodynamics approaches on surfaces for capturing the drift-diffusion dynamics of particles and microstructures immersed within curved fluid interfaces of spherical shape. We take into account the interfacial hydrodynamic coupling, traction coupling with the surrounding bulk fluid, and thermal fluctuations. For fluid-structure interactions, we introduce Immersed Boundary Methods (IBM) and related Stochastic Eulerian-Lagrangian Methods (SELM) for curved surfaces. We use these approaches to investigate the statistics of surface fluctuating hydrodynamics and microstructures. For velocity autocorrelations, we find characteristic power-law scalings $τ^{-1}, τ^{-2}$, and plateaus can emerge. This depends on the physical regime associated with the geometry, surface viscosity, and bulk viscosity. This differs from the characteristic $τ^{-3/2}$ scaling for bulk three dimensional fluids. We develop theory explaining these observed power-laws associated with time-scales for dissipation within the fluid interface and coupling to the surrounding fluid. We then use our introduced methods to investigate a few example systems and roles of hydrodynamic coupling and thermal fluctuations including for the kinetics of passive particles and active microswimmers in curved fluid interfaces.
In this work we consider a relativistic drift-kinetic model for runaway electrons along with a Fokker–Planck operator for small-angle Coulomb collisions, a radiation damping operator, and a secondary knock-on (Boltzmann) collision source. Here, we develop a new scalable fully implicit solver utilizing finite volume and conservative finite difference schemes and dynamic mesh adaptivity. A new data management framework in the PETSc library based on the p4est library is developed to enable simulations with dynamic adaptive mesh refinement (AMR), distributed memory parallelization, and dynamic load balancing of computational work. This framework and the runaway electron solver building on the framework are able to dynamically capture both bulk Maxwellian at the low-energy region and a runaway tail at the high-energy region. To effectively capture features via the AMR algorithm, a new AMR indicator prediction strategy is proposed that is performed alongside the implicit time evolution of the solution. This strategy is complemented by the introduction of computationally cheap feature-based AMR indicators that are analyzed theoretically. Numerical results quantify the advantages of the prediction strategy in better capturing features compared with nonpredictive strategies; and we demonstrate trade-offs regarding computational costs. The robustness with respect to model parameters, algorithmic scalability, and parallel scalability are demonstrated through several benchmark problems including manufactured solutions and solutions of different physics models. We focus on demonstrating the advantages of using implicit time stepping and AMR for runaway electron simulations.
We present the second prototype of a time-to-digital (TDC) ASIC for the upgrade of the ATLAS Monitored Drift Tube (MDT) detector for High-Luminosity LHC operations. Compared to the first prototype, triple modular redundancy has been implemented for the configuration and flow control logic. Here, the total power consumption is increased by less than 10 mW while achieving the same time resolution and channel uniformity. A mini-DAQ system has been built to verify the front-end electronics chain with the new prototype together with other ASICs and boards in triggered mode. Cosmic ray tests with a small-diameter MDT chamber indicate that the configuration and data transmission of the readout electronics perform well. It is expected that this prototype design will be used in the final production.
We report low-frequency (such as 16.77 MHz) RF bunchers are widely used in RF accelerator systems for longitudinal compression of pulses into a single RF bucket, which increases instantaneous beam intensity for time-dependent studies. In this study, the dependency of capture into a 201.25 MHz Drift Tube Linac (DTL) was measured as a function of gap voltage for a 16.77 MHz buncher on chopped H - beam (approximately 25 ns at 750 keV, 10 mA peak current). The multiparticle code PARMILA was used to simulate the phase-space distribution of the 10 mA, 750 keV, H - beam at the entrance to DTL with a wide range of the Low-Frequency Buncher (LFB) field (10 kV to 35 kV). The measurement and simulation indicated that the DTL capture could be dilute (reduced) for a non-optimized buncher field to a pre-configured beamline geometry. The data shows that changing the bunch field while keeping the incoming beam current and energy constant does not significantly alter the beam’s emittance. However, downstream beam capture into the DTL is changed for a non-optimized phase-space bunching distance with the buncher field.
New front-end electronics including ASICs and FPGA boards are under development for the ATLAS Monitored Drift Tube (MDT) detector to handle the large data rates and harsh environment expected at high-luminosity LHC runs. A mobile Data Acquisition (miniDAQ) system is designed to perform integration tests of these front-end electronics. In addition, it will be used for surface commissioning of 96 small-radius MDT (sMDT) chambers and for integration and commissioning of new front-end electronics on the present ATLAS MDT chambers. Details of the miniDAQ hardware and firmware are described in this article. Here, the miniDAQ system is also used to read out new front-end electronics on an sMDT prototype chamber using cosmic muons and results obtained are shown.
The kinetic stability of collisionless, sloshing beam-ion (45° pitch angle) plasma is studied in a three-dimensional (3-D) simple magnetic mirror, mimicking the Wisconsin high-temperature superconductor axisymmetric mirror experiment. The collisional Fokker–Planck code CQL3D-m provides a slowing-down beam-ion distribution to initialize the kinetic-ion/fluid-electron code Hybrid-VPIC, which then simulates free plasma decay without external heating or fuelling. Over 1 – 10 μs, drift-cyclotron loss-cone (DCLC) modes grow and saturate in amplitude. The DCLC scatters ions to a marginally stable distribution with gas-dynamic rather than classical-mirror confinement. Sloshing ions can trap cool (low-energy) ions in an electrostatic potential well to stabilize DCLC, but DCLC itself does not scatter sloshing beam-ions into the said well. Instead, cool ions must come from external sources such as charge-exchange collisions with a low-density neutral population. Manually adding cool ∼1keV ions improves beam-ion confinement several-fold in Hybrid-VPIC simulations, which qualitatively corroborates prior measurements from real mirror devices with sloshing ions.
Lithium hydride (LiH) is a unique, ionic compound with applications in a variety of industries. Unfortunately, LiH is very reactive toward H 2 O even at ppm levels, forming oxide (Li 2 O) and hydroxide (LiOH) corrosion layers while outgassing H 2 . An effective means to eliminate unwanted outgassing is vacuum-heating to convert LiOH into Li 2 O, although subsequent re-exposure to moisture during transport/handling reconverts some Li 2 O back to LiOH. A corrosion growth model for previously vacuum-baked LiH is necessary for long-term prediction of the hydrolysis of LiH. In this work, a para-linear hydroxide corrosion growth model is proposed for the reaction of previously vacuum-baked LiH samples with moisture. This model, composed of two competing diffusion reaction fronts at the LiOH/Li 2 O and Li 2 O/LiH interfaces, is validated experimentally by subjecting a previously vacuum-baked polycrystalline LiH sample to 35 ppm of H 2 O at room temperature while monitoring the corrosion growth as a function of time with diffuse-reflectance infrared Fourier transform (DRIFT) spectroscopy. The para-linear growth model for the hydrolysis of previously vacuum-baked LiH proposed in this report can also serve as a template for the hydrolysis of other hygroscopic oxides grown on metal or metal hydride substrates.
Experiments of in-situ DRIFTS combined with modulation excitation (ME) spectroscopy showed a rich surface chemistry associated with the adsorption of CO 2 on nanocubes and nanospheres of ceria. The nanocubes exposed faces with a (100) orientation, with the edges and corners displaying (110) and (111) orientations, respectively. Here, the nanospheres mainly contained ceria (111) and (110) planes. DFT calculations showed that CO 2 is a multidentate adsorbate on ceria that can undergo changes in its bonding configuration depending on the chemical environment. At 250 °C, a temperature typically used for the conversion of CO 2 into oxygenates, alkanes and olefins, CO 2 reacted with O centers or OH groups present on the nanocubes and nanospheres to yield bi- and tri-dentate carbonates, hydroxycarbonates, and formates. Both nanostructures were highly reactive and a dynamic equilibrium was established: carbonate species were rapidly generated upon the injection of CO 2 and they decomposed upon the removal of CO 2 from the gas phase. In the case of the ceria nanocubes, the adsorption/desorption processes were essentially reversible, opening the door to catalytic transformations. A larger concentration of defects in the ceria nanospheres led to strongly bound carbonates and formates that may be spectators, site blockers, or surface modifiers in catalytic processes. In the ME studies, additional intermediates were detected, and it was clear that the response of surface species to the presence/absence of CO 2 was highly dependent on the morphology of the ceria nanostructures.
Single-site catalysts have drawn broad attention in catalysis because of their maximum atomic utilization and unique catalytic performance. Early work in our group has shown a 40-fold higher activity of methanol decomposition over single-site Pt 1 /CeO 2 catalyst than CeO 2 supported 2.5 nm Pt nanoparticles, while a molecular-level understanding of such enhancement is lacking. Herein, the reaction mechanism of methanol decomposition over Pt 1 /CeO 2 was carefully investigated using in situ DRIFTS, and a reaction pathway was proposed. Methanol molecules were dissociatively adsorbed on nanoceria support first, followed by the diffusion of as-formed methoxy species onto Pt single sites where the dehydrogenation occurs and results in the weakly bonded CO. The ease of methanol dissociative adsorption on nanoceria support and the tailored electronic property of Pt 1 via the metal-support interaction are believed to be strongly correlated with the high activity of Pt 1 /CeO 2 .
Generation and propagation of lower hybrid drift wave (LHDW) near the electron diffusion region (EDR) during guide field reconnection at the magnetopause is studied with data from the Magnetospheric Multiscale mission and a theoretical model. Inside the current sheet, the electron beta (β e ) determines which type of LHDW is excited. Inside the EDR, where the electron beta is high (β e ~5), the long-wavelength electromagnetic LHDW is observed propagating obliquely to the local magnetic field. In contrast, the short-wavelength electrostatic LHDW, propagating nearly perpendicular to the magnetic field, is observed slightly away from the EDR, where β e is small (~0.6). These observed LHDW features are explained by a local theoretical model, including effects from the electron temperature anisotropy, finite electron heat flux, electrostatics, and parallel current. The short-wavelength LHDW is capable of generating significant drag force between electrons and ions.
Abstract In this paper, we derive a new drift‐kinetic method for estimating gradients in the plasma properties through a velocity space distribution at a single point. The gradients are intrinsically related to agyrotropic features of the distribution function. This method predicts the gradients in the magnetized distribution function and can predict gradients of arbitrary moments of the gyrotropic background distribution function. The method allows for estimates on density and pressure gradients on the scale of a Larmor radius, proving to resolve smaller scales than any method currently available to spacecraft. The model is verified with a set of fully kinetic VPIC particle‐in‐cell simulations.
Here, we investigate properties of lower hybrid drift waves (LHDWs) near and inside the electron diffusion regions in 17 magnetopause and 9 magnetotail reconnection events using Magnetospheric MultiScale (MMS) mission observations. Our analysis show that LHDW type depend on the electron beta, as electron beta increases LHDWs become more electromagnetic in nature. The energy transfer from electromagnetic fields to particles is higher in electrostatic LHDWs and it is largely in parallel direction with respect to the local magnetic field. Linear dispersion analysis shows that electrostatic LHDWs are perpendicular propagating while electromagnetic waves propagate in oblique directions and the normalized wavenumber of all LHDW types falls within 0.5–0.8 range. A simple estimate on the LHDW nonlinear saturation suggests a possibly important roles played by these waves in supporting the reconnection electric field.