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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Low-Density Nozzle Flow by the Direct Simulation Monte Carlo and Continuum Methods

Two different approaches, the direct simulation Monte Carlo (DSMC) method based on molecular gasdynamics, and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gasdynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental data, in which density and rotational temperature variations along the centerline and at various locations inside a low-density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data, both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity-temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chang-Hong↗

NASA NTRS2001

An Investigation of the Influence of Urban Areas on Rainfall Using a Cloud-Mesoscale Model and the TRMM Satellite

A recent paper by Shepherd and Pierce (conditionally accepted to Journal of Applied Meteorology) used rainfall data from the Precipitation Radar on NASA's Tropical Rainfall Measuring Mission's (TRMM) satellite to identify warm season rainfall anomalies downwind of major urban areas. A convective-mesoscale model with extensive land-surface processes is employed to (a) determine if an urban heat island (UHI) thermal perturbation can induce a dynamic response to affect rainfall processes and (b) quantify the impact of the following three factors on the evolution of rainfall: (1) urban surface roughness, (2) magnitude of the UHI temperature anomaly, and (3) physical size of the UHI temperature anomaly. The sensitivity experiments are achieved by inserting a slab of land with urban properties (e.g. roughness length, albedo, thermal character) within a rural surface environment and varying the appropriate lower boundary condition parameters. Early analysis suggests that urban surface roughness (through turbulence and low-level convergence) may control timing and initial location of UHI-induced convection. The magnitude of the heat island appears to be closely linked to the total rainfall amount with minor impact on timing and location. The physical size of the city may predominantly impact on the location of UHI-induced rainfall anomaly. The UHI factor parameter space will be thoroughly investigated with respect to their effects on rainfall amount, location, and timing. This study extends prior numerical investigations of the impact of urban surfaces on meteorological processes, particularly rainfall development. The work also contains several novel aspects, including the application of a high-resolution (less than I km) cloud-mesoscale model to investigate urban-induce rainfall process; investigation of thermal magnitude of the UHI on rainfall process; and investigation of UHI physical size on rainfall processes.

Shepherd, J. Marshall↗

NASA NTRS2006

Overview of the Development for a Suite of Low-Thrust Trajectory Analysis Tools

A NASA intercenter team has developed a suite of low-thrust trajectory analysis tools to make a significant improvement in three major facets of low-thrust trajectory and mission analysis. These are: 1) ease of use, 2) ability to more robustly converge to solutions, and 3) higher fidelity modeling and accuracy of results. Due mostly to the short duration of the development, the team concluded that a suite of tools was preferred over having one integrated tool. This tool-suite, their characteristics, and their applicability will be described. Trajectory analysts can read this paper and determine which tool is most appropriate for their problem.

Kos, Larry D.↗

NASA NTRS2018

Effect of Distributed Patch of Smooth Roughness Elements on Transition in a High-Speed Boundary Layer

Surface roughness is known to have a substantial impact on the aerothermodynamic loading of hypersonic vehicles, particularly via its influence on the laminar-turbulent transition process within the boundary layer. Numerical simulations are performed to investigate the effects of a distributed region of densely packed, smooth-shaped roughness elements on the laminar boundary layer over a 7-degree half-angle, circular cone for flow conditions corresponding to a selected trajectory point from the ascent phase of the HIFiRE-1 flight experiment. For peak-to-valley roughness heights of 50 percent or less in comparison with the thickness of the unperturbed boundary layer, the computations converge to a stationary flow, suggesting that the flow is globally stable. Analysis of convective instabilities in the wake of the roughness patch indicates two dominant families of unstable disturbances, namely, a high frequency mode that corresponds to Mack mode waves modified by the wake and a lower frequency mode that corresponds to shear layer instabilities associated with the streaks in the roughness wake. Even though the peak growth rate of the later mode is more than 35 percent greater than the peak growth rates of the Mack modes, the latter modes achieve higher amplification ratios, and hence, are likely to dominate the onset of transition, which is estimated to occur slightly later than that in the unperturbed, i.e., smooth surface boundary layer. Additional computations are performed to investigate the effects of various roughness patch configurations on a Mach 3.5 flat plate boundary layer, to help guide an upcoming experiment in the Mach 3.5 Supersonic Low Disturbance Tunnel at NASA Langley Research Center. In this case, the cumulative reinforcement of basic state distortion over the length of the roughness patch is predicted to yield a significantly earlier transition than that over a smooth plate or a plate with a shorter length roughness patch.

Meelan Choudhari↗

NASA NTRS2022

Anthropogenic Sulfate Aerosol Pollution in South and East Asia Induces Increased Summer Precipitation Over Arid Central Asia

Precipitation has increased across the arid Central Asia region over recent decades. However, the underlying mechanisms of this trend are poorly understood. Here, we analyze multi-model simulations from the Precipitation Driver and Response Model Intercomparison Project (PDRMIP) to investigate potential drivers of the observed precipitation trend. We find that anthropogenic sulfate aerosols over remote polluted regions in South and East Asia lead to increased summer precipitation, especially convective and extreme precipitation, in arid Central Asia. Elevated concentrations of sulfate aerosols over remote polluted Asia cause an equatorward shift of the Asian Westerly Jet Stream through a fast response to cooling of the local atmosphere at mid-latitudes. This shift favours moisture supply from low-latitudes and moisture flux convergence over arid Central Asia, which is confirmed by a moisture budget analysis. High levels of absorbing black carbon lead to opposing changes in the Asian Westerly Jet Stream and reduced local precipitation, which can mask the impact of sulfate aerosols. This teleconnection between arid Central Asia precipitation and anthropogenic aerosols in remote Asian polluted regions highlights long-range impacts of anthropogenic aerosols on atmospheric circulations and the hydrological cycle.

Atmospheric dynamics↗

NASA NTRS1974

Vibration and stress analysis of soft-bonded shuttle insulation tiles. Modal analysis with compact widely space stringers

An efficient iterative procedure is described for the vibration and modal stress analysis of reusable surface insulation (RSI) of multi-tiled space shuttle panels. The method, which is quite general, is rapidly convergent and highly useful for this application. A user-oriented computer program based upon this procedure and titled RESIST (REusable Surface Insulation Stresses) has been prepared for the analysis of compact, widely spaced, stringer-stiffened panels. RESIST, which uses finite element methods, obtains three dimensional tile stresses in the isolator, arrestor (if any) and RSI materials. Two dimensional stresses are obtained in the tile coating and the stringer-stiffened primary structure plate. A special feature of the program is that all the usual detailed finite element grid data is generated internally from a minimum of input data. The program can accommodate tile idealizations with up to 850 nodes (2550 degrees-of-freedom) and primary structure idealizations with a maximum of 10,000 degrees-of-freedom. The primary structure vibration capability is achieved through the development of a new rapid eigenvalue program named ALARM (Automatic LArge Reduction of Matrices to tridiagonal form).

Ojalvo, I. U.↗

NASA NTRS1991

Blade design and analysis using a modified Euler solver

An iterative method for blade design based on Euler solver and described in an earlier paper is used to design compressor and turbine blades providing shock free transonic flows. The method shows a rapid convergence, and indicates how much the flow is sensitive to small modifications of the blade geometry, that the classical iterative use of analysis methods might not be able to define. The relationship between the required Mach number distribution and the resulting geometry is discussed. Examples show how geometrical constraints imposed upon the blade shape can be respected by using free geometrical parameters or by relaxing the required Mach number distribution. The same code is used both for the design of the required geometry and for the off-design calculations. Examples illustrate the difficulty of designing blade shapes with optimal performance also outside of the design point.

Leonard, O.↗

DOE OSTI2025

Bias-Variance Trade-Off in Physics-Informed Neural Networks with Randomized Smoothing for High-Dimensional PDEs

Physics-Informed Neural Networks (PINNs) have triggered a paradigm shift in scientific computing, leveraging mesh-free properties and robust approximation capabilities. While proving effective for low-dimensional partial differential equations (PDEs), the computational cost of PINNs remains a hurdle in high-dimensional scenarios. This is particularly pronounced when computing high-order and high-dimensional derivatives in the physics-informed loss. Randomized Smoothing PINN (RS-PINN) introduces Gaussian noise for stochastic smoothing of the original neural net model, enabling the use of Monte Carlo methods for derivative approximation, which eliminates the need for costly automatic differentiation. Despite its computational efficiency, especially in the approximation of high-dimensional derivatives, RS-PINN introduces biases in both loss and gradients, negatively impacting convergence, especially when coupled with stochastic gradient descent (SGD) algorithms. We present a comprehensive analysis of biases in RS-PINN, attributing them to the nonlinearity of the Mean Squared Error (MSE) loss as well as the intrinsic nonlinearity of the PDE itself. We propose tailored bias correction techniques, delineating their application based on the order of PDE nonlinearity. The derivation of an unbiased RS-PINN allows for a detailed examination of its advantages and disadvantages compared to the biased version. Specifically, the biased version has a lower variance and runs faster than the unbiased version, but it is less accurate due to the bias. To optimize the bias-variance trade-off, we combine the two approaches in a hybrid method that balances the rapid convergence of the biased version with the high accuracy of the unbiased version. In addition to methodological contributions, we present an enhanced implementation of RS-PINN. Extensive experiments on diverse high-dimensional PDEs, including Fokker-Planck, Hamilton-Jacobi-Bellman (HJB), viscous Burgers’, Allen-Cahn, and Sine-Gordon equations, illustrate the bias-variance trade-off and highlight the effectiveness of the hybrid RS-PINN. Empirical guidelines are provided for selecting biased, unbiased, or hybrid versions, depending on the dimensionality and nonlinearity of the specific PDE problem.

97 MATHEMATICS AND COMPUTING↗

NASA NTRSDate not supplied

Boundary-Layer Transition Prediction Through Loose Coupling of OVERFLOW and LASTRAC

Transition prediction based on linear stability theory is expected to more accurately reflect the causality of transition onset than phenomenological transition models based on RANS-like transport equations. To help achieve the CFD vision 2030 aim of building a CFD tool chain with automated prediction of boundary layer transition, a technique to loosely tie the NASA OVERFLOW CFD solver with the LASTRAC stability analysis tool is described. The coupled solver is then used to compute transition over selected over a flat plate in a freestream with sufficiently low levels of turbulence, NLF(1)-0416 airfoil, the 6:1 prolate spheroid at an angle of attack, and a NASA juncture flow model with symmetric wing configuration. The findings show that the loosely coupled approach can reliably predict the transition location accurately in scenarios that are dominated by a single transition mechanism involving Tollmien-Schlichting instabilities, crossflow instabilities, or separation bubble-induced transition, or include a mixture of selected mechanisms. The toolset presents here appears to be robust to the prescription of the initial transition location, and it can lead to a converged solution in four or five rounds of the mean flow calculation and stability analysis, with minimal input from the user.

boundary layer transition↗

NASA NTRS2019

Structural Analysis of Test Flight Vehicles for Application of Multifunctional Energy Storage System

Under the NASA Aeronautics Research Mission Directorate (ARMD) Convergent Aeronautical Solutions (CAS) project, NASA Glenn Research Center has been leading Multifunctional Structures for High Energy Lightweight Load-bearing Storage (M-SHELLS) research efforts. The technology of integrating load-carrying structures with electrical energy storage capacity has the potential to reduce the overall weight of future electric aircraft. The proposed project goals were to develop M-SHELLS in the form of honeycomb coupons and subcomponents, integrate them into the structure, and conduct low-risk flight-tests onboard a remotely piloted small aircraft. Experimental M-SHELLS energy-storing coupons were fabricated and tested in the laboratory for their electrical and mechanical properties. In this report, finite element model development and structural analyses of two small test aircraft candidates are presented. The finite element analysis of the initial two-spar wing is described for strain, deflection, and weight estimation. After a test aircraft Tempest was acquired, a load-deflection test of the wing was conducted. A finite element model of the Tempest was then developed based on the test aircraft dimensions and construction detail. The component weight analyses from the finite element model and test measurements were correlated. Structural analysis results with multifunctional energy storage panels in the fuselage of the test vehicle are presented. Although the flight test was cancelled because of programmatic reasons and time constraints, the structural analysis results indicate that the mid-fuselage floor composite panel could provide structural integrity with minimal weight penalty while supplying electrical energy. To explore potential future applications of the multifunctional structure, analyses of the NASA X-57 Maxwell electric aircraft and a NASA N+3 Technology Conventional Configuration (N3CC) fuselage are presented. Secondary aluminum structures in the fuselage sub-floor and cargo area were partially replaced with reinforced five-layer composite panels with M-SHELLS honeycomb core. The N3CC fuselage weight reduction associated with each design without risking structural integrity are described. The structural analysis and weight estimation with the application of composite M-SHELLS panels to the N3CC fuselage indicate a 3.2% reduction in the fuselage structural weight, prior to accounting for the additional weight of core material required to complete the energy storage functionality.

Mukhopadhyay, Vivek↗

DOE OSTI2024

Ensuring Solution Uniqueness in Three-Phase Power System State Estimation

This paper is concerned with the issue of potential non-unique solutions in three-phase state estimation. Theory of observability analysis for positive sequence power system state estimation is based on certain assumptions that avoid possibility of multiple solutions. Also, it is shown that observability of a positive sequence network remains independent of the network parameters or the operating state. When extending single-phase observability analysis directly to the three-phase case, this paper considers the possibility of converging to multiple solutions, i.e. solution non-uniqueness, even for cases where state estimator successfully converges. The study illustrates via numerical examples the likelihood of converging to entirely different solutions for certain network parameters. It also examines how the operating state, particularly under unbalanced loading, leads to solution non-uniqueness. The paper then describes an alternative approach to ensure a unique solution in three-phase state estimation. This method aims to accurately and uniquely estimate the state of any unbalanced three-phase system, irrespective of load imbalance, network configuration, existence of synchronous generators or transformers.

Power System State Estimation, Three-Phase, Distri↗

NASA NTRS1997

Atomic, Molecular, and Optical Physics: Optical Excitation Function of H(1s-2p) Produced by electron Impact from Threshold to 1.8 keV

The optical excitation function of prompt Lyman-Alpha radiation, produced by electron impact on atomic hydrogen, has been measured over the extended energy range from threshold to 1.8 keV. Measurements were obtained in a crossed-beams experiment using both magnetically confined and electrostatically focused electrons in collision with atomic hydrogen produced by an intense discharge source. A vacuum-ultraviolet mono- chromator system was used to measure the emitted Lyman-Alpha radiation. The absolute H(1s-2p) electron impact excitation cross section was obtained from the experimental optical excitation function by normalizing to the accepted optical oscillator strength, with corrections for polarization and cascade. Statistical and known systematic uncertainties in our data range from +/- 4% near threshold to +/- 2% at 1.8 keV. Multistate coupling affecting the shape of the excitation function up to 1 keV impact energy is apparent in both the present experimental data and present theoretical results obtained with convergent close- coupling (CCC) theory. This shape function effect leads to an uncertainty in absolute cross sections at the 10% level in the analysis of the experimental data. The derived optimized absolute cross sections are within 7% of the CCC calculations over the 14 eV-1.8 keV range. The present CCC calculations converge on the Bethe- Fano profile for H(1s-2p) excitation at high energy. For this reason agreement with the CCC values to within 3% is achieved in a nonoptimal normalization of the experimental data to the Bethe-Fano profile. The fundamental H(1s-2p) electron impact cross section is thereby determined to an unprecedented accuracy over the 14 eV - 1.8 keV energy range.

James, G. K.↗

DOE OSTI2024

Accelerating Multivariate Functional Approximation Computation with Domain Decomposition Techniques⋆

Modeling large datasets through Multivariate Functional Approximations (MFA) provide an elegant way to handle many visualization and scientific analysis workflows. The process necessitates scalable data partitioning methods to compute MFA representations efficiently without compromising the accuracy or continuity of the reconstructed solution. We propose a domain -decomposed method for computing the MFA with B -spline bases, which reduces the total work per task and uses a restricted Additive Schwarz (RAS) method to converge the control point data degrees -of -freedom along subdomain boundaries. We provide an in-depth analysis of the parallel approach with domain decomposition solvers, aiming to minimize local subdomain error residuals and recover high -order continuity at subdomain interfaces with appropriate choices of knot overlaps. The communication cost, determined by the overlap regions in the RAS implementation, is optimized to recover the numerical error profile of the single subdomain case. Our proposed method stands in contrast to previous methods, which typically only recover either C 0 or at best C 1 continuity for arbitrary B -spline degree expansions, or those that require post -processing to blend discontinuities in the reconstructed data. We demonstrate the effectiveness of our approach using analytical and real -world datasets in 1D, 2D, and 3D through both strong and weak scaling studies. The performance results indicate that the overall cost of computing the approximation is directly proportional to the underlying nearest -neighbor communication implementation, and is only weakly dependent on the overlap region size that determines the size of the messages. This finding underscores the efficiency and scalability of our proposed method, making it a promising solution for handling large datasets in scientific workflows.

additive Schwarz solvers↗

DOE OSTI2022

A posteriori superlinear convergence bounds for block conjugate gradient

In this paper, we extend to the block case the a posteriori bound showing superlinear convergence of the conjugate gradient method developed by van der Vorst and Vuik in [J. Comput. Applied Math., 48 (1993), pp. 327–341]. That is, we obtain similar bounds but now for the block conjugate gradient method. We also present a series of computational experiments, illustrating the validity of the bound developed here as well as the bound by Simoncini and Szyld from [SIAM Review, 47 (2005), pp. 247–272] using angles between subspaces. Using these bounds, we make some observations on the onset of superlinearity and how this onset depends on the eigenvalue distribution and the block size.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1996

Relating the Hadamard Variance to MCS Kalman Filter Clock Estimation

The Global Positioning System (GPS) Master Control Station (MCS) currently makes significant use of the Allan Variance. This two-sample variance equation has proven excellent as a handy, understandable tool, both for time domain analysis of GPS cesium frequency standards, and for fine tuning the MCS's state estimation of these atomic clocks. The Allan Variance does not explicitly converge for the nose types of alpha less than or equal to minus 3 and can be greatly affected by frequency drift. Because GPS rubidium frequency standards exhibit non-trivial aging and aging noise characteristics, the basic Allan Variance analysis must be augmented in order to (a) compensate for a dynamic frequency drift, and (b) characterize two additional noise types, specifically alpha = minus 3, and alpha = minus 4. As the GPS program progresses, we will utilize a larger percentage of rubidium frequency standards than ever before. Hence, GPS rubidium clock characterization will require more attention than ever before. The three sample variance, commonly referred to as a renormalized Hadamard Variance, is unaffected by linear frequency drift, converges for alpha is greater than minus 5, and thus has utility for modeling noise in GPS rubidium frequency standards. This paper demonstrates the potential of Hadamard Variance analysis in GPS operations, and presents an equation that relates the Hadamard Variance to the MCS's Kalman filter process noises.

Hutsell, Steven T.↗

DOE OSTI2021

Excited State Intramolecular Proton Transfer with Nuclear-Electronic Orbital Ehrenfest Dynamics

The recent development of the Ehrenfest dynamics approach in the nuclear-electronic orbital (NEO) framework provides a promising way to simulate coupled nuclear-electronic dynamics. Our previous study showed that the NEO-Ehrenfest approach with a semi- classical traveling proton basis method yields accurate predictions of molecular vibrational frequencies. In this work, we provide a more thorough analysis of the semiclassical traveling proton basis method to elucidate its validity and convergence behavior. We also conduct NEO-Ehrenfest dynamics simulations to study an excited state intramolecular proton transfer process. These simulations reveal that nuclear quantum effects influence the predictions of proton transfer reaction rates and kinetic isotope effects due to the intrinsic delocalized nature of the quantum nuclear wave function. This work illustrates the importance of nuclear quantum effects in coupled nuclear-electronic dynamical processes and shows that the NEO-Ehrenfest approach can be a powerful tool to provide insights and predictions for these processes.

Zhao, Luning↗

DOE OSTI2021

Predicting nonresonant pressure-driven MHD modes in equilibria with low magnetic shear

Nonresonant internal modes can be difficult to anticipate as there is no resonant surface in the plasma. However, equilibria that are unstable to multiple nonresonant magnetohydrodynamic (MHD) modes may be more prone to global loss of confinement since these instabilities generate spatially extended linear displacements, potentially enhancing magnetic field line chaos via nonlinear interactions. Here, we successfully predict the unstable nonresonant pressure-driven modes for equilibria with zero shear in the plasma core, irrational q on axis, and a central pressure gradient, which is consistent with pre-crash profiles in sawtoothing tokamak plasmas in the large-aspect-ratio limit. A criterion for identifying nonresonant modes most likely to be unstable is developed from the convergents of the continued fraction representation of q 0 . A higher-order analysis of the standard Energy Principle reveals the conditions under which these modes are expected to dominate. Linear growth rate spectra, as a function of toroidal mode number (up to n = 30), calculated using the initial-value extended-MHD code, M3D-C1, recover the characteristic dependence observed for ideal infernal modes. Nonresonant modes have also been invoked in some ideal sawtooth crash models. Finally, this work provides a mechanism to predict the mode numbers of infernal modes and, potentially, the width of some post-sawtooth-crash profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Provable Convergence of Plug-and-Play Priors With MMSE Denoisers

Plug-and-play priors (PnP) is a methodology for regularized image reconstruction that specifies the prior through an image denoiser. While PnP algorithms are well understood for denoisers performing maximum a posteriori probability (MAP) estimation, they have not been analyzed for the minimum mean squared error (MMSE) denoisers. Here we address this gap by establishing the first theoretical convergence result for the iterative shrinkage/thresholding algorithm (ISTA) variant of PnP for MMSE denoisers. We show that the iterates produced by PnP-ISTA with an MMSE denoiser converge to a stationary point of some global cost function. We validate our analysis on sparse signal recovery in compressive sensing by comparing two types of denoisers, namely the exact MMSE denoiser and the approximate MMSE denoiser obtained by training a deep neural net.

97 MATHEMATICS AND COMPUTING↗