Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear inverse model”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Grail Refinements to Lunar Seismic Structure

To probe a planet's interior, seismology provides the most direct constraints on the variables that govern the dynamic properties of the body. However, the GRAIL (Gravity Recovery and Interior Laboratory) mission's high-resolution measurements of the lunar gravity field provide constraints on crustal thickness, mantle structure, core radius and stratification, and core state (solid vs. molten). These data complement seismic investigations, and joint interpretation permits improved constraints on the Moon's internal structure. Joint interpretation of disparate geophysical datasets helps reduce drawbacks that can result from analyzing them individually. The Apollo seismic network was situated on the lunar nearside surface in a roughly equilateral triangle having sides approximately 1000 km long, with stations 12/14 nearly co-located at one corner. Due to this limited geographical extent, near-surface ray coverage from moonquakes is low, but increase with depth. In comparison, gravity surveys and their resulting gravity anomaly maps have traditionally offered optimal resolution at crustal depths. Gravimetric maps and seismic data sets are therefor well suited to joint inversion, since the complementary information reduces inherent model ambiguity. We will perform a joint inversion of Apollo seismic delay times and gravity data collected by GRAIL lunar gravity mission, in order to recover seismic velocity and density as a function of latitude, longitude and depth within the Moon. We will relate density (rho) to seismic velocity (v) using a linear relationship that is allowed to be depth-dependent. The corresponding coefficient (B) can reflect a variety of material properties that vary with depth, including temperature and composition. The inversion seeks to recover the set of rho, v, and B perturbations that minimize (in a least-squares sense) the difference between the observed and calculated data.

Weber, Renee C.↗
NASA NTRS1994

TOPEX/POSEIDON tides estimated using a global inverse model

Altimetric data from the TOPEX/POSEIDON mission will be used for studies of global ocean circulation and marine geophysics. However, it is first necessary to remove the ocean tides, which are aliased in the raw data. The tides are constrained by the two distinct types of information: the hydrodynamic equations which the tidal fields of elevations and velocities must satisfy, and direct observational data from tide gauges and satellite altimetry. Here we develop and apply a generalized inverse method, which allows us to combine rationally all of this information into global tidal fields best fitting both the data and the dynamics, in a least squares sense. The resulting inverse solution is a sum of the direct solution to the astronomically forced Laplace tidal equations and a linear combination of the representers for the data functionals. The representer functions (one for each datum) are determined by the dynamical equations, and by our prior estimates of the statistics or errors in these equations. Our major task is a direct numerical calculation of these representers. This task is computationally intensive, but well suited to massively parallel processing. By calculating the representers we reduce the full (infinite dimensional) problem to a relatively low-dimensional problem at the outset, allowing full control over the conditioning and hence the stability of the inverse solution. With the representers calculated we can easily update our model as additional TOPEX/POSEIDON data become available. As an initial illustration we invert harmonic constants from a set of 80 open-ocean tide gauges. We then present a practical scheme for direct inversion of TOPEX/POSEIDON crossover data. We apply this method to 38 cycles of geophysical data records (GDR) data, computing preliminary global estimates of the four principal tidal constituents, M(sub 2), S(sub 2), K(sub 1) and O(sub 1). The inverse solution yields tidal fields which are simultaneously smoother, and in better agreement with altimetric and ground truth data, than previously proposed tidal models. Relative to the 'default' tidal corrections provided with the TOPEX/POSEIDON GDR, the inverse solution reduces crossover difference variances significantly (approximately 20-30%), even though only a small number of free parameters (approximately equal to 1000) are actually fit to the crossover data.

Egbert, Gary D.↗
DOE OSTI2022

Translation and rotation equivariant normalizing flow (TRENF) for optimal cosmological analysis

ABSTRACT Our Universe is homogeneous and isotropic, and its perturbations obey translation and rotation symmetry. In this work, we develop translation and rotation equivariant normalizing flow (TRENF), a generative normalizing flow (NF) model which explicitly incorporates these symmetries, defining the data likelihood via a sequence of Fourier space-based convolutions and pixel-wise non-linear transforms. TRENF gives direct access to the high dimensional data likelihood p(x|y) as a function of the labels y, such as cosmological parameters. In contrast to traditional analyses based on summary statistics, the NF approach has no loss of information since it preserves the full dimensionality of the data. On Gaussian random fields, the TRENF likelihood agrees well with the analytical expression and saturates the Fisher information content in the labels y. On non-linear cosmological overdensity fields from N-body simulations, TRENF leads to significant improvements in constraining power over the standard power spectrum summary statistic. TRENF is also a generative model of the data, and we show that TRENF samples agree well with the N-body simulations it trained on, and that the inverse mapping of the data agrees well with a Gaussian white noise both visually and on various summary statistics: when this is perfectly achieved the resulting p(x|y) likelihood analysis becomes optimal. Finally, we develop a generalization of this model that can handle effects that break the symmetry of the data, such as the survey mask, which enables likelihood analysis on data without periodic boundaries.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Elastic robot control - Nonlinear inversion and linear stabilization

An approach to the control of elastic robot systems for space applications using inversion, servocompensation, and feedback stabilization is presented. For simplicity, a robot arm (PUMA type) with three rotational joints is considered. The third link is assumed to be elastic. Using an inversion algorithm, a nonlinear decoupling control law u(d) is derived such that in the closed-loop system independent control of joint angles by the three joint torquers is accomplished. For the stabilization of elastic oscillations, a linear feedback torquer control law u(s) is obtained applying linear quadratic optimization to the linearized arm model augmented with a servocompensator about the terminal state. Simulation results show that in spite of uncertainties in the payload and vehicle angular velocity, good joint angle control and damping of elastic oscillations are obtained with the torquer control law u = u(d) + u(s).

Singh, S. N.↗
DOE OSTI2024

Revealing Decision Conservativeness Through Inverse Distributionally Robust Optimization

This paper introduces Inverse Distributionally Robust Optimization (I-DRO) as a method to infer the conservativeness level of a decision-maker, represented by the size of a Wasserstein metric-based ambiguity set, from the optimal decisions made using Forward Distributionally Robust Optimization (F-DRO). By leveraging the Karush-Kuhn-Tucker (KKT) conditions of the convex F-DRO model, we formulate I-DRO as a bi-linear program, which can be solved using off-the-shelf optimization solvers. Additionally, this formulation exhibits several advantageous properties. We demonstrate that I-DRO not only guarantees the existence and uniqueness of an optimal solution but also establishes the necessary and sufficient conditions for this optimal solution to accurately match the actual conservativeness level in F-DRO. Furthermore, we identify three extreme scenarios that may impact I-DRO effectiveness. Our case study applies F-DRO for power system scheduling under uncertainty and employs I-DRO to recover the conservativeness level of system operators. Numerical experiments based on an IEEE 5-bus system and a realistic NYISO 11-zone system demonstrate I-DRO performance in both normal and extreme scenarios. An extended version of this paper with additional analyses is available at li2024revealing.

distributionally robust optimization↗
NASA NTRS2017

Variability and Quasi-Decadal Changes in the Methane Budget over the Period 2000-2012

Following the recent Global Carbon Project (GCP) synthesis of the decadal methane (CH4) budget over 2000- 2012, we analyse here the same dataset with a focus on quasi-decadal and inter-annual variability in CH4 emissions. The GCP dataset integrates results from top-down studies (exploiting atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up models (including process-based models for estimating land surface emissions and atmospheric chemistry), inventories of anthropogenic emissions, and data-driven approaches. The annual global methane emissions from top-down studies, which by construction match the observed methane growth rate within their uncertainties, all show an increase in total methane emissions over the period 2000-2012, but this increase is not linear over the 13 years. Despite differences between individual studies, the mean emission anomaly of the top-down ensemble shows no significant trend in total methane emissions over the period 2000-2006, during the plateau of atmospheric methane mole fractions, and also over the period 2008-2012, during the renewed atmospheric methane increase. However, the top-down ensemble mean produces an emission shift between 2006 and 2008, leading to 22 [16-32] Tg CH4 yr(exp -1) higher methane emissions over the period 2008-2012 compared to 2002-2006. This emission increase mostly originated from the tropics, with a smaller contribution from mid-latitudes and no significant change from boreal regions. The regional contributions remain uncertain in top-down studies. Tropical South America and South and East Asia seem to contribute the most to the emission increase in the tropics. However, these two regions have only limited atmospheric measurements and remain therefore poorly constrained. The sectorial partitioning of this emission increase between the periods 2002-2006 and 2008-2012 differs from one atmospheric inversion study to another. However, all top-down studies suggest smaller changes in fossil fuel emissions (from oil, gas, and coal industries) compared to the mean of the bottom-up inventories included in this study. This difference is partly driven by a smaller emission change in China from the top-down studies compared to the estimate in the Emission Database for Global Atmospheric Research (EDGARv4.2) inventory, which should be revised to smaller values in a near future. We apply isotopic signatures to the emission changes estimated for individual studies based on five emission sectors and find that for six individual top-down studies (out of eight) the average isotopic signature of the emission changes is not consistent with the observed change in atmospheric 13CH4. However, the partitioning in emission change derived from the ensemble mean is consistent with this isotopic constraint. At the global scale, the top-down ensemble mean suggests that the dominant contribution to the resumed atmospheric CH4 growth after 2006 comes from microbial sources (more from agriculture and waste sectors than from natural wetlands), with an uncertain but smaller contribution from fossil CH4 emissions. In addition, a decrease in biomass burning emissions (in agreement with the biomass burning emission databases) makes the balance of sources consistent with atmospheric 13CH4 observations. In most of the top-down studies included here, OH concentrations are considered constant over the years (seasonal variations but without any inter-annual variability). As a result, the methane loss (in particular through OH oxidation) varies mainly through the change in methane concentrations and not its oxidants. For these reasons, changes in the methane loss could not be properly investigated in this study, although it may play a significant role in the recent atmospheric methane changes as briefly discussed at the end of the paper.

Saunois, Marielle↗
DOE OSTI2024

Development of the tangent linear and adjoint models of the global online chemical transport model MPAS-CO 2 v7.3

We describe the development of the tangent linear (TL) and adjoint models of the Model for Prediction Across Scales (MPAS)-CO 2 transport model, which is a global online chemical transport model developed upon the non-hydrostatic Model for Prediction Across Scales – Atmosphere (MPAS-A). The primary goal is to make the model system a valuable research tool for investigating atmospheric carbon transport and inverse modeling. First, we develop the TL code, encompassing all CO 2 transport processes within the MPAS-CO 2 forward model. Then, we construct the adjoint model using a combined strategy involving re-calculation and storage of the essential meteorological variables needed for CO 2 transport. This strategy allows the adjoint model to undertake a long-period integration with moderate memory demands. To ensure accuracy, the TL and adjoint models undergo vigorous verifications through a series of standard tests. The adjoint model, through backward-in-time integration, calculates the sensitivity of atmospheric CO 2 observations to surface CO 2 fluxes and the initial atmospheric CO 2 mixing ratio. To demonstrate the utility of the newly developed adjoint model, we conduct simulations for two types of atmospheric CO 2 observations, namely the tower-based in situ CO 2 mixing ratio and satellite-derived column-averaged CO 2 mixing ratio (X CO 2 ). A comparison between the sensitivity to surface flux calculated by the MPAS-CO 2 adjoint model with its counterpart from CarbonTracker–Lagrange (CT-L) reveals a spatial agreement but notable magnitude differences. These differences, particularly evident for X CO 2 , might be attributed to the two model systems' differences in the simulation configuration, spatial resolution, and treatment of vertical mixing processes. Moreover, this comparison highlights the substantial loss of information in the atmospheric CO 2 observations due to CT-L's spatial domain limitation. Furthermore, the adjoint sensitivity analysis demonstrates that the sensitivities to both surface flux and initial CO 2 conditions spread out throughout the entire Northern Hemisphere within a month. MPAS-CO 2 forward, TL, and adjoint models stand out for their calculation efficiency and variable-resolution capability, making them competitive in computational cost. In conclusion, the successful development of the MPAS-CO 2 TL and adjoint models, and their integration into the MPAS-CO 2 system, establish the possibility of using MPAS's unique features in atmospheric CO 2 transport sensitivity studies and in inverse modeling with advanced methods such as variational data assimilation.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1981

A model for the acoustic impedance of linear suppressor materials bonded on perforated plate

An analytical model is developed for the acoustic impedance of a composite made of a fine screen bonded to a perforated plate, the resistance having been modeled using steady flow theory and data, and including both linear and nonlinear components. The model of the screen resistance alone is derived first, and it is then extended to include a perforated plate bonded to the screen. Quantities such as effective open ratio and pressure recovery are evaluated. Conclusions are presented, including: (1) when the screen is bonded to a perforated plate, the linear resistance increases as the inverse of the perforate open area, (2) the nonlinearity in the screen resistance depends on the square of the screen wire diameter, and (3) when the resistance is considered as a function of approach velocity, the screen nonlinear term increases as the inverse of the square of the perforated plate open area ratio. The general theory can be applied to any screen structure, but final results are limited to a twilled square weave screen. Results can be applied to the study of sound absorbing materials used for lining the walls of aircraft engine duct suppressors.

Rice, E. J.↗
NASA NTRS1982

Inversion of data from diffraction-limited multiwavelength remote sensors. II - Nonlinear dependence of observables on the geophysical parameters

Linear shift-invariant spatial filtering is applied to inversion of radiometric measurements of the earth in both polarizations at the following frequencies: 6.6, 10.7, 18, 21, and 37 GHz. For purposes of radiative transfer calculations at these frequencies, the state of the ocean-atmosphere system is described by a seven-parameter model. The parameters are near-surface wind speed, sea surface temperature, integrated water vapor mass, scale height of water vapor in an exponential distribution, integrated liquid water mass, height of the liquid water, and characteristic drop radius in a Best (1950) drop size distribution. Also described is a spatial filtering algorithm by means of which geophysical parameters are retrieved from measurements made over the ocean by the scanning multichannel microwave radiometer on the Nimbus 7 satellite. The parameters here are near-surface wind speed, sea surface temperature, rain rate, integrated atmospheric water vapor content, and integrated liquid water content (including rain).

Rosenkranz, P. W.↗
NASA NTRS1985

Robust torque control of an elastic robotic arm based on invertibility and feedback stabilization

An approach to the control of elastic robotic systems for space applications using inversion, servocompensation, and feedback stabilization is presented. For simplicity, a robot arm (PUMA-type) with three rotational joints is considered. The third link is assumed to be elastic. Using an inversion algorithm, a nonlinear decoupling control law, u sub d, is derived such that in the closed loop system, independent control of joint angles by the three joint torquers is accomplished. For the stabilization of elastic oscillations, a linear feedback torquer control law, u sub s, is obtained applying linear quadratic optimization to the linearized arm model augmented with a servocompensator about the terminal state. Simulation results show that, in spite of uncertainties in the payload and vehicle angular velocity, good joint angle control and damping of elastic oscillations are obtained with the torquer control law u = u sub d + u sub s.

Singh, S. N.↗
NASA NTRS1988

Bayesian inference in geomagnetism

The inverse problem in empirical geomagnetic modeling is investigated, with critical examination of recently published studies. Particular attention is given to the use of Bayesian inference (BI) to select the damping parameter lambda in the uniqueness portion of the inverse problem. The mathematical bases of BI and stochastic inversion are explored, with consideration of bound-softening problems and resolution in linear Gaussian BI. The problem of estimating the radial magnetic field B(r) at the earth core-mantle boundary from surface and satellite measurements is then analyzed in detail, with specific attention to the selection of lambda in the studies of Gubbins (1983) and Gubbins and Bloxham (1985). It is argued that the selection method is inappropriate and leads to lambda values much larger than those that would result if a reasonable bound on the heat flow at the CMB were assumed.

Backus, George E.↗
NASA NTRS2007

Physically-Retrieving Cloud and Thermodynamic Parameters from Ultraspectral IR Measurements

A physical inversion scheme has been developed, dealing with cloudy as well as cloud-free radiance observed with ultraspectral infrared sounders, to simultaneously retrieve surface, atmospheric thermodynamic, and cloud microphysical parameters. A fast radiative transfer model, which applies to the clouded atmosphere, is used for atmospheric profile and cloud parameter retrieval. A one-dimensional (1-d) variational multi-variable inversion solution is used to improve an iterative background state defined by an eigenvector-regression-retrieval. The solution is iterated in order to account for non-linearity in the 1-d variational solution. It is shown that relatively accurate temperature and moisture retrievals can be achieved below optically thin clouds. For optically thick clouds, accurate temperature and moisture profiles down to cloud top level are obtained. For both optically thin and thick cloud situations, the cloud top height can be retrieved with relatively high accuracy (i.e., error < 1 km). NPOESS Airborne Sounder Testbed Interferometer (NAST-I) retrievals from the Atlantic-THORPEX Regional Campaign are compared with coincident observations obtained from dropsondes and the nadir-pointing Cloud Physics Lidar (CPL). This work was motivated by the need to obtain solutions for atmospheric soundings from infrared radiances observed for every individual field of view, regardless of cloud cover, from future ultraspectral geostationary satellite sounding instruments, such as the Geosynchronous Imaging Fourier Transform Spectrometer (GIFTS) and the Hyperspectral Environmental Suite (HES). However, this retrieval approach can also be applied to the ultraspectral sounding instruments to fly on Polar satellites, such as the Infrared Atmospheric Sounding Interferometer (IASI) on the European MetOp satellite, the Cross-track Infrared Sounder (CrIS) on the NPOESS Preparatory Project and the following NPOESS series of satellites.

Zhou, Daniel K.↗
NASA NTRS2008

Current Sounding Capability From Satellite Meteorological Observation With Ultraspectral Infrared Instruments

Ultraspectral resolution infrared spectral radiance obtained from near nadir observations provide atmospheric, surface, and cloud property information. The intent of the measurement of tropospheric thermodynamic state and trace abundances is the initialization of climate models and the monitoring of air quality. The NPOESS Airborne Sounder Testbed-Interferometer (NAST-I), designed to support the development of future satellite temperature and moisture sounders, aboard high altitude aircraft has been collecting data throughout many field campaigns. An advanced retrieval algorithm developed with NAST-I is now applied to satellite data collected with the Atmospheric InfraRed Sounder (AIRS) on the Aqua satellite launched on 4 May 2002 and the Infrared Atmospheric Sounding Interferometer (IASI) on the MetOp satellite launched on October 19, 2006. These instruments possess an ultra-spectral resolution, for example, both IASI and NAST-I have ~0.25 cm-1 and a spectral coverage from 645 to 2760 cm-1. The retrieval algorithm with a fast radiative transfer model, including cloud effects, is used for atmospheric profile and cloud parameter retrieval. The physical inversion scheme has been developed, dealing with cloudy as well as cloud-free radiance observed with ultraspectral infrared sounders, to simultaneously retrieve surface, atmospheric thermodynamic, and cloud microphysical parameters. A fast radiative transfer model, which applies to the clouded atmosphere, is used for atmospheric profile and cloud parameter retrieval. A one-dimensional (1-d) variational multi-variable inversion solution is used to improve an iterative background state defined by an eigenvector-regression-retrieval. The solution is iterated in order to account for non-linearity in the 1-d variational solution. It is shown that relatively accurate temperature and moisture retrievals can be achieved below optically thin clouds. For optically thick clouds, accurate temperature and moisture profiles down to cloud top level are obtained. For both optically thin and thick cloud situations, the cloud top height can be retrieved with relatively high accuracy (i.e., error less than 1 km). Retrievals of atmospheric soundings, surface properties, and cloud microphysical properties with the AIRS and IASI observations are obtained and presented. These retrievals are further inter-compared with those obtained from airborne FTS system, such as the NPOESS Airborne Sounder Testbed? Interferometer (NAST I), dedicated dropsondes, radiosondes, and ground based Raman Lidar. The capabilities of satellite ultra-spectral sounder such as the AIRS and IASI are investigated. These advanced satellite ultraspectral infrared instruments are now playing an important role in satellite meteorological observation for numerical weather prediction.

Zhou, Daniel K.↗
DOE OSTI2025

Nuclear Data Adjustment for Nonlinear Applications in the OECD/NEA WPNCS SG14 Benchmark -- A Bayesian Inverse UQ-based Approach for Data Assimilation

The Organization for Economic Cooperation and Development (OECD) Working Party on Nuclear Criticality Safety (WPNCS) proposed a benchmark exercise to assess the performance of current nuclear data adjustment techniques applied to nonlinear applications and experiments with low correlation to applications. This work introduces Bayesian Inverse Uncertainty Quantification (IUQ) as a method for nuclear data adjustments in this benchmark, and compares IUQ to the more traditional methods of Generalized Linear Least Squares (GLLS) and Monte Carlo Bayes (MOCABA). Posterior predictions from IUQ showed agreement with GLLS and MOCABA for linear applications. When comparing GLLS, MOCABA, and IUQ posterior predictions to computed model responses using adjusted parameters, we observe that GLLS predictions fail to replicate computed response distributions for nonlinear applications, while MOCABA shows near agreement, and IUQ uses computed model responses directly. We also discuss observations on why experiments with low correlation to applications can be informative to nuclear data adjustments and identify some properties useful in selecting experiments for inclusion in nuclear data adjustment. Performance in this benchmark indicates potential for Bayesian IUQ in nuclear data adjustments.

FOS: Computer and information sciences↗
NASA NTRS1972

Collisional excitation of carbon monoxide in interstellar clouds.

Investigation of different models for the collisional excitation of carbon monoxide (CO) by solving the rate equations for the ten lowest energy levels, including stimulated emission and absorption by 3 K isotropic radiation, spontaneous emission, and collisions. The results are qualitatively applicable to other linear molecules. They show that a wide range of rotational excitation temperatures, including population inversions, can be produced by the various kinetic temperatures and particle densities which are likely to occur in the interstellar regions containing molecules.

Goldsmith, P. F.↗
NASA NTRS1987

Exact finite element method analysis of viscoelastic tapered structures to transient loads

A general method is presented for determining the dynamic torsional/axial response of linear structures composed of either tapered bars or shafts to transient excitations. The method consists of formulating and solving the dynamic problem in the Laplace transform domain by the finite element method and obtaining the response by a numerical inversion of the transformed solution. The derivation of the torsional and axial stiffness matrices is based on the exact solution of the transformed governing equation of motion, and it consequently leads to the exact solution of the problem. The solution permits treatment of the most practical cases of linear tapered bars and shafts, and employs modeling of structures with only one element per member which reduces the number of degrees of freedom involved. The effects of external viscous or internal viscoelastic damping are also taken into account.

Spyrakos, Constantine Chris↗
NASA NTRS2004

An Analytical Calibration Approach for the Polarimetric Airborne C Band Radiometer

Passive microwave remote sensing is sensitive to the quantity and distribution of water in soil and vegetation. During summer 2000, the Microwave Geophysics Group a t the University of Michigan conducted the seventh Radiobrighness Energy Balance Experiment (REBEX-7) over a corn canopy in Michigan. Long time series of brightness temperatures, soil moisture and micrometeorology on the plot were taken. This paper addresses the calibration of the NASA GSFC polarimetric airborne C band microwave radiometer (ACMR) that participated in REBEX-7. These passive polarimeters are typically calibrated using an end-to-end approach based upon a standard artificial target or a well-known geophysical target. Analyzing the major internal functional subsystems offers a different perspective. The primary goal of this approach is to provide a transfer function that not only describes the system in its entire5 but also accounts for the contributions of each subsystem toward the final modified Stokes parameters. This approach does not assume that the radiometric system is linear as it does not take polarization isolation for granted, and it also serves as a realistic instrument simulator, a useful tool for future designs. The ACMR architecture can be partitioned into functional subsystems. The characteristics of each subsystem was extensively measured and the estimated parameters were imported into the overall dosed form system model. Inversion of the model yields a calibration for the modeled Stokes parameters with uncertainties of 0.2 K for the V and H polarizations and 2.4 K for the 3rd and 4th parameters. Application to the full Stokes parameters over a senescent cornfield is presented.

Pham, Hanh↗
NASA NTRS1995

Inversion Of Dynamical Equations For Control Of Attitude

Method of inverting nonlinear equations of rotational dynamics of rigid body used to design feedback control of orientation of body. Applicable to both direction-cosine and quaternion formulations suitable for large-angle maneuvers. Exploiting some apparently little-known properties of direction cosine and quaternion formulations, method leads to equations for model-follower control system that exhibits exactly linear attitude-error dynamics. Quarternion system more robust in responding to large roll-angle commands.

Bach, Ralph↗