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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
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.↗
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↗
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↗
NASA NTRS1995

An all-at-once reduced Hessian SQP scheme for aerodynamic design optimization

This paper introduces a computational scheme for solving a class of aerodynamic design problems that can be posed as nonlinear equality constrained optimizations. The scheme treats the flow and design variables as independent variables, and solves the constrained optimization problem via reduced Hessian successive quadratic programming. It updates the design and flow variables simultaneously at each iteration and allows flow variables to be infeasible before convergence. The solution of an adjoint flow equation is never needed. In addition, a range space basis is chosen so that in a certain sense the 'cross term' ignored in reduced Hessian SQP methods is minimized. Numerical results for a nozzle design using the quasi-one-dimensional Euler equations show that this scheme is computationally efficient and robust. The computational cost of a typical nozzle design is only a fraction more than that of the corresponding analysis flow calculation. Superlinear convergence is also observed, which agrees with the theoretical properties of this scheme. All optimal solutions are obtained by starting far away from the final solution.

Feng, Dan↗
NASA NTRS2006

High Fidelity System Simulation of Multiple Components in Support of the UEET Program

The High Fidelity System Simulation effort has addressed various important objectives to enable additional capability within the NPSS framework. The scope emphasized High Pressure Turbine and High Pressure Compressor components. Initial effort was directed at developing and validating intermediate fidelity NPSS model using PD geometry and extended to high-fidelity NPSS model by overlaying detailed geometry to validate CFD against rig data. Both "feedforward" and feedback" approaches of analysis zooming was employed to enable system simulation capability in NPSS. These approaches have certain benefits and applicability in terms of specific applications "feedback" zooming allows the flow-up of information from high-fidelity analysis to be used to update the NPSS model results by forcing the NPSS solver to converge to high-fidelity analysis predictions. This apporach is effective in improving the accuracy of the NPSS model; however, it can only be used in circumstances where there is a clear physics-based strategy to flow up the high-fidelity analysis results to update the NPSS system model. "Feed-forward" zooming approach is more broadly useful in terms of enabling detailed analysis at early stages of design for a specified set of critical operating points and using these analysis results to drive design decisions early in the development process.

Plybon, Ronald C.↗
NASA NTRS2001

Modeling and Analysis of Wrinkled Membranes: An Overview

Thin-film membranes are basic elements of a variety of space inflatable/deployable structures. Wrinkling degrades the performance and reliability of these membrane structures, and hence has been a topic of continued interest. Wrinkling analysis of membranes for general geometry and arbitrary boundary conditions is quite challenging. The objective of this presentation is two-fold. Firstly, the existing models of wrinkled membranes and related numerical solution methods are reviewed. The important issues to be discussed are the capability of a membrane model to characterize taut, wrinkled and slack states of membranes in a consistent and physically reasonable manner; the ability of a wrinkling analysis method to predict the formation and growth of wrinkled regions, and to determine out-of-plane deformation and wrinkled waves; the convergence of a numerical solution method for wrinkling analysis; and the compatibility of a wrinkling analysis with general-purpose finite element codes. According to this review, several opening issues in modeling and analysis of wrinkled membranes that are to be addressed in future research are summarized, The second objective of this presentation is to discuss a newly developed membrane model of two viable parameters (2-VP model) and associated parametric finite element method (PFEM) for wrinkling analysis are introduced. The innovations and advantages of the proposed membrane model and PFEM-based wrinkling analysis are: (1) Via a unified stress-strain relation; the 2-VP model treat the taut, wrinkled, and slack states of membranes consistently; (2) The PFEM-based wrinkling analysis has guaranteed convergence; (3) The 2-VP model along with PFEM is capable of predicting membrane out-of-plane deformations; and (4) The PFEM can be integrated into any existing finite element code. Preliminary numerical examples are also included in this presentation to demonstrate the 2-VP model and PFEM-based wrinkling analysis approach.

Yang, B.↗
DOE OSTI2024

Cosmology from cross-correlation of ACT-DR4 CMB lensing and DES-Y3 cosmic shear

Cross-correlation between weak lensing of the Cosmic Microwave Background (CMB) and weak lensing of galaxies offers a way to place robust constraints on cosmological and astrophysical parameters with reduced sensitivity to certain systematic effects affecting individual surveys. We measure the angular cross-power spectrum between the Atacama Cosmology Telescope (ACT) DR4 CMB lensing and the galaxy weak lensing measured by the Dark Energy Survey (DES) Y3 data. Our baseline analysis uses the CMB convergence map derived from ACT-DR4 and Planck data, where most of the contamination due to the thermal Sunyaev Zel’dovich effect is removed, thus avoiding important systematics in the cross-correlation. In our modelling, we consider the nuisance parameters of the photometric uncertainty, multiplicative shear bias and intrinsic alignment of galaxies. The resulting cross-power spectrum has a signal-to-noise ratio = 7.1 and passes a set of null tests. We use it to infer the amplitude of the fluctuations in the matter distribution (S_8 ≡ σ_8(Ω_m/0.3)^0.5 = 0.782 ± 0.059) with informative but well-motivated priors on the nuisance parameters. We also investigate the validity of these priors by significantly relaxing them and checking the consistency of the resulting posteriors, finding them consistent, albeit only with relatively weak constraints. This cross-correlation measurement will improve significantly with the new ACT-DR6 lensing map and form a key component of the joint 6×2pt analysis between DES and ACT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Structures of drug-specific monoclonal antibodies bound to opioids and nicotine reveal a common mode of binding

Opioid-related fatal overdoses have reached epidemic proportions. Because existing treatments for opioid use disorders offer limited long-term protection, accelerating the development of newer approaches is critical. Monoclonal antibodies (mAbs) are an emerging treatment strategy that targets and sequesters selected opioids in the bloodstream, reducing drug distribution across the blood-brain barrier, thus preventing or reversing opioid toxicity. We previously identified a series of murine mAbs with high affinity and selectivity for oxycodone, morphine, fentanyl, and nicotine. To determine their binding mechanism, we used X-ray crystallography to solve the structures of mAbs bound to their respective targets, to 2.2 Å resolution or higher. Structural analysis showed a critical convergent hydrogen bonding mode that is dependent on a glutamic acid residue in the mAbs’ heavy chain and a tertiary amine of the ligand. Further, characterizing drug-mAb complexes represents a significant step toward rational antibody engineering and future manufacturing activities to support clinical evaluation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Error Bounds for Dynamical Spectral Estimation

Dynamical spectral estimation is a well-established numerical approach for estimating eigenvalues and eigenfunctions of the Markov transition operator from trajectory data. Although the approach has been widely applied in biomolecular simulations, its error properties remain poorly understood. Here we analyze the error of a dynamical spectral estimation method called “the variational approach to conformational dynamics" (VAC). We bound the approximation error and estimation error for VAC estimates. Our analysis establishes VAC's convergence properties and suggests new strategies for tuning VAC to improve accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Extratropical Cloud Feedback Constrained by Cloud Sources and Sinks in Cyclones

Constraining cloud feedback in global climate models (GCMs) using observations is important for establishing accurate predictions of future climate. Uncertainty in shortwave cloud feedback (SW FB ) dominates uncertainty in total cloud feedback. Recent studies show a shift toward more positive extratropical SW FB in the latest generations of GCMs leading to the emergence of very high equilibrium climate sensitivity (ECS). In this study, we use precipitation efficiency and albedo susceptibility to constrain liquid water path (LWP) response to warming and SW FB in the Southern Ocean (SO; 50°–80°S). We analyze precipitation in extratropical cyclones (ECs) to learn about extratropical condensed water sink processes, combined with observations of clouds and moisture convergence, and use the analysis to better understand and constrain SW FB . We utilize a perturbed parameter ensemble (PPE) hosted in the Community Atmosphere Model, version 6 (CAM6), to provide a constraint on SW FB based on observations from Clouds and the Earth’s Radiant Energy System (CERES) and Multisensor Advanced Climatology of LWP (MAC-LWP). We apply Gaussian process regression to emulate the model response to all parameters perturbed in the PPE. Confronting the emulator output with observations provides a new estimated response of Earth to global warming. Furthermore, our new estimates of SO LWP reduce the PPE range by 66%–72%, which results in a shortwave cloud radiative effect estimated range that is 27%–34% less than the PPE range. Observations suggest a more positive SO SW FB than the Community Earth System Model, version 2 (CESM2), and consequently do not reject the high climate sensitivity GCMs emerging from the Coupled Model Intercomparison Project phase 6 (CMIP6).

Atmosphere↗
DOE OSTI2022

Simulating Atmospheric Boundary Layer Turbulence with Nek5000/RS

We present large-eddy-simulation (LES) modeling approaches for the simulation of atmospheric boundary layer turbulence that are of direct relevance to wind energy production. In this report, we study a GABLS benchmark problem using high-order spectral element code Nek5000/RS, which is supported under the DOE’s Exascale Computing Project (ECP) Center for Efficient Exascale Discretizations (CEED) project, targeting application simulations on various acceleration-device based exascale computing platforms [1, 2]. We demonstrate our newly developed subgrid-scale (SGS) models based on high-pass filter (HPF), mean-field eddy viscosity (MFEV), and Smagorinsky (SMG) with no-slip and traction boundary conditions, provided with low-order statistics, convergence and turbulent structure analysis. The model fidelity and scaling performance of Nek5000/RS on DOE’s leadership computing platforms in comparison to those of AMR-Wind, a block-structured second-order finite-volume code with adaptive-mesh-refinement capabilities, are discussed in [3].

17 WIND ENERGY↗
NASA NTRS1973

Optimization of a corrugated stiffened composite panel under uniaxial compression

An approach of structural optimization has been used to optimize the weight of a simply supported, corrugated hat stiffened composite panel under uniaxial compression. The approach consists of the employment of nonlinear mathematical programming techniques to reach an optimum solution. Some simplifying assumptions are made in the stress analysis to obtain faster convergence to an optimum solution. With these simplifying assumptions the number of unknown design parameters is reduced to twelve.

Agarwal, B. L.↗
NASA NTRS1981

The stability of pseudospectral-Chebyshev methods

The pseudospectral-Chebyshev methods are shown to be convergent in variable coefficient problems and, in some cases, hyperbolic problems. The analysis demonstrates that the rate of convergence is greater for finite difference methods or the finite element method. For a single first-order hyperbolic equation, the method is seen as remaining stable even when the coefficient changes sign, although in this case it is specified that care must be taken to have adequate spatial resolution. It is noted that this fact, combined with the fact that collocation methods are easy to apply in the nonlinear case, shows that the pseudospectral method is in general preferable to the Galerkin or Tau methods.

Gottlieb, D.↗
NASA NTRS1984

Discrete approximation methods for parameter identification in delay systems

Approximation schemes for parameter identification problems in which the governing state equation is a linear functional differential equation of retarded type are constructed. The basis of the schemes is the replacement of the parameter identification problem having an infinite dimensional state equation by a sequence of approximating parameter identification problems in which the states are given by finite dimensional discrete difference equations. The difference equations are constructed using linear semigroup theory and rational function approximations to the exponential. Sufficient conditions are given for the convergence of solutions to the approximating problems, which can be obtained using conventional methods, to solutions to the original parameter identification problem. Finite difference and spline based schemes using Paderational function approximations to the exponential are constructed, and shown to satisfy the sufficient conditions for convergence. A discussion and analysis of numerical results obtained through the application of the schemes to several examples is included.

Rosen, I. G.↗
NASA NTRS2001

Prospects for Nuclear Electric Propulsion Using Closed-Cycle Magnetohydrodynamic Energy Conversion

Nuclear electric propulsion (NEP) has long been recognized as a major enabling technology for scientific and human exploration of the solar system, and it may conceivably form the basis of a cost-effective space transportation system suitable for space commerce. The chief technical obstacles to realizing this vision are the development of efficient, high-power (megawatt-class) electric thrusters and the development of low specific mass (less than 1 kg/kWe) power plants. Furthermore, comprehensive system analyses of multimegawatt class NEP systems are needed in order to critically assess mission capability and cost attributes. This Technical Publication addresses some of these concerns through a systematic examination of multimegawatt space power installations in which a gas-cooled nuclear reactor is used to drive a magnetohydrodynamic (MHD) generator in a closed-loop Brayton cycle. The primary motivation for considering MHD energy conversion is the ability to transfer energy out of a gas that is simply too hot for contact with any solid material. This has several intrinsic advantages including the ability to achieve high thermal efficiency and power density and the ability to reject heat at elevated temperatures. These attributes lead to a reduction in system specific mass below that obtainable with turbine-based systems, which have definite solid temperature limits for reliable operation. Here, the results of a thermodynamic cycle analysis are placed in context with a preliminary system analysis in order to converge on a design space that optimizes performance while remaining clearly within established bounds of engineering feasibility. MHD technology issues are discussed including the conceptual design of a nonequilibrium disk generator and opportunities for exploiting neutron-induced ionization mechanisms as a means of increasing electrical conductivity and enhancing performance and reliability. The results are then used to make a cursory examination of piloted Mars missions during the 2018 opportunity.

Litchford, R. J.↗