Engineering PapersSearch

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 37 records · Page 2

Comparison of linear inversion methods by examination of the duality between iterative and inverse matrix methods

Linear numerical inversion methods applied to atmospheric remote sounding generally can be categorized in two ways: (1) iterative, and (2) inverse matrix methods. However, these two categories are not unrelated; a duality exists between them. In other words, given an iterative scheme, a corresponding inverse matrix method exists, and conversely. This duality concept is developed for the more familiar linear methods. The iterative duals are compared with the classical linear iterative approaches and their differences analyzed. The importance of the initial profile in all methods is stressed. Calculations using simulated data are made to compare accuracies and to examine the dependence of the solution on the initial profile.

Fleming, H. E.

A finite element method for transmission in non-uniform ducts without flow Comparison with the method of weighted residuals

A finite element method is derived for studying the transmission of sound in nonuniform ducts without flow. The method is based on a weighted residuals technique and eight noded isoparametric elements. The method employs two computational schemes: the Helmholtz equation found through conservation equations, and the conservation equations themselves. Solutions are found for transmission and reflection coefficients for both two-dimensional and circular geometries. Numerical solutions are compared with those for a method of weighted residuals in the form of a modified Galerkin method, showing good agreement.

Astley, R. J.

A comparison between progressive extension method (PEM) and iterative method (IM) for magnetic field extrapolations in the solar atmosphere

This paper presents a comparison between two numerical methods for the extrapolation of nonlinear force-free magnetic fields, viz the Iterative Method (IM) and the Progressive Extension Method (PEM). The advantages and disadvantages of these two methods are summarized, and the accuracy and numerical instability are discussed. On the basis of this investigation, it is claimed that the two methods do resemble each other qualitatively.

Wu, S. T.

Gel stretch method: a new method to measure constitutive properties of cardiac muscle cells

Diastolic dysfunction is an important cause of congestive heart failure; however, the basic mechanisms causing diastolic congestive heart failure are not fully understood, especially the role of the cardiac muscle cell, or cardiocyte, in this process. Before the role of the cardiocyte in this pathophysiology can be defined, methods for measuring cardiocyte constitutive properties must be developed and validated. Thus this study was designed to evaluate a new method to characterize cardiocyte constitutive properties, the gel stretch method. Cardiocytes were isolated enzymatically from normal feline hearts and embedded in a 2% agarose gel containing HEPES-Krebs buffer and laminin. This gel was cast in a shape that allowed it to be placed in a stretching device. The ends of the gel were held between a movable roller and fixed plates that acted as mandibles. Distance between the right and left mandibles was increased using a stepper motor system. The force applied to the gel was measured by a force transducer. The resultant cardiocyte strain was determined by imaging the cells with a microscope, capturing the images with a CCD camera, and measuring cardiocyte and sarcomere length changes. Cardiocyte stress was characterized with a finite-element method. These measurements of cardiocyte stress and strain were used to determine cardiocyte stiffness. Two variables affecting cardiocyte stiffness were measured, the passive elastic spring and viscous damping. The passive spring was assessed by increasing the force on the gel at 1 g/min, modeling the resultant stress vs. strain relationship as an exponential [sigma = A/k(ekepsilon - 1)]. In normal cardiocytes, A = 23.0 kN/m2 and k = 16. Viscous damping was assessed by examining the loop area between the stress vs. strain relationship during 1 g/min increases and decreases in force. Normal cardiocytes had a finite loop area = 1.39 kN/m2, indicating the presence of viscous damping. Thus the gel stretch method provided accurate measurements of cardiocyte constitutive properties. These measurements have allowed the first quantitative assessment of passive elastic spring properties and viscous damping in normal mammalian cardiocytes.

Non-NASA Center

A Simplified FANG Cell-Centered Finite-Volume Method and Comparison with Other Methods for Trouble-Prone Grids

We propose a simplication of the face-averaged nodal-gradient (FANG) method for a cell-centered finite-volume Euler/Navier-Stokes solver on arbitrary grids, and compare it with other gradient methods for trouble-prone grids in two dimensions. The implementation of the FANG method is simplified by adding the face-neighbor cells of the cells around a node to a least-squares gradient stencil. The resulting method is stable for both triangular and quadrilateral grids. Although it increases the residual stencil for triangular grids, it allows the solver to work seamlessly for mixed grids and greatly simplifies the implementation, especially in three dimensions. For comparison, only explicit weighted/unweighted least-squares cell-centered and nodal gradient methods are considered. These gradients are used in both inviscid and viscous schemes, and we investigate their impact on the iterative convergence of an implicit defect-correction solver on difficult grids such as highly-curved-and-thin grids and highly distorted anisotropic grids. Finally, we will also consider a face-stencil-based limiter and compare it with a conventional cell-stencil-based limiter.

Hiroaki Nishikawa

Spectral methods for the Euler equations: Fourier methods and shock-capturing

Spectral methods for compressible flows are introduced in relation to finite difference and finite element techniques within the framework of the method of weighted residuals. Current spectral collocation methods are put in historical context. The basic concepts of Fourier spectral collocation methods are provided. Filtering strategies for shock-capturing approaches are also presented. Fourier shock capturing techniques are evaluated using a one dimensional, periodic astrophysical ""nozzle'' problem.

Hussaini, M. Y.

A method for computing the leading-edge suction in a higher-order panel method

Experimental data show that the phenomenon of a separation induced leading edge vortex is influenced by the wing thickness and the shape of the leading edge. Both thickness and leading edge shape (rounded rather than point) delay the formation of a vortex. Existing computer programs used to predict the effect of a leading edge vortex do not include a procedure for determining whether or not a vortex actually exists. Studies under NASA Contract NAS1-15678 have shown that the vortex development can be predicted by using the relationship between the leading edge suction coefficient and the parabolic nose drag. The linear theory FLEXSTAB was used to calculate the leading edge suction coefficient. This report describes the development of a method for calculating leading edge suction using the capabilities of the higher order panel methods (exact boundary conditions). For a two dimensional case, numerical methods were developed using the double strength and downwash distribution along the chord. A Gaussian quadrature formula that directly incorporates the logarithmic singularity in the downwash distribution, at all panel edges, was found to be the best method.

Ehlers, F. E.

Spectral methods for the Euler equations. I - Fourier methods and shock capturing

Spectral methods for compressible flows are introduced in relation to finite difference and finite element techniques within the framework of the method of weighted residuals. Current spectral collociation methods are put into historical context. The basic concepts of Fourier spectral collocation methods are provided. Filtering strategies for shock-capturing approaches are also presented. Fourier shock-capturing techniques are evaluated using a one-dimensional, periodic astrophysical 'nozzle' problem.

Hussaini, M. Y.

Optimization methods, flux conserving methods for steady state Navier-Stokes equation

Navier-Stokes equation as discretized by new flux conserving method proposed by Chang and Scott results in the system: vector F(vector x) = 0, where F is a vector valued function. The Optimization method we use is based on Quasi-Newton methods: given a nonlinear function vector F(vector x) = 0, we solve, Delta(vector x) = -BF(vector x), where Delta(vector x) is the correction term and B is the inverse Jacobian of F(x). Then, iteratively, vector(x(sub (i+1))) = vector(x (sub i)) + alpha.Delta(vector x(sub i)), where alpha is a line search correction term determined by a line search routine. We use the BFCG's update the Jacobian matrix B(sub k) at each iteration. It is well known that B(sub k) approaches B(*) at the solution X(*). This algorithm has several advantages over the Newton-Raphson method. For example, we do not need to calculate the Jacobian matrix at each iteration which is computationally very expensive.

Adeyeye, John

Analysis of corner cracks at hole by a 3-D weight function method with stresses from finite element method

Stress intensity factors for quarter-elliptical corner cracks emanating from a circular hole are determined using a 3-D weight function method combined with a 3-D finite element method. The 3-D finite element method is used to analyze uncracked configuration and provide stress distribution in the region where crack is to occur. Using this stress distribution as input, the 3-D weight function method is used to determine stress intensity factors. Three different loading conditions, i.e. remote tension, remote bending and wedge loading, are considered for a wide range in geometrical parameters. The significance in using 3-D uncracked stress distribution and the difference between single and double corner cracks are studied. Typical crack opening displacements are also provided. Comparisons are made with solutions available in the literature.

Zhao, W.

A New Method for Nonlinear and Nonstationary Time Series Analysis: The Empirical Mode Decomposition Method

A new method for analyzing nonlinear and nonstationary data has been developed. The key part of the method is the Empirical Mode Decomposition method with which any complicated data set can be decomposed into a finite and often small number of Intrinsic Mode Functions (IMF). An IMF is defined as any function having the same numbers of zero-crossing and extrema, and also having symmetric envelopes defined by the local maxima and minima respectively. The IMF also admits well-behaved Hilbert transform. This decomposition method is adaptive, and, therefore, highly efficient. Since the decomposition is based on the local characteristic time scale of the data, it is applicable to nonlinear and nonstationary processes. With the Hilbert transform, the Intrinsic Mode Functions yield instantaneous frequencies as functions of time that give sharp identifications of imbedded structures. The final presentation of the results is an energy-frequency-time distribution, designated as the Hilbert Spectrum. Classical nonlinear system models are used to illustrate the roles played by the nonlinear and nonstationary effects in the energy-frequency-time distribution.

Huang, Norden E.

Pre-coding method and apparatus for multiple source or time-shifted single source data and corresponding inverse post-decoding method and apparatus

A pre-coding method and device for improving data compression performance by removing correlation between a first original data set and a second original data set, each having M members, respectively. The pre-coding method produces a compression-efficiency-enhancing double-difference data set. The method and device produce a double-difference data set, i.e., an adjacent-delta calculation performed on a cross-delta data set or a cross-delta calculation performed on two adjacent-delta data sets, from either one of (1) two adjacent spectral bands coming from two discrete sources, respectively, or (2) two time-shifted data sets coming from a single source. The resulting double-difference data set is then coded using either a distortionless data encoding scheme (entropy encoding) or a lossy data compression scheme. Also, a post-decoding method and device for recovering a second original data set having been represented by such a double-difference data set.

Yeh, Pen-Shu

Pre-coding method and apparatus for multiple source or time-shifted single source data and corresponding inverse post-decoding method and apparatus

A pre-coding method and device for improving data compression performance by removing correlation between a first original data set and a second original data set, each having M members, respectively. The pre-coding method produces a compression-efficiency-enhancing double-difference data set. The method and device produce a double-difference data set, i.e., an adjacent-delta calculation performed on a cross-delta data set or a cross-delta calculation performed on two adjacent-delta data sets, from either one of (1) two adjacent spectral bands coming from two discrete sources, respectively, or (2) two time-shifted data sets coming from a single source. The resulting double-difference data set is then coded using either a distortionless data encoding scheme (entropy encoding) or a lossy data compression scheme. Also, a post-decoding method and device for recovering a second original data set having been represented by such a double-difference data set.

Yeh, Pen-Shu

Robust and Accurate Shock Capturing Method for High-Order Discontinuous Galerkin Methods

A simple yet robust and accurate approach for capturing shock waves using a high-order discontinuous Galerkin (DG) method is presented. The method uses the physical viscous terms of the Navier-Stokes equations as suggested by others; however, the proposed formulation of the numerical viscosity is continuous and compact by construction, and does not require the solution of an auxiliary diffusion equation. This work also presents two analyses that guided the formulation of the numerical viscosity and certain aspects of the DG implementation. A local eigenvalue analysis of the DG discretization applied to a shock containing element is used to evaluate the robustness of several Riemann flux functions, and to evaluate algorithm choices that exist within the underlying DG discretization. A second analysis examines exact solutions to the DG discretization in a shock containing element, and identifies a "model" instability that will inevitably arise when solving the Euler equations using the DG method. This analysis identifies the minimum viscosity required for stability. The shock capturing method is demonstrated for high-speed flow over an inviscid cylinder and for an unsteady disturbance in a hypersonic boundary layer. Numerical tests are presented that evaluate several aspects of the shock detection terms. The sensitivity of the results to model parameters is examined with grid and order refinement studies.

Atkins, Harold L.

A Numerical Combination of Extended Boundary Condition Method and Invariant Imbedding Method Applied to Light Scattering by Large Spheroids and Cylinders

The extended boundary condition method (EBCM) and invariant imbedding method (IIM) are two fundamentally different T-matrix methods for the solution of light scattering by nonspherical particles. The standard EBCM is very efficient but encounters a loss of precision when the particle size is large, the maximum size being sensitive to the particle aspect ratio. The IIM can be applied to particles in a relatively large size parameter range but requires extensive computational time due to the number of spherical layers in the particle volume discretization. A numerical combination of the EBCM and the IIM (hereafter, the EBCM+IIM) is proposed to overcome the aforementioned disadvantages of each method. Even though the EBCM can fail to obtain the T-matrix of a considered particle, it is valuable for decreasing the computational domain (i.e., the number of spherical layers) of the IIM by providing the initial T-matrix associated with an iterative procedure in the IIM. The EBCM+IIM is demonstrated to be more efficient than the IIM in obtaining the optical properties of large size parameter particles beyond the convergence limit of the EBCM. The numerical performance of the EBCM+IIM is illustrated through representative calculations in spheroidal and cylindrical particle cases.

invariant imbeddings

A method to estimate weight and dimensions of aircraft gas turbine engines. Volume 1: Method of analysis

Weight and envelope dimensions of aircraft gas turbine engines are estimated within plus or minus 5% to 10% using a computer method based on correlations of component weight and design features of 29 data base engines. Rotating components are estimated by a preliminary design procedure where blade geometry, operating conditions, material properties, shaft speed, hub-tip ratio, etc., are the primary independent variables used. The development and justification of the method selected, the various methods of analysis, the use of the program, and a description of the input/output data are discussed.

Pera, R. J.

Comparison of results obtained by solving the radiative transfer equation with an iterative method and a spherical harmonics method

Fluxes and intensities of light scattered by a model atmosphere are computed by a spherical harmonics approximation and by an iterative method of solving the radiative transfer equation and are compared. The large differences in the net fluxes and intensities reported by Dave and Armstrong (1974) for the two methods are reduced here by making a few changes in the iterative routine. Decreasing the polar angle increment from 2 to 1 deg in the iterative method of computing the source function does not improve the results as suggested by Dave and Armstrong.

Bahethi, O. P.

An artificial compression method for ENO schemes: The slope modification method

A simple and effective method of artificial compression is introduced. This method is based on a modification of the slopes of the ENO (essentially nonoscillatory) reconstruction and, with the help of suitable chosen parameters, greatly improves the resolution of the contact discontinuities. Numerical examples are provided to test the performance of the method and to give some suggestions as to the choice of the parameters.

Yang, Huanan