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,495 records · Page 83

Non-autoregressive time-series methods for stable parametric reduced-order models

Advection-dominated dynamical systems, characterized by partial differential equations, are found in applications ranging from weather forecasting to engineering design where accuracy and robustness are crucial. There has been significant interest in the use of techniques borrowed from machine learning to reduce the computational expense and/or improve the accuracy of predictions for these systems. These rely on the identification of a basis that reduces the dimensionality of the problem and the subsequent use of time series and sequential learning methods to forecast the evolution of the reduced state. Often, however, machine-learned predictions after reduced-basis projection are plagued by issues of stability stemming from incomplete capture of multiscale processes as well as due to error growth for long forecast durations. To address these issues, we have developed a non-autoregressive time series approach for predicting linear reduced-basis time histories of forward models. In particular, we demonstrate that non-autoregressive counterparts of sequential learning methods such as long short-term memory (LSTM) considerably improve the stability of machine-learned reduced-order models. Further, we evaluate our approach on the inviscid shallow water equations and show that a non-autoregressive variant of the standard LSTM approach that is bidirectional in the principal component directions obtains the best accuracy for recreating the nonlinear dynamics of partial observations. Moreover-and critical for many applications of these surrogates-inference times are reduced by three orders of magnitude using our approach, compared with both the equation-based Galerkin projection method and the standard LSTM approach.

97 MATHEMATICS AND COMPUTING↗

A discontinuous piecewise polynomial generalized moving least squares scheme for robust finite element analysis on arbitrary grids

A variational approach is developed with a meshless discretization to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. Here, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality, which we call discontinuous piecewise polynomial generalized moving least squares (DPP-GMLS). At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework, providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5 x improvement in accuracy and allows for taking an 8 x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method.

97 MATHEMATICS AND COMPUTING↗

WARP3D-Release 10.8: Dynamic Nonlinear Analysis of Solids using a Preconditioned Conjugate Gradient Software Architecture

This report describes theoretical background material and commands necessary to use the WARP3D finite element code. WARP3D is under continuing development as a research code for the solution of very large-scale, 3-D solid models subjected to static and dynamic loads. Specific features in the code oriented toward the investigation of ductile fracture in metals include a robust finite strain formulation, a general J-integral computation facility (with inertia, face loading), an element extinction facility to model crack growth, nonlinear material models including viscoplastic effects, and the Gurson-Tver-gaard dilatant plasticity model for void growth. The nonlinear, dynamic equilibrium equations are solved using an incremental-iterative, implicit formulation with full Newton iterations to eliminate residual nodal forces. The history integration of the nonlinear equations of motion is accomplished with Newmarks Beta method. A central feature of WARP3D involves the use of a linear-preconditioned conjugate gradient (LPCG) solver implemented in an element-by-element format to replace a conventional direct linear equation solver. This software architecture dramatically reduces both the memory requirements and CPU time for very large, nonlinear solid models since formation of the assembled (dynamic) stiffness matrix is avoided. Analyses thus exhibit the numerical stability for large time (load) steps provided by the implicit formulation coupled with the low memory requirements characteristic of an explicit code. In addition to the much lower memory requirements of the LPCG solver, the CPU time required for solution of the linear equations during each Newton iteration is generally one-half or less of the CPU time required for a traditional direct solver. All other computational aspects of the code (element stiffnesses, element strains, stress updating, element internal forces) are implemented in the element-by- element, blocked architecture. This greatly improves vectorization of the code on uni-processor hardware and enables straightforward parallel-vector processing of element blocks on multi-processor hardware.

Koppenhoefer, Kyle C.↗

Observability Analysis of a Power System Stochastic Dynamic Model Using a Derivative-Free Approach

Serving as a prerequisite to power system dynamic state estimation, the observability analysis of a power system dynamic model has recently attracted the attention of many power engineers. However, because this model is typically nonlinear and large-scale, the analysis of its observability is a challenge to the traditional derivative-based methods. Indeed, the linear-approximation-based approach may provide unreliable results while the nonlinear-technique-based approach inevitably faces extremely complicated derivations. Furthermore, because power systems are intrinsically stochastic, the traditional deterministic approaches may lead to inaccurate observability analyses. In this work, facing these challenges, we propose a novel polynomial-chaos-based derivative-free observability analysis approach that not only is free of any linear approximations, but also accounts for the stochasticity of the dynamic model while bringing a low implementation complexity. Furthermore, this approach enables us to quantify the degree of observability of a stochastic model, what conventional deterministic methods cannot do. The excellent performance of the proposed method has been demonstrated by performing extensive simulations using a synchronous generator model with IEEE-DC1A exciter and the TGOV1 turbine governor.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Describing Function Analysis of Transformer Magnetizing Inductance for Direct Power Control of Back-to-Back Modular Multilevel Converters with Advanced Grid Support

This paper provides a detailed investigation into the application of describing function-based analysis for assessing transformer magnetizing inductance and its impact on system performance. The focus is on a back-to-back modular multilevel converter architecture, designed to interconnect systems operating at different frequencies. The study explores the implementation of a Direct Power Control strategy, examining its effects on transformer magnetizing inductance saturation and offering effective mitigation techniques. Furthermore, the integration of advanced grid support functionalities is highlighted, demonstrating how these enhancements bolster the converter's ability to improve grid stability and power quality, positioning it as a robust solution for modern power systems. The proposed approach is validated through extensive computer simulations based MAT LAB/Simulink domain, supported by significant case study results, confirming its practical effectiveness.

back-to-back modular multilevel converters (B2B- M↗

High-order gradients with the shifted boundary method: An embedded enriched mixed formulation for elliptic PDEs

Here, we propose an extension of the embedded boundary method known as “shifted boundary method” to elliptic diffusion equations in mixed form (e.g., Darcy flow, heat diffusion problems with rough coefficients, etc.). Our aim is to obtain an improved formulation that, for linear finite elements, is at least second-order accurate for both flux and primary variable, when either Dirichlet or Neumann boundary conditions are applied. Following previous work of Nishikawa and Mazaheri in the context of residual distribution methods, we consider the mixed form of the diffusion equation (i.e., with Darcy-type operators), and introduce an enrichment of the primary variable. This enrichment is obtained exploiting the relation between the primary variable and the flux variable, which is explicitly available at nodes in the mixed formulation. The proposed enrichment mimics a formally quadratic pressure approximation, although only nodal unknowns are stored, similar to a linear finite element approximation. We consider both continuous and discontinuous finite element approximations and present two approaches: a non-symmetric enrichment, which, as in the original references, only improves the consistency of the overall method; and a symmetric enrichment, which enables a full error analysis in the classical finite element context. Combined with the shifted boundary method, these two approaches are extended to high-order embedded computations, and enable the approximation of both primary and flux (gradient) variables with second-order accuracy, independently on the type of boundary conditions applied. We also show that the primary variable is third-order accurate, when pure Dirichlet boundary conditions are embedded.

97 MATHEMATICS AND COMPUTING↗

Hydrogen bond arrangements in (H 2 O) 20, 24, 28 clathrate hydrate cages: Optimization and many-body analysis

Here we provide a detailed study of hydrogen bonding arrangements, relative stability, residual entropy, and an analysis of the many-body effects in the (H 2 O) 20 (D-cage), (H 2 O) 24 (T-cage), and (H 2 O) 28 (H-cage) hollow cages making up structures I (sI) and II (sII) of clathrate hydrate lattices. Based on the enumeration of the possible hydrogen bonding networks for a fixed oxygen atom scaffold, the residual entropy (S 0 ) of these three gas phase cages was estimated at 0.754 82, 0.754 44, and 0.754 17 · Nk b , where N is the number of molecules and k b is Boltzmann’s constant. A previously identified descriptor of enhanced stability based on the relative arrangement and connectivity of nearest-neighbor fragments on the polyhedral water cluster [strong-weak-effective-bond model] also applies to the larger hollow cages. The three cages contain a maximum of 7, 9, and 11 such preferable arrangements of trans nearest dimer pairs with one “free” OH bond on the donor molecule (t1d dimers). The Many-Body Expansion (MBE) up to the 4-body suggests that the many-body terms vary nearly linearly with the cluster binding energy. Using a hierarchical approach of screening the relative stability of networks starting from optimizations with the TIP4P, TTM2.1-F, and MB-pol classical potentials, subsequently refining at more accurate levels of electronic structure theory (DFT and MP2), and finally correcting for zero-point energy, we were able to identify a group of four low-lying isomers of the (H 2 O) 24 T-cage, two of which are antisymmetric and the other two form a pair of antipode configurations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonlinear fishbone dynamics in spherical tokamaks

Linear and nonlinear kinetic-MHD hybrid simulations have been carried out to investigate linear stability and nonlinear dynamics of beam-driven fishbone instability in spherical tokamak plasmas. Realistic NSTX parameters with finite toroidal rotation were used. The results show that the fishbone is driven by both trapped and passing particles. The instability drive of passing particles is comparable to that of trapped particles in the linear regime. The effects of rotation are destabilizing and a new region of instability appears at higher q min (>1.5) values, q min being the minimum of safety factor profile. In the nonlinear regime, the mode saturates due to flattening of beam ion distribution, and this persists after initial saturation while mode frequency chirps down in such a way that the resonant trapped particles move out radially and keep in resonance with the mode. Correspondingly, the flattening region of beam ion distribution expands radially outward. A substantial fraction of initially non-resonant trapped particles become resonant around the time of mode saturation and keep in resonance with the mode as frequency chirps down. On the other hand, the fraction of resonant passing particles is significantly smaller than that of trapped particles. Our analysis shows that trapped particles provide the main drive to the mode in the nonlinear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structural and Spectroscopic Characterization of Plutonium and Other Tetravalent Metals Complexed to a Keggin Ion

Here, we report the isolation of the first plutonium(IV) complex with a Keggin ion chelator: Cs 20 [Pu(PW 11 O 39 ) 2 ] 2 ·13H 2 O. Single crystal XRD and solid-state UV–vis absorbance analysis demonstrate the stabilization of Pu 4+ by the Keggin ligand. The unit cell contains two [Pu(PW 11 O 39 ) 2 ] 10– complexes (Pu(PW 11 ) 2 ) bridged by Cs + . Raman and 31 P NMR spectra of Pu(PW 11 ) 2 are consistent with the analogous Zr 4+ , Hf 4+ , Ce 4+ , and Th 4+ complexes. The Pu–O bond distances at the two Pu sites are 2.35(3) and 2.34(3) Å, matching the value extrapolated from the bonding trend built with the other 8-coordinated tetravalent cations. However, the long-range arrangement of the Pu(PW 11 ) 2 complexes within the lattice is unique in the series of M IV (PW 11 ) 2 compounds: pairs of Pu(PW 11 ) 2 are organized perpendicular to each other. Based on solution-state UV–visible absorbance, small-angle X-ray scattering (SAXS), and 31 P NMR, the tetravalent cations quantitatively form the 1:2 species in solution ([Pu(PW 11 O 39 ) 2 ] 10– (aq) ) and no 1:1 species ([Pu(PW 11 O 39 )(H 2 O) x ] 3– (aq) ). Finally, a linear correlation exists between the metal–oxygen distances in the M IV (PW 11 ) 2 compounds and the corresponding metal dioxides, allowing for extrapolation for Pa 4+ , Am 4+ , and Bk 4+ . The results indicate that our microscale POM approach represents a viable pathway to probe properties of rare actinide ions in discrete molecules, beyond the traditional oxide extended solids.

and nuclear chemistry↗

Nonlinear fishbone dynamics in spherical tokamaks

Linear and nonlinear kinetic-MHD hybrid simulations have been carried out to investigate linear stability and nonlinear dynamics of beam-driven fishbone instability in spherical tokamak plasmas. Realistic NSTX parameters with finite toroidal rotation were used. Our results show that the fishbone is driven by both trapped and passing particles. The instability drive of passing particles is comparable to that of trapped particles in the linear regime. The effects of rotation are destabilizing and a new region of instability appears at higher q min (>1.5) values, q min being the minimum of safety factor profile. In the nonlinear regime, the mode saturates due to flattening of beam ion distribution, and this persists after initial saturation while mode frequency chirps down in such a way that the resonant trapped particles move out radially and keep in resonance with the mode. Correspondingly, the flattening region of beam ion distribution expands radially outward. Furthermore, a substantial fraction of initially non-resonant trapped particles become resonant around the time of mode saturation and keep in resonance with the mode as frequency chirps down. On the other hand, the fraction of resonant passing particles is significantly smaller than that of trapped particles. Finally, our analysis shows that trapped particles provide the main drive to the mode in the nonlinear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

New Material for Surface-Enhanced Raman Spectroscopy

A chemical method of synthesis and application of coating materials that are especially suitable for surface-enhanced Raman spectroscopy (SERS) has been developed. The purpose of this development is to facilitate the utilization of the inherently high sensitivity of SERS to detect chemicals of interest (analytes) in trace amounts, without need for lengthy sample preparation. Up to now, the use of SERS has not become routine because the methods available have not been able to reproduce sampling conditions and provide quantitative measurements. In contrast, the coating materials of the present method enable analysis with minimum preparation of samples, and SERS measurements made using these materials are reproducible and reversible. Moreover, unlike in methods investigated in prior efforts to implement SERS, sampling is not restricted to such specific environments as electrolytes or specific solvents. The coating materials of this method are porous glasses, formed in sol-gel processes, that contain small particles of gold or silver metal. Materials of this type can be applied to the sample-contact surfaces of a variety of sampling and sensing devices, including glass slides, glass vials, fiber-optic probes, and glass tubes. Glass vials with their insides coated according to this method are particularly convenient for SERS to detect trace chemicals in solutions: One simply puts a sample solution containing the analyte(s) into a vial, then puts the vial into a Raman spectrometer for analysis. The chemical ingredients and the physical conditions of the sol-gel process have been selected so that the porous glass formed incorporates particles of the desired metal with size(s) to match the wavelength(s) of the SERS excitation laser in order to optimize the generation of surface plasmons. The ingredients and processing conditions have further been chosen to tailor the porosity and polarity of the glass to optimize the sample flow and the interaction between the analyte(s) and the plasmon field that generates Raman photons. The porous silica network of a sol-gel glass creates a unique environment for stabilizing SERS-active metal particles. Relative to other material structures that could be considered for SERS, the porous silica network offers higher specific surface area and thus greater interaction between analyte molecules and metal particles. Efforts to perform SERS measurements with the help of sampling devices coated by this method have been successful. In tests, numerous organic and inorganic chemicals were analyzed in several solvents, including water. The results of the tests indicate that the SERS measurements were reproducible within 10 percent and linear over five orders of magnitude. One measure of the limits of detectability of chemicals in these tests was found to be a concentration of 300 parts per billion. Further development may eventually make it possible to realize the full potential sensitivity of SERS for detecting some analytes in quantities as small as a single molecule.

Farquharson, Stuart↗

Multi-scale, Multi-disciplinary, and Multi-agent Explainable AI with Koopman-Undergirded Learning, Prediction, and Analysis (M3EA KULPA) (Project Closeout Report)

The goal of this project was to develop and use domain-aware machine learning formulations, based on the Koopman Operator (KO), for modelling multi-scale, multi-disciplinary (e.g., multi-physics), and/or multi-agent systems. The project developed these formulations for the following cases: • Systems with dynamics at two separate time scales, • Systems with a bi-level hierarchical control structure, • Systems with bi-level hierarchical control and dynamics at two separate time scales (the lower level controls operating at the faster time scale), and • Systems with n separate but interacting agents/disciplines (with/without control, respectively); the controls for each agent could include bi-level hierarchical control and dynamics at two separate time scales as described above. The project then defined a set of dynamical systems consisting of different nonlinear oscillators that could be used to test these different formulations and then subsequently learned the KO models for those systems. With the KO models, we were able to do the following: • Quantify system stability, including both long-term and transient behavior, • Quantify the effects of feedbacks between the different time scales and agents/disciplines in terms of those feedbacks’ effects on system stability, • Replace a standard Proportional-Integral (PI) control in the hierarchical control structure with a KO-based Linear-Quadratic Regular (LQR), a form of optimal control, • Calculate optimal supervisory control policies a) with and without time scale separated dynamics at the lower level control levels and b) with both PI and KO-based LQR lower level control policies, and • Calculate dynamic Nash equilibria for multi-agent systems where each agent makes its own control decisions.

97 MATHEMATICS AND COMPUTING↗

SU-F-I-70: Investigation of Gafchromic EBT3 Film Energy Dependence Using Proton, Photon, and Electron Beams

Purpose: To investigate the energy dependence of Gafchromic EBT3 film over a range of clinically used proton, photon and electron energies. Methods: Proton beam energies of 117 and 204 MeV, corresponding respectively to ranges in water of 10 cm and 27 cm from a Mevion S250 double scatter system unit were used. Electron energies of 6 and 20 MeV and photon energies of 6 and 18 MV from a Varian Clinac 21EX Linac were used. Two pieces of film (5×5 cm{sup 2}) were irradiated sequentially for doses of 100, 500, and 1000 cGy for all energies and modalities. Films were placed on the central beam axis for a 10×10 cm{sup 2} field size in the middle of spread out Bragg peak (SOBP) for proton and in respective dmax for photon and electron energies. Films were scanned on a flatbed Epson Expression 10000 XL scanner on the central region of the scanning window using 48-bit, 300 dpi, and landscape orientation after 48 hours post-irradiation of film to account for optical density (OD) stabilization. Film analysis of the red channel was performed using ImageJ 1.48v (National Institutes of Health). Results: The energy dependency of EBT3 among all energies and modalities for all doses studied was small within measurement uncertainties (1σ = ± 4.1%). The mean net OD in red channel for films receiving the same dose in the same energy modality had standard deviations within 0.9% for photons, 4.9% for electrons and 1.8% for protons. It was observed that film pieces were activated during proton irradiation, e.g., 7 mR/hr at surface after 30 minutes of irradiation, lasting for 2 hours post irradiation. Conclusion: EBT3 energy dependency was evaluated for clinically used proton, photon, and electron energies. The film self-activation may have contributed to fog and negligible dose.

60 APPLIED LIFE SCIENCES↗

Discontinuous diffusion synthetic acceleration for S{sub n} transport on 2D arbitrary polygonal meshes

In this paper, a Diffusion Synthetic Acceleration (DSA) technique applied to the S{sub n} radiation transport equation is developed using Piece-Wise Linear Discontinuous (PWLD) finite elements on arbitrary polygonal grids. The discretization of the DSA equations employs an Interior Penalty technique, as is classically done for the stabilization of the diffusion equation using discontinuous finite element approximations. The penalty method yields a system of linear equations that is Symmetric Positive Definite (SPD). Thus, solution techniques such as Preconditioned Conjugate Gradient (PCG) can be effectively employed. Algebraic MultiGrid (AMG) and Symmetric Gauss–Seidel (SGS) are employed as conjugate gradient preconditioners for the DSA system. AMG is shown to be significantly more efficient than SGS. Fourier analyses are carried out and we show that this discontinuous finite element DSA scheme is always stable and effective at reducing the spectral radius for iterative transport solves, even for grids with high-aspect ratio cells. Numerical results are presented for different grid types: quadrilateral, hexagonal, and polygonal grids as well as grids with local mesh adaptivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Energetics of oxidation and formation of uranium mononitride

Uranium mononitride (UN) is an advanced nuclear fuel currently being considered for use in several generation IV fast and thermal neutron spectrum core designs, with additional applications to thermal and electric nuclear propulsion reactors. Here, to better understand the thermal behavior and thermodynamic stability of UN, we investigated the bulk thermal oxidation process and thermochemical reactions, including the enthalpy of oxidation and standard enthalpy of formation, by conducting thermalgravimetric analysis – differential scanning calorimetry coupled with mass spectrometry (TGA-DSC-MS), and high temperature transposed temperature drop and oxide melt drop solution calorimetry. The bulk oxidation of UN (containing a small amount of α-UN 1.5+x ) in air was found to follow a step-wise process characterized by consecutive oxidative reactions UN-UN 1.5+x -UO 2 → UO 2 -UO 3 —N k → UO 3 —N k → UO 3 → U 3 O 8 . TGA results support that the UO 2 – U 2 N 3+x passivating layer delays the onset of rapid bulk oxidation of UN in air up to 662 K. Synchrotron X-ray diffraction (XRD) and extended X-ray absorption fine structure (EXAFS) analyses were performed to characterize UN and its final oxidized product. The standard enthalpy of formation (ΔH° f ) of UN was determined to be –144.4 ± 5.9 kJ/mol·atom, in good agreement with previously determined values from Pt encapsulation and bomb calorimetric experiments. Lastly, a negative linear correlation between ΔH° f and the N/U molar ratio was established based on the thermochemical data obtained in this work and previously reported enthalpies of formations of β-UN 1.5-x and α-UN 1.5+x .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Domain Structures of Ni and NiFe (Oxy)Hydroxide Oxygen-Evolution Catalysts from X-ray Pair Distribution Function Analysis

Ni-Fe (oxy)hydroxides, Ni( i-z )Fe z O x H y , are among the fastest-known water oxidation catalysts in alkaline media on a per-cation basis. At current densities relevant for electrolysis (e.g., >0.5 A/cm( -2 )), mass and electron transport through catalyst films with high mass loading are critical and depend substantially on the extended and intermediate catalyst architecture. We use X-ray pair distribution function (PDF) analysis to determine the intermediate nanostructures of electrodeposited Ni( i-z )Fe z O x H y films. We report the effects of electrodeposition technique (pulsed versus continuous), electrochemical cycling, and Fe content on the structure of the catalyst film. The PDF patterns for Ni( i-z )Fe z O x H y films are best simulated by model structures consisting of brucite-like β-Ni(OH) 2 fragments 1 to 3 layers in thickness. Only the oxidation state of the film significantly affects the intralayer scattering behavior (i.e., metal-oxygen bond distance). The interlayer interactions, however, are affected by Fe content and deposition conditions. The domain size of many of the systems are similar, extending to similar to ~5 nm, which are best modeled by sheets containing upward of similar to ~250 metal cations. Smaller domains were found for films deposited through a larger number of electrochemical cathodic current pulses. Films can be cycled between as-deposited, oxidized, and reduced states, with minimal loss of intrasheet coherence, indicating a degree of structural stability. We estimate heterogeneity in the domain structures by modeling the PDF data to linear combinations of oxyhydroxide fragments with different sizes and numbers of layers.

36 MATERIALS SCIENCE↗

Stability analysis and future singularity of the m{sup 2} R □{sup -2} R model of non-local gravity

We analyse the classical stability of the model proposed by Maggiore and Mancarella, where gravity is modified by a term ∼ m{sup 2} R □{sup -2} R to produce the late-time acceleration of the expansion of the universe. Our study takes into account all excitations of the metric that can potentially drive an instability. There are some subtleties in identifying these modes, as a non-local field theory contains dynamical fields which yet do not correspond to degrees of freedom. Since some of them are ghost-like, we clarify the impact of such modes on the stability of the solutions of interest that are the flat space-time and cosmological solutions. We then find that flat space-time is unstable under scalar perturbations, but the instability manifests itself only at cosmological scales, i.e. out of the region of validity of this solution. It is therefore the stability of the FLRW solution which is relevant there, in which case the scalar perturbations are known to be well-behaved by numerical studies. By finding the analytic solution for the late-time behaviour of the scale factor, which leads to a big rip singularity, we argue that the linear perturbations are bounded in the future because of the domination of Hubble friction. In particular, this effect damps the scalar ghost perturbations which were responsible for destabilizing Minkowski space-time. Thus, the model remains phenomenologically viable.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of An Intelligent Flight Propulsion Control System

The initial design and demonstration of an Intelligent Flight Propulsion and Control System (IFPCS) is documented. The design is based on the implementation of a nonlinear adaptive flight control architecture. This initial design of the IFPCS enhances flight safety by using propulsion sources to provide redundancy in flight control. The IFPCS enhances the conventional gain scheduled approach in significant ways: (1) The IFPCS provides a back up flight control system that results in consistent responses over a wide range of unanticipated failures. (2) The IFPCS is applicable to a variety of aircraft models without redesign and,(3) significantly reduces the laborious research and design necessary in a gain scheduled approach. The control augmentation is detailed within an approximate Input-Output Linearization setting. The availability of propulsion only provides two control inputs, symmetric and differential thrust. Earlier Propulsion Control Augmentation (PCA) work performed by NASA provided for a trajectory controller with pilot command input of glidepath and heading. This work is aimed at demonstrating the flexibility of the IFPCS in providing consistency in flying qualities under a variety of failure scenarios. This report documents the initial design phase where propulsion only is used. Results confirm that the engine dynamics and associated hard nonlineaaities result in poor handling qualities at best. However, as demonstrated in simulation, the IFPCS is capable of results similar to the gain scheduled designs of the NASA PCA work. The IFPCS design uses crude estimates of aircraft behaviour. The adaptive control architecture demonstrates robust stability and provides robust performance. In this work, robust stability means that all states, errors, and adaptive parameters remain bounded under a wide class of uncertainties and input and output disturbances. Robust performance is measured in the quality of the tracking. The results demonstrate the flexibility of the IFPCS architecture and the ability to provide robust performance under a broad range of uncertainty. Robust stability is proved using Lyapunov like analysis. Future development of the IFPCS will include integration of conventional control surfaces with the use of propulsion augmentation, and utilization of available lift and drag devices, to demonstrate adaptive control capability under a greater variety of failure scenarios. Further work will specifically address the effects of actuator saturation.

Calise, A. J.↗