Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical stability”

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 811 records · Page 45

Design of feedforward controllers for multivariable plants

Simple methods for the design of feedforward controllers to achieve steady-state disturbance rejection and command tracking in stable multivariable plants are developed in this paper. The controllers are represented by simple and low-order transfer functions and are not based on reconstruction of the states of the commands and disturbances. For unstable plants, it is shown that the present method can be applied directly when an additional feedback controller is employed to stabilize the plant. The feedback and feedforward controllers do not affect each other and can be designed independently based on the open-loop plant to achieve stability, disturbance rejection and command tracking, respectivley. Numerical examples are given for illustration.

Seraji, H.↗

Progress on a Taylor weak statement finite element algorithm for high-speed aerodynamic flows

A new finite element numerical Computational Fluid Dynamics (CFD) algorithm has matured to the point of efficiently solving two-dimensional high speed real-gas compressible flow problems in generalized coordinates on modern vector computer systems. The algorithm employs a Taylor Weak Statement classical Galerkin formulation, a variably implicit Newton iteration, and a tensor matrix product factorization of the linear algebra Jacobian under a generalized coordinate transformation. Allowing for a general two-dimensional conservation law system, the algorithm has been exercised on the Euler and laminar forms of the Navier-Stokes equations. Real-gas fluid properties are admitted, and numerical results verify solution accuracy, efficiency, and stability over a range of test problem parameters.

Baker, A. J.↗

Nonlinear interactions in mixing layers and compressible heated round jets

The nonlinear interactions between a fundamental instability mode and both its harmonics and the changing mean flow are studied using the weakly nonlinear stability theory of Stuart and Watson, and numerical solutions of coupled nonlinear partial differential equations. The first part focuses on incompressible cold (or isothermal; constant temperature throughout) mixing layers, and for these, the first and second Landau constants are calculated as functions of wavenumber and Reynolds number. It is found that the dominant contribution to the Landau constants arises from the mean flow changes and not from the higher harmonics. In order to establish the range of validity of the weakly nonlinear theory, the weakly nonlinear and numerical solutions are compared and the limitation of each is discussed. At small amplitudes and at low-to-moderate Reynolds numbers, the two results compare well in describing the saturation of the fundamental, the distortion of the mean flow, and the initial stages of vorticity roll-up. At larger amplitudes, the interaction between the fundamental, second harmonic, and the mean flow is strongly nonlinear and the numerical solution predicts flow oscillations, whereas the weakly nonlinear theory yields saturation. In the second part, the weakly nonlinear theory is extended to heated (or nonisothermal; mean temperature distribution) subsonic round jets where quadratic and cubic nonlinear interactions are present, and the Landau constants also depend on jet temperature ratio, Mach number and azimuthal mode number. Under exponential growth and nonlinear saturation, it is found that heating and compressibility suppress the growth of instability waves, that the first azimuthal mode is the dominant instability mode, and that the weakly nonlinear solution describes the early stages of the roll-up of an axisymmetric shear layer. The receptivity of a typical jet flow to pulse type input disturbance is also studied by solving the initial value problem and then examining the behavior of the long-time solution.

Jarrah, Yousef Mohd↗

Multivariable flight control synthesis and literal robustness analysis for an aeroelastic vehicle

An integrated flight/aeroelastic control law is developed analytically for a hypothetical large supersonic transport aircraft in which the first aeroelastic mode frequency of the fuselage (6 rad/sec) is near the short-period mode (2 rad/sec). The approach employed is based on a linear-quadratic-regulator (LQR) formulation (yielding model-following state-feedback gains), followed by asymptotic loop-transfer recovery of LQR robustness (to produce an output-feedback control law). The derivation is outlined, and numerical results comparing the performance and multivariate stability robustness of the present controller with those of a classical controller are presented in graphs. The two controllers are shown to have similar characteristics, even with respect to the sources of limitations on robustness.

Schmidt, David K.↗

On the design of decoupling controllers for advanced rotorcraft in the hover case

A methodology for design of helicopter control systems is proposed that can account for various types of concurrent specifications: stability, decoupling between longitudinal and lateral motions, handling qualities, and physical limitations of the swashplate motions. This is achieved by synergistic use of analytical techniques (Q-parameterization of all stabilizing controllers, transfer function interpolation) and advanced numerical optimization techniques. The methodology is used to design a controller for the UH-60 helicopter in hover. Good results are achieved for decoupling and handling quality specifications.

Fan, M. K. H.↗

Analytical multiple scattering correction to the Mie theory: Application to the analysis of the lidar signal

The contribution of the multiple scattering to the lidar signal is dependent on the optical depth tau. Therefore, the radar analysis, based on the assumption that the multiple scattering can be neglected is limited to cases characterized by low values of the optical depth (tau less than or equal to 0.1) and hence it exclude scattering from most clouds. Moreover, all inversion methods relating lidar signal to number densities and particle size must be modified since the multiple scattering affects the direct analysis. The essential requests of a realistic model for lidar measurements which include the multiple scattering and which can be applied to practical situations follow. (1) Requested are not only a correction term or a rough approximation describing results of a certain experiment, but a general theory of multiple scattering tying together the relevant physical parameter we seek to measure. (2) An analytical generalization of the lidar equation which can be applied in the case of a realistic aerosol is requested. A pure analytical formulation is important in order to avoid the convergency and stability problems which, in the case of numerical approach, are due to the large number of events that have to be taken into account in the presence of large depth and/or a strong experimental noise.

Flesia, C.↗

Interface behavior of a multi-layer fluid configuration subject to acceleration in a microgravity environment, supplement 1

With the increasing opportunities for research in a microgravity environment, there arises a need for understanding fluid mechanics under such conditions. In particular, a number of material processing configurations involve fluid-fluid interfaces which may experience instabilities in the presence of external forcing. In a microgravity environment, these accelerations may be periodic or impulse-type in nature. This research investigates the behavior of a multi-layer idealized fluid configuration which is infinite in extent. The analysis is linear, and each fluid region is considered inviscid, incompressible, and immiscible. An initial parametric study of confiquration stability in the presence of a constant acceleration field is performed. The zero mean gravity limit case serves as the base state for the subsequent time-dependent forcing cases. A stability analysis of the multi-layer fluid system in the presence of periodic forcing is investigated. Floquet theory is utilized. A parameter study is performed, and regions of stability are identified. For the impulse-type forcing case, asymptotic stability is established for the configuration. Using numerical integration, the time response of the interfaces is determined.

Lyell, M. J.↗

Studies of Premixed Laminar and Turbulent Flames at Microgravity

Several topics related to the experimental and theoretical study of the combustion limits in premixed flames at reduced gravity have been investigated. The topics studied include: (1) numerical simulation of flame ball structure and stability; (2) interferometry for micro-g combustion experiments; (3) simulation of buoyancy effects in premixed flames using aqueous autocatalytic reactions; (4) premixed flame propagation in Hele-Shaw cells; and (5) experimental characterization of flame balls.

Ronney, Paul D.↗

Annual Research Briefs - 2000: Center for Turbulence Research

This report contains the 2000 annual progress reports of the postdoctoral Fellows and visiting scholars of the Center for Turbulence Research (CTR). It summarizes the research efforts undertaken under the core CTR program. Last year, CTR sponsored sixteen resident Postdoctoral Fellows, nine Research Associates, and two Senior Research Fellows, hosted seven short term visitors, and supported four doctoral students. The Research Associates are supported by the Departments of Defense and Energy. The reports in this volume are divided into five groups. The first group largely consists of the new areas of interest at CTR. It includes efficient algorithms for molecular dynamics, stability in protoplanetary disks, and experimental and numerical applications of evolutionary optimization algorithms for jet flow control. The next group of reports is in experimental, theoretical, and numerical modeling efforts in turbulent combustion. As more challenging computations are attempted, the need for additional theoretical and experimental studies in combustion has emerged. A pacing item for computation of nonpremixed combustion is the prediction of extinction and re-ignition phenomena, which is currently being addressed at CTR. The third group of reports is in the development of accurate and efficient numerical methods, which has always been an important part of CTR's work. This is the tool development part of the program which supports our high fidelity numerical simulations in such areas as turbulence in complex geometries, hypersonics, and acoustics. The final two groups of reports are concerned with LES and RANS prediction methods. There has been significant progress in wall modeling for LES of high Reynolds number turbulence and in validation of the v(exp 2) - f model for industrial applications.

Source record↗

Ultra-Light Precision Membrane Optics

SRS Technologies and NASA Marshall Space Flight Center have conducted a research effort to explore the possibility of developing ultra-lightweight membrane optics for future imaging applications. High precision optical flats and spherical mirrors were produced under this research effort. The thin film mirrors were manufactured using surface replication casting of CPI(Trademark), a polyimide material developed specifically for UV hardness and thermal stability. In the course of this program, numerous polyimide films were cast with surface finishes better than 1.5 nanometers rms and thickness variation of less than 63 nanometers. Precision membrane optical flats were manufactured demonstrating better than 1/13 wave figure error when measured at 633 nanometers. The aerial density of these films is 0.037 kilograms per square meter. Several 0.5-meter spherical mirrors were also manufactured. These mirrors had excellent surface finish (1.5 nanometers rms) and figure error on the order of tens of microns. This places their figure error within the demonstrated correctability of advanced wavefront correction technologies such as real time holography.

Moore, Jim↗

Robust Control Design for Systems With Probabilistic Uncertainty

This paper presents a reliability- and robustness-based formulation for robust control synthesis for systems with probabilistic uncertainty. In a reliability-based formulation, the probability of violating design requirements prescribed by inequality constraints is minimized. In a robustness-based formulation, a metric which measures the tendency of a random variable/process to cluster close to a target scalar/function is minimized. A multi-objective optimization procedure, which combines stability and performance requirements in time and frequency domains, is used to search for robustly optimal compensators. Some of the fundamental differences between the proposed strategy and conventional robust control methods are: (i) unnecessary conservatism is eliminated since there is not need for convex supports, (ii) the most likely plants are favored during synthesis allowing for probabilistic robust optimality, (iii) the tradeoff between robust stability and robust performance can be explored numerically, (iv) the uncertainty set is closely related to parameters with clear physical meaning, and (v) compensators with improved robust characteristics for a given control structure can be synthesized.

Crespo, Luis G.↗

Enabling the Next Generation of Spaceborne Quadrupole Mass Spectrometers

The quadrupole mass spectrometer (QMS) has over 30 years of spaceflight heritage in making important neutral gas and low energy ion observations. Given their geometrical constraints, these instruments are currently operated at the extreme limit of their capabilities. However, a technique called higher order auxiliary excitation provides a set of novel, robust, electronics-based solutions for improving the performance of these sensors. By driving the quadrupole rods with an additional frequency nearly twice that of the normal RF operating frequency, substantially increased abundance sensitivity, maximum attainable mass resolution, and peak stability can be achieved through operation of voltage scan lines through the center of formed upper stability islands. Such improvements are modeled using numerical simulations of ion trajectories in a quadrupole field with and without applied higher order auxiliary excitation. When compared to a traditional QMS with a mass range up to 500Da, sensors can be designed with the same precision electronics to have expected mass ranges beyond 1500Da with a power increase of less than twice that of its heritage implementations.

Gershman, D. J.↗

Time Domain Approach for Acoustic Liner Analysis

- Procedures for implementing the Ingard-Myers impedance condition in time domain - Instability issue of the Ingard-Myers Condition - Proposed stabilization of the Ingard-Myers condition - Numerical examples

time domain impedance boundary condition↗

Learning Local Volt/VAR Controllers Toward Efficient Network Operation with Stability Guarantees: Preprint

This paper considers the problem of voltage regulation in distribution network. The primary motivation is to keep voltages within pre-assigned operating limits by commanding the reactive power output of distributed energy resources (DERs) deployed in the grid. We develop a framework for developing local Volt/Var control that comprises of two main steps. In the first, exploiting historical data and for each DER, we learn a function representing desirable equilibrium points for the power network. These points approximate solutions of an Optimal Power Flow problem. In the second, we propose a control scheme for steering the network towards these favorable configurations. Theoretical conditions are derived to formally guarantee the stability of the developed control scheme and numerical simulations illustrate the effectiveness of the proposed approach.

data-driven control↗

Theoretical analyses of baroclinic flows

Completed and ongoing research activities are discussed briefly, including a three-dimensional, linear stability analysis of the baroclinic Hadley cell and a numerical model of the baroclinic flow between two rotating concentric spheres. This model simulates axisymmetric flow in the Atmospheric General Circulation Experiment configuration. A computer code designed to solve the strongly nonlinear stability problem for the Eady basic state is mentioned.

Antar, B. N.↗

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.↗

The structure and stability of coronal magnetic fields

The high degree of symmetry often assumed in studies of the structure and stability of coronal magnetic field configurations is restrictive and can yield misleading results. Fully three-dimensional numerical methods for constructing force-free equilibria and for examining their stability properties are developed which make no assumptions about symmetry. A test of the stability analysis has been performed by applying it to the Gold-Hoyle twisted flux tube, which is known to be kink-unstable if the helical field makes more than about one turn between the line-tying end-plates. The preliminary result is that the critical number of turns is about 1.1, in good agreement with the previous best estimate. However, it is found that the growth rate, which has not been discussed previously, is orders of magnitude smaller than expected, even when the flux tube is twisted far beyond the stability limit.

Mcclymont, A. N.↗

Numerical Analysis of Flow Evolution in a Helium Jet Injected into Ambient Air

A computational model to study the stability characteristics of an evolving buoyant helium gas jet in ambient air environment is presented. Numerical formulation incorporates a segregated approach to solve for the transport equations of helium mass fraction coupled with the conservation equations of mixture mass and momentum using a staggered grid method. The operating parameters correspond to the Reynolds number varying from 30 to 300 to demarcate the flow dynamics in oscillating and non-oscillating regimes. Computed velocity and concentration fields were used to analyze the flow structure in the evolving jet. For Re=300 case, results showed that an instability mode that sets in during the evolution process in Earth gravity is absent in zero gravity, signifying the importance of buoyancy. Though buoyancy initiates the instability, below a certain jet exit velocity, diffusion dominates the entrainment process to make the jet non-oscillatory as observed for the Re=30 case. Initiation of the instability was found to be dependent on the interaction of buoyancy and momentum forces along the jet shear layer.

Satti, Rajani P.↗