Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

A method for solving the transonic full-potential equation for general configurations

A method is developed for solving the full-potential equation for two-dimensional and axisymmetric flow which retains the grid and boundary condition simplicity of the transonic small-disturbance codes. The method is based on a finite-volume formulation of the mass conservation equation in a Cartesian coordinate system, and is an extension of the method of Purvis and Burkhalter (1979). This finite-volume approach, combined with the simple boundary treatment, is shown to result in a highly robust method applicable to a wide range of geometries and flow conditions. The accuracy of the method is demonstrated for general geometries in two-dimensional and axisymmetric flows. The use of this method results in significant gains in convergence rate over the vertical-line over-relaxation scheme by incorporating an AF2-type algorithm (Ballhaus et al., 1978). It is suggested that the simplicity of this method shold allow a relatively easy extension to complex geometries in three-dimensional flows, and complex two-dimensional configurations such as multielement airfoils should be amenable to this method.

Wedan, B.↗
NASA NTRS1984

Improved methods of vibration analysis of pretwisted, airfoil blades

Vibration analysis of pretwisted blades of asymmetric airfoil cross section is performed by using two mixed variational approaches. Numerical results obtained from these two methods are compared to those obtained from an improved finite difference method and also to those given by the ordinary finite difference method. The relative merits, convergence properties and accuracies of all four methods are studied and discussed. The effects of asymmetry and pretwist on natural frequencies and mode shapes are investigated. The improved finite difference method is shown to be far superior to the conventional finite difference method in several respects. Close lower bound solutions are provided by the improved finite difference method for untwisted blades with a relatively coarse mesh while the mixed methods have not indicated any specific bound.

Subrahmanyam, K. B.↗
NASA NTRS1985

An evaluation of fracture analysis methods

The results of an experimental and predictive round robin on the applications of fracture analysis methods are presented. The objective of the round robin was to verify whether fracture analysis methods currently in use can or cannot predict failure loads on complex structural components containing cracks. Fracture results from tests on a number of compact specimens were used to make the predictions. The accuracy of the prediction methods was evaluated in terms of the variation in the ratio of predicted to experimental failure loads, and the predictions methods are ranked in order of minimum standard error. A range of applicability of the different methods was also considered in assessing their usefulness. For 7075-T651 aluminum alloy, the best methods were: the effective K sub R curve; the critical crack-tip opening displacement (CTOD) criterion using a finite element analysis; and the K sub R curve with the Dugdale model. For the 2024-T351 aluminum alloy, the best methods included: the two-parameter fracture criterion (TPFC); the CTOD parameter using finite element analysis; the K-curve with the Dugdale model; the deformation plasticity failure assessment diagram (DPFAD); and the effective K sub R curve with a limit load condition. For 304 stainless steel, the best methods were the limit load analysis; the CTOD criterion using finite-element analysis TPFC and DPFAD. Some sample experimental results are given in an appendix.

Newman, J. C., Jr.↗
NASA NTRS1987

A finite-difference outer layer and integral inner layer method for the solution of the turbulent boundary layer equations

A new turbulent boundary-layer method is developed which models the inner region with the law of the wall while the outer region uses Clauser's eddy viscosity in Matsuno's finite-difference method. The match point between the inner and outer regions as well as the wall shear stress are determined at each marching step during the computation. Results obtained for incompressible, two-dimensional flow over flat plates and ellipses are compared with solutions from a baseline method which uses a finite-difference method for the entire boundary layer. Since the present method used the finite-difference method in the outer region only, the number of grid points required was about half that needed for the baseline method. Accurate displacement and momentum thicknesses were predicted for all cases. Skin friction was predicted well for the flat plate, but the accuracy decreased significantly for the ellipses. Adding a wake functions to the law of the wall allows some of the pressure gradient effect to be taken into account thereby increasing the accuracy of the method.

Barnwell, R. W.↗
NASA NTRS1986

A comparison of the bounded derivative and the normal-mode initialization methods using real data

Application of the bounded-derivative and normal-mode methods to a simple linear barotropic model at a typical middle latitude shows that the two methods lead to identical constraints up to a certain degree of approximation. Beyond this accuracy the two methods may differ from each other. When applied to a global nonlinear barotropic model using real data, again the two methods lead to similar balanced initial states. The gravity oscillations in the unbalanced height field, which have amplitudes of up to 60 m with a dominant periodicity of about 5 to 6 h, are practically eliminated by both initialization methods. The rotational wind component is smooth even for the unbalanced initial state. The small-scale spatial features of the irrotational wind component are drastically reduced by initialization. Both the nonlinear normal-mode and the bounded-derivative initialization methods yield similar divergence fields centered around the areas of highest orography. The comparison shows that there is no significant loss of information in the mass and momentum fields, despite the fact that the bounded-derivative method employs only the original, rotational wind component to construct a balanced initial state compared to the normal-mode method, which, in addition, makes use of the unbalanced divergent wind and height fields.

Semazzi, F. H. M.↗
NASA NTRS1988

Improved method for stress and compatibility analysis of multicomponent rotating systems

An improved method of analyzing multicomponent rotating assemblies for the determination of operating stresses and component compatibility has been developed. In this method, a single finite element model is developed which contains all of the separate components in the rotating assembly. This is made possible by using gap elements to simulate the contact surfaces between components. The MARC finite element computer program is then used to perform the analysis. This improved method is less time consuming and more reliable than the conventional method of analyzing such systems. Results are presented for two different stress-compatibility analyses of a six-component axial flow compressor rotor. The results for the previously used flexibility analysis method are compared with those for the improved analysis method. The stresses predicted by each method compare quite well with each other. The predictions of the component compatibility, as well as the magnitude of the forces at the contact surfaces, also compare well for these two analysis procedures. It is therefore recommended that the improved analysis method be used to determine the stress-compatibility characterstics of multicomponent rotating systems.

Carek, Gerald A.↗
NASA NTRS1988

Velocity-pressure integrated versus penalty finite element methods for high Reynolds number flows

Velocity-pressure integrated and consistent penalty finite element computations of high Reynolds number, laminar flows are presented. In both of the methods, the pressure has been interpolated using linear shape functions for a triangular element. The triangular element is contained inside the bi-quadratic isoparametric element. It has been reported previously that the pressure interpolation method, when used in the velocity-pressure integrated method, yielded accurate computational results for high Reynolds number flows. It is shown that use of the same pressure interpolation method in the consistent penalty finite element method yielded accurate velocity and pressure fields which were comparable to those obtained using the velocity-pressure integrated method. Accuracy of the two finite element methods has been demonstrated by comparing the computational results with available experimental data and/or fine-grid finite difference computational results. Advantages and disadvantages of the two methods are discussed on the basis of accuracy and convergence nature. Example problems considered include a lid-driven cavity flow for Reynolds number of 10,000, a laminar backward-facing step flow, a laminar flow through a nest of cylinders, and a channel flow with an internal blockage. A finite element computer program (NSFLOW/P) for the 2-D, incompressible Navier-Stokes equations is also presented.

Kim, Sang-Wook↗
NASA NTRS1989

An investigation of new methods for estimating parameter sensitivities

The method proposed for estimating sensitivity derivatives is based on the Recursive Quadratic Programming (RQP) method and in conjunction a differencing formula to produce estimates of the sensitivities. This method is compared to existing methods and is shown to be very competitive in terms of the number of function evaluations required. In terms of accuracy, the method is shown to be equivalent to a modified version of the Kuhn-Tucker method, where the Hessian of the Lagrangian is estimated using the BFS method employed by the RQP algorithm. Initial testing on a test set with known sensitivities demonstrates that the method can accurately calculate the parameter sensitivity.

Beltracchi, Todd J.↗
NASA NTRS1989

Velocity-pressure integrated versus penalty finite element methods for high Reynolds number flows

Velocity-pressure integrated and consistent penalty finite element computations of high Reynolds number laminar flows are presented. In both methods the pressure has been interpolated using linear shape functions for a triangular element which is contained inside the biquadratic flow element. It has been shown previously that the pressure interpolation method, when used in conjunction with the velocity-pressure integrated method, yields accurate computational results for high-Reynolds-number flows. It is shown in this paper that use of the same pressure interpolation method in the consistent penalty finite element method yields computational results which are comparable to those of the velocity-pressure integrated method for both the velocity and the pressure fields. Accuracy of the two finite element methods has been demonstrated by comparing the computational results with available experimental data and/or fine grid finite difference computational results. Advantages and disadvantages of the two finite element methods are discussed on the basis of accuracy and convergence nature. Example problems considered include a lid-driven cavity flow of Reynolds number 10000, a laminar backward-facing step flow and a laminar flow through a nest of cylinders.

Kim, S.-W.↗
NASA NTRS1989

An advanced probabilistic structural analysis method for implicit performance functions

In probabilistic structural analysis, the performance or response functions usually are implicitly defined and must be solved by numerical analysis methods such as finite element methods. In such cases, the most commonly used probabilistic analysis tool is the mean-based, second-moment method which provides only the first two statistical moments. This paper presents a generalized advanced mean value (AMV) method which is capable of establishing the distributions to provide additional information for reliability design. The method requires slightly more computations than the second-moment method but is highly efficient relative to the other alternative methods. In particular, the examples show that the AMV method can be used to solve problems involving non-monotonic functions that result in truncated distributions.

Wu, Y.-T.↗
NASA NTRS1989

Experimental tests of methods for the measurement of rainfall rate using an airborne dual-wavelength radar

Several attenuation-based methods for estimating the rainfall rate were applied to measurements made by an airborne dual-wavelength radar operating at 0.87 cm, the K(a)-band, and at 3 cm, the X-band. These methods included the traditional Z-R methods, designated Z(X)-R and Z(K)-R for the X- and K(a) band wavelengths, respectively; single- and dual-wavelength surface reference techniques (SRT and DSRT, respectively); and standard dual-wavelength methods with and without range-averaging. As the primary sources of error for these methods are nearly independent, agreement among the rain rates obtained with these methods would lend confidence in the results. Correlation coefficients obtained between the rainfall rates with the Z(X)-R and DSRT methods were generally between 0.7 and 0.9. Good agreement among the methods occurred most often in stratiform rain for rain rates betwen a few mm/hr to about 15 mm/hr, i.e., where attenuation at the shorter wavelength is significant but not so severe as to result in a loss of signal.

Meneghini, R.↗
NASA NTRS1989

Performance of FFT methods in local gravity field modelling

Fast Fourier transform (FFT) methods provide a fast and efficient means of processing large amounts of gravity or geoid data in local gravity field modelling. The FFT methods, however, has a number of theoretical and practical limitations, especially the use of flat-earth approximation, and the requirements for gridded data. In spite of this the method often yields excellent results in practice when compared to other more rigorous (and computationally expensive) methods, such as least-squares collocation. The good performance of the FFT methods illustrate that the theoretical approximations are offset by the capability of taking into account more data in larger areas, especially important for geoid predictions. For best results good data gridding algorithms are essential. In practice truncated collocation approaches may be used. For large areas at high latitudes the gridding must be done using suitable map projections such as UTM, to avoid trivial errors caused by the meridian convergence. The FFT methods are compared to ground truth data in New Mexico (xi, eta from delta g), Scandinavia (N from delta g, the geoid fits to 15 cm over 2000 km), and areas of the Atlantic (delta g from satellite altimetry using Wiener filtering). In all cases the FFT methods yields results comparable or superior to other methods.

Forsberg, Rene↗
NASA NTRS1989

Fourier transform methods in local gravity modeling

New algorithms were derived for computing terrain corrections, all components of the attraction of the topography at the topographic surface and the gradients of these attractions. These algoriithms utilize fast Fourier transforms, but, in contrast to methods currently in use, all divergences of the integrals are removed during the analysis. Sequential methods employing a smooth intermediate reference surface were developed to avoid the very large transforms necessary when making computations at high resolution over a wide area. A new method for the numerical solution of Molodensky's problem was developed to mitigate the convergence difficulties that occur at short wavelengths with methods based on a Taylor series expansion. A trial field on a level surface is continued analytically to the topographic surface, and compared with that predicted from gravity observations. The difference is used to compute a correction to the trial field and the process iterated. Special techniques are employed to speed convergence and prevent oscillations. Three different spectral methods for fitting a point-mass set to a gravity field given on a regular grid at constant elevation are described. Two of the methods differ in the way that the spectrum of the point-mass set, which extends to infinite wave number, is matched to that of the gravity field which is band-limited. The third method is essentially a space-domain technique in which Fourier methods are used to solve a set of simultaneous equations.

Harrison, J. C.↗
NASA NTRS1990

Experimental implementation of the modified independent modal space control method

An experimental realization of a modified independent modal space control (MIMSC) method to control the vibration of a flexible cantilevered beam is presented. In its operation the method relies on the use of one piezoelectric actuator to control several vibration modes through a time-sharing strategy. The effectiveness of the MIMSC method in damping out the beam vibration is demonstrated by comparing the results with those obtained by other modal control methods. Two methods are considered, the independent modal space control method and the pseudo-inverse method. The feasibility of the MIMSC method as a viable alternative for controlling large flexible structures with a very small number of actuators is emphasized.

Baz, A.↗
NASA NTRS1990

A first-order Lyapunov robustness method for linear systems with uncertain parameters

A method for stability-robustness analysis based on a quadratic Liapunov function that varies linearly with uncertainty parameters is derived. Linear time-invariant systems with structured uncertainties are discussed. The Liapunov function is optimized numerically to maximize the robustness region in parameter space. Numerical results are given for four examples in which the first-order method is compared to previous Liapunov methods. While the zero-order method is slightly better than the first-order method for one example, the first-order method is clearly superior in the other three (more realistic) examples. The first-order method is especially superior for the active control of flexible structures, where robustness with respect to (1) unmodeled coupling between modeled modes and (2) unmodeled modes is important. For such applications, the first-order method is much better at detecting the increased robustness associated with increased separation between frequencies.

Leal, M. A.↗
NASA NTRS1991

A method for the geometrically nonlinear analysis of compressively loaded prismatic composite structures

A method was developed for the geometrically nonlinear analysis of the static response of thin-walled stiffened composite structures loaded in uniaxial or biaxial compression. The method is applicable to arbitrary prismatic configurations composed of linked plate strips, such as stiffened panels and thin-walled columns. The longitudinal ends of the structure are assumed to be simply supported, and geometric shape imperfections can be modeled. The method can predict the nonlinear phenomena of postbuckling strength and imperfection sensitivity which are exhibited by some buckling-dominated structures. The method is computer-based and is semi-analytic in nature, making it computationally economical in comparison to finite element methods. The method uses a perturbation approach based on the use of a series of buckling mode shapes to represent displacement contributions associated with nonlinear response. Displacement contributions which are of second order in the model amplitudes are incorported in addition to the buckling mode shapes. The principle of virtual work is applied using a finite basis of buckling modes, and terms through the third order in the model amplitudes are retained. A set of cubic nonlinear algebraic equations are obtained, from which approximate equilibrium solutions are determined. Buckling mode shapes for the general class of structure are obtained using the VIPASA analysis code within the PASCO stiffened-panel design code. Thus, subject to some additional restrictions in loading and plate anisotropy, structures which can be modeled with respect to buckling behavior by VIPASA can be analyzed with respect to nonlinear response using the new method. Results obtained using the method are compared with both experimental and analytical results in the literature. The configurations investigated include several different unstiffened and blade-stiffening panel configurations, featuring both homogeneous, isotropic materials, and laminated composite material.

Stoll, Frederick↗
NASA NTRS1989

Projection methods for the numerical solution of Markov chain models

Projection methods for computing stationary probability distributions for Markov chain models are presented. A general projection method is a method which seeks an approximation from a subspace of small dimension to the original problem. Thus, the original matrix problem of size N is approximated by one of dimension m, typically much smaller than N. A particularly successful class of methods based on this principle is that of Krylov subspace methods which utilize subspaces of the form span(v,av,...,A(exp m-1)v). These methods are effective in solving linear systems and eigenvalue problems (Lanczos, Arnoldi,...) as well as nonlinear equations. They can be combined with more traditional iterative methods such as successive overrelaxation, symmetric successive overrelaxation, or with incomplete factorization methods to enhance convergence.

Saad, Youcef↗
NASA NTRS1992

Improved method for solving the viscous shock layer equations

An improved method for solving the viscous shock layer equations for supersonic/hypersonic flows past blunt-nosed bodies is presented. The method is capable of handling slender to thick bodies. The solution is obtained by solving a coupled set of five equations, built of the four basic viscous shock layer equations and an additional equation for the standoff distance. The coupling of the equations prevents the local iterations divergence problems encountered by previous methods of solution far downstream on slender bodies. It also eliminates the need for local iterations, which were required by previous methods of solution, for a first-order scheme in the streamwise direction. A new global iteration procedure is employed to impose the shock boundary conditions. The procedure prevents the global iteration instability encountered by the basic method of solution and improves the convergence rate of the global iteration procedure of later methods devised to overcome this difficulty. The new technique reduces the computation time by 65-95 percent as compared to previous methods of solution. The method can efficiently be implemented in vector/parallel computers.

Gordon, Rachel↗