Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence behaviors”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Nonlinear dynamics and numerical uncertainties in CFD

The application of nonlinear dynamics to improve the understanding of numerical uncertainties in computational fluid dynamics (CFD) is reviewed. Elementary examples in the use of dynamics to explain the nonlinear phenomena and spurious behavior that occur in numerics are given. The role of dynamics in the understanding of long time behavior of numerical integrations and the nonlinear stability, convergence, and reliability of using time-marching, approaches for obtaining steady-state numerical solutions in CFD is explained. The study is complemented with spurious behavior observed in CFD computations.

Yee, H. C.↗
DOE OSTI2023

Multigrid Reduction in Time for Chaotic Dynamical Systems

As CPU clock speeds have stagnated and high performance computers continue to have ever higher core counts, increased parallelism is needed to take advantage of these new architectures. Traditional serial time-marching schemes can be a significant bottleneck, as many types of simulations require large numbers of time-steps which must be computed sequentially. Parallel-in-time schemes, such as the Multigrid Reduction in Time (MGRIT) method, remedy this by parallelizing across time-steps and have shown promising results for parabolic problems. However, chaotic problems have proved more difficult, since chaotic initial value problems (IVPs) are inherently ill-conditioned. MGRIT relies on a hierarchy of successively coarser time-grids to iteratively correct the solution on the finest time-grid, but due to the nature of chaotic systems, small inaccuracies on the coarser levels can be greatly magnified and lead to poor coarse-grid corrections. Here we introduce a modified MGRIT algorithm based on an existing quadratically converging nonlinear extension to the multigrid Full Approximation Scheme (FAS), as well as a novel time-coarsening scheme. Together, these approaches better capture long-term chaotic behavior on coarse-grids and greatly improve convergence of MGRIT for chaotic IVPs. Further, we introduce a novel low-memory variant of the algorithm for solving chaotic PDEs with MGRIT which not only solves the IVP, but also provides estimates for the unstable Lyapunov vectors of the system. Finally, we provide supporting numerical results for the Lorenz system and demonstrate parallel speedup for the chaotic Kuramoto–Sivashinsky PDE over a significantly longer time-domain than in previous works.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1996

Some Aspects of Nonlinear Dynamics and CFD

The application of nonlinear dynamics to improve the understanding of numerical uncertainties in computational fluid dynamics (CFD) is reviewed. Elementary examples in the use of dynamics to explain the nonlinear phenomena and spurious behavior that occur in numerics are given. The role of dynamics in the understanding of long time behavior of numerical integrations and the nonlinear stability, convergence, and reliability of using time-marching approaches for obtaining steady-state numerical solutions in CFD is explained. The study is complemented with examples of spurious behavior observed in CFD computations.

Yee, Helen C.↗
NASA NTRS1997

Dynamics of Numerics & Spurious Behaviors in CFD Computations

The global nonlinear behavior of finite discretizations for constant time steps and fixed or adaptive grid spacings is studied using tools from dynamical systems theory. Detailed analysis of commonly used temporal and spatial discretizations for simple model problems is presented. The role of dynamics in the understanding of long time behavior of numerical integration and the nonlinear stability, convergence, and reliability of using time-marching approaches for obtaining steady-state numerical solutions in computational fluid dynamics (CFD) is explored. The study is complemented with examples of spurious behavior observed in steady and unsteady CFD computations. The CFD examples were chosen to illustrate non-apparent spurious behavior that was difficult to detect without extensive grid and temporal refinement studies and some knowledge from dynamical systems theory. Studies revealed the various possible dangers of misinterpreting numerical simulation of realistic complex flows that are constrained by available computing power. In large scale computations where the physics of the problem under study is not well understood and numerical simulations are the only viable means of solution, extreme care must be taken in both computation and interpretation of the numerical data. The goal of this paper is to explore the important role that dynamical systems theory can play in the understanding of the global nonlinear behavior of numerical algorithms and to aid the identification of the sources of numerical uncertainties in CFD.

Yee, Helen C.↗
NASA NTRS2016

Supplement of: The Influence of Volcanic Eruptions on the Climate of Tropical South America During the Last Millennium in an Isotope-Enabled General Circulation Model

Currently, little is known on how volcanic eruptions impact large-scale climate phenomena such as South American paleo-intertropical Convergence Zone (ITCZ) position and summer monsoon behavior. In this paper, an analysis of observations and model simulations is employed to assess the influence of large volcanic eruptions on the climate of tropical South America. This problem is first considered for historically recent volcanic episodes for which more observations are available but where fewer events exist and the confounding effects of El NioSouthern Oscillation (ENSO) lead to inconclusive interpretation of the impact of volcanic eruptions at the continental scale. Therefore, we also examine a greater number of reconstructed volcanic events for the period 850CE to present that are incorporated into the NASA GISS ModelE2-R simulation of the last millennium.An advantage of this model is its ability to explicitly track water isotopologues throughout the hydrologic cycle and simulating the isotopic imprint following a large eruption. This effectively removes a degree of uncertainty associated with error-prone conversion of isotopic signals into climate variables, and allows for a direct comparison between GISS simulations and paleoclimate proxy records.Our analysis reveals that both precipitation and oxygen isotope variability respond with a distinct seasonal and spatial structure across tropical South America following an eruption. During austral winter, the heavy oxygen isotope in precipitation is enriched, likely due to reduced moisture convergence in the ITCZ domain and reduced rainfall over northern South America. During austral summer, however, more negative values of the precipitation isotopic composition are simulated over Amazonia, despite reductions in rainfall, suggesting that the isotopic response is not a simple function of the amount effect. During the South American monsoon season, the amplitude of the temperature response to volcanic forcing is larger than the rather weak and spatially less coherent precipitation signal, complicating the isotopic response to changes in the hydrologic cycle.

Colose, Christopher↗
NASA NTRS2016

The Influence of Volcanic Eruptions on the Climate of Tropical South America During the Last Millennium in an Isotope-Enabled General Circulation Model

Currently, little is known on how volcanic eruptions impact large-scale climate phenomena such as South American paleo-intertropical Convergence Zone (ITCZ) position and summer monsoon behavior. In this paper, an analysis of observations and model simulations is employed to assess the influence of large volcanic eruptions on the climate of tropical South America. This problem is first considered for historically recent volcanic episodes for which more observations are available but where fewer events exist and the confounding effects of El Niño-Southern Oscillation (ENSO) lead to inconclusive interpretation of the impact of volcanic eruptions at the continental scale. Therefore, we also examine a greater number of reconstructed volcanic events for the period 850 CE to present that are incorporated into the NASA GISS ModelE2-R simulation of the last millennium. An advantage of this model is its ability to explicitly track water isotopologues throughout the hydrologic cycle and simulating the isotopic imprint following a large eruption. This effectively removes a degree of uncertainty associated with error-prone conversion of isotopic signals into climate variables, and allows for a direct comparison between GISS simulations and paleoclimate proxy records. Our analysis reveals that both precipitation and oxygen isotope variability respond with a distinct seasonal and spatial structure across tropical South America following an eruption. During austral winter, the heavy oxygen isotope in precipitation is enriched, likely due to reduced moisture convergence in the ITCZ domain and reduced rainfall over northern South America. During austral summer, however, more negative values of the precipitation isotopic composition are simulated over Amazonia, despite reductions in rainfall, suggesting that the isotopic response is not a simple function of the "amount effect". During the South American monsoon season, the amplitude of the temperature response to volcanic forcing is larger than the rather weak and spatially less coherent precipitation signal, complicating the isotopic response to changes in the hydrologic cycle.

Colose, Christopher M.↗
DOE OSTI2021

Evolution of Chain Dynamics and Oxidation States with Increasing Chain Length for a Donor–Acceptor-Conjugated Oligomer Series

While it is known that the chain length strongly affects the properties of π-conjugated polymers, the effects of chain length on the molecular structure, chain conformation, and oxidation state properties in donor–acceptor-type conjugated structures remain unclear. This limits our understanding of how the polymer molecular weight impacts material properties. Here, a discrete and monodisperse oligomer series (n = 3, 5, 7, 9, 15, and 21) and polymers (nPB), composed of the donor 3,4-propylenedioxythiophene (ProDOT), the acceptor benzothiadiazole (BTD), and methylthio end-capping groups, are synthesized by C–H-activated cross-coupling. The molecular structure, molecular weight, and dispersity of each oligomer/polymer are thoroughly characterized by nuclear magnetic resonance spectroscopy, mass spectrometry, and gel permeation chromatography. This series reveals a rod-to-coil transition at n = 15 and coil formation at polymer length scales of ~28 units via solution small-angle neutron scattering characterization. The oxidation states are deciphered via cyclic voltammetry, differential pulse voltammetry, spectroelectrochemistry, and density functional theory calculations. Oligomers 3–9 undergo successive one-electron oxidation steps, while 15 and higher undergo multielectron oxidations per step in CH 2 Cl 2 –TBAPF6 at a Pt ultramicroelectrode. The electronic transition of each oxidation state (1+, 2+, 3+, etc.) is tracked by absorption spectroscopy, revealing a “bipolaron to di-polaron” transition at n = 7 at which shorter oligomers prefer bipolaron formation and longer oligomers prefer di-polaron formation in their dication states. Furthermore, oxidized 15 has a lower electronic transition energy compared to its polymer homologue, as shown by spectroelectrochemistry, revealing a synergy between the chain length and the oxidation state properties. This study proves that the convergence limit between small molecule and polymer behavior occurs at approximately 15 units and highlights the property transitions that occur as a function of chain length for a donor–acceptor class of conjugated organic materials.

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

A physics-constrained neural network for multiphase flows

The present study develops a physics-constrained neural network (PCNN) to predict sequential patterns and motions of multiphase flows (MPFs), which includes strong interactions among various fluid phases. To predict the order parameters, which locate individual phases in the future time, a neural network (NN) is applied to quickly infer the dynamics of the phases by encoding observations. The multiphase consistent and conservative boundedness mapping algorithm (MCBOM) is next implemented to correct the predicted order parameters. This enforces the predicted order parameters to strictly satisfy the mass conservation, the summation of the volume fractions of the phases to be unity, the consistency of reduction, and the boundedness of the order parameters. Then, the density of the fluid mixture is updated from the corrected order parameters. Finally, the velocity in the future time is predicted by another NN with the same network structure, but the conservation of momentum is included in the loss function to shrink the parameter space. The proposed PCNN for MPFs sequentially performs (NN)-(MCBOM)-(NN), which avoids nonphysical behaviors of the order parameters, accelerates the convergence, and requires fewer data to make predictions. Numerical experiments demonstrate that the proposed PCNN is capable of predicting MPFs effectively.

Mechanics↗
NASA NTRS1956

Extrapolation techniques applied to matrix methods in neutron diffusion problems

A general matrix method is developed for the solution of characteristic-value problems of the type arising in many physical applications. The scheme employed is essentially that of Gauss and Seidel with appropriate modifications needed to make it applicable to characteristic-value problems. An iterative procedure produces a sequence of estimates to the answer; and extrapolation techniques, based upon previous behavior of iterants, are utilized in speeding convergence. Theoretically sound limits are placed on the magnitude of the extrapolation that may be tolerated. This matrix method is applied to the problem of finding criticality and neutron fluxes in a nuclear reactor with control rods. The two-dimensional finite-difference approximation to the two-group neutron fluxes in a nuclear reactor with control rods. The two-dimensional finite-difference approximation to the two-group neutron-diffusion equations is treated. Results for this example are indicated.

Mccready, Robert R↗
NASA NTRS2016

Closing the Contrast Gap Between Testbed and Model Prediction with WFIRST-CGI Shaped Pupil Coronagraph

JPL has recently passed an important milestone in its technology development for a proposed NASA WFIRST mission: demonstration of a better than 1x10-8 contrast over broad bandwidth (10%) on both shaped pupil coronagraph (SPC) and hybrid Lyot coronagraph (HLC) testbeds with WFIRST telescope optics. Challenges remain, however, in the technology readiness for the proposed mission. One such is the discrepancy between the achieved contrast on testbeds and their corresponding model predictions. A series of testbed diagnosis and modeling activity were planned and (some) carried out on SPC testbed in order to close the gap. A very useful tool we developed was to analyze the measured testbed Jacobian and compare with its model version that was used to control image plane speckle pattern. The difference between these two is the error in the control Jacobian. When the error, which includes both amplitude and phase, was inserted into a separate prediction model with nominal Jacobian, the model prediction became closely matching the SPC testbed behavior in both contrast floor and contrast convergence speed. This offers a new perspective and is a step closer toward model validation for high contrast coronagraph. Further Jacobian analysis and modeling provided clues to the possible sources for the mismatch: the DM misregistration and the underrepresented testbed optics wavefront error (WFE) and the deformable mirror (DM) setting for correcting this WFE. These analyses suggested that high contrast coronagraph has a tight tolerance in the accuracy of its control Jacobian. Modifications to both testbed control model as well as prediction model are being implemented and future works discussed.

diffraction modeling↗
DOE OSTI2021

When band convergence is not beneficial for thermoelectrics

Abstract Band convergence is considered a clear benefit to thermoelectric performance because it increases the charge carrier concentration for a given Fermi level, which typically enhances charge conductivity while preserving the Seebeck coefficient. However, this advantage hinges on the assumption that interband scattering of carriers is weak or insignificant. With first-principles treatment of electron-phonon scattering in the CaMg 2 Sb 2 -CaZn 2 Sb 2 Zintl system and full Heusler Sr 2 SbAu, we demonstrate that the benefit of band convergence can be intrinsically negated by interband scattering depending on the manner in which bands converge. In the Zintl alloy, band convergence does not improve weighted mobility or the density-of-states effective mass. We trace the underlying reason to the fact that the bands converge at a one k-point, which induces strong interband scattering of both the deformation-potential and the polar-optical kinds. The case contrasts with band convergence at distant k-points (as in the full Heusler), which better preserves the single-band scattering behavior thereby successfully leading to improved performance. Therefore, we suggest that band convergence as thermoelectric design principle is best suited to cases in which it occurs at distant k-points.

42 ENGINEERING↗
DOE OSTI2021

Optimizing multigrid reduction-in-time and Parareal coarse-grid operators for linear advection

Parallel-in-time methods, such as multigrid reduction-in-time (MGRIT) and Parareal, provide an attractive option for increasing concurrency when simulating time-dependent partial differential equations (PDEs) in modern high-performance computing environments. While these techniques have been very successful for parabolic equations, it has often been observed that their performance suffers dramatically when applied to advection-dominated problems or purely hyperbolic PDEs using standard rediscretization approaches on coarse grids. In this paper, we apply MGRIT or Parareal to the constant-coefficient linear advection equation, appealing to existing convergence theory to provide insight into the typically nonscalable or even divergent behavior of these solvers for this problem. To overcome these failings, we replace rediscretization on coarse grids with improved coarse-grid operators that are computed by applying optimization techniques to approximately minimize error estimates from the convergence theory. Therefore, one of our main findings is that, in order to obtain fast convergence as for parabolic problems, coarse-grid operators should take into account the behavior of the hyperbolic problem by tracking the characteristic curves. Our approach is tested for schemes of various orders using explicit or implicit Runge–Kutta methods combined with upwind-finite-difference spatial discretizations. In all cases, we obtain scalable convergence in just a handful of iterations, with parallel tests also showing significant speed-ups over sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Superconducting Quantum Metamaterials from Convergence of Soft and Hard Condensed Matter Science

Superconducting quantum metamaterials are expected to exhibit a variety of novel properties, but have been a major challenge to prepare as a result of the lack of appropriate synthetic routes to high-quality materials. Here, the discovery of synthesis routes to block copolymer (BCP) self-assembly-directed niobium nitrides and carbonitrides is described. The resulting materials exhibit unusual structure retention even at temperatures as high as 1000 °C and resulting critical temperature, T c , values comparable to their bulk analogues. Applying the concepts of soft matter self-assembly, it is demonstrated that a series of four different BCP-directed mesostructured superconductors are accessible from a single triblock terpolymer. Resulting materials display a mesostructure-dependent T c without substantial variation of the XRDmeasured lattice parameters. Lastly, field-dependent magnetization measurements of a sample with double-gyroid morphology show abrupt jumps comparable in overall behavior to flux avalanches. Results suggest a fruitful convergence of soft and hard condensed matter science.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Convergence of Public-Private Benefits in Denver: Surveys and Analyses to Inform Urban Mobility-, Energy-, Infrastructure- and Behavior-Related Innovation: Preprint

Cities, public transit agencies, and new private ride hailing services seek to understand emerging traveler dynamics, the shifting demographics of urban travelers,and new energy-efficient mobility opportunities. This includes exploring how new infrastructure investments, public and private mobility services, and smart-phone mobility apps are reshaping behaviors, demands (e.g. mobility-on-demand services), travel experiences and energy-efficient urban travel preferences. Currently, cities and metropolitan regions are providing and experimenting with many new mobility options, technologies, and personalized information services at the intersection of urban mobility, energy, and infrastructure systems (e.g., new commuter rail). To date, technology alone has not been able to crack the nut of 'creating faster trip times, less congestion, safer streets, and cleaner air for its citizens through fewer cars on the road'. This paper focuses on this gap by offering new concepts and potential for integrated approaches. Accommodating more vehicles miles traveled in cities, without increases in person miles traveled (PMT), could be costly, generating: 1) tremendous demands for new infrastructure, land, road space, materials, and energy; 2) higher traffic fatality risks; and 3) worsening air quality. Therefore, this study focuses on reducing single occupancy vehicle use by enhancing integrated mobility, helping transit and ridehailing increase occupancy in ways that also reduce energy use, and improve quality of life for urban travelers and communities. This study focuses on a survey of urban travelers in Denver, as a representative case study for metropolitan regions experiencing rapid growth, ageing populations, increased urban sprawl, traffic-related delays, and inefficient energy use per PMT.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Comprehensive Basis-Set Testing of Extended Symmetry-Adapted Perturbation Theory and Assessment of Mixed-Basis Combinations to Reduce Cost

Hybrid or “extended” symmetry-adapted perturbation theory (XSAPT) replaces traditional SAPT’s treatment of dispersion with better performing alternatives while at the same time extending two-body (dimer) SAPT to a many-body treatment of polarization using a self-consistent charge embedding procedure. The present work presents a systematic study of how XSAPT interaction energies and energy components converge with respect to the choice of Gaussian basis set. Errors can be reduced in a systematic way using correlation-consistent basis sets, with aug-cc-pVTZ results converged within <0.1 kcal/mol. Similar (if slightly less systematic) behavior is obtained using Karlsruhe basis sets at much lower cost, and we introduce new versions with limited augmentation that are even more efficient. Pople-style basis sets, which are more efficient still, often afford good results if a large number of polarization functions are included. The dispersion models used in XSAPT afford much faster basis-set convergence as compared to the perturbative description of dispersion in conventional SAPT, meaning that “compromise” basis sets (such as jun-cc-pVDZ) are no longer required and benchmark-quality results can be obtained using triple-ζ basis sets. The use of diffuse functions proves to be essential, especially for the description of hydrogen bonds. As a result, the “δ(Hartree–Fock)” correction for high-order induction can be performed in double-ζ basis sets without significant loss of accuracy, leading to a mixed-basis approach that offers 4× speedup over the existing (cubic scaling) XSAPT approach.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1977

On the singular behavior of linear acoustic theory in near-sonic duct flows

The propagation of sound in a converging-diverging duct containing a quasi-one-dimensional steady flow with a high subsonic throat Mach number is studied. The behavior of linearized acoustic theory at the throat of the duct is shown to be singular, and the explicit form of the singularity is determined for two special types of area variation. Numerically computed results showing the development of the singularity are presented. The singularity implies that linearized acoustic theory is invalid for sound propagating through a throat carrying a near-sonic flow.

Myers, M. K.↗
DOE OSTI2025

Self-regulating behavior of hybrid membrane systems as demonstrated in an element-scale forward osmosis-reverse osmosis hybrid system

Hybrid membrane systems can be difficult to design due to the requisite flow rate matching between up- and downstream unit operations. In this work, we use a forward osmosis-reverse osmosis (FO-RO) hybrid system to demonstrate how some membrane systems can exhibit self-regulating behavior due to osmotic coupling. This can reduce the need for complex control systems for flow balancing. We show this behavior using a module-scale test bed that can mimic the behavior of larger scale operations. The system shows permeate flow rate near-convergence between the FO and RO modules after startup or when perturbed by a change in RO module pressure. The behavior of this hybrid system demonstrates that some membrane operations can exploit osmotic interdependence, rather than expensive control systems, to achieve steady state operation.

Debottlenecking↗
DOE OSTI2020

Dirty dishes or dirty laundry? Comparing two methods for quantifying American consumers' preferences for load management in a smart home

One challenge of transitioning to renewable energy is that household electricity use and renewable generation are often misaligned. Smart home energy management systems hold promise for shifting usage to match generation, but these systems need to be designed with the occupants’ preferences in mind. The purpose of the present research is to compare two approaches for collecting and modeling consumers’ load management preferences, both of which are amenable to use in a home energy management system. Specifically, we examine the performance of Simple Multi-Attribute Rating Technique Exploiting Ranks (SMARTER) and Analytic Hierarchy Process (AHP) in quantifying consumers’ preferences regarding air temperature (air conditioning and heating), water heating, dishwashing, clothes washing and drying, monetary costs, environmental impacts, and comfort/convenience. Two studies are presented: Study 1 examines the SMARTER approach, and Study 2 focuses on the AHP approach. In both studies, online surveys (N SMARTER = 956 and N AHP = 1023) were conducted to elicit preferences from participants across the United States. The preferences modeled by both approaches were validated based on (a) their ability to predict participants’ choices in a Discrete Choice Experiment and (b) their convergence with previous research on load-shifting behavior. The validation procedure suggests that the SMARTER approach is superior in modeling consumers’ preferences for load management. Overall, this research lays the groundwork for designing a smart home interface capable of collecting occupants’ preferences and using those preferences to deliver improved occupant comfort, lower operating costs, reduced environmental impact, and more significant demand response than exists today.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗