Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transient problem”

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 397 records · Page 22

Large-angle slewing maneuvers for flexible spacecraft

A new class of closed-form solutions for finite-time linear-quadratic optimal control problems is presented. The solutions involve Potter's solution for the differential matrix Riccati equation, which assumes the form of a steady-state plus transient term. Illustrative examples are presented which show that the new solutions are more computationally efficient than alternative solutions based on the state transition matrix. As an application of the closed-form solutions, the neighboring extremal path problem is presented for a spacecraft retargeting maneuver where a perturbed plant with off-nominal boundary conditions now follows a neighboring optimal trajectory. The perturbation feedback approach is further applied to three-dimensional slewing maneuvers of large flexible spacecraft. For this problem, the nominal solution is the optimal three-dimensional rigid body slew. The perturbation feedback then limits the deviations from this nominal solution due to the flexible body effects. The use of frequency shaping in both the nominal and perturbation feedback formulations reduces the excitation of high-frequency unmodeled modes. A modified Kalman filter is presented for estimating the plant states.

Chun, Hon M.↗

Single Event Transients in Linear Integrated Circuits

On November 5, 2001, a processor reset occurred on board the Microwave Anisotropy Probe (MAP), a NASA mission to measure the anisotropy of the microwave radiation left over from the Big Bang. The reset caused the spacecraft to enter a safehold mode from which it took several days to recover. Were that to happen regularly, the entire mission would be compromised, so it was important to find the cause of the reset and, if possible, to mitigate it. NASA assembled a team of engineers that included experts in radiation effects to tackle the problem. The first clue was the observation that the processor reset occurred during a solar event characterized by large increases in the proton and heavy ion fluxes emitted by the sun. To the radiation effects engineers on the team, this strongly suggested that particle radiation might be the culprit, particularly when it was discovered that the reset circuit contained three voltage comparators (LM139). Previous testing revealed that large voltage transients, or glitches appeared at the output of the LM139 when it was exposed to a beam of heavy ions [NI96]. The function of the reset circuit was to monitor the supply voltage and to issue a reset command to the processor should the voltage fall below a reference of 2.5 V [PO02]. Eventually, the team of engineers concluded that ionizing particle radiation from the solar event produced a negative voltage transient on the output of one of the LM139s sufficiently large to reset the processor on MAP. Fortunately, as of the end of 2004, only two such resets have occurred. The reset on MAP was not the first malfunction on a spacecraft attributed to a transient. That occurred shortly after the launch of NASA s TOPEX/Poseidon satellite in 1992. It was suspected, and later confirmed, that an anomaly in the Earth Sensor was caused by a transient in an operational amplifier (OP-15) [KO93]. Over the next few years, problems on TDRS, CASSINI, [PR02] SOHO [HA99,HA01] and TERRA were also attributed to transients. In some cases, such events produced resets by falsely triggering circuits designed to protect against over- voltage or over-current. On at least three occasions, transients caused satellites to switch into "safe mode" in which most of the systems on board the satellites were powered down for an extended period. By the time the satellites were reconfigured and returned to full operational state, much scientific data had been lost. Fortunately, no permanent damage occurred in any of the systems and they were all successfully re-activated.

Buchner, Stephen↗

The initial-value problem for viscous channel flows

Plane viscous channel flows are perturbed and the ensuing initial-value problems are investigated in detail. Unlike traditional methods where traveling wave normal modes are assumed for solution, this works offers a means whereby completely arbitrary initial input can be specified without having to resort to eigenfunction expansions. The full temporal behavior, including both early time transients and the long time asymptotics, can be determined for any initial disturbance. Effects of three-dimensionality can be assessed. The bases for the analysis are: (a) linearization of the governing equations; (b) Fourier decomposition in the spanwise and streamwise directions of the flow; and (c) direct numerical integration of the resulting partial differential equations. All of the stability data that are known for such flows can be reproduced. Also, the optimal initial condition can be determined in a straight forward manner and such optimal conditions clearly reflect transient growth data that is easily determined by a rational choice of a basis for the initial conditions. Although there can be significant transient growth for subcritical values of the Reynolds number using this approach it does not appear possible that arbitrary initial conditions will lead to the exceptionally large transient amplitudes that have been determined by optimization of normal modes. The approach is general and can be applied to other classes of problems where only a finite discrete spectrum exists, such as the boundary layer for example.

Criminale, W. O.↗

Structural/control system optimization with variable actuator masses

A method is presented to integrate the design space for structural/control system optimization problems in the case of linear state feedback control. Nonstructural lumped masses and control system design variables as well as structural sizing variables are all treated equally as independent design variables in the optimization process. Structural and control design variable linking schemes are used in order to avoid a prohibitively large increase in the total number of independent design variables. When actuator masses are treated as nonstructural lumped mass design variables, special consideration is given to the relation between the transient peak responses and the required actuator masses which is formulated as a behavior constraint form. The original nonlinear mathematical programming problem based on a finite element formulation and linear state feedback is replaced by a sequence of explicit approximate problems exploiting various approximation concepts such as design variable linkings, temporary constraint deletion and first order Taylor series expansion of nonlinear behavior constraints in terms of intermediate design variables. Examples which involve a variety of dynamic behavior constraints (including constraints on closed-loop eigenvalues, peak transient displacements, peak actuator forces, and relations between the peak responses and the actuator masses) are effectively solved by using the method presented.

Jin, Ik M.↗

Supersonic wind tunnel optimization

Attention is given to major problems arising in the course of definition studies for an optimum supersonic wind tunnel facility that incorporates advanced technologies. The issues identified encompass (1) large amplitude stream disturbances due to acoustic radiation from nozzle wall boundary layers; (2) flow field and model shape distortion caused by sting supports and their installation; (3) lack of sufficient three-dimensional flow visualizaton and diagnostic capabilities that allow ready identification of 'Reynolds number effects'; (4) lack of aerodynamic/propulsion facilities for investigation of the requisite range of vehicle attitudes and transient propulsion behavior; and (5) inordinately large energy usage, primarily due to diffuser inefficiency. Remedies are offered for these problems.

Bushnell, D. M.↗

Special Test Methods for Batteries

Various methods are described for measuring heat generation in primary and secondary batteries as well as the specific heat of batteries and cell thermal conductance. Problems associated with determining heat generation in large batteries are examined. Special attention is given to monitoring temperature gradients in nickel cadmium cells, the use of auxiliary electrodes for conducting tests on battery charge control, evaluating the linear sweep of current from charge to discharge, and determining zero current voltage. The fast transient behavior of batteries in the microsecond range, and the electrical conductance of nickel sinters in the thickness direction are also considered. Mechanical problems experienced in the vibration of Ni-Cd batteries and tests to simulate cyclic fatigue of the steel table connecting the plates to the comb are considered. Methods of defining the distribution of forces when cells are compressed during battery packaging are also explored.

Gross, S.↗

Application of digital computer APU modeling techniques to control system design.

Study of the required controls for a H2-O2 auxiliary power unit (APU) technology program for the Space Shuttle. A steady-state system digital computer program was prepared and used to optimize initial system design. Analytical models of each system component were included. The program was used to solve a nineteen-dimensional problem, and then time-dependent differential equations were added to the computer program to simulate transient APU system and control. Some system parameters were considered quasi-steady-state, and others were treated as differential variables. The dynamic control analysis proceeded from initial ideal control modeling (which considered one control function and assumed the others to be ideal), stepwise through the system (adding control functions), until all of the control functions and their interactions were considered. In this way, the adequacy of the final control design over the required wide range of APU operating conditions was established.

Bailey, D. A.↗

Determination of rotor wake induced empennage airloads

A computer program has been developed that predicts the unsteady aerodynamic forces that impinge on the empennage surface due to its interactions with the main rotor wake. The program was utilized to determine the vibration airloads acting on the Black Hawk horizontal stabilizer and on the CH-53A stabilizer under high speed forward flight conditions. The results demonstrate that it is possible to compute the high frequency empennage vibration airloads efficiently by utilizing suitable analysis techniques; a CH-53A case required only 1.3 minutes of computer time. Numerical problems associated with the unsteady Kutta condition have been minimized. Even though the analysis is of transient nature, the results for harmonic cases show a very fast convergence to periodicity. In general, good correlation was shown between the predicted vibratory airloads and the measured stabilizer airloads for a CH-53A helicopter operating at high speed flight conditions.

Gangwani, S. T.↗

An adaptive finite element scheme for the Euler and Navier-Stokes equations

The implementation of two explicit finite element schemes for the solution of the compressible Euler and Navier-Stokes equations is presented. The schemes can be employed with general unstructured triangular meshes in two dimensions. Either scheme can therefore be used as the basic solver in a solution adaptive mesh procedure in a direct manner. The particular adaptive approach which is advocated here is intended for the solution of steady state problems only and involves an adaptive regeneration of the grid at prescribed stages during the false transient. The grid regeneration is accomplished by a mesh generator which has the capability of generating triangular grids over computational domains of arbitrary shape. The procedure is illustrated by solving transonic flows over multi-airfoil configurations and high speed flows, involving shock interactions, past circular cylinders.

Morgan, K.↗

Application of optical distributed sensing and computation to control of large space structures

A real time holographic sensing technique is introduced and its advantages are investigated from the filtering and control point of view. A feature of holographic sensing is its capability to make distributed measurements of the position and velocity of moving objects, such as a vibrating flexible space structure. This work is based upon the distributed parameter models of linear time invariant systems, particularly including the linear oscillator equations describing the vibration of large flexible space structures. The general conclusion is that application of optical distributed sensors bring gains in the situation where Kalman filtering is necessary for state estimation. In this case, both steady state and transient filtering error covariance become smaller. This in turn results in smaller cost in the LQG problem.

Balakrishnan, A. V.↗

An extended laser flash technique for thermal diffusivity measurement of high-temperature materials

Knowledge of thermal diffusivity data for high-temperature materials (solids and liquids) is very important in analyzing a number of processes, among them solidification, crystal growth, and welding. However, reliable thermal diffusivity versus temperature data, particularly those for high-temperature liquids, are still far from complete. The main measurement difficulties are due to the presence of convection and the requirement for a container. Fortunately, the availability of levitation techniques has made it possible to solve the containment problem. Based on the feasibility of the levitation technology, a new laser flash technique which is applicable to both levitated liquid and solid samples is being developed. At this point, the analysis for solid samples is near completion and highlights of the technique are presented here. The levitated solid sample which is assumed to be a sphere is subjected to a very short burst of high power radiant energy. The temperature of the irradiated surface area is elevated and a transient heat transfer process takes place within the sample. This containerless process is a two-dimensional unsteady heat conduction problem. Due to the nonlinearity of the radiative plus convective boundary condition, an analytic solution cannot be obtained. Two options are available at this point. Firstly, the radiation boundary condition can be linearized, which then accommodates a closed-form analytic solution. Comparison of the analytic curves for the temperature rise at different points to the experimentally-measured values will then provide the thermal diffusivity values. Secondly, one may set up an inverse conduction problem whereby experimentally obtained surface temperature history is used as the boundary conditions. The thermal diffusivity can then be elevated by minimizing the difference between the real heat flux boundary condition (radiation plus convection) and the measurements. Status of an experimental study directed at measuring the thermal diffusivity of high-temperature solid samples of pure Nickel and Inconel 718 superalloys are presented. Preliminary measurements showing surface temperature histories are discussed.

Shen, F.↗

Investigation of Rim Cracking in Turbine Wheels with Welded Blades

Rim cracking in turbine wheels with welded blades was evaluated. The problem is explained on the basis of the occurrence of plastic flow in the rim during transient starting conditions when thermal compressive stresses resulting from high-temperature gradients exceed the proportional elastic limit of the material.

Millenson, M. B.↗

An Analytical Solution for Transient Thermal Response of an Insulated Structure

An analytical solution was derived for the transient response of an insulated aerospace vehicle structure subjected to a simplified heat pulse. This simplified problem approximates the thermal response of a thermal protection system of an atmospheric entry vehicle. The exact analytical solution is solely a function of two non-dimensional parameters. A simpler function of these two parameters was developed to approximate the maximum structural temperature over a wide range of parameter values. Techniques were developed to choose constant, effective properties to represent the relevant temperature and pressure-dependent properties for the insulator and structure. A technique was also developed to map a time-varying surface temperature history to an equivalent square heat pulse. Using these techniques, the maximum structural temperature rise was calculated using the analytical solutions and shown to typically agree with finite element simulations within 10 to 20 percent over the relevant range of parameters studied.

Blosser, Max L.↗

SPAR thermal analysis processors reference manual, system level 16. Volume 1: Program executive. Volume 2: Theory. Volume 3: Demonstration problems. Volume 4: Experimental thermal element capability. Volume 5: Programmer reference

User instructions are given for performing linear and nonlinear steady state and transient thermal analyses with SPAR thermal analysis processors TGEO, SSTA, and TRTA. It is assumed that the user is familiar with basic SPAR operations and basic heat transfer theory.

Marlowe, M. B.↗

Development and Characterization of a Rate-Dependent Three-Dimensional Macroscopic Plasticity Model Suitable for Use in Composite Impact Problems

Several key capabilities have been identified by the aerospace community as lacking in the material/models for composite materials currently available within commercial transient dynamic finite element codes such as LS-DYNA. Some of the specific desired features that have been identified include the incorporation of both plasticity and damage within the material model, the capability of using the material model to analyze the response of both three-dimensional solid elements and two dimensional shell elements, and the ability to simulate the response of composites composed with a variety of composite architectures, including laminates, weaves and braids. In addition, a need has been expressed to have a material model that utilizes tabulated experimentally based input to define the evolution of plasticity and damage as opposed to utilizing discrete input parameters (such as modulus and strength) and analytical functions based on curve fitting. To begin to address these needs, an orthotropic macroscopic plasticity based model suitable for implementation within LS-DYNA has been developed. Specifically, the Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The coefficients in the yield function are determined based on tabulated stress-strain curves in the various normal and shear directions, along with selected off-axis curves. Incorporating rate dependence into the yield function is achieved by using a series of tabluated input curves, each at a different constant strain rate. The non-associative flow-rule is used to compute the evolution of the effective plastic strain. Systematic procedures have been developed to determine the values of the various coefficients in the yield function and the flow rule based on the tabulated input data. An algorithm based on the radial return method has been developed to facilitate the numerical implementation of the material model. The presented paper will present in detail the development of the orthotropic plasticity model and the procedures used to obtain the required material parameters. Methods in which a combination of actual testing and selective numerical testing can be combined to yield the appropriate input data for the model will be described. A specific laminated polymer matrix composite will be examined to demonstrate the application of the model.

Impact↗

Optimal symmetric flight with an intermediate vehicle model

Optimal flight in the vertical plane with a vehicle model intermediate in complexity between the point-mass and energy models is studied. Flight-path angle takes on the role of a control variable. Range-open problems feature subarcs of vertical flight and singular subarcs. The class of altitude-speed-range-time optimization problems with fuel expenditure unspecified is investigated and some interesting phenomena uncovered. The maximum-lift-to-drag glide appears as part of the family, final-time-open, with appropriate initial and terminal transient exceeding level-flight drag, some members exhibiting oscillations. Oscillatory paths generally fail the Jacobi test for durations exceeding a period and furnish a minimum only for short-duration problems.

Menon, P. K. A.↗

Structural dynamic analysis on a parallel computer - The finite element machine

Recent and prospective advances in parallel multiple introduction multiple data (MIMD) computers offer significant improvements in the range of structural problems that can be solved, as well as the speed of solutions. These improvements have their basis in the effective selection and implementation of algorithms exploiting parallel computation. Attention is presently given to the solution of the transient response calculations of an experimental MIMD computer designated the 'Finite Element Machine', including its algorithm and its results for representative one- and two-dimensional dynamic response test problems. A factor of 6.5 is noted for the computational speed improvement.

Storaasli, O.↗

A three-dimensional transport model for the middle atmosphere

In this paper we describe fundamental properties of an 'off-line' three-dimensional transport model, that is, a model which uses prescribed rather than predicted winds. The model is currently used primarily for studying problems of the middle atmosphere because we have not (yet) incorporated a formulation for the convective transport of trace species, a prerequisite for many tropospheric problems. The off-line model is simpler and less expensive than a model which predicts the wind and mass evolution (an 'on-line' model), but it is more complex than the two-dimensional (2-D) zonally averaged transport models often used in the study of chemistry and transport in the middle atmosphere. It thus serves as a model of intermediate complexity and can fill a useful niche for the study of transport and chemistry. We compare simulations of four tracers, released in the lower stratosphere, in both the on- and off-line models to document the difference resulting from differences in modeling the same problem with this intermediate model. These differences identify the price to be paid in going to a cheaper and simpler calculation. The off-line model transports a tracer in three dimensions. For this reason, it requires fewer approximations than 2-D transport model, which must parameterize the effects of mixing by transient and zonally asymmetric wind features. We compare simulations of the off-line model with simulations of a 2-D model for two problems. First, we compare 2-D and three-dimensional (3-D) models by simulating the emission of an NO(x)-like tracer by a fleet of high-speed aircraft. The off-line model is then used to simulate the transport of C-14 and to contrast its simulation properties to that of the host of 2-D models which participated in an identical simulation in a recent NASA model intercomparison. The off-line model is shown to be somewhat sensitive to the sampling strategy for off-line winds. Simulations with daily averaged winds are in very good qualitative agreement but are less diffusive than when driven with instantaneous winds sampled at half-hour intervals. Simulations with the off-line and 2-D models are quite similar in the middle and upper stratosphere but behave quite differently in the lower stratosphere, where the 3-D model has a substantially more vigorous circulation. The off-line model is quite realistic in its simulation of C-14. While there are still systematic differences between the 3-D calculation and the observations, the differences seem to be substantially reduced when compared with the body of 2-D simulations documented in the above mentioned NASA intercomparison, particularly at 31 deg N.

Rasch, Philip J.↗