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 solution of the one dimensional planar interface solidification problem is presented. Numerical solutions are applied to experimental solute profiles to determine the effective diffusion constant for Hg Co Te. Progress on the design and construction of a high gradient furnace is reported.
The indirect synthesis method is developed and shown to be capable of leading a near-optimal design of multidegree-of-freedom and multidesign-element transient nonlinear dynamical systems. The basis of the approach is to select the open design parameters such that the response of the portion of the system being designed approximates the limiting performances solution. The limiting performance problem can be formulated as one of linear programming by replacing all portions of the system subject to transient disturbances by control forces and supposing that the remaining portions are linear as are the overall kinematic constraints. One then selects the design parameters that respond most closely to the limiting performance solution, which can be achieved by unconstrained curve-fitting techniques.
Methods for calculating sensitivity derivatives for discrete structural systems are surveyed, primarily covering literature published during the past two decades. Methods are described for calculating derivatives of static displacements and stresses, eigenvalues and eigenvectors, transient structural response, and derivatives of optimum structural designs with respect to problem parameters. The survey is focused on publications addressed to structural analysis, but also includes a number of methods developed in nonstructural fields such as electronics, controls, and physical chemistry which are directly applicable to structural problems. Most notable among the nonstructural-based methods are the adjoint variable technique from control theory, and the Green's function and FAST methods from physical chemistry.
First-order, second-order, and implicit total variation diminishing (TVD) schemes are reviewed using the modified flux approach. Some transient and steady-state calculations are then carried out to illustrate the applicability of these schemes to the Euler equations. It is shown that the second-order explicit TVD schemes generate good shock resolution for both transient and steady-state one-dimensional and two-dimensional problems. Numerical experiments for a quasi-one-dimensional nozzle problem show that the second-order implicit TVD scheme produces a fairly rapid convergence rate and remains stable even when running with a Courant number of 10 to the 6th.
Methods for calculating sensitivity derivatives for discrete structural systems are surveyed, primarily covering literature published during the past two decades. Methods are described for calculating derivatives of static displacements and stresses, eigenvalues and eigenvectors, transient structural response, and derivatives of optimum structural designs with respect to problem parameters. The survey is focused on publications developed in nonstructural fields such as electronics, controls, and physical chemistry which are directly applicable to structural problems. Most notable among the nonstructural-based methods are the adjoint variable technique from control theory, and the Green's function and FAST methods from physical chemistry.
A general three-field variational principle is obtained for the motion of an acoustic fluid enclosed in a rigid or flexible container by the method of canonical decomposition applied to a modified form of the wave equation in the displacement potential. The general principle is specialized to a mixed two-field principle that contains the fluid displacement potential and pressure as independent fields. This principle contains a free parameter alpha. Semidiscrete finite-element equations of motion based on this principle are displayed and applied to the transient response and free-vibrations of the coupled fluid-structure problem. It is shown that a particular setting of alpha yields a rich set of formulations that can be customized to fit physical and computational requirements. The variational principle is then extended to handle slosh motions in a uniform gravity field, and used to derived semidiscrete equations of motion that account for such effects.
A general three-field variational principle is obtained for the motion of an acoustic fluid enclosed in a rigid or flexible container by the method of canonical decomposition applied to a modified form of the wave equation in the displacement potential. The general principle is specialized to a mixed two-field principle that contains the fluid displacement potential and pressure as independent fields. This principle contains a free parameter alpha. Semidiscrete finite-element equations of motion based on this principle are displayed and applied to the transient response and free-vibrations of the coupled fluid-structure problem. It is shown that a particular setting of alpha yields a rich set of formulations that can be customized to fit physical and computational requirements. The variational principle is then extended to handle slosh motions in a uniform gravity field, and used to derive semidiscrete equations of motion that account for such effects.
The phenomenon of coronal mass ejections (CME) is examined with reference to recent work in this field. In particular, attention is given to CMEs close to the sun, CMEs in interplanetary space, solar cycle variations of CMEs, and magnetic field topology and the problem of magnetic flux balance in interplanetary space. It is demonstrated that transient ejections of material from the sun are the cause of almost all large nonrecurrent geomagnetic substorms, as originally suggested by Chapman and his colleagues.
A fundamental problem in Solar-Terrestrial Physics is the origin of the solar transient plasma output, which includes the coronal mass ejection and its interplanetary manifestation, e.g. the magnetic cloud. The traditional blast wave model resulted from solar thermal pressure impulse has faced with challenge during recent years. In the MHD numerical simulation study of CME, the authors find that the basic feature of the asymmetrical event on 18 August 1980 can be reproduced neither by a thermal pressure nor by a speed increment. Also, the thermal pressure model fails in simulating the interplanetary structure with low thermal pressure and strong magnetic field strength, representative of a typical magnetic cloud. Instead, the numerical simulation results are in favor of the magnetic field expansion as the likely mechanism for both the asymmetrical CME event and magnetic cloud.
In a continuing program of research and development, a system has been demonstrated that makes high-speed measurements of thermal infrared radiance from gas-turbine engine exhaust streams. When a gas-turbine engine is operated under conditions that minimize the emission of pollutants, there is a risk of crossing the boundary from stable to unstable combustion. Combustion instability can lead to engine damage and even catastrophic failure. Sensor systems of the type under development could provide valuable data during the development testing of gas-turbine engines or of engine components. A system of the type under development makes high-speed measurements of thermal infrared radiance from the engine exhaust stream. The sensors of this system can be mounted outside the engine, which eliminates the need for engine case penetrations typical with other engine dynamics monitors. This is an important advantage in that turbine-engine manufacturers consider such penetrations to be very undesirable. A prototype infrared sensor system has been built and demonstrated on a turbine engine. This system includes rugged and inexpensive near-infrared sensors and filters that select wavelengths of infrared radiation for high sensitivity. In experiments, low-frequency signatures were consistently observed in the detector outputs. Under some conditions, the signatures also included frequency components having one or two radiance cycles per engine revolution. Although it has yet to be verified, it is thought that the low-frequency signatures may be associated with bulk-mode combustion instabilities or flow instabilities in the compressor section of the engine, while the engine- revolution-related signatures may be indicative of mechanical problems in the engine. The system also demonstrated the ability to detect transient high-radiance events. These events indicate hot spots in the exhaust stream and were found to increase in frequency during engine acceleration.
Two methods for calculating linear frequency domain aerodynamic coefficients from a time marching Full Potential cascade solver are developed and verified. In the first method, the Influence Coefficient, solutions to elemental problems are superposed to obtain the solutions for a cascade in which all blades are vibrating with a constant interblade phase angle. The elemental problem consists of a single blade in the cascade oscillating while the other blades remain stationary. In the second method, the Pulse Response, the response to the transient motion of a blade is used to calculate influence coefficients. This is done by calculating the Fourier Transforms of the blade motion and the response. Both methods are validated by comparison with the Harmonic Oscillation method and give accurate results. The aerodynamic coefficients obtained from these methods are used for frequency domain flutter calculations involving a typical section blade structural model. An eigenvalue problem is solved for each interblade phase angle mode and the eigenvalues are used to determine aeroelastic stability. Flutter calculations are performed for two examples over a range of subsonic Mach numbers.
Two methods for calculating linear frequency domain aerodynamic coefficients from a time marching Full Potential cascade solver are developed and verified. In the first method, the Influence Coefficient, solutions to elemental problems are superposed to obtain the solutions for a cascade in which all blades are vibrating with a constant interblade phase angle. The elemental problem consists of a single blade in the cascade oscillating while the other blades remain stationary. In the second method, the Pulse Response, the response to the transient motion of a blade is used to calculate influence coefficients. This is done by calculating the Fourier Transforms of the blade motion and the response. Both methods are validated by comparison with the Harmonic Oscillation method and give accurate results. The aerodynamic coefficients obtained from these methods are used for frequency domain flutter calculations involving a typical section blade structural model. An eigenvalue problem is solved for each interblade phase angle mode and the eigenvalues are used to determine aeroelastic stability. Flutter calculations are performed for two examples over a range of subsonic Mach numbers.
A fault model and a system model are necessary to assess the tolerance or resilience of a digital system to lightning-induced transients. It is noted that these models are usually developed separately. A new approach is outlined here for this assessment problem which combines the fault and system models into an overall model. With this approach, referred to as the containment set approach, an assessment of the effects of lightning-induced transients can be made in terms of a state transition matrix. This matrix can be generated by means of fault injection experiments. In addition, certain nonredundancy-oriented design alternatives to the achievement of lightning-induced transient tolerance in digital systems are indicated by the containment set approach.
The extent to which suitable solutions may be obtained for one physics problem and two engineering type problems is traced. NASTRAN appears to be a practical tool to solve one-group steady-state neutron diffusion equations. Transient diffusion analysis may be performed after new levels that allow time-dependent temperature calculations are developed. NASTRAN piecewise linear anlaysis may be applied to solve those plasticity problems for which a smooth stress-strain curve can be used to describe the nonlinear material behavior. The accuracy decreases when sharp transitions in the stress-strain relations are involved. Improved NASTRAN usefulness will be obtained when nonlinear material capabilities are extended to axisymmetric elements and to include provisions for time-dependent material properties and creep analysis. Rigid formats 3 and 5 proved to be very convenient for the buckling and normal-mode analysis of a nuclear fuel element.
Implicit solution methods are important in applications modeled by PDEs with disparate temporal and spatial scales. Because such applications require high resolution with reasonable turnaround, "routine" parallelization is essential. The pseudo-transient matrix-free Newton-Krylov-Schwarz (Psi-NKS) algorithmic framework is presented as an answer. We show that, for the classical problem of three-dimensional transonic Euler flow about an M6 wing, Psi-NKS can simultaneously deliver: globalized, asymptotically rapid convergence through adaptive pseudo- transient continuation and Newton's method-, reasonable parallelizability for an implicit method through deferred synchronization and favorable communication-to-computation scaling in the Krylov linear solver; and high per- processor performance through attention to distributed memory and cache locality, especially through the Schwarz preconditioner. Two discouraging features of Psi-NKS methods are their sensitivity to the coding of the underlying PDE discretization and the large number of parameters that must be selected to govern convergence. We therefore distill several recommendations from our experience and from our reading of the literature on various algorithmic components of Psi-NKS, and we describe a freely available, MPI-based portable parallel software implementation of the solver employed here.
A theoretical overview of the response coefficient method for heat conduction transients with time-dependent input forcing functions is presented with a number of illustrative applications. The method may be the most convenient and economical if the same problem is to be solved many times with different input-time histories or if the solution time is relatively long. The method is applicable to a wide variety of problems, including irregular geometries, position-dependent boundary conditions, position-dependent physical properties, and nonperiodic irregular input histories. Nonuniform internal energy generation rates within the structure can also be handled by the method. The area of interest is long-time solutions, in which initial condition is unimportant, and not the early transient period. The method can be applied to one dimensional problems in cartesian, cylindrical, and spherical coordinates as well as to two dimensional problems in cartesian and cylindrical coordinates.
One of the propulsion system concepts to be considered for the High-Speed Civil Transport (HSCT) is an underwing, dual-propulsion, pod-per-wing installation. Adverse transient phenomena such as engine compressor stall and inlet unstart could severely degrade the performance of one of these propulsion pods. The subsequent loss of thrust and increased drag could cause aircraft stability and control problems that could lead to a catastrophic accident if countermeasures are not in place to anticipate and control these detrimental transient events. Aircraft system engineers must understand what happens during an engine compressor stall and inlet unstart so that they can design effective control systems to avoid and/or alleviate the effects of a propulsion pod engine compressor stall and inlet unstart. The objective of the Inlet Unstart Propulsion Airframe Integration test program was to assess the underwing flow field of a High-Speed Civil Transport propulsion system during an engine compressor stall and subsequent inlet unstart. Experimental research testing was conducted in the 10- by 10-Foot Supersonic Wind Tunnel at the NASA Glenn Research Center at Lewis Field. The representative propulsion pod consisted of a two-dimensional, bifurcated inlet mated to a live turbojet engine. The propulsion pod was mounted below a large flat plate that acted as a wing simulator. Because of the plate s long length (nominally 10-ft wide by 18-ft long), realistic boundary layers could form at the inlet cowl plane. Transient instrumentation was used to document the aerodynamic flow-field conditions during an unstart sequence. Acquiring these data was a significant technical challenge because a typical unstart sequence disrupts the local flow field for about only 50 msec. Flow surface information was acquired via static pressure taps installed in the wing simulator, and intrusive pressure probes were used to acquire flow-field information. These data were extensively analyzed to determine the impact of the unstart transient on the surrounding flow field. This wind tunnel test program was a success, and for the first time, researchers acquired flow-field aerodynamic data during a supersonic propulsion system engine compressor stall and inlet unstart sequence. In addition to obtaining flow-field pressure data, Glenn researchers determined other properties such as the transient flow angle and Mach number. Data are still being reduced, and a comprehensive final report will be released during calendar year 2000.
Performability analysis is a technique developed for evaluating the effectiveness of fault-tolerant computing systems in multiphase missions. Performability was evaluated for its accuracy, practical usefulness, and relative cost. The evaluation was performed by applying performability and the fault tree method to a set of sample problems ranging from simple to moderately complex. The problems involved as many as five outcomes, two to five mission phases, permanent faults, and some functional dependencies. Transient faults and software errors were not considered. A different analyst was responsible for each technique. Significantly more time and effort were required to learn performability analysis than the fault tree method. Performability is inherently as accurate as fault tree analysis. For the sample problems, fault trees were more practical and less time consuming to apply, while performability required less ingenuity and was more checkable. Performability offers some advantages for evaluating very complex problems.