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,117 records · Page 62

Ballistic Leverage for Conjunction Avoidance During the Lunar Transit Trajectory of NASA's Co-Manifested Vehicle

Analysis is completed to investigate the sensitivity of the Lunar Transit Trajectory of NASA’s Co-Manifested Vehicle (CMV) to the insertion of intentional coasting into an otherwise nearly continuous low thrust maneuver. A reference trajectory is developed which explicitly models low altitude thruster shutdowns in a level of detail not included in previous Design Reference Missions. This reference trajectory is used to investigate the sensitivity of the Lunar Transit to additional inserted coasting, presented in the form of separation of the perturbed trajectory from the reference as a function of mission elapsed time and duration and location of the coast segment. These results are shown to inform potential conjunction avoidance strategies during the CMV’s long low thrust spiral from a highly elliptic Earth parking orbit to the target Gateway Near Rectilinear HALO Orbit (NRHO).

low thrust↗

Effect of resonance-oblateness coupling on a satellite orbit

Second order effects of the coupling between geopotential resonance and oblateness on a satellite orbit are calculated. Results show that: (1) these effects arise from the interaction of resonance with the secular changes of the orbit's node, perigee, and mean anomaly; (2) they have the same period and phase as first order resonance perturbations; and (3) their amplitudes are proportional to the square of the period and dominate the first order effects as the orbit becomes commensurate. A striking example of this coupling is seen in the 18 day resonance variation of the node of the orbit of the first earth resources technology satellite. Analysis of this one arc second (31m) variation yielded a strong 14th order constraint to the geopotential for odd degree terms. This constraint is poorly predicted by current models.

Wagner, C. A.↗

Effect of resonance-oblateness coupling on a satellite orbit

Second-order effects of the coupling between geopotential resonance and oblateness on a satellite orbit are calculated. These effects arise from the interaction of resonance with the secular changes of the orbit's node, perigee, and mean anomaly. They have the same period and phase as first-order resonance perturbations. But their amplitudes are proportional to the square of the period and dominate the first-order effects as the orbit becomes commensurate. A striking example of this coupling is seen in the 18-day resonance variation of the node of the orbit of the first earth resources technology satellite. Analysis of this 1 arc sec (31 m) variation yields a strong fourteenth-order constraint to the geopotential field for odd degree terms. This constraint (lumped coefficient) is poorly predicted by most current models.

Wagner, C. A.↗

Validation of zero-order feedback strategies for medium range air-to-air interception in a horizontal plane

A zero order feedback solution of a variable speed interception game between two aircraft in the horizontal plane, obtained by using the method of forced singular perturbation (FSP), is compared with the exact open loop solution. The comparison indicates that for initial distances of separation larger than eight turning radii of the evader, the accuracy of the feedback approximation is better than one percent. The result validates the zero order FSP approximation for medium range air combat analysis.

Shinar, J.↗

Rotordynamic coefficients for compressible flow in tapered annular seals

Derivation of the governing equations for compressible flow in a tapered annular seal is based on Hirs' turbulent bulk-flow model. Zeroth and first-order perturbation equations are developed by an expansion in the eccentricity ratio. These equations are numerically integrated to obtain the leakage, and the direct and cross-coupled stiffness and damping coefficients. Seal parameters similar to the Space Shuttle Main Engine High Pressure Oxidizer Turbopump are used to demonstrate output from the analysis procedure. The effects of preswirl and seal taper are shown for three different length-to-diameter ratios. Generally the results indicate that prerotating the fluid significantly increases the cross-coupled stiffness but has little effect on the other coefficients, and increasing the convergent taper increases the direct stiffness while decreasing the direct damping and cross-coupled stiffness.

Nelson, C. C.↗

Earth's external magnetic fields at low orbital altitudes

Under our Jun. 1987 proposal, Magnetic Signatures of Near-Earth Distributed Currents, we proposed to render operational a modeling procedure that had been previously developed to compute the magnetic effects of distributed currents flowing in the magnetosphere-ionosphere system. After adaptation of the software to our computing environment we would apply the model to low altitude satellite orbits and would utilize the MAGSAT data suite to guide the analysis. During the first year, basic computer codes to run model systems of Birkeland and ionospheric currents and several graphical output routines were made operational on a VAX 780 in our research facility. Software performance was evaluated using an input matchstick ionospheric current array, field aligned currents were calculated and magnetic perturbations along hypothetical satellite orbits were calculated. The basic operation of the model was verified. Software routines to analyze and display MAGSAT satellite data in terms of deviations with respect to the earth's internal field were also made operational during the first year effort. The complete set of MAGSAT data to be used for evaluation of the models was received at the end of the first year. A detailed annual report in May 1989 described these first year activities completely. That first annual report is included by reference in this final report. This document summarizes our additional activities during the second year of effort and describes the modeling software, its operation, and includes as an attachment the deliverable computer software specified under the contract.

Klumpar, D. M.↗

Persistence and predictability in a perfect model

A realistic two-level GCM is used to examine the relationship between predictability and persistence. Predictability is measured by the average divergence of ensembles of solutions starting from perturbed initial conditions, and persistence is defined in terms of the autocorrelation function based on a single long-term model integration. The average skill of the dynamical forecasts is compared with the skill of simple persistence-based statistical forecasts. For initial errors comparable in magnitude to present-day analysis errors, the statistical forecast loses all skill after about one week, reflecting the lifetime of the lowest frequency fluctuations in the model. Large ensemble mean dynamical forecasts would be expected to remain skillful for about 3 wk. The disparity between the skill of the statistical and dynamical forecasts is greater for the higher frequency modes, which have little memory beyond 1 d, yet remain predictable for about 2 wk. The results are analyzed in terms of two characteristic time scales.

Schubert, Siegfried D.↗

Linearization of manipulator dynamics using spatial operators

Linearized dynamics models for manipulators are useful in robot analysis, motion planning, and control applications. Techniques from the spatial operator algebra are used to obtain closed form operator expressions for two types of linearized dynamics models, the linearized inverse and forward dynamics models. Spatially recursive algorithms of O(n) and O(n-squared) complexity for the computation of the perturbation vector and coefficient matrices for the linearized inverse dynamics model are developed first. Subsequently, operator factorization and inversion identities are used to develop corresponding closed-form expressions for the linearized forward dynamics model (LFDM). Once again, these are used to develop algorithms of O(n) and O(n-squared) complexity for the computation of the perturbation vector and the coefficient matrices. The algorithms for the LFDM do not require the explicit computation of the mass matrix nor its numerical inversion and are also of lower complexity than the conventional O(n-cubed) algorithms.

Jain, A.↗

Protein changes associated with reprotonation of the Schiff base in the photocycle of Asp96-->Asn bacteriorhodopsin. The MN intermediate with unprotonated Schiff base but N-like protein structure

The difference Fourier transform infrared spectrum for the N intermediate in the photoreaction of the light-adapted form of bacteriorhodopsin can be recorded at pH 10 at 274 K (Pfefferle, J.-M., Maeda, A., Sasaki, J., and Yoshizawa, T. (1991) Biochemistry 30, 6548-6556). Under these conditions, Asp96-->Asn bacteriorhodopsin gives a photoproduct which shows changes in protein structure similar to those observed in N of wild-type bacteriorhodopsin. However, decreased intensity of the chromophore bands and the single absorbance maximum at about 400 nm indicate that the Schiff base is unprotonated, as in the M intermediate. This photoproduct was named MN. At pH 7, where the supply of proton is not as restricted as at pH 10, Asp96-->Asn bacteriorhodopsin yields N with a protonated Schiff base. The Asn96 residue, which cannot deprotonate as Asp96 in wild-type bacteriorhodopsin, is perturbed upon formation of both MN at pH 10 and N at pH 7. We suggest that the reprotonation of the Schiff base is preceded by a large change in the protein structure including perturbation of the residue at position 96.

Non-NASA Center↗

Quantitative Assessment of Fatigue Damage Accumulation in Wavy Slip Metals from Acoustic Harmonic Generation

A comprehensive, analytical treatment is presented of the microelastic-plastic nonlinearities resulting from the interaction of a stress perturbation with dislocation substructures (veins and persistent slip bands) and cracks that evolve during high-cycle fatigue of wavy slip metals. The nonlinear interaction is quantified by a material (acoustic) nonlinearity parameter beta extracted from acoustic harmonic generation measurements. The contribution to beta from the substructures is obtained from the analysis of Cantrell [Cantrell, J. H., 2004, Proc. R. Soc. London A, 460, 757]. The contribution to beta from cracks is obtained by applying the Paris law for crack propagation to the Nazarov-Sutin crack nonlinearity equation [Nazarov, V. E., and Sutin, A. M., 1997, J. Acoust. Soc. Am. 102, 3349]. The nonlinearity parameter resulting from the two contributions is predicted to increase monotonically by hundreds of percent during fatigue from the virgin state to fracture. The increase in beta during the first 80-90 percent of fatigue life is dominated by the evolution of dislocation substructures, while the last 10-20 percent is dominated by crack growth. The model is applied to the fatigue of aluminium alloy 2024-T4 in stress-controlled loading at 276MPa for which experimental data are reported. The agreement between theory and experiment is excellent.

Cantrell, John H.↗

How Does Mediterranean Basin's Atmosphere Become Weak Moisture Source During Negative Phase of NAO: Use of AIRS, AMSR, TOVS, & TRMM Satellite Datasets Over Last Two NAO Cycles to Examine Governing Controls on E-P

The Mediterranean Sea is a noted 'concentration" basin in that it almost continuously exhibits positive evaporation minus precipitation (E - P ) properties -- throughout the four seasons and from one year to the next. Nonetheless, according to the ECMWF Era-40 48-year (1958-2005) climate reanalysis dataset, for various phases of the North Atlantic Oscillation (NAO) when the pressure gradient between Portugal and Iceland becomes either very relaxed (large negative NAO-Index) or in transition (small positive or negative NAO-Index), the atmospheric moisture source properties of the basin become weak, at times even reversed for several months (i.e., negative E - P). This behavior poses numerous questions concerning how and why these events occur. Moreover, it begs the question of what it would take for the basin to reach its tipping point in which P would exceed E throughout the rainy season (some six months) on an annually persistent basis -- and the sea would possibly transform to a recurring "dilution" basin. This talk investigates these questions by: (1) establishing over a period from 1979 to present, based on detailed analyses of satellite retrieval products from a combination of NASA-AQUA, NOAA-LEO, NASA/JAXA Scatterometer, and NASA-TRMM platforms, plus additional specialized satellite data products and ancillary meteorological datasets, the actual observation-based behavior of E - P, (2) diagnosing the salient physical and meteorological mechanisms that lead to the weaker E - P events during the analysis period, partly based on analyzing surface and upper air data at discrete stations in the western and eastern Mediterranean -- while at the same time evaluating the quality of the ERA-40 data over this same time period, (3) conducting GCM and high-resolution regional modeling experiments to determine if perturbed but realistic meteorological background conditions could maintain Mediterranean as a "dilution" basin through the October to March rainy season on annually recurring basis, and (4) investigating how such conditions might modify important internal and external climatic processes known to be closely related to the dynamical, thermal, and hydrological properties of the basin (e.g., drought frequency over Iberian peninsula, rainfall accumulation within Sahel, alteration of Levantine branch of east-west aligned open thermohaline cell, and modification of warm-salty intermediate flow through Gibraltar straight into North Atlantic).

Smith, Eric A.↗

Time-Distance Helioseismology Data-Analysis Pipeline for Helioseismic and Magnetic Imager Onboard Solar Dynamics Observatory (SDO-HMI) and Its Initial Results

The Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) provides continuous full-disk observations of solar oscillations. We develop a data-analysis pipeline based on the time-distance helioseismology method to measure acoustic travel times using HMI Doppler-shift observations, and infer solar interior properties by inverting these measurements. The pipeline is used for routine production of near-real-time full-disk maps of subsurface wave-speed perturbations and horizontal flow velocities for depths ranging from 0 to 20 Mm, every eight hours. In addition, Carrington synoptic maps for the subsurface properties are made from these full-disk maps. The pipeline can also be used for selected target areas and time periods. We explain details of the pipeline organization and procedures, including processing of the HMI Doppler observations, measurements of the travel times, inversions, and constructions of the full-disk and synoptic maps. Some initial results from the pipeline, including full-disk flow maps, sunspot subsurface flow fields, and the interior rotation and meridional flow speeds, are presented.

Sun: helioseismology↗

Adiabatic Gravitational Waveform Model for Compact Objects Undergoing Quasi-Circular Inspirals Into Rotating Massive Black Holes

We present bhpwave: a new python-based, open-source tool for generating the gravitational waveforms of stellar-mass compact objects undergoing quasicircular inspirals into rotating massive black holes. These binaries, known as extreme-mass-ratio inspirals (EMRIs), are exciting mHz gravitational wave sources for future space-based detectors such as the Laser Interferometer Space Antenna (LISA). Relativistic models of EMRI gravitational wave signals are necessary to unlock the full scientific potential of mHz detectors, yet few open-source EMRI waveform models exist. Thus we built bhpwave, which uses the adiabatic approximation from black hole perturbation theory to rapidly construct gravitational waveforms based on the leading-order inspiral dynamics of the binary. In this work, we present the theoretical and numerical foundations underpinning bhpwave. We also demonstrate how bhpwave can be used to assess the impact of EMRI modeling errors on LISA gravitational wave data analysis. In particular, we find that for retrograde orbits and slowly spinning black holes we can mismodel the gravitational wave phasing by as much as ∼10 radians without significantly biasing EMRI parameter estimation.

Zachary Nasipak↗

A linear thermal stability analysis for the vertical structure of alpha model accretion disks

The linear equations that describe the thermal stability of the convective and radiative vertical structure of an alpha model accretion disk are derived. For a specific model chosen to be representative of a dwarf nova disk, it is found that the eigenfunction for the temperature shows large excursions in the zone of partial ionization, but that the eigenfunctions for the unstable fundamental mode perturbations of the vertical distance and flux are nearly constant over the vertical extent of the disk. The growth of any instability is thus nearly homologous and the stability criteria derived for vertically averaged structures represent an adequate approximation.

Saio, Hideyuki↗

Development of separation due to interaction between a shock wave and a turbulent boundary layer perturbed by rarefaction waves

The interaction of a turbulent boundary layer with rarefaction and shock waves in flows past inclined steps was investigated experimentally and theoretically. The experiments were carried out in supersonic wind tunnels with test sections of 0.6 x 0.6 m and 0.2 x 0.2 m, with adiabatic conditions on the surface model. Based on detailed measurements of pressure fields, velocity, and surface friction, combined with results of optical visualization and analysis of limiting flow lines, three characteristic flow regimes are identified. These are (1) nonseparated flow at small step angles, (2) formation of a local separation zone with a free separation point at moderate step angles, and (3) formation of a large-scale separated flow with a fixed separation point at sufficiently large step angles. The quantitative gasdynamic schemes and surface pressure distributions for these flow regimes are shown for Mach 2.85.

Zheltovodov, A. A.↗

A method for the geometrically nonlinear analysis of compressively loaded prismatic composite structures

A method was developed for the geometrically nonlinear analysis of the static response of thin-walled stiffened composite structures loaded in uniaxial or biaxial compression. The method is applicable to arbitrary prismatic configurations composed of linked plate strips, such as stiffened panels and thin-walled columns. The longitudinal ends of the structure are assumed to be simply supported, and geometric shape imperfections can be modeled. The method can predict the nonlinear phenomena of postbuckling strength and imperfection sensitivity which are exhibited by some buckling-dominated structures. The method is computer-based and is semi-analytic in nature, making it computationally economical in comparison to finite element methods. The method uses a perturbation approach based on the use of a series of buckling mode shapes to represent displacement contributions associated with nonlinear response. Displacement contributions which are of second order in the model amplitudes are incorported in addition to the buckling mode shapes. The principle of virtual work is applied using a finite basis of buckling modes, and terms through the third order in the model amplitudes are retained. A set of cubic nonlinear algebraic equations are obtained, from which approximate equilibrium solutions are determined. Buckling mode shapes for the general class of structure are obtained using the VIPASA analysis code within the PASCO stiffened-panel design code. Thus, subject to some additional restrictions in loading and plate anisotropy, structures which can be modeled with respect to buckling behavior by VIPASA can be analyzed with respect to nonlinear response using the new method. Results obtained using the method are compared with both experimental and analytical results in the literature. The configurations investigated include several different unstiffened and blade-stiffening panel configurations, featuring both homogeneous, isotropic materials, and laminated composite material.

Stoll, Frederick↗

M2-SCREAM: A Stratospheric Composition Reanalysis of Aura MLS

Stratospheric composition and transport patterns are changing because of increasing greenhouse gas concentrations and declining levels of ozone depleting substances. Superimposed on forced trends are: dynamical variability at a range of time scales and unforeseen perturbations such as those resulting from injections of volcanic material and smoke from PyroCb events. Chemical reanalyses can be powerful tools for studying the composition of the stratosphere and disentangling forced responses from internal variability. This presentation introduces a new chemical reanalysis of stratospheric constituents developed and produced at NASA’s Global Modeling and Assimilation Office (GMAO) using the recently developed GMAO Constituent Data Assimilation System. Called the MERRA-2 Stratospheric Composition Reanalysis with Aura MLS (M2-SCREAM), this reanalysis consists of assimilated global three-dimensional fields of stratospheric ozone, water vapor, hydrogen chloride (HCl), nitric acid (HNO3) and nitrous oxide (N2O) mixing ratios and covers the period since the beginning of MLS observations in September 2004 to the present. M2-SCREAM assimilates version 4.2 MLS profiles of the five constituents alongside total ozone column from the Ozone Monitoring Instrument. The dynamics and tropospheric water vapor are constrained by assimilated meteorological fields from MERRA-2. The assimilated constituent fields and meteorology are provided to users at a 50-km horizontal resolution and a three-hourly frequency. Monthly analysis uncertainties and flags are also provided to the users. M2-SCREAM is an accurate and dynamically consistent high-resolution data record of the five constituents, all of which are of primary importance to stratospheric chemistry and transport studies. We will present a description of M2-SCREAM and selected results of a process-based evaluation of this product using independent data. As an illustration, we will discuss the perturbation to stratospheric composition from the eruption of Hunga Tonga–Hunga Ha'apai in January 2022 as seen in the reanalysis. We will also propose potential scientific applications of M2-SCREAM and outline plans for an upcoming comprehensive composition reanalysis that is being developed at NASA GMAO.

MERRA-2↗

LightForce Photon-Pressure Collision Avoidance: Efficiency Assessment on an Entire Catalogue of Space Debris

The potential to perturb debris orbits using photon pressure from ground-based lasers has been confirmed by independent research teams. Two useful applications of this scheme are protecting space assets from impacts with debris and stabilizing the orbital debris environment, both relying on collision avoidance rather than de-orbiting debris. This paper presents the results of a new assessment method to analyze the efficiency of the concept for collision avoidance. Earlier research concluded that one ground based system consisting of a 10 kW class laser, directed by a 1.5 m telescope with adaptive optics, can prevent a significant fraction of debris-debris collisions in low Earth orbit. That research used in-track displacement to measure efficiency and restricted itself to an analysis of a limited number of objects. As orbit prediction error is dependent on debris object properties, a static displacement threshold should be complemented with another measure to assess the efficiency of the scheme. In this paper we present the results of an approach using probability of collision. Using a least-squares fitting method, we improve the quality of the original TLE catalogue in terms of state and co-state accuracy. We then calculate collision probabilities for all the objects in the catalogue. The conjunctions with the highest risk of collision are then engaged by a simulated network of laser ground stations. After those engagements, the perturbed orbits are used to re-assess the collision probability in a 20 minute window around the original conjunction. We then use different criteria to evaluate the utility of the laser-based collision avoidance scheme and assess the number of base-line ground stations needed to mitigate a significant number of high probability conjunctions. Finally, we also give an account how a laser ground station can be used for both orbit deflection and debris tracking.

laser↗