SEARCH · Engineering Papers
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.
The finite element method for fluid dynamics was used to develop a two dimensional mesh generation scheme. The method consists of shearing and conformal maps with upper and lower surfaces handled independently to allow sharp leading edges. The method also generates meshes of triangular or quadrilateral elements.
The study of time-dependent flows generated by surface-tension gradients in thin liquid layers is reported. Numerical techniques for dealing with the complete float-zone process were investigated. In particular, finite element methods with fast iterative procedures, using multipoint constraints and nonlinear optimization were examined.
An assessment of the potential impact of adaptive finite element technology on the analysis part of the aircraft structural synthesis process is presented. The main conclusion is that adaptive application of the p-version of the finite element method based on indirect error estimation procedures results in substantial cost reduction and increased reliability of the computed data. Adaptivity based on direct a posteriori error estimation has the potential for additional savings.
Reduced models and reduced controllers for systems governed by matrix-second-order differential equations are obtained by retaining those modes which make the largest contributions to quadratic control objectives. Such contributions, expressed in terms of modal data, used as mode truncation criteria, allow the statement of the specific control objectives to influence the early model reduction from very high order models which are available, for example, from finite element methods. The relative importance of damping, frequency, and eigenvector in the mode truncation decisions are made explicit for each of these control objectives: attitude control, vibration suppression and figure control. The paper also shows that using modal cost analysis (MCA) on the closed loop modes of the optimally controlled system allows the construction of reduced control policies which feedback only those closed loop modal coordinates which are most critical to the quadratic control performance criterion. In this way, the modes which should be controlled (and hence the modes which must be observable by choice of measurements), are deduced from truncations of the optimal controller.
A general mapping procedure is described and applied to the study of noise propagation in variable area ducts. The mapping provides a boundary fitted co-ordinate system which is ideal for the finite difference solution of acoustic fields with irregular boundaries, without the burden of large matrices required by finite element methods. The procedure is first described in general and then applied to a particular two-dimensional geometry under current experimental investigation. This method should be ideally suited to the study of high frequency noise propagation in variable area ducts and in cases where the far field is included in the calculation procedure. Moreover, the current approach can be directly extended to three-dimensions, resulting in numerical calculation over a rectangular parallelepiped in the transformed plane.
Finite element method is employed in the solution of full-potential equation for analyzing steady, transonic flows through a series of highly staggered cascades. Convergence characteristics and efficiency of the developed numerical procedure are discussed. Accuracy of the mathematical model is investigated in terms of the employed computational grid. In particular, the importance of accurate representation of the trailing edge together with the shock on the obtained flow results are studied. Numerical results include choked flows, for a cascade of NACA 0012 airfoils with pitch/chord length = 1, and stagger angle of 45 deg.
A constitutive law by which structural damping is modeled as a relationship between stress, strain, and strain rate in a material is used in conjunction with the finite element method to develop general integral expressions for viscous and nonviscous damping matrices. To solve the set of nonlinear equations resulting from the presence of nonviscous damping, a solution technique is developed by modifying the Newmark method to accommodate an iterative solution and treat the nonviscous damping as a pseudo-force. The technique is then checked for accuracy and convergence in single- and multi-degree-of-freedom problems, and is found to be accurate and efficient for initial-condition problems with small nonviscous damping.
The problem of forced fluid vibrations in a partially filled spinning spherical tank is solved numerically by using the finite element method. The governing equations include Coriolis acceleration and spatially homogeneous vorticity. An exponential instability is detected in the present simulation for fill ratios below 0.5 and centrifugal acceleration to thrust ratios less than 1.7. This fictitious instability appears in the model as a result of the homogeneous vortex assumption since the free slosh equations are neutrally stable in the Liapunov sense.
An algorithm is developed for the numerical simulation of large motions of an unrestrained, nonuniform beam in orbit. Explicit provisions are made for the accommodation of general, physically significant initial conditions and for the specification of control forces in practical terms. The evaluation of certain constants that appear in the equations of motion is considered, the finite element method being used for this purpose. The importance of initial conditions in the study of spacecraft motions is stressed.
An investigation is conducted regarding the feasibility of approaches for improving the efficiency and stability of existing finite element method (FEM) schemes, taking into account both analytical and numerical studies. Of the four schemes considered for the 'steady' problem, the Hermitian Galerkin formulation appears to be the most efficient and therefore the most suitable scheme for futher full scale implementation. The Hermitian residual least squares (RLS) scheme although comparable in accuracy for the cases considered exhibits a slight tendency to cumulative errors. The performance of both the Lagrangian element schemes considered compares poorly with that of their Hermitian element counterparts. This is particularly true of the Lagrangian RLS scheme. The presence of internal oscillatory components is an inevitable consequence of all Galerkin schemes irrespective of element type.
Interlaminar crack growth behavior in fiber-reinforced composites subjected to fatigue loading was investigated experimentally and theoretically. In the experimental phase, inter-laminar crack propagation rates and mechanisms were determined for the cases of various geometries, laminate parameters and cyclic stress levels. A singular hybrid-stress finite element method was used in conjuction with the experimental results to examine the local crack-tip behavior and to characterize the crack propagation during fatigue. Results elucidate the basic nature of the cyclic delamination damage, and relate the interlaminar crack growth rate to the range of mixed-mode crack-tip stress intensity factors. The results show that crack growth rates are directly related to the range of the mixed-mode cyclic stress intensity factors by a power law relationship.
The acoustic field within a duct containing a nonuniform steady flow was predicted. This analysis used the finite element method to calculate the velocity potential within the duct.
An iterative finite element integral technique is used to predict the sound field radiated from the JT15D turbofan inlet. The sound field is divided into two regions: the sound field within and near the inlet which is computed using the finite element method and the radiation field beyond the inlet which is calculated using an integral solution technique. The velocity potential formulation of the acoustic wave equation was employed in the program. For some single mode JT15D data, the theory and experiment are in good agreement for the far field radiation pattern as well as suppressor attenuation. Also, the computer program is used to simulate flight effects that cannot be performed on a ground static test stand.
Numerical solution of nonlinear equilibrium problems of structures by means of Newton-Raphson type iterations is reviewed. Each step of the iteration is shown to correspond to the solution of a linear problem, therefore the feasibility of the finite element method for nonlinear analysis is established. Organization and flow of data for various types of digital computers, such as single-processor/single-level memory, single-processor/two-level-memory, vector-processor/two-level-memory, and parallel-processors, with and without sub-structuring (i.e. partitioning) are given. The effect of the relative costs of computation, memory and data transfer on substructuring is shown. The idea of assigning comparable size substructures to parallel processors is exploited. Under Cholesky type factorization schemes, the efficiency of parallel processing is shown to decrease due to the occasional shared data, just as that due to the shared facilities.
General nonlinear discretized governing equations of motion of spinning elastic solids and structures are adjusted for the case of a spinning cantilever with initial geometric imperfections. Consideration is given to second degree nonlinearities in the strain-displacement and velocity-displacement relationships. Parameters of the discretization are developed to include the type and number of the coordinate functions used in the admissible trial solution in order to unify the discretization approaches associated with stationarity principles. The coordinate functions comprise both sets of continuous and piecewise continuous functions employed in the Rayleigh-Ritz and the finite element methods, respectively. Coefficient matrices are provided which contain the energy density expressions and which are adaptable to computer programming.
Steady state crack propagation is investigated numerically using a model consisting of 236 free atoms in two (010) planes of bcc alpha iron. The continuum region is modeled using the finite element method with 175 nodes and 288 elements. The model shows clear (010) plane fracture to the edge of the discrete region at moderate loads. Analysis of the results obtained indicates that models of this type can provide realistic simulation of steady state crack propagation.
An application of the finite element method is presented in order to study the flow through a two-dimensional channel including the choked and nearly-choked flow conditions. The mathematical procedure provides a combined treatment of potential and Euler equations, where the steady Euler equations are solved through the integration of a pseudo-time system which is equivalent to a relaxation scheme; the isentropic, potential flow equations are embedded in this formulation. To analyze the transonic flow through a channel with near-choked flow conditions, the flow through a parallel channel with a 10 percent circular arc bump is calculated as a sample problem.
Gear tooth average and instantaneous surface temperatures were measured with a fast response infrared radiometric microscope, while operating at arious speeds, loads and oil jet pressures. Increased oil jet pressure had a significant effect on both average and peak surface temperatures at all test conditions, increasing the speed at constant load and increasing the load at constant speed causes a significant rise in average and peak surface temperatures of gear teeth. A gear tooth temperature analysis was conducted by a finite element method combined with a calculated heat input and oil jet impingment depth with estimated heat transfer coefficients based on the experimental data. It is concluded that oil jet pressures required for adequate cooling at high load and speed conditions must be high enough to get full penetration depth of the teeth.