On the stability of Poiseuille flow in a pipe
Poiseuille pipe flow stability from finite difference equations approximation to nonlinear axisymmetric Navier-Stokes equations under stream function perturbation
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.
Poiseuille pipe flow stability from finite difference equations approximation to nonlinear axisymmetric Navier-Stokes equations under stream function perturbation
The optical response of metal-barrier-metal (MBM) junctions is discussed, with emphasis on the nonlinearity of response associated with reversal of polarity as optical intensity is increased. This nonlinear behavior can be used to measure correlation functions of optical pulses or to generate ultrashort electrical pulses. Expressions for the current flowing through the insulating barrier (junction), and for the rectification current, are given. Pulse signal and polarity of the electrical pulse evolves as the junction bias is varied about the null-signal bias. The apparatus used to obtain rectification pulse signals is described; the calculated evolution of the rectified pulse shape is found to correspond to that of pulse signals obtained experimentally. The generation of pulses six orders of magnitude smaller than a previously-obtained instrument-limited 1.6 nanosecond pulse has been observed.
Wind tunnel measurements show that the wake of a rotor, except at near-hovering speeds, is not like that of a propeller. The wake is more like that of a wing except that, because of the slow speeds, the wake velocities may be much greater. The helicopter can produce a wake hazard to following light aircraft that is disproportionately great compared to an equivalent fixed-wing aircraft. This hazard should be recognized by both pilots and airport controllers when operating in congested areas. Even simple momentum theory shows that, in autorotation and partial-power descent, the required power is a complex function of both airspeed and descent angle. The nonlinear characteristic, together with an almost total lack of usable instrumentation at low airspeeds, has led to numerous power-settling accidents. The same theory shows that there is a minimum forward speed at which a rotor can autorotate. Neglect of, or inadequate appraisal of this minimum speed has also led to numerous accidents. Ground effect and the problems it creates is discussed.
An operational N-ray/X-ray computer analysis system for nondestructive testing problems was realized by developing: (1) more efficient data transfer routines for use in existing image processing systems; (2) the capability to store and analyze multiple N-ray/X-ray images; (3) an efficient 2-D Fast Fourier Transform; and (4) transfer functions compensation capability to correct for distortions and nonlinearities of N-rays/X-rays. Applications of the technique are summarized for the following: thermal protection system tile with hydraulic fluid; thermal protection system tile with corrosion; thermal protection system tile debonding; stainless steel tubing corrosion; defective electrical devices; fast film restoration and deblur; corrosion tracking over time; reduction of source vignetting, and identification of materials (0 rings, Teflon parts, soft goods) via table look up.
The structural behavior of an elastic beam-column placed with a gap between two nonlinearity elastic layers each resting on a rigid foundation was examined. The beam-column was laterally supported at both ends and subjected to a uniform transverse load and axial compression. Its slenderness was such that the axial compressive force exceeds the amount that would be necessary to buckle it as a simple supported column. The elastic layers were represented by an elastic foundation with a strongly nonlinear specific reaction taken as a rapidly increasing function of the layer compression. The analytical model developed simulated the entire pattern of the deflection and stress state including layer and end support reactions, under gradually increasing axial force.
The experimental determination of the subsonic roll derivatives of the Standard Dynamics Model, which is representative of a current fighter aircraft configuration, is described. The direct, cross and cross-coupling derivatives are presented for angles of attack up to 41 deg and sideslip angles in the range from -5 deg to 5 deg, as functions of oscillation frequency. The derivatives exhibited significant nonlinear trends at high incidences and were found to be extremely sensitive to sideslip angle at angles of attack near 36 deg. The roll damping and dynamic cross derivatives were highly frequency dependent at angles of attack above 30 deg. The highest values measured for the dynamic cross and cross-coupling derivatives were comparable in magnitude with the maximum roll damping. The effects of oscillation amplitude and Mach number were also investigated, and the direct derivatives were correlated with data from another facility.
A scheme is developed for solving constrained optimization problems in which the objective function and the constraint function are dependent on the solution of the nonlinear flow equations. The scheme updates the design parameter iterative solutions and the flow variable iterative solutions simultaneously. It is applied to an advanced propeller design problem with the Euler equations used as the flow governing equations. The scheme's accuracy, efficiency and sensitivity to the computational parameters are tested.
A scheme is developed for solving constrained optimization problems in which the objective function and the constraint function are dependent on the solution of the nonlinear flow equations. The scheme updates the design parameter iterative solutions and the flow variable iterative solutions simultaneously. It is applied to an advanced propeller design problem with the Euler equations used as the flow governing equations. The scheme's accuracy, efficiency and sensitivity to the computational parameters are tested.
This user's manual is presented for an aerodynamic optimization program that updates flow variables and design parameters simultaneously. The program was developed for solving constrained optimization problems in which the objective function and the constraint function are dependent on the solution of the nonlinear flow equations. The program was tested by applying it to the problem of optimizing propeller designs. Some reference to this particular application is therefore made in the manual. However, the optimization scheme is suitable for application to general aerodynamic design problems. A description of the approach used in the optimization scheme is first presented, followed by a description of the use of the program.
The three dimensional marginal separation of a boundary layer along a line of symmetry is considered. The key equation governing the displacement function is derived, and found to be a nonlinear integral equation in two space variables. This is solved iteratively using a pseudo-spectral approach, based partly in double Fourier space, and partly in physical space. Qualitatively, the results are similar to previously reported two dimensional results (which are also computed to test the accuracy of the numerical scheme); however quantitatively the three dimensional results are much different.
A parallel algorithm has been developed for rapidly solving trajectory optimization problems. The goal of the work has been to develop an algorithm that is suitable to do real-time, on-line optimal guidance through repeated solution of a trajectory optimization problem. The algorithm has been developed on an INTEL iPSC/860 message passing parallel processor. It uses a zero-order-hold discretization of a continuous-time problem and solves the resulting nonlinear programming problem using a custom-designed augmented Lagrangian nonlinear programming algorithm. The algorithm achieves parallelism of function, derivative, and search direction calculations through the principle of domain decomposition applied along the time axis. It has been encoded and tested on 3 example problems, the Goddard problem, the acceleration-limited, planar minimum-time to the origin problem, and a National Aerospace Plane minimum-fuel ascent guidance problem. Execution times as fast as 118 sec of wall clock time have been achieved for a 128-stage Goddard problem solved on 32 processors. A 32-stage minimum-time problem has been solved in 151 sec on 32 processors. A 32-stage National Aerospace Plane problem required 2 hours when solved on 32 processors. A speed-up factor of 7.2 has been achieved by using 32-nodes instead of 1-node to solve a 64-stage Goddard problem.
An examination is conducted of the steady and unsteady behavior of isothermal and nonisothermal liquid bridge systems under both normal and low-gravity conditions, recasting the full nonlinear equations in terms of a stream-function vorticity formulation in conjunction a nonorthogonal coordinate transformation. For the cases examined, the shape of the bridge is more sensitive to typical spacecraft accelerations than the buoyancy-driven flow. The effect of thermocapillary flow on the surface shape is found to be small, for the range of capillary and Reynolds numbers considered.
The design optimization of air-breathing propulsion engine concepts has been accomplished by soft-coupling the NASA Engine Performance Program (NEPP) analyzer with the NASA Lewis multidisciplinary optimization tool COMETBOARDS. Engine problems, with their associated design variables and constraints, were cast as nonlinear optimization problems with thrust as the merit function. Because of the large number of mission points in the flight envelope, the diversity of constraint types, and the overall distortion of the design space; the most reliable optimization algorithm available in COMETBOARDS, when used by itself, could not produce satisfactory, feasible, optimum solutions. However, COMETBOARDS' unique features-which include a cascade strategy, variable and constraint formulations, and scaling devised especially for difficult multidisciplinary applications-successfully optimized the performance of subsonic and supersonic engine concepts. Even when started from different design points, the combined COMETBOARDS and NEPP results converged to the same global optimum solution. This reliable and robust design tool eliminates manual intervention in the design of air-breathing propulsion engines and eases the cycle analysis procedures. It is also much easier to use than other codes, which is an added benefit. This paper describes COMETBOARDS and its cascade strategy and illustrates the capabilities of the combined design tool through the optimization of a high-bypass- turbofan wave-rotor-topped subsonic engine and a mixed-flow-turbofan supersonic engine.
This paper presents the results of applying two different types of neural networks in two different approaches to the sensor validation problem. The first approach uses a functional approximation neural network as part of a nonlinear observer in a model-based approach to analytical redundancy. The second approach uses an auto-associative neural network to perform nonlinear principal component analysis on a set of redundant sensors to provide an estimate for a single failed sensor. The approaches are demonstrated using a nonlinear simulation of a turbofan engine. The fault detection and sensor estimation results are presented and the training of the auto-associative neural network to provide sensor estimates is discussed.
Experimental aerodynamic researchers require estimated precision and bias uncertainties of measured physical quantities, typically at 95 percent confidence levels. Uncertainties of final computed aerodynamic parameters are obtained by propagation of individual measurement uncertainties through the defining functional expressions. In this paper, rigorous mathematical techniques are extended to determine precision and bias uncertainties of any instrument-sensor system. Through this analysis, instrument uncertainties determined through calibration are now expressed as functions of the corresponding measurement for linear and nonlinear univariate and multivariate processes. Treatment of correlated measurement precision error is developed. During laboratory calibration, calibration standard uncertainties are assumed to be an order of magnitude less than those of the instrument being calibrated. Often calibration standards do not satisfy this assumption. This paper applies rigorous statistical methods for inclusion of calibration standard uncertainty and covariance due to the order of their application. The effects of mathematical modeling error on calibration bias uncertainty are quantified. The effects of experimental design on uncertainty are analyzed. The importance of replication is emphasized, techniques for estimation of both bias and precision uncertainties using replication are developed. Statistical tests for stationarity of calibration parameters over time are obtained.
In recent years, there has been a growing interest in the development of passive optical power limiters for the protection of the human eye and solid-state sensors from damage caused by energetic light pulses and also for other switching applications. One of the key issues involved is the search for appropriate materials that show effective reverse saturable absorption. Phthalocyanines seem to be good candidates for such applications because of their higher third order nonlinearity and the unique electronic absorption characteristics. A series of 1,4,8,11,15, 18,22,25-octa-alkoxy metallophthalocyanines containing various central metal atoms such as zinc, copper, palladium, cobalt and nickel were characterized for their third order nonlinearity and for their nonlinear absorptive properties to evaluate their suitability to function as reverse saturable absorbers.
The Monte Carlo simulation method in conjunction with the finite element large deflection modal formulation are used to estimate fatigue life of aircraft panels subjected to stationary Gaussian band-limited white-noise excitations. Ten loading cases varying from 106 dB to 160 dB OASPL with bandwidth 1024 Hz are considered. For each load case, response statistics are obtained from an ensemble of 10 response time histories. The finite element nonlinear modal procedure yields time histories, probability density functions (PDF), power spectral densities and higher statistical moments of the maximum deflection and stress/strain. The method of moments of PSD with Dirlik's approach is employed to estimate the panel fatigue life.
We present simulations of the interannual variability of martian global dust storms (GDSs) with a simplified low-order model (LOM) of the general circulation. The simplified model allows one to conduct computationally fast long-term simulations of the martian climate system. The LOM is constructed by Galerkin projection of a 2D (zonally averaged) general circulation model (GCM) onto a truncated set of basis functions. The resulting LOM consists of 12 coupled nonlinear ordinary differential equations describing atmospheric dynamics and dust transport within the Hadley cell. The forcing of the model is described by simplified physics based on Newtonian cooling and Rayleigh friction. The atmosphere and surface are coupled: atmospheric heating depends on the dustiness of the atmosphere, and the surface dust source depends on the strength of the atmospheric winds. Parameters of the model are tuned to fit the output of the NASA AMES GCM and the fit is generally very good. Interannual variability of GDSs is possible in the IBM, but only when stochastic forcing is added to the model. The stochastic forcing could be provided by transient weather systems or some surface process such as redistribution of the sand particles in storm generating zones on the surface. The results are sensitive to the value of the saltation threshold, which hints at a possible feedback between saltation threshold and dust storm activity. According to this hypothesis, erodable material builds up its a result of a local process, whose effect is to lower the saltation threshold until a GDS occurs. The saltation threshold adjusts its value so that dust storms are barely able to occur.