Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence analysis”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An Evaluation of Tropical Cyclone Rainfall Structures in the HighResMIP Simulations against Satellite Observations

This study evaluates tropical cyclone (TC) rainfall structures in the CMIP6 HighResMIP global climate model (GCM) simulations against satellite rainfall retrievals. Here, we specifically focus on TCs within the deep tropics (25°S–25°N). Analysis of TC rain rate composites indicates that in comparison to the satellite observations at the same intensity, many HighResMIP simulations tend to overproduce rain rates around TCs, in terms of both maximum rain rate magnitude and area-averaged rain rates. In addition, as model horizontal resolution increases, the magnitude of the peak rain rate appears to increase. However, the area-averaged rain rates decrease with increasing horizontal resolution, partly due to the TC eyewall being located closer to the TC center, thus occupying a smaller area and contributing less to the area-averaged rain rates. The effect of ocean coupling is to lower the TC rain rates, bringing them closer to the satellite observations, due to reduced horizontal moisture flux convergence and surface latent heat flux beneath TCs. Examination of horizontal rain rate distributions indicates that vertical wind shear–induced rainfall asymmetries in HighResMIP-simulated TCs are qualitatively consistent with the observations. In addition, a positive relationship is observed between the area-averaged inner-core rainfall and TC intensification likelihoods across the HighResMIP simulations, as GCM simulations producing stronger TCs more frequently have the greater rainfall close to the center, in agreement with previous theoretical and GCM simulation results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

What Is the Limit of Quantification for the Minor Phase in Time-of-Flight Neutron Diffraction? A Case Study on Fe and Ni Powder Mixtures at VULCAN

A phase present in small quantities within materials may not simply serve as a secondary component; it can play a crucial role in determining the integrity, properties, and performance of the material. These minor but important phases usually draw attention in material design and processing for fundamental understanding as well as material quality control. Accurately quantifying a minor phase amid a majority phase, especially at extremely low fractions, remains a challenging task. Time-of-flight neutron diffraction, coupled with advanced pattern analysis techniques like Rietveld refinement, is a powerful tool for crystal structure identification and phase quantification. The deep penetrating capability of neutrons enables the detection and quantification of trace phases within materials. In this study, the quantification limits of time-of-flight neutron diffraction were explored using the VULCAN diffractometer at the Spallation Neutron Source, using Fe–Ni powder mixtures as a sample system. By comparing the refinement results to the known weighed values, it was determined that the reliable quantification of a minor Ni phase is achievable down to about 0.1 wt% while a Ni fraction as low as 0.02 wt% is difficult to trace. Effective control of the refinement parameters, especially the profile function parameters, are found to significantly influence the convergence of fittings and the accuracy of phase quantification.

Rietveld refinement↗
DOE OSTI2020

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Phenomena-based graph representations and applications to chemical process simulation

Rapid and robust simulation of chemical processes is critical to conduct process design, optimization, techno-economic analysis, and sustainability analysis. Yet, efficiently solving simulation models remains a challenge due to the highly coupled and nonlinear nature of the underlying algebraic equations that capture the physical phenomena taking place in the process (e.g., material and energy conservation, phase equilibrium, reactions). In this work, we show that graph-theoretic representations of the physical phenomena within unit operations can help navigate and decompose equations to systematically identify alternative approaches for fast and robust numerical solutions. Specifically, we present a graph-theoretic abstraction that captures the connectivity between the model variables/equations and use this abstraction to group variables/equations into fundamental phenomena. We show that phenomena-based decomposition of the underlying equations can help decouple nonlinearities and enforce material/energy conservation at the process level to accelerate convergence. The proposed decomposition approach differs from the more traditional sequential modular simulation approach, in which equations are grouped and decomposed by unit operations. We implemented the phenomena-based decomposition in BioSTEAM—an open-source process simulation platform in Python—and demonstrated that this approach can converge a variety of separation process models. Compared to sequential modular simulation, the phenomena-based approach can converge idealized systems faster, but it can be slower for (or even fail to converge) highly coupled and nonideal process systems.

Convergence↗
DOE OSTI2022

Adaptive elasticity policies for staging-based in situ visualization

In situ processing aims to alleviate the growing gap between computation and I/O capabilities by performing data processing close to the data source. In situ processing is widely used to process data generated by multiple data sources, including observation data from edge devices or scientific observational facilities and the simulation data generated by scientific computation on a high-performance computing (HPC) platform. For a scientific workflow that is run on an HPC platform and composed of a simulation program and an in situ data analytics or visualization (abbreviated as ana/vis) task, there is an implicit assumption that the computing resources assigned to the workflow keep static during the workflow execution. However, with the converging trend between the HPC and cloud computing platform, running the in situ ana/vis task in an elastic way is promising to decrease its overhead and improve its resource utilization rate. Resource elasticity represents the ability to change resource configurations such as the number of computing nodes/processes during workflow execution. An elastic job may dynamically adjust resource configurations; it may use a few resources at the beginning and more resources toward the end of the job when interesting data appear. However, it is hard to predict a priori how many computing nodes/processes need to be added/removed during the workflow execution to adapt to changing workflow needs. How to efficiently guide elasticity operations, such as growing or shrinking the number of processes used for in situ analysis during workflow execution, is an open-ended research question. In this article, we present adaptive elasticity policies that adopt workflow runtime information collected during workflow execution to predict how to trigger the addition/removal of processes in order to minimize in situ processing overhead. Taking in situ visualization tasks as an example, we integrate the presented elasticity policies into a staging-based elastic workflow and evaluate its efficiency in multiple elasticity scenarios. Compared with the situation without elasticity or with a static elasticity policy that uses a fixed number of processes for each rescaling operation, the adaptive elasticity policy can save overhead in finding a proper resource configuration and improve resource utilization efficiency. Furthermore, one experiment illustrates that the adaptive elasticity policy saves 41% of core-hours compared with the situation without the resource elasticity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Proper orthogonal descriptors for efficient and accurate interatomic potentials

Here, we present the proper orthogonal descriptors for efficient and accuracy representation of the potential energy surface. The potential energy surface is represented as a many-body expansion of parametrized potentials in which the potentials are functions of atom positions and parameters. The proper orthogonal decomposition is employed to decompose the parametrized potentials into a set of proper orthogonal descriptors (PODs). Because of the rapid convergence of the proper orthogonal decomposition, relevant snapshots can be sampled exhaustively to represent the atomic neighborhood environment accurately with a small number of descriptors. The proper orthogonal descriptors are used to develop interatomic potentials by using a linear expansion of the descriptors and determining the expansion coefficients from a weighted least-squares regression against a density functional theory (DFT) training set. We present a comprehensive evaluation of the POD potentials on previously published DFT data sets comprising Li, Mo, Cu, Ni, Si, Ge, and Ta elements. The data sets represent a diverse pool of metals, transition metals, and semiconductors. The accuracy of the POD potentials are comparable to that of state-of-the-art machine learning potentials such as the spectral neighbor analysis potential (SNAP) and the atomic cluster expansion (ACE).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Glancing Angle Deposition in Gas Sensing: Bridging Morphological Innovations and Sensor Performances

Glancing Angle Deposition (GLAD) has emerged as a versatile and powerful nanofabrication technique for developing next-generation gas sensors by enabling precise control over nanostructure geometry, porosity, and material composition. Through dynamic substrate tilting and rotation, GLAD facilitates the fabrication of highly porous, anisotropic nanostructures, such as aligned, tilted, zigzag, helical, and multilayered nanorods, with tunable surface area and diffusion pathways optimized for gas detection. This review provides a comprehensive synthesis of recent advances in GLAD-based gas sensor design, focusing on how structural engineering and material integration converge to enhance sensor performance. Key materials strategies include the construction of heterojunctions and core–shell architectures, controlled doping, and nanoparticle decoration using noble metals or metal oxides to amplify charge transfer, catalytic activity, and redox responsiveness. GLAD-fabricated nanostructures have been effectively deployed across multiple gas sensing modalities, including resistive, capacitive, piezoelectric, and optical platforms, where their high aspect ratios, tailored porosity, and defect-rich surfaces facilitate enhanced gas adsorption kinetics and efficient signal transduction. These devices exhibit high sensitivity and selectivity toward a range of analytes, including NO2, CO, H2S, and volatile organic compounds (VOCs), with detection limits often reaching the parts-per-billion level. Emerging innovations, such as photo-assisted sensing and integration with artificial intelligence for data analysis and pattern recognition, further extend the capabilities of GLAD-based systems for multifunctional, real-time, and adaptive sensing. Finally, current challenges and future research directions are discussed, emphasizing the promise of GLAD as a scalable platform for next-generation gas sensing technologies.

Chemistry↗
DOE OSTI2024

Theoretical analysis of the OH-initiated atmospheric oxidation reactions of imidazole

Imidazoles are present in Earth's atmosphere in both the gas-phase and in aerosol particles, and have been implicated in the formation of brown carbon aerosols. The gas-phase oxidation of imidazole (C 3 N 2 H 4 ) by hydroxyl radicals has been shown to be preferentially initiated via OH-addition to position C5, producing the 5-hydroxyimidazolyl radical adduct. However, the fate of this adduct upon reaction with O 2 in the atmospheric gas-phase is currently unknown. We employed an automated approach to investigate the reaction mechanism and kinetics of imidazole's OH-initiated gas-phase oxidation, in the presence of O 2 and NO x . The explored mechanism included reactions available to first-generation RO 2 radicals, as well as alkoxyl radicals produced from RO 2 + NO reactions. Product distributions were obtained by assembling and solving a master equation, under conditions relevant to the Earth's atmosphere. Our calculations show a complex, branched reaction mechanism, which nevertheless converges to yield two major closed-shell products: 4H-imidazol-4-ol (4H-4ol) and N,N′-diformylformamidine (FMF). At 298 K and 1 atm, we estimate the yields of 4H-4ol and FMF from imidazole oxidation initiated via OH-addition to position C5 to be 34 : 66, 12 : 85 and 2 : 95 under 10 ppt, 100 ppt and 1 ppb of NO respectively. This work also revealed O 2 -migration pathways between the α-N-imino peroxyl radical isomers. This reaction channel is fast for the first-generation RO 2 radicals, and may be important during the atmospheric oxidation of other unsaturated organic nitrogen compounds as well.

Almeida, Thomas Golin↗
DOE OSTI2024

Effects of Injected Current Streams on MHD Equilibrium Reconstruction of Local Helicity Injection Plasmas in a Spherical Tokamak

Open field line currents are intrinsic to DC helicity injection plasma startup and pose a challenge for inferring the plasma equilibrium with standard reconstruction analysis. Local helicity injection (LHI) is a type of DC helicity injection which uses small, modular current sources to drive force-free current along helical field lines to produce tokamak plasmas. MHD modeling and magnetic measurements during LHI indicate the injected current streams remain coherent as helical structures on the outboard edge of a core toroidal plasma that is tokamak-like in a toroidally averaged sense. To extract core plasma equilibrium properties, external magnetic diagnostics corrected for contributions from the injected current streams are fitted by a standard Grad-Shafranov equilibrium code. An iterative approach for estimating and subtracting the stream contributions from the diagnostic signals is described and applied to a model equilibrium database to reduce systematic errors introduced by the streams. Convergence is usually attained with 2 to 4 iterations, with derived equilibrium parameters matching the prescribed axisymmetric core values to within estimated experimental uncertainties. Accurate recovery of core parameters occurs when the ratio of the net toroidal windup current from the streams to the core plasma current is less than 0.2, which is typically satisfied in most experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effects of Substance Use and Antisocial Personality on Neuroimaging-Based Machine Learning Prediction of Schizophrenia

Abstract Background and hypothesis Neuroimaging-based machine learning (ML) algorithms have the potential to aid the clinical diagnosis of schizophrenia. However, literature on the effect of prevalent comorbidities such as substance use disorder (SUD) and antisocial personality (ASPD) on these models’ performance has remained unexplored. We investigated whether the presence of SUD or ASPD affects the performance of neuroimaging-based ML models trained to discern patients with schizophrenia (SCH) from controls. Study design We trained an ML model on structural MRI data from public datasets to distinguish between SCH and controls (SCH = 347, controls = 341). We then investigated the model’s performance in two independent samples of individuals undergoing forensic psychiatric examination: sample 1 was used for sensitivity analysis to discern ASPD (N = 52) from SCH (N = 66), and sample 2 was used for specificity analysis to discern ASPD (N = 26) from controls (N = 25). Both samples included individuals with SUD. Study results In sample 1, 94.4% of SCH with comorbid ASPD and SUD were classified as SCH, followed by patients with SCH + SUD (78.8% classified as SCH) and patients with SCH (60.0% classified as SCH). The model failed to discern SCH without comorbidities from ASPD + SUD (AUC = 0.562, 95%CI = 0.400–0.723). In sample 2, the model’s specificity to predict controls was 84.0%. In both samples, about half of the ASPD + SUD were misclassified as SCH. Data-driven functional characterization revealed associations between the classification as SCH and cognition-related brain regions. Conclusion Altogether, ASPD and SUD appear to have effects on ML prediction performance, which potentially results from converging cognition-related brain abnormalities between SCH, ASPD, and SUD.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Overlapping Domain Coupling of Multidimensional and System Codes in NEAMS - Pronghorn and SAM

This report describes the implementation and testing of domain-overlapping coupling of the Pronghorn and SAM codes. Both Pronghorn and SAM are codes developed by the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program for analysis of advanced nuclear reactors. Pronghorn focuses on analyzing multi-dimensional core flow conditions, while SAM focuses on analysis of the entire system including, among other components, piping, pumps, and heat exchangers. It is desirable for many advanced reactor thermal-hydraulics simulations to couple Pronghorn and SAM to obtain self-consistent solutions in the core and the system. The domain-overlapping coupling approach provides a robust and numerically efficient candidate for this coupling. We describe the algorithm for obtaining self-consistent solutions between a SAM network of flow channels and an overlapped multi-dimensional Pronghorn domain. The algorithm is tested for a variety of simple canonical test problems and a model of the Molten Salt Reactor Experiment (MSRE). Results demonstrate that the domain overlapping coupling approach provides consistent pressure, mass flow rate, enthalpy, and passive scalar distributions between Pronghorn and SAM and converges reliably and efficiently (at most 15 iterations, usually much less than 10) for the considered test problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

As a Matter of Dynamical Range – Scale Dependent Energy Dynamics in MHD Turbulence

Abstract Magnetized turbulence is ubiquitous in many astrophysical and terrestrial plasmas but no universal theory exists. Even the detailed energy dynamics in magnetohydrodynamic (MHD) turbulence are still not well understood. We present a suite of subsonic, super-Alfvénic, high plasma beta MHD turbulence simulations that only vary in their dynamical range, i.e., in their separation between the large-scale forcing and dissipation scales, and their dissipation mechanism (implicit large eddy simulation, ILES, and direct numerical simulation (DNS)). Using an energy transfer analysis framework we calculate the effective numerical viscosities and resistivities, and demonstrate that all ILES calculations of MHD turbulence are resolved and correspond to an equivalent visco-resistive MHD turbulence calculation. Increasing the number of grid points used in an ILES corresponds to lowering the dissipation coefficients, i.e., larger (kinetic and magnetic) Reynolds numbers for a constant forcing scale. Independently, we use this same framework to demonstrate that—contrary to hydrodynamic turbulence—the cross-scale energy fluxes are not constant in MHD turbulence. This applies both to different mediators (such as cascade processes or magnetic tension) for a given dynamical range as well as to a dependence on the dynamical range itself, which determines the physical properties of the flow. We do not observe any indication of convergence even at the highest resolution (largest Reynolds numbers) simulation at 2048 3 cells, calling into question whether an asymptotic regime in MHD turbulence exists, and, if so, what it looks like.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

pyEF: A Python Framework for QM and QM/MM Atom-Wise Electric Field Analysis

We introduce pyEF, a software package for computing molecular electric fields, electrostatic interaction energies, and electrostatic potentials from quantum mechanical (QM) atom-centered multipole expansions with atom-wise decomposable contributions. We demonstrate the computational efficiency and accuracy of this QM-derived electric field evaluation tool through several tests. To assess the influence of the underlying QM method and charge partitioning scheme on these electrostatic quantities, we analyze over 250 configurations of an acetone solute molecule in five solvents of variable polarity. We find that electric field calculations are highly sensitive to the choice of charge partitioning method. Even among real-space charge schemes, acetone Stark tuning rates differ by up to a factor of 2. Benchmarking computed solvent dipole moments against experimental bulk values, we conclude that the CM5, ADCH, and Hirshfeld-I charge schemes most reliably capture solvent electrostatics and therefore provide a more faithful foundation for computing electric fields. When constructed from these real-space charges, electric fields are nearly insensitive to basis set size and monotonically increase in magnitude with higher Fock exchange. We also demonstrate efficient convergence of QM electrostatics when more distant molecules are represented solely by MM point charges, reducing computational overhead. Leveraging these findings, we demonstrate the use of pyEF to deduce environmental effects on a transition metal complex from a Ga 4 L 6 12– nanocage and quantify the dominant role of organic linkers in orchestrating electrostatic preorganization.

electric fields↗
DOE OSTI2024

A convergent genus expansion for the plateau

We conjecture a formula for the spectral form factor of a double-scaled matrix integral in the limit of large time, large density of states, and fixed temperature. The formula has a genus expansion with a nonzero radius of convergence. To understand the origin of this series, we compare to the semiclassical theory of “encounters” in periodic orbits. In Jackiw-Teitelboim (JT) gravity, encounters correspond to portions of the moduli space integral that mutually cancel (in the orientable case) but individually grow at low energies. At genus one we show how the full moduli space integral resolves the low energy region and gives a finite nonzero answer.

2D Gravity↗
DOE OSTI2020

Stereoscopic high-speed microscopy to understand transient internal flow processes in high-pressure nozzles

The flow and cavitation behavior inside fuel injectors is known to affect spray development, mixing and combustion characteristics. We discuss that while diesel fuel injectors with converging and hydro-eroded holes are generally known to limit cavitation and feature higher discharge coefficients during the steady period of injection, less is known about the flow during transient periods corresponding to needle opening and closing. Multiple injection strategies involve short injections, multiplying these aspects and giving them a growing importance as part of the fuel delivery process. In this study, single-hole transparent nozzles were manufactured with the same hole inlet radius and diameter as the Engine Combustion Network Spray D nozzle, mounted to a modified version of a common-rail Spray A injector body and needle. Needle opening and closing periods were visualized with stereoscopic high-speed microscopy at injection pressures relevant to modern diesel engines. Time-resolved sac pressure was extracted via elastic deformation analysis of the transparent nozzles. Sources of cavitation were observed and tracked, enabling the identification of a gas exchange process after the end of injection with ingestion of chamber gas into the sac and orifice. We observed that the gas exchange contributed widely to disrupting the start of injection and outlet flow during the subsequent injection event.

42 ENGINEERING↗
DOE OSTI2021

Kinetic equilibrium reconstructions of plasmas in the MAST database and preparation for reconstruction of the first plasmas in MAST upgrade

Reconstructions of plasma equilibria using magnetic sensors were routine during operation of the Mega Ampere Spherical Tokamak (MAST) device, but reconstructions using kinetic profiles were not. These are necessary for stability and disruption analysis of the MAST database, as well as for operation in the upgrade to the device, MAST-U. The three-dimensional (3D) code VALEN is used to determine eddy currents in the 3D vessel structures for vacuum coil test shots, which are then mapped to effective resistances in the two-dimensional vessel groupings in the EFIT equilibrium reconstruction code to be used in conjunction with nearby loop voltage measurements for estimated currents in the structures during reconstruction. Kinetic equilibrium reconstructions with EFIT, using all available magnetic sensors as well as Thomson scattering measurements of electron temperature and density, charge exchange recombination spectroscopy measurements of ion temperature, and internal magnetic field pitch angle measurements from a motional Stark effect (MSE) diagnostic are performed for a large database of MAST discharges. Excellent convergence errors are obtained for the portions of the discharges where the stored energy was not too low, and it is found that reconstructions performed with temperature and density measurements but without MSE data usually already match the pitch angle measurements well. A database of 275 kinetic equilibria is used to test the ideal MHD stability calculation capability for MAST. In conclusion, the necessary changes to conducting structure in VALEN, and diagnostic setup in EFIT have been completed for the upgrade from MAST to MAST-U, enabling kinetic reconstructions to commence from the first plasma discharges of the upgraded device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Air Blast Mesh Sensitivity and Pressure Mapping Study

Nose cone structural and aerodynamic characteristics are essential for intelligent design of aircraft, spacecraft, and ballistic systems. Finite element analysis can be used to help understand the structural integrity and flight characteristics of different nose cones. A mesh sensitivity study was undertaken for a particular nose cone geometry that was used in tests at LANL facilities in order to confirm the integrity of the meshed geometry. A simple cone that best matched closed-form theoretical solutions was modeled, and received good correlation to the theory. Complexity was then added back to the nose cone. Parameters applied to the simple cone were then implemented in the nose cone geometry giving assurance of accuracy after the geometry was changed. Nose cone results averaged 6.3% error for radial displacement when compared with the theoretical. Hoop stress averaged 6.0% error and meridional stress averaged 5.7% error at the finest mesh level. Meshes showed signs of convergence when compared to all three theoretical solutions. Finally, pressure time-history data from LANL computational fluid dynamics simulations was applied to the surface of the final nose cone geometry. The pressure data was interpolated from pressure gauge locations onto nearby meshed elements, which allowed for FEA software to run simulations on the cone with the pressure data as a loading condition. The pressure mapping resulted in the ability to understand the nose cone’s rigid body motion that in turn can inform design of future nose cones.

42 ENGINEERING↗
DOE OSTI2025

Using ARM Observations to Evaluate Process-Interactions in MCS Simulations Across Scales (Final Progress Report)

This project, funded by DOE Atmospheric System Research (DE-SC0020050), focused on improving the representation of mesoscale convective systems (MCSs) in numerical weather and climate models by leveraging high-resolution observations from the DOE Atmospheric Radiation Measurement (ARM) program. The research aimed to evaluate model sensitivities to grid spacing, microphysics, and planetary boundary layer (PBL) schemes, with a particular emphasis on improving convection parameterization for high-resolution modeling. Findings from this work highlight several key advancements. Model validation against ARM radar wind profiler data from the Southern Great Plains (SGP) and Manaus (MAO) sites revealed systematic biases in simulated convective mass flux profiles, leading to the development of an observationally constrained evaluation framework for diagnosing and improving model performance. Sensitivity analyses demonstrated that the representation of Amazonian MCSs was highly dependent on PBL scheme selection, while mid-latitude MCSs were more strongly influenced by microphysics parameterizations. A series of high-resolution WRF simulations, ranging from 4 km to 125 m grid spacing, provided insight into the behavior of convective drafts across scales. While updraft properties converged at sub-kilometer resolutions, biases in downdraft intensity persisted even at the finest resolution tested, emphasizing the need for further refinements in model physics. Additionally, comparisons of MCS vertical structures between mid-latitude and tropical environments revealed stronger updrafts and larger mass flux in mid-latitude MCSs, providing critical insights for improving climate model representations of storm-scale dynamics. The project’s findings have already contributed to advancing numerical modeling capabilities, particularly in WRF, MPAS, ICON, and DOE’s SCREAM model, by refining how convective processes are represented in high-resolution climate simulations. Results were disseminated through peer-reviewed publications, conference presentations, and ARM/ASR Research Highlights, engaging the broader scientific community. The project also provided valuable training opportunities for two postdoctoral researchers, who played central roles in model development, analysis, and dissemination of results. Their work contributed to several publications and conference presentations, helping prepare them for careers in atmospheric modeling. By improving the simulation of MCSs, this research directly supports the development of next-generation climate models capable of more accurately representing extreme precipitation and convective processes. The insights gained will inform future improvements in convective parameterization and guide the design of high-resolution weather and climate simulations, ultimately enhancing the reliability of climate projections and weather forecasts.

54 ENVIRONMENTAL SCIENCES↗