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 901 records · Page 50

All sources NASA NTRS DOE OSTI

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.↗

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.↗

NASA NTRS1987

Modeling and control of beam-like structures

The most popular finite element codes are based upon appealing theories of convergence of modal frequencies. For example, the popularity of cubic elements for beam-like structures is due to the rapid convergence of modal frequencies and stiffness properties. However, for those problems in which the primary consideration is the accuracy of response of the structure at specified locations it is more important to obtain accuracy in the modal costs than in the modal frequencies. The modal cost represents the contribution of a mode in the norm of the response vector. This paper provides a complete modal cost analysis for beam-like continua. Upper bounds are developed for mode truncation errors in the model reduction process and modal cost analysis dictates which modes to retain in order to reduce the model for control design purposes.

Hu, A.↗

NASA NTRS2014

3D Space Radiation Transport in a Shielded ICRU Tissue Sphere

A computationally efficient 3DHZETRN code capable of simulating High Charge (Z) and Energy (HZE) and light ions (including neutrons) under space-like boundary conditions with enhanced neutron and light ion propagation was recently developed for a simple homogeneous shield object. Monte Carlo benchmarks were used to verify the methodology in slab and spherical geometry, and the 3D corrections were shown to provide significant improvement over the straight-ahead approximation in some cases. In the present report, the new algorithms with well-defined convergence criteria are extended to inhomogeneous media within a shielded tissue slab and a shielded tissue sphere and tested against Monte Carlo simulation to verify the solution methods. The 3D corrections are again found to more accurately describe the neutron and light ion fluence spectra as compared to the straight-ahead approximation. These computationally efficient methods provide a basis for software capable of space shield analysis and optimization.

Wilson, John W.↗

NASA NTRS1983

An improved finite-difference analysis of uncoupled vibrations of tapered cantilever beams

An improved finite difference procedure for determining the natural frequencies and mode shapes of tapered cantilever beams undergoing uncoupled vibrations is presented. Boundary conditions are derived in the form of simple recursive relations involving the second order central differences. Results obtained by using the conventional first order central differences and the present second order central differences are compared, and it is observed that the present second order scheme is more efficient than the conventional approach. An important advantage offered by the present approach is that the results converge to exact values rapidly, and thus the extrapolation of the results is not necessary. Consequently, the basic handicap with the classical finite difference method of solution that requires the Richardson's extrapolation procedure is eliminated. Furthermore, for the cases considered herein, the present approach produces consistent lower bound solutions.

Subrahmanyam, K. B.↗

NASA NTRS1986

Analysis of the SOR iteration for the 9-point Laplacian

The SOR iteration for solving linear systems of equations depends upon an overrelaxation factor omega. A theory for determining omega was given by Young (1950) for consistently ordered matrices. Here the optimal omega is determined for the 9-point stencil for the model problem of Laplace's equation on a square. Several orderings of the equations are considered, including the natural rowwise and multicolor orderings, all of which lead to non-consistently ordered matrices, and two equivalence classes of orderings are found with different convergence behavior and optimal omega's. The results for the natural rowwise ordering are compared to those of Garabedian (1956) and it is explained why both results are, in a sense, correct, even though they differ. Also analyzed is a pseudo SOR method for the model problem and it is shown that it is not as effective as the SOR methods. Finally, the point SOR methods are compared to known results for line SOR methods for this problem.

Loyce M. Adams↗

NASA NTRS1980

High-Pressure Hot-Gas Self-Acting Floating Ring Shaft Seal for Liquid Rocket Turbopumps

Design analysis, detail design, fabrication, and experimental evaluation was performed on two self acting floating ring shaft seals for a rocket engine turbopump high pressure 24132500 n/sq m (3500 psig) hot gas 533 K 9500 F) high speed 3142 rad/sec (30000 rmp) turbine. The initial design used Rayleigh step hydrodynamic lift pads to assist in centering the seal ring with minimum rubbing contact. The final design used a convergent tapered bore to provide hydrostatic centering force. The Rayleigh step design was tested for 107 starts and 4.52 hours total. The leakage was satisfactory; however, the design was not acceptable due to excessive wear caused by inadequate centering force and failure of the sealing dam caused by erosion damage. The tapered bore design was tested for 370 starts and 15.93 hours total. Satisfactory performance for the required life of 7.5 hours per seal was successfully demonstrated.

Burcham, R. E.↗

NASA NTRS1993

Preconditioning for the Navier-Stokes equations with finite-rate chemistry

The preconditioning procedure for generalized finite-rate chemistry and the proper preconditioning for the one-dimensional Navier-Stokes equations are presented. Eigenvalue stiffness is resolved and convergence-rate acceleration is demonstrated over the entire Mach-number range from the incompressible to the hypersonic. Specific benefits are realized at low and transonic flow speeds. The extended preconditioning matrix accounts for thermal and chemical non-equilibrium and its implementation is explained for both explicit and implicit time marching. The effect of higher-order spatial accuracy and various flux splittings is investigated. Numerical analysis reveals the possible theoretical improvements from using proconditioning at all Mach numbers. Numerical results confirm the expectations from the numerical analysis. Representative test cases include flows with previously troublesome embedded high-condition-number regions.

Godfrey, Andrew G.↗

NASA NTRS1999

Engineering Overview of a Multidisciplinary HSCT Design Framework Using Medium-Fidelity Analysis Codes

An objective of the HPCC Program at NASA Langley has been to promote the use of advanced computing techniques to more rapidly solve the problem of multidisciplinary optimization of a supersonic transport configuration. As a result, a software system has been designed and is being implemented to integrate a set of existing discipline analysis codes, some of them CPU-intensive, into a distributed computational framework for the design of a High Speed Civil Transport (HSCT) configuration. The proposed paper will describe the engineering aspects of integrating these analysis codes and additional interface codes into an automated design system. The objective of the design problem is to optimize the aircraft weight for given mission conditions, range, and payload requirements, subject to aerodynamic, structural, and performance constraints. The design variables include both thicknesses of structural elements and geometric parameters that define the external aircraft shape. An optimization model has been adopted that uses the multidisciplinary analysis results and the derivatives of the solution with respect to the design variables to formulate a linearized model that provides input to the CONMIN optimization code, which outputs new values for the design variables. The analysis process begins by deriving the updated geometries and grids from the baseline geometries and grids using the new values for the design variables. This free-form deformation approach provides internal FEM (finite element method) grids that are consistent with aerodynamic surface grids. The next step involves using the derived FEM and section properties in a weights process to calculate detailed weights and the center of gravity location for specified flight conditions. The weights process computes the as-built weight, weight distribution, and weight sensitivities for given aircraft configurations at various mass cases. Currently, two mass cases are considered: cruise and gross take-off weight (GTOW). Weights information is obtained from correlations of data from three sources: 1) as-built initial structural and non-structural weights from an existing database, 2) theoretical FEM structural weights and sensitivities from Genesis, and 3) empirical as-built weight increments, non-structural weights, and weight sensitivities from FLOPS. For the aeroelastic analysis, a variable-fidelity aerodynamic analysis has been adopted. This approach uses infrequent CPU-intensive non-linear CFD to calculate a non-linear correction relative to a linear aero calculation for the same aerodynamic surface at an angle of attack that results in the same configuration lift. For efficiency, this nonlinear correction is applied after each subsequent linear aero solution during the iterations between the aerodynamic and structural analyses. Convergence is achieved when the vehicle shape being used for the aerodynamic calculations is consistent with the structural deformations caused by the aerodynamic loads. To make the structural analyses more efficient, a linearized structural deformation model has been adopted, in which a single stiffness matrix can be used to solve for the deformations under all the load conditions. Using the converged aerodynamic loads, a final set of structural analyses are performed to determine the stress distributions and the buckling conditions for constraint calculation. Performance constraints are obtained by running FLOPS using drag polars that are computed using results from non-linear corrections to the linear aero code plus several codes to provide drag increments due to skin friction, wave drag, and other miscellaneous drag contributions. The status of the integration effort will be presented in the proposed paper, and results will be provided that illustrate the degree of accuracy in the linearizations that have been employed.

Weston, R. P.↗

NASA NTRS2007

Direct Calculation of the Scattering Amplitude Without Partial Wave Decomposition: Inclusion of Correlation Effects - III

In the first two papers in this series, we developed a method for studying electron-hydrogen scattering that does not use partial wave analysis. We constructed an ansatz for the wave function in both the static and static exchange approximations and calculated the full scattering amplitude. Here we go beyond the static exchange approximation, and include correlation in the wave function via a modified polarized orbital. This correlation function provides a significant improvement over the static exchange approximation: the resultant elastic scattering amplitudes are in very good agreement with fully converged partial wave calculations for electron-hydrogen scattering. A fully variational modification of this approach is discussed in the conclusion of the article Popular summary of Direct calculation of the scattering amplitude without partial wave expansion. III ....." by J. Shertzer and A. Temkin. In this paper we continue the development of In this paper we continue the development of a new approach to the way in which researchers have traditionally used to calculate the scattering cross section of (low-energy) electrons from atoms. The basic mathematical problem is to solve the Schroedinger Equation (SE) corresponding the above physical process. Traditionally it was always the case that the SE was reduced to a sequence of one-dimensional (ordinary) differential equations - called partial waves which were solved and from the solutions "phase shifts" were extracted, from which the scattering cross section was calculated.

Shertzer, Janine↗

NASA NTRS2019

Regularization and Error Characterization of GRACE Mascons

We present a new global time-variable gravity mascon solution derived from Gravity Recovery and Climate Experiment (GRACE) Level 1B data. The new product from the NASA Goddard Space Flight Center (GSFC) results from a novel approach that combines an iterative solution strategy with geographical binning of inter-satellite range-acceleration residuals in the construction of time-dependent regularization matrices applied in the inversion of mascon parameters. This estimation strategy is intentionally conservative as it seeks to maximize the role of the GRACE measurements on the final solution while minimizing the influence of the regularization design process. We fully reprocess the Level 1B data in the presence of the final mascon solution to generate true post-fit inter-satellite residuals, which are utilized to confirm solution convergence and to validate the mascon noise uncertainties. We also present the mathematical case that regularized mascon solutions are biased, and that this bias, or leakage, must be combined with the estimated noise variance to accurately assess total mascon uncertainties. The estimated leakage errors are determined from the monthly resolution operators. We present a simple approach to compute the total uncertainty for both individual mascon and regional analysis of the GSFC mascon product, and validate the results in comparison with independent mascon solutions and calibrated Stokes uncertainties. Lastly, we present the new solution and uncertainties with global analyses of the mass trends and annual amplitudes, and compute updated trends for the global ocean, and the respective contributions of the Greenland Ice Sheet, Antarctic Ice Sheet, Gulf of Alaska, and terrestrial water storage. This analysis highlights the successful closure of the global mean sea level budget, that is, the sum of global ocean mass from the GSFC mascons and the steric component from Argo floats agrees well with the total determined from sea surface altimetry.

B. D. Loomis↗

NASA NTRS2019

Regularization and Error Characterization of GRACE Mascons

Loomis, B. D.↗

NASA NTRS1980

Electromagnetic scattering from two identical pseudospheres

In an attempt to gain some insight into the problem of how close particles have to be on the average before phase effects become important, the cross section and phase matrices for two identical particles as a function of the distance between their centers for both broadside and end-on illumination have been calculated. The point dipole approximation is employed to build two pseudospherical particles each consisting of 32 dipoles with an effective size parameter of 0.9283 and refractive index of 1.54-0.0i. In regard to cross sections, it is found that in certain cases several percent deviation from the converged value, which is twice the single particle cross section, could be found up to distances of separation of 25 single particle diameters. The maxima and minima in the cross section curves can be explained in terms of simple interference effects. Also the element P11 of the phase matrix showed a great deal of structure, and the position of maxima and minima could be explained by phase shift analysis.

Kattawar, G. W.↗