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 steady discharge of a high-voltage corona is analyzed as a singular perturbation problem. The small parameter is the ratio of the length of the ionization region to the total gap length. By this method, current versus voltage characteristics can be calculated analytically.
This lecture is introductory to the subject of unsteady subsonic and supersonic flows. The primary objective is to present fundamental concepts in order to promote an understanding of the relations between the basic physical problems and their mathematical formulation as well as to establish a common foundation for the more detailed presentations of subsequent lectures in this session. Linearized (small-perturbation) potential flow is emphasized, although needs beyond that limit are indicated. The basic equations, concepts, and procedures common to all the methods are reviewed first, followed by the development, discussion, and status of methods for creating two-dimensional incompressible flow, strip theory, subsonic lifting-surface theory, subsonic/supersonic surface-panel methods, and supersonic lifting-surface theory.
Most published work relating to dynamically loaded journal bearings are directed to determining the minimum film thickness from the predicted journal trajectories. These do not give any information about the subsynchronous whirl stability of journal bearing systems since they do not consider the equations of motion. It is, however, necessary to know whether the bearing system operation is stable or not under such an operating condition. The stability characteristics of the system are analyzed. A linearized perturbation theory about the equilibrium point can predict the threshold of stability; however it does not indicate postwhirl orbit detail. The linearized method may indicate that a bearing is unstable for a given operating condition whereas the nonlinear analysis may indicate that it forms a stable limit cycle. For this reason, a nonlinear transient analysis of a rigid rotor supported on oil journal bearings under: (1) a unidirectional constant load, (2) a unidirectional periodic load, and (3) variable rotating load are performed. The hydrodynamic forces are calculated after solving the time-dependent Reynolds equation by a finite difference method with a successive overrelaxation scheme. Using these forces, equations of motion are solved by the fourth-order Runge-Kutta method to predict the transient behavior of the rotor. With the aid of a high-speed digital computer and graphics, the journal trajectories are obtained for several different operating conditions.
An analytical method for studying characteristics of frozen orbits using mean orbital elements is presented. By removing some of the approximations of Brouwer and Kosai (1959) and using the method of integration by parts, an analytic solution is obtained accurate to the second order in the sense of expansion of the perturbing potential of an aspherical approximation. The solution applies to orbits at and near the critical inclination as well as to frozen orbits and, of course, to unfrozen orbits. The analysis reveals that a frozen orbit can be interpreted as an orbit in which movements of the long-periodic part of the rate of change of periapsis cancels the secular part exactly at any time.
On large scales, the higher order moments of the mass distribution, S(sub J) = bar-zeta(sub J)/bar-zeta(sup J-1)(sub 2), e.g., the skewness S(sub 3) and kurtosis S(sub 4), can be predicted using nonlinear perturbation theory. Comparison of these predictions with moments of the observed galaxy distribution probes the bias between galaxies and mass. Applying this method to models with initially Gaussian fluctuations and power spectra P(k) similar to that of galaxies in the Automatic Plate Measuring (APM) survey, we find that the predicted higher order moments S(sub J)(R) are in good agreement with those directly inferred from the APM survey in the absence of bias. We use this result to place limits on the linear and nonlinear bias parameters. Models in which the extra power observed on large scales (with respect to the standard cold dark matter (CDM) model) is produced by scale-dependent bias match the APM higher order amplitudes only if nonlinear bias (rather than nonlinear gravity) generates the observed higher order moments. When normalized to Cosmic Background Explorer Differential Microwave Radiometer (COBE DMR), these models are siginificantly ruled out by the S(sub 3) observations. The cold plus hot dark matter model normalized to COBE can reproduce the APM higher order correlations if one introduces nonlinear bias terms, while the low-density CDM model with a cosmological constant does not require any bias to fit the large-scale amplitudes.
This paper describes the formation of nonlinear structure in flat (zero curvature) Friedmann cosmological models. We consider models with two components: the usual nonrelativistic component that evolves under gravity and eventually forms the large-scale structure of the universe, and a uniform dark matter component that does not clump under gravity, and whose energy density varies with the scale factor a(t) like a(t)(sup -n), where n is a free parameter. Each model is characterized by two parameters: the exponent n and the present density parameter Omega(sub 0) of the nonrelativistic component. The linear perturbation equations are derived and solved for these models, for the three different cases n = 3, n is greater than 3, and n is less than 3. The case n = 3 is relevant to model with massive neutrinos. The presence of the uniform component strongly reduces the growth of the perturbation compared with the Einstein-de Sitter model. We show that the Meszaros effect (suppression of growth at high redshift) holds not only for n = 4, radiation-dominated models, but for all models with n is greater than 3. This essentially rules out any such model. For the case n is less than 3, we numerically integrate the perturbation equations from the big bang to the present, for 620 different models with various values of Omega(sub 0) and n. Using these solutions, we show that the function f(Omega(sub 0), n) = (a/delta(sub +))d(delta)(sub +)/da, which enters in the relationship between the present density contrast delta(sub 0) and peculiar velocity field u(sub 0) is essentially independent of n. We derive approximate solutions for the second-order perturbation equations. These second-order solutions are tested against the exact solutions and the Zel'dovich approximation for spherically symmetric perturbations in the marginally nonlinear regime (the absolute value of delta is less than or approximately 1). The second-order and Zel'dovich solutions have comparable accuracy, significantly higher than the accuracy of the linear solutions. We then investigate the dependence of the delta(sub 0) - u(sub 0) relationship upon the value of n in the nonlinear regime using the second-order solutions for marginally nonlinear, general perturbations, and the exact solutions for strongly nonlinear, spherically symmetric perturbations. In both cases, we find that the delta(sub 0) - u(sub 0) relationship remains independent of n. We speculate that this result extends to strongly nonlinear, general perturbations as well. This eliminates any hope to determine the presence of the uniform component or the value of n using dynamical methods. Finally, we compute the nonlinear evolution of the skewness of the distribution of values of delta, assuming Gaussian initial conditions. We find that the skewness is not only independent of n, but also of Omega(sub 0). Thus the skewness cannot be used to discriminate among various models with Gaussian initial conditions. However, it can be used for testing the Gaussianity of the initial conditions themselves. We conclude that the uniform component leaves no observable signature in the present large-scale structure of the universe. To determine its presence and nature, we must investigate the relationship between the past and present universe, using redshift-dependent tests.
A new non-perturbing technique for feedback control of plasma instabilities has been developed in the Columbia Linear Machine (CLM). The feedback control scheme involves the injection of a feedback modulated ion beam as a remote suppressor. The ion beam was obtained from a compact ion beam source which was developed for this purpose. A Langmuir probe was used as the feedback sensor. The feedback controller consisted of a phase-shifter and amplifiers. This technique was demonstrated by stabilizing various plasma instabilities to the background noise level, like the trapped particle instability, the ExB instability and the ion-temperature-gradient (ITG) driven instability. An important feature of this scheme is that the injected ion beam is non-perturbing to the plasma equilibrium parameters. The robustness of this feedback stabilization scheme was also investigated. The principal result is that the scheme is fairly robust, tolerating about 100% variation about the nominal parameter values. Next, this scheme is extended to the unsolved general problem of controlling multimode plasma instabilities simultaneously with a single sensor-suppressor pair. A single sensor-suppressor pair of feedback probes is desirable to reduce the perturbation caused by the probes. Two plasma instabilities the ExB and the ITG modes, were simultaneously stabilized. A simple 'state' feedback type method was used where more state information was generated from the single sensor Langmuir probe by appropriate signal processing, in this case, by differentiation. This proof-of-principle experiment demonstrated for the first time that by designing a more sophisticated electronic feedback controller, many plasma instabilities may be simultaneously controlled. Simple theoretical models showed generally good agreement with the feedback experimental results. On a parallel research front, a better understanding of the saturated state of a plasma instability was sought partly with the help of feedback. A plasma instability is usually observed in its saturated state and appears as a single feature in the frequency spectrum with a single azimuthal and parallel wavenumbers. The physics of the non-zero spectral width was investigated in detail because the finite spectral width can cause "turbulent" transport. One aspect of the "turbulence" was investigated by obtaining the scaling of the linear growth rate of the instabilities with the fluctuation levels. The linear growth rates were measured with the established gated feedback technique. The research showed that the ExB instability evolves into a quasi-coherent state when the fluctuation level is high. The coherent aspects were studied with a bispectral analysis. Moreover, the single spectral feature was discovered to be actually composed of a few radial harmonics. The radial harmonics play a role in the nonlinear saturation of the instability via three-wave coupling.
In order to address the problem of sources and sinks of CO2 measurements are needed on a global scale. Clearly a satellite is a promising approach to meeting this requirement. Unfortunately, most methods for making a CO2 measurement from space involve the whole column. Since sources and sinks at the surface represent a small perturbation to the total column one is faced with the need to measure the column with a precision better than 1%. No species has ever been measured from space at this level. We have developed over the last 3 years a small instrument based upon a Fabry-Perot interferometer that is very sensitive to atmospheric CO2 and has a high signal to noise ratio. We have tested this instrument in a ground based configuration and from aircraft platforms simulating operation from a satellite. We will present results from these tests and discuss ways that this promising new instrument could be used to improve our understanding of the global carbon budget.
The Electron Electrostatic Analyzer (EEA) is part of the HERMES instrument suite on Lunar Gateway and is designed to measure electrons in the energy range from 1 eV to 18 keV. Electrons in this range can be deflected or blocked from entering a detector by electric fields from a charged spacecraft platform, so many heliophysics missions (MMS, Cluster, etc.) use active potential control to minimize stray electric fields which would corrupt low-energy plasma measurements. Gateway, however, does not have such capabilities so the surface charging effects will influence the electron data. Further, some elements of Gateway partially block the instrument’s field of view to space. The effects of the Gateway platform on the low-energy electron population must be understood to maximize the accuracy of the EEA data. To investigate the electrostatic environment near Gateway and its effects on the EEA measurements, we use Nascap-2K to model the expected charging of Gateway under representative solar wind, magnetosheathic, and magnetospheric conditions. Next, we use the potential and electric field outputs from Nascap-2K along with in-house particle tracing codes to investigate how low-energy electrons are perturbed, blocked, or deflected by the spacecraft. We discuss how these perturbations impact the plasma moments computed from the EEA data and possible approaches for correcting the charging effects including application of advanced methods such as neural networks. In addition to electrons from the space environment being influenced by the electrostatic environment around Gateway, photoelectrons generated on Gateway surfaces will also be measured by EEA. We are investigating how knowledge of the potentials of various Gateway surfaces can be extracted from these measured photoelectrons using advanced numerical modeling and analysis. Such information would be valuable for validating charging models and planning mission operations around the platform. Finally, we discuss how the specific workflow developed for Gateway can be extended to a generalized process that is applicable to other future missions. As the space industry prepares for human exploration beyond low-Earth orbit, the need to obtain accurate space weather data from rideshare science payloads is critical. Our process for removing platform effects from the science data serves as a pathfinder for leveraging advanced simulations to maximize science returns and ensure the safety and success of future crewed missions.
A novel method is presented for the investigation of Fabry-Perot etalon defects, incorporating the effects of microsurface imperfection, bowing, and departure from parallelism, as well as the effect of finite aperture on the etalon, into the perturbation function. The tolerance imperfections of reflective layers can be incorporated into the transfer coefficients of a noncorrugated thin-film system, and any defect can be investigated either separately or in combination with others through the use of the appropriate analytical or numerical simulation of the corrugated function. In the context of numerical computation, the present method permits the simulation and evaluation of different specific situations. The expressions derived for the transfer coefficients of a corrugated, equidistant thin film system may be applicable to other branches of applied optics.
A method for computing the flow around a hovering helicopter rotor is presented. The solution procedure is divided into two parts. The first part is a finite volume solution of the Euler equations for the near field flow around a rotor blade. The second part is a free wake approach for determining the wake geometry and induced velocities. These two parts are solved in a coupled fashion. The rolled-up vortex wake of the rotor is included in the Euler computational domain through a nonlinear perturbation technique. This eliminates the numerical diffusion of the wake vorticity due to truncation error and artificial viscosity without the need for excessive grid resolution near the vortex core. The method is used to compute both hovering rotor and wing/vortex interaction flows, and the results are compared to experiment.
Small perturbation theory is applied to compute the deflection of the wind blowing across land that has an irregular topography. As an illustration, the method is applied first to the flow around a single hill of Gaussian profile. Then calculations are made for the irregular topography on the east side of Cleveland where the elevation changes by several hundred feet. It was found that the topography produced small wind deflections that would not be of practical importance in air pollution dispersion studies. The calculations were for a neutrally stable atmosphere. Although they are not investigated here, other factors such as thermal stratification of the atmosphere, diurnal variations, and convection currents resulting from the proximity of Lake Erie and the city heat island effect are expected to be more significant than the influence of topography.
Derivation of the sensitivity equations that yield the sensitivity derivatives directly, which avoids the costly and inaccurate perturb-and-reoptimize approach, is discussed and solvability of the equations is examined. The equations apply to optimum solutions obtained by direct search methods as well as those generated by procedures of the sequential unconstrained minimization technique class. Applications are discussed for the use of the sensitivity derivatives in extrapolation of the optimal objective function and design variable values for incremented parameters, optimization with multiple objectives, and decomposition of large optimization problems.
Measurement of the gradient of the gravitational acceleration from a satellite platform is likely to provide the next improvement in knowledge of the earth's gravity field after the upcoming Geopotential Research Mission (GRM). Observations from the subsatellite of a Tethered Satellite System (TSS) would increase sensitivity and resolution due to the low altitude possible. However, the TSS is a dynamically 'noisy' system and would be perturbed by atmospheric drag fluctuations. The dynamic noise is being modeled in order to evaluate the feasibility of TSS gradiometry and to design methods of abating the error caused by this noise. The demonstration flights of the TSS will provide an opportunity to directly observe the dynamical environment and refine modeling techniques.
An analytical method of obtaining the natural frequencies and mode shapes of clamped, rectangular plates having a small initial curvature is presented. Specifically, the singular perturbation technique is used to reduce the fourth-order plate vibration problem to the simpler membrane problem with modified boundary conditions that account for the bending effects. The eigenfrequencies for plates with inverse aspect ratios varying between 0.1 and 1.0 and for the dimensionless normal prestress between 0.1 and 1.0 have been presented for values of epsilon, the normalized bending rigidity, ranging between 0.0010 and 0.2500. It is established that a small initial curvature has no effect on the frequency of vibration of the plate. However, its effect is manifested in the eigenmodes.
A spectral method which employs trigonometric functions and Chebyshev polynomials is used to compute the steady, incompressible laminar flow past a circular cylinder. Linear stability methods are used to formulate a pair of decoupled generalized eigenvalue problems for the growth of symmetric and asymmetric (about the dividing streamline) perturbations. It is shown that, while the symmetric disturbances are stable, the asymmetric perturbations become unstable at a Reynolds number about 40 with a Strouhal number about 0.12. The critical conditions are found to depend on the size of the computational domain in a manner similar to that observed in the laboratory.
The magnetospheric ULF waves observed during the magnetospheric compression event of November 1, 1984, were investigated using the magnetic field and medium-energy ion data obtained by the AMPTE Charge Composition Explorer. The electric field of ULF waves was inferred from the particle data, using a method developed on the basis of the gyration acceleration mechanism. Three types of waves were found: (1) a 10-15-min wave with perturbations in the magnetic field intensity and in the flux of ions; (2) a 3-5-min standing Alfven wave with perturbations in the azimuthal magnetic field component and in the radial component of the electric field; and (3) a 2-5-min irregular disturbance near the magnetopause, which involves all components of the magnetic field and the intensity of the ion flux. The origins of these waves are discussed.
The evolution of a small perturbation introduced periodically into a boundary layer with a moderately strong adverse pressure gradient (APG) was investigated experimentally using the hot-wire method. The magnitude of the disturbance is found to decay with the streamwise distance in the APG before undergoing rapid growth and triggering transition in a very repetitive manner. The evidence suggests that the laminar layer separates and that the transition mechanism is an inviscid shear-layer type of instability. Contours of spanwise vorticity in the plane of the disturbance show the formation of roll-ups in the initial stages. Subsequent interaction with the wall leads to the formation of large-scale vortex loops which merge into a single vortex loop further downstream.