Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical model”

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,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multi-scale dynamics of magnetic flux tubes and inverse magnetic energy transfer

Here, we report on an analytical and numerical study of the dynamics of a three-dimensional array of identical magnetic flux tubes in the reduced-magnetohydrodynamic description of the plasma. We propose that the long-time evolution of this system is dictated by flux-tube mergers, and that such mergers are dynamically constrained by the conservation of the pertinent (ideal) invariants, viz. the magnetic potential and axial fluxes of each tube. We also propose that in the direction perpendicular to the merging plane, flux tubes evolve in a critically balanced fashion. These notions allow us to construct an analytical model for how quantities such as the magnetic energy and the energy-containing scale evolve as functions of time. Of particular importance is the conclusion that, like its two-dimensional counterpart, this system exhibits an inverse transfer of magnetic energy that terminates only at the system scale. We perform direct numerical simulations that confirm these predictions and reveal other interesting aspects of the evolution of the system. We find, for example, that the early time evolution is characterized by a sharp decay of the initial magnetic energy, which we attribute to the ubiquitous formation of current sheets. We also show that a quantitatively similar inverse transfer of magnetic energy is observed when the initial condition is a random, small-scale magnetic seed field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Fluctuation dynamo in a weakly collisional plasma

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Vacuum magnetic fields with exact quasisymmetry near a flux surface. Part 1. Solutions near an axisymmetric surface

While several results have pointed to the existence of exactly quasisymmetric fields on a surface (Garren & Boozer, Phys. Fluids B, vol. 3, 1991, pp. 2805–2821; 2822–2834; Plunk & Helander, J. Plasma Phys. , vol. 84, 2018, 905840205), we have obtained the first such solutions using a vacuum surface expansion formalism. We obtain a single nonlinear parabolic partial differential equation for a function $\eta$ such the field strength satisfies $B = B(\eta )$ . Closed-form solutions are obtained in cylindrical, slab and isodynamic geometries. Numerical solutions of the full nonlinear equations in general axisymmetric toroidal geometry are obtained, resulting in a class of quasihelical local vacuum equilibria near an axisymmetric surface. The analytic models provide additional insight into general features of the nonlinear solutions, such as localization of the surface perturbations on the inboard side. The local solutions thus obtained can be continued globally only for special initial surfaces.

Physics↗
DOE OSTI2022

Relativistic non-thermal particle acceleration in two-dimensional collisionless magnetic reconnection

Magnetic reconnection, especially in the relativistic regime, provides an efficient mechanism for accelerating relativistic particles and thus offers an attractive physical explanation for non-thermal high-energy emission from various astrophysical sources. I present a simple analytical model that elucidates key physical processes responsible for reconnection-driven relativistic non-thermal particle acceleration in the large-system, plasmoid-dominated regime in two dimensions. The model aims to explain the numerically observed dependencies of the power-law index $p$ and high-energy cutoff $\gamma _c$ of the resulting non-thermal particle energy spectrum $f(\gamma )$ on the ambient plasma magnetization $\sigma$ , and (for $\gamma _c$ ) on the system size $L$ . In this self-similar model, energetic particles are continuously accelerated by the out-of-plane reconnection electric field $E_{\rm rec}$ until they become magnetized by the reconnected magnetic field and eventually trapped in plasmoids large enough to confine them. The model also includes diffusive Fermi acceleration by particle bouncing off rapidly moving plasmoids. I argue that the balance between electric acceleration and magnetization controls the power-law index, while trapping in plasmoids governs the cutoff, thus tying the particle energy spectrum to the plasmoid distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effect of insulator end cap thickness on time-dependent Hartmann flow in a rotating mirror

We present a framework for analysing plasma flow in a rotating mirror. By making a series of physical assumptions, we reduce the magnetohydrodynamic (MHD) equations in a three-dimensional cylindrical system to a one-dimensional system in a shallow, cuboidal channel within a transverse magnetic field, similar to the Hartmann flow in ducts. We then solve the system both numerically and analytically for a range of values of the Hartmann number and calculate the dependence of the plasma flow speed on the thickness of the insulating end cap. We observe that the mean flow overshoots and decelerates before achieving a steady-state value, a phenomenon that the analytical model cannot capture. This overshoot is directly proportional to the thickness of the insulating end cap and the external electric field, with a weak dependence on the external magnetic field. Our simplified model can act as a benchmark for future simulations of the supersonic mirror device CMFX (centrifugal magnetic fusion experiment), which will employ more sophisticated physics and realistic magnetic field geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Approximating the particle distribution in rotating and tandem mirror traps

Steady-state distribution functions can be used to calculate stability conditions for modes, radiation energy losses and particle loss rates. Heuristic analytic approximations to these distributions can capture key behaviors of the true distributions such as the relative speeds of different transport processes while possessing computational advantages over their numerical counterparts. In this paper, we motivate and present a closed-form analytic model for a distribution of particles in a centrifugal or tandem mirror. We find that our model outperforms other known models in approximating numerical steady-state simulations outside of a narrow range of low confining potentials. We demonstrate the model’s suitability in the high confining potential regime for applications such as loss-cone stability thresholds, fusion yields and available energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Effects of Quantum Dot Loading on the Radioluminescence Efficiency in Quantum-Dot-Embedded Composites

Nanoparticle-embedded plastic scintillators are an emerging technology for fast, large-area, high-resolution radiation detection and imaging. Here, this study investigates the properties of such composites, focusing on the effects of the quantum dot (QD) concentration on the radioluminescence (RL) intensity, spectra, and dynamics. Experiments using CdSe/CdS QDs in a polymer reveal a superlinear increase in RL with the QD concentration despite optical losses from inner filtering and interparticle interactions. When corrected for inner filtering, RL shows a quadratic concentration dependence, consistent with simple analytical models of improving the secondary electron capture. Practically, the benefits of high QD concentrations are muted by optical losses, but the findings apply to other systems with insulating hosts. In addition to manipulating emission for large effective Stokes shifts, future improvements may come from hosts with higher stopping power and better charge transport, which enable more effective funneling of excitations but without concomitant optical losses associated with high nanoparticle concentrations.

Auger recombination↗
DOE OSTI2019

Native Void Space for Maximum Volumetric Capacity in Silicon-Based Anodes

Volumetric energy density is considered a primary factor in developing high-energy batteries. Despite its significance, less efforts have been devoted to its improvement. Silicon-based materials have emerged as next-generation anodes for lithium-ion batteries due to their high specific capacity. However, their volumetric capacities are limited by the volume expansion rate of silicon, which restricts mass loading in the electrodes. To address this challenge, we introduce porous silicon templated from earth-abundant minerals with native internal voids, capable of alleviating volumetric expansion during repeated cycles. In situ transmission electron microscopy analysis allows the precise determination of the expansion rate of silicon, thus presenting an analytical model for finding the optimal content in silicon/graphite composites. The inner pores in silicon reduce problems associated with its expansion and allow higher silicon loading of 42% beyond the conventional limitations of 13–14%. Consequently, the anode designed in this work can deliver a volumetric capacity of 978 mAh cc–1. Thus, suppressing volume expansion with natural abundant template-assisted materials opens new avenues for cost-effective fabrication of high volumetric capacity batteries.

Yeom, Su Jeong↗
DOE OSTI2021

3D Printed Nickel–Molybdenum-Based Electrocatalysts for Hydrogen Evolution at Low Overpotentials in a Flow-Through Configuration

Three-dimensional (3D) printed, hierarchically porous nickel molybdenum (NiMo) electrocatalysts were synthesized and evaluated in a flow-through configuration for the hydrogen evolution reaction (HER) in 1.0 M KOH(aq) in a simple electrochemical H-cell. 3D NiMo electrodes possess hierarchically porous structures because of the resol-based aerogel precursor, which generates superporous carbon aerogel as a catalyst support. Relative to a traditional planar electrode configuration, the flow-through configuration allowed efficient removal of the hydrogen bubbles from the catalyst surface, especially at high operating current densities, and significantly decreased the overpotentials required for HER. An analytical model that accounted for the electrokinetics of HER as well as the mass transport with or without the flow-through configuration was developed to quantitatively evaluate voltage losses associated with kinetic overpotentials and ohmic resistance due to bubble formation in the porous electrodes. Furthermore, the chemical composition, electrochemical surface area (ECSA), and roughness factor (RF) were also systematically studied to assess the electrocatalytic performance of the 3D printed, hierarchically porous NiMo electrodes. An ECSA of 25163 cm2 was obtained with the highly porous structures, and an average overpotential of 45 mV at 10 mA cm–2 was achieved over 24 h by using the flow-through configuration. The flow-through configuration evaluated in the simple H-cell achieved high electrochemical accessible surface areas for electrochemical reactions and provided useful information for adaption of the porous electrodes in flow cells.

08 HYDROGEN↗
DOE OSTI2025

Off-Equilibrium Reactivity of Boron-Enriched Metal Diboride Surfaces in Electroreduction Conditions

Boron-based materials, featuring B-dependent reactivity and diverse phases, are emerging as promising catalyst systems. However, the catalytic mechanism on many borides remains poorly understood due to complex surface reconstructions under reaction conditions. Here, we investigate the MoB 2 surface in conditions of hydrogen evolution reaction in acidic media, using grand canonical global optimization, grand canonical density functional theory, ab initio molecular dynamics, free energy surface sampling, and an analytical model for electrochemical barrier evaluation. We propose a boron-enrichment strategy to tune the surface reactivity of the hexagonal face of MoB 2 . We reveal the dynamic nature of the B-enriched surface under H coverage and kinetic trapping of the system in the metastable regime with an extensive examination of the deactivation pathways. The metastable center B site on B-enriched surfaces, featuring buckled-up configuration and a usual relaxation effect, is found to be highly active toward HER via the Volmer–Heyrovsky mechanism. In conclusion, this work demonstrates how off-equilibrium behaviors can arise from the interplay between adsorbate coverage and surface reconstruction on a seemingly simple surface, and we present a theoretical framework and computational workflows to address these behaviors, along with other realistic complexities, in kinetics simulations.

Adsorption↗
DOE OSTI2025

Tutorial: Machine-Learning-Based CREASE-2D Analysis of 2D SAXS Profiles to Characterize Anisotropic Nanostructures in Soft Materials

We present a tutorial to guide users on how to extend the Computational Reverse Engineering Analysis of Scattering Experiments-2D (CREASE-2D) framework to interpret their experimental two-dimensional small-angle scattering (SAS) data from soft materials (e.g., polymers, peptide amphiphiles, biomolecular fibrils). Unlike most traditional SAS analysis approaches, which typically rely on azimuthally averaged onedimensional (1D) profiles, CREASE-2D utilizes the complete 2D scattering profile to reveal information about anisotropy in the structure. In past applications, CREASE has provided insights into complex structural features, including the cross-sectional shapes of assembled nanostructures and dispersity in these features, which are difficult to discern with existing analytical models. While (1D- ) CREASE has been applied to SANS and SAXS data, this tutorial shares the steps for implementing CREASE-2D using an example of a dipeptide solution system, for which we have SAXS data. We present details for these steps involved in using CREASE-2D to interpret SAXS profiles: how to preprocess SAXS data, define relevant structural features, generate three-dimensional real-space structures for specific values of these features, train a machine learning (ML) surrogate model to predict scattering profiles for given structural features, and optimize these features using genetic algorithms (GA). Then, we use these steps to interpret complex 2DSAXS data collected from dipeptide solutions that, in microscopy images, exhibit nanoscale structures that could be elliptical tubes/ flat tapes/cylinders or a combination of these cross sections. Open-source codes, computational hardware, and software requirements, as well as the strengths and limitations of this protocol, are also presented. We expect researchers working with (soft) biomaterials, peptide amphiphiles, amphiphilic polymer solutions, polymer nanocomposites, and blends of particles/polymers will find this CREASE-2D method and this tutorial of use.

CREASE↗
DOE OSTI2022

Polariton Dynamics in Two-Dimensional Ruddlesden–Popper Perovskites Strongly Coupled with Plasmonic Lattices

Strong coupling between light and matter can produce hybrid eigenstates known as exciton-polaritons. Although polariton dynamics are important photophysical properties, the relaxation pathways of polaritons in different coupling regimes have seen limited attention. This paper reports the dynamics of hybridized states from 2D Ruddlesden-Popper perovskites coupled to plasmonic nanoparticle lattices. The open cavity architecture of Al lattices enables the coupling strength to be modulated by varying either the lead halide perovskite film thickness or the superstrate refractive index. Both experiments and finite-difference time-domain simulations of the optical dispersion diagrams showed avoided crossings that are a signature of strong coupling. Furthermore, our analytical model also elucidated the correlation between the exciton/plasmon mixing ratio and polariton coupling strength. Using fs-transient absorption spectroscopy, we found that both the upper and lower polaritons have shorter lifetimes than the excitons and that polaritons can show faster excited state dynamics when they have access to additional energy transfer channels.

2D perovskite↗
DOE OSTI2025

Chromium versus Aluminum: Impact of Nickel Alloy Composition and Interfacial Kinetics on High-Temperature Passivating Oxide Formation

High-temperature corrosion resistance depends critically on the formation of a passivating surface oxide, which is highly sensitive to alloy composition and structure. Such details often elude experimental investigation, and simplified analytical models fail to provide a truly chemical view of passivating oxide evolution. Here, we explicitly compare the fundamental chemistry of Cr and Al as prototypical passivating elements in Ni alloys by directly simulating competing reaction and diffusion processes within the oxide film using kinetic Monte Carlo and density functional theory. We find that the origin and expression of passivating behavior during early-stage thermal oxidation are qualitatively different between the two alloy systems. Ni–Cr alloys feature a sudden onset of passivation associated with a sharp phase transition upon Cr enrichment that directly couples oxidation kinetics to phase transformation behavior. In contrast, Ni–Al alloys display more continuous oxide phase variation with Al enrichment, ultimately resulting in a lower composition threshold for passivation and a thinner passivating layer. In addition, we elucidate the nonobvious role of metal exchange within the alloy near the oxide boundary, which fundamentally alters film composition and passivating behavior. Furthermore, our results have key implications for engineering improved corrosion-resistant alloys, both in terms of compositional variation and processing.

Alloys↗
DOE OSTI2025

Native Top-Down Mass Spectrometry Characterization of Model Integral Membrane Protein Bacteriorhodopsin

Bacteriorhodopsin (bR) from Halobacterium salinarum has been a model system for structural biology and is a structural template for the characterization of membrane G-protein couple receptors (GPCRs) in particular. Here, in this study, wild-type bacteriorhodopsin and two single-residue mutants were characterized by native top-down mass spectrometry (nTD-MS) with Orbitrap-based high-energy collision dissociation (HCD) and electron capture dissociation (ECD). After in-source dissociation ejected the membrane protein from detergent micelles, high-resolution native MS measurement allowed for identification of multiple proteoforms as well as lipid-bound forms. Further top-down MS measurements by HCD produced a large number of product ions for in-depth sequencing and unambiguous localization of post-translational modifications. For the first time, native TD-MS with ECD was used to characterize an integral membrane protein. ECD yielded fragments originating from all helices and loop regions, even accessing a sequence stretch that HCD could not. Combining HCD and ECD fragmentation patterns significantly enhanced the sequence coverage of bR. We propose bR to be a model analyte for testing nTD-MS performance for membrane proteins.

crystal cleavage↗
DOE OSTI2019

Helium in Stream Water as a Volcanic Monitoring Tool

Abstract We show that synoptic sampling of streams can be used to characterize volcanic volatiles in groundwater over large spatial scales. Synoptic sampling of dissolved noble gases, 222 Rn, major ions, and stream discharge was carried out along a 30 km reach of the Gibbon River, near Norris Geyser Basin in Yellowstone National Park, USA. Groundwater discharge location, volume, and composition were estimated by constrained calibration of a stream flow and solute transport model. Estimated groundwater composition from stream modeling was compared to shallow groundwater concentrations measured in nearby springs. 3 He, 222 Rn, and Cl signatures in the Gibbon River are indicative of groundwater discharge with a volcanic signature along the study reach. Stream water noble gas isotopic composition has similar isotopic mixing patterns to springs. The model‐estimated composition of groundwater discharging to the Gibbon agrees well with observed groundwater composition from nearby springs for all modeled analytes. We present the first observations of elevated mantle helium in stream water and show that stream water can be used as a convenient collection point to estimate spatially distributed groundwater composition and to monitor changes in volatile flux over large spatial areas. These results offer the possibility that stream surveys in volcanic terrain could be a new method for distributed volcanic monitoring at the catchment scale and beyond.

Gardner, W. Payton↗
DOE OSTI2021

Enabling the Application of Large Footprint Open-Bottom Permeameters through New Shape Factors

Open-bottom permeameters (OBP) are commonly used for estimating the hydraulic properties of sediment beds under water bodies and fluxes across the sediment-water interface. The specific shape factor for an OBP flow situation is necessary for data interpretation. It is determined by the ratio of the permeameter tube radius to the tube's penetration depth within the sediment. Previous permeameter designs favored small radii which allowed for deeper penetration and stability of the tube. The shape factor is provided here for large footprint OBPs that allows for shallow penetration and greater radius. Through a numerical and analytical modeling study of the OBP hydrodynamics, previous knowledge on the OBP shape factor is generalized to the entire range of radius-to-penetration ratios that might be practically applied. The results reveal the advantages of the large footprint design, including better sensitivity due to the ability to impose larger discharge over a shorter test duration.

42 ENGINEERING↗
DOE OSTI2021

Injection Parameters That Promote Branching of Hydraulic Cracks

Fluid injection into rock formations can either produce complex branched hydraulic fractures, create simple planar fractures, or be dominated by porous diffusion. Currently, the optimum injection parameters to create branched fractures are unknown. We conducted repeatable hydraulic fracturing experiments using analog-rock samples with controlled heterogeneity to quantify the fluid parameters that promote fracture branching. A large range of injection rates and fluid viscosities were used to investigate their effects on induced fracture patterns. Paired with a simple analytical model, our results identify the threshold at which fracture transitions from an isolated planar crack to branched cracks when closed natural fractures exist. These results demonstrate that this transition can be controlled by injection rate and fluid viscosity. In relation to the field practices, the present model predicts slickwater and lower viscosity fluid injections promote fracture branching, with the Marcellus shale used as an example.

58 GEOSCIENCES↗
DOE OSTI2022

Leakage From Coexisting Geologic Forcing and Injection‐Induced Pressurization: A Semi‐Analytical Solution for Multilayered Aquifers With Multiple Wells

Abstract Abnormal fluid pressures (above or below hydrostatic pressure) can develop and persist in sedimentary basins. The common occurrence of abnormal pressures may cause challenges for project permitting of geological carbon sequestration (GCS), particularly in reservoirs with pre‐injection overpressure. The leaky wells that may exist in some sedimentary basins can provide flow paths between deep brine aquifers and shallower freshwater aquifers. Pre‐injection relative overpressures can cause brine leakage through leaky wells even before any injection occurs. The tendency for flow through leaky wells is coupled with the process of pressure dissipation that occurs through aquitards. Specifically, with non‐zero permeability, aquitards can dissipate pressure over large areal extents and thereby reduce leakage rates through leaky wells. This study presents development of a semi‐analytical solution for hydraulic head and brine leakage in multilayered aquifer–aquitard systems with geologic pressure forcing. The geologic forcing that causes abnormal pressures in the multilayered system can coexist with any number of injection, extraction, and leaky wells that also affect fluid pressure. The semi‐analytical model is applied to explore how leakage through leaky wells varies as functions of pressurization rate, along with aquitard and leaky well properties in an overpressured multilayered system. The results show that although injection‐induced pressures can dissipate rapidly through suitably permeable aquitards, coexisting geologic forcing may create sustained rates of brine leakage into freshwater aquifers through leaky wells. In GCS, a very low‐permeability aquitard with high capillary entry pressure to free‐phase CO 2 is desired to serve as the caprock to prevent leakage of CO 2 from the storage reservoir. Nevertheless, the results from this study show that the brine leakage impact to shallow freshwater aquifers through leaky wells might substantially decrease with increasing aquitard permeability values, as long as small aquitard permeability and high capillary entry pressure serve to prevent CO 2 leakage.

Cihan, A.↗