Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,351 records · Page 75

Precision beam diagnostics at the NuMI facility using muon monitor observations

The Neutrinos at the Main Injector (NuMI) facility at Fermilab delivers an intense neutrino beam for multiple experiments by producing pions that decay into neutrinos, muons, and other particles. Magnetic horns—the primary pion focusing elements in the NuMI beamline—exhibit predominantly linear optics, enabling a predictable relationship between the proton beam and the resulting pion and muon phase spaces. This study has two primary objectives: first, to evaluate and confirm the linearity of the horn focusing mechanism using analytical models and numerical simulations; and second, to demonstrate that key beam parameters—such as proton beam intensity, beam position on target, and horn current—can be extracted from muon monitor observations within this linear optics framework. Using a machine learning model trained on spill-by-spill muon monitor data, we infer the horn current with a precision of ±0.05%, the beam intensity with ±0.1%, and the beam position on target with ±0.018⁢ mm horizontally and ±0.013⁢ mm vertically. This approach provides a reliable cross-check of beam parameters, helping to reduce systematic uncertainties that are critical for future experiments such as the Deep Underground Neutrino Experiment, which will rely on the neutrino beam produced by the Long-Baseline Neutrino Facility.

Beam control↗

A Review of Modeling Approaches for Predicting Frost Growth and Defrosting on Tube-Fin Heat Exchangers: Preprint

Frost formation and growth on the evaporator surface is a common process that deteriorates the air-refrigerant heat transfer and restricts airflow. This degrades the performance of the vapor compression system by increasing temperature lift and air-side pressure drop. To accurately predict these effects during coil frosting, as well as the energy use and duration of the defrost process, there is a need to estimate the heat and mass transfer, momentum transport, and solid-liquid and liquid-vapor phase change. Therefore, in the past few decades, continuous effort has been made to model frosting and defrosting processes using approaches ranging from empirical correlations to computational fluid dynamic models. To provide a clearer overview for researchers, engineers, and manufacturers in this field, this paper provides a comprehensive literature review for frosting and defrosting models. The paper begins with theoretical background of frost formation and defrost processes, and then reviews the common modeling approaches in literature and their underlying assumptions when trying to account for various physical phenomenon. Based on the literature review, the most critical modeling effort for frost formation is the determination of frost densification rate and frost growth rate. Various methods to predict these two parameters are reviewed. Empirical correlations commonly used for frost density and thermal conductivity are presented and compared. For the defrost process, various multi-stage models have been proposed with different assumptions. Some assume the presence of air gap between the tube wall and the frost, while others consider the melted frost flow due to gravity. We also review physics-based and empirical approaches to integrate defrost models into heat pump models. We conclude by identifying research gaps and providing recommendations.

defrost↗

Data-Driven RANS Turbulence Closures for Forced Convection Flow in Reactor Downcomer Geometry

Recent progress in data-driven turbulence modeling has shown its potential to enhance or replace traditional equation-based Reynolds-averaged Navier-Stokes (RANS) turbulence models. Here, this work utilizes invariant neural network (NN) architectures to model Reynolds stresses and turbulent heat fluxes in forced convection flows (when the models can be decoupled). As the considered flow is statistically one dimensional, the invariant NN architecture for the Reynolds stress model reduces to the linear eddy viscosity model. To develop the data-driven models, direct numerical and RANS simulations in vertical planar channel geometry mimicking a part of the reactor downcomer are performed. Different conditions and fluids relevant to advanced reactors (sodium, lead, unitary-Prandtl-number fluid, and molten salt) constitute the training database. The models enabled accurate predictions of velocity and temperature, and compared to the baseline k–τ turbulence model with the simple gradient diffusion hypothesis, do not require tuning of the turbulent Prandtl number. The data-driven framework is implemented in the open-source graphics processing unit–accelerated spectral element solver nekRS and has shown the potential for future developments and consideration of more complex mixed convection flows.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

sup3ruhi (Super Resolution for Renewable Resource Data and Urban Heat Islands) [SWR-25-05]

Urban heat is a growing concern, particularly in dense metropolitan areas where high temperatures increase the risk of heat-related illness and drive energy expenses for cooling. Estimating the effects of urban heat remains a challenge due to limitations in describing the built environment, computational constraints, and the need for high-resolution data. This software presents open-source, computationally efficient machine learning methods that enhance the accuracy of urban temperature estimates compared to historical reanalysis data. Models trained using this software have been applied to urban microclimates in Los Angeles and Seattle showing greater accuracy and less bias when compared to low-resolution reanalysis datasets like ERA5 and even when compared to high-resolution mesoscale numerical weather models like WRF with an urban canopy model. Initial findings highlight how machine learning can support urban heat resilience planning by enabling improved assessments of local heat islands, mitigation strategies, and their energy implications. This software is an extension of (sup3r). This software supports the following publication: Buster, Grant, et al. Tackling Extreme Urban Heat: A Machine Learning Approach to Assess the Impacts of Climate Change and the Efficacy of Climate Adaptation Strategies in Urban Microclimates. arXiv:2411.05952, arXiv, 8 Nov. 2024. arXiv.org, https://doi.org/10.48550/arXiv.2411.05952. And has related public data records available at: Buster, Grant, Cox, Jordan, Benton, Brandon, and King, Ryan. Super-Resolution for Renewable Resource Data and Urban Heat Islands (Sup3rUHI). United States: N.p., 16 Oct, 2024. Web. https://data.openei.org/submissions/6220.

Buster, Grant [National Renewable Energy Laborator↗

The three-dimensional structure of black hole accretion flows within the plunging region

ABSTRACT We analyse, using new analytical models and numerical general relativistic magnetohydrodynamic simulations, the three-dimensional properties of accretion flows inside the plunging region of black hole spacetimes (i.e. at radii smaller than the innermost stable circular orbit). These simulations are of thick discs, with aspect ratios of order unity $h/r \sim 1$, and with a magnetic field geometry given by the standard low-magnetization ‘SANE’ configuration. This work represents the first step in a wider analysis of this highly relativistic region. We show that analytical expressions derived in the ‘thin disc’ limit describe the numerical results remarkably well, despite the large aspect ratio of the flow. We further demonstrate that accretion within this region is typically mediated via spiral arms, and that the geometric properties of these spiral structures can be understood with a simple analytical model. These results highlight how accretion within the plunging region is fundamentally two-dimensional in character, which may have a number of observational implications. We derive a modified theoretical description of the pressure within the plunging region which accounts for turbulent heating and may be of use to black hole image modelling.

Astronomy & Astrophysics↗

Spark Channel Dynamics of Electrostatic Discharges

When two differently-charged objects are brought in close proximity to each other, the resulting high electric fields can cause electron avalanche breakdown of the air gap separating the objects, a process known as electrostatic discharge (ESD). If enough initial charge is stored on the objects, the electrical breakdown can proceed to ionize the air to such a degree that a highly conductive filament of plasma forms in the gap, known as a spark channel. The spark electrically bridges the air gap, resulting in a rapid pulse of current that neutralizes the charge difference. The current pulse produces significant heating of the gas in the spark, resulting in dissociation, ionization, thermal radiation, and hydrodynamic expansion. ESD presents a hazard to electrically-sensitive devices, with consequences such as economic losses (e.g. damaged electronics) or unsafe response (e.g. unintended ignition of flammable gas mixtures, initiation of detonators, etc.). For this thesis, the ESD spark is taken to occur between two conducting electrodes, with the spark channel being axisymmetric in a cylindrical coordinate system centered on the channel. An RLC-type circuit is used for the discharge model of the ESD event. The spark is treated as a time-dependent resistance that is in series with a capacitance, an inductance, and (optionally) a load resistance representing a “victim” component under threat from the ESD event. The primary motivation of this work is to use a numerical hydrodynamic model to understand the energy dissipation and transport processes in the spark. The model consists of the compressible Euler equations of mass, momentum, and energy conservation together with an Eddington/P1 approximation for thermal radiation transport. To close the hydrodynamic system, an equation of state (EOS) was fitted from tabular data for air that accounts for the dissociation and ionization of air species. The hydrodynamic equations are solved using a conservative Lagrangian finite volume method. These partial differential equations are coupled to the circuit equations by calculation of the spark resistance via numerical integration of the electrical conductivity of the channel. Computational results are compared against experimental measurements of discharge current and radial density of the spark channel.

42 ENGINEERING↗

A stochastic solver based on the residence time algorithm for crystal plasticity models

Abstract The deformation of crystalline materials by dislocation motion takes place in discrete amounts determined by the Burgers vector. Dislocations may move individually or in bundles, potentially giving rise to intermittent slip. This confers plastic deformation with a certain degree of variability that can be interpreted as being caused by stochastic fluctuations in dislocation behavior. However, crystal plasticity (CP) models are almost always formulated in a continuum sense, assuming that fluctuations average out over large material volumes and/or cancel out due to multi-slip contributions. Nevertheless, plastic fluctuations are known to be important in confined volumes at or below the micron scale, at high temperatures, and under low strain rate/stress deformation conditions. Here, we develop a stochastic solver for CP models based on the residence-time algorithm that naturally captures plastic fluctuations by sampling among the set of active slip systems in the crystal. The method solves the evolution equations of explicit CP formulations, which are recast as stochastic ordinary differential equations and integrated discretely in time. The stochastic CP model is numerically stable by design and naturally breaks the symmetry of plastic slip by sampling among the active plastic shear rates with the correct probability. This can lead to phenomena such as intermittent slip or plastic localization without adding external symmetry-breaking operations to the model. The method is applied to body-centered cubic tungsten single crystals under a variety of temperatures, loading orientations, and imposed strain rates.

42 ENGINEERING↗

Hydrodynamic theory of scrambling in chaotic long-range interacting systems

The Fisher-Kolmogorov-Petrovsky-Piskunov (FKPP) equation provides a mean-field theory of out-of-time-ordered commutators in locally interacting quantum chaotic systems at high energy density. In systems with power-law interactions, the corresponding fractional-derivative FKPP equation provides an analogous mean-field theory. However, the fractional FKPP description is potentially subject to strong quantum fluctuation effects, so it is not clear a priori if it provides a suitable effective description for generic chaotic systems with power-law interactions. Here, in this work, we study this problem using a model of coupled quantum dots with interactions decaying as 1/r α , where each dot hosts N degrees of freedom. The large-N limit corresponds to the mean-field description, while quantum fluctuations contributing to the OTOC can be modeled by 1/N corrections consisting of a cutoff function and noise. Within this framework, we show that the parameters of the effective theory can be chosen to reproduce the butterfly light cone scalings previously found for N=1 and generic finite N. In order to reproduce these scalings, the fractional index μ in the FKPP equation needs to be shifted from the naïve value of μ=2⁢α–1 to a renormalized value μ=2⁢α–2. We provide supporting analytic evidence for the cutoff model and numerical confirmation for the full fractional FKPP equation with cutoff and noise.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ion Trapping Studies and Mitigation Strategies for the EIC ERL-Based Strong Hadron Cooler

An Energy Recovery Linac based strong hadron cooler was previously considered for the Electron-Ion Collider. The required electron beam parameters for variable-energy strong hadron cooling place significant constraints on ion trapping and collective effects. This paper presents initial studies of these constraints through a combination of analytical modelling and numerical simulations of ion production, trapping behaviour, and mitigation strategies. A multi-bunch tracking framework based on ELEGANT with the ionEffects module is used to simulate machine operation over millisecond time scales, corresponding to more than 3 × 10^5 electron bunches. The simulations include modelling of ionisation processes together with transverse electron–ion dynamics, allowing the evolution and accumulation of ions to be investigated. Analytical expressions based on Gaussian beam distributions are used to estimate ion trapping conditions and benchmark the simulation results. A bi-periodic bunch spacing scheme is also investigated as a possible mitigation method by detuning the ion oscillation frequency. These studies provide an initial assessment of ion trapping in the strong hadron cooler and demonstrate possible approaches for reducing beam–ion effects.

Bi, R. [Lancaster University, Cockcroft Institute]↗

Feedbacks Between Surface Deformation and Permafrost Degradation in Ice Wedge Polygons, Arctic Coastal Plain, Alaska

In the past three decades, an abrupt, pan-Arctic acceleration of ice wedge melting has transformed tundra landscapes, spurring the formation of hummock-like features known as high-centered polygons (HCPs). This rapid geomorphic transition profoundly alters regional hydrology and influences surface emissions of CO 2 and CH 4 . In Arctic Alaska, most recent instances of ice wedge degradation have arrested within 15–20 years of inception, stabilizing HCP microtopography. However, feedbacks between ground surface deformation and permafrost stability are incompletely understood, limiting our capacity to predict trajectories of landscape evolution in a still warmer future. In this work, we use field data from a site near Prudhoe Bay, Alaska, to develop a modeling-based framework for assessing the strength of positive (i.e., exacerbating) feedbacks on ice wedge degradation, focusing on the importance of heterogeneity in surface drainage and microtopographic conditions. Our simulations suggest that, when troughs are narrow, positive feedbacks on ice wedge melting (associated with thermokarst pool formation) are relatively weak. Positive feedbacks are markedly stronger beneath wide troughs, such as those that form above older, larger ice wedges. Seasonal thaw abruptly accelerates once a talik begins to form beneath wide and deep thermokarst pools. Once a talik initiates, winter severity and snowpack thickness increase in importance as predictors of thaw intensity in summer. Our results indicate that meter-scale heterogeneity in polygonal microtopography potentially exerts strong, nonlinear controls on thermokarst trajectories. These findings are useful for predicting future thermokarst dynamics and for interpreting the results from coarser-resolution land surface models operating at greater spatial and temporal scales.

54 ENVIRONMENTAL SCIENCES↗

An algorithm for physics informed scan path optimization in additive manufacturing

Site specific microstructure control is a critical research area within the field of additive manufacturing due to its potential to revolutionize part performance. One way to achieve site specific microstructure control is through control of the solidification conditions via the construction of intricate scan paths; however, the search space for such a problem is large. Previous attempts only considered the solidification conditions at the top surface while also requiring either lots of manual-fine tuning or large amounts of computational resources. This paper introduces a general method for scan path optimization which considers the solidification conditions in the bulk of the material without an increase in computational expense. This method consists of three core components:1. A heat transfer model for simulating the temperature field at a given time.2. A surrogate model which takes scan pattern information and temperature data and predicts the solidification conditions of the bulk as well as the meltpool depths for a spot melt.3. A decision algorithm to decide which spot melt should be printed next based on the outputs of the surrogate model.Each of these components can be changed without changing the overall method. Within this work, this method is applied in the creation of an algorithm containing a semi-analytic heat transfer model to simulate the temperature field, a fully convolutional neural network (FCNN) as the surrogate model, and a greedy decision algorithm. The resulting algorithm produced complex scan patterns which gave strong results for simulated microstructure control.

36 MATERIALS SCIENCE↗

Concrete Structure Health Monitoring Using Vibro-acoustic Testing and Machine Learning

Assessment and management of aging concrete structures in nuclear power plants require a more systematic approach than simple reliance on margins of safety in existing codes. Research is required to deploy online monitoring capability to assess the health of concrete structures in nuclear power plants and perform predictive analytics based on heterogeneous measurements, in order to reduce the operation and maintenance costs. The research effort described in this report focuses on the health monitoring of concrete slabs containing reactive aggregates using the vibro-acoustic modulation (VAM) technique and predictive modeling. The research is developing an approach to enable detection, localization, and estimation of the extent of degradation in concrete slabs due to alkali-silica reaction (ASR). The project outcomes are generic and extendable to other degradation modes in concrete. A concrete slab without reinforcements was cast at Vanderbilt University and cured in aggressive conditions to accelerate ASR degradation. In this report, we develop a methodology for the automation of VAM testing, and the construction of machine learning models to facilitate damage localization in concrete structures. We investigate how machine learning models based on physics simulation data can be effectively utilized for damage localization (damage location and depth) and estimation in a concrete slab specimen. Our study indicates that these models could potentially support damage localization and estimation for NPP structures.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Fracture Network Influence on Rock Damage and Gas Transport following an Underground Explosion

Simulations of rock damage and gas transport following underground explosions that omit preexisting fracture networks in the subsurface cannot fully characterize the influence of geo-structural variability on gas transport. Previous studies do not consider the impact that fracture network structure and variability have on gas seepage. In this study, we develop a sequentially coupled, axi-symmetric model to look at the damage pattern and resulting gas breakthrough curves following an underground explosion given different fracture network realizations. We simulate 0.327 and 0.164 kT chemical explosives with burial depths of 100 m for 90 stochastically generated fracture networks. Gases quickly reach the surface in 30% of the higher yield simulations and 5% of the lower yield simulations. The fast breakthrough can be attributed to the formation of connected pathways between fractures to the surface. The formation of a connected damage pathway to the surface is not clearly correlated with the fracture intensity (P32) in our simulations. Breakthrough curves with slower transport are highly variable depending on the fracture network sample. The variability in the breakthrough behavior indicates that ignoring the influence of fracture networks on rock damage, which strongly influences the hydraulic properties following an underground explosion, will likely lead to a large underestimation of the uncertainty in the gas transport to the surface. This work highlights the need for incorporation of fracture networks into models for accurately predicting gas seepage following underground explosions.

58 GEOSCIENCES↗

Material Discovery and Design Principles of Perovskite Oxides for Reversible Solid Oxide Cells (R-SOC)

Reversible solid oxide cells (R-SOCs) are highly efficient devices for energy conversion and storage, capable of operating for both hydrogen utilization and production. In fuel cell mode, an R-SOC consumes hydrogen or natural gas to generate electricity, while in electrolysis mode, it produces hydrogen from steam. The discover of new materials with rapid oxygen surface exchange kinetics and enduring stability is crucial for the economically viable commercialization of R-SOCs. To facilitate this pursuit, we conducted extensive Density Functional Theory (DFT) calculations and developed Machine Learning (ML) models to predict critical catalytic properties essential for R-SOCs, such as oxygen surface exchange/diffusivity, and area-specific resistance (ASR). BaCoxFeyZrzO3-d(BFCZ)(x+y+z=1) emerged as a promising family of electrode materials with high activity and stability, validated through systematic experimental study. Moreover, a robust numerical multiphysics model was developed to optimize materials and microstructure parameters, providing the ability to predict the performance of functional R-SOCs.

Liu, Jian↗

Advances in Clayff Molecular Simulation of Layered and Nanoporous Materials and Their Aqueous Interfaces

As a general-purpose force field for molecular simulations of layered materials and their fluid interfaces, Clayff continues to see broad usage in atomistic computational modeling for numerous geoscience and materials science applications due to its (1) success in predicting properties of bulk nanoporous materials and their interfaces, (2) transferability to a range of layered and nanoporous materials, and (3) simple functional form which facilitates incorporation into a variety of simulation codes. Here, we review applications of Clayff to model bulk phases and interfaces not included in the original parameter set and recent modifications for modeling surface terminations such as hydroxylated nanoparticle edges. Overall, we conclude with a discussion of expectations for future developments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Provably accurate simulation of gauge theories and bosonic systems

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

Tong, Yu↗