Search NASA⌕ Search

SEARCH · Search NASA

Results for “methods”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Data-parallel lower-upper relaxation method for reacting flows

The implicit lower-upper symmetric Gauss-Seidel (LU-SGS) method of Yoon and Jameson is modified for use on massively parallel computers. The method has been implemented on the Thinking Machines CM-5 and the MasPar MP-1 and MP-2, where large percentages of the theoretical peak floating point performance are obtained. It is shown that the new data-parallel LU relaxation method has better convergence properties than the original method for two different inviscid compressible flow simulations. The convergence is also improved for five-species reacting air computations. The performance of the method on various partitions of the CM-5 and on the MasPar computers is discussed. The new method shows promise for the efficient simulation of very large perfect gas and reacting flows.

Candler, Graham V.↗
NASA NTRS1995

Investigation of dispersion-relation-preserving scheme and spectral analysis methods for acoustic waves

Important characteristics of the aeroacoustic wave propagation are mostly encoded in their dispersion relations. Hence, a computational aeroacoustic (CAA) algorithm, which reasonably preserves these relations, was investigated. It was derived using an optimization procedure to ensure, that the numerical derivatives preserved the wave number and angular frequency of the differential terms in the linearized, 2-D Euler equations. Then, simulations were performed to validate the scheme and a compatible set of discretized boundary conditions. The computational results were found to agree favorably with the exact solutions. The boundary conditions were transparent to the outgoing waves, except when the disturbance source was close to a boundary. The time-domain data generated by such CAA solutions were often intractable until their spectra was analyzed. Therefore, the relative merits of three different methods were included in the study. For simple, periodic waves, the periodogram method produced better estimates of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window was more effective when used with the weighted-overlapped-segment-averaging and Blackman-Tukey methods gave better results than the periodogram method. Finally, it was demonstrated that the representation of time domain-data was significantly dependent on the particular spectral analysis method employed.

Vanel, Florence O.↗
NASA NTRS1995

A deterministic particle method for one-dimensional reaction-diffusion equations

We derive a deterministic particle method for the solution of nonlinear reaction-diffusion equations in one spatial dimension. This deterministic method is an analog of a Monte Carlo method for the solution of these problems that has been previously investigated by the author. The deterministic method leads to the consideration of a system of ordinary differential equations for the positions of suitably defined particles. We then consider the time explicit and implicit methods for this system of ordinary differential equations and we study a Picard and Newton iteration for the solution of the implicit system. Next we solve numerically this system and study the discretization error both analytically and numerically. Numerical computation shows that this deterministic method is automatically adaptive to large gradients in the solution.

Mascagni, Michael↗
NASA NTRS1970

A Study of the Electrotransport in Thin Films: An Advanced Method for the Study of Electrotransport in Thin Films Using Electrical Resistance Measurements

To date, very little is understood about electrotransport in thin films. One reason for this is the lack of methods in which electromigration can be adequately monitored dynamically. In this thesis, a new nondestructive method for the study of the electrotransport phenomenon in thin aluminum films is presented. This method makes use of electrical resistance measurements of various regions along the test sample to monitor changes in resistance resulting from electromigration. The advantage of this method is that resistance changes can be observed long before void formation can be seen using standard microscopy methods. In order to achieve the best efficiency using this method, special consideration must be given to the design of the sample. Design problems and their respective solutions are given. The results of the electrical resistance measurements are interpreted by means of photographs taken with an scanning electron microscope. The increase in resistance of the aluminum samples is interpreted to be due to void formation. The expected decrease in resistance near the anode was never observed. It is shown that mass accumulation takes on the form of hillocks and whiskers on the surface of the aluminum film, which does not contribute appreciably to conductivity.

Hummel, R. E.↗
NASA NTRS1996

Element Library for Three-Dimensional Stress Analysis by the Integrated Force Method

The Integrated Force Method, a recently developed method for analyzing structures, is extended in this paper to three-dimensional structural analysis. First, a general formulation is developed to generate the stress interpolation matrix in terms of complete polynomials of the required order. The formulation is based on definitions of the stress tensor components in term of stress functions. The stress functions are written as complete polynomials and substituted into expressions for stress components. Then elimination of the dependent coefficients leaves the stress components expressed as complete polynomials whose coefficients are defined as generalized independent forces. Such derived components of the stress tensor identically satisfy homogenous Navier equations of equilibrium. The resulting element matrices are invariant with respect to coordinate transformation and are free of spurious zero-energy modes. The formulation provides a rational way to calculate the exact number of independent forces necessary to arrive at an approximation of the required order for complete polynomials. The influence of reducing the number of independent forces on the accuracy of the response is also analyzed. The stress fields derived are used to develop a comprehensive finite element library for three-dimensional structural analysis by the Integrated Force Method. Both tetrahedral- and hexahedral-shaped elements capable of modeling arbitrary geometric configurations are developed. A number of examples with known analytical solutions are solved by using the developments presented herein. The results are in good agreement with the analytical solutions. The responses obtained with the Integrated Force Method are also compared with those generated by the standard displacement method. In most cases, the performance of the Integrated Force Method is better overall.

Kaljevic, Igor↗
NASA NTRS1995

A test of a vortex method for the computation of flap side edge noise

Upon approach to landing, a major source location of airframe noise occurs at the side edges of the part span, trailing edge flaps. In the vicinity of these flaps, a complex arrangement of spanwise flow with primary and secondary tip vortices may form. Each of these vortices is observed to become fully three-dimensional. In the present study, a numerical model is developed to investigate the noise radiated from the side edge of a flap. The inherent three-dimensionality of this flow forces us to carefully consider a numerical scheme which will be both accurate in its prediction of the flow acoustics and also computationally efficient. Vortex methods have offered a fast and efficient means of simulating many two and three-dimensional, vortex dominated flows. In vortex methods, the time development of the flow is tracked by following exclusively the vorticity containing regions. Through the Biot-Savart law, knowledge of the vorticity field enables one to obtain flow quantities at any desired location during the flow evolution. In the present study, a numerical procedure has been developed which incorporates the Lagrangian approach of vortex methods into a calculation for the noise radiated by a flow-surface interaction. In particular, the noise generated by a vortex in the presence of a flat half plane is considered. This problem serves as a basic model of flap edge flow. It also permits the direct comparison between our computed results and previous acoustic analyses performed for this problem. In our numerical simulations, the mean flow is represented by the complex potential W(z) = Aiz(exp l/2), which is obtained through conformal mapping techniques. The magnitude of the mean flow is controlled by the parameter A. This mean flow has been used in the acoustic analysis by Hardin and is considered a reasonable model of the flow field in the vicinity of the edge and away from the leading and trailing edges of the flap. To represent the primary vortex which occurs near the flap, a point vortex is introduced just below the flat half plane. Using a technique from panel methods, boundary conditions on the flap surface are satisfied by the introduction of a row of stationary point vortices along the extent of the flap. At each time step in the calculation, the strength of these vortices is chosen to eliminate the normal velocity at intermediary collocation points. The time development of the overall flow field is then tracked using standard techniques from vortex methods. Vortex trajectories obtained through this computation are in good agreement with those predicted by the analytical solution given by Hardin, thus verifying the viability of this procedure for more complex flow arrangements. For the flow acoustics, the Ffowcs Williams-Hawkings equation is numerically integrated. This equation supplies the far field acoustic pressure based upon pressures occurring along the flap surface. With our vortex method solution, surface pressures may be obtained with exceptional resolution. The Ffowcs Williams-Hawkings equation is integrated using a spatially fourth order accurate Simpson's rule. Rational function interpolation is used to obtain the surface pressures at the appropriate retarded times. Comparisons between our numerical results for the acoustic pressure and those predicted by the Hardin analysis have been made. Preliminary results indicate the need for an improved integration technique. In the future, the numerical procedure developed in this study will be applied to the case of a rectangular flap of finite thickness and ultimately modified for application to the fully three-dimensional problem.

Martin, James E.↗
NASA NTRS1995

Standard Methods for Bolt-Bearing Testing of Textile Composites

The response of three 2-D braided materials to bolt bearing loading was evaluated using data generated by Boeing Defense and Space Group in Philadelphia, PA. Three test methods, stabilized single shear, unstabilized single shear, and double shear, were compared. In general, these textile composites were found to be sensitive to bolt bearing test methods. The stabilized single shear method yielded higher strengths than the unstabilized single shear method in all cases. The double shear test method always produced the highest strengths but these results may be somewhat misleading. It is therefore recommended that standard material comparisons be made using the stabilized single shear test method. The effects of two geometric parameters, W/D and e/D, were also studied. An evaluation of the effect of the specimen width (W) to hole diameter (D) ratio concluded that bolt bearing responses were consistent with open hole tension results. A W/D ratio of 6 or greater should be maintained. The proximity of the hole to the specimen edge significantly affected strength. In all cases, strength was improved by increasing the ratio of the distance from the hole center to the specimen edge (e) to the hole diameter (D) above 2. An e/D ratio of 3 or greater is recommended.

Portanova, M. A.↗
NASA NTRS1996

Artificial Boundary Conditions Based on the Difference Potentials Method

While numerically solving a problem initially formulated on an unbounded domain, one typically truncates this domain, which necessitates setting the artificial boundary conditions (ABC's) at the newly formed external boundary. The issue of setting the ABC's appears to be most significant in many areas of scientific computing, for example, in problems originating from acoustics, electrodynamics, solid mechanics, and fluid dynamics. In particular, in computational fluid dynamics (where external problems present a wide class of practically important formulations) the proper treatment of external boundaries may have a profound impact on the overall quality and performance of numerical algorithms. Most of the currently used techniques for setting the ABC's can basically be classified into two groups. The methods from the first group (global ABC's) usually provide high accuracy and robustness of the numerical procedure but often appear to be fairly cumbersome and (computationally) expensive. The methods from the second group (local ABC's) are, as a rule, algorithmically simple, numerically cheap, and geometrically universal; however, they usually lack accuracy of computations. In this paper we first present a survey and provide a comparative assessment of different existing methods for constructing the ABC's. Then, we describe a relatively new ABC's technique of ours and review the corresponding results. This new technique, in our opinion, is currently one of the most promising in the field. It enables one to construct such ABC's that combine the advantages relevant to the two aforementioned classes of existing methods. Our approach is based on application of the difference potentials method attributable to V. S. Ryaben'kii. This approach allows us to obtain highly accurate ABC's in the form of certain (nonlocal) boundary operator equations. The operators involved are analogous to the pseudodifferential boundary projections first introduced by A. P. Calderon and then also studied by R. T. Seeley. The apparatus of the boundary pseudodifferential equations, which has formerly been used mostly in the qualitative theory of integral equations and PDE'S, is now effectively employed for developing numerical methods in the different fields of scientific computing.

Tsynkov, Semyon V.↗
NASA NTRS1995

Development of a Probabilistic Component Mode Synthesis Method for the Analysis of Non-Deterministic Substructures

Standard methods of structural dynamic analysis assume that the structural characteristics are deterministic. Recognizing that these characteristics are actually statistical in nature, researchers have recently developed a variety of methods that use this information to determine probabilities of a desired response characteristic, such as natural frequency, without using expensive Monte Carlo simulations. One of the problems in these methods is correctly identifying the statistical properties of primitive variables such as geometry, stiffness, and mass. This paper presents a method where the measured dynamic properties of substructures are used instead as the random variables. The residual flexibility method of component mode synthesis is combined with the probabilistic methods to determine the cumulative distribution function of the system eigenvalues. A simple cantilever beam test problem is presented that illustrates the theory.

Brown, Andrew M.↗
NASA NTRS1996

The Multigrid-Mask Numerical Method for Solution of Incompressible Navier-Stokes Equations

A multigrid-mask method for solution of incompressible Navier-Stokes equations in primitive variable form has been developed. The main objective is to apply this method in conjunction with the pseudospectral element method solving flow past multiple objects. There are two key steps involved in calculating flow past multiple objects. The first step utilizes only Cartesian grid points. This homogeneous or mask method step permits flow into the interior rectangular elements contained in objects, but with the restriction that the velocity for those Cartesian elements within and on the surface of an object should be small or zero. This step easily produces an approximate flow field on Cartesian grid points covering the entire flow field. The second or heterogeneous step corrects the approximate flow field to account for the actual shape of the objects by solving the flow field based on the local coordinates surrounding each object and adapted to it. The noise occurring in data communication between the global (low frequency) coordinates and the local (high frequency) coordinates is eliminated by the multigrid method when the Schwarz Alternating Procedure (SAP) is implemented. Two dimensional flow past circular and elliptic cylinders will be presented to demonstrate the versatility of the proposed method. An interesting phenomenon is found that when the second elliptic cylinder is placed in the wake of the first elliptic cylinder a traction force results in a negative drag coefficient.

Ku, Hwar-Ching↗
NASA NTRS1996

Krylov Subspace and Multigrid Methods Applied to the Incompressible Navier-Stokes Equations

We consider numerical solution methods for the incompressible Navier-Stokes equations discretized by a finite volume method on staggered grids in general coordinates. We use Krylov subspace and multigrid methods as well as their combinations. Numerical experiments are carried out on a scalar and a vector computer. Robustness and efficiency of these methods are studied. It appears that good methods result from suitable combinations of GCR and multigrid methods.

Vuik, C.↗
NASA NTRS1997

Estimation of Slow Crack Growth Parameters for Constant Stress-Rate Test Data of Advanced Ceramics and Glass by the Individual Data and Arithmetic Mean Methods

The two estimation methods, individual data and arithmetic mean methods, were used to determine the slow crack growth (SCG) parameters (n and D) of advanced ceramics and glass from a large number of room- and elevated-temperature constant stress-rate ('dynamic fatigue') test data. For ceramic materials with Weibull modulus greater than 10, the difference in the SCG parameters between the two estimation methods was negligible; whereas, for glass specimens exhibiting Weibull modulus of about 3, the difference was amplified, resulting in a maximum difference of 16 and 13 %, respectively, in n and D. Of the two SCG parameters, the parameter n was more sensitive to the estimation method than the other. The coefficient of variation in n was found to be somewhat greater in the individual data method than in the arithmetic mean method.

Choi, Sung R.↗
NASA NTRS1996

A General Interface Method for Aeroelastic Analysis of Aircraft

The aeroelastic analysis of an aircraft requires an accurate and efficient procedure to couple aerodynamics and structures. The procedure needs an interface method to bridge the gap between the aerodynamic and structural models in order to transform loads and displacements. Such an interface method is described in this report. This interface method transforms loads computed by any aerodynamic code to a structural finite element (FE) model and converts the displacements from the FE model to the aerodynamic model. The approach is based on FE technology in which virtual work is employed to transform the aerodynamic pressures into FE nodal forces. The displacements at the FE nodes are then converted back to aerodynamic grid points on the aircraft surface through the reciprocal theorem in structural engineering. The method allows both high and crude fidelities of both models and does not require an intermediate modeling. In addition, the method performs the conversion of loads and displacements directly between individual aerodynamic grid point and its corresponding structural finite element and, hence, is very efficient for large aircraft models. This report also describes the application of this aero-structure interface method to a simple wing and an MD-90 wing. The results show that the aeroelastic effect is very important. For the simple wing, both linear and nonlinear approaches are used. In the linear approach, the deformation of the structural model is considered small, and the loads from the deformed aerodynamic model are applied to the original geometry of the structure. In the nonlinear approach, the geometry of the structure and its stiffness matrix are updated in every iteration and the increments of loads from the previous iteration are applied to the new structural geometry in order to compute the displacement increments. Additional studies to apply the aero-structure interaction procedure to more complicated geometry will be conducted in the second phase of the present contract.

Tzong, T.↗
NASA NTRS1998

Interfacing with Legacy using Remote Method Invocation

The assignment described was enough to make a neophyte Java developer bolt for the door: provide a remote method for use by an applet which invokes a native method that wraps a function in an existing legacy library. The purpose of the remote method is to return an instance of a class object whose contents reflect the data structure returned by the legacy function. While embroiled in implementation, I would have spent the time wading through their JNI use group archive as well, but I couldn't seem to locate one. Subsequently, I made the decision to try to document my findings in order to assist others. Before we start on the class design, let's look at what the existing legacy code does. The C function to be called, Get-Legacy-Data, consists of two steps: an ASII file is read from the local disk and its contents are parsed into a Legacy_Type structure whose address is passed as an argument by the caller. The legacy code was compiled into a shared object library, legacy. so, using the IRIX 6.2 compiler and then loaded onto the Web server, a Silicon Graphics Indy station loaded with the IRIX 6.4 operating system. As far as the class design is concerned, the first thing required is a class to act as a template for the data structure returned by the legacy function. This class, JLegacy, declares a series of public instance variables which correspond to the members of Legacy_Type and provides a parameterless constructor. This constructor is never called, not even by the native method which allocates the object for return to the remote method. Next, the remote interface declaration for the remote object must be defined. In order for JLegacyRO to implement getJLegacy, JLegacyRO must interface with the existing legacy code through a native method, getn. getn is declared in the JLegacyRO class but implemented in C, just like the legacy code. getn returns a JLegacy instance and is declared static since its implementation is the same for all instances of the JLegacyRO class.

Howard, Scott M.↗
NASA NTRS1998

Application of the Probabilistic Dynamic Synthesis Method to the Analysis of a Realistic Structure

The Probabilistic Dynamic Synthesis method is a new technique for obtaining the statistics of a desired response engineering quantity for a structure with non-deterministic parameters. The method uses measured data from modal testing of the structure as the input random variables, rather than more "primitive" quantities like geometry or material variation. This modal information is much more comprehensive and easily measured than the "primitive" information. The probabilistic analysis is carried out using either response surface reliability methods or Monte Carlo simulation. A previous work verified the feasibility of the PDS method on a simple seven degree-of-freedom spring-mass system. In this paper, extensive issues involved with applying the method to a realistic three-substructure system are examined, and free and forced response analyses are performed. The results from using the method are promising, especially when the lack of alternatives for obtaining quantitative output for probabilistic structures is considered.

Brown, Andrew M.↗
NASA NTRS1997

A Modified Linear-Mixing Method for Calculating Atmospheric Path Radiances of Aerosol Mixtures

The top-of-atmosphere (TOA) path radiance generated by an aerosol mixture can be synthesized by linearly adding the contributions of the individual aerosol components, weighted by their fractional optical depths. The method, known as linear mixing, is exact in the single-scattering limit. When multiple scattering is significant, the method reproduces the atmospheric path radiance of the mixture with less than 3% errors for weakly absorbing aerosols up to optical thickness of 0.5. However, when strongly absorbing aerosols are included in the mixture, the errors are much larger. This is due to neglecting the effect of multiple interactions between the aerosol components, especially when the values of the single-scattering albedos of these components are so different that the parameter epsilon = (Sigma)f(sub i) absolute value of bar omega(sub i) - bar omega(sub mix)/bar omega(sub i), is larger than approx. 0.1, where bar omega(sub i) and f(sub i) are the single-scattering albedo and the fractional abundance of the i th component, and bar omega(sub i) is the effective single-scattering albedo of the mixture. We describe an empirical, modified linear-mixing method which effectively accounts for the multiple interactions between aerosol components. The modified and standard methods are identical when epsilon = 0.0 and give similar results when epsilon less than or equal to 0.05. For optical depths larger than approx. 0.5, or when epsilon greater than 0.05, only the modified method can reproduce the radiances within 5% error for common aerosol types up to optical thickness of 2.0. Because this method facilitates efficient and accurate atmospheric path radiance calculations for mixtures of a wide variety of aerosol types, it will be used as part of the aerosol retrieval methodology for the Earth Observing System (EOS) multiangle imaging spectroradiometer (MISR), scheduled for launch into polar orbit in 1998.

Abdou, W. A.↗
NASA NTRS1997

A Modified Linear-Mixing Method for Calculating Atmospheric Path Radiances of Aerosol Mixtures

The top-of-atmosphere (TOA) path radiance generated by an aerosol mixture can be synthesized by linearly adding the contributions of the individual aerosol components, weighted by their fractional optical depths. The method, known as linear mixing, is exact in the single-scattering limit. When multiple scattering is significant, the method reproduces the atmospheric path radiance of the mixture with less than 3% errors for weakly absorbing aerosols up to optical thickness of 0.5. However, when strongly absorbing aerosols are included in the mixture, the errors are much larger. This is due to neglecting the effect of multiple interactions between the aerosol components, especially when the values of the single-scattering albedos of these components are so different that the parameter e = the sum of f(sub i)[(bar)omega(sub i) - (bar)omega(sub mix)]/(bar)omega(sub i) is larger than approximately 0.1, where (bar)omega(sub i)and f(sub i) are the single-scattering albedo and the fractional abundance of the ith component, and (bar)omega(sub mix) is the effective single-scattering albedo of the Mixture. We describe an empirical, modified linear-mixing method which effectively accounts for the multiple interactions between aerosol components. The modified and standard methods are identical when epsilon = 0.0 and give similar results when epsilon is less than or equal to 0.05. For optical depths larger than approximately 0.5, or when epsilon is greater than 0.05, only the modified method can reproduce the radiances within 5% error for common aerosol types up to optical thickness of 2.0. Because this method facilitates efficient and accurate atmospheric path radiance calculations for mixtures of a wide variety of aerosol types, it will be used as part of the aerosol retrieval methodology for the Earth Observing System (EOS) multiangle imaging spectroradiometer (MISR), scheduled for launch into polar orbit in 1998.

Abdou, W. A.↗
NASA NTRS1998

Application of the Probabilistic Dynamic Synthesis Method to Realistic Structures

The Probabilistic Dynamic Synthesis method is a technique for obtaining the statistics of a desired response engineering quantity for a structure with non-deterministic parameters. The method uses measured data from modal testing of the structure as the input random variables, rather than more "primitive" quantities like geometry or material variation. This modal information is much more comprehensive and easily measured than the "primitive" information. The probabilistic analysis is carried out using either response surface reliability methods or Monte Carlo simulation. In previous work, the feasibility of the PDS method applied to a simple seven degree-of-freedom spring-mass system was verified. In this paper, extensive issues involved with applying the method to a realistic three-substructure system are examined, and free and forced response analyses are performed. The results from using the method are promising, especially when the lack of alternatives for obtaining quantitative output for probabilistic structures is considered.

Brown, Andrew M.↗