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 793 records · Page 44

Robust control of systems with real parameter uncertainty and unmodelled dynamics

During this research period we have made significant progress in the four proposed areas: (1) design of robust controllers via H infinity optimization; (2) design of robust controllers via mixed H2/H infinity optimization; (3) M-delta structure and robust stability analysis for structured uncertainties; and (4) a study on controllability and observability of perturbed plant. It is well known now that the two-Riccati-equation solution to the H infinity control problem can be used to characterize all possible stabilizing optimal or suboptimal H infinity controllers if the optimal H infinity norm or gamma, an upper bound of a suboptimal H infinity norm, is given. In this research, we discovered some useful properties of these H infinity Riccati solutions. Among them, the most prominent one is that the spectral radius of the product of these two Riccati solutions is a continuous, nonincreasing, convex function of gamma in the domain of interest. Based on these properties, quadratically convergent algorithms are developed to compute the optimal H infinity norm. We also set up a detailed procedure for applying the H infinity theory to robust control systems design. The desire to design controllers with H infinity robustness but H(exp 2) performance has recently resulted in mixed H(exp 2) and H infinity control problem formulation. The mixed H(exp 2)/H infinity problem have drawn the attention of many investigators. However, solution is only available for special cases of this problem. We formulated a relatively realistic control problem with H(exp 2) performance index and H infinity robustness constraint into a more general mixed H(exp 2)/H infinity problem. No optimal solution yet is available for this more general mixed H(exp 2)/H infinity problem. Although the optimal solution for this mixed H(exp 2)/H infinity control has not yet been found, we proposed a design approach which can be used through proper choice of the available design parameters to influence both robustness and performance. For a large class of linear time-invariant systems with real parametric perturbations, the coefficient vector of the characteristic polynomial is a multilinear function of the real parameter vector. Based on this multilinear mapping relationship together with the recent developments for polytopic polynomials and parameter domain partition technique, we proposed an iterative algorithm for coupling the real structured singular value.

Chang, Bor-Chin↗

A microflare-related activation of a filament observed in H-alpha and C IV lines

A filament in active region AR 2717 was observed in two lines formed at different temperatures (H-alpha at 10 exp 4 K and C IV at 10 exp 5 K) with the multichannel-double-pass (MSDP) spectrograph and the ultraviolet spectrometer and polarimeter (UVSP). The partial disparition brusque (DB) of the filament as observed in H-alpha was due to the heating of a filament section in the vicinity of a bright point. A propagating disturbance followed this event. A detailed analysis of C IV rasters shows that this disturbance was not a passive perturbation, but was itself triggering active phenomena at various locations along its path, resulting in energy releases. It is suggested that this propagation of brightness was due to fast successive reconnections between fine looplike structures of the filament.

Schmieder, B.↗

A COMPACT model reduction methodology for articulated, multi-flexible bodies structures

To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one need appropriate reduced-order models for all the flexible components involved. The component modes projection and assembly model reduction methodology is one way to construct these reduced-order component models. However, the generated reduced-order system models typically contain high frequency extraneous modes. The presence of high frequency extraneous modes in a reduced-order system model forces us to use a small stepsize in the numerical integration of that model, resulting in an increase in the elapsed time of the integration process. To eliminate these high-frequency modes, we first identify and discard the component modes that 'contributed' to the generations of these modes. A sensitivity analysis is then used to determine how the remaining component modes should be perturbed in order to still closely capture the selected system modes. The effectiveness of the proposed methodology has been successfully verified using a finite-element model of the Galileo spacecraft.

Lee, Allan Y.↗

Benefit Assessment of the Precision Departure Release Capability Concept

A Precision Departure Release Capability concept is being evaluated by both the National Aeronautics and Space Administration and the Federal Aviation Administration as part of a larger goal of improving throughput, efficiency and capacity in integrated departure, arrival and surface operations. The concept is believed to have the potential of increasing flight efficiency and throughput by avoiding missing assigned slots and minimizing speed increase or path stretch to recover the slot. The main thrust of the paper is determining the impact of early and late departures from the departure runway when an aircraft has a slot assigned either at a meter fix or at the arrival airport. Results reported in the paper are for two scenarios. The first scenario considers flights out of Dallas/Fort Worth destined for Hartsfield-Jackson International Airport in Atlanta flying through the Meridian meter-fix in the Memphis Center with miles-in-trail constraints. The second scenario considers flights destined to George Bush Intercontinental/Houston Airport with specified airport arrival rate constraint. Results show that delay reduction can be achieved by allowing reasonable speed changes in scheduling. It was determined that the traffic volume between Dallas/Fort Worth and Atlanta via the Meridian fix is low and the departures times are spread enough that large departure schedule uncertainty can be tolerated. Flights can depart early or late within 90 minutes without accruing much more delay due to miles-in-trail constraint at the Meridian fix. In the Houston scenario, 808 arrivals from 174 airports were considered. Results show that delay experienced by the 16 Dallas/Fort Worth departures is higher if initial schedules of the remaining 792 flights are kept unaltered while they are rescheduled. Analysis shows that the probability of getting the initially assigned slot back after perturbation and rescheduling decreases with increasing standard deviation of the departure delay distributions. Results show that most Houston arrivals can be expected to be on time based on the assumed zero-mean Normal departure delay distributions achievable by Precision Departure Release Capability. In the current system, airport-departure delay, which is the sum of gate-departure delay and taxi-out delay, is observed at the airports. This delay acts as a bias, which can be reduced by Precision Departure Release Capability.

Palopo, Kee↗

Data Processing of Miniaturized Laser Heterodyne Radiometer (mini-LHR) Ground Instrument Retrievals

The Miniaturized Laser Heterodyne Radiometer (mini-LHR) is a passive ground instrument that observes the mole fraction of carbon dioxide (CO) and methane (CH)in the atmospheric column by measuring their absorption of sunlight at 1.6 microns. A laser heterodyne radiometer is similar in design to the super heterodyne radio receiver that is well known by ham radio enthusiasts. While not previously a commercial technique, laser heterodyne radiometers have a history of measuring atmospheric trace gases that started in the 1960s. With the commercial availability of inexpensive, low-power, thumbnail-sized lasers developed for the telecommunications industry, it was possible to miniaturize this technique and ultimately commercialize it. The mini-LHR has been under development at NASA GSFC since 2009. During that time, in addition to signicant technical improvements, processing has also evolved and been streamlined. Here we present details of the processing approach for raw mini-LHR data to produce 30- and 60-minute data products of CH and CO column mole fractions. Processing occurs in two general stages: a python-based pre-processing of raw data, followed by ingestion into a Planetary Spectrum Generator (PSG) retrieval algorithm. Raw data processing involves removal of outliers, correcting for changes in air mass throughout the day, averaging scans, and ultimately converting averaged scans into transmittance vs. wavelength. The PSG retrieval simulates a spectra for the time/day/location of the scan with meteorological inputs from Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) data set and then perturbs the concentrations of CO and CH to obtain a t based on an iterative least-squares curve fitting procedure.

Giancarlo Roberto Zambrano↗

Shock-ramp analysis test problem

Quasi-isentropic (ramp) compression is now a well-established experimental method and so are the analysis techniques to give Lagrangian sound speed, pressure, and density along the sample material's isentrope. A shock followed by ramp compression is a natural extension to investigate, for example, shock melt and refreeze on compression, or isentropes of states off the Hugoniot or principal isentrope. In practice, graded-density impactors produce initial shocks, compression by shaped laser pulses may be unable to produce a smooth pressure increase from zero, and incidental perturbations on the drive pulse may also give rise to shocks, so robust shock-ramp analysis methods will be needed. Appropriate analysis methods are needed for shock-ramp experiments, based on those for quasi-isentropic compression, and these require validation. This paper describes three different analyses of a shock-ramp test problem, including an assessment of their estimated errors. The methods tested were based on hydrodynamic characteristics or integration backward in space. All methods gave the known Lagrangian sound speed to within ~1%, and pressure and volume to within less than 2% and 1%, demonstrating that the analysis methods of isentropic compression experiments can be confidently extended to the analysis of shock and ramp compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

General framework for quantifying dissipation pathways in open quantum systems. III. Off-diagonal subsystem–bath couplings

This paper extends the previously reported theory of dissipation pathways [C. W. Kim and I. Franco, J. Chem. Phys. 160, 214111 (2024)] to incorporate off-diagonal subsystem–bath coupling, which is often required to model molecular systems where the environment directly influences transitions and couplings between subsystem states. We systematically derive master equations for both population transfer and dissipation into individual bath components, for which we also rigorously prove energy conservation and detailed balance. The approach is based on second-order perturbation theory with respect to the subsystem–bath couplings, whose form is not limited to any specific model. The accuracy of the developed method is tested by applying it to diverse model Hamiltonians involving linearly coupled harmonic oscillator baths and comparing the outcomes against the hierarchical equations of motion (HEOM) method. Overall, our method accurately quantifies the contributions of specific bath components to the overall dissipation while significantly reducing the computational cost compared to numerically exact methods such as HEOM, thus offering a path to examine how vibronic interactions steer non-adiabatic processes in realistic chemical systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical predictions of latitude dependencies in the solar wind

Results are presented which were obtained with the Winge-Coleman model for theoretical predictions of latitudinal dependencies in the solar wind. A first-order expansion is described which allows analysis of first-order latitudinal variations in the coronal boundary conditions and results in a second-order partial differential equation for the perturbation stream function. Latitudinal dependencies are analytically separated out in the form of Legendre polynomials and their derivative, and are reduced to the solution of radial differential equations. This analysis is shown to supply an estimate of how large the coronal variation in latitude must be to produce an 11 km/sec/deg gradient in the radial velocity of the solar wind, assuming steady-state processes.

Winge, C. R., Jr.↗

The use of atmospheric measurements to constrain model predictions of ozone change from chlorine perturbations

Atmospheric photochemistry models have been used to predict the sensitivity of the ozone layer to various perturbations. These same models also predict concentrations of chemical species in the present day atmosphere which can be compared to observations. Model results for both present day values and sensitivity to perturbation depend upon input data for reaction rates, photodissociation rates, and boundary conditions. A method of combining the results of a Monte Carlo uncertainty analysis with the existing set of present atmospheric species measurements is developed. The method is used to examine the range of values for the sensitivity of ozone to chlorine perturbations that is possible within the currently accepted ranges for input data. It is found that model runs which predict ozone column losses much greater than 10 percent as a result of present fluorocarbon fluxes produce concentrations and column amounts in the present atmosphere which are inconsistent with the measurements for ClO, HCl, NO, NO2, and HNO3.

Douglass, Anne R.↗

Score-based denoising for atomic structure identification

We propose an effective method for removing thermal vibrations that complicate the task of analyzing complex dynamics in atomistic simulation of condensed matter. Our method iteratively subtracts thermal noises or perturbations in atomic positions using a denoising score function trained on synthetically noised but otherwise perfect crystal lattices. The resulting denoised structures clearly reveal underlying crystal order while retaining disorder associated with crystal defects. Purely geometric, agnostic to interatomic potentials, and trained without inputs from explicit simulations, our denoiser can be applied to simulation data generated from vastly different interatomic interactions. The denoiser is shown to improve existing classification methods, such as common neighbor analysis and polyhedral template matching, reaching perfect classification accuracy on a recent benchmark dataset of thermally perturbed structures up to the melting point. Demonstrated here in a wide variety of atomistic simulation contexts, the denoiser is general, robust, and readily extendable to delineate order from disorder in structurally and chemically complex materials.

36 MATERIALS SCIENCE↗

An adjoint sensitivity study of blocking in a two-layer isentropic model

This paper presents a new methodology for adjoint sensitivity analysis, previously developed in general terms by Cacuci, into a form directly applicable 10 meteorological problems. This technique is illustrated by examining the sensitivity of a blocking index in a two-layer isentropic model. The index represents a response function for the sensitivity analysis that unlike previous meteorological applications, is an operator and not a functional, and thus extends the scope of adjoint sensitivity to general operator-type responses depending on time and/or space. The sensitivity of the blocking index to perturbations introduced into the model atmosphere, as well as to model parameters, is discussed. The methodology of generalized adjoint sensitivity analysis described in this paper constitutes a prototype for further applications in the atmospheric and/or oceanic sciences.

Zou, X.↗

A Design Study of Onboard Navigation and Guidance During Aerocapture at Mars

The navigation and guidance of a high lift-to-drag ratio sample return vehicle during aerocapture at Mars are investigated. Emphasis is placed on integrated systems design, with guidance algorithm synthesis and analysis based on vehicle state and atmospheric density uncertainty estimates provided by the navigation system. The latter utilizes a Kalman filter for state vector estimation, with useful update information obtained through radar altimeter measurements and density altitude measurements based on IMU-measured drag acceleration. A three-phase guidance algorithm, featuring constant bank numeric predictor/corrector atmospheric capture and exit phases and an extended constant altitude cruise phase, is developed to provide controlled capture and depletion of orbital energy, orbital plane control, and exit apoapsis control. Integrated navigation and guidance systems performance are analyzed using a four degree-of-freedom computer simulation. The simulation environment includes an atmospheric density model with spatially correlated perturbations to provide realistic variations over the vehicle trajectory. Navigation filter initial conditions for the analysis are based on planetary approach optical navigation results. Results from a selection of test cases are presented to give insight into systems performance.

Fuhry, Douglas Paul↗

Eigenmode analysis of unsteady one-dimensional Euler equations

The initial boundary value problem describing the evolution of unsteady linearized perturbations of a steady, uniform subsonic flow is analyzed. The eigenmodes and eigenfrequencies of the system are derived and several examples are presented to illustrate the effect of different boundary conditions on the exponential decay rate of the eigenmodes. The resultant implications for the stability and convergence rates of finite difference computations are discussed.

Giles, M.↗

Adaptation to transient postural perturbations

This research was first proposed in May, 1986, to focus on some of the problems encountered in the analysis of postural responses gathered from crewmembers. The ultimate driving force behind this line of research was the desire to treat, predict, or explain 'Space Adaptation Syndrome' (SAS) and hence circumvent any adverse effects of space motion sickness on crewmember performance. The aim of this project was to develop an easily implemented analysis of the transient responses to platform translation that can be elicited with a protocol designed to force sensorimotor reorganization, utilizing statistically reliable criterion measures. This report will present: (1) a summary of the activity that took place in each of the three funded years of the project; (2) discussion of experimental results and their implications for future research; and (3) a list of presentations and publications resulting from this project.

Andres, Robert O.↗

Analysis of the effects of mean local node-crossing time on the evolution of Sun-synchronous orbits

An investigation of the effect of mean local node-crossing time on the evolution of Sun-synchronous orbits was undertaken during Phase-A orbit analysis for the National Oceanic and Atmospheric Administration (NOAA) O,P,Q environmental spacecraft. That analysis added to the growing body of evidence that individual Sun-synchronous missions, at differing node-crossing times, experience nodal drift rates that can differ in both magnitude and direction. A Sun-synchronous orbit is obtained by means of a nodal drift rate approximating the 0.9856-degree-per-day apparent precession of the position of the mean Sun. This drift rate is achieved through the interaction of the orbital semimajor axis and inclination in Earth's geopotential field. Influencing perturbations include atmospheric drag and, most important, the effects of solar gravitation on inclination. The present analysis examines a series of Sun-synchronous orbits with mean local node-crossing times at 1-hour intervals from 6 a.m. to 6 p.m. It considers the fixed geometry of each orbital plane with respect to both the Sun and the diurnal atmospheric bulge, then analyzes the influence of these features upon the evolution of the semimajor axix and inclination and thus upon the rate of the nodal drift in the course of 1 year.

Jordan, P.↗

Physics-Based Methods of Failure Analysis and Diagnostics in Human Space Flight

The Integrated Health Management (IHM) for the future aerospace systems requires to interface models of multiple subsystems in an efficient and accurate information environment at the earlier stages of system design. The complexity of modern aeronautic and aircraft systems (including e.g. the power distribution, flight control, solid and liquid motors) dictates employment of hybrid models and high-level reasoners for analysing mixed continuous and discrete information flow involving multiple modes of operation in uncertain environments, unknown state variables, heterogeneous software and hardware components. To provide the information link between key design/performance parameters and high-level reasoners we rely on development of multi-physics performance models, distributed sensors networks, and fault diagnostic and prognostic (FD&P) technologies in close collaboration with system designers. The main challenges of our research are related to the in-flight assessment of the structural stability, engine performance, and trajectory control. The main goal is to develop an intelligent IHM that not only enhances components and system reliability, but also provides a post-flight feedback helping to optimize design of the next generation of aerospace systems. Our efforts are concentrated on several directions of the research. One of the key components of our strategy is an innovative approach to the diagnostics/prognostics based on the real time dynamical inference (DI) technologies extended to encompass hybrid systems with hidden state trajectories. The major investments are into the multiphysics performance modelling that provides an access of the FD&P technologies to the main performance parameters of e.g. solid and liquid rocket motors and composite materials of the nozzle and case. Some of the recent results of our research are discussed in this chapter. We begin by introducing the problem of dynamical inference of stochastic nonlinear models and reviewing earlier results. Next, we present our analytical approach to the solution of this problem based on the path integral formulation. The resulting algorithm does not require an extensive global search for the model parameters, provides optimal compensation for the effects of dynamical noise, and is robust for a broad range of dynamical models. In the following Section the strengths of the algorithm are illustrated illustrated by inferring the parameters of the stochastic Lorenz system and comparing the results with those of earlier research. Next, we discuss a number of recent results in application to the development of the IHM for aerospace system. Firstly, we apply dynamical inference approach to a solution of classical three tank problems with mixed unknown continuous and binary parameters. The problem is considered in the context of ground support system for filling fuel tanks of liquid rocket motors. It is shown that the DI algorithm is well suited for successful solution of a hybrid version of this benchmark problem even in the presence of additional periodic and stochastic perturbation of unknown strength. Secondly, we illustrate our approach by its application to an analysis of the nozzle fault in a solid rocket motor (SRM). The internal ballistics of the SRM is modelled as a set of one-dimensional partial differential equations coupled to the dynamics of the propellant regression. In this example we are specifically focussed on the inference of discrete and continuous parameters of the nozzle blocking fault and on the possibility of an application of the DI algorithm to reducing the probability of "misses" of an on-board FD&P for SRM. In the next section re-contact problem caused by first stage/upper stage separation failure is discussed. The reaction forces imposed on the nozzle of the upper stage during the re-contact and their connection to the nozzle damage and to the thrust vector control (TVC) signal are obtained. It is shown that transient impact induced torquean be modelled as a response of an effective damped oscillator. A possible application of the DI algorithm to the inference of damage parameters and predicting fault dynamics ahead of time using the actuator signal is discussed. Finally, we formulate Bayesian inferential framework for development of the IHM system for in-flight structural health monitoring (SHM) of composite materials. We consider the signal generated by piezoelectric actuator mounted on composite structure generating elastic waves in it. The signal received by the sensor is than compared with the baseline signal. The possibility of damage inference is discussed in the context of development of the SHM.

Smelyanskiy, Vadim N.↗

Investigating the impact of intermediate-mode perturbations on diagnosing plasma conditions in DT cryogenic implosions via synthetic x-ray Thomson scattering

The pursuit of inertial confinement fusion ignition target designs requires precise experimental validation of the conditions within imploding capsules, in particular the density and temperature of the compressed shell. Previous work has identified x-ray Thomson scattering (XRTS) as a viable diagnostic tool for inferring the in-flight compressed deuterium-tritium shell conditions during capsule implosions (Poole et al 2022 Phys. Plasmas 29 072703). However, this study focused on one-dimensional simulations, which do not account for the growth of hydrodynamic instabilities. In this work, two-dimensional DRACO simulations incorporating intermediate-mode perturbations up to Legendre mode l = 50 were used to generate synthetic XRTS spectra with the SPECT3D code. The analysis employed Markov-Chain Monte Carlo techniques to infer plasma conditions from these spectra. The results demonstrate that the XRTS diagnostic platform can effectively discern the in-flight compressed shell conditions for targets with varying adiabats, even in the presence of intermediate-mode perturbations. This work underscores the potential of XRTS for realistic inertial confinement fusion experiments, providing a robust method for probing the complex dynamics of fusion implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗