Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear stability 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,315 records · Page 73

Semi-Implicit Computation of Fast Modes in a Scheme Integrating Slow Modes by a Leapfrog Method Based on a Selective Implicit Time Filter

Abstract A scheme for integration of atmospheric equations containing terms with differing time scales is developed. The method employs a filtered leapfrog scheme utilizing a fourth-order implicit time filter with one function evaluation per time step to compute slow-propagating phenomena such as advection and rotation. The terms involving fast-propagating modes are handled implicitly with an unconditionally stable method that permits application of larger time steps and faster computations compared to fully explicit treatment. Implementation using explicit and recurrent formulation is provided. Stability analysis demonstrates that the method is conditionally stable for any combination of frequencies involved in the slow and fast terms as they approach the origin. The implicit filter used in the method damps the computational modes without noticeably sacrificing the accuracy of the physical mode. TheO[(Δt 4 )] accuracy for amplitude errors achieved by the implicitly filtered leapfrog is preserved in applications where terms responsible for fast propagation are integrated with a semi-implicit method. Detailed formulation of the method for soundproof nonhydrostatic anelastic equations is provided. Procedures for implementation in global spectral shallow-water models are also given. Examples comparing numerical and analytical solutions for linear gravity waves demonstrate the accuracy of the scheme. The performance is also shown in more practical nonlinear applications, where numerical solutions accomplished by the method are evaluated against those computed from a scheme where the slow terms are handled by the third-order Runge–Kutta scheme. It demonstrates that the method is able to accurately resolve fine-scale dynamics of Kelvin–Helmholtz shear instabilities, the evolution of density current, and nonlinear drifts of twin tropical cyclones.

Meteorology & Atmospheric Sciences↗

The Aeroacoustics of Supersonic Coaxial Jets

Instability waves have been established as the dominant source of mixing noise radiating into the downstream arc of a supersonic jet when the waves have phase velocities that are supersonic relative to ambient conditions. Recent theories for supersonic jet noise have used the concepts of growing and decaying linear instability waves for predicting radiated noise. This analysis is extended to the prediction of noise radiation from supersonic coaxial jets. Since the analysis requires a known mean flow and the coaxial jet mean flow is not described easily in terms of analytic functions, a numerical prediction is made for its development. The Reynolds averaged, compressible, boundary layer equations are solved using a mixing length turbulence model. Empirical correlations are developed for the effects of velocity and temperature ratios and Mach number. Both normal and inverted velocity profile coaxial jets are considered. Comparisons with measurements for both single and coaxial jets show good agreement. The results from mean flow and stability calculations are used to predict the noise radiation from coaxial jets with different operating conditions. Comparisons are made between different coaxial jets and a single equivalent jet with the same total thrust, mass flow, and exit area. Results indicate that normal velocity profile jets can have noise reductions compared to the single equivalent jet. No noise reductions are found for inverted velocity profile jets operated at the minimum noise condition compared to the single equivalent jet. However, it is inferred that changes in area ratio may provide noise reduction benefits for inverted velocity profile jets.

Dahl, Milo D.↗

Robust multi‐loop control of a floating wind turbine

Abstract A principal challenge facing the control of floating offshore wind turbines (FOWTs) is the problem of instability, or “negative damping,” when using blade pitch feedback to control generator speed. This closed‐loop instability can be attributed to non‐minimum phase zeros in the transfer function from blade pitch to generator speed. Standard approaches to improving stability and performance include robust tuning of control gains and introducing multiple feedback loops to respond to platform motion. Combining these approaches is nontrivial because multiple control loops complicate the impact of coupling in the system dynamics. The single‐loop approach to analyzing stability robustness neglects inter‐loop coupling, while a simplistic multi‐loop approach is highly sensitive to dimensional scaling and overestimates the robustness of the single‐loop controller. This work proposes a sensitivity representation that separates some of the natural FOWT dynamic coupling into a parallel feedback loop in the sensitivity function loop to address both of these concerns. The modified robustness measure is used with a simplified linear FOWT model to optimize scheduled multi‐loop control parameters in an automated tuning procedure. This controller is implemented for the 10‐MW Ultraflexible Smart FLoating Offshore Wind Turbine (USFLOWT) and compared against conventional single‐ and multi‐loop controllers tuned using frequency‐domain analysis and high‐fidelity OpenFAST simulations. The multi‐loop robust controller shows the highest overall performance in generator speed regulation and tower load reduction, though consideration of power quality, actuator usage, and other structural loading leads to additional trade‐offs.

Stockhouse, David↗

Nonlinear axisymmetric liquid currents in spherical annuli

A numerical analysis of non-linear axisymmetric viscous flows in spherical annuli of different gap sizes is presented. Only inner sphere was supposed to rotate at a constant angular velocity. The streamlines, lines of constant angular velocity, kinetic energy spectra, and spectra of velocity components are obtained. A total kinetic energy and torque needed to rotate the inner sphere are calculated as functions of Re for different gap sizes. In small-gap annulus nonuniqueness of steady solutions of Navier-Stokes equations is established and regions of different flow regime existences are found. Numerical solutions in a wide-gap annulus and experimental results are used in conclusions about flow stability in the considered range of Re. The comparison of experimental and numerical results shows close qualitative and quantitative agreement.

Astafyeva, N. M.↗

Eigenmode analysis of the sheared-flow Z-pinch

Experiments have demonstrated that a Z-pinch can persist for thousands of times longer than the growth time of global magnetohydrodynamic (MHD) instabilities such as the m=0 sausage and m=1 kink modes. These modes have growth times on the order of ta=a/vi, where vi is the ion thermal speed and a is the pinch radius. Axial flows with duz/dr ≲ vi/a have been measured during the stable period, and the commonly accepted theory is that this amount of shear is sufficient to stabilize these modes as predicted by numerical studies using the ideal MHD equations. However, these studies only consider specific equilibrium profiles that typically have a modest magnitude for the logarithmic pressure gradient, qP≡d ln P/d ln r, and may not represent experimental conditions. Linear stability of the sheared-flow Z-pinch is studied here via a direct eigen-decomposition of the matrix operator obtained from the linear ideal MHD equations. Several equilibrium profiles with a large variation of qP are examined. Considering a practical range of k, 1/3 ≲ ka ≲ 10, it is shown that the shear required to stabilize m=0 modes can be expressed as duz/dr≥Cγ0/(ka)α. Here, γ0=γ0(ka) is the profile-specific growth rate in the absence of shear, which scales approximately with |qP|. Both C and α are profile-specific constants, but C is order unity and α≈1. It is further demonstrated that even a large value of shear, duz/dr=3vi/a, is not sufficient to provide linear stabilization of the m=1 kink mode for all profiles considered. This result is in contrast to the currently accepted theory predicting stabilization at much lower shear, duz/dr=0.1vi/a, and suggests that the experimentally observed stability cannot be explained within the linear ideal-MHD model.

Angus, J. R. (ORCID:0000000314740002)↗

First-principles investigation of thermodynamic, elastic and electronic properties of Al{sub 3}V and Al{sub 3}Nb intermetallics under pressures

The thermodynamic, elastic, and electronic properties of D0{sub 22}-type Al{sub 3}V and Al{sub 3}Nb intermetallics were studied using the first-principle method. The results showed the pressure has profound effects on the structural, mechanical and electronic properties in both Al{sub 3}V and Al{sub 3}Nb. Thermodynamically, the formation enthalpies for Al{sub 3}V and Al{sub 3}Nb were derived, which agreed well with available experimental and theoretical values. Comparably, Al{sub 3}Nb was a more stable phase with the more negative H{sub f} than Al{sub 3}V. Mechanically, the calculated elastic constants showed linearly increasing tendencies, and satisfied the Born's criteria from 0–20 GPa, indicating the mechanically stability of Al{sub 3}V and Al{sub 3}Nb under this pressure range. Further, the mechanical parameters (i.e., bulk modulus (B), shear modulus (G), and Young's modulus (E)) were derived using the Voigt-Reuss-Hill (VRH) method, and in good agreement with available experimental results at the ground state. All these parameters presented the linearly increasing dependences on the external pressure. The B/G ratios and Poisson's ratio indicated that the Al{sub 3}V and Al{sub 3}Nb crystals should exhibit brittle behavior at 0–20 GPa. Additionally, the bulk modulus can be obtained through fitting the Birch-Murnaghan equation (B{sub 0}), computing by VRH method (B{sub H}), and deriving from the elastic theory (B{sub relax}) in both intermetallics. The uniformity of these calculated bulk moduli in each compound exhibited the excellent reliability and self-consistency. In addition, Debye temperature was estimated from the average sound velocity. The Debye temperature showed an increasing dependence on the pressures. Finally, through density of states analysis, Al{sub 3}V and Al{sub 3}Nb were suggested to possess naturally metallic behavior. Under pressures, it was noted that the shapes of peaks and pseudogaps exhibited relative few changes, suggesting Al{sub 3}V and Al{sub 3}Nb has kept structurally stable up to 20 GPa. At zero pressure, Al{sub 3}Nb was considered as a more structurally stable phase with the more number of bonding electrons per atom than Al{sub 3}V. This conclusion was in consistent with the one drawn from the thermodynamic analysis.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Transition In Supersonic Flows With Corners

It is proposed to investigate the effect of sudden change in the slope of the surface on the stability and transition onset point in supersonic and hypersonic boundary layers. Since we can not use the linear stability theory in flows where there are discontinuities we have to solve the full Navier-Stokes equations to determine the amplification of the disturbances across the oblique shock which exists at the corner. We will use the Navier-Stokes code which was developed by Dr. Harold Atkins in the aerodynamic and acoustic methods branch. The code is developed using the higher order ENO scheme and is currently set up for two-dimensional flows. In the first year of our investigation, we will investigate the evolution of small amplitude two-dimensional disturbances and determine the N-factors for the most amplified disturbances. We will perform the analysis for the parameter similar to X-31 model and flight experiment conditions. From these computation we will infer what is the maximum amplification rate possible in quite environments without any tripping devices. In the second part of the investigation, we will perform nonlinear computations to determine what is the minimum amplitudes necessary to cause the transition at a designed location. In the third part, we have to determine what kind of trip devices necessary to excite the required disturbances. The principal investigator of this project is Dr. P. Balakumar, Associate Professor of Department of Aerospace Engineering, Old Dominion University. We will actively collaborate with research scientists in the Aero-Thermodynamics Branch at NASA/Langley Research Center. A half-time graduate student will also participate in the project.

Balakumar, P.↗

Microbial extracellular polysaccharide production and aggregate stability controlled by switchgrass ( Panicum virgatum ) root biomass and soil water potential

Deep-rooting perennial grasses are promising feedstocks for biofuel production, especially in marginal soils lacking organic material, nutrients, and/or that experience significant water stress. Perennial grass roots influence surrounding soil conditions and microbial activities, and produce extracellular polymeric substances (EPS) composed primarily of extracellular polysaccharides (EPSac). These polymers can alleviate microbial moisture and nutrient stress, and enhance soil characteristics through improved water retention and aggregate stability—which may in turn enhance carbon persistence. Here, we used a 13 CO 2 greenhouse tracer experiment to examine the effect of switchgrass cultivation on EPSac production and origin in a marginal soil with five fertilization/water treatments (control, +N, +NP, +P, low water), and compared these results with measurements of field soils collected after long-term switchgrass cultivation. Soils with added nitrogen and phosphorus (+NP) had the highest root biomass, EPSac and percentage of water-stable soil aggregates. Multiple linear regression analyses revealed that root biomass and soil water potential were important determinants of soil EPSac production, potentially by controlling carbon supply and diurnal changes in moisture stress. Path analysis showed that soil aggregation was positively correlated with bulk soil EPSac content and also regulated by soil water potential. High mannose content indicated the majority of EPSac was of microbial origin and 13 CO 2 labeling indicated that 0.18% of newly fixed plant carbon was incorporated into EPSac. Analysis of field soils suggests that EPSac is significantly enhanced after long-term switchgrass cultivation. Taken as a whole, our greenhouse and field results demonstrate that switchgrass cultivation can promote microbial production of EPSac, providing a mechanism to enhance aggregation in marginal soils.

13c labeling↗

Insights into Tetravalent Np Speciation in HNO 3 through Spectroelectrochemistry and Multivariate Analysis

In situ optical spectroscopy, spectropotentiometry, and multivariate analysis were applied to the Np(IV) nitrate system to better understand speciation and quantify HNO 3 concentration. Thin-layer spectropotentiometry, or spectroelectrochemistry, was leveraged to isolate and stabilize Np(IV) without compromising the solution conditions and generate representative Vis-NIR absorption spectra from 0.5 to 10 M HNO 3 and benchmark the corresponding Np(IV) molar absorptivity coefficients. Spectra were described with principal component analysis (PCA) to identify the purest Np(IV) absorbance spectra among other oxidation states [e.g., Np(V/VI)] at each acid concentration and then to identify the primary sources of variance within each Np(IV) spectrum with respect to Np(IV) nitrate complexes. Then, partial least-squares regression (PLSR) and support vector regression (SVR) models were built to predict HNO 3 concentration from the Np(IV) spectral data. The nonlinear SVR model outperformed the linear PLSR model for the HNO 3 concentration predictions. Finally, the inclusion of spectra collected in edge and center point HNO 3 concentrations in the calibration set was determined to be crucial for producing models with strong predictive capabilities. The multivariate approach used in this study makes it possible to quantify HNO 3 concentration solely based on Np(IV) absorption spectra, which is essential to quantifying processing streams in various online monitoring applications.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Impact of impurities on peeling–ballooning modes and turbulence in tokamak plasmas

This study investigates the impact of various impurity species on peeling–ballooning (PB) modes and microturbulence in tokamak plasmas through the extension of traditional two-fluid and gyro-landau-fluid (GLF) models. By incorporating finite Larmor radius (FLR) effects, the analysis provides a comprehensive understanding of impurity-driven impact and its interaction with plasma turbulence. Depending on charge state and local plasma conditions, heavy impurities may exhibit gyro-radii larger than those of main ions, which are captured in the extended GLF model presented. Following the presentation of modified two-fluid equations incorporating impurity effects, we systematically analyze the distinctions between impurity and main ion dynamics and their resultant feedback mechanisms on plasma behavior. Derivation of the linear dispersion relation enables quantification of impurity-mediated modifications to: plasma vorticity, diamagnetic drift and gyroviscous effects, electron Hall physics, and FLR effects. BOUT++ – based linear simulations corroborate this formalism, demonstrating systematic stabilization of PB modes upon impurity seeding. And then operational implications for practical impurity control strategies in tokamak devices are proposed. The results underscore the necessity of impurity management to maintain stability and optimize plasma confinement, with specific focus on how FLR effects contribute to transport dynamics. This work paves the way for enhanced modeling and simulation efforts, supporting the development of strategies to control impurity-induced turbulence and improve overall reactor performance.

BOUT++ simulation↗

Reduced Order Modeling for Prediction and Control of Large-Scale Systems.

This report describes work performed from June 2012 through May 2014 as a part of a Sandia Early Career Laboratory Directed Research and Development (LDRD) project led by the first author. The objective of the project is to investigate methods for building stable and efficient proper orthogonal decomposition (POD)/Galerkin reduced order models (ROMs): models derived from a sequence of high-fidelity simulations but having a much lower computational cost. Since they are, by construction, small and fast, ROMs can enable real-time simulations of complex systems for onthe- spot analysis, control and decision-making in the presence of uncertainty. Of particular interest to Sandia is the use of ROMs for the quantification of the compressible captive-carry environment, simulated for the design and qualification of nuclear weapons systems. It is an unfortunate reality that many ROM techniques are computationally intractable or lack an a priori stability guarantee for compressible flows. For this reason, this LDRD project focuses on the development of techniques for building provably stable projection-based ROMs. Model reduction approaches based on continuous as well as discrete projection are considered. In the first part of this report, an approach for building energy-stable Galerkin ROMs for linear hyperbolic or incompletely parabolic systems of partial differential equations (PDEs) using continuous projection is developed. The key idea is to apply a transformation induced by the Lyapunov function for the system, and to build the ROM in the transformed variables. It is shown that, for many PDE systems including the linearized compressible Euler and linearized compressible Navier-Stokes equations, the desired transformation is induced by a special inner product, termed the “symmetry inner product”. Attention is then turned to nonlinear conservation laws. A new transformation and corresponding energy-based inner product for the full nonlinear compressible Navier-Stokes equations is derived, and it is demonstrated that if a Galerkin ROM is constructed in this inner product, the ROM system energy will be bounded in a way that is consistent with the behavior of the exact solution to these PDEs, i.e., the ROM will be energy-stable. The viability of the linear as well as nonlinear continuous projection model reduction approaches developed as a part of this project is evaluated on several test cases, including the cavity configuration of interest in the targeted application area. In the second part of this report, some POD/Galerkin approaches for building stable ROMs using discrete projection are explored. It is shown that, for generic linear time-invariant (LTI) systems, a discrete counterpart of the continuous symmetry inner product is a weighted L2 inner product obtained by solving a Lyapunov equation. This inner product was first proposed by Rowley et al., and is termed herein the “Lyapunov inner product“. Comparisons between the symmetry inner product and the Lyapunov inner product are made, and the performance of ROMs constructed using these inner products is evaluated on several benchmark test cases. Also in the second part of this report, a new ROM stabilization approach, termed “ROM stabilization via optimization-based eigenvalue reassignment“, is developed for generic LTI systems. At the heart of this method is a constrained nonlinear least-squares optimization problem that is formulated and solved numerically to ensure accuracy of the stabilized ROM. Numerical studies reveal that the optimization problem is computationally inexpensive to solve, and that the new stabilization approach delivers ROMs that are stable as well as accurate. Summaries of “lessons learned“ and perspectives for future work motivated by this LDRD project are provided at the end of each of the two main chapters.

97 MATHEMATICS AND COMPUTING↗

Aeromechanical stability of helicopters with a bearingless main rotor. Part 1: Equations of motion

Equations of motion for a coupled rotor-body system were derived for the purpose of studying air and ground resonance characteristics of helicopters that have bearingless main rotors. For the fuselage, only four rigid body degrees of freedom are considered; longitudinal and lateral translations, pitch, and roll. The rotor is assumed to consist of three or more rigid blades. Each blade is joined to the hub by means of a flexible beam segment (flexbeam or strap). Pitch change is accomplished by twisting the flexbeam with the pitch-control system, the characteristics of which are variable. Thus, the analysis is capable of implicitly treating aeroelastic couplings generated by the flexbeam elastic deflections, the pitch-control system, and the angular offsets of the blade and flexbeam. The linearized equations are written in the nonrotating system retaining only the cyclic rotor modes; thus, they comprise a system of homogeneous ordinary differential equations with constant coefficients. All contributions to the linearized perturbation equations from inertia, gravity, quasi-steady aerodynamics, and the flexbeam equilibrium deflections are retained exactly.

Hodges, D. H.↗

Dynamic analysis of a magnetic bearing system with flux control

Using measured values of two-dimensional forces in a magnetic actuator, equations of motion for an active magnetic bearing are presented. The presence of geometric coupling between coordinate directions causes the equations of motion to be nonlinear. Two methods are used to examine the unbalance response of the system: simulation by direct integration in time; and determination of approximate steady state solutions by harmonic balance. For relatively large values of the derivative control coefficient, the system behaves in an essentially linear manner, but for lower values of this parameter, or for higher values of the coupling coefficient, the response shows a split of amplitudes in the two principal directions. This bifurcation is sensitive to initial conditions. The harmonic balance solution shows that the separation of amplitudes actually corresponds to a change in stability of multiple coexisting solutions.

Knight, Josiah↗

Dehydration of Methyl Lactate on Alkali Cation-Exchanged Faujasite: Effects of Metal Cation Identity and Water Pressure

Turnover rates for the catalytic dehydration of methyl lactate (ML) over ion-exchanged faujasite (FAU) catalysts depend on the identity of alkali metal cations (Na + , K + , Cs + ) and local solvation effects. Analysis of rate measurements and in situ infrared spectroscopy gives evidence that the reaction involves kinetically relevant dissociation of adsorbed methyl lactate upon alkali metal cations. This process involves concerted methyl transfer to the surface and dissociation of the alkali metal from the framework, which occurs at cationic active sites that remain predominantly unoccupied under relevant conditions (0.5–10 kPa ML, 0.5–15 kPa H 2 O, 563–583 K). Despite the mechanistic similarities, apparent activation enthalpies (ΔH app ‡ ) decrease linearly (47 kJ mol –1 from Na + to Cs + ) with ionization energy and cationic radius, and apparent activation entropies (ΔS app ‡ ) decrease 74 J mol –1 K –1 . These trends reflect electrostatic interactions that stabilize the cations to the anionic sites on the zeolite: stronger association between these charges leads to increasingly endothermic processes to displace the alkali metal to form a cationic methoxy and an intrapore metal lactate intermediate. Water physisorption measurements suggest alkali metal ions bind superstoichiometric quantities of water within FAU pores, and in situ infrared spectra suggest the concerted adsorption of ML requires reorganization of this water. Consequently, these processes introduce entropic gains that partially offset entropy losses associated with ML adsorption. Hence, turnover rates differ only by a factor of 2 among Na-, K-, and Cs-FAU at 573 K (ΔΔG app ‡ = 5 kJ mol –1 ). These findings demonstrate the interplay of alkali metal ions with zeolite active sites and intrapore water clusters for ML dehydration, indicating that these interactions can be leveraged to deliver optimal performance under different reaction conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of tearing instability in a current sheet forming by sheared incompressible flow

Sweet–Parker current sheets in high Lundquist number plasmas are unstable to tearing, suggesting they will not form in physical systems. Understanding magnetic reconnection thus requires study of the stability of a current sheet as it forms. Formation can occur due to sheared, sub-Alfvénic incompressible flows which narrow the sheet. Standard tearing theory is not immediately applicable to such forming sheets for two reasons: first, because the flow introduces terms not present in the standard calculation; second, because the changing equilibrium introduces time dependence to terms which are constant in the standard calculation, complicating the formulation of an eigenvalue problem. This paper adapts standard tearing mode analysis to confront these challenges. In an initial phase when any perturbations are primarily governed by ideal magnetohydrodynamics, a coordinate transformation reveals that the flow compresses and stretches perturbations. Here, a multiple scale formulation describes how linear tearing mode theory can be applied to an equilibrium changing under flow, showing that the flow affects the separable exponential growth only implicitly, by making the standard scalings time dependent. In the nonlinear Rutherford stage, the coordinate transformation shows that standard theory can be adapted by adding to the stationary rates time dependence and an additional term due to the strengthening equilibrium magnetic field. Overall, this understanding supports the use of flow-free scalings with slight modifications to study tearing in a forming sheet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

TU-F-CAMPUS-I-05: Semi-Automated, Open Source MRI Quality Assurance and Quality Control Program for Multi-Unit Institution

Purpose: Phantom measurements allow for the performance of magnetic resonance (MR) systems to be evaluated. Association of Physicists in Medicine (AAPM) Report No. 100 Acceptance Testing and Quality Assurance Procedures for MR Imaging Facilities, American College of Radiology (ACR) MR Accreditation Program MR phantom testing, and ACR MRI quality control (QC) program documents help to outline specific tests for establishing system performance baselines as well as system stability over time. Analyzing and processing tests from multiple systems can be time-consuming for medical physicists. Besides determining whether tests are within predetermined limits or criteria, monitoring longitudinal trends can also help prevent costly downtime of systems during clinical operation. In this work, a semi-automated QC program was developed to analyze and record measurements in a database that allowed for easy access to historical data. Methods: Image analysis was performed on 27 different MR systems of 1.5T and 3.0T field strengths from GE and Siemens manufacturers. Recommended measurements involved the ACR MRI Accreditation Phantom, spherical homogenous phantoms, and a phantom with an uniform hole pattern. Measurements assessed geometric accuracy and linearity, position accuracy, image uniformity, signal, noise, ghosting, transmit gain, center frequency, and magnetic field drift. The program was designed with open source tools, employing Linux, Apache, MySQL database and Python programming language for the front and backend. Results: Processing time for each image is <2 seconds. Figures are produced to show regions of interests (ROIs) for analysis. Historical data can be reviewed to compare previous year data and to inspect for trends. Conclusion: A MRI quality assurance and QC program is necessary for maintaining high quality, ACR MRI Accredited MR programs. A reviewable database of phantom measurements assists medical physicists with processing and monitoring of large datasets. Longitudinal data can reveal trends that although are within passing criteria indicate underlying system issues.

60 APPLIED LIFE SCIENCES↗

General Rotorcraft Aeromechanical Stability Program (GRASP): Theory manual

The general rotorcraft aeromechanical stability program (GRASP) was developed to calculate aeroelastic stability for rotorcraft in hovering flight, vertical flight, and ground contact conditions. GRASP is described in terms of its capabilities and its philosophy of modeling. The equations of motion that govern the physical system are described, as well as the analytical approximations used to derive them. The equations include the kinematical equation, the element equations, and the constraint equations. In addition, the solution procedures used by GRASP are described. GRASP is capable of treating the nonlinear static and linearized dynamic behavior of structures represented by arbitrary collections of rigid-body and beam elements. These elements may be connected in an arbitrary fashion, and are permitted to have large relative motions. The main limitation of this analysis is that periodic coefficient effects are not treated, restricting rotorcraft flight conditions to hover, axial flight, and ground contact. Instead of following the methods employed in other rotorcraft programs. GRASP is designed to be a hybrid of the finite-element method and the multibody methods used in spacecraft analysis. GRASP differs from traditional finite-element programs by allowing multiple levels of substructure in which the substructures can move and/or rotate relative to others with no small-angle approximations. This capability facilitates the modeling of rotorcraft structures, including the rotating/nonrotating interface and the details of the blade/root kinematics for various types. GRASP differs from traditional multibody programs by considering aeroelastic effects, including inflow dynamics (simple unsteady aerodynamics) and nonlinear aerodynamic coefficients.

Hodges, Dewey H.↗

Aeroservoelastic Model Validation and Test Data Analysis of the F/A-18 Active Aeroelastic Wing

Model validation and flight test data analysis require careful consideration of the effects of uncertainty, noise, and nonlinearity. Uncertainty prevails in the data analysis techniques and results in a composite model uncertainty from unmodeled dynamics, assumptions and mechanics of the estimation procedures, noise, and nonlinearity. A fundamental requirement for reliable and robust model development is an attempt to account for each of these sources of error, in particular, for model validation, robust stability prediction, and flight control system development. This paper is concerned with data processing procedures for uncertainty reduction in model validation for stability estimation and nonlinear identification. F/A-18 Active Aeroelastic Wing (AAW) aircraft data is used to demonstrate signal representation effects on uncertain model development, stability estimation, and nonlinear identification. Data is decomposed using adaptive orthonormal best-basis and wavelet-basis signal decompositions for signal denoising into linear and nonlinear identification algorithms. Nonlinear identification from a wavelet-based Volterra kernel procedure is used to extract nonlinear dynamics from aeroelastic responses, and to assist model development and uncertainty reduction for model validation and stability prediction by removing a class of nonlinearity from the uncertainty.

Brenner, Martin J.↗