Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Convergence analysis”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Severe storm initiation and development from satellite infrared imagery and Rawinsonde data

The geographical distribution of potential temperatures, mixing ratio, and streamlines of flow patterns at 850, 700, and 500 mb heights are used to understand the prestorm convection and the horizontal convergence of moisture. From the analysis of 21 tornadoes the following conclusions are reached: (1) Strong horizontal convergence of moisture appeared at the 850, 700, and 500 mb levels in the area 12 hours before the storm formation; (2) An abundantly moist atmosphere below 3 km (700 mb) becomes convectively unstable during the time period between 12 and 24 hours before the initiation of the severe storms; (3) Strong winds veering with height with direction parallel to the movement of a dryline, surface fronts, etc; (4) During a 36-hour period, a tropopause height in the areas of interest is lowest at the time of tornadic cloud formation; (5) A train of gravity waves is detected before and during the cloud formation period. Rapid-scan infrared imagery provides near real-time information on the life cycle of the storm which can be summarized as follows: (1) Enhanced convection produced an overshooting cloud top penetrating above the tropopause, making the mass density of the overshooting cloud much greater than the mass density of the surrounding air; (2) The overshooting cloud top collapsed at the end of the mature stage of the cloud development; (3) The tornado touchdown followed the collapse of the overshooting cloud top.

Hung, R. J.↗
DOE OSTI2024

Benchmarking Correlation-Consistent Basis Sets for Frequency-Dependent Polarizabilities with Multiresolution Analysis

This paper presents the first converged frequency-dependent HF polarizability results for general molecules, on a set of 89 closed-shell atoms and molecules. The solver employs multiresolution analysis (MRA) in a multiwavelet basis to compute both ground and response states to a guaranteed precision, which are validated against independent numerical grid calculations on atoms and linear molecules. The MRA ground-state energies and response properties are used to evaluate results in correlation-consistent basis sets up to 5Z augmented with either single or double diffuse functions and core-polarization functions. Systematic trends are revealed through consideration of chemical composition as well as the use of machine learning to cluster convergence trends, the latter suggesting the possibility of learning and correcting basis-set error.

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

Thermodynamic and dynamic contributions to future changes in summer precipitation over Northeast Asia and Korea: a multi-RCM study

This study examines future changes in precipitation over Northeast Asia and Korea using five regional climate model (RCM) simulations driven by single global climate model (GCM) under two representative concentration pathway (RCP) emission scenarios. Focusing on summer season (June–July–August) when heavy rains dominate in this region, future changes in precipitation and associated variables including temperature, moisture, and winds are analyzed by comparing future conditions (2071–2100) with a present climate (1981–2005). Here, physical mechanisms are examined by analyzing moisture flux convergence at 850 hPa level, which is found to have a close relationship to precipitation and by assessing contribution of thermodynamic effect (TH, moisture increase due to warming) and dynamic effect (DY, atmospheric circulation change) to changes in the moisture flux convergence. Overall background warming and moistening are projected over the Northeast Asia with a good inter-RCM agreement, indicating dominant influence of the driving GCM. Also, RCMs consistently project increases in the frequency of heavy rains and the intensification of extreme precipitation over South Korea. Analysis of moisture flux convergence reveals competing impacts between TH and DY. The TH effect contributes to the overall increases in mean precipitation over Northeast Asia and in extreme precipitation over South Korea, irrespective of models and scenarios. However, DY effect is found to induce local-scale precipitation decreases over the central part of the Korean Peninsula with large inter-RCM and inter-scenario differences. Composite analysis of daily anomaly synoptic patterns indicates that extreme precipitation events are mainly associated with the southwest to northeast evolution of large-scale low-pressure system in both present and future climates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Pitfalls in the n -mode representation of vibrational potentials

Simulations of anharmonic vibrational motion rely on computationally expedient representations of the governing potential energy surface. The n-mode representation (n-MR)—effectively a many-body expansion in the space of molecular vibrations—is a general and efficient approach that is often used for this purpose in vibrational self-consistent field (VSCF) calculations and correlated analogues thereof. In the present analysis, a lack of convergence in many VSCF calculations is shown to originate from negative and unbound potentials at truncated orders of the n-MR expansion. For cases of strong anharmonic coupling between modes, the n-MR can both dip below the true global minimum of the potential surface and lead to effective single-mode potentials in VSCF that do not correspond to bound vibrational problems, even for bound total potentials. The present analysis serves mainly as a pathology report of this issue. Furthermore, this insight into the origin of VSCF non-convergence provides a simple, albeit ad hoc, route to correct the problem by “painting in” the full representation of groups of modes that exhibit these negative potentials at little additional computational cost. Somewhat surprisingly, this approach also reasonably approximates the results of the next-higher n-MR order and identifies groups of modes with particularly strong coupling. The method is shown to identify and correct problematic triples of modes—and restore SCF convergence—in two-mode representations of challenging test systems, including the water dimer and trimer, as well as protonated tropine.

Chemistry↗
DOE OSTI2023

Unified analysis of finite-size error for periodic Hartree-Fock and second order Møller-Plesset perturbation theory

Despite decades of practice, finite-size errors in many widely used electronic structure theories for periodic systems remain poorly understood. For periodic systems using a general Monkhorst-Pack grid, there has been no comprehensive and rigorous analysis of the finite-size error in the Hartree-Fock theory (HF) and the second order Møller-Plesset perturbation theory (MP2), which are the simplest wavefunction based method, and the simplest post-Hartree-Fock method, respectively. Such calculations can be viewed as a multi-dimensional integral discretized with certain trapezoidal rules. Due to the Coulomb singularity, the integrand has many points of discontinuity in general, and standard error analysis based on the Euler-Maclaurin formula gives overly pessimistic results. The lack of analytic understanding of finite-size errors also impedes the development of effective finite-size correction schemes. We propose a unified analysis to obtain sharp convergence rates of finite-size errors for the periodic HF and MP2 theories. Our main technical advancement is a generalization of the result of Lyness [Math. Comp. 30 (1976), pp. 1–23] for obtaining sharp convergence rates of the trapezoidal rule for a class of non-smooth integrands. Our result is applicable to three-dimensional bulk systems as well as low dimensional systems (such as nanowires and 2D materials). Our unified analysis also allows us to prove the effectiveness of the Madelung-constant correction to the Fock exchange energy, and the effectiveness of a recently proposed staggered mesh method for periodic MP2 calculations (see X. Xing, X. Li, and L. Lin [J. Chem. Theory Comput. 17 (2021), pp. 4733–4745]). In conclusion, our analysis connects the effectiveness of the staggered mesh method with integrands with removable singularities, and suggests a new staggered mesh method for reducing finite-size errors of periodic HF calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Augmenting the one-shot framework by additional constraints

The (multistep) one-shot method for design optimization problems has been successfully implemented for various applications. To this end, a slowly convergent primal fixed-point iteration of the state equation is augmented by an adjoint iteration and a corresponding preconditioned design update. In this paper we present a modification of the method that allows for additional equality constraints besides the usual state equation. Finally, a retardation analysis and the local convergence of the method in terms of necessary and sufficient conditions are given, which depend on key characteristics of the underlying problem and the quality of the utilized preconditioner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A multigrid waveform relaxation method for solving the poroelasticity equations

In this work, a multigrid waveform relaxation method is proposed for solving a collocated finite difference discretization of the linear Biot’s model. This gives rise to the first space–time multigrid solver for poroelasticity equations in the literature. The waveform relaxation iteration is based on a point-wise Vanka smoother that couples the pressure variable at a grid-point with the displacements around it. A semi-algebraic mode analysis is proposed to theoretically analyze the convergence of the multigrid waveform relaxation algorithm. This analysis is novel since it combines the semi-algebraic analysis, suitable for parabolic problems, with the non-standard analysis for overlapping smoothers. The practical utility of the method is illustrated through several numerical experiments in one and two dimensions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2016

Wavenumber-frequency deconvolution of aeroacoustic microphone phased array data of arbitrary coherence

Deconvolution of aeroacoustic data acquired with microphone phased arrays is a computationally-challenging task for distributed sources with arbitrary coherence. A new technique for performing such deconvolution is proposed. This technique relies on analysis of the array data in the wavenumber-frequency domain, allowing for fast convolution and reduced storage requirements when compared to traditional coherent deconvolution. A positive semidefinite constraint for the iterative deconvolution procedure is implemented and shows improved behavior in terms of quantifiable convergence metrics when compared to a standalone covariance inequality constraint. A series of simulations validates the method's ability to resolve coherence and phase angle relationships between partially-coherent sources, as well as determines convergence criteria for deconvolution analysis. Simulations for point sources near the microphone phased array show potential for handling such data in the wavenumber-frequency domain. In particular, a physics-based integration boundary calculation is described, and can successfully isolate sources and track the appropriate integration bounds with and without the presence of flow. Magnitude and phase relationships between multiple sources are successfully extracted. Limitations of the deconvolution technique are determined from the simulations, particularly in the context of a simulated acoustic field in a closed test section wind tunnel with strong boundary layer contamination. A final application to a trailing edge noise experiment conducted in an open-jet wind tunnel matches best estimates of acoustic levels from traditional calculation methods and qualitatively assesses the coherence characteristics of the trailing edge noise source.

phased array↗
NASA NTRS1993

Tectonics and volcanism of Eastern Aphrodite Terra, Venus - No subduction, no spreading

Eastern Aphrodite Terra, a deformed region with high topographic relief on Venus, has been interpreted as analogous to a terrestrial extensional or convergent plate boundary. However, analysis of geological and structural relations indicates that the tectonics of eastern Aphrodite Terra is dominated by blistering of the crust by magma diapirs. The findings imply that, within this region, vertical tectonism dominates over horizontal tectonism and, consequently, that this region is neither a divergent nor a convergent plate boundary.

Hansen, Vicki L.↗
NASA NTRS1997

Gradient-Based Aerodynamic Shape Optimization Using ADI Method for Large-Scale Problems

A gradient-based shape optimization methodology, that is intended for practical three-dimensional aerodynamic applications, has been developed. It is based on the quasi-analytical sensitivities. The flow analysis is rendered by a fully implicit, finite volume formulation of the Euler equations.The aerodynamic sensitivity equation is solved using the alternating-direction-implicit (ADI) algorithm for memory efficiency. A flexible wing geometry model, that is based on surface parameterization and platform schedules, is utilized. The present methodology and its components have been tested via several comparisons. Initially, the flow analysis for for a wing is compared with those obtained using an unfactored, preconditioned conjugate gradient approach (PCG), and an extensively validated CFD code. Then, the sensitivities computed with the present method have been compared with those obtained using the finite-difference and the PCG approaches. Effects of grid refinement and convergence tolerance on the analysis and shape optimization have been explored. Finally the new procedure has been demonstrated in the design of a cranked arrow wing at Mach 2.4. Despite the expected increase in the computational time, the results indicate that shape optimization, which require large numbers of grid points can be resolved with a gradient-based approach.

Pandya, Mohagna J.↗
DOE OSTI2017

The shifted boundary method for embedded domain computations. Part I: Poisson and Stokes problems

We propose a new finite element method for embedded domain computations, which falls in the category of surrogate/approximate boundary algorithms. The key feature of the proposed approach is the idea of shifting the location where boundary conditions are applied from the true to the surrogate boundary, and to appropriately modify the shifted boundary conditions, enforced weakly, in order to preserve optimal convergence rates of the numerical solution. This process yields a method which, in our view, is simple, efficient, and also robust, since it is not affected by the small-cut-cell problem. Although general in nature, here we apply this new concept to the Poisson and Stokes problems. We present in particular the full analysis of stability and convergence for the case of the Poisson operator, and numerical tests for both the Poisson and Stokes equations, for geometries of progressively higher complexity.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1987

Thermal analysis and design of a cooling system for a Mach 14 nozzle

The analysis and design of a Mach 14 converging diverging nozzle wall liner is provided. The analysis indicates that: no fin on the coolant side of the nozzle wall is optimum, the thermal stresses are dominant, and the critical area is very near the throat. The molybdenum alloy TZM, with a wall thickness of 2.0 mm in the throat area, appears to be the only material capable of meeting design requirements. Additionally, cooling water at 2000 psia with a flow velocity of 25 m/s in the coolant passages is required.

Mullisen, Ronald↗
DOE OSTI2022

Anderson acceleration stability in NDA-accelerated k-eigenvalue problems

Anderson acceleration (AA) has been used to improve the stability and convergence rate of multiphysics iterative methods for reactor analysis. Most applications studied assume a tightly converged solution for the different physics problems, and AA is usually applied to state variables like temperature, density, and heat generation rate. In this paper, we study the theoretical performance of AA in NDA-accelerated k-eigenvalue problems. The problems and algorithms studied are simplified from the coupled iteration scheme adopted by MPACT and many other high-fidelity whole-core reactor codes. Compared to previous analyses of AA for these iteration schemes, we study the case with a partially converged neutronics solution and possibly partially converged nonlinear diffusion acceleration (NDA)/coarse mesh finite difference (CMFD) solutions. We observe that the performance of the iteration scheme with AA is very sensitive to the initial guess and is affected by the partially converged CMFD solutions. When the NDA solution is fully converged, using AA cannot achieve the optimal convergence rate in large-sized problems. Conversely, if the NDA solution is partially converged, the iteration scheme with AA can diverge or converge extremely slowly. It is found that the loss of robustness for AA is due to the fact that it is applied to the iterative subspace of state variables rather than the fundamental unknowns of the governing equations. To improve the robustness, the scalar flux should also be considered in the implementation of AA. After considering the residuals of flux, we observe that the stability is regardless of the partial convergence of NDA solutions. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Cosmological constraints from the convergence 1-point probability distribution

Here, we examine the cosmological information available from the 1-point probability density function (PDF) of the weak-lensing convergence field, utilizing fast l-picola simulations and a Fisher analysis. We find competitive constraints in the Ω m–σ8 plane from the convergence PDF with 188 arcmin 2 pixels compared to the cosmic shear power spectrum with an equivalent number of modes (ℓ < 886). The convergence PDF also partially breaks the degeneracy cosmic shear exhibits in that parameter space. A joint analysis of the convergence PDF and shear 2-point function also reduces the impact of shape measurement systematics, to which the PDF is less susceptible, and improves the total figure of merit by a factor of 2–3, depending on the level of systematics. Finally, we present a correction factor necessary for calculating the unbiased Fisher information from finite differences using a limited number of cosmological simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

NUMERICAL CONVERGENCE IN SMOOTHED PARTICLE HYDRODYNAMICS

We study the convergence properties of smoothed particle hydrodynamics (SPH) using numerical tests and simple analytic considerations. Our analysis shows that formal numerical convergence is possible in SPH only in the joint limit N → ∞, h → 0, and N{sub nb} → ∞, where N is the total number of particles, h is the smoothing length, and N{sub nb} is the number of neighbor particles within the smoothing volume used to compute smoothed estimates. Previous work has generally assumed that the conditions N → ∞ and h → 0 are sufficient to achieve convergence, while holding N{sub nb} fixed. We demonstrate that if N{sub nb} is held fixed as the resolution is increased, there will be a residual source of error that does not vanish as N → ∞ and h → 0. Formal numerical convergence in SPH is possible only if N{sub nb} is increased systematically as the resolution is improved. Using analytic arguments, we derive an optimal compromise scaling for N{sub nb} by requiring that this source of error balance that present in the smoothing procedure. For typical choices of the smoothing kernel, we find N{sub nb} ∝N {sup 0.5}. This means that if SPH is to be used as a numerically convergent method, the required computational cost does not scale with particle number as O(N), but rather as O(N {sup 1} {sup +} {sup δ}), where δ ≈ 0.5, with a weak dependence on the form of the smoothing kernel.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2018

Quantitative validation of Ames IR intensity and new line lists for (32/33/34)S(16)O2, (32)S(18)O2 and (16)O(32)S(18)O

The quality of Ames-296K SO2 Infrared (IR) line list intensities is first validated by quantitative exploration of several dipole moment surfaces (DMSs) and partition sum convergence. The DMSs are computed with several of Dunning’s correlation-consistent basis sets and their vibrational dependence are compared to the empirical model derived from Stark effect experiments reported by D. Patel, D. Margolese, and T.R. Dykea [J.Chem.Phys. 70, 2740 (1979)]. The effective dipole deviations from the DMS adopted in the Ames IR lists is 0.2-0.4% for vibrational states up to 3v3. The vibrational dependence of the dipole moment is also in good agreement, except for nv1. Partition sum convergence at 296K is confirmed by new calculations with rotational quantum number J up to 150 and upper state E’ up to 8000 cm-1. The isotopologue consistency of the Ames IR line lists is superior relative to the regular Effective Hamiltonian (EH) models and Effective Dipole Moment (EDM) models. The v1+v2 and v2+v3 intensity consistency check reveals the recently reported experimental intensities need significant improvement or re-analysis. After the accuracy, convergence, and isotopologue consistency have been confirmed, the theoretical Ames-296K intensities are combined with the experimental line positions or EH models that experimental spectroscopists published after 2009. Three high-resolution IR line sets are reported for the (32/33/34)S(16)O2, (32)S(18)O2 and (16)O(32)S(18)O isotopologues: (1) the “New Lines Sets” include experimentally measured line positions; (2) the “Expanded Line Sets” include possible transitions among new rovibrational levels assigned in experiments and ground state (GS) levels predicted by reliable EH models; (3) the “Ames + MARVEL Sets” include possible transitions among all those levels reported in a recent MARVEL analysis. [Tóbiás et al, JQSRT 208, 152 (2018)]. Compared to the limited data in High-resolution TRANsmission molecular absorption database (HITRAN), these line sets have significantly improved the data coverage up to 4000 cm-1. Some missing bands can be traced to the unpublished experimental data. The isotopologue consistency of these line sets will help identify the uncertainties and defects in the experimental EH and EDM models. These line sets are good candidates for the next HITRAN update, if line shape parameters are available. The line sets can be downloaded from supplementary files or from the Ames Molecular Spectroscopic Database at http://huang.seti.org.

Xinchuan Huang↗