Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multi-”

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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Multiscale simulations for multi-continuum Richards equations

In this paper, we study a multiscale method for simulating a dual-continuum unsaturated flow problem within complex heterogeneous fractured porous media. Mathematically, each of the dual continua is modeled by a multiscale Richards equation (for pressure head), and these equations are coupled to one another by transfer terms. On its own, Richards equation is already a nonlinear partial differential equation, and it is exceedingly difficult to solve numerically due to the extra nonlinear dependencies involving the soil water. To deal with multiple scales, our strategy is that starting from a microscopic scale, we upscale the coupled system of dual-continuum Richards equations via homogenization by the two-scale asymptotic expansion, to obtain a homogenized system, at an intermediate scale (level). Based on a hierarchical approach, the homogenization’s effective coefficients are computed through solving the arising cell problems. Furthermore, to tackle the nonlinearity, after time discretization, we use Picard iteration procedure for linearization of the homogenized Richards equations. At each Picard iteration, some degree of multiscale still remains from the intermediate level, so we utilize the generalized multiscale finite element method (GMsFEM) combining with a multi-continuum approach, to upscale the homogenized system to a macroscopic (coarse-grid) level. This scheme involves building uncoupled and coupled multiscale basis functions, which are used not only to construct coarse-grid solution approximation with high accuracy but also (with the coupled multiscale basis) to capture the interactions among continua. These prospects and convergence are demonstrated by several numerical results for the proposed method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Prediction of local concentration fields in porous media with chemical reaction using a multi scale convolutional neural network

The study of solute transport in porous media is of interest in many chemical engineering systems. Some example applications include packed bed catalytic reactors, filtration devices, and batteries. The pore scale modeling of these systems is time consuming and may require large computing resources, for this reason computational fluid dynamics (CFD) simulations are not practical if a large number of simulations is required, like in multiscale modeling, where a model at a large scale calls for pore scale simulations. It has been shown that neural networks can be trained with a dataset of flow simulations and then predict fields orders of magnitude faster, and with less computational resources, in new domains. However, it is crucial to provide the neural network with an effective description of the domain and the undergoing operating conditions to be able to train models that generalize accurately in unseen samples. Therefore, research is needed to employ neural networks in new complex systems. The appropriate training of a network for predicting coupled flow and solute transport processes is an outstanding problem due to the complex interplay between geometry and operating conditions. In this work, we train a multi scale convolutional neural network (MSNet) with a diverse dataset of simulations of transport and chemical reaction in porous media to predict the local concentration fields in images of porous media. Our dataset contains a wide diversity of sphere pack arrangements under different operating conditions (Péclet and Reynolds numbers). Further, we train a robust model by employing different input descriptors that represent the medium and the different operating conditions of each system. Our trained model is able to provide nearly instantaneous predictions, compared to around twenty hours of the CFD workflow, with less than 3.5% error on new geometries and transport conditions. Thus the model could be easily integrated in a multiscale workflow where fast response is needed.

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

A multi-region approach for the analysis of porous materials and gas-solid reactions using USAXS-SAXS-WAXS: application to CaO carbonation

The reaction rate of gas-solid non-catalytic reactions is typically investigated using reactant conversion data over time and ex-situ measurements of the porous solid reactant textural properties; these data enable the experimental estimation of the initial intrinsic reaction rate, whereas the evolution of the textural properties and the reaction rate over time are evaluated theoretically using reaction models. In this work, a different methodology is presented, based on: a) in-situ time-resolved USAXS-SAXS-WAXS measurements of the solid sample, and b) a multi-region modeling approach; this methodology allows for the estimation of the textural properties and the intrinsic reaction rate at any time during the reaction. Referring to a reaction in which a solid product is obtained from a solid reactant and a gaseous reactant, the porous particle is described as consisting of several distinct regions with different microstructural properties, and both a gas-solid reaction model and a SAXS model are derived and applied to the carbonation of porous CaO. The proposed reaction model highlights the role of the inaccessible reactant and accurately predicts the solid reactant conversion versus time profile. The SAXS model accurately predicts the measured linear trends of the Porod invariant and of the pre-factor of the power law scattering profile versus CaO mass fraction; in the absence of chemical reactions, the proposed equations extend classical SAXS theory to porous materials containing macropores and nonporous solid phases in addition to standard nanoscale inhomogeneous regions.

CO2 solid sorbents↗
DOE OSTI2021

A fast multi-fidelity method with uncertainty quantification for complex data correlations: Application to vortex-induced vibrations of marine risers

Here we develop a fast multi-fidelity modeling method for very complex correlations between high- and low-fidelity data by working in modal space to extract the proper correlation function. We apply this method to infer the amplitude of motion of a flexible marine riser in cross-flow, subject to vortex-induced vibrations (VIV). VIV are driven by an absolute instability in the flow, which imposes a frequency (Strouhal) law that requires a matching with the impedance of the structure; this matching is easily achieved because of the rapid parametric variation of the added mass force. As a result, the wavenumber of the riser spatial response is within narrow bands of uncertainty. Hence, an error in wavenumber prediction can cause significant phase-related errors in the shape of the amplitude of response along the riser, rendering correlation between low- and high-fidelity data very complex. Working in modal space as outlined herein, dense data from low-fidelity data, provided by the semi-empirical computer code VIVA, can correlate in modal space with few high-fidelity data, obtained from experiments or fully-resolved CFD simulations, to correct both phase and amplitude and provide predictions that agree very well overall with the correct shape of the amplitude response. We also quantify the uncertainty in the prediction using Bayesian modeling and exploit this uncertainty to formulate an active learning strategy for the best possible location of the sensors providing the high fidelity measurements.

42 ENGINEERING↗
DOE OSTI2021

Kink-antikink collisions and multi-bounce resonance windows in higher-order field theories

In this work, we study collisions of coherent structures in higher-order field-theoretic models, such as the Φ 8 , Φ 10 and Φ 12 ones. The main distinguishing feature, of the example models considered herein, is that the collision arises due to the long-range interacting algebraic tails of these solitary waves. We extend the approach to suitably initialize the relevant kinks, in the additional presence of finite initial velocity, in order to minimize the dispersive wave radiation potentially created by their slow spatial decay. We find that, when suitably initialized, these models still feature the multi-bounce resonance windows earlier found in models in which the kinks bear exponential tails, such as the Φ 4 and Φ 6 field theories among others. Also present is the self-similar structure of the associated windows with three- and more-bounce windows at the edges of two- and lower-bounce ones. Moreover, phenomenological but highly accurate (and predictive), scaling relations are derived for the dependence of the time between consecutive collisions and, e.g., the difference in kinetic energy between the incoming one and the critical one for one-bounces. Such scalings are traced extensively over two-bounce collision windows throughout the three models, hinting at the possibility of an analytical theory in this direction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Experimental measurement of the rapid mixing of fuel and air in a multi-element diffusion (Hencken) burner

Flat flame burners are valuable tools for studying complex combustion phenomena like particle synthesis due to the well-controlled and well-behaved initial and boundary conditions. In this study, the methane/air flame temperature and tracer density profiles of a multi-element diffusion burner equipped with a central tube were characterized using krypton (Kr) x-ray fluorescence at the Advanced Photon Source at Argonne National Laboratory. The experiments were performed using a fuel-to-oxygen equivalence ratio of 0.8 and 3.3% Kr by volume was added to the reactant gases. One- and two-dimensional scans of non-combusting (i.e., cold flow) and combusting conditions, both with and without flow through the central tube were performed. Spatially-resolved Kr number density profiles were obtained at different planar locations and heights above the surface of the burner. The results show rapid mixing occurred near the exit-plane of the burner and increased spatial uniformity was achieved from 0.1 mm to 0.5 mm above the surface of the burner. Further, the data also showed the impact of the central tube flow on the main methane/air flame, including evidence of the central tube flow mixing with the main flow. Flame temperatures were calculated from the fluorescence measurements at various heights above the burner surface and were in good agreement with model predictions for a burner stabilized flame. These high-fidelity measurements provide accurate baseline data that are particularly valuable for understanding the effects of fuel and air mixing on downstream phenomena, such as the formation of condensed-phase products.

42 ENGINEERING↗
DOE OSTI2021

A multi-phase-field model of topological pattern formation during electrochemical dealloying of binary alloys

Porous metal structures produced by electrochemical dealloying have been extensively studied for use in various applications that require high surface areas. A new comprehensive multi-phase-field (MPF) model is proposed to study topological porous patterns formed by spontaneously etching a bulk binary alloy that involves electrochemical reactions, bulk and surface diffusion, ion transport, applied electrode potential, and charge conservation. The governing equations for the alloy-porous cluster-electrolyte system account for a generalized Butler-Volmer electrochemical reaction and are in accordance with the classical nucleation theory. Based on a quantitative examination of the effects of electrode potential and precursor composition, the simulation results reproduce typical phenomena including passive surface dealloying, active porosity evolution, critical potential, and characteristic length scale in two-dimension (2-D) and three-dimension (3-D). To consider more complicated dealloying systems, the evolutions of a bimodal porous metal structure, nanocomposite, and nested porous network with a structured hierarchy are investigated. Finally, the proposed model can be a useful tool for understanding and predicting the morphology evolution of diverse porous structures during electrochemical dealloying.

Hierarchical porous structure↗
DOE OSTI2026

Effects of interstitial oxygen on ω transformations and twin formation in bcc NbTaTiHf multi-principal element alloy from first-principles

Transformation- and twinning-induced plasticity (TRIP and TWIP) have been reported to contribute to the low-temperature deformation of some body-centered cubic (bcc) multi-principal element alloys (MPEAs) containing large fractions of group IV transition metals. The influence of interstitial solutes on the mechanisms underlying these forms of plasticity, however, remains unclear. Using first-principles calculations, we study the effects of interstitial O atoms on the relative stability of bcc and ω phases and on unstable and twin boundary stacking fault energy profiles in a representative bcc MPEA with high group-IV elemental fraction: NbTaTiHf. We find that O additions generally promote the relaxation of ω configurations back to their parent bcc structure, therefore inhibiting ω transformation. Calculations of the Rice parameter for bulk bcc and phases, as well as bcc-ω interfaces, further show that ω formation is a potent embrittlement factor, an effect that is enhanced by O additions, suggesting that the formation of bcc-ω interfaces is energetically preferred over the formation of the bulk ω phase. By contrast, the Rice parameter for twin boundaries indicates that these interfaces do not embrittle the material, even with O atoms at twin boundaries, providing a more favorable pathway for plastic deformation compared to ω transformation.

Density functional theory↗
DOE OSTI2025

Tree-level carbon stock estimations across diverse species using multi-source remote sensing integration

Forests are critical carbon sinks, and remote sensing has been increasingly widely used for forest monitoring and biomass estimations. However, species-specific tree-level studies remain limited. In this study, we demonstrated the feasibility of integrating UAV-based LiDAR with high-resolution optical satellite imagery (0.5 m) to estimate biomass for individual trees across different species. The proposed method accurately estimated biomass for 53 trees (R² = 0.82, rRMSE = 0.44), with species-specific datasets, showing an average 25.2% increase in R² and a 14.8% reduction in rRMSE. A novel vegetation index combining forest structure parameters with vegetation indices (VIs) was developed using high-resolution multispectral satellite data (3 m) to explore its relationship with individual tree biomass. Combining forest structural parameters with VIs further improved estimation accuracy, achieving an R²of 0.89 and an rRMSE of 0.34. Species-specific datasets show an 11.6% increase in R²compared to methods without VIs, and a 22.2% improvement over methods using only VIs. SHapley Additive exPlanations (SHAP) analysis shows that the volume feature played a key role in model performance and remained stable throughout the training process. Altogether, the proposed approach enhances individual tree biomass and carbon sink estimations, showing great potential for large-scale precise forest carbon monitoring using multi-source remote sensing data.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Natural gas maximal load delivery for multi-contingency analysis

An increasing dependence on natural gas has amplified existing vulnerabilities to the power grid, including disruptions to gas transmission networks from natural and man-made disasters. To address the operational challenges arising from these disruptions, we, in this study, consider the problem of estimating the steady-state operating capacity of a damaged gas pipeline network while ensuring the maximal delivery of load. Specifically, we formulate the mixed-integer nonconvex maximal load delivery (MLD) problem, which proves difficult to solve on large-scale networks. To address this challenge, we present a relaxation of the MLD problem and use it to determine bounds on the transport capacity of a gas pipeline system. A rigorous computational evaluation over network models ranging in size from 11 to 4,197 junctions shows that the relaxation-based method is suitable for analyzing the impacts of multi-contingency network disruptions, often converging to the optimal solution of the relaxation in less than ten seconds.

03 NATURAL GAS↗
DOE OSTI2020

An interface-aware sub-scale dynamics multi-material cell model for solids with void closure and opening at all speeds

Here, we present a multi-material cell model (closure model) for demanding arbitrary Lagrangian-Eulerian (ALE) simulations of fluids and solids. It is based on the interface-aware sub-scale dynamics (IASSD) approach which utilizes the exact material interface geometry within the computational cell to calculate internal material interactions. Our formulation of the closure model also aims to improve the accuracy in low-speed impact events. Voids are used to represent ambient vacuum and internal free boundaries of the distinct materials. Void regions can close and open at contact surfaces, allowing a transition from contact physics to free motion in vacuum. The coupling of void closure and opening with a new formulation of the IASSD model for solids is tested on several one- and two-dimensional numerical examples, ranging from gas expansion in vacuum to planar and round object impacts at various speeds.

42 ENGINEERING↗
DOE OSTI2025

Direct Discontinuous Galerkin methods for the reacting multi-component flow equations

The Direct Discontinuous Galerkin (DDG (Liu and Yan, 2008)) method and a counterpart with Interface Correction (DDGIC (Danis and Yan, 2022)) are extended to compute diffusion terms that arise when solving the compressible multi-component flow equations in thermochemical nonequilibrium. Thermodynamic properties, transport properties, chemical reaction rates, and energy exchange terms are computed using Mutation++ (Scoggins et al., 2020). The DG method is applied on unstructured grids, where the accuracy and convergence rates can be sensitive to the numerical method chosen for parabolic terms. A method for determining the homogeneity tensor of the flow equations required for DDGIC is shown. The convergence properties of the DDG methods are studied and compared to the Interior Penalty (IP) method. A number of numerical experiments are conducted to assess the accuracy and performance of the method. The numerical results and convergence studies indicate that DDG and DDGIC provide accurate solutions and perform well for general flows in thermochemical nonequilibrium.

Diffusion↗
DOE OSTI2021

Caravel: A C++ framework for the computation of multi-loop amplitudes with numerical unitarity

We present the first public version of Caravel, a C++17 framework for the computation of multi-loop scattering amplitudes in quantum field theory, based on the numerical unitarity method. Caravel is composed of modules for the D-dimensional decomposition of integrands of scattering amplitudes into master and surface terms, the computation of tree-level amplitudes in floating point or finite-field arithmetic, the numerical computation of one- and two-loop amplitudes in QCD and Einstein gravity, and functional reconstruction tools. Here, we provide programs that showcase Caravel's main functionalities and allow to compute selected one- and two-loop amplitudes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

$\mathrm{DRAGON}$: A multi-GPU orbital-free density functional theory molecular dynamics simulation package for modeling of warm dense matter

As progress in electronic structure theoretical methods is made, ab initio molecular dynamics (MD) based on orbital-free density functional theory (OF-DFT) is becoming increasingly more successful at substituting the traditional, very accurate but computationally costly Kohn–Sham (KS) approach for simulations of matter at the challenging warm dense matter (WDM) regime. However, despite the significant cost alleviation of eliminating the dependence on the KS orbitals, OF-DFT MD runs require ~10 2 to 10 3 CPU cores running for days, or even weeks, for simulations of systems comprised of 10 2 to 10 3 atoms, depending on thermodynamic conditions. Here we present DRAGON, a multi-GPU OF-DFT MD code for fast and efficient simulations of WDM. With a relatively small allocation of resources (4 to 8 GPU devices) it can provide an order of magnitude speedup for simulations containing $\mathscr{O}$(10 4 ) atoms and target systems composed of $\mathscr{O}$(10 5 ) atoms at conditions within the WDM regime, which is currently outside the capabilities of CPU codes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Towards automated and real-time multi-object detection of anguilliform fishes from sonar data using YOLOv8 deep learning algorithm

Eels (Anguilla spp.), including American eels (Anguilla rostrata), European eels (Anguilla anguilla), and Japanese eels (Anguilla japonica), are species of critical management and regulatory concern due to their vulnerability to various stressors during downstream migrations. Accurate and efficient detection of migrating eels can improve our understanding of fish behaviors and fish-hydraulic structure interactions, thereby facilitating the design, operation, and optimization of more effective downstream passage facilities from both biological and economic perspectives. However, a real-time, automated framework for detecting migrating eels in real-world applications is currently lacking. Leveraging imaging sonar as a reliable technology for fish passage monitoring, field data are acquired using imaging sonar and then converted to single sonar frames/images for subsequent analysis. In this study, a framework based on the You Only Look Once Version 8 (YOLOv8)-based convolutional neural network is proposed for multi-object detection of eels and non-eel fish using the sonar images after image subtraction and additional wavelet denoising. The results from both training and testing phases demonstrate that the framework's ability can successfully detect both eels and non-eel fish in preprocessed sonar images, achieving F1-scores and mAP@0.50 exceeding 0.84. Additionally, the incorporation of wavelet denoising during preprocessing slightly improve detection performance. Furthermore, the transferability of this framework from eel to lamprey detection is demonstrated to be feasible given the similar morphological characteristics of these two species. Overall, the proposed framework achieves accurate and efficient detection of migrating eels, providing reliable and real-time information that can help conserve vulnerable eel and eel-like populations.

Deep learning↗
DOE OSTI2022

A multi-criteria land suitability assessment of field allocation decisions for switchgrass

With the goal of evaluating biomass resource availability in certain regions in the United States, the Department of Energy, Department of Agriculture, Environmental Protection Agency, national laboratories, and Forest Service research laboratories, together with academic and industry collaborators, performed a resource analysis study, 2016 Billion-Ton Report: Advancing Domestic Resources for a Thriving Bioeconomy, Volume 1: Economic Availability of Feedstocks, using POLYSYS model to estimate biomass potential. The study provided county- and state-level biomass potential projections, with biomass feedstock production estimates and suitability assessment of lands acreage likely to transition to energy crops. For this analysis, we estimate biomass supply by year 2028, more specifically, lands acreage transition to switchgrass in the Midwest of the United States, across 50 counties in 3 states (Nebraska, Kansas, and Colorado). We propose a multi-criteria site suitability framework, to determine the potential acreage and production for switchgrass, at the field level. Site selection criteria span agronomic and ecosystem services such as previous land cover, crop rotation, soil water holding capacity, crop productivity, field efficiency, topography, soil organic matter, soil leaching propensity, distance to surface water bodies, and soil erosivity. Efficient agricultural field selection for dedicated energy crops field allocation is critical to evaluate biomass resource availability, and it is our position that estimation at the field level would provide a more accurate and useful estimation to support an emerging bioeconomy. Using field suitability scores and county-level assessments, we model land allocation decision for energy crop in the future, estimating the spatial distribution of dynamics of energy crop adoption with priority to fields scoring the highest. Additionally, compared with POLYSYS results from the report mentioned above, our method shows good performance overall, with about 2% and 1.2% of difference in switchgrass acreage allocation and crop production, respectively.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

New insights into the flow dynamics of a deep freshwater aquifer in the semi-arid and saline Cuvelai-Etosha Basin, Northern Namibia: Results of a multi-environmental tracer study

Study region A paleo-megafan system of the Cubango River in the northern parts of the semi-arid Cuvelai-Etosha Basin, shared by Angola and Namibia. It hosts a deep freshwater aquifer, the so-called Kalahari-Ohangwena 2 (KOH-2), with the potential to resolve the imminent regional water supply shortages. Study focus Hydrogeochemical and multi-environmental tracer studies incorporating the use of age tracers 14 C, 36 Cl, 81 Kr and 4 He to determine the age of groundwater and provide insights into the flow dynamics of the KOH-2. New hydrological insights for the region Stable water isotopes and noble gas thermometry show that in a period with higher rainfall and recharge, temperatures were at least 3 – 4 °C lower than today. Several arguments led to the conclusion that younger groundwater, possibly of an age of 35,000 years, is mixed with ancient saline pore water. These include: 1) the correlation of measured 36 Cl and 81 Kr ratios, as well as 4 He concentrations, using a binary mixing model, and 2) the substantial variation in 81 Kr ages, ranging from 40,000 to 170,000 years, over relatively short distances—a phenomenon challenging to explain by advective groundwater flow equations. Consequently, the ages derived from 81 Kr measurements serve as indicators of the extent of freshening and therefore describe mixing ages rather than absolute travel times.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Localized Corrosion Behavior of Non-Equiatomic NiFeCrMnCo Multi-Principal Element Alloys

A series of single phase Ni 38 Fe 20 Cr x Mn 21–0.5x Co 21–0.5x multi-principal element alloys (MPEAs) were fabricated and their localized corrosion behavior was tested in 0.6 M NaCl solution. These MPEAs were found to be susceptible to localized corrosion, exhibiting pitting at room temperature, as indicated by positive hysteresis loops and low repassivation potentials observed in cyclic potentiodynamic polarization experiments. The pitting resistance increased with increasing Cr content. Interestingly, even the MPEAs with low chromium content, Ni 38 Fe 20 Cr 6 Mn 18 Co 18 and Ni 38 Fe 20 Cr 10 Mn 16 Co 16 , exhibited passivity and localized corrosion. Here, the passivity observed in these low chromium MPEAs is likely derived from their relatively high nickel content. Crystallographic pit morphology was observed in most of these MPEAs, suggesting the absence of a salt film on the dissolving pit surface. Such crystallographic pit morphology indicates that the pit growth in most of these MPEAs was under charge-transfer/ohmic control.

36 MATERIALS SCIENCE↗