Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spreading models”

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 559 records · Page 31

Effect of microstructure and temperature on nuclear graphite oxidation using the 3D Random Pore Model

This work compared the 3D Random Pore Model (3D-RPM) with experimental and characterization data to systematically study the effect of nuclear graphite microstructure and air oxidation temperature. It is well known that oxidation-induced weight loss at elevated temperatures degrades graphite structure and properties; however, a fundamental understanding of the role of graphite microstructure is still unclear. In this work, three diverse grades of nuclear graphite—IG-110, NBG-18, and PCEA—were examined and tested at air oxidation temperatures of 600, 650, 700, and 750 °C following ASTM D7542. Here, the 3D-RPM reproduced the microstructure and temperature dependence of mass loss curves for these grades. For the first time, measurements and modeling were combined to show how bulk density and the amount and spatial distribution of open and closed porosity affects oxidation. IG-110 was less dense and had an open pore network that is finer and more uniformly spread, and showed fastest oxidation. PCEA porosity was less uniform and showed less oxidation, while NBG-18 was more dense, had the least fine and uniform porosity, and showed the slowest oxidation. In general, oxidation proceeds faster if open porosity is more uniformly distributed and can incrementally access closed pore surface area with more ease.

36 MATERIALS SCIENCE↗

Hyperdiffusion of dust particles in a turbulent tokamak plasma

In this work, the effect of plasma turbulence on the trajectories of dust particles is investigated for the first time. The dynamics of dust particles is computed using the ad hoc developed Dust Injection Simulator code, using a 3D turbulent plasma background computed with the TOKAM3X code. As a result, the evolution of the particle trajectories is governed by the ion drag force, and the shape of the trajectory is set by the Stokes number St∝a d /n 0 , with a d the dust radius and n 0 the density at the separatrix. The plasma turbulence is observed to scatter the dust particles, exhibiting a hyperdiffusive regime in all cases. The amplitude of the turbulent spread of the trajectories Δr 2 is shown to depend on the ratio Ku/St, with Ku∝u rms the Kubo number and urms the fluctuation level of the plasma flow. These results are compared with a simple analytical model, predicting Δr 2 ∝(Ku/St) 2 t 3 , or Δr 2 ∝(u rms n 0 /a d ) 2 t 3 . As the dust is heated by the plasma fluxes, thermionic emission sets the dust charge, originally negative, to slightly positive values. This results in a substantial reduction of the ion drag force through the suppression of its Coulomb scattering component. The dust grain inertia is then no longer negligible and drives the transition from a hyperdiffusive regime toward a ballistic one.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

PRIME

SAND2021-0565 O PRIME is a modeling framework designed for the real-time characterization and forecasting of partially observed epidemics. The method is designed to help guide medical resource allocation in the early epoch of the outbreak. Characterization is the estimation of infection spread parameters using daily counts of symptomatic patients. The estimation problem is posed as one of Bayesian inference and solved using a Markov Chain Monte Carlo technique. The framework can accommodate multiple epidemic waves and can help identify different disease dynamics at the regional, state, and country levels. Examples are provided using publicly available COVID-19 data. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Safta, Cosmin↗

Remote Sensing-Based Estimation of Advanced Perennial Grass Biomass Yields for Bioenergy

A sustainable bioeconomy would require growing high-yielding bioenergy crops on marginal agricultural areas with minimal inputs. To determine the cost competitiveness and environmental sustainability of such production systems, reliably estimating biomass yield is critical. However, because marginal areas are often small and spread across the landscape, yield estimation using traditional approaches is costly and time-consuming. This paper demonstrates the (1) initial investigation of optical remote sensing for predicting perennial bioenergy grass yields at harvest using a linear regression model with the green normalized difference vegetation index (GNDVI) derived from Sentinel-2 imagery and (2) evaluation of the model’s performance using data from five U.S. Midwest field sites. The linear regression model using midsummer GNDVI predicted yields at harvest with R2 as high as 0.879 and a mean absolute error and root mean squared error as low as 0.539 Mg/ha and 0.616 Mg/ha, respectively, except for the establishment year. Perennial bioenergy grass yields may be predicted 152 days before the harvest date on average, except for the establishment year. The green spectral band showed a greater contribution for predicting yields than the red band, which is indicative of increased chlorophyll content during the early growing season. Although additional testing is warranted, this study showed a great promise for a remote sensing approach for forecasting perennial bioenergy grass yields to support critical economic and logistical decisions of bioeconomy stakeholders.

09 BIOMASS FUELS↗

Bound state dynamics in the long-range spin- 1 2 XXZ model

Experimental platforms based on trapped ions, cold molecules, and Rydberg atoms have made possible the investigation of highly nonlocal spin-1 / 2 Hamiltonians with long-range couplings. Here, we study the effects of such nonlocal couplings in the long-range spin-1 / 2 XXZ Heisenberg Hamiltonian. Here, we calculate explicitly the two-spin energy spectrum, which describes all possible energetic configurations of two spins pointing in a specific direction embedded in a background of spins with opposite orientation. For fast decay of the spin-spin couplings, we find that the two-spin energy spectrum is characterized by well-defined discrete values, corresponding to bound states, separated by a set of continuum states describing the scattering region. In the deep long-range regime instead, the bound states disappear as they get incorporated by the scattering region. The presence of two-spin bound states results to be crucial to determine both two- and many-spin dynamics. On one hand, radically different two-spin spreadings can be observed by tuning the decay of the spin couplings. On the other hand, two-spin bound states enable the dynamical stabilization of effective antiferromagnetic states in the presence of ferromagnetic couplings. Finally, we propose a novel scheme based on a trapped-ion quantum simulator to experimentally realize the long-range XXZ model and to study its out-of-equilibrium properties.

74 ATOMIC AND MOLECULAR PHYSICS↗

A semi-automated algorithm for designing stellarator divertor and limiter plates and application to HSX

We present a semi-automated algorithm for designing three-dimensional divertor or limiter plates targeting low heat loads. The algorithm designs the plates in two stages: firstly, the parallel heat flux distribution is caught on vertically-inclined plates at one or several toroidal locations. Secondly, the power per unit area is reduced by stretching, tilting and bending the plates toroidally. Heat transport is modelled using the EMC3-Lite code, which uses an anisotropic diffusion model. We apply this scheme to HSX, a medium-sized stellarator located at the University of Wisconsin–Madison. Starting from the current machine with an extended vessel wall, we construct plates which are able to effectively catch and spread the heat for three different magnetic configurations. The scheme has a computational cost in the order of tens of CPU-minutes, making it a powerful tool for semi-automated plasma-facing component design in three-dimensional environments.

anisotropic diffusion↗

Gambusia holbrooki Survive Shear Stress, Pressurization and Avoid Blade Strike in a Simulated Pumped Hydroelectric Scheme

Pumped hydroelectric energy storage (PHES) projects are being considered worldwide as a means of achieving political renewable energy targets in a way that stabilises baseload energy supply from often intermittent renewable energy sources. Unlike a conventional hydroelectric system that only pass water in a downstream direction, a feature of PHES is that it relies on the bi-directional flow of water. In some cases this flow can be across different waterbodies or catchments, posing a risk of inadvertently expanding the range of aquatic animals like fish. The risk of this happening depends on the likelihood of survival of individuals, which remains poorly understood for turbines that are pumping rather than generating. This study quantified the survival of a globally widespread and invasive poeciliid fish, Eastern gambusia (Gambusia holbrooki) when exposed hydraulic stresses characteristic of what would be experienced through a PHES during the pumping phase. A shear flume and hyperbaric chamber were used to expose fish to different strain rates and rapid and sustained pressurisation. A blade strike model was also used to predict survival of fish passing through a Francis dual turbine / pump. The ranges simulated were based on design and operational conditions provided for a PHES scheme being proposed in south-eastern Australia. All gambusia tested survived extremely high (up to 7600 kPa gauge pressure) pressurisation, high levels of shear stress (up to 1853 s -1 ), and the majority (> 93 %) were unlikely to be struck by a turbine blade. Given their tolerance to these extreme simulated stresses, we conclude that gambusia will likely survive passage through a PHES scheme. Therefore, where a new PHES poses the risk of inadvertently expanding the range of gambusia or similar poeciliid species, measures to minimise their spread or mitigate their ecosystem impacts should be considered.

25 ENERGY STORAGE↗

Comparison of Numerical Methods for Modeling the Wave Field Effects Generated by Individual Wave Energy Converters and Multiple Converter Wave Farms

This numerical study compares the wave field generated by the spectral wave action balance code, SNL-SWAN, to the linear-wave boundary-element method (BEM) code, WAMIT. The objective of this study is to assess the performance of SNL-SWAN for modeling wave field effects produced by individual wave energy converters (WECs) and wave farms comprising multiple WECs by comparing results from SNL-SWAN with those produced by the BEM code WAMIT. BEM codes better model the physics of wave-body interactions and thus simulate a more accurate near-field wave field than spectral codes. In SNL-SWAN, the wave field’s energy extraction is modeled parametrically based on the WEC’s power curve. The comparison between SNL-SWAN and WAMIT is made over a range of incident wave conditions, including short-, medium-, and long-wavelength waves with various amounts of directional spreading, and for three WEC archetypes: a point absorber (PA), a pitching flap (PF) terminator, and a hinged raft (HR) attenuator. Individual WECs and wave farms of five WECs in various configuration were studied with qualitative comparisons made of wave height and spectra at specific locations, and quantitative comparisons of the wave fields over circular arcs around the WECs as a function of radial distance. Results from this numerical study demonstrate that in the near-field, the difference between SNL-SWAN and WAMIT is relatively large (between 20% and 50%), but in the far-field from the array the differences are minimal (between 1% and 5%). The resultant wave field generated by the two different numerical approaches is highly dependent on parameters such as: directional wave spreading, wave reflection or scattering, and the WEC’s power curve.

16 TIDAL AND WAVE POWER↗

Stochastic social behavior coupled to COVID-19 dynamics leads to waves, plateaus, and an endemic state

It is well recognized that population heterogeneity plays an important role in the spread of epidemics. While individual variations in social activity are often assumed to be persistent, that is, constant in time, here we discuss the consequences of dynamic heterogeneity. By integrating the stochastic dynamics of social activity into traditional epidemiological models, we demonstrate the emergence of a new long timescale governing the epidemic, in broad agreement with empirical data. Our stochastic social activity model captures multiple features of real-life epidemics such as COVID-19, including prolonged plateaus and multiple waves, which are transiently suppressed due to the dynamic nature of social activity. The existence of a long timescale due to the interplay between epidemic and social dynamics provides a unifying picture of how a fast-paced epidemic typically will transition to an endemic state.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Lessons Learned from AskGDR: Usage and Impact Analysis of the Geothermal Data Repository's AI Research Assistant: Preprint

In October of 2024, the Department of Energy's (DOE) Geothermal Data Repository (GDR) team officially launched AskGDR, an AI research assistant resulting from the integration of a Large Language Model (LLM) with the metadata and supporting documents associated with GDR datasets. AskGDR allows GDR users to ask deeper questions about the origin of datasets, the methods used to collect them, and the findings they help support. Using Retrieval Augmented Generation (RAG), AskGDR can be used to summarize findings spread across dozens of papers and technical reports or to extract relevant information describing a single data field. However, generative AI is experimental. The National Renewable Energy Laboratory (NREL) has been collecting metrics on AskGDR and documenting lessons learned during its deployment. This paper will outline the efficacy and impact of AskGDR through analysis of its use, operating costs, number and types of questions asked, and the quality of answers provided.

15 GEOTHERMAL ENERGY↗

Chromatic effects on the PSF and shear measurement for the Roman Space Telescope High-Latitude Wide Area Survey

Weak gravitational lensing (WL) is a key cosmological probe that requires precise measurement of galaxy images to infer shape distortions, or shear, and constrain cosmology. Accurate estimation of the point spread function (PSF) is crucial for shear measurement, but the wavelength dependence of the PSF introduces chromatic biases that can systematically impact shear inference. We focus on biases arising from spectral energy distribution (SED) differences between stars, used for PSF modelling, and galaxies, used for shear measurement. We investigate these effects in Roman’s four design reference mission WL bands (Y106, J129, H158, F184) and wide filter (W146). Using Roman-like image simulations, we quantify the induced shear biases and compare them to requirements on those biases. Multiplicative biases over all galaxies hover around ~0.2 per cent in the WL bands and 2 per cent in the wide filter, exceeding the mission requirement of |m| < 0.032 per cent and relaxed requirement of |m| < 0.1 per cent⁠. In individual redshift bins, biases can reach 0.4–0.9 per cent for the WL bands and 3–6 per cent for the wide filter. Additive biases remain acceptable in the WL bands but exceed systematic limits in the wide filter. We develop and test PSF-level corrections, showing that a first-order correction reduces biases within survey requirements for the WL bands; however, higher-order terms are necessary for the wide filter. Our results highlight the necessity of chromatic corrections for precision WL with Roman and provide a framework for mitigating these biases. Finally, we compare analytical colour-based corrections to self-organizing maps and find that both methods effectively reduce biases.

79 ASTRONOMY AND ASTROPHYSICS↗

Local Electricity Markets

This book provides a thorough, yet concise, review on the current status of LEM development all around the world, including the most promising research models and opportunities that are being proposed; an overview of the regulatory issues on the different Continents a discussion on the current infrastructure (both hardware and software) that is implemented and that is expected to facilitate the implementation and wide spread application of local energy markets (e.g., resulting from the investments already made in SG deployment); a review of current applications and practical implementations of LEMs; and a wrap-up and discussion on the most relevant paths for future research and development in this field of study.

day-ahead market, peer-to-peer, bilateral market, ↗

Joint Modeling of GD-1 and C-19 as Old Streams

DESI observational data for the GD-1 and C-19 streams are compared to stream simulations in a common evolving multi-halo potential of a Milky Way-like galaxy based on a cosmological simulation. The goal is to find the best match of the stream velocity spread and the density power spectrum stream density to simulations having either CDM or WDM subhalos. The cocoon velocity width integrated over the length of the stream is independent of orbital blurring along the stream and the power spectrum integrates over the width of the stream, sidestepping the geometric details of the streams. Streams develop from star clusters inserted at $\simeq$1 Gyr after the Big Bang and evolved for 13 Gyr to their current orbital positions. Streams in a CDM subhalo population provide the best match to the velocity width, with streams younger than 10 Gyr ruled out as insufficiently hot. The progenitor star cluster masses, which determine the fraction of stars released at late times which comprise the stream core, are found to be $\simeq 8\times 10^4 M_\odot$ for GD-1 and $\simeq 4\times 10^4 M_\odot$ for C-19, although the mass depends on the star cluster half mass radius. Stream heating leads to stream lumpiness which is measurable for the relatively large and clean GD-1 dataset. The stream density power spectrum measured along the length of the DESI GD-1 sample is in good agreement with CDM simulations, with 1.7 to 1.9 times more power than WDM 7 keV and 5.5 keV simulations.

Carlberg, Raymond G. [Toronto U.] (ORCID:000000027↗

CFD Modeling of an 81" Single Slope Barrier Wall Drainage Window for on Grade and Sag Conditions

The main goal of this study is to design a drainage window in an 81" single slope barrier wall for on-grade and sag conditions that could replace the standard ODOT I-3 barrier catch basin. The ODOT design guidelines are as follows: a) Examine configurations of the opening to the flume behind the wall under the wall section to maximize conveyance of flow to a roadside ditch section. b) Examine window lengths and configurations to maximize hydraulic efficiency on a grade up to an absolute maximum window length of 20 ft. c) Use a minimum opening height of 4 in. and consider tapering to an increased height at the outlet to the flume to reduce the potential for debris clogging beneath the barrier section. On-grade locations are analyzed for longitudinal slope up to 5% and spread up to 10 ft. Sag locations, represented by pavement at zero longitudinal slope, are analyzed for water depths up to 18 in. The roadway cross section consists of a 10-ft shoulder at a 4.0% cross slope and pavement at a 1.6% cross slope. All surfaces are represented with Manning’s n = 0.015. Computational fluid dynamics (CFD) is used to perform the analysis. The modeling approach follows previous studies by the authors. The new design is similar to a slotted inlet: a rectangular opening cut into the bottom of a barrier (e.g., a temporary barrier or single slope barrier wall), which typically has a uniform cross section. Findings from laboratory tests of flow under Florida DOT barrier walls were used for additional validation of the modeling approach used in this study.

42 ENGINEERING↗

3-D Radiological Data Acquisition, Visualization and Modeling - 20211

The U.S. Army Corps of Engineers (USACE) was tasked to investigate and remediate low activity radiological contamination from research and production of the nation's first nuclear weapons at the former DuPont Chambers Works Formerly Utilized Sites Remedial Action Program (FUSRAP) site (DuPont). The DuPont site had several buildings used for the Manhattan project that were demolished in the 1940's and 1950's apparently using heavy earthmoving equipment. Some of the contaminated rubble from the demolition appears to have been spread out by this equipment resulting in somewhat random scattering of radiologically contaminated soil and debris along with aqueous spills. Traditional investigative methods such as soil borings, test pits and 2-dimensional gamma walkovers were only partially successful in delineating the radiological contamination at the site. It was feared that even 'chasing' the contamination during remediation would miss contamination if the demolition resulted in discontinuous trails of radiologically contaminated soils. In evaluating the data generated over the interceding decades, the USACE determined that a better method to collect and process the remedial action radiological data was needed to enable the project team to optimize predictive planning and meet documentation expectations. The purpose of this paper is to provide an overview of the effort and progress to combine and organize radiological survey methods into a highly flexible sampling, modeling, and decision analysis approach that emphasizes the quality of decision-making during remediation. This innovative system blends multiple tools to develop a methodology that can extend MARRSIM [1] into the subsurface and provide tools that can be applied to other sites. (authors)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Asymmetry of the Ferroelectric Phase Transition in BaTiO 3

Symmetry changes during phase transformations fundamentally determine the behavior of thermodynamic systems, governing phenomena as diverse as water evaporation, fermion condensation, and epidemic spreading. Phase transitions are conventionally divided into two classes, and this classification is typically assumed to remain invariant upon reversing the transition. Here, an asymmetric phase transformation is uncovered in the ferroelectric–paraelectric transition of single-crystal BaTiO 3 , a model system for first-order transitions. Under slow temperature variation (≤0.1 °C min −1 ), thermodynamic, dielectric, and domain-structure measurements reveal that the ferroelectric-to-paraelectric transition exhibits latent heat, phase coexistence, and a discontinuous order parameter, while these signatures are absent upon cooling. Complementary phase-field simulations demonstrate similar behavior, attributing it to distinct elastic strain energy accumulation and release during heating and cooling. These findings reveal a first-order character upon heating but a second-order-like behavior upon cooling, challenging the conventional paradigm of symmetric phase-transition classification and suggesting new possibilities for ferroelectric-based energy storage.

Hershkovitz, Asaf [Technion-Israel Institute of Te↗

Spatiotemporal dynamics of ionic reorganization near biological membrane interfaces

Electrical signals in excitable cells involve spatially localized ionic fluxes through ion channels and pumps on cellular lipid membranes. Common approaches to understand how these fluxes spread assume that the membrane and the surrounding electrolyte comprise an equivalent circuit of capacitors and resistors, which ignores the localized nature of transmembrane ion transport, the resulting ionic gradients and electric fields, and their spatiotemporal relaxation. Here, we consider a model of localized ion pumping across a lipid membrane, and use theory and simulation to investigate how the electrochemical signal propagates spatiotemporally in and out of plane along the membrane. The localized pumping generates long-ranged electric fields with three distinct scaling regimes along the membrane: a constant potential near-field region, an intermediate monopolar region, and a far-field dipolar region. Upon sustained pumping, the monopolar region expands radially in plane with a steady speed that is enhanced by the dielectric mismatch and the finite thickness of the lipid membrane. For unmyelinated lipid membranes in physiological settings, we find remarkably fast propagation speeds of ∼ 40 m / s , allowing faster ionic reorganization compared to bare diffusion. Together, our paper shows that transmembrane ionic fluxes induce transient long-ranged electric fields in electrolyte solutions, which may play hitherto unappreciated roles in biological signaling. Published by the American Physical Society 2025

Row, Hyeongjoo (ORCID:000000033623512X)↗

Hyperdiffusion of dust particles in a turbulent tokamak plasma

The effect of plasma turbulence on the trajectories of dust particles is investigated for the first time. The dynamics of dust particles is computed using the ad-hoc developed Dust Injection Simulator code, using a 3D turbulent plasma background computed with the TOKAM3X code. As a result, the evolution of the particle trajectories is governed by the ion drag force, and the shape of the trajectory is set by the Stokes number $St\propto a_d/n_0$, with $a_d$ the dust radius and $n_0$ the density at the separatrix. The plasma turbulence is observed to scatter the dust particles, exhibiting a hyperdiffusive regime in all cases. The amplitude of the turbulent spread of the trajectories $\Delta r^2$ is shown to depend on the ratio $Ku/St$, with $Ku\propto u_{rms}$ the Kubo number and $u_{rms}$ the fluctuation level of the plasma flow. These results are compared with a simple analytical model, predicting $\Delta r^2\propto (Ku/St)^2t^3$, or $\Delta r^2\propto (u_{rms}n_0/a_d)^2t^3$. As the dust is heated by the plasma fluxes, thermionic emission sets the dust charge, originally negative, to slightly positive values. This results in a substantial reduction of the ion drag force through the suppression of its Coulomb scattering component. The dust grain inertia is then no longer negligible, and drives the transition from a hyperdiffusive regime towards a ballistic one.

dust↗