Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 NTRS1987

An evaluation of higher-order modal methods for calculating transient structural response

The present study evaluates a higher-order modal method proposed by Leung for transient structural analysis entitled the force-derivative method. This method repeatedly integrates by parts with respect to time the convolution-integral form of the structural response to produce successively better approximations to the contribution of the higher modes which are neglected in the modal summation. Comparisons are made of the force-derivative, the mode-displacement, and the mode-acceleration methods for several numerical example problems for various times, levels of damping, and forcing functions. The example problems include a tip-loaded cantilevered beam and a simply-supported multispan beam. The force-derivative method is shown to converge to an accurate solution in fewer modes than either the mode-displacement or the mode-acceleration methods. In addition, for problems in which there are a large number of closely-spaced frequencies whose mode shapes have a negligible contribution to the response, the force derivative method is very effective in representing the effect of the important, but otherwise neglected, higher modes.

Camarda, Charles J.↗
NASA NTRS1988

Application of a nonisentropic full potential method to AGARD standard airfoils

An entropy-correction method for the unsteady full potential equation is presented. The unsteady potential equation is modified to model the entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. A flux-biasing differencing method, which generates the proper amounts of artificial viscosity in supersonic regions, is used to discretize the flow equations in space. Calculated results are presented for the NLR 7301, NACA 0012, and NACA 64A010A airfoils. Comparisons of the present method and solutions of the Euler equations are presented for the NLR 7301 airfoil, and comparisons of the present method and experimental data are presented for all three airfoils. The comparisons show that the present method more accurately models solutions of the Euler equations and experiment than does the isentropic potential formulation. In addition, it is shown that modeling shock-generated entropy extends the range of validity of the full potential method.

Whitlow, Woodrow, Jr.↗
NASA NTRS1988

Accuracy analysis of the semi-analytical method for shape sensitivity calculation

The semianalytical method, widely used for calculating derivatives of static response with respect to design variables for structures modeled by finite elements, is studied in this paper. The paper shows that the method can have serious accuracy problems for shape design variables in structures modeled by beam, plate, truss, frame, and solid elements. The errors are shown to be associated with the structural model. An error index is developed to test the accuracy of the semianalytical method. It characterizes the difference in errors between a general finite-difference method and the semianalytical method. Moreover, a method for improving the accuracy of the semianalytical method (when possible) is provided. Examples are presented to demonstrate the use of the error index.

Barthelemy, Bruno↗
NASA NTRS1988

A comparison between onestep and other multiconic trajectory propagation methods

The paper presents an optimized version of the one-step multiconic method of spacecraft trajectory propagation in the presence of two attracting bodies. This method is compared to eleven other methods ranging from simple conic propagation to full integration of the equations of motion by applying the methods to a wide variety of test cases taken from planned space missions. This optimize one-step method is shown to be an order of magnitude better than other single-step methods for interplanetary and most lunar trajectories. Multistep multiconics methods and integration of the state offer another order of magnitude improvement, but at a proportionally higher cost and calculation time.

Kiedron, Krystyna↗
NASA NTRS1988

An overview of very high level software design methods

Very High Level design methods emphasize automatic transfer of requirements to formal design specifications, and/or may concentrate on automatic transformation of formal design specifications that include some semantic information of the system into machine executable form. Very high level design methods range from general domain independent methods to approaches implementable for specific applications or domains. Applying AI techniques, abstract programming methods, domain heuristics, software engineering tools, library-based programming and other methods different approaches for higher level software design are being developed. Though one finds that a given approach does not always fall exactly in any specific class, this paper provides a classification for very high level design methods including examples for each class. These methods are analyzed and compared based on their basic approaches, strengths and feasibility for future expansion toward automatic development of software systems.

Asdjodi, Maryam↗
NASA NTRS1989

Computational Methods for Structural Mechanics and Dynamics, part 1

The structural analysis methods research has several goals. One goal is to develop analysis methods that are general. This goal of generality leads naturally to finite-element methods, but the research will also include other structural analysis methods. Another goal is that the methods be amenable to error analysis; that is, given a physical problem and a mathematical model of that problem, an analyst would like to know the probable error in predicting a given response quantity. The ultimate objective is to specify the error tolerances and to use automated logic to adjust the mathematical model or solution strategy to obtain that accuracy. A third goal is to develop structural analysis methods that can exploit parallel processing computers. The structural analysis methods research will focus initially on three types of problems: local/global nonlinear stress analysis, nonlinear transient dynamics, and tire modeling.

Stroud, W. Jefferson↗
NASA NTRS1990

Differential equation based method for accurate approximations in optimization

A method to efficiently and accurately approximate the effect of design changes on structural response is described. The key to this method is to interpret sensitivity equations as differential equations that may be solved explicitly for closed form approximations, hence, the method is denoted the Differential Equation Based (DEB) method. Approximations were developed for vibration frequencies, mode shapes and static displacements. The DEB approximation method was applied to a cantilever beam and results compared with the commonly-used linear Taylor series approximations and exact solutions. The test calculations involved perturbing the height, width, cross-sectional area, tip mass, and bending inertia of the beam. The DEB method proved to be very accurate, and in most cases, was more accurate than the linear Taylor series approximation. The method is applicable to simultaneous perturbation of several design variables. Also, the approximations may be used to calculate other system response quantities. For example, the approximations for displacements are used to approximate bending stresses.

Pritchard, Jocelyn I.↗
NASA NTRS1990

Differential equation based method for accurate approximations in optimization

This paper describes a method to efficiently and accurately approximate the effect of design changes on structural response. The key to this new method is to interpret sensitivity equations as differential equations that may be solved explicitly for closed form approximations, hence, the method is denoted the Differential Equation Based (DEB) method. Approximations were developed for vibration frequencies, mode shapes and static displacements. The DEB approximation method was applied to a cantilever beam and results compared with the commonly-used linear Taylor series approximations and exact solutions. The test calculations involved perturbing the height, width, cross-sectional area, tip mass, and bending inertia of the beam. The DEB method proved to be very accurate, and in msot cases, was more accurate than the linear Taylor series approximation. The method is applicable to simultaneous perturbation of several design variables. Also, the approximations may be used to calculate other system response quantities. For example, the approximations for displacement are used to approximate bending stresses.

Pritchard, Jocelyn I.↗
NASA NTRS1990

Computation of maximum gust loads in nonlinear aircraft using a new method based on the matched filter approach and numerical optimization

Time-correlated gust loads are time histories of two or more load quantities due to the same disturbance time history. Time correlation provides knowledge of the value (magnitude and sign) of one load when another is maximum. At least two analysis methods have been identified that are capable of computing maximized time-correlated gust loads for linear aircraft. Both methods solve for the unit-energy gust profile (gust velocity as a function of time) that produces the maximum load at a given location on a linear airplane. Time-correlated gust loads are obtained by re-applying this gust profile to the airplane and computing multiple simultaneous load responses. Such time histories are physically realizable and may be applied to aircraft structures. Within the past several years there has been much interest in obtaining a practical analysis method which is capable of solving the analogous problem for nonlinear aircraft. Such an analysis method has been the focus of an international committee of gust loads specialists formed by the U.S. Federal Aviation Administration and was the topic of a panel discussion at the Gust and Buffet Loads session at the 1989 SDM Conference in Mobile, Alabama. The kinds of nonlinearities common on modern transport aircraft are indicated. The Statical Discrete Gust method is capable of being, but so far has not been, applied to nonlinear aircraft. To make the method practical for nonlinear applications, a search procedure is essential. Another method is based on Matched Filter Theory and, in its current form, is applicable to linear systems only. The purpose here is to present the status of an attempt to extend the matched filter approach to nonlinear systems. The extension uses Matched Filter Theory as a starting point and then employs a constrained optimization algorithm to attack the nonlinear problem.

Pototzky, Anthony S.↗
NASA NTRS1990

A decentralized linear quadratic control design method for flexible structures

A decentralized suboptimal linear quadratic control design procedure which combines substructural synthesis, model reduction, decentralized control design, subcontroller synthesis, and controller reduction is proposed for the design of reduced-order controllers for flexible structures. The procedure starts with a definition of the continuum structure to be controlled. An evaluation model of finite dimension is obtained by the finite element method. Then, the finite element model is decomposed into several substructures by using a natural decomposition called substructuring decomposition. Each substructure, at this point, still has too large a dimension and must be reduced to a size that is Riccati-solvable. Model reduction of each substructure can be performed by using any existing model reduction method, e.g., modal truncation, balanced reduction, Krylov model reduction, or mixed-mode method. Then, based on the reduced substructure model, a subcontroller is designed by an LQ optimal control method for each substructure independently. After all subcontrollers are designed, a controller synthesis method called substructural controller synthesis is employed to synthesize all subcontrollers into a global controller. The assembling scheme used is the same as that employed for the structure matrices. Finally, a controller reduction scheme, called the equivalent impulse response energy controller (EIREC) reduction algorithm, is used to reduce the global controller to a reasonable size for implementation. The EIREC reduced controller preserves the impulse response energy of the full-order controller and has the property of matching low-frequency moments and low-frequency power moments. An advantage of the substructural controller synthesis method is that it relieves the computational burden associated with dimensionality. Besides that, the SCS design scheme is also a highly adaptable controller synthesis method for structures with varying configuration, or varying mass and stiffness properties.

Su, Tzu-Jeng↗
NASA NTRS1991

The L sub 1 finite element method for pure convection problems

The least squares (L sub 2) finite element method is introduced for 2-D steady state pure convection problems with smooth solutions. It is proven that the L sub 2 method has the same stability estimate as the original equation, i.e., the L sub 2 method has better control of the streamline derivative. Numerical convergence rates are given to show that the L sub 2 method is almost optimal. This L sub 2 method was then used as a framework to develop an iteratively reweighted L sub 2 finite element method to obtain a least absolute residual (L sub 1) solution for problems with discontinuous solutions. This L sub 1 finite element method produces a nonoscillatory, nondiffusive and highly accurate numerical solution that has a sharp discontinuity in one element on both coarse and fine meshes. A robust reweighting strategy was also devised to obtain the L sub 1 solution in a few iterations. A number of examples solved by using triangle and bilinear elements are presented.

Jiang, Bo-Nan↗
NASA NTRS1990

Accuracy analysis of the semi-analytical method for shape sensitivity calculation

The semianalytical method of design sensitivity analysis that is widely used for calculating derivatives of static response with respect to design variables for structures modeled by finite elements is studied. It is shown that the method can have serious accuracy problems for shape design variables in structures modeled by beam, plate, truss, frame, and solid elements. Errors are shown to be associated with an incompatibility of the sensitivity field with the structure. An error index is developed to test the accuracy of the semianalytical method. It characterizes the difference in errors between a general finite-difference method and the semianalytical method. A method for improving the accuracy of the semianalytical method (when possible) is provided. Examples are presented to demonstrate the use of the error index.

Barthelemy, Bruno↗
NASA NTRS1992

Divergence preserving discrete surface integral methods for Maxwell's curl equations using non-orthogonal unstructured grids

Several new discrete surface integral (DSI) methods for solving Maxwell's equations in the time-domain are presented. These methods, which allow the use of general nonorthogonal mixed-polyhedral unstructured grids, are direct generalizations of the canonical staggered-grid finite difference method. These methods are conservative in that they locally preserve divergence or charge. Employing mixed polyhedral cells, (hexahedral, tetrahedral, etc.) these methods allow more accurate modeling of non-rectangular structures and objects because the traditional stair-stepped boundary approximations associated with the orthogonal grid based finite difference methods can be avoided. Numerical results demonstrating the accuracy of these new methods are presented.

Madsen, Niel K.↗
NASA NTRS1992

Comparison of epifluorescent viable bacterial count methods

Two methods, the 2-(4-Iodophenyl) 3-(4-nitrophenyl) 5-phenyltetrazolium chloride (INT) method and the direct viable count (DVC), were tested and compared for their efficiency for the determination of the viability of bacterial populations. Use of the INT method results in the formation of a dark spot within each respiring cell. The DVC method results in elongation or swelling of growing cells that are rendered incapable of cell division. Although both methods are subjective and can result in false positive results, the DVC method is best suited to analysis of waters in which the number of different types of organisms present in the same sample is assumed to be small, such as processed waters. The advantages and disadvantages of each method are discussed.

Rodgers, E. B.↗
NASA NTRS1992

Applications of a direct/iterative design method to complex transonic configurations

The current study explores the use of an automated direct/iterative design method for the reduction of drag in transport configurations, including configurations with engine nacelles. The method requires the user to choose a proper target-pressure distribution and then develops a corresponding airfoil section. The method can be applied to two-dimensional airfoil sections or to three-dimensional wings. The three cases that are presented show successful application of the method for reducing drag from various sources. The first two cases demonstrate the use of the method to reduce induced drag by designing to an elliptic span-load distribution and to reduce wave drag by decreasing the shock strength for a given lift. In the second case, a body-mounted nacelle is added and the method is successfully used to eliminate increases in wing drag associated with the nacelle addition by designing to an arbitrary pressure distribution as a result of the redesigning of a wing in combination with a given underwing nacelle to clean-wing, target-pressure distributions. These cases illustrate several possible uses of the method for reducing different types of drag. The magnitude of the obtainable drag reduction varies with the constraints of the problem and the configuration to be modified.

Smith, Leigh Ann↗
NASA NTRS1992

The solution of coupled Schroedinger equations using an extrapolation method

In this paper, extrapolation to the limit in a finite-difference method is applied to solve a system of coupled Schroedinger equations. This combination results in a method that only requires knowledge of the potential energy functions for the system. This numerical procedure has several distinct advantages over the more conventional methods. Namely, initial guesses for the term values are not needed; assumptions need be made about the behavior of the wavefunctions, such as the slope or magnitude in the nonclassical region; and the algorithm is easy to implement, has a firm mathematical foundation, and provides error estimates. Moreover, the method is less sensitive to round-off error than other methods since a small number of mesh points is used and it can be implemented on small computers. A comparison of the method with another numerical method shows results agreeing within 1 part in 10 exp 4.

Goorvitch, D.↗
NASA NTRS1992

The structures, binding energies and vibrational frequencies of Ca3 and Ca4: An application of the CCSD(T) method

The Ca3 and Ca4 metallic clusters have been investigated using state-of-the-art ab initio quantum mechanical methods. Large atomic natural orbital basis sets have been used in conjunction with the singles and doubles coupled-cluster (CCSD) method, a coupled-cluster method that includes a perturbational estimate of connected triple excitations, denoted CCSD(T), and the multireference configuration interaction (MRCI) method. The equilibrium geometries, binding energies and harmonic vibrational frequencies have been determined with each of the methods so that the accuracy of the coupled-cluster methods may be assessed. Since the CCSD(T) method reproduces the MRCI results very well, cubic and quartic force fields of Ca3 and Ca4 have been determined using this approach and used to evaluate the fundamental vibrational frequencies. The infrared intensities of both the e' mode of Ca3 and the t2 mode of Ca4 are found to be small. The results obtained in this study are compared and contrasted with those from our earlier studies on small Be and Mg clusters.

Lee, Timothy J.↗
NASA NTRS1992

An alternate method for achieving temperature control in the -130 C to 75 C range

Thermal vacuum testing often requires temperature control of chamber shrouds and heat exchangers within the -130 C to 75 C range. There are two conventional methods which are normally employed to achieve control through this intermediate temperature range: (1) single-pass flow where control is achieved by alternately pulsing hot gaseous nitrogen (GN2) and cold LN2 into the feed line to yield the setpoint temperature; and (2) closed-loop circulation where control is achieved by either electrically heating or LN2 cooling the circulating GN2 to yield the setpoint temperature. A third method, using a mass flow ratio controller along with modulating control valves on GN2 and LN2 lines, provides excellent control but equipment for this method is expensive and cost-prohibitive for all but long-term continuous processes. The single-pass method provides marginal control and can result in unexpected overcooling of the test article from even a short pulse of LN2. The closed-loop circulation method provides excellent control but requires an expensive blower capable of operating at elevated pressures and cryogenic temperatures. Where precise control is needed (plus or minus 2 C), single-pass flow systems typically have not provided the precision required, primarily because of overcooling temperature excursions. Where several individual circuits are to be controlled at different temperatures, the use of expensive cryogenic blowers for each circuit is also cost-prohibitive, especially for short duration of one-of-a-kind tests. At JPL, a variant of the single-pass method was developed that was shown to provide precise temperature control in the -130 C to 75 C range while exhibiting minimal setpoint overshoot during temperature transitions. This alternate method uses a commercially available temperature controller along with a GN2/LN2 mixer to dampen the amplitude of cold temperature spikes caused by LN2 pulsing. The design of the GN2/LN2 mixer, the overall control system configuration, the operational procedure, and the prototype system test results are described.

Johnson, Kenneth R.↗