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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

PIV Measurements of Supersonic Internally-Mixed Dual-Stream Jets

While externally mixed, or separate flow, nozzle systems are most common in high bypass-ratio aircraft, they are not as attractive for use in lower bypass-ratio systems and on aircraft that will fly supersonically. The noise of such propulsion systems is also dominated by jet noise, making the study and noise reduction of these exhaust systems very important, both for military aircraft and future civilian supersonic aircraft. This paper presents particle image velocimetry of internally mixed nozzle with different area ratios between core and bypass, and nozzles that are ideally expanded and convergent. Such configurations independently control the geometry of the internal mixing layer and of the external shock structure. These allow exploration of the impact of shocks on the turbulent mixing layers, the impact of bypass ratio on broadband shock noise and mixing noise, and the impact of temperature on the turbulent flow field. At the 2009 AIAA/CEAS Aeroacoustics Conference the authors presented data and analysis from a series of tests that looked at the acoustics of supersonic jets from internally mixed nozzles. In that paper the broadband shock and mixing noise components of the jet noise were independently manipulated by holding Mach number constant while varying bypass ratio and jet temperature. Significant portions of that analysis was predicated on assumptions regarding the flow fields of these jets, both shock structure and turbulence. In this paper we add to that analysis by presenting particle image velocimetry measurements of the flow fields of many of those jets. In addition, the turbulent velocity data documented here will be very useful for validation of computational flow codes that are being developed to design advanced nozzles for future aircraft.

Bridges, James E.↗
NASA NTRSDate not supplied

Efficient Parametric Uncertainty Analysis of an Earth Entry Vehicle Concept Using Least Angle Regression

The objective of this work was to outline and apply an efficient and accurate parametric un-certainty propagation approach to the analysis of convective heating on an Earth entry vehicle concept. The described approach was based on Least Angle Regression used to solve a sparse and underdetermined linear system in the point-collocation non-intrusive polynomial chaos surrogate method. This approach involved an iterative process to computing the non-zero terms of the underlying polynomial chaos model using only enough samples to converge uncertainty interval predictions and Sobol index values based global nonlinear sensitivity estimates. The Earth entry vehicle was analyzed at three points along a representative trajectory for a Mars return mission. 329 sources of uncertainty were identified in the computational fluid dynamics model used to predict the forebody convective heating. These included uncertainty in flow field chemical rates, collision integrals, heats of formation, surface finite rate char model reaction rates, wall roughness height, and the turbulent Schmidt number. Results from this study showed that convective heating uncertainty as high as 50% of the nominal was predicted with only about 50 evaluations of the computational model. This was far fewer than would be required for a sampling-based approach or a full basis polynomial chaos model, which would have required over 50,000 samples. Additionally, results showed that over 90% of the total convective heating uncertainty was due to uncertainty in the N2catalytic rate on the surface, while the remainder of the uncertainty was attributed to the turbulent Schmidt number and the wall roughness uncertainties.

Thomas K West IV↗
NASA NTRS1981

Numerical analysis of spectral properties of coupled oscillator Schroedinger operators. I - Single and double well anharmonic oscillators

Several methods for computing many eigenvalues and eigenfunctions of a single anharmonic oscillator Schroedinger operator whose potential may have one or two minima are described. One of the methods requires the solution of an ill-conditioned generalized eigenvalue problem. This method has the virtue of using a bounded amount of work to achieve a given accuracy in both the single and double well regions. Rigorous bounds are given, and it is proved that the approximations converge faster than any inverse power of the size of the matrices needed to compute them. The results of computations for the g:phi(4):1 theory are presented. These results indicate that the methods actually converge exponentially fast.

Isaacson, D.↗
DOE OSTI2021

The halo light-cone catalogues of ABACUSSUMMIT

We describe a method for generating halo catalogues on the light-cone using the AbacusSummit suite of N-body simulations. The main application of these catalogues is the construction of realistic mock galaxy catalogues and weak lensing maps on the sky. Our algorithm associates the haloes from a set of coarsely spaced snapshots with their positions at the time of light-cone crossing by matching halo particles to on-the-fly light-cone particles. It then records the halo and particle information into an easily accessible product, which we call the AbacusSummit halo light-cone catalogues. Our recommended use of this product is in the halo mass regime of M halo > 2.1 x 10 11 M ⊙ h -1 for the base resolution simulations, i.e. haloes containing at least 100 particles, where the interpolated halo properties are most reliable. To test the validity of the obtained catalogues, we perform various visual inspections and consistency checks. In particular, we construct galaxy mock catalogues of emission-line galaxies (ELGs) at z ~ 1 by adopting a modified version of the AbacusHOD script, which builds on the standard halo occupation distribution (HOD) method by including various extensions. We find that the multipoles of the autocorrelation function are consistent with the predictions from the full-box snapshot, implicitly validating our algorithm. In addition, we compute and output CMB convergence maps and find that the auto- and cross-power spectrum agrees with the theoretical prediction at the sub-per-cent level.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Transitioning from Stage-III to Stage-IV: cosmology from galaxy×CMB lensing and shear×CMB lensing

ABSTRACT We examine the cosmological constraining power from two cross-correlation probes between galaxy and cosmic microwave background (CMB) surveys: the cross-correlation of lens galaxy density with CMB lensing convergence 〈δgκCMB〉, and source galaxy weak lensing shear with CMB lensing convergence 〈γκCMB〉. These two cross-correlation probes provide an independent cross-check of other large-scale structure constraints and are insensitive to galaxy-only or CMB-only systematic effects. In addition, when combined with other large-scale structure probes, the cross-correlations can break degeneracies in cosmological and nuisance parameters, improving both the precision and robustness of the analysis. In this work, we study how the constraining power of 〈δgκCMB〉 + 〈γκCMB〉 changes from Stage-III (ongoing) to Stage-IV (future) surveys. Given the flexibility in selecting the lens galaxy sample, we also explore systematically the impact on cosmological constraints when we vary the redshift range and magnitude limit of the lens galaxies using mock galaxy catalogs. We find that in our setup, the contribution to cosmological constraints from 〈δgκCMB〉 and 〈γκCMB〉 are comparable in the Stage-III data sets; but in Stage-IV surveys, the noise in 〈δgκCMB〉 becomes subdominant to cosmic variance, preventing 〈δgκCMB〉 to further improve the constraints. This implies that to maximize the cosmological constraints from future 〈δgκCMB〉 + 〈γκCMB〉 analyses, we should focus more on the requirements on 〈γκCMB〉 instead of 〈δgκCMB〉. Furthermore, the selection of the lens sample should be optimized in terms of our ability to characterize its redshift or galaxy bias instead of its number density.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1976

Electrostatic analysis of charge-coupled structures

An analysis is presented which is based upon a numerical solution of Gauss's law for the multiple dielectric semiconductor and electrode structure. The formulation resulted in a large set of equations, usually nonlinear, which were solved by an iterative scheme based upon the Gauss-Seidel method employing a relaxation parameter. Proper consideration of the nonlinearity of the space-charge in the semiconductor allowed convergent solutions which can give reasonable approximations to the channel potential profiles though not as good for estimating the space charge itself. The results of this work, which are three computer programs, are listed. None of the programs give the electric field explicitly as output data. However, when this information is desired, the programs can be expanded to include field estimations based on a polynominal fit to the potential field, or they may serve as subprograms for main programs dealing with general aspects of charge transfer.

Gassaway, J. D.↗
NASA NTRS1981

Software models: A Bayesian approach to parameter estimation in the Jelenski-Moranda software reliability model

Maximum likelihood estimation procedures for the Jelinski-Moranda software reliability model often give misleading answers. A reparameterization and a Bayesian analysis eliminate some of the problems incurred by MLE methods and often give better predictions on sets of real and simulated data. Practical difficulties in estimating the initial number of errors N and the failure rate of each error phi by the method of maximum likelihood are: N, the MLE of N, is occasionally infinite (i.e., the routines for calculating N and phi do not converge). It is shown that N is finite sub i only if the regression line of the interevent times t sub i vs. i has positive slope. A serious problem is that often N approximates n, the sample size, and sometimes N = n. Thus the MLE predicts that the program is perfect even when it is far from being so. Only when almost all failures have been removed can N and phi be trusted near the end of debugging.

Littlewood, B.↗
NASA NTRS1986

Optimization methods and silicon solar cell numerical models

The goal of this project is the development of an optimization algorithm for use with a solar cell model. It is possible to simultaneously vary design variables such as impurity concentrations, front junction depth, back junctions depth, and cell thickness to maximize the predicted cell efficiency. An optimization algorithm has been developed and interfaced with the Solar Cell Analysis Program in 1 Dimension (SCAPID). SCAPID uses finite difference methods to solve the differential equations which, along with several relations from the physics of semiconductors, describe mathematically the operation of a solar cell. A major obstacle is that the numerical methods used in SCAPID require a significant amount of computer time, and during an optimization the model is called iteratively until the design variables converge to the value associated with the maximum efficiency. This problem has been alleviated by designing an optimization code specifically for use with numerically intensive simulations, to reduce the number of times the efficiency has to be calculated to achieve convergence to the optimal solution. Adapting SCAPID so that it could be called iteratively by the optimization code provided another means of reducing the cpu time required to complete an optimization. Instead of calculating the entire I-V curve, as is usually done in SCAPID, only the efficiency is calculated (maximum power voltage and current) and the solution from previous calculations is used to initiate the next solution.

Girardini, K.↗
NASA NTRS1986

Optimization methods and silicon solar cell numerical models

An optimization algorithm for use with numerical silicon solar cell models was developed. By coupling an optimization algorithm with a solar cell model, it is possible to simultaneously vary design variables such as impurity concentrations, front junction depth, back junction depth, and cell thickness to maximize the predicted cell efficiency. An optimization algorithm was developed and interfaced with the Solar Cell Analysis Program in 1 Dimension (SCAP1D). SCAP1D uses finite difference methods to solve the differential equations which, along with several relations from the physics of semiconductors, describe mathematically the performance of a solar cell. A major obstacle is that the numerical methods used in SCAP1D require a significant amount of computer time, and during an optimization the model is called iteratively until the design variables converge to the values associated with the maximum efficiency. This problem was alleviated by designing an optimization code specifically for use with numerically intensive simulations, to reduce the number of times the efficiency has to be calculated to achieve convergence to the optimal solution.

Girardini, K.↗
NASA NTRS1986

Cloud and moisture fields derived from the GLA retrievals of HIRS2/MSU data

The GLA retrieval scheme for the analysis of HIRS and MSU radiances is applied to derive cloud and humidity fields from the HIRS2/MSU data for June 1979. For the retrieval of cloud fraction and cloud top pressure, the original algorithm of Susskind et al. (1983) and Susskind et al. (1984) was improved. The derived profiles of the monthly mean fields of cloud fraction and cloud top pressure clearly show the intertropical convergence zone, with the most intense convection in the monsoonal region of the southern Asia and over Central America, which show up as containing the highest cloud top levels and largest cloud amount. For the retrieval of humidity profiles, which are not one of the products of the original processing system, a new algorithm was derived.

Reuter, D.↗
NASA NTRS1993

Scattering by a groove in an impedance plane

An analysis of two-dimensional scattering from a narrow groove in an impedance plane is presented. The groove is represented by a impedance surface and the problem reduces to that of scattering from an impedance strip in an otherwise uniform impedance plane. On the basis of this model, appropriate integral equations are constructed using a form of the impedance plane Green's functions involving rapidly convergent integrals. The integral equations are solved by introducing a single basis representation of the equivalent current on the narrow impedance insert. Both transverse electric (TE) and transverse magnetic (TM) polarizations are treated. The resulting solution is validated by comparison with results from the standard boundary integral method (BIM) and a high frequency solution. It is found that the presented solution for narrow impedance inserts can be used in conjunction with the high frequency solution for the characterization of impedance inserts of any given width.

Bindiganavale, Sunil↗
NASA NTRS2010

Mineralogy of SNC Meteorite EET79001 by Simultaneous Fitting of Moessbauer Backscatter Spectra

We have acquired M ssbauer spectra for SNC meteorite EET79001 with a MIMOS II backscatter M ssbauer spectrometer [1] similar to those now operating on Mars as part of the Mars Exploration Rover (MER) missions. We are working to compare the Fe mineralogical composition of martian meteorites with in-situ measurements on Mars. Our samples were hand picked from the >1 mm size fraction of saw fines on the basis of lithology, color, and grain size (Table 1). The chips were individually analyzed at approx.300K by placing them on a piece of plastic that was in turn supported by the contact ring of the instrument (oriented vertically). Tungsten foil was used to mask certain areas from analysis. As shown in Figure 1, a variety of spectra was obtained, each resulting from different relative contributions of the Fe-bearing minerals present in the sample. Because the nine samples are reasonably mixtures of the same Fe-bearing phases in variable proportions, the nine spectra were fit simultaneously (simfit) with a common model, adjusting parameters to a single minimum chi-squared convergence criterion [2]. The starting point for the fitting model and values of hyperfine parameters was the work of Solberg and Burns [3], who identified olivine, pyroxene, and ferrous glass as major, and ilmenite and a ferric phase as minor (<5%), Fe-bearing phases in EET79001.

Morris, Richard V.↗
NASA NTRS2010

Predicting Tropical Cyclogenesis with a Global Mesoscale Model: Hierarchical Multiscale Interactions During the Formation of Tropical Cyclone Nargis(2008)

Very severe cyclonic storm Nargis devastated Burma (Myanmar) in May 2008, caused tremendous damage and numerous fatalities, and became one of the 10 deadliest tropical cyclones (TCs) of all time. To increase the warning time in order to save lives and reduce economic damage, it is important to extend the lead time in the prediction of TCs like Nargis. As recent advances in high-resolution global models and supercomputing technology have shown the potential for improving TC track and intensity forecasts, the ability of a global mesoscale model to predict TC genesis in the Indian Ocean is examined in this study with the aim of improving simulations of TC climate. High-resolution global simulations with real data show that the initial formation and intensity variations of TC Nargis can be realistically predicted up to 5 days in advance. Preliminary analysis suggests that improved representations of the following environmental conditions and their hierarchical multiscale interactions were the key to achieving this lead time: (1) a westerly wind burst and equatorial trough, (2) an enhanced monsoon circulation with a zero wind shear line, (3) good upper-level outflow with anti-cyclonic wind shear between 200 and 850 hPa, and (4) low-level moisture convergence.

Shen, B.-W.↗
NASA NTRS1984

Radar and surface data analysis of a microburst in JAWS

Two microbursts that occurred simultaneously and were produced by two neighboring storms on June 22, 1982 have been examined using Doppler radar analyses and PAM surface network wind observations. The development of the microbursts is described, and the probable importance of forcing mechanisms is considered. It is found that a sounding conductive to deep convection with a very dry, unstable boundary layer gives rise to hail storms which produce microbursts associated with high-reflectivity cores. The two microbursts were associated with intensifying shafts of very high reflectivities. Convergence into the area of downdraft within the dry layer just below the cloud base, evidence of low reflectivity notches in the neighborhood of the microburst downdrafts, and subsaturated surface conditions in the microburst outflows suggest active evaporation processes may be an important source of energy for these two microbursts.

Hjelmfelt, M. R.↗
DOE OSTI2021

Development of the reactor core isolation cooling system model for the extended station black-out accident

In order to accurately predict the behavior of the Reactor Core Isolation Cooling System (RCIC) in the context of the FLEX strategies under the extended station black-out accident (SBO), we present a mechanistic model of the RCIC in this work. The control volume formulation of angular momentum conservation is employed to model the Terry turbine in the RCIC system. The two-phase nozzle flow model developed in the previous work is utilized to simulate the convergent-divergent nozzles in the Terry turbine. A nozzle flow table is generated to cover the operation range of the RCIC during beyond design basis accidents (BDBAs). The mechanistic RCIC turbine-pump model and the nozzle flow table are incorporated in a TRACE BWR model. Simulations of short-term and long-term SBO transients are performed to demonstrate dynamic RCIC responses. A system analysis code with the proposed mechanistic model of the RCIC system will have a better prediction of system safety response during SBO accidents. With the addition of such more realistic models, the system code will be capable of providing expanded understanding of safety margin, which can help help to form the technical basis for the successful implementation of the FLEX and SAMG measures under Extended Loss of AC Power (ELAP) conditions for BWR or PWR reactor designs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Mathematical modelling of the concave front in the adjacent high explosive detonation problem

This study presents an analysis of the transition-zone in adjacent high explosive (HE) detonation problems which uses a $D, 𝜅, \dot{D}$ relationship, where $D$ is the detonation front-normal velocity, 𝜅 is the detonation front curvature and $\dot{D}$ is the time derivative of detonation front-normal velocity. Our approach extends the traditional $(D, 𝜅)$ model to accurately predict the behaviour of both diverging and converging detonation shock fronts. Our findings affirm that a hyperbolic type of front evolution equation, enhanced with wave acceleration, provides a robust framework for modelling complex shock front dynamics in HE materials. This approach not only captures the natural effects of straightness and boundary slope jumps in the transition-zone but also bridges the gap between mathematical predictions and experimental observations, offering insights into the behaviour of both diverging and converging detonation propagations in a homogeneous HE.

acceleration↗
DOE OSTI2020

Atmospheric condition identification in multivariate data through a metric for total variation

Identification of atmospheric conditions within a multivariable atmospheric data set is a necessary step in the validation of emerging and existing high-fidelity models used to simulate wind plant flows and operation.Atmospheric conditions relevant for wind energy research include stationary conditions, given the need for well-converged statistics for model validation, as well as conditions observed less frequently, such as extreme atmospheric events, which are used in wind turbine and wind plant design.Aggregation of observations without regard to covariance between time series discounts the dynamical nature of the atmosphere and is not sufficiently representative of atmospheric conditions.Identification and characterization of continuous time periods with atmospheric conditions that have a high value for analysis or simulation set the stage for more advanced model validation and the development of real-time control and operational strategies.The current work explores a single metric for variation in a multivariate data sample that quantifies variability within each channel as well as covariance between channels.The total variation is used to identify conditions of interest that conform to desired objective functions, such as stationary conditions, ramps or waves of wind speed, and changes in wind direction.Total variation is somewhat sensitive to the presence of outliers in the input data, and the method is best complemented by quality-control procedures to ensure reliable results.The direct detection and classification of events or conditions of interest within atmospheric data sets is vital to developing our understanding of wind plant response and to the formulation of forecasting and control models.

17 WIND ENERGY↗
NASA NTRS1971

The pressure and deformation profile between two colliding lubricated cylinders

The pressure and deformation profiles between two colliding lubricated cylinders are obtained by solving the coupled, time-dependent elastohydrodynamic equations with an iterative procedure. The analysis includes effects which were not considered in a previous solution, namely, the effect of the lubricant compressibility and the effect of a lubricant with composite pressure-viscosity coefficients. It is found that the local approach velocity plays an important role during final stages of normal approach. It causes the lubricant to be entrapped within the contact region; neither the pressure nor the deformation profile converges to the Hertzian profile for a dry contact. The use of a smaller pressure-viscosity coefficient at high pressures reduces the sharp pressure gradient at the center of the contact and produces a much milder variation of load with respect to the film thickness. The effect of compressibility of the lubricant is found to be relatively small.

Lee, K.↗