Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “perturbation analysis”

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 1,153 records · Page 64

The Influence of Particle Charge on Heterogeneous Reaction Rate Coefficients

The effects of particle charge on heterogeneous reaction rates are presented. Many atmospheric particles, whether liquid or solid are charged. This surface charge causes a redistribution of charge within a liquid particle and as a consequence a perturbation in the gaseous uptake coefficient. The amount of perturbation is proportional to the external potential and the square of the ratio of debye length in the liquid to the particle radius. Previous modeling has shown how surface charge affects the uptake coefficient of charged aerosols. This effect is now included in the heterogeneous reaction rate of an aerosol ensemble. Extension of this analysis to ice particles will be discussed and examples presented.

Aikin, A. C.↗

Sub-Micrometer Scale Minor Element Mapping in Interplanetary Dust Particles: A Test for Stratospheric Contamination

Combined X-ray microprobe (XRM), energy dispersive x-ray fluorescence using a Transmission Electron Microscope (TEM), and electron microprobe measurements have determined that the average bulk chemical composition of the interplanetary dust particles (IDPs) collected from the Earth s stratosphere is enriched relative to the CI meteorite composition by a factor of 2 to 4 for carbon and for the moderately volatile elements Na, K, P, Mn, Cu, Zn, Ga, Ge, and Se, and enriched to approximately 30 times CI for Br. However, Jessberger et al., who have reported similar bulk enrichments using Proton Induced X-ray Emission (PIXE), attribute the enrichments to contamination by meteor-derived atmospheric aerosols during the several weeks these IDPs reside in the Earth s atmosphere prior to collection. Using scanning Auger spectroscopy, a very sensitive surface analysis technique, Mackinnon and Mogk have observed S contamination on the surface of IDPs, presumably due to the accretion of sulfate aerosols during stratospheric residence. But the S-rich layer they detected was so thin (approximately 100 angstroms thick) that the total amount of S on the surface was too small to significantly perturb the bulk S-content of a chondritic IDP. Stephan et al. provide support for the contamination hypothesis by reporting the enrichment of Br on the edges of the IDPs using Time-of-Flight Secondary-Ion Mass-Spectrometry (TOFSIMS), but TOF-SIMS is notorious for producing false edge-effects, particularly on irregularly-shaped samples like IDPs. Sutton et al. mapped the spatial distribution of Fe, Ni, Zn, Br, and Sr, at the approximately 2 m scale, in four IDPs using element-specific x-ray fluorescence (XRF) computed microtomography. They found the moderately volatile elements Zn and Br, although spatially inhomogeneous, were not concentrated on the surface of any of the IDPs they examined, suggesting that the Zn and the Br enrichments in the IDPs are not due to contamination during stratospheric residence.

Flynn, G. J.↗

Motion of the Mantle in the Translational Modes of the Earth and Mercury

Slichter modes refer to the translational motion of the inner core with respect to the outer core and the mantle. The polar Slichter mode is the motion of the inner core along the axis of rotation. Busse presented an analysis of the polar mode which yielded an expression for its period. Busse's analysis included the assumption that the mantle was stationary. This approximation is valid for planets with small inner cores, such as the Earth whose inner core is about 1/60 of the total planet mass. On the other hand, many believe that Mercury's core may be enormous. If so, the motion of the mantle should be expected to produce a significant effect. We present a formal framework for including the motion of the mantle in the analysis of the translational motion of the inner core. We analyze the effect of the motion of the mantle on the Slichter modes for a non-rotating planet with an inner core of arbitrary size. We omit the effects of viscosity in the outer core, magnetic effects, and solid tides. Our approach is perturbative and is based on a linearization of Euler's equations for the motion of the fluid and Newton's second law for the motion of the inner core. We find an analytical expression for the period of the Slichter mode. Our result agrees with Busse's in the limiting case of small inner core. We present the unexpected result that even for Mercury the motion of the mantle does not significantly change the period of oscillation.

Grinfeld, Pavel↗

Nonlinear Evolution of Instabilities in a Laminar Separation Bubble at a Hypersonic Mach Number

The development of both convective stationary perturbation as well as global instabilities in the vicinity of a laminar separation bubble above an axisymmetric compression corner in a hypersonic flow is investigated using numerical simulations. The flow configuration of primary interest corresponds to the cone-cylinder-flare model used in experimental measurements in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. For a flare angle of 10 degrees and a unit Reynolds number of 11.5 x 10 6 m -1 , their surface flow visualizations identified the presence of streamwise elongated thermal streaks near the reattachment position and the dominant azimuthal spacing between the streaks was determined to be approximately 10 degrees (i.e., an azimuthal mode number of 36). Previous linear stability analyses predicted that the amplification characteristics of small amplitude, unsteady, convective instabilities within this flow were consistent with the surface pressure fluctuations measured in the experiment. However, their accompanying analysis of global instabilities had found the separation bubble to be weakly unstable at the 10 degrees flare angle, with the most unstable global mode corresponding to a stationary disturbance with an azimuthal wavenumber of m ≈ 5, which was well below the measured wavenumber of m = 36. Besides confirming the presence of the global instability at these flow conditions, the present numerical simulations quantify the details of the stationary equilibrium state associated with the supercritical bifurcation resulting from the nonlinear saturation of the unstable global mode. Although velocity perturbations associated with the saturated global mode are dominated by the fundamental spanwise wavelength associated with the linear global instability, the surface heat flux downstream of reattachment is dominated by m = 36 in agreement with experimental measurements. However, the shorter wavelength peak ( m = 36) in the azimuthal spectrum is not manifested as streaks with an azimuthal spacing of 10 degrees. An additional simulation of the combined evolution of both unsteady convective and global instabilities indicates that the nonlinear interactions between a broadband spectrum of oblique Mack’s first mode instabilities generate stationary vortex modes with a higher dominant wavenumber ( m = 60) as that observed in the experiments. Additional analyses are performed to predict the global mode behavior at other flare angles and how it differs from the convective centrifugal instability due to concave streamline curvature over the compression corner.

Hypersonic flow↗

Nonlinear Evolution of Instability in a Laminar Separation Bubble at Hypersonic Mach Number

The development of both convective stationary perturbation as well as global instabilities in the vicinity of a laminar separation bubble above an axisymmetric compression corner in a hypersonic flow is investigated using numerical simulations. The flow configuration of primary interest corresponds to the cone-cylinder-flare model used in experimental measurements in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. For a flare angle of 10 degrees and a unit Reynolds number of 11.5 x 10 6 m -1 , their surface flow visualizations identified the presence of streamwise elongated thermal streaks near the reattachment position and the dominant azimuthal spacing between the streaks was determined to be approximately 10 degrees (i.e., an azimuthal mode number of 36). Previous linear stability analyses predicted that the amplification characteristics of small amplitude, unsteady, convective instabilities within this flow were consistent with the surface pressure fluctuations measured in the experiment. However, their accompanying analysis of global instabilities had found the separation bubble to be weakly unstable at the 10 degrees flare angle, with the most unstable global mode corresponding to a stationary disturbance with an azimuthal wavenumber of m ≈ 5, which was well below the measured wavenumber of m = 36. Besides confirming the presence of the global instability at these flow conditions, the present numerical simulations quantify the details of the stationary equilibrium state associated with the supercritical bifurcation resulting from the nonlinear saturation of the unstable global mode. Although velocity perturbations associated with the saturated global mode are dominated by the fundamental spanwise wavelength associated with the linear global instability, the surface heat flux downstream of reattachment is dominated by m = 36 in agreement with experimental measurements. However, the shorter wavelength peak ( m = 36) in the azimuthal spectrum is not manifested as streaks with an azimuthal spacing of 10 degrees. An additional simulation of the combined evolution of both unsteady convective and global instabilities indicates that the nonlinear interactions between a broadband spectrum of oblique Mack’s first mode instabilities generate stationary vortex modes with a higher dominant wavenumber ( m = 60) as that observed in the experiments. Additional analyses are performed to predict the global mode behavior at other flare angles and how it differs from the convective centrifugal instability due to concave streamline curvature over the compression corner.

High speed flow↗

Lunar surface solar wind observations at the Apollo 12 and Apollo 15 sites

The solar wind-lunar magnetic field interaction is examined by comparing Apollo 12 and Apollo 15 solar wind spectrometric data. The instrumentation and data analysis methods are described along with the observations at the Apollo 15 site, and the data from the two sites are compared. The results show no noticeable differences between the properties of the upstream solar wind and the plasma observed at the Apollo 15 site (where the local magnetic field is relatively weak) and strong perturbations in the solar wind at the Apollo 12 site (where the field is relatively strong), which include deceleration, deflection, and heating of solar wind protons, focusing or defocusing of the ion flux, and an increased level of plasma parameter fluctuations. These effects are shown to require a charge-separated electric field above the lunar surface and a scale size of about 5 km for the local magnetic field at the Apollo 12 site and to suggest that local lunar magnetic field regions cause lunar limb compression waves which should be more noticeable for low solar wind dynamic pressures.

Clay, D. R.↗

Investigation of the effects of external current systems on the MAGSAT data utilizing grid cell modeling techniques

The feasibility of modeling magnetic fields due to certain electrical currents flowing in the Earth's ionosphere and magnetosphere was investigated. A method was devised to carry out forward modeling of the magnetic perturbations that arise from space currents. The procedure utilizes a linear current element representation of the distributed electrical currents. The finite thickness elements are combined into loops which are in turn combined into cells having their base in the ionosphere. In addition to the extensive field modeling, additional software was developed for the reduction and analysis of the MAGSAT data in terms of the external current effects. Direct comparisons between the models and the MAGSAT data are possible.

Klumpar, D. M.↗

Computational Investigation of Experimental Interaction Impedance Obtained by Perturbation for Helical Traveling-Wave Tube Structures

Conventional methods used to measure the cold-test interaction impedance of helical slow-wave structures involve perturbing a helical circuit with a cylindrical dielectric rod placed on the central axis of the circuit. It has been shown that the difference in resonant frequency or axial phase shift between the perturbed and unperturbed circuits can be related to the interaction impedance. However, because of the complex configuration of the helical circuit, deriving this relationship involves-several approximations. With the advent of accurate three- dimensional helical circuit models, these standard approximations can be fully investigated. This paper addresses the most prominent approximations made in the analysis for measured interaction impedance by Lagerstrom and investigates their accuracy using the three-dimensional simulation code MAFIA. It is shown that a more accurate value of interaction impedance can be obtained by using three-dimensional computational methods rather than performing costly and time consuming experimental cold-test measurements.

Kory, Carol L.↗

A Flight Dynamics Model for a Multi-Actuated Flexible Rocket Vehicle

A comprehensive set of motion equations for a multi-actuated flight vehicle is presented. The dynamics are derived from a vector approach that generalizes the classical linear perturbation equations for flexible launch vehicles into a coupled three-dimensional model. The effects of nozzle and aerosurface inertial coupling, sloshing propellant, and elasticity are incorporated without restrictions on the position, orientation, or number of model elements. The present formulation is well suited to matrix implementation for large-scale linear stability and sensitivity analysis and is also shown to be extensible to nonlinear time-domain simulation through the application of a special form of Lagrange s equations in quasi-coordinates. The model is validated through frequency-domain response comparison with a high-fidelity planar implementation.

Orr, Jeb S.↗

Transients and Dayside Magnetospheric Interactions

Observations at and deep within the Earth’s dayside magnetosphere provide considerable evidence indicating that the solar wind-magnetosphere interaction is often unsteady. Bursty reconnection resulting in the formation of flux transfer events and transient auroral brightenings is common. Fluctuations in solar wind parameters may trigger some bursts of reconnection, while other bursts may result from intrinsic magnetopause instabilities. Large amplitude solitons often travel along the dayside magnetopause. Deeper within the dayside magnetosphere, quasi-periodic transient compressions energize radiation belt electrons and ring current ions. Most of these perturbations occur behind the quasi-parallel bow shock, indicating that they are driven by density structures generated by kinetic processes within the foreshock. This presentation compares and contrasts the amplitudes and occurrence patterns of density and dynamic pressure perturbations intrinsic to the solar wind with those generated within the foreshock, specifies how they move across the magnetosphere, examines recent global MHD simulation results for their interaction with the bow shock, magnetosheath, and magnetopause. It then outlines new paths toward quantifying the significance of bursty reconnection to the overall interaction, via new imaging and analysis techniques.

D. G. Sibeck↗

Stabilization of Hypersonic Boundary Layers by Linear and Nonlinear Optimal Perturbations

The effect of stationary, finite-amplitude, linear and nonlinear optimal perturbations on the modal disturbance growth in a Mach 6 axisymmetric flow over a 7 deg. half-angle cone with 0:126 mm nose radius and 0:305 m length is investigated. The freestream parameters (M = 6, Re(exp 1) = 18 x 10(exp. 6) /m) are selected to match the flow conditions of a previous experiment in the VKI H3 hypersonic tunnel. Plane-marching parabolized stability equations are used in conjunction with a partial-differential equation based planar eigenvalue analysis to characterize the boundary layer instability in the presence of azimuthally periodic streaks. The streaks are observed to stabilize nominally planar Mack mode instabilities, although oblique Mack mode and first-mode disturbances are destabilized. Experimentally measured transition onset in the absence of any streaks correlates with an amplification factor of N = 6 for the planar Mack modes. For high enough streak amplitudes, the transition threshold of N = 6 is not reached by the Mack mode instabilities within the length of the cone; however, subharmonic first-mode instabilities, which are destabilized by the presence of the streaks, do reach N = 6 near the end of the cone. The highest stabilization is observed at streak amplitudes of approximately 20 percent of the freestream velocity. Because the use of initial disturbance profiles based on linear optimal growth theory may yield suboptimal control in the context of nonlinear streaks, the computational predictions are extended to nonlinear optimal growth theory. Results show that by using nonlinearly optimal perturbation leads to slightly enhanced stabilization of plane Mack mode disturbances as well as reduced destabilization of subharmonic first-mode disturbances.

Incompressible boundary layer↗

On the determination and investigation of the terrestrial ionospheric refractive indices using GEOS-3/ATS-6 satellite-to-satellite tracking data

An analysis of the shortening and lengthening of the phase of satellite-to-satellite (SST) data that passed within 40-700 km above the earth surface during its ATS-6 to GEOS-3 to ATS-6 path resulted in refractivity vs height profiles. The SST Doppler data were used directly to adjust the GEOS-3 orbit. Perturbations from the moon, sun, and a 15th-order/degree earth gravity field were included in the orbit solution. This orbit was continued through the occulation period and a model ionosphere was estimated by a least-squares adjustment of the Chapman ionosphere parameters from the SST data residuals. The refractivity profile obtained by this model ionosphere was compared to a refractivity profile obtained by a direct integral inversion of the SST data residuals. Systematic differences between the two methods were caused by orbital errors, which propagated into the solution.

Liu, A. S.↗

Biological extinction in earth history

Virtually all plant and animal species that have ever lived on the earth are extinct. For this reason alone, extinction must play an important role in the evolution of life. The five largest mass extinctions of the past 600 million years are of greatest interest, but there is also a spectrum of smaller events, many of which indicate biological systems in profound stress. Extinction may be episodic at all scales, with relatively long periods of stability alternating with short-lived extinction events. Most extinction episodes are biologically selective, and further analysis of the victims and survivors offers the greatest chance of deducing the proximal causes of extinction. A drop in sea level and climatic change are most frequently invoked to explain mass extinctions, but new theories of collisions with extraterrestrial bodies are gaining favor. Extinction may be constructive in a Darwinian sense or it may only perturb the system by eliminating those organisms that happen to be susceptible to geologically rare stresses.

Raup, D. M.↗

Origins of hydrodynamic forces on centrifugal pump impellers

Hydrodynamic interactions that occur between a centrifugal pump impeller and volute are experimentally and theoretically investigated. The theoretical analysis considers the inability of the blades to perfectly guide the flow through the impeller, and also includes a quasi-one dimensional treatment of the flow in the volute. The disturbance at the impeller discharge and the resulting forces are determined by the theoretical model. The model is then extended to obtain the hydrodynamic force perturbations that are caused by the impeller whirling eccentrically in the volute. Under many operating conditions, these force perturbations were found to be destablizing. Comparisons are made between the theoretical model and the experimental measurements of pressure distributions and radial forces on the impeller. The theoretical model yields fairly accurate predictions of the radial forces caused by the flow through the impeller. However, it was found that the pressure acting on the front shroud of the impeller has a substantial effect on the destablizing hydrodynamic forces.

Adkins, Douglas R.↗

Computational aspects of sensitivity calculations in linear transient structural analysis

The calculation of sensitivities of displacements, velocities, accelerations, and stresses in linear, structural, and transient response problems is studied. Several existing sensitivity calculation methods and two new methods are compared for three example problems. Approximation vectors such as vibration mode shapes are used to reduce the dimensionality of the finite model. To be computationally practical in large-order problems, the overall finite difference methods must use the approximation vectors from the original design in the analyses of the perturbed models. This was found to result in poor convergence of stress sensitivities in several cases. Two semianalytical techniques are developed to overcome this poor convergence. Both new methods result in very good convergence of the stress sensitivities; the computational cost is much less than would result if the vibration modes were recalculated and then used in an overall finite difference method.

Greene, W. H.↗

Analysis of earth albedo effect on sun sensor measurements based on theoretical model and mission experience

Analysis of flight data from previous missions indicates that anomalous Sun sensor readings could be caused by Earth albedo interference. A previous Sun sensor study presented a detailed mathematical model of this effect. The model can be used to study the effect of both diffusive and specular reflections and to improve Sun angle determination based on perturbed Sun sensor measurements, satellite position, and an approximate knowledge of attitude. The model predicts that diffuse reflected light can cause errors of up to 10 degrees in Coarse Sun Sensor (CSS) measurements and 5 to 10 arc sec in Fine Sun Sensor (FSS) measurements, depending on spacecraft orbit and attitude. The accuracy of these sensors is affected as long as part of the illuminated Earth surface is present in the sensor field of view. Digital Sun Sensors (DSS) respond in a different manner to the Earth albedo interference. Most of the time DSS measurements are not affected, but for brief periods of time the Earth albedo can cause errors which are a multiple of the sensor least significant bit and may exceed one degree. This paper compares model predictions with Tropical Rainfall Measuring Mission (TRMM) CSS measurements in order to validate and refine the model. Methods of reducing and mitigating the impact of Earth albedo are discussed. ne CSS sensor errors are roughly proportional to the Earth albedo coefficient. Photocells that are sensitive only to ultraviolet emissions would reduce the effective Earth albedo by up to a thousand times, virtually eliminating all errors caused by Earth albedo interference.

Brasoveanu, Dan↗

Reassessment of net radiation measurements in the atmosphere of Venus

An analysis was performed of the possible error mechanisms which degraded the IR net flux measurements made by the three small atmospheric probes dispatched from the Pioneer Venus spacecraft. The larger errors began below 30 km, and caused the data to be inconsistent with previous estimates of the atmospheric opacity. Evaluations were made of the possible radiation field perturbations behind each probe, cloud particle deposition on the sensor windows, and thermal disturbances within the sensors because of gas flow through the window retainers. The gas flow through the retainers was identified as the most likely error source, and was demonstrated in laboratory tests. A strong Reynolds number dependence was also found. Radiative transfer calculations were performed to account for the errors, using the constraints defined by the tests. Upper and lower bounds were calculated for the true net flux for both day and night conditions.

Revercomb, H. E.↗

Effect of anisotropic thermal conductivity on the morphological stability of a binary alloy

A linear morphological stability analysis is presented for a planar interface during unidirectional solidification of a binary alloy, in the case of a crystal having an anisotropic thermal conductivity. A dispersion relation shows that the onset of instability depends on the orientation of the growth direction with respect to principal crystallographic axes and on the orientation of the wavevector of the perturbation. The onset of instability can be either oscillatory (traveling waves) or nonoscillatory in time. For growth along a principal axis of the crystal there is an exchange of stabilities, and the onset of instability is nonoscillatory.

Coriell, S. R.↗