Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical stability”

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

Understanding the effect of density functional choice and van der Waals treatment on predicting the binding configuration, loading, and stability of amine-grafted metal organic frameworks

Metal organic frameworks (MOFs) are crystalline, three-dimensional structures with high surface areas and tunable porosities. Made from metal nodes connected by organic linkers, the exact properties of a given MOF are determined by node and linker choice. MOFs hold promise for numerous applications, including gas capture and storage. M2(4,4′-dioxidobiphenyl-3,3′-dicarboxylate)—henceforth simply M2(dobpdc), with M = Mg, Mn, Fe, Co, Ni, Cu, or Zn—is regarded as one of the most promising structures for CO2 capture applications. Further modification of the MOF with diamines or tetramines can significantly boost gas species selectivity, a necessity for the ultra-dilute CO2 concentrations in the direct-air capture of CO2. There are countless potential diamines and tetramines, paving the way for a vast number of potential sorbents to be probed for CO2 adsorption properties. The number of amines and their configuration in the MOF pore are key drivers of CO2 adsorption capacity and kinetics, and so a validation of computational prediction of these quantities is required to suitably use computational methods in the discovery and screening of amine-functionalized sorbents. In this work, we study the predictive accuracy of density functional theory and related calculations on amine loading and configuration for one diamine and two tetramines. In particular, we explore the Perdew–Burke–Ernzerhof (PBE) functional and its formulation for solids (PBEsol) with and without the Grimme-D2 and Grimme-D3 pairwise corrections (PBE+D2/3 and PBEsol+D2/3), two revised PBE functionals with the Grimme-D2 and Grimme-D3 pairwise corrections (RPBE+D2/3 and revPBE+D2/3), and the nonlocal van der Waals correlation (vdW-DF2) functional. We also investigate a universal graph deep learning interatomic potential’s (M3GNet) predictive accuracy for loading and configuration. These results allow us to identify a useful screening procedure for configuration prediction that has a coarse component for quick evaluation and a higher accuracy component for detailed analysis. Our general observation is that the neural network-based potential can be used as a high-level and rapid screening tool, whereas PBEsol+D3 gives a completely qualitatively predictive picture across all systems studied, and can thus be used for high accuracy motif predictions. We close by briefly exploring the predictions of relative thermal stability for the different functionals and dispersion corrections.

Chemistry↗
DOE OSTI2024

Low-Level Liquid-Bearing Clouds Contribute to Seasonal Lower Atmosphere Stability and Surface Energy Forcing over a High-Mountain Watershed Environment

Measurements of atmospheric structure and surface energy budgets distributed along a high-altitude moun tain watershed environment near Crested Butte, Colorado, from two separate, but coordinated, field campaigns, Surface Atmosphere Integrated field Laboratory (SAIL) and Study of Precipitation, the Lower Atmosphere, and Surface for Hydrometeorology (SPLASH), are analyzed. This study identifies similarities and differences in how clouds influence the radiative budget over one snow-free summer season (2022) and two snow-covered seasons (2021/22; 2022/23) for this alpine location. A relationship between lower-tropospheric stability stratification and longwave radiative flux from the presence or absence of clouds is identified. When low clouds persisted, often with signatures of supercooled liquid in win ter, the lower troposphere experienced weaker stability, while radiatively clear skies that are less likely to be influenced by liquid droplets were associated with appreciably stronger lower-tropospheric stratification. Corresponding surface turbu lent heat fluxes partitioned differently based upon the cloud–stability stratification regime derived from early morning radiosounding profiles. Combined with the differences in the radiative budget largely resulting from dramatic seasonal dif ferences in surface albedo, the lower atmosphere stratification, surface energy budget, and near-surface thermodynamics are shown to be modified by the effective longwave radiative forcing of clouds. The diurnal evolution of thermodynamics and surface energy components varied depending on the early morning stratification state. Thus, the importance of quies cent versus synoptically active large-scale meteorology is hypothesized as a critical forcing for cloud properties and associ ated surface energy budget variations. The physical relationships between clouds, radiation, and stratification can provide a useful suite of metrics for process understanding and to evaluate numerical models in such an undersampled, highly com plex terrain environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Quantum stabilization of unexpected ordered phases on the honeycomb lattice

In this work, I discuss and showcase the utility of the method called minimally-augmented spin-wave theory (MAGSWT) as a relatively simple semi-analytical approach to the phase diagrams of quantum magnets. It complements numerical methods by providing physical insight into which states are competitive and by yielding approximate phase boundaries that agree well with much more numerically intensive calculations. This approach enabled me to construct the faithful phase diagrams of two paradigmatic honeycomb-lattice models in the quantum S=1/2 limit, the J 1 –J 3 FM-AF and J 1 –J 2 AF models, for the collinear quantum phases that replace or extend the classical ones. The results are in good qualitative and semi-quantitative agreement with state-of-the-art numerical studies for these models, correctly capturing the emergence of unexpected quantum phases and the suppression of classically favored spiral orders by quantum fluctuations. This study provides a much-needed important guidance to the ongoing theoretical and experimental searches of the unconventional quantum states. This work will be of significant and timely interest to both theorists and experimentalists in the field of quantum magnetism, broadly defined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A discontinuous piecewise polynomial generalized moving least squares scheme for robust finite element analysis on arbitrary grids

A variational approach is developed with a meshless discretization to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. Here, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality, which we call discontinuous piecewise polynomial generalized moving least squares (DPP-GMLS). At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework, providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5 x improvement in accuracy and allows for taking an 8 x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Instabilities of a vortex-ring-bright soliton in trapped binary three-dimensional Bose-Einstein condensates

Instabilities of vortex-ring-bright coherent structures in harmonically trapped two-component three-dimensional Bose-Einstein condensates are studied numerically within the coupled Gross-Pitaevskii equations and interpreted analytically. Interestingly, the filled vortex core with a sufficiently large amount of the bright component is observed to reduce the parametric interval of stability of the vortex ring. Here, we have identified the mechanisms of several linear instabilities and one nonlinear parametric instability in this connection. Two of the linear instabilities are qualitatively different from ones reported earlier and are associated with azimuthal modes of m = 0 and m = 1, i.e., deviations of the vortex from the stationary ring shape. Our nonlinear parametric resonance instability occurs between the m = 0 and m = 2 modes and signals the exchange of energy between them.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Critical role of magnetic moments in the lattice dynamics of YBa 2 Cu 3 O 6

The role of lattice dynamics in unconventional high-temperature superconductivity is still vigorously debated. Theoretical insights into this problem have long been prevented by the absence of an accurate first-principles description of the combined electronic, magnetic, and lattice degrees of freedom. Utilizing the recently constructed r 2 SCAN density functional that stabilizes the antiferromagnetic (AFM) state of the pristine oxide YBa 2 ⁢Cu 3 ⁢O 6 , we faithfully reproduce the experimental dispersion of key phonon modes. Here we further find significant magnetoelastic coupling in numerous high-energy Cu-O bond stretching optical branches, where the AFM results improve over the soft nonmagnetic phonon bands.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Enhanced radiative and thermal properties from surface encapsulation of InAs nanowires

In this work, photoluminescence (PL), quantum efficiency and carrier dynamics are investigated in indium arsenide (InAs) nanowires (NWs) with various surface treatments, including a molecular beam epitaxy (MBE)-grown semiconductor shell passivation, sulfur-passivation, alumina (Al 2 O 3 ) coating by atomic layer deposition (ALD) and polydimethylsiloxane (PDMS) spin-coating. The ALD-dielectric layer-coated InAs core-shell NWs show a maximum 13-fold increase in PL intensity. In contrast to the previous reports, this enhancement is found to be due to increased radiative rate from an enhanced Purcell factor, better thermal conductance and higher carrier injection within the NWs instead of improved surface quality. Numeric simulations confirm the experimentally observed increased radiative rate. Further improvements are suggested with even thicker capped InAs NWs. Carrier lifetime in surface-treated NWs is extended and shows long-term stability, critical for practical devices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids. LDRD Final Report

This report summarizes the work performed under a one-year LDRD project aiming to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. In this project, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality. At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5x improvement in accuracy and allows for taking an 8x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method. The report concludes with a brief summary of ongoing projects and collaborations that utilize or extend the products of this work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Impact of particle size and particle-flow-wall coupling on pressurized oxy-combustion in the down-fired burner

Concerns over climate change have led to numerous efforts in developing low-carbon energy technologies. Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in power generation. Designing an effective burner plays a vital role in developing new coal combustion technologies. Because of the high pressure, the volume fraction of coal particles at the fuel inlet of the burner of a pressurized oxy-combustor is close or even higher than the maximum limit that commercial CFD codes (e.g., ANSYS FLUENT) can handle. At this high particle volume fraction, the interactions among particles, fluid flow, and wall need to be re-evaluated. The present computational work is the first step of a systematic analysis of the particle influence, like releasing method, releasing locating, and particle size, in a pilot-scale POC combustor, developed at Washington University in St. Louis (WUSTL). In a POC process, pulverized coal is burned under elevated pressure and O2-CO2 environment. Specifically, a 15-bar, 100 kWth, POC combustor is modeled employing ANSYS FLUENT, using Reynolds-averaged Navier-Stokes (RANS) modeling. It is revealed that for this pilot-scale, pressurized burner, velocity profiles of the near-wall region in this POC facility exhibit some discrepancy against the near-wall turbulent flow velocity profile. In order to investigate the particle influence in the near-wall region, particle releasing location will be investigated. Numerical simulation results exhibit the coupling effect of turbulence flow and particles in this case. The particle size also demonstrates a great effect on particle trajectory, then further has an impact on flame stability and temperature.

Li, Lei↗
DOE OSTI2026

Towards High-Throughput Computation of Phase-and Defect Diagrams

The past decade has seen immense advances in our understanding of defect thermodynamics, and the use of machine learning and data science approaches has played a critical role in these advances [1–14]. In the area of grain boundaries (GBs), a particular focus has been placed on the effects of alloying – namely, GB solute segregation or more broadly, GB alloying [15–25], which has been observed and catalogued across a vast range of systems [26–50]. The impacts of solute segregation to GBs are numerous, and can range from negative effects such as embrittlement – for example, due to impurities [51–53], during irradiation [54–61], or during heat treatment [62–65] – to positive effects such as the stabilization against grain growth [66–69], thus enabling the design of nanocrystalline alloys with access to an enhanced range of functional and mechanical properties, and the reduction of embrittlement through the segregation of GB strengthening solutes [49,70–79].

36 MATERIALS SCIENCE↗
DOE OSTI2020

STRUCTURAL STABILITY OF MAGNETIC TUNNEL JUNCTION BASED MOLECULAR SPINTRONICS DEVICES (MTJMSD)

Harnessing the exotic properties of molecular level nanostructures to produce novel sensors, metamaterials, and futuristic computer devices can be technologically transformative. In addition, connecting the molecular nanostructures to ferromagnetic electrodes bring the unprecedented opportunity of making spin property based molecular devices. We have demonstrated that magnetic tunnel junction based molecular spintronics device (MTJMSD) approach to address numerous technological hurdles that have been inhibiting this field for decades. MTJMSD approach is based on producing a capacitor like a testbed where two metal electrodes are separated by an ultrathin insulator and subsequently bridging the molecule nanostructure across the insulator to transform a capacitor into a molecular device. Our prior work showed that MTJMSDs produced extremely intriguing phenomenon such as room temperature current suppression by six orders, spin photovoltaic effect, and evolution of new forms of magnetic metamaterials arising due to the interaction of the magnetic a molecule with two ferromagnetic thin films. However, making robust and reproducible electrical connections with exotic molecules with ferromagnetic electrodes is full of challenges and requires attention to MTJMSD structural stability. This paper focuses on MTJMSD stability by describing the overall fabrication protocol and the associated potential threat to reliability. MTJMSD is based on microfabrication methods such as (a) photolithography for patterning the ferromagnetic electrodes, (b) sputtering of metallic thin films and insulator, and (c) at the end electrochemical process for bridging the molecules between two ferromagnetic films separated by ~ 2nm insulating gap.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Direct Numerical Simulation of Flame-Wall Interaction for Low-Carbon Gas Turbine Combustion

Green hydrogen (H2) and ammonia (NH3) are emerging as carbon free alternatives to hydrocarbon fuels. In gas turbines, flames interact with the combustor liner (wall), which affects the pollutant emissions, the burning efficiency, and the thermal load on the liner. We lack understanding of this important flame-wall interaction (FWI) for alternative fuels. FWI occurs at the scale of the flame thickness, and such scales are only resolved by direct numerical simulation. Here, we resolve the FWI for NH3 and H2 flames. Preliminary results of two distinct projects are presented: 2D laminar anchored V-flame, and 3D turbulent swirling flame. For the laminar flame, we show flame quenching at the wall and resulting leakage of NH3 and increased N2O, but decreased NO emissions. For the swirling flame, we show CO emissions, and comparison of flame shapes and location to experimental measurements. Ultimately we show that quenching and pollutant emissions strongly affected by fuel/air ratio, and that an inclined wall can stabilize very lean turbulent CH4/H2 flames, but prone to large CO emissions.

adaptive mesh refinement↗
DOE OSTI2021

Efficient Topology Design Algorithms for Power Grid Stability

The dynamic response of power grids to small disturbances influences their overall stability. This letter examines the effect of network topology on the linearized time-invariant dynamics of electric power systems. The proposed framework utilizes H 2 -norm based stability metrics to study the optimal placement of lines on existing networks as well as the topology design of new networks. The design task is first posed as an NP-hard mixed-integer nonlinear program (MINLP) that is exactly reformulated as a mixed-integer linear program (MILP) using McCormick linearization. To improve computation time, graph-theoretic properties are exploited to derive valid inequalities (cuts) and tighten bounds on the continuous optimization variables. Moreover, a cutting plane generation procedure is put forth that is able to interject the MILP solver and augment additional constraints to the problem on-the-fly. Finally, the efficacy of our approach in designing optimal grid topologies is demonstrated through numerical tests on the IEEE 39-bus network.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Aerosols consistently suppress the convective boundary layer development

Aerosols with different vertical distribution and various optical properties induce diverse heating rates and thereby affecting convective boundary layer (CBL) development. Our results showed consistent CBL-suppression of aerosols during daytime with numerical experiments, in which aerosols were specified at different heights with synthesized single scattering albedo from 64 studies and asymmetry factor from 20 studies globally. Absorbing aerosols concentrated below but close to the CBL top had the strongest suppression effect on CBL development relative to that concentrated near surface or above CBL. Aerosol cooling effect by attenuating incident solar radiation and surface heat flux exceeded its warming effect by reheating the atmosphere layer with absorbed shortwave radiation, and eventually declined net heating rate, which inhibited CBL development, lowered mixed-layer potential temperature and stabilized atmospheric stratification. Stove effect of absorbing aerosols (CBL enhancement) under a zero background aerosol extinction coefficient is negligible for dominant dome effect (CBL suppression) which consistently suppresses CBL development regardless of aerosol vertical height and background aerosol extinction coefficient. Our study also highlighted the importance of specifying background aerosol extinction coefficient in numerical experiments for accurate assessment of aerosol radiative forcing and CBL-aerosol interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Microstructure and creep properties of cast near-eutectic Al–Ce–Ni alloys

This study investigates the as-cast and aged microstructures, thermal stability, ambient temperature strengthening, and creep resistance of three ternary Al–Ce–Ni alloys (wt%): near-eutectic Al–10Ce–5Ni (with both eutectic and hypoeutectic regions), hypoeutectic Al-7.5Ce-3.75Ni (with numerous primary Al dendrites), and hypereutectic Al-12.5Ce-6.25Ni (with coarse, blocky primary Al3Ni and Al11Ce3 precipitates and some primary Al dendrites). Depending on the alloy composition and local solidification conditions, the following eutectic morphologies are found: (i) coarse Al–Ce eutectic colonies where Al 11 Ce 3 is in the form of “Chinese script”, (ii) intermingled regions of binary Al–Ce and Al–Ni eutectic colonies, with finer Al 3 Ni and Al 11 Ce 3 fibers, (iii) large ternary eutectic colonies, where the binary Al 3 Ni and Al 11 Ce 3 phases are alternating or intertwining within the individual, fine fibers (diameters of ~60–170 nm, depending on solidification rates), and (iv) ternary eutectic zones (between primary Al dendrites), where fine Al 3 Ni and Al 11 Ce 3 build up a 3D-interconnected network. The high volume fraction of intermetallic phases and extremely fine eutectic spacing/fiber diameter both contribute to high ambient strengthening (higher as-cast microhardness than binary Al–Ce or Al–Ni), and also provide enhanced creep resistance at 300 and 350 °C. Additionally, the alloys are coarsening-resistant up to 425 °C for extended periods, with a gradual decrease in microhardness. The alloys aged at 400 °C to 1050 h show fiber fragmentation and coarsening of the resulting particles, with the faster-diffusing Ni driving more rapid coarsening of the Al 3 Ni particles which engulf finer, more stable Al 11 Ce 3 particles. Severe overaging (performed at 590 °C for 24 h) leads to Al 3 Ni and Al 11 Ce 3 spheroids which remain submicron-sized in eutectic colonies, but micron-sized at colony boundary and at Al dendrite-eutectic interface. Creep resistance at 300 °C of overaged Al–10Ce–5Ni remains substantial, consistent with load-transfer based composite strengthening being an important strengthening mechanism in these alloys, making them excellent candidates for replacement of heavier steel or titanium parts operating under stress up to 300 °C.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Examination of Semi-Analytical Solution Methods in the Coarse Operator of Parareal Algorithm for Power System Simulation

With continuing advances in high-performance parallel computing platforms, parallel algorithms have become powerful tools for development of faster than real-time power system dynamic simulations. In particular, it has been demonstrated in recent years that parallel-in-time (Parareal) algorithms have the potential to achieve such an ambitious goal. Here, the selection of a fast and reasonably accurate coarse operator of the Parareal algorithm is crucial for its effective utilization and performance. This paper examines semi-analytical solution (SAS) methods as the coarse operators of the Parareal algorithm and explores performance of the SAS methods to the standard numerical time integration methods. Two promising time-power series-based SAS methods were considered; Adomian decomposition method and Homotopy analysis method with a windowing approach for improving the convergence. Numerical performance case studies on 10-generator 39-bus system and 327-generator 2383-bus system were performed for these coarse operators over different disturbances, evaluating the number of Parareal iterations, computational time, and stability of convergence. All the coarse operators tested with different scenarios have converged to the same corresponding true solution (if they are convergent) and the SAS methods provide comparable computational speed, while having more stable convergence to the true solution in many cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

On the Training and Generalization of Deep Operator Networks

Here, we present a novel training method for deep operator networks (DeepONets), one of the most popular neural network models for operators. DeepONets are constructed by two subnetworks, namely the branch and trunk networks. Typically, the two subnetworks are trained simultaneously, which amounts to solving a complex optimization problem in a high dimensional space. In addition, the nonconvex and nonlinear nature makes training very challenging. To tackle such a challenge, we propose a two-step training method that trains the trunk network first and then sequentially trains the branch network. The core mechanism is motivated by the divide-and-conquer paradigm and is the decomposition of the entire complex training task into two subtasks with reduced complexity. Therein the Gram–Schmidt orthonormalization process is introduced which significantly improves stability and generalization ability. On the theoretical side, we establish a generalization error estimate in terms of the number of training data, the width of DeepONets, and the number of input and output sensors. Numerical examples are presented to demonstrate the effectiveness of the two-step training method, including Darcy flow in heterogeneous porous media.

deep operator networks↗
DOE OSTI2022

The DPG Method for the Convection-Reaction Problem, Revisited

In this work, we study both conforming and non-conforming versions of the practical DPG method for the convection-reaction problem. We determine that the most common approach for DPG stability analysis – construction of a local Fortin operator – is infeasible for the convection-reaction problem. We then develop a line of argument based on a direct proof of discrete stability; we find that employing a polynomial enrichment for the test space does not suffice for this purpose, motivating the introduction of a (two-element) subgrid mesh. The argument combines mathematical analysis with numerical experiments.

97 MATHEMATICS AND COMPUTING↗