Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Scaling Studies”

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 883 records · Page 49

Simulations of flaming combustion and flaming-to-smoldering transition in wildland fire spread at flame scale

Our objective in the present study is to provide basic insights into the coupling between external-gas and solid biomass vegetation processes that control the dynamics of flame spread in wildland fire problems. Here we focus on a modeling approach that resolves processes occurring at vegetation and flame scales, i.e., the formation of flammable vapors due to the thermal degradation of the solid biomass, the subsequent combustion in ambient air, the thermal feedback to the biomass through radiative and convective heat transfer, and the possible transition from flaming combustion (taking place outside of the solid biomass) to smoldering combustion (taking place inside the solid biomass). The capability uses a multiphase combustion framework and treats external-gas processes through a Large Eddy Simulation solver and solid biomass processes through a discrete particle model. The discrete particle model adopts a one-dimensional porous medium formulation, includes descriptions of drying, thermal pyrolysis, oxidative pyrolysis, and char oxidation, as well as a description of the external-gas-to-solid-biomass diffusion of oxygen mass; the discrete particle model thereby provides a treatment of in-depth oxidative processes and allows the simulation of smoldering combustion. The modeling capability is applied to the simulation of fire spread across a surrogate biomass vegetation bed corresponding to a discrete array of cylindrical-shaped, vertically-oriented, pine wood sticks, characterized by a monomodal size distribution, in horizontal flat terrain and under wind-aided conditions. The numerical results demonstrate that the model can simulate successful flaming-to-smoldering transition followed by complete biomass consumption.

01 COAL, LIGNITE, AND PEAT↗

Quantitative principles of microbial metabolism shared across scales

Metabolism is the complex network of chemical reactions occurring within every cell and organism, maintaining life, mediating ecosystem processes and affecting Earth’s climate. Experiments and models of microbial metabolism often focus on one specific scale, overlooking the connectivity between molecules, cells and ecosystems. Here we highlight quantitative metabolic principles that exhibit commonalities across scales, which we argue could help to achieve an integrated perspective on microbial life. Mass, electron and energy balance provide quantitative constraints on their flow within metabolic networks, organisms and ecosystems, shaping how each responds to its environment. The mechanisms underlying these flows, such as enzyme–substrate interactions, often involve encounter and handling stages that are represented by equations similar to those for cells and resources, or predators and prey. Here we propose that these formal similarities reflect shared principles and discuss how their investigation through experiments and models may contribute to a common language for studying microbial metabolism across scales.

59 BASIC BIOLOGICAL SCIENCES↗

Wedge-Shaped Hopper Design for Milled Woody Biomass Flow

The bioenergy industry has been challenged by unstable flow and transport of milled biomass in material handling operations. Handling issues such as hopper clogging and auger jamming are attributed to knowledge gaps between existing handling units designed for bulk solids and their suitability for milled biomass with high compressibility. This work investigates various flow behaviors of granular woody biomass in wedge-shaped hoppers. Herein, hopper flow physical experiments and numerical simulations are conducted to study the influence of the critical material attributes and critical processing parameters on the flow pattern, arching, and throughput. The results show that (1) the preferred flow pattern, mass flow, can be achieved by controlling the material’s internal friction angle, hopper inclination, and hopper wall friction; (2) hopper arching, governed by the competing gravity-driven force against flow resistance from material internal friction and material–wall friction, can be controlled by the hopper wall friction angle and the inclination angle; and (3) flow throughput can be accurately estimated from our empirical equation with inputs of hopper outlet geometry and particle-scale to bulk-scale material attributes. This study elucidates woody biomass flow physics and provides guidance for industrial equipment design.

09 BIOMASS FUELS↗

Relativistically transparent magnetic filaments: scaling laws, initial results and prospects for strong-field QED studies

Abstract Relativistic transparency enables volumetric laser interaction with overdense plasmas and direct laser acceleration of electrons to relativistic velocities. The dense electron current generates a magnetic filament with field strength of the order of the laser amplitude (>10 5 T). The magnetic filament traps the electrons radially, enabling efficient acceleration and conversion of laser energy into MeV photons by electron oscillations in the filament. The use of microstructured targets stabilizes the hosing instabilities associated with relativistically transparent interactions, resulting in robust and repeatable production of this phenomenon. Analytical scaling laws are derived to describe the radiated photon spectrum and energy from the magnetic filament phenomenon in terms of the laser intensity, focal radius, pulse duration, and the plasma density. These scaling laws are compared to 3D particle-in-cell (PIC) simulations, demonstrating agreement over two regimes of focal radius. Preliminary experiments to study this phenomenon at moderate intensity ( a 0 ∼ 30) were performed on the Texas Petawatt Laser. Experimental signatures of the magnetic filament phenomenon are observed in the electron and photon spectra recorded in a subset of these experiments that is consistent with the experimental design, analytical scaling and 3D PIC simulations. Implications for future experimental campaigns are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Comparative Analysis of Cosmological Simulations via Angular Power Spectra

Cosmology is the study of the universe's origin, large-scale structure, and dynamics. One such area of study is understanding dark energy and dark matter. Dark matter has not been observed to interact with luminous matter other than through gravitational attractions, forming large scale structures observed throughout the universe. In this area, the primary model is the Lambda Cold Dark Matter (ΛCDM) model. This standard model consists of cosmological dark energy (Λ), cold dark matter (CDM), and normal luminous matter. We use angular power spectra (APS) to compare simulated ΛCDM models with observed large-scale galaxy distributions from the DECam Local Volume Exploration Survey (DELVE) [1]. After ensuring that the simulations can produce APS that match observational results, we vary cosmological parameters to constrain the allowed ranges of parameter values that are consistent with the observed data.

79 ASTRONOMY AND ASTROPHYSICS↗

Comparative analysis of cosmological simulations via angular power spectra

Cosmology is the study of the universe's origin, large-scale structure, and dynamics. One such area of study is understanding dark energy and dark matter. Dark matter has not been observed to interact with luminous matter other than through gravitational attractions, forming large scale structures observed throughout the universe. In this area, the primary model is the Lambda Cold Dark Matter (ΛCDM) model. This standard model consists of cosmological dark energy (Λ), cold dark matter (CDM), and normal luminous matter. We use Angular Power Spectra (APS) to compare simulated ΛCDM models with observed large-scale galaxy distributions from the DECam Local Volume Exploration Survey (DELVE). After ensuring that the simulations can produce APS which match observational results, we vary cosmological parameters to constrain the allowed ranges of parameter values that are consistent with the observed data.

Feithen, William↗

Formulation and solution approach for calibrating activity-based travel demand model-system via microsimulation

This study addresses the problem of calibrating utility-maximizing nested logit activity-based travel demand model-systems. After estimation, it is common practice to use aggregate measurements to calibrate the estimated model-system’s parameters prior to their application in transportation planning, policy making, and operations. However, calibration of activity-based model-systems has received much less attention. Existing calibration approaches are myopic heuristics in the sense that they do not consider the fundamental inter-dependencies among choice-models and do not have a systematic way to adjust model parameters. Also, other purely simulation-based approaches do not perform well in large-scale applications. In this study, we focus on utility-maximizing nested logit activity-based model-systems and calibrating aggregate statistics such as activity shares, mode shares, time-dependent & mode-specific OD flows, and time-dependent & mode-specific sensor counts. We formulate the calibration problem as a simulation-based optimization problem and propose a stochastic gradient-based solution procedure to solve it. The solution procedure relies on microsimulation to calculate expectations of the aggregate statistics of interest to the calibration problem. Additionally, we derive approximate analytical expressions for the gradient of the objective function —that are evaluated through microsimulation on mini-batches of the population. The proposed solution procedure is sensitive to the fundamental structure of the activity-based model-system and is non-myopic in considering the dependencies across its model components. The formulated optimization problem is non-convex, highly nonlinear, and potentially has multiple-minima. Lastly, we show —through a real-world application— that the proposed solution procedure outperforms other state-of-the-art purely simulation-based optimization approaches in terms of computational efficiency, stability, and convergence. We also compare various gradient-based solution algorithms to determine the best algorithm to update the parameters. This work has the potential to facilitate wider and easier application of activity-based model-systems.

97 MATHEMATICS AND COMPUTING↗

Final Report for CSP Tower Public Opinion and Education Project

As part of the CSP Plant Optimization Study for the California Power Market (DE-EE0009809) the project wanted to understand the public’s opinion of the technology and and explore the types of community engagement that is needed to support such development. The initial objectives for the public perception activities were twofold. The primary objective was to gather public opinion and feedback from the communities living near the two operating solar power tower plants, Ivanpah and Crescent Dunes, and to distill lessons learned from the community engagement conducted before, during, and after the plants were developed to inform future development. This included outreach to nearby airports. The other objective was to gauge public opinion about large-scale solar, specifically CSP towers, to start educating the public on the benefits and to begin building relationships with communities of interest, initially targeting the Kingman, Arizona area. The objectives shifted after the first exploratory trip to Ivanpah and Kingman, however, as it became apparent that it would be challenging to gather public opinion from the community surrounding Ivanpah that could be useful for other community profiles, and the company wasn’t ready to address the concerns in Kingman. Another area was chosen, therefore, to represent those where future development is possible. The recently published Lawrence Berkeley National Lab Perceptions of Large-Scale Solar Project Neighbors Study exemplified public perception polling based on social science and served as the foundation for the survey questions taken to the field. The intent was to ensure that people knew their input was valued and that the time they spent was valuable for the participant as well. It has been noted in the literature that in-person interaction has greater benefits than activities online or via mail, as well as limits the expense.

14 SOLAR ENERGY↗

Assessing Cloud and Precipitation Properties on Temporal and Spatial Scales Using LASSO Simulations over ENA

Low clouds and precipitation representation remain a major source of uncertainty in Earth System Models (ESMs), particularly due to challenges in representing their sub-grid variability and scale-dependent sampling. This study evaluates the performance of preliminary simulations from the Large-Eddy Simulation (LES) ARM Symbiotic Simulation and Observation (LASSO) project over the Eastern North Atlantic (ENA), with a focus on liquid water path (LWP), ice water path (IWP), cloud fraction (CF), and surface precipitation simulated across closed-cell, open-cell, and transitional cloud regimes. Using LES (100 m horizontal grid spacing) driven by ERA5 and MERRA-2 reanalyses, we assess the representativeness of ground-based point observations by analyzing their correspondence to model-resolved spatial and temporal means. Results suggest that observational sampling of at least 6 hours is required to achieve consistency with domain-scale averages, in particular for observations that exhibit pronounced sub-grid heterogeneity, such as precipitation. ERA5-forced simulations exhibit improved spatial coherence and agreement with domain-averaged quantities when compared to MERRA-2 runs, with performance discrepancies largest for convective cloud conditions due to differences in forcing fidelity and temporal resolution. These findings highlight the importance of regime-aware model evaluation strategies and potentially demonstrate how LES can inform observation-model comparison practices and the development of cloud and precipitation parameterizations in ESMs.

Liang, Jiakun [University of Hawai'i at Manoa] (OR↗

Seismic Tomographic Modeling of the Crust and Upper Mantle beneath Israel and the Middle East: Improved Resolution through Optimized Model Parameterization

Accurate regional seismic travel-time (RSTT) predictions rely on regional phases (e.g., Pg, Lg, Pn, Sn ) to account for 3D effects in the crust and upper mantle that are not captured by 1D models traditionally used for real-time location. The RSTT prediction model accounts for regional-scale crust and upper mantle structure globally by incorporating regional seismic phases into its travel-time calculations. Previous versions of the RSTT model have used a constant grid cell size of 1°. To improve the tomographic accuracy of recovering velocity structure at regional scales, we, in this study, perform data-driven grid refinement on the RSTT model down to a 0.125° grid (~14 km) in pursuit of two main goals: (1) to test the limits of RSTT capability and accuracy of determined velocity structure through variable grid refinement and (2) to image smaller structures in Israel and the Middle East and illuminate upper mantle dynamics operating in this complex tectonic area. We investigate the effects of model parameterization as grid cell size decreases and the trade-offs between recovered velocity structures. Our final dataset includes 4751 events and 499 stations that recorded 79,344 Pn and 7489 Pg . The variable grid refinement method allows recovery of finer-scale velocity structures and reduces travel-time residuals in areas with the highest data coverage. At smaller grid cell sizes, longer paths need to be upweighted to stabilize the inversion. Results illuminate tectonic features undefined in coarser grid-size models; in particular, we observe mantle perturbations related to the subduction zone around the Cyprian arc and crustal anomalies near the Dead Sea fault and throughout the Anatolian plate.

58 GEOSCIENCES↗

Orographic Cirrus and Its Radiative Forcing in NCAR CAM6

Cirrus clouds play an important role in the Earth's radiative energy budget, thereby affecting the climate state and climate change. Orographic gravity wave (OGW)-induced sub-grid scale vertical velocity (i.e., cooling rate) is not resolved by large-scale models and its impact on ice formation in cirrus clouds is not well quantified. In this study, one sub-grid scale OGW scheme (e.g., McFarlane) is used in the Community Atmosphere Model version 6 (CAM6) to generate vertical velocity variance (σw) for cirrus formation. Results from the default model and simulations with the OGW-induced σw are evaluated against the DOE ARM Small Particles in Cirrus (SPARTICUS) campaign observations. The OGW based on the McFarlane scheme increases the sub-grid scale σw over mountains compared to the default model and improves the model agreement with the SPARTICUS observations. Larger σ w due to OGWs can trigger more frequent homogeneous nucleation in orographic cirrus and generates a higher number concentration of ice crystals observed during the SPARTICUS campaign. Moreover, our evaluation of the model simulations against satellite observations indicates that the McFarlane scheme generates high in-cloud ice number concentrations (>200 L -1 ) in the upper troposphere over mountains and high plateaus at mid- and high-latitudes of the winter hemisphere as shown in the observations. More ice crystals with smaller sizes absorb more infrared radiation (+0.523 ± 0.125 W m -2 ). The net radiative cloud forcing change at the top of the atmosphere is +0.330 W m -2 due to the orographic cirrus.

54 ENVIRONMENTAL SCIENCES↗

Multiscale Concurrent Atomistic-Continuum (CAC) modeling of multicomponent alloys

We report strengthening in complex multicomponent systems such as solid solution alloys is controlled primarily by the dynamic interactions between dislocation lines and heterogeneously distributed solute species. Modeling of extended defect length scales in such multicomponent systems becomes prohibitively expensive, motivating the development of reduced order approaches. This work explores the application of the Concurrent Atomistic-Continuum (CAC) method to model dislocation mobility in random alloys at extended length scales. By employing recently developed average-atom interatomic potentials, the average “bulk” material response in coarse-grained regions interacts with true random solute species in the atomistic-scale domain. We demonstrate that spurious stresses in domain resolution transition regions are eliminated entirely due to the CAC formulation. Simultaneously, the key details of local stress fluctuation due to randomness in the dislocation core region are captured, and fluctuating stress smoothly decays to the long-range dislocation stress field response. Dislocation mobility calculations, for line lengths over 400 nm, are computed as a function of alloy composition in the model FeNiCr system and compared to full molecular dynamics (MD). The results capture the composition-dependent trends, while reducing degrees of freedom by nearly 40%. This approach can be readily extended to any system described by an EAM potential and facilitates the study of large-scale defect dynamics in complex solute environments to support computational alloy design.

36 MATERIALS SCIENCE↗

First measurement of electron neutrino scattering cross section on argon

We report the first electron neutrino cross section measurements on argon, based on data collected by the ArgoNeuT experiment running in the GeV-scale NuMI beamline at Fermilab. A flux-averaged $\nu_e + \overline{\nu}_e$ total and a lepton angle differential cross section are extracted using 13 $\nu_e$ and $\overline{\nu}_e$ events identified with fully-automated selection and reconstruction. We employ electromagnetic-induced shower characterization and analysis tools developed to identify $\nu_e/\overline{\nu}_e$-like events among complex interaction topologies present in ArgoNeuT data ($\langle E_{\bar{\nu}_e} \rangle = 4.3$ GeV and $\langle E_{\nu_e} \rangle = 10.5$ GeV). The techniques are widely applicable to searches for electron-flavor appearance at short- and long-baseline using liquid argon time projection chamber technology. Notably, the data-driven studies of GeV-scale $\nu_e/\overline{\nu}_e$ interactions presented in this Letter probe an energy regime relevant for future DUNE oscillation physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Convolutional Neural Network for Segmenting Micro-X-ray Computed Tomography Images of Wood Cellular Structures

To further enhance the performance of wood products, improved tools are needed to study in situ cellular scale phenomena like mechanical deformations and moisture swelling. Micro-X-ray computed tomography (μXCT) using brilliant synchrotron light sources now has the spatial and temporal resolution for real-time visualization of phenomena in three-dimensional cellular structures. However, the tradeoff for speed includes the loss of intensity contrast between different types of materials within the imaged structure, such as cell wall and air in wood. This loss of contrast prevents traditional histogram-based segmentation methods from being used effectively. A new convolutional neural network (CNN) approach was therefore developed to segment fast μXCT images of wood into cell wall and air volumes. The fast μXCT and segmentation were demonstrated in the study of moisture swelling in loblolly pine (Pinus taeda) earlywood and latewood cellular structures conditioned at 0%, 33%, 75%, and 95% relative humidity (RH). The CNN segmentation results had a mean intersection over union (IoU) metric accuracy of 96%. Initial analysis of the swelling in the latewood revealed cell walls swelled about 25% when conditioned from 0% to 95% RH. Additionally, the widths of ray cell lumina in the transverse plane of latewood could be observed to increase at higher RH. The segmentation method presented here will facilitate future quantitative analyses in in situ μXCT studies of wood and other similar cellular materials.

Arzola-Villegas, Xavier (ORCID:0000000305369766)↗

Transferable dispersion-aware machine learning interatomic potentials for multilayer transition metal dichalcogenide heterostructures

Stacking atomically thin transition metal dichalcogenides (TMDs) into heterostructures enables exploration of exotic quantum phases, particularly through twist-angle-controlled moiré superlattices. These structures exhibit novel electronic and optical behaviors driven by atomic-scale structural reconstruction. However, studying such systems with DFT is computationally demanding due to their large unit cells and van der Waals (vdW) interactions between layers. To address this, we develop a transferable neural network potential (NNP) that includes long-range vdW corrections up to 12Å with minimal overhead. Trained on vdW-corrected DFT data for Mo- and W-based TMDs with S, Se, and Te, the NNP accurately models monolayers, bilayers, heterostructures, and their interaction with h-BN substrates. It reproduces equilibrium structures, energy landscapes, phonon dispersions, and matches experimental atomic reconstructions in twisted WS2 and MoS2/WSe2 systems. We demonstrate that our NNP achieves DFT-level accuracy and high computational efficiency, enabling large-scale simulations of TMD-based moiré superlattices both with and without substrates.

materials for energy and catalysis↗

Development of Data-Driven Models for Performance Prediction and Chemical Dosing of a Full-Scale Controlled Phosphorus Precipitation Reactor

This study evaluated the use of data-driven models to improve control of a struvite precipitation reactor that removes phosphorus from wastewater while producing a fertilizer product. The researchers developed predictive models for influent orthophosphate concentration, effluent orthophosphate concentration, and phosphorus removal using operational data from a full-scale MagPrex™ reactor at a water resource recovery facility in Denver, Colorado. Model predictions were used to recommend magnesium chloride dosing adjustments needed to achieve a target effluent phosphorus concentration. Several machine learning approaches were tested, with ridge regression providing the best predictions for influent orthophosphate concentration and phosphorus removal, and XGBoost providing the best predictions for effluent orthophosphate concentration. Simulation results indicated that the decision-support approach could correctly identify dosing adjustments in most cases and reduce chemical use. Full-scale implementation achieved lower accuracy due to changing operating conditions and limited historical data in some operating ranges. Here, the results demonstrate the potential of data-driven tools to support phosphorus recovery process control while also identifying practical limitations that affect deployment in full-scale systems.

42 ENGINEERING↗

An Open-Source Parallel EMT Simulation Framework

As the integration level of inverter-based resources (IBRs) increases, ensuring the reliable operation of the bulk power systems requires the use of electromagnetic transient (EMT) simulation tools to identify and mitigate system-wide stability risks. Conducting EMT studies for large-scale, IBR-rich grids, however, is challenging due to the inherent computational bottleneck caused by the underlying high-fidelity models and required small time steps. This paper introduces ParaEMT: an open-source, generic EMT simulation framework designed to accelerate simulations by leveraging advanced parallel computational technologies, such as high-performance computers. This paper presents a comprehensive exposition of ParaEMT, covering its modeling library, simulation strategy, framework structure, operational procedures, and auxiliary features, alongside its extensible parallel computational architecture. Notably, ParaEMT is a publicly accessible and modularized framework written in Python, thereby facilitating future development and the integration of new models and algorithms. The accuracy and efficiency of ParaEMT are demonstrated by rigorous validations via multiple case studies.

electromagnetic transient simulation↗

An Open-Source Parallel EMT Simulation Framework: Preprint

As the integration level of inverter-based resources (IBR) increases, ensuring the reliable operation of the bulk power systems requires the use of electromagnetic transient (EMT) simulation tools to identify and mitigate system-wide stability risks. Conducting EMT studies for large-scale, IBR-rich grids, however, is challenging due to the inherent computational bottleneck caused by the underlying high-fidelity models and required small time steps. This paper introduces ParaEMT: an open-source, generic EMT simulation framework designed to accelerate simulations by leveraging advanced parallel computational technologies, such as high-performance computers. This paper presents a comprehensive exposition of ParaEMT, covering its modeling library, simulation strategy, framework structure, operational procedures, and auxiliary features, alongside its extensible parallel computational architecture. Notably, ParaEMT is a publicly accessible and modularized framework written in Python, thereby facilitating future development and the integration of new models and algorithms. The accuracy and efficiency of ParaEMT are demonstrated by rigorous validations via multiple case studies.

electromagnetic transient simulation↗