Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Simulation”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Analysis of heat transfer and AuNPs-mediated photo-thermal inactivation of E. coli at varying laser powers using single-phase CFD modeling

In the wake of the COVID-19 pandemics, the demand for innovative and effective methods of bacterial inactivation has become a critical area of research, providing the impetus for this study. The purpose of this research is to analyze the AuNPs-mediated photothermal inactivation of E. coli. Gold nanoparticles irradiated by laser represent a promising technique for combating bacterial infection that combines high-tech and scientific progress. The intermediate aim of the work was to present the calibration of the model with respect to the gold nanorods experiment. The purpose of this work is to study the effect of initial concentration of E. coli bacteria, the design of the chamber and the laser power on heat transfer and inactivation of E. coli bacteria. Using the CFD simulation, the work combines three main concepts. 1. The conversion of laser light to heat has been described by a combination of three distinctive approximations: a- Discrete particle integration to take into account every nanoparticle within the system, b- Rayleigh-Drude approximation to determine the scattering and extinction coefficients and c- Lambert–Beer–Bourger law to describe the decrease in laser intensity across the AuNPs. 2. The contribution of the presence of E. coli bacteria to the thermal and fluid-dynamic fields in the microdevice was modeled by single-phase approach by determining the effective thermophysical properties of the water-bacteria mixture. 3. An approach based on a temperature threshold attained at which bacteria will be inactivated, has been used to predict bacterial response to temperature increases. The comparison of the thermal fields and temporal temperature changes obtained by the CFD simulation with those obtained experimentally confirms the accuracy of the light-heat conversion model derived from the aforementioned approximations. The results show a linear relationship between maximum temperature and variation in laser power over the range studied, which is in line with previous experimental results. It was also found that the temperature inside the microchamber can exceed 55 °C only when a laser power higher than 0.8 W is used, so bacterial inactivation begins. The experimental data allows to determinate the concentration of nanoparticles. This parameter is introduced into the mathematical model obtaining the same number of AuNPs. However, this assumption introduces a certain simplification, as in the mathematical model the distribution of nanoparticles is uniform. This work is directly connected to the use of gold nanoparticles for energy conversion, as well as the field of bacterial inactivation in microfluidic systems such as lab-on-a-chip. Presented mathematical and numerical models can be extended to the entire spectrum of wavelengths with particular use of white light in the inactivation of bacteria. This work represents a significant advancement in the field, as to the best of the authors’ knowledge, it is the first to employ a single-phase computational fluid dynamics (CFD) approach specifically combined with the thermal inactivation of bacteria. Moreover, this research pioneers the use of a numerical simulation to analyze the temperature threshold of photothermal inactivation of E. coli mediated by gold nanorods (AuNRs). The integration of these methodologies offers a new perspective on optimizing bacterial inactivation techniques, making this study a valuable contribution to both computational modeling and biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hardware-Efficient Quantum Phase Estimation via Local Control

Quantum phase estimation plays a central role in quantum simulation as it enables the study of spectral properties of many-body quantum systems. Most variants of the phase estimation algorithm require the application of the global unitary evolution conditioned on the state of one or more auxiliary qubits, posing a significant challenge for current quantum devices. In this work, we present an approach to quantum phase estimation that uses only locally controlled operations, resulting in a significantly reduced circuit depth. At the heart of our approach are efficient routines to measure the complex phase of the expectation value of the time-evolution operator, the so-called Loschmidt echo, for both circuit dynamics and Hamiltonian dynamics. By tracking changes in the phase during the dynamics, the routines trade circuit depth for increased sampling cost and classical postprocessing. Our approach does not rely on reference states and is applicable to any efficiently preparable state, regardless of its correlations. We provide a comprehensive analysis of the sample complexity and illustrate the results with numerical simulations. Our methods offer a practical pathway for measuring spectral properties in large many-body quantum systems using current quantum devices.

Schiffer, Benjamin F. [Max Planck Institute of Qua↗
DOE OSTI2024

A discontinuous piecewise polynomial generalized moving least squares scheme for robust finite element analysis on arbitrary grids

A variational approach is developed with a meshless discretization to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. Here, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality, which we call discontinuous piecewise polynomial generalized moving least squares (DPP-GMLS). At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework, providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5 x improvement in accuracy and allows for taking an 8 x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Toward a quantitative description of solvation structure: a framework for differential solution scattering measurements

Appreciating that the role of the solute–solvent and other outer-sphere interactions is essential for understanding chemistry and chemical dynamics in solution, experimental approaches are needed to address the structural consequences of these interactions, complementing condensed-matter simulations and coarse-grained theories. High-energy X-ray scattering (HEXS) combined with pair distribution function analysis presents the opportunity to probe these structures directly and to develop quantitative, atomistic models of molecular systems in situ in the solution phase. However, at concentrations relevant to solution-phase chemistry, the total scattering signal is dominated by the bulk solvent, prompting researchers to adopt a differential approach to eliminate this unwanted background. Though similar approaches are well established in quantitative structural studies of macromolecules in solution by small- and wide-angle X-ray scattering (SAXS/WAXS), analogous studies in the HEXS regime—where sub-ångström spatial resolution is achieved—remain underdeveloped, in part due to the lack of a rigorous theoretical description of the experiment. To address this, herein we develop a framework for differential solution scattering experiments conducted at high energies, which includes concepts of the solvent-excluded volume introduced to describe SAXS/WAXS data, as well as concepts from the time-resolved X-ray scattering community. Our theory is supported by numerical simulations and experiment and paves the way for establishing quantitative methods to determine the atomic structures of small molecules in solution with resolution approaching that of crystallography.

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

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗
DOE OSTI2026

Mesoscale Organization in Cumulus-Coupled Stratocumulus

Marine cloud systems cover a substantial portion of the world’s oceans. Most of these clouds form relatively close to the ocean surface, typically within one to two kilometers, a region referred to by meteorologists as the marine boundary layer. They are composed predominantly of liquid water, although ice particles can occur in mid- and high-latitude marine clouds during winter. In satellite imagery, these clouds appear bright against the darker ocean surface below, reflecting a large fraction of incoming sunlight back into space that would otherwise warm the ocean. Because marine boundary layer clouds cover such an extensive area of the ocean, they exert a significant influence on Earth’s overall transfer of solar energy absorbed by the surface and thermal energy emitted to space, a balance known as the planetary radiation budget. Marine boundary layer clouds are typically thin, and their formation and dissipation depend on a delicate balance between processes acting at the ocean surface below and the warm, dry air above. They are notoriously difficult to simulate accurately in weather forecast models, which often produce too few marine low clouds in midlatitudes and clouds in tropical regions that are excessively bright, meaning they reflect too much solar radiation. The marine boundary layer is frequently characterized by widespread overcast cloud cover that often transitions from a continuous, single-layer deck to more broken cloud fields toward the tropics. These transitions typically proceed through an intermediate stage in which shallow, broken clouds form beneath the overlying stratiform cloud deck. Once broken clouds develop below the overcast, they frequently self-organize into cloud clusters known as marine boundary layer convective complexes (MBLCCs), although the mechanisms governing the formation and organization of MBLCCs remain poorly understood. Accurately representing these transitions in long-range weather forecast models is essential because they influence the properties of air masses advected over the continental United States and Europe, and they become increasingly important for forecasts on seasonal and longer timescales. We employed two complementary approaches to investigate the processes controlling MBLCCs and their impact on marine cloud cover. Long-term observations from the U.S. Department of Energy’s Eastern North Atlantic (ENA) Observatory provided a unique dataset that allowed us to characterize fundamental properties of MBLCCs, including their typical size and frequency of occurrence. These observations were combined with high-resolution numerical simulations performed on supercomputers to examine the evolution of MBLCCs during cold-air outbreaks over the ENA region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Comparison of interlaminar damage modeling strategies for hybrid composite/aluminum laminates subjected to low-velocity impact

Low-velocity impact of hybrid metal-composite structures was investigated experimentally and computationally. Composite laminates consisting of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) were joined with a 6061-T6 aluminum plate using an epoxy adhesive. Two variations of the structure were studied; one consisting of all plies oriented at 0° and one consisting of all plies oriented at 45°. A drop tower was used to impact structures at a range of energies, including energies above and below the threshold at which the aluminum layer was perforated. Numerical simulations were implemented using Sierra/SM, an in-house transient dynamics finite element code developed at Sandia National Laboratories. A Hosford plasticity model was used to describe the response of the aluminum layer. A newly implemented orthotropic continuum damage mechanics (CDM) constitutive model was used to represent the composite laminate. This 3D-CDM model was compared to a cohesive zone model (2D-CDM/CZM) to investigate efficacy of aluminum perforation energy prediction, delamination prediction, and computational cost. Accuracy of each model was evaluated using the experimental results. Each showed good agreement with the tests for both the force and velocity histories, as well as the observed damage mechanisms. The 2D-CDM/CZM model was marginally more accurate in capturing both the composite and aluminum behavior — this model averaged error percentages of -11.2% and 10.8% for residual velocity and peak force, respectively. Meanwhile, the 3D-CDM model predictions yielded average error percentages of -35.5% (velocity) and 22.6% (force). However, the 3D-CDM model generally resulted in a decreased computational cost; the average run time was 14% shorter than the 2D-CDM/CZM model and 3x as many timesteps per hour were computed using the same computational resources. In conclusion, new experimental data on the impact and perforation resistance of metal-composite laminates is presented in addition to numerical predictions of the impact behavior.

Carbon fiber↗
DOE OSTI2025

Taking control of compressible modes: bulk viscosity and the turbulent dynamo

Many polyatomic astrophysical plasmas are compressible and out of chemical and thermal equilibrium, introducing a bulk viscosity into the plasma via the internal degrees of freedom of the molecular composition, directly impacting the decay of compressible modes, $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$. This is especially important for small-scale, turbulent dynamo processes in the interstellar medium (ISM), which are known to be sensitive to the effects of compression. To control the viscous properties of $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$, we perform trans-sonic, visco-resistive dynamo simulations with additional bulk viscosity $\nu _{\text{bulk}}$, deriving a new $\nu _{\text{bulk}}$ Reynolds number $\text{Re}_{\text{bulk}}$, and viscous Prandtl number $\text{P}\nu \equiv \text{Re}_{\text{bulk}}/ \text{Re}_{\text{shear}}$, where $\text{Re}_{\text{shear}}$ is the shear viscosity Reynolds number. We derive a framework for decomposing $E_{\rm mag}$ growth rates into incompressible and compressible terms via orthogonal tensor decompositions of $\boldsymbol {\nabla }\otimes \mathrm{{\boldsymbol {\mathit {v}}}}$, where $\mathrm{{\boldsymbol {\mathit {v}}}}$ is the fluid velocity. We find that $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ play a dual role, growing and decaying $E_{\rm mag}$, and that field-line stretching is the main driver of growth, even in compressible dynamos. In the absence of $\nu _{\text{bulk}}$ ($\text{P}\nu \rightarrow \infty$), $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ pile up on small-scales, creating a spectral bottleneck, which disappears for $\text{P}\nu \approx 1$. As $\text{P}\nu$ decreases, $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ are dissipated at increasingly larger scales, in turn suppressing incompressible modes through a coupling between high-k modes. We emphasize the importance of further understanding the role of $\nu _{\text{bulk}}$ in compressible astrophysical plasmas, which we estimate could be as strong as the shear viscosity in the cold ISM, and highlight that compressible direct numerical simulations without bulk viscosity have unresolved compressible mode dissipation scales.

MHD↗
DOE OSTI2024

Numerical modeling and experimental validation of low velocity impact of woven GFRP/CFRP composites

Low-velocity impact of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) composite laminates was studied experimentally and numerically. Hybrid laminates containing blocked layers of GFRP/CFRP/GFRP with all plies oriented at 0° were investigated. Relatively high impact energies were used to obtain full perforation of the laminate in a low-velocity impact setup. Numerical simulations were carried out using the in-house transient dynamics finite element code, Sierra/SM, developed at Sandia National Laboratories. A three-dimensional continuum damage model was used to describe the response of a woven composite ply. Two methods for handling delamination were considered and compared: (1) cohesive zone modeling and (2) continuum damage mechanics. The reduced model size achieved by omission of the cohesive zone elements produced acceptable results at reduced computational cost. Further, the comparison between different modeling techniques can be used to inform modeling decisions relevant to low velocity impact scenarios. The modeling was validated by comparing with the experimental results and showed good agreement in terms of predicted damage mechanisms and impactor velocity and force histories.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Modeling and Experimental Demonstration of Flux Spreading in Light Trapping Planar-Cavity Solar Enclosed Particle Receivers

This study experimentally validates and numerically models the flux-spreading effect in a light-trapping planar-cavity solar receiver) for particle-based concentrating solar power systems. The receiver's shallow cavity with vertical planar walls redistributes concentrated solar flux, reducing peak intensity and achieving uniform heat flux. On-sun tests at National Renewable Energy Laboratory's High-Flux Solar Furnace under flux up to 1500 kW/m2 measured cavity wall temperatures, which were compared with Monte Carlo ray-tracing (SolTrace) and computational numerical simulations. Three angular absorptance models were evaluated: constant absorptance, a Pyromark-based directional model, and a Fresnel-based Cr2O3 model. The Fresnel-derived model showed the best agreement with experiments, achieving high correlation (PC > 0.85), structural similarity (SSIM > 0.98), and signal-to-noise ratios (PSNR > 40 dB), with temperature prediction errors of 1-11%. Results confirm that flux spreading mitigates local overheating and validate the integrated modeling approach, supporting the solar receiver scalability for high-efficiency, high-temperature concentrating solar power applications.

14 SOLAR ENERGY↗
DOE OSTI2026

Best practices in software development for robust and reproducible geoscientific models based on insights from the Global Carbon Budget's dynamic vegetation models

Computational models play an increasingly vital role in scientific research by enabling the numerical simulation of complex processes. Such models are also fundamental in geosciences. For instance, they offer critical insights into the impacts of global change on the Earth system today and in the future. Beyond their value as research tools, models are also software products and should therefore adhere to certain established software engineering standards. However, scientists are rarely trained as software developers, which can lead to potential deficiencies in software quality like unreadable, inefficient, or erroneous code. The complexity of models, coupled with their integration into broader workflows, also often makes it challenging to reproduce results, evaluate processes, and build upon them. In this paper, we review the state and current practices of the development processes of the state-of-the-art land surface models used by the Global Carbon Budget. We combine the experience of modelers from the respective research groups with the expertise of software engineers from tech companies to outline key principles and tools for improving software quality in research. We explore four main areas: (1) model testing and validation, (2) scientific, technical, and user documentation, (3) version control, continuous integration, and code review, and (4) the portability and reproducibility of workflows. Our review reveals that while modeling communities are incorporating many best practices, significant room for improvement remains in areas such as automated testing, automated documentation, and reproducibility. Therefore, we here identify and promote essential software engineering practices, including numerous examples of practices from within the community that can serve as guidelines for other models and could help streamline processes across the entire community. We conclude with an open-source example implementation of these principles, demonstrating portable and reproducible data flows, a continuous integration setup, and web-based visualizations. This example may serve as a practical resource for model developers, users, and all scientists engaged in scientific programming.

Gregor, Konstantin [Technical Univ. of Munich (Ger↗
DOE OSTI2024

A New Simplified Discrete Fracture Model for Shearing of Intersecting Fractures and Faults

Abstract Shearing of fractures and faults is important because it can result in permeability change or even induce seismicity—both are keys for efficient and safe energy recovery and storage in Earth systems. Quantitative analysis of shearing of intersecting fractures and faults is challenging because it can involve dynamic frictional contacts that are complicated by deformation of the rock matrix. To predict the shearing of intersecting fractures/faults, we attempt to answer the question of how intersections impact the shearing of a fracture network and whether we can simplify the description as compared to classical discrete fracture network (DFN) models. To answer these questions, we conducted a series of numerical simulations on scenarios for variable numbers of intersecting fractures. All these examples yield consistent results: the results of using DFNs are consistent with those of using hypothetical major paths. This leads to a new model, which we name simplified discrete fracture network model , to analyze shearing of intersecting fractures/faults using major path(s). We found that the intersections of fractures do not fundamentally change the shearing of two intersecting fractures if the intersecting angles are small. Furthermore, increasing the number of fractures/faults may relax the stress as more fractures/faults become available for shearing and distributing the stress. The simplified DFN model, which can capture efficiently the shearing behavior of each major paths from a large number of intersecting fractures/faults, will be a promising conceptual model that is complementary to existing equivalent continuum and discrete fracture models to analyze shearing of intersecting fractures/faults.

58 GEOSCIENCES↗
DOE OSTI2025

On the viability of stimulated hydrogen generation from iron-rich formations

Hydrogen-based technologies present a promising solution for the global energy transition. In addition to electrolytic production, subsurface geological formations provide a potential natural source of hydrogen. Iron-rich ultramafic rocks, in particular, are favorable for hydrogen generation through natural processes such as serpentinization. Naturally occurring reactions and migration can be enhanced through various types of stimulation, including thermal, hydraulic, and chemical treatment. Through numerical simulations, we analyzed the complex interplay of factors influencing the production and migration within the subsurface, emphasizing the importance of different stimulation techniques, catalysts, and conditions. Our findings indicate that key parameters, such as damage zone permeability and width, significantly impact producible hydrogen mass. Our results indicate that a combination of large damage zone widths, high permeability, and a stimulated reaction rate of 1 × 10 -9 can yield economically viable production rates of up to 1 kg s -1 at the wellhead. Moreover, the availability of ferrous iron, rather than the serpentinization rate itself, has been identified as the primary limiting factor in achieving economically sustainable hydrogen production. In conclusion, while an unstimulated rock volume of 0.165 km 3 yields only 45t of hydrogen in two years, various stimulation techniques can increase production to 18500t.

08 - HYDROGEN↗
DOE OSTI2026

Oscillatory and Collective Dynamics of Gold‐Nanoparticle‐Laden Droplets Driven by Photothermal‐Induced Thermocapillarity

Droplets have long intrigued researchers due to their ability to exhibit complex and fascinating behavior when subjected to external stimuli. Here, a coupled oscillatory behavior of gold-nanoparticle-surfactant-laden aqueous droplets is investigated at an oil-oil interface stimulated by light. This study shows that the interaction between light and the droplets gives rise to a range of oscillatory modes, including bouncing and clustering, where droplets exhibit collective movement. From experiment and numerical simulations, this study elucidates the underlying mechanism: upon laser irradiation, gold nanoparticles convert light into heat, generating asymmetric thermal gradients that drive upward thermocapillary flows and a hydrodynamic force from photothermal convection. These forces compete with gravity and buoyancy to induce droplet bouncing, while the resulting asymmetric flow fields bias neighboring droplets toward the illuminated droplet, leading to clustering. These findings not only expand the library of dynamic droplet behavior but also offer insights into the potential applications of light-driven systems in materials science, soft robotics, and microfluidics.

Marangoni effect↗
DOE OSTI2026

Terahertz Chiral Metamaterials Enabling Broadband Polarization Conversion Using Polarization Guiding Effect

Polarization control plays a vital role in terahertz (THz) photonics, enabling a wide range of applications from imaging and spectroscopy to sensing and wireless communication. However, conventional polarization control methods at THz frequencies are limited by narrow operational bandwidths and excessive absorption losses. In this work, we present a broadband THz polarization control based on the polarization guiding effect. We design and fabricate a THz polarization rotator capable of rotating incident linear polarization by 90°. The device consists of 24 twistingly stacked silicon–air metagrating layers that introduce a large degree of form birefringence, with each metagrating layer being 280 µm thick and sequentially twisted by 3.75°. Numerical simulations using the Berreman 4 × 4 method predict broadband operation of the device. Experimental validation using THz time-domain spectroscopy confirms efficient and broadband 90° polarization rotation from 0.2 to 1.25 THz. Based on this concept, we further propose THz Q-plates for the generation of broadband structured vector beams with radial and azimuthal polarization states. This work offers a scalable, material-agnostic platform for advanced THz polarization manipulation and vector beam engineering, opening new pathways for next-generation THz photonic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Laterally Modulating Carrier Concentration by Ion Irradiation in CdO Thin Films for Mid‐IR Plasmonics

This report demonstrates tunable carrier densities in CdO thin films through local ion irradiation, providing lateral control of mid‐IR optical properties. Ion‐solid interactions produce donor‐like defects that boost electron concentrations from the practical minimum of 2.5 × 10 19 cm −3 to a maximum of 2.5 × 10 20 cm −3 by metered ion exposure. This range is achieved using He, N, Ar, or Au ions at 1–2.8 MeV; when normalized by displacements per atom, all ion species produce comparable results. Since CdO is well‐described by the Drude model, irradiation‐tuned carrier densities directly alter the infrared dielectric function, and in turn, mid‐infrared optical properties. Further, it is demonstrated that by combining irradiation with traditional lithography, CdO films expose to ions in the presence of 3‐µm thick, patterned photoresist exhibit lateral carrier density profiles with ≈400‐nm resolution. Scanning near‐field optical microscopy reveals sharp optical interfaces with almost no companion contrast in surface morphology, microstructure, or crystallinity. Finally, CdO lateral homostructures supporting surface plasmon polaritons (SPPs) are demonstrated whose dispersion relation can be tuned through periodic patterning in a monolithic platform by simple nanofabrication. Numerical simulations show these polaritons result from strong coupling between excitations at CdO plasma frequencies and SPPs supported by the platinum substrate.

cadmium oxide↗
DOE OSTI2024

Dynamics of Intra-Cell Thermal Front Propagation in Lithium-Ion Battery Safety Issues

Thermal runaway (TR), a critical failure mode in lithium-ion batteries (LIBs), poses significant safety risks and hinders wider application of LIBs. TR typically begins at a localized heat source and spreads across the cell. Understanding thermal front propagation (TFP) characteristics, such as front and velocity, is crucial for assessing energy release and temperature distribution for battery hazardous estimation. Recent studies assume that TR within cells propagates at a near-constant velocity, based on the reaction kinetics and thermal properties. Here, in this study, an intra-battery TR model is further proposed and it indicates that TFP velocity stabilizes when the front is distanced from the heat source. Theoretical estimates for propagation velocity and front are developed and validated through numerical simulations and experimental tests from the NREL Battery Failure Databank. The energy release rate during TFP and the impact of preheating based on a point heat source are explored. This work clarifies the long-standing clouds of the thermal font propagation behaviors within the single cell, highlights the power and beauty of mathematics modeling to describe the complicated thermal behaviors, and provides important guidelines for thermal hazardous understanding for next-generation batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

Anthropogenic effects on flood hazards in a hyper-arid watershed: The 2015 Atacama floods

An unprecedented precipitation event in the hyper-arid Atacama Desert of Northern Chile occurred in March 2015. Geomorphic alterations to the river channel and the coastal zone, coupled with the exceptional magnitude of the rainfall, caused catastrophic damage and loss of life. On the coast of the El Salado watershed, legacy mine tailings infilled the watershed-ocean connection, while the river channel was altered both by tailings and urbanization. The consequences of this event resulted from the coupling of anthropogenic geomorphic changes with an unusual climate event. Lack of field data, complex geomorphology and sediment loads influenced by human activity make analysing floods in these regions especially challenging. The objective of this work is to improve our understanding of the factors that control flood hazards by using numerical simulations to reconstruct the 2015 flood in El Salado. We carry out unsteady two-dimensional simulations fully coupled with the sediment concentration to identify the influence of tailing deposits, considering high-resolution data of the pre- and post-2015 flood topography. In conclusion, the results highlight the importance of specific event-based studies, using models that can help designing better strategies for climate change adaptation and risk mitigation, while providing information for risk reduction and channel restoration.

54 ENVIRONMENTAL SCIENCES↗