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.
Element-by-element solution strategies are developed for transient heat conduction problems. Results of numerical tests indicate the effectiveness of the procedures proposed. The small database requirements and attractive architectural features of the algorithms suggest considerable potential for solving large scale problems.
The controls analysis based on a truncated finite element model of the 122m. Hoop/Column Antenna System focuses on an analysis of the controllability as well as the synthesis of control laws. Graph theoretic techniques are employed to consider controllability for different combinations of number and locations of actuators. Control law synthesis is based on an application of the linear regulator theory as well as pole placement techniques. Placement of an actuator on the hoop can result in a noticeable improvement in the transient characteristics. The problem of orientation and shape control of an orbiting flexible beam, previously examined, is now extended to include the influence of solar radiation environmental forces. For extremely flexible thin structures modification of control laws may be required and techniques for accomplishing this are explained. Effects of environmental torques are also included in previously developed models of orbiting flexible thin platforms.
The specifications of the flywheel system for momentum storage and vehicle torquing are somewhat dependent upon the attitude control requirements of the space station in orbit. As a ground rule, the flywheel system will be sized large enough to provide all attitude maneuvers, if practical, to avoid or minimize turning on the reaction control system (RCS). The RCS, whenever used, expels expensive mass and tends to contaminate optical surfaces of the vehicle. The vehicle rate and acceleration specifications of 0.10 deg/sec and 0.01 deg/square sec are tentative, and may be reduced if lesser values are more practical for flywheel design. For local vertical attitude hold, the average attitude error should be zero, and not the classical 1 degree, since control moment gyro (CMG) gimbal angles provide an exact reference feedback for gravity gradient momentum. Docking presents a problem for docking transients and attitude alignment which will require use of the RCS.
During final development testing of the Space Shuttle Reaction Control Subsystem propellant tanks, a problem with pressure transients in the system was uncovered. Due to the nature of the tanks, performance tests to determine the impact of the transients on the expulsion efficiency of the tanks could not directly simulate the actual conditions which would be present in a low-gravity environment. However, by masking or covering various segments of the propellant acquisition device, a good simulation of a low-gravity environment was achieved in ground testing.
The convergence properties of three modal methods for solving structural transient response (STR) problems are examined by comparing the accuracy of the total number of modes used to represent the responses. The goal is to identify a method which involves as few degrees of freedom as possible. Two mode superposition methods, the mode displacement method (MDM) and the mode acceleration method (MAM), are fitted with equivalent expressions for damped response. An error norm is employed to compare the accuracies of the methods in calculating displacements, moments, and shear forces for a viscously-damped cantilevered beam experiencing several dynamic loading conditions and levels of damping. The MAM converges quicker and uses fewer dof, mainly because it has a pseudo-static response term which accounts for some of the flexibility of higher modes. A third method, which integrates the convolution integral iteratively over time, generates a higher order approximation to the STR with a reduced number of modes. It converges fastest among the three, but requires definition of an appropriate forcing function. Finally, the MAM and MDM models are applied to study the effects of the forcing frequency and the initial conditions.
Techniques for eliminating time lags in the implicit finite-difference solution of partial differential equations are investigated analytically, with a focus on transient fluid dynamics problems on overlapping multicomponent grids. The fundamental principles of the approach are explained, and the method is shown to be applicable to both rectangular and curvilinear grids. Numerical results for sample problems are compared with exact solutions in graphs, and good agreement is demonstrated.
Report discusses use of explicit exponential finite-difference technique to solve various diffusion-type partial differential equations. Study extends technique to transient-heat-transfer problems in one dimensional cylindrical coordinates and two and three dimensional Cartesian coordinates and to some nonlinear problems in one or two Cartesian coordinates.
The aim was to study the solution of transient finite element problems on the Transputer system of parallel processors. The central difference time integration rule was used so that no equation solving was necessary. Also investigated was subcycling time integration which uses different time steps in different subdomains of the finite element mesh. A one-dimensional bar problem was analyzed using the parallel time integration algorithm. This involves subdividing the bar into subproblems which are assigned to different processors. Results show that the significant speed-up can be obtained through parallel processing. Also subcycling can give an additional speed-up in certain classes of problems. A two-dimensional problem was also examined to evaluate the effect of the communication to computation ratio on solution time.
Most high speed cryogenic turbopumps for liquid propulsion rocket engines currently use ball or roller contact bearings for rotor support. The operating speeds, loads, clearances, and environments of these pumps combine to make bearing wear a limiting factor on turbopump life. An example is the high pressure oxygen turbopump (HPOTP) used in the Space Shuttle Main Engine (SSME). Although the HPOTP design life is 27,000 seconds at 30,000 rpms, or approximately 50 missions, bearings must currently be replaced after 2 missions. One solution to the bearing wear problem in the HPOTP, as well as in future turbopump designs, is the utilization of fluid film bearings in lieu of continuous contact bearings. Hydrostatic, hydrodynamic, and damping seal bearings are all replacement candidates for contact bearings in rocket engine high speed turbomachinery. These three types of fluid film bearings have different operating characteristics, but they share a common set of mechanical design opportunities and difficulties. Results of research to define some of the mechanical design issues are given. Problems considered include transient strat/stop rub, non-operational rotor support, bearing wear inspection and measurement, and bearing fluid supply route. Emphasis is given to the HPOTP preburner pump (PBP) bearing, but the results are pertinent to high-speed cryogenic turbomachinery in general.
A parallel program to analyze transient finite element problems was written and implemented on a system of transputer processors. The program uses the explicit time integration algorithm which eliminates the need for equation solving, making it more suitable for parallel computations. An interprocessor communication scheme was developed for arbitrary two dimensional grid processor configurations. Several 3-D problems were analyzed on a system with a small number of processors.
The useful and powerful features of NASTRAN and three real world problems for the testing of the capabilities of different NASTRAN versions are discussed. The test problems involve direct transient analysis, nonlinear analysis, and static analysis. The experiences in using graphics software packages are also discussed. It was found that MSC/XL can be more useful if it can be improved to generate picture files of the analysis results and to extend its capabilities to support finite element codes other than MSC/NASTRAN. It was found that the current version of SDRC/I-DEAS (version VI) may have bugs in the module 'Data Loader'.
The application of adaptive finite element methods to the solution of transient heat conduction problems in two dimensions is investigated. The computational domain is represented by an unstructured assembly of linear triangular elements and the mesh adaptation is achieved by local regeneration of the grid, using an error estimation procedure coupled to an automatic triangular mesh generator. Two alternative solution procedures are considered. In the first procedure, the solution is advanced by explicit timestepping, with domain decomposition being used to improve the computational efficiency of the method. In the second procedure, an algorithm for constructing continuous lines which pass only once through each node of the mesh is employed. The lines are used as the basis of a fully implicit method, in which the equation system is solved by line relaxation using a block tridiagonal equation solver. The numerical performance of the two procedures is compared for the analysis of a problem involving a moving heat source applied to a convectively cooled cylindrical leading edge.
New developments describing the theoretical basis towards an effective virtual-pulse time integral procedure are presented for general nonlinear transient heat transfer problems. To validate the proposed methodology of computation, numerical test examples are given and comparisons are drawn with the implicit second-order accurate Crank-Nicolson technique. The proposed method shows superior or improved accuracy and stability characteristics for the models tested.
The behavior of a turbulent boundary layer which has been subjected to a local ramp-like deceleration in the external velocity field, which leads to forced separation, has been studied experimentally. The data of this study are re-interpreted in light of more recent findings concerning the temporal nature of boundary layer turbulence in the presence of forced unsteady shear. In particular, the robustness of the near-wall turbulent motions to organized deformation is recognized. Their resilence during unsteady shearing action promotes continued efficient turbulent mixing and rapid redistribution of turbulent kinetic energy during forced transients. In aerodynamic problems, the rapid nature of the adjustment of the turbulence field to a new temporal boundary condition necessitates equally rapid remedial measures to be taken if means of control/prevention of forced unsteady separation are to be deployed to maximum effect. This requirement suggests exploration of the use of simple, real-time statistical forecasting techniques, based upon time-series analysis of easily-measurable features of the flow, to help assure timely deployment of mechanisms of boundary-layer control. This paper focuses upon the nature of turbulence in boundary layers undergoing forced deceleration which would lead to separation. A preliminary form of a forecasting model is presented and evaluated. Using observations of the previous two large eddies passing a detector, it forecasts the behavior of the future large eddy rather well.
The focus of this research is on the development of analysis and sensitivity analysis equations for nonlinear, transient heat transfer problems modeled by p-version, time discontinuous finite element approximation. The resulting matrix equation of the state equation is simply in the form ofA(x)x = c, representing a single step, time marching scheme. The Newton-Raphson's method is used to solve the nonlinear equation. Examples are first provided to demonstrate the accuracy characteristics of the resultant finite element approximation. A direct differentiation approach is then used to compute the thermal sensitivities of a nonlinear heat transfer problem. The report shows that only minimal coding effort is required to enhance the analysis code with the sensitivity analysis capability.
The heat transfer of a water spray impinging upon a surface in a very low pressure environment is of interest to cooling of space vehicles during launch and re-entry, and to industrial processes where flash evaporation occurs. At very low pressure, the process occurs near the triple point of water, and there exists a transient multiphase transport problem of ice, water and water vapor. At the impingement location, there are three heat transfer mechanisms: evaporation, freezing and sublimation. A preliminary heat transfer model was developed to explore the interaction of these mechanisms at the surface and within the spray.
A current research topic in coupled-field problems is the development of effective transient algorithms that permit different time integration methods with different time steps to be used simultaneously in various regions of the problems. The implicit-explicit approach seems to be very successful in structural, fluid, and fluid-structure problems. This paper summarizes this research direction. A family of mixed time integration schemes, with the capabilities mentioned above, is also introduced for transient thermal analysis. A stability analysis and the computer implementation of this technique are also presented. In particular, it is shown that the mixed time implicit-explicit methods provide a natural framework for the further development of efficient, clean, modularized computer codes.
The classical EHL (elastohydrodynamic lubrication) point contact problem is solved using a new system-approach, similar to that introduced by Houpert and Hamrock for the line-contact problem. Introducing a body-fitted coordinate system, the troublesome free-boundary is transformed to a fixed domain. The Newton-Raphson method can then be used to determine the pressure distribution and the cavitation boundary subject to the Reynolds boundary condition. This method provides an efficient and rigorous way of solving the EHL point contact problem with the aid of a supercomputer and a promising method to deal with the transient EHL point contact problem. A typical pressure distribution and film thickness profile are presented and the minimum film thicknesses are compared with the solution of Hamrock and Dowson. The details of the cavitation boundaries for various operating parameters are discussed.