Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FACTS”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Pressure-stabilized fixed-stress iterative solutions of compositional poromechanics

We consider the numerical behavior of the fixed-stress splitting method for coupled poromechanics as undrained regimes are approached. We explain that pressure stability is related to the splitting error of the scheme, not the fact that the discrete saddle point matrix never appears in the fixed-stress approach. This observation reconciles previous results regarding the pressure stability of the splitting method. Using examples of compositional poromechanics with application to geological CO sequestration, we see that solutions obtained using the fixed-stress scheme with a low order finite element-finite volume discretization which is not inherently inf-sup stable can exhibit the same pressure oscillations obtained with the corresponding fully implicit scheme. Moreover, pressure jump stabilization can effectively remove these spurious oscillations in the fixed-stress setting, while also improving the efficiency of the scheme in terms of the number of iterations required at every time step to reach convergence.

42 ENGINEERING↗
DOE OSTI2024

Equivariant graph convolutional neural networks for the representation of homogenized anisotropic microstructural mechanical response

Composite materials with different microstructural material symmetries are common in engineering applications where grain structure, alloying and particle/fiber packing are optimized via controlled manufacturing. In fact these microstructural tunings can be done throughout a part to achieve functional gradation and optimization at a structural level. To predict the performance of particular microstructural configuration and thereby overall performance, constitutive models of materials with microstructure are needed. In this work we provide neural network architectures that provide effective homogenization models of materials with anisotropic components. These models satisfy equivariance and material symmetry principles inherently through a combination of equivariant and tensor basis operations. We demonstrate them on datasets of stochastic volume elements with different textures and phases where the material undergoes elastic and plastic deformation, and show that the these network architectures provide significant performance improvements.

anisotropy↗
DOE OSTI2025

An end-to-end deep learning method for solving nonlocal Allen–Cahn and Cahn–Hilliard phase-field models

Here, we propose an efficient end-to-end deep learning method for solving nonlocal Allen–Cahn (AC) and Cahn–Hilliard (CH) phase-field models. One motivation for this effort emanates from the fact that discretized partial differential equation-based AC or CH phase-field models result in diffuse interfaces between phases, with the only recourse for remediation is to severely refine the spatial grids in the vicinity of the true moving sharp interface whose width is determined by a grid-independent parameter that is substantially larger than the local grid size. In this work, we introduce non-mass conserving nonlocal AC or CH phase-field models with regular, logarithmic, or obstacle double-well potentials. Because of non-locality, some of these models feature totally sharp interfaces separating phases. The discretization of such models can lead to a transition between phases whose width is only a single grid cell wide. Another motivation is to use deep learning approaches to ameliorate the otherwise high cost of solving discretized nonlocal phase-field models. To this end, loss functions of the customized neural networks are defined using the residual of the fully discrete approximations of the AC or CH models, which results from applying a Fourier collocation method and a temporal semi-implicit approximation. To address the long-range interactions in the models, we tailor the architecture of the neural network by incorporating a nonlocal kernel as an input channel to the neural network model. We then provide the results of extensive computational experiments to illustrate the accuracy, predictive capabilities, and cost reductions of the proposed method.

42 ENGINEERING↗
DOE OSTI2021

Platinum group metal-free Fe-based (Fe-N-C) oxygen reduction electrocatalysts for direct alcohol fuel cells

In direct alcohol fuel cells (DAFCs), the oxidation of alcohols happens at the expenses of the oxygen reduction reaction at the cathode. DAFCs' cathodes use a significant amount of platinum that is an expensive critical raw material. Moreover, platinum oxidizes alcohols, a fact that combined with alcohol crossover, decreases significantly the performance of the cells. The use of Fe-based (Fe-N-C) platinum group metal–free (PGM-free) cathodes is a convenient strategy to overcome these limitations. This review analyzes the application of PGM-free cathodes to DAFCs. The discussion focuses on acidic systems and covers the following subjects: (i) the breakdown of DAFC potential in its components, (ii) the analysis of the advantages from the use of the PGM-free cathode, and (iii) a review of the performance and durability of DAFCs. Finally, the review closes with a view of the authors of the future perspective for the research.

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

Theoretical kinetics predictions for NH 2 + HO 2

Recent modeling studies of NH 3 oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies between existing literature mechanisms. In this study high level theoretical kinetics predictions have been obtained for the reaction of NH 2 with HO 2 , which has previously been highlighted as an important reaction with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations including large basis sets and high-level corrections to yield high accuracy (~0.2 kcal/mol) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature and pressure dependent kinetics. For the radical-radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constant as well as the branching to the NH 3 + O 2 , H 2 NO + OH, and HNO + H 2 O channels. Rate constants are also reported for the H 2 NO + OH reaction as they arise naturally from the analysis. The rate constant and branching fraction determined in this work for the NH 2 + HO 2 reaction deviate significantly from values used in most previous modeling studies. The fact that the main product channel is chain terminating, rather than propagating, has strong implications for modeling NH 3 ignition and oxidation, in particular at intermediate temperatures and elevated pressure.

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

Co-optimized machine-learned manifold models for large eddy simulation of turbulent combustion

Many modeling approaches in large eddy simulation (LES) of turbulent combustion employ a projection of the thermochemical state onto a low-dimensional manifold within state space to reduce the number of transported variables and hence computational cost. Flamelet-generated manifolds (FGM) is an example of a well-established, physics-based approach, but increasingly, principal component analysis (PCA) is being used as a data-driven method for generating manifold models. For both approaches, the nonlinear relationship between the location on the predefined manifold and the outputs of interest, such as reaction rates, can be tabulated or encoded in a neural network. This work proposes a new approach for manifold modeling that extends these existing approaches. A modified neural network structure simultaneously encodes the definition of the manifold variables, the nonlinear mapping, and the subfilter closure for LES. This allows all three of these aspects of the model to be co-optimized, generating a model from any source of combustion thermochemical state data. The manifold parameterizing variables are constrained to be linear combinations of species, as in FGM and PCA-based models, to aid in interpretability and implementation. For LES, subfilter variances of the manifold variables are also included as inputs. Two types of a priori analysis are performed to evaluate the new approach. In the first, the model is trained on data from one-dimensional premixed flames. In this case, the approach recovers the behavior of flamelet-based manifold approaches, and in fact slightly improves performance by identifying an optimized progress variable. The approach is also applied to data from direct numerical simulations of spherical ignition kernels in isotropic turbulence. For any specified manifold dimensionality, the new approach provides substantially lower prediction errors than a PCA-based model developed from the same data set. Additionally, the LES formulation of the new approach can provide accurate predictions for filtered reaction rates across a variety of filter widths.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Direct numerical simulation of a turbulent methane/air flame impinged by a sub-breakdown electric field

Here, the effects of an external electric field on a turbulent methane/air diffusion flame are analyzed in this work using direct numerical simulations. The analyzed configuration consists of a temporally evolving mixing layer of air and a mixture of methane and nitrogen at 1 atm that is impinged by an electric field in the direction normal to the mean mixing plane. The combustion and chemi-ionization reactions involved in the flow are modeled using finite-rate chemistry and a reduced reaction scheme consisting of 26 species and 134 reactions. The mass diffusion and ion-wind effects are modeled using a detailed description of the diffusion coefficients and electric mobilities based on kinetic theory. The presented calculations show that the turbulence generated within the mixing layer is mostly unaffected by the applied electric field for the configuration under exam. In fact, the electric body force is developed away from the mixing region where the flow is uniform. Conversely, the turbulence is able to introduce very high intermittency in the electrically charged species concentration and, consequently, in the electric body force. Such intermittency will constitute a challenge in the future formulation of combustion models that take into account ion-wind effects.

09 BIOMASS FUELS↗
DOE OSTI2022

A phase-field model for void and gas bubble superlattice formation in irradiated solids

A phase-field model to simulate the formation of both void and gas bubble superlattices is derived from a grand potential functional, assuming 1D diffusion of self-interstitial atoms. The model is capable of accounting for superlattice formation by either a nucleation and growth or spinodal decomposition mechanism; in this work, we focus on the nucleation and growth mechanism, using a discrete nucleation approach in the phase-field model. In simulations of void formation, short aligned rows of voids were initially formed, followed by growth in the size of the aligned rows, and finally leading to superlattice formation, consistent with experimental observations. Void superlattice spacing as a function of model parameters was quantified. Increased nucleation rates led to smaller superlattice spacing, while increased vacancy diffusivity led to larger superlattice spacing. In simulations of gas bubble superlattice formation, superlattice spacing decreased with increasing gas atom flux, consistent with experimental observations. Here, simulations with increasing ratio of gas atom flux to vacancy-interstitial pair production led to increased superlattice spacing, which is inconsistent with experimental observations at high gas atom flux rates; this discrepancy is attributed to the fact that the model does not account for bound gas atom-vacancy pairs that may dominate at high gas flux rates.

36 MATERIALS SCIENCE↗
DOE OSTI2025

On the numerical sensitivity of cellular automata grain structure predictions to large thermal gradients and cooling rates

Cellular automata (CA) models of as-solidified grain structure, originally developed and applied to casting, have become a common means of predicting grain structure resulting from Additive Manufacturing (AM) processes. The majority of these models are based on the decentered octahedron approach, which attempts to correct for the effect of grid anisotropy on the prediction of competitive solidification of dendritic grains. However, AM solidification occurs under cooling rates ($\dot{T}$) and thermal gradients (G) that are orders of magnitude larger than those encountered in casting, and no systematic investigation on the effect of the CA model cell size (Δx) and time step (Δt) on AM microstructure predictions has been performed. Here, in this study, such an investigation is first performed via simulation of individual grains of various crystallographic orientations with a fixed, unidirectional G, showing that CA prediction of the steady-state undercooling matched the expected values based on the interfacial response function at small G and deviated from the expected values at large G. Simulation of competitive growth of multiple grains showed a weakening of the predicted texture as G and Δx became large. Simulation of solidification under AM conditions, where G and $\dot{T}$ vary spatially across the melt pools, showed that not only does grain selection weaken and deviate from expectations at large Δx, but grains with crystallographic $\langle$100$\rangle$ aligned with the grid directions are more adversely affected by the temperature field discontinuities than grains with other crystallographic orientations. Despite the fact that the exact grain competition results depended on Δt, the overall texture development was notably less sensitive to Δt than Δx, provided that a reasonable value of Δt is selected based on the ratio of Δx to the maximum local solidification velocity in the simulation domain. Finally, from the directional solidification and AM simulation results, an analysis of computational cost compared to simulation resolution is performed based on an equation derived to quantify the relatively inaccuracy in grain selection based on the model and temperature field inputs. From this analysis, it is concluded that there is a need for algorithmic improvements to improve CA grain competition accuracy for large G processing conditions as sufficiently small Δx to resolve the necessary competition is intractable for many AM processing conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine learning for domain transfer between simulated and experimental 2D X-ray diffraction patterns using generative adversarial networks

X-ray diffraction (XRD) is a well-established technique for analyzing materials at an atomic level. Dynamic compression experiments (DCE), in which materials are subject to extreme pressures, can provide fundamental understanding to pressure-induced phase transitions and compression of the crystal lattice. The analysis of XRD patterns from highly compressed samples is non-trivial given the sparsity of data, high experimental costs, and the fact that the data is often marred with X-ray background and other artifacts. While accurate computational frameworks exist, they solve the forward problem—from structures and orientations to XRD patterns. Solving the inverse problem for 2D experimental diffraction patterns is currently a complex manual process of matching and comparing experimentally observed patterns to computationally generated ones. Machine learning is a promising tool for automating the matching process but often requires data-intensive architectures. Here, in this study, we use a CycleGAN to translate the domain of limited experimental data to a domain in which there is readily available simulated data. This domain shift allows data-intensive machine learning models that have only been trained on simulated XRD patterns to be used in the analysis of experiments.

Brozak, Samantha Jean [Sandia National Laboratorie↗
DOE OSTI2021

A nudged hybrid analysis and modeling approach for realtime wake-vortex transport and decay prediction

Here, we put forth a long short-term memory (LSTM) nudging framework for the enhancement of reduced order models (ROMs) of fluid flows utilizing noisy measurements for air traffic improvements. Toward emerging applications of digital twins in aviation, the proposed approach allows for constructing a realtime predictive tool for wake-vortex transport and decay systems. We build on the fact that in realistic application, there are uncertainties in initial and boundary conditions, model parameters, as well as measurements. Moreover, conventional nonlinear ROMs based on Galerkin projection (GROMs) suffer from imperfection and solution instabilities, especially for advection-dominated flows with slow decay in the Kolmogorov n-width. In the presented LSTM nudging (LSTM-N) approach, we fuse forecasts from a combination of imperfect GROM and uncertain state estimates, with sparse Eulerian sensor measurements to provide more reliable predictions in a dynamical data assimilation framework. We illustrate our concept by solving the two-dimensional vorticity transport equation. We investigate the effects of measurements noise and state estimate uncertainty on the performance of the LSTM-N behavior. We also demonstrate that it can sufficiently handle different levels of temporal and spatial measurement sparsity, and offer a huge potential in developing next-generation digital twin technologies for aerospace applications.

42 ENGINEERING↗
DOE OSTI2023

Moving-horizon false data injection attack design against cyber–physical systems

Systematic attack design is essential to understanding the vulnerabilities of cyber–physical systems (CPSs), to better design for resiliency. In particular, false data injection attacks (FDIAs) are well-known and have been shown to be capable of bypassing bad data detection (BDD) while causing targeted biases in resulting state estimates. However, their effectiveness against moving horizon estimators (MHE) is not well understood. In fact, this paper shows that conventional FDIAs are generally ineffective against MHE. One of the main reasons is that the moving window renders the static FDIA recursively infeasible. Here, this paper proposes a new attack methodology, moving-horizon FDIA (MH-FDIA), by considering both the performance of historical attacks and the current system’s status. Theoretical guarantees for successful attack generation and recursive feasibility are given. Numerical simulations on the IEEE-14 bus system further validate the theoretical claims and show that the proposed MH-FDIA outperforms state-of-the-art counterparts in both stealthiness and effectiveness. In addition, an experiment on a path-tracking control system of an autonomous vehicle shows the feasibility of the MH-FDIA in real-world nonlinear systems.

42 ENGINEERING↗
DOE OSTI2021

Modeling of High-Temperature Corrosion of Zirconium Alloys Using the eXtended Finite Element Method (X-FEM)

Oxidation modeling in modern nuclear fuel performance codes is currently limited by the lack of coupling with mechanics, thus preventing proper description of how high-temperature oxidation impacts mechanical properties. This is mostly due to the fact that the finite difference formalism adopted in corrosion models is incompatible with the direct coupling with mechanics in the finite element modeling employed in modern nuclear fuel performance codes. In this study, a physically based zirconium alloy corrosion model called the Coupled-Current Charge Compensation (C4) model, which was initially developed for operating temperature conditions, has been updated to include high-temperature corrosion in order to provide additional critical information (e.g., oxygen concentration profile) under loss-of-coolant accident (LOCA) conditions—information lacking in existing empirical models. The C4 model was implemented in the MOOSE finite-element framework developed at Idaho National Laboratory, enabling it to be used in the BISON nuclear fuel performance code based on the MOOSE framework. To precisely track the different interfaces at a relatively low computational cost, the eXtended Finite Element Method (X-FEM) was applied in MOOSE. The model’s results were compared to those of existing empirical models as well as metallographic analysis of high-temperature oxidized Zircaloy-4 coupons. Oxygen diffusivities in the a and ß phases resulting from this comparison closely agree with those found in the literature. The C4 model implemented with X-FEM in MOOSE now has the capability to accurately predict oxide, oxygen-stabilized a, and prior ß phase layer growth kinetics under isothermal exposure at high temperature (1000–1500°C). Furthermore, in contrast with the empirical models, the C4 model accounts for the finite thickness of the fuel cladding. It can predict the oxygen concentration profile evolution through the whole cladding, enabling evaluation of the remaining ductile thickness—a crucial variable for modeling the mechanical behavior of the fuel cladding under LOCA. Furthermore, this implementation allows direct coupling with mechanics, at a low computing cost, using finite-element-based nuclear fuel performance codes such as BISON.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A preference for cold dark matter over Superfluid Dark Matter in local Milky Way data

There are many well-known correlations between dark matter and baryons that exist on galactic scales. These correlations can essentially be encompassed by a simple scaling relation between observed and baryonic accelerations, historically known as the Mass Discrepancy Acceleration Relation (MDAR). The existence of such a relation has prompted many theories that attempt to explain the correlations by invoking additional fundamental forces on baryons. The standard lore has been that a theory that reduces to the MDAR on galaxy scales but behaves like cold dark matter (CDM) on larger scales provides an excellent fit to data, since CDM is desirable on scales of clusters and above. However, this statement should be revised in light of recent results showing that a fundamental force that reproduces the MDAR is challenged by local Milky Way dynamics and rotation curve data between 5–18 kpc. In this study, we test this claim on the example of Superfluid Dark Matter. We find that a standard CDM model is preferred over a static superfluid profile assuming a steady-state Galactic disk and discuss the robustness of this conclusion to disequilibrium effects. This preference is due to the fact that the superfluid model over-predicts vertical accelerations, even while reproducing galactic rotation curves. Our results establish an important criterion that any dark matter model must satisfy within the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Optimizing PCM-integrated walls for potential energy savings in U.S. Buildings

Buildings in the United States account for nearly half of total U.S. energy use. The energy used for space conditioning can be reduced by utilizing thermal energy storage, such as phase change materials (PCMs), into building envelopes; however, the energy savings of PCM-integrated building envelopes reported in the literature vary widely. In the absence of established guidelines, thermophysical requirements of an optimal PCM, its method of application into the building envelope, and the corresponding energy savings under various climates remain unknown. Here, we perform an extensive numerical investigation on the integration of PCM into building walls to establish the key conditions required for effective utilization of PCM in reducing heat gains in the cooling season and heat losses in the heating season. We also determine the optimal transition temperature, optimal PCM location in the wall, and the energy-saving potential of the PCM-integrated building walls in five U.S. cities located in different International Energy Conservation Code climate zones. Results show that employing PCMs in building walls does not always lead to an improvement; in fact, incorrect applications of PCMs can substantially increase energy use in the buildings. In the climates we studied, PCMs were found effective in reducing heat gains during the cooling season while mostly ineffective in managing heat losses during heating season. Depending on the climate, optimized PCMs in U.S. building walls can provide reduction in the annual heat gain in the range of 3.5% to 47.2% and the annual heat loss in the range of -2.8% to 8.3%. Future consideration of buildings with substantial solar gains in winter may lead to more reduction in heat losses by PCMs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Energy impact of human health and wellness lighting recommendations for office and classroom applications

The goal of this investigation was to evaluate potential energy impacts of circadian lighting design recommendations gaining attention in a variety of common applications such as offices and classrooms. The renewed focus on health along with advances in SSL technology capabilities has underscored that there is still much to learn regarding the relationship between light and human physiology. The energy implications of designing to address these possible physiological effects are not yet fully understood. Beyond the fact that the basic metric of luminous efficacy (lumens per watt) does not cover these other effects, the emerging science seems to indicate that addressing a holistic view of the human needs in most applications may mean a need for increased light and associated energy use by electric lighting systems. Within the two applications considered, lumen output, spectral characteristics, surface reflectance distribution and desk orientation were varied to explore the magnitude of potential effects. Meeting current IES illuminance recommendations did not satisfy existing EML and CS recommendations for any of the simulations. In some cases, satisfying circadian metric recommendations required average illuminance that was more than double IES recommendations, which may negatively impact lighting quality as well. Using results from 45 unique simulation conditions, it was estimated that lighting energy use may increase between 10% and 100% due to increased luminaire light levels used to meet circadian lighting design recommendations listed in current building standards such as WELL v2 Q2 2019, UL Design Guideline 24480, and CHPS Core Criteria 3.0.

circadian lighting, lighting simulation, solid sta↗
DOE OSTI2022

Assessment of lignite upgrade and hydrogen evolution via electrolysis

Lignite, considered as the lowest-rank coal, accounts for 45% of the total coal resources worldwide, and its low carbon and high ash contents limit its large-scale application. In order to upgrade lignite effectively, a modified coal electrolytic cell (MCEC) was designed, built, and then evaluated for lignite upgrade and hydrogen production through electrolysis, and the operating conditions were optimized. The results show that the performance of the MCEC outperformed the previous coal electrolytic cell (CEC). The operating temperature and iron concentration significantly impact the average current density per milligram of catalyst (ACDMC). The content of dry fixed carbon and heating value increased by 20.9% (wt%) and 8.15%, respectively, in electrolyzed lignite compared with raw lignite. Besides, 71.5% ash and 39.6% mercury were removed from raw lignite after electrolysis, which is the first report. The Faradaic efficiency (FE) for hydrogen evolution reached 99% using MCEC. In addition, the analysis of spectroscopy methods (SEM, EDX, FTIR, and XRD) revealed that the lignite particles were cracked; the oxygen-containing functional groups increased; the minerals of quartz and hematite formed on the surface of lignite particles after electrolysis. Furthermore, the facts indicate that electrolysis of lignite is a competitive approach for lignite upgrade and hydrogen production.

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

Mapping thermal energy storage technologies with advanced nuclear reactors

Advanced nuclear power plants (NPPs) will potentially need to operate in environments where power generation flexibility is more highly valued than the stability or baseload generation capability for conventional demand curves. Thermal energy storage (TES) systems would enable NPPs to respond nimbly to market variability and could also position advanced NPPs to participate differently in restructured markets, thus further enhancing their economic competitiveness. TES systems could also benefit the electric grid by eliminating the need for peaking plants, as well as by improving the economic performance of baseload NPPs. While TES technologies afford a unique opportunity to address many of these challenges, the applicability of these systems is also complicated by the fact that various advanced NPPs are designed differently, each with its own temperature range, size, operating fluids, and operating conditions. Hence, TES systems face significant barriers to investment, as more information on their compatibility and performance metrics is needed to quantify the advantages provided by each, as well as the challenges these technologies might face if coupled with a particular type of advanced NPP. This study explores the possibility of integrating a wide variety of TES technologies with various categories of advanced NPPs, based on their operating characteristics. To help decision makers, users and developers decide which TES technology is best suited to a particular category of advanced NPPs, this research present a Phenomena Identification and Ranking Table (PIRT) analysis of 10 TES systems that could potentially be coupled with advanced NPPs, which themselves are divided into nine categories based on their operating conditions. Then, each advanced NPP category is evaluated for compatibility with the 10 TES systems by assembling and discussing a database of information concerning 10 engineering questions, defined herein in as figures of merit (FOMs), such as: technology readiness level (TRL), temperature compatibility, energy density, size, cycle frequency, ramp time, realignment frequency, geographic needs, environmental impact, and interventions. By assembling a database of information concerning the TES technologies' compatibility with various advanced NPP systems, this study can help developers acquaint themselves with a particular TES technology before choosing to build a new integrated installation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗