Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic parameter”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Multiscale electro-chemo-mechanical model of high-capacity amorphous silicon anode material in sodium-ion batteries

The growing interest in sodium-ion batteries (SIBs) is fueled by their abundant resources and environmentally friendly nature, with amorphous silicon (a-Si) emerging as a promising anode material for enhancing capacity. However, the key challenge lies in sustaining reversible capacity during cycling. Here, in this work, we developed a multiscale electrochemical model incorporating an a-Si anode to elucidate the performance parameters of SIBs. Additionally, we integrated an electro-chemo-mechanical model at the particle level to capture stress generation, an essential factor in the degradation of high-capacity electrodes. Unlike existing models, our approach accounts for large-deformation chemo-mechanics at the particle scale and includes simulations under varying charge rates to explore multiscale behavior. The results reveal that coupled sodiation significantly prolongs complete cycling times and influences discharge dynamics, indicating that neglecting this coupling leads to an underestimation of actual capacity. Furthermore, we observed pronounced polarization effects at higher charge rates, resulting in heterogeneous stress distributions across the anode. With the identification of critical failure parameters for both active particles and binder materials, offering novel insights for mitigating degradation in high-capacity electrode systems.

Amorphous silicon anode↗
DOE OSTI2023

Exploiting electricity market dynamics using flexible electrolysis units for retrofitting methanol synthesis

Here we investigate the economic viability of integrating flexible electrolysis units to produce hydrogen in methanol synthesis processes. Specifically, we investigate whether this approach can help reduce methanol production costs by strategically exploiting dynamics of electricity markets. Our study integrates high-fidelity process simulations, optimization tools, and microkinetic modeling (informed by density functional theory) to conduct detailed techno-economic analyses and to compare performance against traditional processes that use hydrogen produced via steam-methane reforming (SMR). We also use this approach to estimate the levelized cost of hydrogen (LCOH) as a function of time-varying electricity prices (from day-ahead and real-time prices) and of key techno-economic parameters. Our results show that the proposed electrification framework is cost-competitive under certain electricity market conditions. Specifically, we find that, when the electrolysis system is operated in flexible mode (and can respond to dynamics of electricity markets), the associated electricity cost nearly collapses to zero. Conversely, when the unit is not flexible (and cannot respond to markets), the electricity cost comprises 60% of the total cost. Our results also reveal that the LCOH of the flexible electrolysis system participating in real-time electricity markets is 31% lower than the LCOH obtained from SMR. Overall, this indicates that exploiting the dynamics of electricity markets can make hydrogen production cost-competitive and this can lead to viable alternatives to electrify methanol production and other hydrogen-based processes.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Cross-Cutting Metrology Tools for In Operando Characterization of Carrier Dynamics in Photovoltaic Devices (Final Technical Report)

Understanding the nature of recombination and its dependence on defects and interfaces is essential for engineering materials and contacts for higher V oc and photovoltaic (PV) efficiency. Time-resolved photoluminescence (TRPL) has conventionally been used to evaluate recombination, but not all materials are strongly emissive or otherwise suitable. Time-resolved terahertz spectroscopy is a pump – probe method that presents valuable complementary information, wherein an optical pump pulse photoexcites carriers within the absorber and the transient photoconductivity is probed with a terahertz pulse. Until now, experimental constraints have prohibited the use of terahertz probes to interrogate full device stacks, and measurements were instead made on exfoliated films or films grown on unconventional substrates. However, interfaces are critical to the behavior of photoexcited carriers in solar cells, and the substrates themselves often influence the film growth and bulk properties. Therefore, it is important to probe the behavior of PV absorbers as close as possible to their normal operating conditions. Here we have developed cross-cutting metrology that enable in operando characterization of carrier dynamics and recombination mechanisms in working PV devices. Through a combination of complementary non-contact TRTS and TRPL experiments and modeling, we have obtained key parameters including photoexcited carrier bulk lifetime and interface and back surface recombination velocities with greater precision and accuracy than can be achieved with conventional TRPL alone. By varying pulsed photoexcitation conditions, the contributions of interface and bulk recombination mechanisms can be determined. We have developed and refined this characterization approach by investigating CdTe PVs, which is a well-established technology but with significant margin for further improvements in efficiency. Results were published in Journal of Applied Physics (DOI: 10.1063/5.0064730) and Proceedings of the IEEE PVSC (DOI: 10.1109/PVSC43889.2021.9518559). After validating the approach, we have applied it to correlate recombination rates to processing conditions and material properties, including composition and defects as well as interfaces and surface treatments. Such feedback can inform processing and design choices to enable higher PV efficiency, not only for CdTe, but also CdSe x Te 1-x , CIGS, perovskites, kesterites, and future technologies. A perspective on predicting solar cell performance from terahertz spectroscopy was published in Advanced Energy Materials (DOI: 10.1002/aenm.202102776). The cross-cutting metrology tools developed here will enable determination of the locus and mechanism of performance-limiting recombination in thin film PV devices. This feedback can guide the engineering of PV devices with higher V oc and efficiency, thereby leading to reductions in levelized cost of electricity to meet the SunShot 2030 target of $0.03/kWh for utility-scale PV. Lower costs will enable rapid expansion of renewable electricity generation that is nearly free of carbon emissions. Reducing carbon emissions of the electricity sector is one of the most important solutions to climate change. The project’s specific focus on CdTe PV can benefit First Solar, the US company that is the global leader in CdTe technology, through both technology and workforce development.

14 SOLAR ENERGY↗
DOE OSTI2022

Benchmarking and Exploring Parameter Space of the 2-Phase Bubble Tracking Model for Liquid Mercury Target Simulation

High intensity proton pulses strike the Spallation Neutron Source (SNS)’s mercury target to provide bright neutron beams. These strikes deposit extensive energy into the mercury and its steel vessel. Prediction of the resultant loading on the target is difficult when helium gas is intentionally injected into the mercury to reduce the loading and to mitigate the pitting damage on the vessel. A 2-phase material model that incorporates the Rayleigh-Plesset (R-P) model is expected to address this complex multi-physics dynamics problem by including the bubble dynamics in the liquid mercury. We present a study comparing the measured target strains in the SNS target station with the simulation results of the solid mechanics simulation framework. We investigate a wide range of various physical model parameters, including the number of bubble families, bubble size distribution, viscosity, surface tension, etc. to understand their impact on simulation accuracy. Our initial findings reveal that using 8-10 bubble families in the model renders a simulation strain envelope that covers the experimental ones. Further optimization studies are planned to predict the strain response more accurately.

Lin, Lianshan↗
DOE OSTI2024

A reactive electrochemomechanical theory for growth and remodeling of polyelectrolyte hydrogels and application to dynamic polymerization of $\text{DNA}$ hydrogels

Here this study develops a framework for growth and remodeling of active polyelectrolyte hydrogels that accounts for effects of compositional changes on the mechanical response. By developing a reactive electrochemomechanical theory, thermodynamical constraints upon reactive and remodeling processes are elucidated within a general framework that allows any number of chemical reactions to evolve the response of the gel and transfer mass and charge between constituents. Fully coupled, nonlinear constitutive relations are adopted for molar fluxes, allowing exploration of effects including cross-diffusion, electrophoresis, and electro-osmosis. A robust finite element implementation is developed in the open source FEBio software (febio.org) by exploiting an equivalence between electrochemomechanics and mixture theory. The implementation is verified against analytical solutions for free swelling, and a proper reduction to a prior chemomechanical theory is demonstrated for neutral gels swollen only by a solvent with no solutes. The theory and implementation are then applied to model the tunable large swelling achieved through dynamic polymerization of DNA crosslinkers seen in our recently developed experimental hydrogel system (Cangialosi et al., 2017). A novel constitutive model for reaction-driven evolution of the locking stretch λ L in a non-Gaussian mechanical free energy was developed, where the increasing concentration of DNA crosslinkers makes further swelling energetically favorable. With a single free parameter, excellent agreement was found between measured and predicted equilibrium swelling ratios. This study demonstrated the ability to extend the electrochemomechanical framework to include chemical reactions and composition-aware constitutive models, and showed that development of reactive models allows simulation of complex dynamic polymerization phenomena not treated before. The theoretical frame here can be further expanded in scope to incorporate additional non-ideal and nonlinear phenomena.

42 ENGINEERING↗
DOE OSTI2025

Lattice Structure and Dynamics of Sparse Molecular Crystals: OsO 4 and RuO 4

OsO 4 and RuO 4 are molecular oxides with unique tetrameric structures and rare +8 oxidation states. Accurately modeling their properties remains challenging for density functional theory (DFT) due to weak intertetramer interactions, which standard functionals fail to capture. Here, in this work, we show that the van der Waals (vdW)-corrected density functional (vdW-DF-optB86b) provides structural parameters that are much closer to experimental values than the standard generalized gradient approximation, with volume predictions that fall within the experimentally observed range. Phonon band structure analysis shows that the inclusion of vdW interactions stabilizes soft phonon modes, highlighting the importance of dispersion corrections for accurate predictions of lattice dynamics. Experimental measurements of the phonon density of states for OsO 4 , obtained via inelastic neutron scattering, demonstrate good agreement with our vdW-DF-optB86b calculations. These results validate OsO 4 and RuO 4 as valuable benchmarks for structural and vibrational calculations via vdW-corrected DFT methods and offer insights for studying the broader class of sparse molecular materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Dynamic Rupture Modeling of the M7.1, 2019 Ridgecrest, California, Earthquake

An important part of the simulation technique implemented in the strong motion prediction recipe developed in Japan, originally proposed by Irikura and Miyake (2011) is the earthquake rupture characterization that is based on established empirical scaling relationships between source parameters and seismic moment. These scaling relationships and rupture characterization are derived from recorded data and are tested by simulations. Recent analysis of results obtained with the Recipe for the 2016 M7.1 Kumamoto earthquake performed by several researchers, indicate that in order to reproduce near-fault long-period ground motions with permanent displacements, typical for earthquakes with surface rupturing, the recipe needs modifications. In particular, improvements are required in characterizing the rupture kinematics at shallow depths, above the seismogenic zone (Irikura et al., 2019). Modeling of rupture dynamics for recorded earthquakes, and numerical simulations for different rupture scenarios, provides physical constrains on slip evolution in time and space, which can be used in adjustments of empirical source scaling relationships for large crustal earthquakes. In this study we mainly focused on rupture dynamics modeling of the 2019 M7 Ridgecrest, California earthquake. The characterization of the kinematic rupture model derived from the dynamic rupture models of this earthquake can inform further refinements of rupture modeling proposed by the Recipe. Through different rupture scenarios we investigated potential relationships between characteristics of Long Period Ground Motion Generation Area (LMGA), used in the Recipe, and asperity depth and rupture initiation.

58 GEOSCIENCES↗
DOE OSTI2023

Robotic pendant drop: containerless liquid for μs-resolved, AI-executable XPCS

Abstract The dynamics and structure of mixed phases in a complex fluid can significantly impact its material properties, such as viscoelasticity. Small-angle X-ray Photon Correlation Spectroscopy (SA-XPCS) can probe the spontaneous spatial fluctuations of the mixed phases under various in situ environments over wide spatiotemporal ranges (10 −6 –10 3 s /10 −10 –10 −6 m). Tailored material design, however, requires searching through a massive number of sample compositions and experimental parameters, which is beyond the bandwidth of the current coherent X-ray beamline. Using 3.7-μs-resolved XPCS synchronized with the clock frequency at the Advanced Photon Source, we demonstrated the consistency between the Brownian dynamics of ~100 nm diameter colloidal silica nanoparticles measured from an enclosed pendant drop and a sealed capillary. The electronic pipette can also be mounted on a robotic arm to access different stock solutions and create complex fluids with highly-repeatable and precisely controlled composition profiles. This closed-loop, AI-executable protocol is applicable to light scattering techniques regardless of the light wavelength and optical coherence, and is a first step towards high-throughput, autonomous material discovery.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Enhancing Even Gas Distribution in Porous Media with Radial Flow

Packed beds, such as those used in thermal energy storage (TES) systems, typically use flow from one end to the other. This axial flow configuration leads to inefficiencies due to dispersion effects and high pressure drop. Instead of axial flow, this work proposes radial flow in packed beds, where a central tube provides flow that transports fluid from the center to the bed’s wall. Radial flow could be a promising solution to increase the efficiency of charging/discharging processes in TES systems. For instance, studies have been conducted on the thermal behavior during the charging process using radial flow and axial flow, and it was found that radial flow is better than axial flow in terms of thermal performance, where more energy can be stored in the storage tank during the charging process. However, the design of the radial pipe should be optimized to improve even flow distribution. In this work, numerical investigations on different designs were analyzed to enhance even flow gas distribution in the porous media when using this radial technique. This work shows the impact of different parameters on the even flow distribution into the packed bed for two different designs: 1) one radial tube at the center to provide the radial flow along with an annular tube at the wall to receive the flow, and 2) one radial tube at the center and four radial tubes at the wall. Air was used as fluid and 6 mm alumina beads as packing materials. Computational fluid dynamics (CFD) models in COMSOL Multiphysics were used to simulate the behavior of air flow through the piping and packed bed. It was found that the flow into the bed from the radial tube can be affected by different parameters: space between the holes, size/diameter of the holes, number of segments/zones, number of the holes in each zone and the length of each zone.

Beck, David M.↗
DOE OSTI2022

Provably accurate simulation of gauge theories and bosonic systems

Quantum many-body systems involving bosonic modes or gauge fields have infinite-dimensional local Hilbert spaces which must be truncated to perform simulations of real-time dynamics on classical or quantum computers. To analyze the truncation error, we develop methods for bounding the rate of growth of local quantum numbers such as the occupation number of a mode at a lattice site, or the electric field at a lattice link. Our approach applies to various models of bosons interacting with spins or fermions, and also to both abelian and non-abelian gauge theories. We show that if states in these models are truncated by imposing an upper limit Λ on each local quantum number, and if the initial state has low local quantum numbers, then an error at most ϵ can be achieved by choosing Λ to scale polylogarithmically with ϵ − 1 , an exponential improvement over previous bounds based on energy conservation. For the Hubbard-Holstein model, we numerically compute a bound on Λ that achieves accuracy ϵ , obtaining significantly improved estimates in various parameter regimes. We also establish a criterion for truncating the Hamiltonian with a provable guarantee on the accuracy of time evolution. Building on that result, we formulate quantum algorithms for dynamical simulation of lattice gauge theories and of models with bosonic modes; the gate complexity depends almost linearly on spacetime volume in the former case, and almost quadratically on time in the latter case. We establish a lower bound showing that there are systems involving bosons for which this quadratic scaling with time cannot be improved. By applying our result on the truncation error in time evolution, we also prove that spectrally isolated energy eigenstates can be approximated with accuracy ϵ by truncating local quantum numbers at Λ = polylog ( ϵ − 1 ) .

Tong, Yu↗
DOE OSTI2020

Discrete Fracture Network Modeling to Estimate Upscaled Parameters for the Topopah Spring, Lava Flow, and Tiva Canyon Aquifers at Pahute Mesa, Nevada National Security Site

This report describes the results of Discrete Fracture Network (DFN) simulations for the Topopah Spring Aquifer (TSA), Lava Flow Aquifer, and Tiva Canyon Aquifer (TCA), at Pahute Mesa on the Nevada National Security Site (NNSS), formerly the Nevada Test Site. The research focuses on calculating upscaled groundwater flow and contaminant transport parameters using DFNs generated according to fracture characteristics observed in the TSA, LFA and TCA at Pahute Mesa. The highly fractured and heterogeneous nature of these aquifers makes them candidates for stochastic DFN modeling of radionuclide transport on a small scale with subsequent upscaling. One hundred independent DFN realizations are generated for each aquifer, and the upscaled parameters for continuum simulations of subsurface flow and transport in fractured media at Pahute Mesa are calculated. Our goal is to implement a modeling approach that can translate parameters to larger-scale models that account for local-scale flow and transport processes, such as channelization of flow and transport along a few well connected, large fractures. Additionally, to simulate advective and advective-diffusive transport through the fracture networks, the Time Domain Random Walk (TDRW) approach is applied to account for matrix diffusion into a finite half-space. Moreover, a novel approach to calculate dynamic (active) fracture surface area to reflect flow channeling is implemented. This work will improve the representation of radionuclide transport processes in largescale, regulatory-focused models by providing estimates of hard-to-measure flow and contaminant transport parameters at large scales. In this report, we (1) show recent results of flow and transport simulations on multiple DFN realizations of the TSA, LFA, TCA; (2) discuss the resulting distributions of estimated upscaled parameters; (3) describe the estimation of upscaled parameters for an equivalent parallel-plate continuum model and (4) present a comparison between simulated transport from the equivalent continuum model and an actual DFN.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Optimal operational controls for power grid

Power grid operation is facing increased penetration of renewables and distributed energy resources at both transmission and distribution levels. This presents an increased level of randomness and uncertainty for its operational quality, where the challenges are mainly how controls and optimizations can be performed so that the operational quality of the power grid can be maintained. An interesting feature is also that we need to make a good use of uncertainties embedded in the swing equation from the renewable contributions rather than simply trying to minimize the impact of uncertainties. Furthermore, this requires effective solutions on system estimation, control, and optimization. Indeed, in past decades, work has been carried out on the development of advanced controls and optimization for such complicated dynamic systems. Examples are optimal operational control for both performance and stability enhancement considering random renewables penetration, frequency, and active power control for microgrids and parameter estimation for various parts of the power grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Affleck-Dine cogenesis of baryon and dark matter

We propose a mechanism for cogenesis of baryon and dark matter (DM) in the universe via the Affleck-Dine (AD) route. An AD field which breaks the lepton number symmetry, leads to the generation of lepton asymmetry by virtue of its cosmic evolution, which then gets transferred into lepton and dark sectors. While the lepton asymmetry gets converted into baryon asymmetry via sphalerons, the dark sector asymmetry leads to the final DM abundance with the symmetric part being annihilated away due to resonantly enhanced annihilation, which we choose to be provided by a gauged B – L portal. Stringent constraints from DM direct detection forces DM and B – L gauge boson masses to be light, in the few GeV ballpark. While a large portion of the model parameter space is already ruled out, the remaining parameter space is within sensitivity of laboratory as well as cosmology based experiments. The AD field also plays the role of inflaton with the required dynamics by virtue of its non-minimal coupling to gravity, consistent with observations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Cutting force and stability for inserted cutters using structured light metrology

This paper describes force and stability prediction for inserted end and face mills using a reverse engineering approach. Structured light scanning is used to identify the spatial coordinates of the points that define the multiple insert cutting edges. These points are used to extract the cutting edge radius and angle and are then incorporated in a time domain simulation that predicts cutting force and tool/workpiece deflection for user-selected operating parameters. Good agreement between predicted and measured cutting forces is demonstrated. The process stability is then predicted using the validated force model and measured structural dynamics.

Machining↗
DOE OSTI2022

Laser-driven, ion-scale magnetospheres in laboratory plasmas. I. Experimental platform and first results

We report magnetospheres are a ubiquitous feature of magnetized bodies embedded in a plasma flow. While large planetary magnetospheres have been studied for decades by spacecraft, ion-scale “mini” magnetospheres can provide a unique environment to study kinetic-scale, collisionless plasma physics in the laboratory to help validate models of larger systems. In this work, we present preliminary experiments of ion-scale magnetospheres performed on a unique high-repetition-rate platform developed for the Large Plasma Device at the University of California, Los Angeles. The experiments utilize a high-repetition-rate laser to drive a fast plasma flow into a pulsed dipole magnetic field embedded in a uniform magnetized background plasma. 2D maps of the magnetic field with high spatial and temporal resolution are measured with magnetic flux probes to examine the evolution of magnetosphere and current density structures for a range of dipole and upstream parameters. The results are further compared to 2D particle-in-cell simulations to identify key observational signatures of the kinetic-scale structures and dynamics of the laser-driven plasma. We find that distinct 2D kinetic-scale magnetopause and diamagnetic current structures are formed at higher dipole moments, and their locations are consistent with predictions based on pressure balances and energy conservation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Device Voltage Stress from Ground Leakage Current in Medium-Voltage Solid-State Transformer

Grounding related issues are critical for safe and reliable operation of solid-state transformer (SST) in medium-voltage (MV) applications, e.g., solar photovoltaic and energy storage integration, date center, electric vehicle fast charging, etc. This article presents for the first time the issue of additional device voltage stress due to grounding-loop current in current-source SST, using the soft-switching solid-state transformer (S4T) as an example. The S4T features single-stage isolated AC-AC, AC-DC, or DC-DC conversion with full-range ZVS, derived from flyback converter or current-source converter (CSC). However, the flyback operation for CSC-based SST means that magnetizing current flows through the reverse-blocking devices on only one side of the medium-frequency transformer (MFT) at a time. Then, the voltages across the devices, especially those on the other side of the MFT, can be influenced by parasitic current. A parasitic model of a modular S4T (M-S4T) prototype is developed from direct measurements and datasheets. Using the developed parasitic model and equivalent circuits, the causes of the voltage stress are analyzed. A voltage-stress mitigation scheme of connecting additional grounding capacitors is proposed. Damping resistors are also installed to damp out the grounding-loop resonance. A robust parameter design of the proposed scheme is given. The existence of the voltage stress issue and the effectiveness of the proposed scheme are verified experimentally on an MV SiC M-S4T prototype with inherent parameter variations among the five modules in the prototype. Here, both single-module and stacked-module operation are demonstrated during steady state and dynamic conditions up to 4 kV peak.

42 ENGINEERING↗
DOE OSTI2025

Detecting outbreaks using a spatial latent field

In this paper, we present a method for estimating the infection-rate of a disease as a spatial-temporal field. Our data comprises time-series case-counts of symptomatic patients in various areal units of a region. We extend an epidemiological model, originally designed for a single areal unit, to accommodate multiple units. The field estimation is framed within a Bayesian context, utilizing a parameterized Gaussian random field as a spatial prior. We apply an adaptive Markov chain Monte Carlo method to sample the posterior distribution of the model parameters condition on COVID-19 case-count data from three adjacent counties in New Mexico, USA. Our results suggest that the correlation between epidemiological dynamics in neighboring regions helps regularize estimations in areas with high variance (i.e., poor quality) data. Using the calibrated epidemic model, we forecast the infection-rate over each areal unit and develop a simple anomaly detector to signal new epidemic waves. Our findings show that anomaly detector based on estimated infection-rates outperforms a conventional algorithm that relies solely on case-counts.

Safta, Cosmin [Sandia National Laboratories (SNL-C↗
DOE OSTI2022

Computational Fluid Dynamics and Heat Transfer Modeling of a Dimpled Heat Exchanger

Multiphysics and analytical calculations were conducted for a heat exchanger with passive, natural circulation flow. A glycol/water working fluid convects the heat to a dimpled heat exchanger shell, which subsequently transfers the heat to the soil, which acts as the ultimate heat sink. Because the system is fully-passive, it is not subject to the expenses, maintenance, and mechanical breakdowns associated with moving parts. Density, heat capacity, and thermal conductivity material properties were measured for various soil samples, and subsequently included as input for the soil heat conduction model. The soil model was coupled to a computational fluid dynamics (CFD) heat exchanger model that included the dynamic Smagorinsky large eddy simulation and k- omega turbulence models. The analysis showed that the fluid dynamics and heat transfer models worked properly, albeit at a slow pace. Nevertheless, the coupled CFD/heat conduction simulation ran long enough to determine a key parameter—the amount of heat conducted from the heat exchanger to the ground. This unique performance value, along with experimental data, was used as input for stand-alone, fast-running CFD models, as well as boundaries to obtain solutions to partial differential equations for soil heat conduction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗