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 937 records · Page 52

Seasat orbit refinement for altimetry application

This paper describes the use of stochastic differential correction models in refining the Seasat orbit based on post-flight analysis of tracking data. The objective is to obtain orbital-height precision that is commensurate with the inherent Seasat altimetry data precision level of 10 cms. Local corrections to a mean ballistic arc, perturbed principally by atmospheric drag variations and local gravitational anomalies, are obtained by the introduction of stochastic dynamical models in conjunction with optimal estimation/smoothing techniques. Assessment of the resulting orbit with 'ground truth' provided by Seasat altimetry data shows that the orbital height precision is improved by 32% when compared to a conventional least-squares solution using the same data set. The orbital height precision realized by employing stochastic differential correction models is in the range of 73 cms to 208 cms rms.

Mohan, S. N.↗

Optimal heading change with minimum energy loss for a hypersonic gliding vehicle

A three state model is presented for analyzing the problem of optimal changes in heading with minimum energy loss for a hypersonic gliding vehicle. A further model order reduction to a single state model is examined using singular perturbation theory. The optimal solution for the reduced problem defines an optimal altitude profile dependent on the current energy of the vehicle, and the corresponding optimal lift and bank angle. A separate boundary layer analysis, based on an expansion of the necessary conditions about the reduced solution, is used to account for altitude and flight path angle dynamics and to derive a guidance law in feedback form. The guidance law is evaluated for a hypothetical vehicle.

Calise, Anthony J.↗

Buckling and Damage Resistance of Transversely-Loaded Composite Shells

Experimental and numerical work was conducted to better understand composite shell response to transverse loadings which simulate damage-causing impact events. The quasi-static, centered, transverse loading response of laminated graphite/epoxy shells in a [+/-45(sub n)/O(sub n)](sub s) layup having geometric characteristics of a commercial fuselage are studied. The singly-curved composite shell structures are hinged along the straight circumferential edges and are either free or simply supported along the curved axial edges. Key components of the shell response are response instabilities due to limit-point and/or bifurcation buckling. Experimentally, deflection-controlled shell response is characterized via load-deflection data, deformation-shape evolutions, and the resulting damage state. Finite element models are used to study the kinematically nonlinear shell response, including bifurcation, limit-points, and postbuckling. A novel technique is developed for evaluating bifurcation from nonlinear prebuckling states utilizing asymmetric spatial discretization to introduce numerical perturbations. Advantages of the asymmetric meshing technique (AMT) over traditional techniques include efficiency, robustness, ease of application, and solution of the actual (not modified) problems. The AMT is validated by comparison to traditional numerical analysis of a benchmark problem and verified by comparison to experimental data. Applying the technique, bifurcation in a benchmark shell-buckling problem is correctly identified. Excellent agreement between the numerical and experimental results are obtained for a number of composite shells although predictive capability decreases for stiffer (thicker) specimens which is attributed to compliance of the test fixture. Restraining the axial edge (simple support) has the effect of creating a more complex response which involves unstable bifurcation, limit-point buckling, and dynamic collapse. Such shells were noted to bifurcate into asymmetric deformation modes but were undamaged during testing. Shells in this study which were damaged were not observed to bifurcate. Thus, a direct link between bifurcation and atypical damage could not be established although the mechanism (bifurcation) was identified. Recommendations for further work in these related areas are provided and include extensions of the AMT to other shell geometries and structural problems.

Wardle, Brian L.↗

Detection and Characterization of Instrumental Transients in LISA Pathfinder and Their Projection to LISA

The LISA Pathfinder (LPF) mission succeeded outstandingly in demonstrating key technologicalaspects of future space-borne gravitational-wave detectors, such as the Laser Interferometer SpaceAntenna (LISA). Specifically, LPF demonstrated with unprecedented sensitivity the measurementof the relative acceleration of two free-falling cubic test masses. Although most disruptive non-gravitational forces have been identified and their effects mitigated through a series of calibrationprocesses, some faint transient signals of yet unexplained origin remain in the measurements. If theyappear in the LISA data, these perturbations (also called glitches) could skew the characterizationof gravitational-wave sources or even be confused with gravitational-wave bursts. For the firsttime, we provide a comprehensive census of LPF transient events. Our analysis is based on aphenomenological shapelet model allowing us to derive simple statistics about the physical featuresof the glitch population. We then implement a generator of synthetic glitches designed to be usedfor subsequent LISA studies, and perform a preliminary evaluation of the effect of the glitches onfuture LISA data analyses.

Quentin Baghi↗

Lapse-Rate Feeback: The Key to Reconciling TOA and Surface Attributions of Surface Warming

The impact of climate feedbacks on surface warming is conventionally evaluated using a decomposition of the top-of-atmosphere (TOA) energy budget. Alternatively, the climate feedback analysis can also be carried out using the surface energy budget. However, the two perspectives do not provide the same interpretation of process contributions to surface warming, particularly when executing a spatial analysis. The TOA energy budget is equal to the sum of the surface and atmospheric energy budgets. Using the CMIP5 RCP 8.5 model projections, we show that the major discrepancies between the surface and TOA climate feedback attributions of surface warming are due to non-negligible changes in the atmospheric energy budget that differ from their counterparts at the surface. Individual radiative and non-radiative processes cause vertically non-uniform energy flux perturbations, in response to an external forcing, that naturally lead to a vertically non-uniform temperature change. The TOA lapse-rate feedback is the manifestation of these multiple processes that produce a vertically non-uniform warming response such that it accounts for the asymmetry between the changes in the atmospheric and surface energy budgets. Using the climate feedback-response analysis method, we can decompose the lapse-rate feedback into contributions by individual processes. The negative lapse-rate feedback in the tropics is due to greater moist convection and condensational heating along with ocean heat storage. On the other hand, the positive lapse-rate feedback in polar regions is primarily due to the surface albedo and water vapor feedbacks. Combining the process contributions that are hidden within the lapse-rate feedback with their respective direct impacts on the TOA energy budget allows for a very consistent picture of process contributions to surface warming and its inter-model spread as that given by the surface energy budget approach. With the modified approach, both perspectives indicate water vapor feedback is the largest contributor in the tropics, while surface albedo feedback is the greatest contributor in polar regions. Dynamics and ocean heat storage are the main suppressors of surface warming, except over Antarctica. Both perspectives show there is large inter-model uncertainty in the contributions of water vapor, clouds, albedo, and dynamics plus ocean heat storage.

Sergio A Sejas↗

Analysis of Drift Errors in the JPL/UCLA Micromachined Gyroscope

The model of the JPL/UCLA micromachined vibratory gyroscope will be enhanced to include time varying effects. First, they will be shown to exist through trends in the experimental results. Causes of mechanical error will be further explained by analyzing possible perturbations to the physical model.

coriolois force gyroscope gyroscope design gyrosco↗

Development of Monte Carlo Capability for Orion Parachute Simulations

Parachute test programs employ Monte Carlo simulation techniques to plan testing and make critical decisions related to parachute loads, rate-of-descent, or other parameters. This paper describes the development and use of a MATLAB-based Monte Carlo tool for three parachute drop test simulations currently used by NASA. The Decelerator System Simulation (DSS) is a legacy 6 Degree-of-Freedom (DOF) simulation used to predict parachute loads and descent trajectories. The Decelerator System Simulation Application (DSSA) is a 6-DOF simulation that is well suited for modeling aircraft extraction and descent of pallet-like test vehicles. The Drop Test Vehicle Simulation (DTVSim) is a 2-DOF trajectory simulation that is convenient for quick turn-around analysis tasks. These three tools have significantly different software architectures and do not share common input files or output data structures. Separate Monte Carlo tools were initially developed for each simulation. A recently-developed simulation output structure enables the use of the more sophisticated DSSA Monte Carlo tool with any of the core-simulations. The task of configuring the inputs for the nominal simulation is left to the existing tools. Once the nominal simulation is configured, the Monte Carlo tool perturbs the input set according to dispersion rules created by the analyst. These rules define the statistical distribution and parameters to be applied to each simulation input. Individual dispersed parameters are combined to create a dispersed set of simulation inputs. The Monte Carlo tool repeatedly executes the core-simulation with the dispersed inputs and stores the results for analysis. The analyst may define conditions on one or more output parameters at which to collect data slices. The tool provides a versatile interface for reviewing output of large Monte Carlo data sets while preserving the capability for detailed examination of individual dispersed trajectories. The Monte Carlo tool described in this paper has proven useful in planning several Crew Exploration Vehicle parachute tests.

Moore, James W.↗

Galaxy formation by dust

It has been known since the early 1940's that radiation can cause an instability in the interstellar medium. Absorbing dust particles in an isotropic radiation field shadow each other by a solid angle which is inversely proportional to the square of the distance between the two particles, leading to an inverse-square attractive force - mock gravity. The effect is largest in an optically thin medium. Recently Hogan and White (HW, hereafter) proposed that if the pre-galactic universe contained suitable sources of radiation and dust, instability in the dust distribution caused by mock gravity may have led to the formation of galaxies and galaxy clusters. In their picture of a well-coupled dust-gas medium, HW show that mock gravity begins to dominate gravitational instability when the perturbation becomes optically thin, provided that the radiation field at the time is strong enough. The recent rocket observation of the microwave background at submillimeter wavelengths by Matsumoto et al. might be from pre-galactic stars, the consequence of the absorption of ultraviolet radiation by dust, and infrared reemission which is subsequently redshifted. HW's analysis omits radiative drag, incomplete collisional coupling of gas and dust, finite dust albedo, and finite matter pressure. These effects could be important. In a preliminary calculation including them, the authors have confirmed that mock gravitational instability is effective if there is a strong ultraviolet radiation at the time, but any galaxies that form would be substantially enriched in heavy elements because the contraction of the dust is more rapid than that of the gas. Moreover, since the dust moves with supersonic velocity through the gas soon after the perturbation becomes optically thin, the sputtering of dust particles by gas is significant, so the dust could disappear before the instability develops significantly. They conclude that the mock gravity by dust is not important in galaxy formations.

Wang, Boqi↗

Computational analysis of flow in 3D propulsive transition ducts

A numerical analysis of fully three dimensional, statistically steady flows in propulsive transition ducts being considered for use in future aircraft of higher maneuverability is investigated. The purpose of the transition duct is to convert axisymmetric flow from conventional propulsion systems to that of a rectangular geometry of high aspect ratio. In an optimal design, the transition duct would be of minimal length in order to reduce the weight penalty, while the geometrical change would be gradual enough to avoid detrimental flow perturbations. Recent experiments conducted at the Propulsion Aerodynamics Branch have indicated that thrust losses in ducts of superelliptic cross-section can be surprisingly low, even if flow separation occurs near the divergent walls. In order to address the objective of developing a rational design procedure for optimal transition ducts, it is necessary to have available a reliable computational tool for the analysis of flows achieved in a sequence of configurations. Current CFD efforts involving complicated geometries usually must contend with two separate but interactive aspects: namely, grid generation and flow solution. The first two avenues of the present investigation were comprised of suitable grid generation for a class of transition ducts of superelliptic cross-section, and the subsequent application of the flow solver PAB3D to this geometry. The code, PAB3D, was developed as a comprehensive tool for the solution of both internal and external high speed flows. The third avenue of investigation has involved analytical formulations to aid in the understanding of the nature of duct flows, and also to provide a basis of comparison for subsequent numerical solutions. Numerical results to date include the generation of two preliminary grid systems for duct flows, and the initial application of PAB3D to the corresponding geometries, which are of the class tested experimentally.

Sepri, Paavo↗

Numerical experiments in supersonic boundary-layer stability

The three-dimensional (3-D) time-dependent compressible Navier-Stokes equations are numerically solved by a Fourier-Chebyshev collocation method to study the stability of supersonic flows over a flat plate. Several direct simulations carried out in this study suggest the existence of a secondary instability that might provide a route to transition. The interaction of the modes involved in the secondary instability is possibly amenable to a Floquet-type analysis. Pertinent differences between this instability and the analogous incompressible K-type secondary instability are pointed out. Some preliminary results of a 2-D direct simulation of the nonlinear evolution of a second mode perturbation wave are also discussed.

Erlebacher, Gordon↗

Low order dynamical models for mixing layers

The goal of this project was to investigate the feasibility of developing useful low order dynamical systems to allow one to determine cheaply the effect of various initial conditions on the evolution of mixing layers. A mathematically rigorous approach for a time evolving base flow as found in a mixing layer seems to be the singular value decomposition of the linear operator transforming some initial perturbations into the solution at a given later time. Unfortunately, this is very costly, if at all feasible, and only deals with the linear evolution of a disturbance. A different approach based on weakly nonlinear analysis was attempted, but numerous difficulties were encountered. Several formulations were proposed which may be ill-posed. The reasons for this possible failure of the attempted approach is related to the fundamental limitations of weakly nonlinear analysis. The clarification of these issues is attempted in this report.

Riahi, D. N.↗

Pressure oscillation in the leakage annulus between a shrouded impeller and its housing due to impeller-discharge-pressure disturbances

The perturbed flow in the leakage path between a shrouded-pump impeller and its housing is analyzed using experiences with the Space Shuttle Main Engine (SSME), high pressure fuel turbopump (HPFTP) wearing-ring seals. Analysis is based on a bulk-flow model which consists of the path-momentum, circumferential momentum, and continuity equations. The pressure oscillations in the leakage annulus are driven by a circumferential variation of the impeller discharge pressure. It is shown that the occurrence and nature of the pressure oscillations depend on the tangential-velocity ratio of the fluid entering the seal, the order of the Fourier coefficient, the closeness of the precessional frequency of the rotating pressure field to the first natural frequency of the fluid annulus, and the clearance of the wearing-ring seal. The results obtained may explain the internal melting observed on SSME HPFTP seal parts.

Childs, D. W.↗

Modulational stability of periodic solutions of the Kuramoto-Sivaskinsky equation

We study the long-wave, modulational, stability of steady periodic solutions of the Kuramoto-Sivashinsky equation. The analysis is fully nonlinear at first, and can in principle be carried out to all orders in the small parameter, which is the ratio of the spatial period to a characteristic length of the envelope perturbations. In the linearized regime, we recover a high-order version of the results of Frisch, She, and Thual, which shows that the periodic waves are much more stable than previously expected.

Papageorgiou, Demetrios T.↗

The Effect of Three-Dimensional Freestream Disturbances on the Supersonic Flow Past a Wedge

The interaction between a shock wave (attached to a wedge) and small amplitude, three-dimensional disturbances of a uniform, supersonic, freestream flow are investigated. The paper extends the two-dimensional study of Duck et al, through the use of vector potentials, which render the problem tractable by the same techniques as in the two-dimensional case, in particular by expansion of the solution by means of a Fourier-Bessel series, in appropriately chosen coordinates. Results are presented for specific classes of freestream disturbances, and the study shows conclusively that the shock is stable to all classes of disturbances (i.e. time periodic perturbations to the shock do not grow downstream), provided the flow downstream of the shock is supersonic (loosely corresponding to the weak shock solution). This is shown from our numerical results and also by asymptotic analysis of the Fourier-Bessel series, valid far downstream of the shock.

Duck, Peter W.↗

Atmospheric Models for Aerocapture

There are eight destinations in the solar System with sufficient atmosphere for aerocapture to be a viable aeroassist option - Venus, Earth, Mars, Jupiter, Saturn and its moon Titan, Uranus, and Neptune. Engineering-level atmospheric models for four of these targets (Earth, Mars, Titan, and Neptune) have been developed for NASA to support systems analysis studies of potential future aerocapture missions. Development of a similar atmospheric model for Venus has recently commenced. An important capability of all of these models is their ability to simulate quasi-random density perturbations for Monte Carlo analyses in developing guidance, navigation and control algorithm, and for thermal systems design. Similarities and differences among these atmospheric models are presented, with emphasis on the recently developed Neptune model and on planned characteristics of the Venus model. Example applications for aerocapture are also presented and illustrated. Recent updates to the Titan atmospheric model are discussed, in anticipation of applications for trajectory and atmospheric reconstruct of Huygens Probe entry at Titan.

Justus, C. G.↗

Computational Analysis of the Flow and Acoustic Effects of Jet-Pylon Interaction

Computational simulation and prediction tools were used to understand the jet-pylon interaction effect in a set of bypass-ratio five core/fan nozzles. Results suggest that the pylon acts as a large scale mixing vane that perturbs the jet flow and jump starts the jet mixing process. The enhanced mixing and associated secondary flows from the pylon result in a net increase of noise in the first 10 diameters of the jet s development, but there is a sustained reduction in noise from that point downstream. This is likely the reason the pylon nozzle is quieter overall than the baseline round nozzle in this case. The present work suggests that focused pylon design could lead to advanced pylon shapes and nozzle configurations that take advantage of propulsion-airframe integration to provide additional noise reduction capabilities.

Hunter, Craig A.↗

Stability analysis of a liquid fuel annular combustion chamber

The problems of combustion instability in an annular combustion chamber are investigated. A modified Galerkin method was used to produce a set of modal amplitude equations from the general nonlinear partial differential acoustic wave equation. From these modal amplitude equations, the two variable perturbation method was used to develop a set of approximate equations of a given order of magnitude. These equations were modeled to show the effects of velocity sensitive combustion instabilities by evaluating the effects of certain parameters in the given set of equations. By evaluating these effects, parameters which cause instabilities to occur in the combustion chamber can be ascertained. It is assumed that in the annular combustion chamber, the liquid propellants are injected uniformly across the injector face, the combustion processes are distributed throughout the combustion chamber, and that no time delay occurs in the combustion processes.

Mcdonald, G. H.↗

Analysis of crustal structure of Venus utilizing residual Line-of-Sight (LOS) gravity acceleration and surface topography data. A trial of global modeling of Venus gravity field using harmonic spline method

To construct Venus' gravity disturbance field (or gravity anomaly) with the spacecraft-observer line of site (LOS) acceleration perturbation data, both a global and a local approach can be used. The global approach, e.g., spherical harmonic coefficients, and the local approach, e.g., the integral operator method, based on geodetic techniques are generally not the same, so that they must be used separately for mapping long wavelength features and short wavelength features. Harmonic spline, as an interpolation and extrapolation technique, is intrinsically flexible to both global and local mapping of a potential field. Theoretically, it preserves the information of the potential field up to the bound by sampling theorem regardless of whether it is global or local mapping, and is never bothered with truncation errors. The improvement of harmonic spline methodology for global mapping is reported. New basis functions, a singular value decomposition (SVD) based modification to Parker & Shure's numerical procedure, and preliminary results are presented.

Fang, Ming↗