Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Instability control”

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 343 records · Page 19

Self-Contained Automated Methodology for Optimal Flow Control

This paper describes a self-contained, automated methodology for active flow control which couples the time-dependent Navier-Stokes system with an adjoint Navier-Stokes system and optimality conditions from which optimal states, i.e., unsteady flow fields and controls (e.g., actuators), may be determined. The problem of boundary layer instability suppression through wave cancellation is used as the initial validation case to test the methodology. Here, the objective of control is to match the stress vector along a portion of the boundary to a given vector; instability suppression is achieved by choosing the given vector to be that of a steady base flow. Control is effected through the injection or suction of fluid through a single orifice on the boundary. The results demonstrate that instability suppression can be achieved without any a priori knowledge of the disturbance, which is significant because other control techniques have required some knowledge of the flow unsteadiness such as frequencies, instability type, etc. The present methodology has been extended to three dimensions and may potentially be applied to separation control, re-laminarization, and turbulence control applications using one to many sensors and actuators.

Joslin, Ronald D.↗

Instabilities caused by floating-point arithmetic quantization.

It is shown that an otherwise stable digital control system can be made unstable by signal quantization when the controller operates on floating-point arithmetic. Sufficient conditions of instability are determined, and an example of loss of stability is treated when only one quantizer is operated.

Phillips, C. L.↗

Transonic adaptive flutter suppression using approximate unsteady time domain aerodynamics

A digital adaptive controller is applied to the active flutter suppression problem of a wing under time varying flight conditions in subsonic and transonic flow. Linear quadratic controller gain at each time step is obtained using an iterative Riccati solver. The digital adaptive optimal controller is robust with respect to the unknown external loads. Flutter and divergence instabilities are simultaneously suppressed using a trailing-edge control surface and displacement sensing. A new transonic unsteady aerodynamic approximation methodology is developed which enables one to carry out the rapid calculation required for transonic aeroservoelastic applications. This approximation is based on a combination of unsteady subsonic aerodynamics combined with a transonic correction procedure. Aeroservoelastic transient time response is obtained using Roger's approximation, state transition matrices and an iterative time marching algorithm. The aeroservoelastic system in the time domain is modeled using a deterministic ARMA model together with a parameter estimator. Transonic flutter boundaries of a wing structure are computed, in the time domain, using an estimated aeroelastic system matrix and are in good agreement with experimental data for the low transonic Mach number range.

Pak, Chan-Gi↗

Adaptive control of large space structures using recursive lattice filters

The use of recursive lattice filters for identification and adaptive control of large space structures is studied. Lattice filters were used to identify the structural dynamics model of the flexible structures. This identification model is then used for adaptive control. Before the identified model and control laws are integrated, the identified model is passed through a series of validation procedures and only when the model passes these validation procedures is control engaged. This type of validation scheme prevents instability when the overall loop is closed. Another important area of research, namely that of robust controller synthesis, was investigated using frequency domain multivariable controller synthesis methods. The method uses the Linear Quadratic Guassian/Loop Transfer Recovery (LQG/LTR) approach to ensure stability against unmodeled higher frequency modes and achieves the desired performance.

Sundararajan, N.↗

Control of Large-Amplitude Combustion Oscillations Using Nanosecond Repetitively Pulsed Plasmas

This paper details the use of nanosecond repetitively pulsed discharges to attenuate combustion instabilities in a 14 kW swirl-stabilized methane/air combustor. The combustor exhibits large-amplitude pressure oscillations ranging from 1 to 4% of the mean pressure during which the flame exhibits bulk motion in each instability cycle, upstream and downstream, as revealed by high-speed chemiluminescence. Control is accomplished with an electrode comprising a pin anode at the centerline of the combustor, allowing a nanosecond spark to be generated in a region spanning close to the flame base, through the shear layers of the swirling flow and ending at the metallic combustor wall. The discharges are generated using 20 kV, 9 kHz pulses; and they correspond to about 120 W of mean power. This results in a suppression of the peak amplitude of the pressure oscillations by a factor of two to four, and 5 dB in the rms value. Using phase-averaged visualizations of the flame with and without plasma, we detail the sequence of flame motion in the course of the instability. With the plasma active, this reveals significant interactions between the flame and the plasma during the suppression. Finally, we present a state-space model of the thermoacoustic system, and we demonstrate open-loop control of the instabilities.

Engineering↗

Generation of large-volume high-pressure plasma by spatiotemporal control of space charge

Any attempt to scale pressure and volume of nonthermal plasma usually leads to instabilities due to the formation of localized space charge. The control of the plasma is limited by the discharge geometry, type of excitation, and gas composition. This article explores the possibility of controlling the space charge in a discharge with a spatially and temporally varying electric field. It is shown that a phase-staggered sinusoidal excitation to a set of conformal azimuthal electrodes in a cylindrical geometry leads to a traveling electric field. Simulations show that in space charge dominated transport, the charged species are dispersed both in the radial and azimuthal directions. This will lead to better control of the space charge and stable discharges near atmospheric pressures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Charge-Based Droop Control Addressing Control Saturation for Low-Inertia Converters

With bulky dc-link capacitors/inductors, traditional voltage-source converters (VSCs) and current-source converters (CSCs) feature a large inertia to facilitate the converter control under large transients. To achieve high power density with reduced cost, low-inertia converters (LICs) featuring significantly reduced dc-link capacitors/inductors have attracted growing attention. However, without bulky dc-link energy buffer, LICs are prone to control saturation under large transients, resulting in undesired oscillation and instability. This issue cannot be managed by traditional proportional-integral based control due to the low inertia. And it will deteriorate when several LICs are connected in series. To address this challenge in LICs, a model-predictive control (MPC) with a computation-inexpensive feed-forward compensation method has been proposed to provide fast dynamic responses. But it would suffer control saturation under larger transients that degrades the control performance. In this article, a charge-based droop control (CDC) is proposed to address this remaining challenge. This paper firstly analyzes the control saturation challenge in LICs by using a tri-port soft-switching solid-state transformer as an example. Next, the operating principle of the proposed CDC are introduced. Two different implementation approaches are discussed in detail. Lastly, the proposed scheme is validated in simulation with a high-fidelity model of hardware prototype. As a result, the proposed CDC eliminated the 2 kHz oscillation and reduced the dc-link ripple and overshoot due to control saturation by 75% and 50%, respectively.

30 DIRECT ENERGY CONVERSION↗

Subspace Iteration Method for Complex Eigenvalue Problems with Nonsymmetric Matrices in Aeroelastic System

Modern airplane design is a multidisciplinary task which combines several disciplines such as structures, aerodynamics, flight controls, and sometimes heat transfer. Historically, analytical and experimental investigations concerning the interaction of the elastic airframe with aerodynamic and in retia loads have been conducted during the design phase to determine the existence of aeroelastic instabilities, so called flutter .With the advent and increased usage of flight control systems, there is also a likelihood of instabilities caused by the interaction of the flight control system and the aeroelastic response of the airplane, known as aeroservoelastic instabilities. An in -house code MPASES (Ref. 1), modified from PASES (Ref. 2), is a general purpose digital computer program for the analysis of the closed-loop stability problem. This program used subroutines given in the International Mathematical and Statistical Library (IMSL) (Ref. 3) to compute all of the real and/or complex conjugate pairs of eigenvalues of the Hessenberg matrix. For high fidelity configuration, these aeroelastic system matrices are large and compute all eigenvalues will be time consuming. A subspace iteration method (Ref. 4) for complex eigenvalues problems with nonsymmetric matrices has been formulated and incorporated into the modified program for aeroservoelastic stability (MPASES code). Subspace iteration method only solve for the lowest p eigenvalues and corresponding eigenvectors for aeroelastic and aeroservoelastic analysis. In general, the selection of p is ranging from 10 for wing flutter analysis to 50 for an entire aircraft flutter analysis. The application of this newly incorporated code is an experiment known as the Aerostructures Test Wing (ATW) which was designed by the National Aeronautic and Space Administration (NASA) Dryden Flight Research Center, Edwards, California to research aeroelastic instabilities. Specifically, this experiment was used to study an instability known as flutter. ATW was a small-scale airplane wing comprised of an airfoil and wing tip boom. This wing was formulated based on a NACA-65A004 airfoil shape with a 3.28 aspect ratio. The wing had a span of 18 inch with root chord length of 13.2 inch and tip chord length of 8.7 inch. The total area of this wing was 197 square inch. The wing tip boom was a 1 inch diameter hollow tube of length 21.5 inch. The total weight of the wing was 2.66 lbs.

Pak, Chan-gi↗

Magnetized Target Fusion Driven by Plasma Liners

Magnetized target fusion (MTF) attempts to combine the favorable attributes of magnetic confinement fusion (MCF) for energy confinement with the attributes of inertial confinement fusion (ICF) for efficient compression heating and wall-free containment of the fusing plasma. It uses a material liner to compress and contain a magnetized plasma. For practical applications, standoff drivers to deliver the imploding momentum flux to the target plasma remotely are required. Spherically converging plasma jets have been proposed as standoff drivers for this purpose. The concept involves the dynamic formation of a spherical plasma liner by the merging of plasma jets, and the use of the liner so formed to compress a spheromak or a field reversed configuration (FRC). For the successful implementation of the scheme, plasma jets of the requisite momentum flux density need to be produced. Their transport over sufficiently large distances (a few meters) needs to be assured. When they collide and merge into a liner, relative differences in velocity, density and temperature of the jets could give rise to instabilities in the development of the liner. Variation in the jet properties must be controlled to ensure that the growth rate of the instabilities are not significant over the time scale of the liner formation before engaging with the target plasma. On impact with the target plasma, some plasma interpenetration might occur between the liner and the target. The operating parameter space needs to be identified to ensure that a reasonably robust and conducting contact surface is formed between the liner and the target. A mismatch in the "impedance" between the liner and the target plasma could give rise to undesirable shock heating of the liner leading to increased entropy (thermal losses) in the liner. Any irregularities in the liner will accentuate the Rayleigh-Taylor instabilities during the compression of the target plasma by the liner.

Thio, Y. C. Francis↗

Laser-plasma instabilities in long scale-length plasmas relevant to shock-ignition

We present 2D particle-in-cell simulations of laser-plasma instabilities (LPI) performed using conditions relevant to experiments on the OMEGA and NIF laser facilities. The laser intensity used is 2×10 15 W cm -2 , which allows comparison with previous experiments and simulations investigating shock ignition, though this is a relatively low intensity for the scheme. We find that the large convective gain of the instabilities leads to dynamics that are largely controlled by pump depletion and to scattering or absorption by LPIs predominantly occurring away from the quarter-critical density. The two-plasmon decay (TPD) instability is dominant for OMEGA-type conditions and remains important at the NIF-scale. In both cases, most absorption from TPD occurs near its Landau cutoff density, and we examine its nonlinear dynamics in this region, focusing on the ion-acoustic waves responsible for saturation. For NIF-type conditions, stimulated Raman scattering (SRS) and stimulated Brillouin scattering divert significant fractions of laser energy, and we find that SRS backscatter can occur with high gain at low densities due to kinetic effects. An in-depth analysis of hot-electron production is performed to identify the main sources and their characteristic electron temperature. For both simulations presented, the overall hot-electron temperature lies between 30 and 35 keV. This is relatively low and will likely increase shock pressure; however, the distributions also contain a significant number of high-energy (E k >100 keV) electrons that would likely cause unacceptable preheat. We, therefore, suggest some strategies that may be utilized to minimize this high-energy component. I. INTROD

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Identifying the microtearing modes in the pedestal of DIII-D H-modes using gyrokinetic simulations

Recent evidence points toward the microtearing mode (MTM) as an important fluctuation in the H-mode pedestal for anomalous electron heat transport. A study of the instabilities in the pedestal region carried out using gyrokinetic simulations to model an ELMy H-mode DIII-D discharge (USN configuration, 1.4 MA plasma current, and 3 MW heating power) is presented. The simulations produce MTMs, identified by predominantly electromagnetic heat flux, small particle flux, and a substantial degree of tearing parity. The magnetic spectrogram from Mirnov coils exhibits three distinct frequency bands—two narrow bands at lower frequency (~35–55 kHz and ~70–105 kHz) and a broader band at higher frequency (~300–500 kHz). Global linear GENE simulations produce MTMs that are centered at the peak of the ω* profile and correspond closely with the bands in the spectrogram. The three distinctive frequency bands can be understood from the basic physical mechanisms underlying the instabilities. For example (i) instability of certain toroidal mode numbers (n) is controlled by the alignment of their rational surfaces with the peak in the ω* profile, and (ii) MTM instabilities in the lower n bands are the conventional collisional slab MTM, whereas the higher n band depends on curvature drive. While many features of the modes can be captured with the local approximation, a global treatment is necessary to quantitatively reproduce the detailed band gaps of the low-n fluctuations. Notably, the transport signatures of the MTM are consistent with careful edge modeling by SOLPS.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical simulations of thermal instabilities in stratified gases. II - Exploration of the parameter space

The temporal evolution of density perturbations in an initially hydrostatic isothermal atmosphere consisting of an optically thin radiating compressible plasma is studied. Numerical techniques are used to describe the nonlinear evolution of the perturbations, and the relative equilibrium between dynamic and thermal instabilities as governed by three independent control parameters are examined, namely, the initial density contrast of the perturbation, the ratio of the local buoyancy oscillation period to the local radiative cooling time, and the ratio of the perturbation radius to the local scaleheight. Four orders of magnitude of initial density contrasts and ratios of buoyancy and cooling times, and one order of magnitude of the bubble dimensions are explored. Well-defined oscillations were found to occur in a limited parameter range, and thermal instability to occur even within secondary condensations deriving from the bubble fragmentation.

Reale, F.↗

Phase development and its role on subharmonic control

The conditions for resonance interaction between two instability waves in an axisymmetric jet were investigated. Considerations of the energy equation of the wave resulting from the interaction indicate that the phase angle between the wave-induced stresses and the wave-induced strains plays a crucial role in the resonance interaction. This fact is demonstrated experimentally by exciting a jet at fundamental and subharmonic frequencies. The phase angle between the waves stresses and strains was varied by varying the initial phase-difference between the two excitation waves. The subharmonic resonance was found to be highly dependent on this angle. Favorable agreement was found between the phase angles predicted by a nonlinear theory and the measured ones. The theory is used to explain the subharmonic's resonance in terms of the phase-angles.

Mankbadi, Reda R.↗

Phase development and its role on subharmonic control

The conditions for resonance interaction between two instability waves in an axisymmetric jet were investigated. Considerations of the energy equation of the wave resulting from the interaction indicate that the phase angle between the wave-induced stresses and the wave-induced strains plays a crucial role in the resonance interaction. This fact is demonstrated experimentally by exciting a jet at fundamental and subharmonic frequencies. The phase angle between the waves' stresses and strains was varied by varying the initial phase-difference between the two excitation waves. The subharmonic resonance was found to be highly dependent on this angle. Favorable agreement was found between the phase angles predicted by a nonlinear theory and the measured ones. The theory is used to explain the subharmonic's resonance in terms of the phase-angles.

Mankbadi, Reda R.↗

Roughness Based Crossflow Transition Control for a Swept Airfoil Design Relevant to Subsonic Transports

A high fidelity transition prediction methodology has been applied to a swept airfoil design at a Mach number of 0.75 and chord Reynolds number of approximately 17 million, with the dual goal of an assessment of the design for the implementation and testing of roughness based crossflow transition control and continued maturation of such methodology in the context of realistic aerodynamic configurations. Roughness based transition control involves controlled seeding of suitable, subdominant crossflow modes in order to weaken the growth of naturally occurring, linearly more unstable instability modes via a nonlinear modification of the mean boundary layer profiles. Therefore, a synthesis of receptivity, linear and nonlinear growth of crossflow disturbances, and high-frequency secondary instabilities becomes desirable to model this form of control. Because experimental data is currently unavailable for passive crossflow transition control for such high Reynolds number configurations, a holistic computational approach is used to assess the feasibility of roughness based control methodology. Potential challenges inherent to this control application as well as associated difficulties in modeling this form of control in a computational setting are highlighted. At high Reynolds numbers, a broad spectrum of stationary crossflow disturbances amplify and, while it may be possible to control a specific target mode using Discrete Roughness Elements (DREs), nonlinear interaction between the control and target modes may yield strong amplification of the difference mode that could have an adverse impact on the transition delay using spanwise periodic roughness elements.

Li, Fei↗

Avoiding fusion plasma tearing instability with deep reinforcement learning

For stable and efficient fusion energy production using a tokamak reactor, it is essential to maintain a high-pressure hydrogenic plasma without plasma disruption. Therefore, it is necessary to actively control the tokamak based on the observed plasma state, to manoeuvre high-pressure plasma while avoiding tearing instability, the leading cause of disruptions. This presents an obstacle-avoidance problem for which artificial intelligence based on reinforcement learning has recently shown remarkable performance. However, the obstacle here, the tearing instability, is difficult to forecast and is highly prone to terminating plasma operations, especially in the ITER baseline scenario. Previously, we developed a multimodal dynamic model that estimates the likelihood of future tearing instability based on signals from multiple diagnostics and actuators. Here we harness this dynamic model as a training environment for reinforcement-learning artificial intelligence, facilitating automated instability prevention. We demonstrate artificial intelligence control to lower the possibility of disruptive tearing instabilities in DIII-D, the largest magnetic fusion facility in the United States. The controller maintained the tearing likelihood under a given threshold, even under relatively unfavourable conditions of low safety factor and low torque. In particular, it allowed the plasma to actively track the stable path within the time-varying operational space while maintaining H-mode performance, which was challenging with traditional preprogrammed control. This controller paves the path to developing stable high-performance operational scenarios for future use in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An analytical and experimental investigation of output feedback vs. linear quadratic regulator

This paper presents analytical and experimental comparison of three control laws for a laboratory structure designed to simulate large space structures. The first control law is the standard time invariant linear quadratic regulator with state estimation, which requires model reduction for practical implementation. This model reduction is associated with the so-called spillover instability. Two new simple direct output feedback control laws guaranteeing stability are proposed. One minimizes the maximum control force, and the other minimizes the same quadratic performance index as the linear quadratic regulator. The three control laws are found to give comparable performance for the modes retained in the reduced model. However, the standard linear quadratic regulator with state estimation provides almost no margin of stability for the unmodeled modes, while the simpler direct feedback laws provide significant stability margins to all modes. The analytical results were verified experimentally using a digital implementation of the control laws. Good agreement between the analytical predictions and experimental measurements was observed.

Martinovic, Zoran N.↗

Aeroservoelastic Stability Analysis of the X-43A Stack

The first air launch attempt of an X-43A stack, consisting of the booster, adapter and Hyper-X research vehicle, ended in failure shortly after the successful drop from the National Aeronautics and Space Administration Dryden Flight Research Center (Edwards, California) B-52B airplane and ignition of the booster. The stack was observed to begin rolling and yawing violently upon reaching transonic speeds, and the grossly oscillating fins of the booster separated shortly thereafter. The flight then had to be terminated with the stack out of control. Very careful linear flutter and aeroservoelastic analyses were subsequently performed as reported herein to numerically duplicate the observed instability. These analyses properly identified the instability mechanism and demonstrated the importance of including the flight control laws, rigid-body modes, structural flexible modes and control surface flexible modes. In spite of these efforts, however, the predicted instability speed remained more than 25 percent higher than that observed in flight. It is concluded that transonic shock phenomena, which linear analyses cannot take into account, are also important for accurate prediction of this mishap instability.

Pak, Chan-gi↗