Search NASA⌕ Search

SEARCH · Search NASA

Results for “STABLE OSCILLATION”

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 469 records · Page 26

Optimal Control Modification for Time-Scale Separated Systems

Recently a new optimal control modification has been introduced that can achieve robust adaptation with a large adaptive gain without incurring high-frequency oscillations as with the standard model-reference adaptive control. This modification is based on an optimal control formulation to minimize the L2 norm of the tracking error. The optimal control modification adaptive law results in a stable adaptation in the presence of a large adaptive gain. This study examines the optimal control modification adaptive law in the context of a system with a time scale separation resulting from a fast plant with a slow actuator. A singular perturbation analysis is performed to derive a modification to the adaptive law by transforming the original system into a reduced-order system in slow time. A model matching conditions in the transformed time coordinate results in an increase in the actuator command that effectively compensate for the slow actuator dynamics. Simulations demonstrate effectiveness of the method.

Adaptive Control↗

Radiofrequency synthesiser with an intrinsic instability of 5 × 10{sup −15} at the averaging time of 1 s based on a femtosecond erbium-doped fibre laser

A radio-optical synthesiser intended for operation in a radio-frequency master oscillator with an optical He – Ne/CH{sub 4} frequency standard (λ = 3.39 μm) is developed on the basis of a femtosecond erbium-doped fibre laser. The synthesiser generates equidistant harmonics in the frequency range of 1 – 10 GHz with stability determined by the optical frequency standard. A stable supercontinuum spectrum is formed in the range around 1.06 μm, which provides stable 24-hour operation of the synthesiser and is important for practical applications in off-laboratory conditions. A direct comparison of the output frequencies of two synthesisers shows that the up-grade of the fibre laser and detection system of femtosecond pulses results in the synthesiser intrinsic instability of 5 × 10{sup −15} at the averaging time of 1 s. Such a value is by an order of magnitude less than that obtained in our earlier works. (laser frequency standards)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Generating Orbital Angular Momentum Beams in an FEL Oscillator

With the storage ring FEL at Duke University, we have established fundamental harmonic FEL lasing in a variety of coherently mixed orbital angular momentum (OAM) modes, with photons in the |l> and |-l> superposition states (l = 1, 2, 3) and with each eigenstate carrying l_h OAM. Spatially organized in Laguerre Gaussian modes, stable FEL lasing of these OAM beams has been achieved with reasonable power. Optical techniques have been developed to characterize these beams. The operation of such an OAM FEL paves the way for the generation of OAM gamma-ray beams via Compton scattering.

Wu, Y.K.↗

Galileo Probe Doppler Residuals as the Wave-Dynamical Signature of Weakly Stable, Downward-Increasing Stratification in Jupiter's Deep Wind Layer

Doppler radio tracking of the Galileo probe-to-orbiter relay, previously analyzed for its in situ measure of Jupiter's zonal wind at the equatorial entry site, also shows a record of significant residual fluctuations apparently indicative of varying vertical motions. Regular oscillations over pressure depth in the residual Doppler measurements of roughly 1-8 Hz (increasing upward), as filtered over a 134 sec window, are most plausibly interpreted as gravity waves, and imply a weak, but downward increasing static stability within the 5 - 20 bar region of Jupiter's atmosphere. A matched extension to deeper levels of an independent inertial stability constraint from the measured vertical wind shear at 1 - 4 bars is roughly consistent with a static stability of approximately 0.5 K/km near the 20 bar level, as independently detected by the probe Atmospheric Structure Instrument.

Allison, Michael↗

ELM Suppression by Boron Powder Injection and Comparison with Lithium Powder Injection on EAST

Type I edge-localized modes (ELMs) in the Experimental Advanced Superconducting Tokamak (EAST) were completely suppressed via boron powder injection into the X-point region of an upper-single null configuration over a wide range of operating conditions (2.8 < P aux < 7.5 MW, 3.8 × 10 19 < n e < 6 × 10 19 m -3 , RF-only and RF + NBI heating scenarios, both grad-B drift directions, and even He ion majority plasmas) (Sun et al. in Nucl. Fusion, 2020). A window of edge B concentration for stable long pulse operation was identified: too low and ELMs return, too high and the discharge suffers radiative collapse. The injection of boron powder above the minimum for ELM suppression coincided with the occurrence of an edge harmonic oscillation detected in magnetics (both on the high-field side and low-field side), in AXUV diodes near the upper X-point, divertor Dα emission, and in a range of other diagnostics (Diallo et al. in: Proceedings of 2020 IAEA fusion energy conference, 2021). Furthermore, no harmonic oscillation was observed when ELMs were present, and stored energy was slightly increased at constant density during ELM suppression. Core tungsten emission during ELM suppression either increased or decreased relative to ELMy H-mode, but the W emission was maintained at acceptable levels. The threshold B injection rate was measured for several conditions, and found to increase with heating power. Li powder injection into comparable discharges also resulted in a short phase of ELM suppression, but density and stored energy both decreased due to the strong pumping effect of lithium; no edge harmonic oscillation was observed with Li injection, indicating that the ELM suppression mechanisms differ. The new set of B-seeded, ELM-suppressed discharges exhibited certain characteristics of quiescent H-mode (Burrell et al. in Phys Plasmas 8:2153, 2001), but did not require high shear, counter beams, etc. The wide operating window and compatibility with RF-only discharges paves the way for future experiments targeting long pulse H-mode discharges with complete ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Non-isobaric Thermal Instability

Multiphase media have very complex structure and evolution. Accurate numerical simulations are necessary to make advances in our understanding of this rich physics. Because simulations can capture both the linear and nonlinear evolution of perturbations with a relatively wide range of sizes, it is important to thoroughly understand the stability of condensation and acoustic modes between the two extreme wavelength limits of isobaric and isochoric instability as identified by Field. Partially motivated by a recent suggestion that large non-isobaric clouds can "shatter" into tiny cloudlets, we revisit the linear theory to survey all possible regimes of thermal instability. We uncover seven regimes in total, one of which allows three unstable condensation modes. Using the code Athena++, we determine the numerical requirements to properly evolve small amplitude perturbations of the entropy mode into the nonlinear regime. Here, our 1D numerical simulations demonstrate that for a typical AGN cooling function, the nonlinear evolution of a single eigenmode in an isobarically unstable plasma involves increasingly larger amplitude oscillations in cloud size, temperature, and density as the wavelength increases. Such oscillations are the hallmark behavior of non-isobaric multiphase gas dynamics and may be observable as correlations between changes in brightness and the associated periodic redshifts and blueshifts in systems that can be spatially resolved. Intriguingly, we discuss regimes and derive characteristic cloud sizes for which the saturation process giving rise to these oscillations can be so energetic that the cloud may indeed break apart. However, we dub this process "splattering" instead of "shattering," as it is a different fragmentation mechanism that is triggered when the cloud suddenly "lands" on the stable cold branch of the equilibrium curve.

79 ASTRONOMY AND ASTROPHYSICS↗

Trajectory Design for the Transiting Exoplanet Survey Satellite

The Transiting Exoplanet Survey Satellite (TESS) is a National Aeronautics and Space Administration (NASA) mission, scheduled to be launched in 2017. TESS will travel in a highly eccentric orbit around Earth, with initial perigee radius near 17 Earth radii (Re) and apogee radius near 59 Re. The orbit period is near 2:1 resonance with the Moon, with apogee nearly 90 degrees out-of-phase with the Moon, in a configuration that has been shown to be operationally stable. TESS will execute phasing loops followed by a lunar flyby, with a final maneuver to achieve 2:1 resonance with the Moon. The goals of a resonant orbit with long-term stability, short eclipses and limited oscillations of perigee present significant challenges to the trajectory design. To rapidly assess launch opportunities, we adapted the Schematics Window Methodology (SWM76) launch window analysis tool to assess the TESS mission constraints. To understand the long-term dynamics of such a resonant orbit in the Earth-Moon system we employed Dynamical Systems Theory in the Circular Restricted 3-Body Problem (CR3BP). For precise trajectory analysis we use a high-fidelity model and multiple shooting in the General Mission Analysis Tool (GMAT) to optimize the maneuver delta-V and meet mission constraints. Finally we describe how the techniques we have developed can be applied to missions with similar requirements. Keywords: resonant orbit, stability, lunar flyby, phasing loops, trajectory optimization

Stationkeepting↗

Trajectory Design for the Transiting Exoplanet Survey Satellite (TESS)

The Transiting Exoplanet Survey Satellite (TESS) is a National Aeronautics and Space Administration (NASA) mission launching in 2017. TESS will travel in a highly eccentric orbit around Earth, with initial perigee radius near 17 Earth radii (Re) and apogee radius near 59 Re. The orbit period is near 2:1 resonance with the Moon, with apogee nearly 90 degrees out-of-phase with the Moon, in a configuration that has been shown to be operationally stable. TESS will execute phasing loops followed by a lunar flyby, with a final maneuver to achieve 2:1 resonance with the Moon. The goals of a resonant orbit with long-term stability, short eclipses and limited oscillations of perigee present significant challenges to the trajectory design. To rapidly assess launch opportunities, we adapted the SWM76 launch window tool to assess the TESS mission constraints. To understand the long-term dynamics of such a resonant orbit in the Earth-Moon system we employed Dynamical Systems Theory in the Circular Restricted 3-Body Problem (CR3BP). For precise trajectory analysis we use a high-fidelity model and multiple shooting in the General Mission Analysis Tool (GMAT) to optimize the maneuver delta-V and meet mission constraints. Finally we describe how the techniques we have developed can be applied to missions with similar requirements.

Monte Carlo↗

Effects of Extreme Obliquity Variations on the Habitability of Exoplanets

We explore the impact of obliquity variations on planetary habitability in hypothetical systems with high mutual inclination. We show that large-amplitude, high-frequency obliquity oscillations on Earth-like exoplanets can suppress the ice-albedo feedback, increasing the outer edge of the habitable zone. We restricted our exploration to hypothetical systems consisting of a solar-mass star, an Earth-mass planet at 1 AU, and 1 or 2 larger planets. We verified that these systems are stable for 108 years with N-body simulations and calculated the obliquity variations induced by the orbital evolution of the Earth-mass planet and a torque from the host star. We ran a simplified energy balance model on the terrestrial planet to assess surface temperature and ice coverage on the planet's surface, and we calculated differences in the outer edge of the habitable zone for planets with rapid obliquity variations. For each hypothetical system, we calculated the outer edge of habitability for two conditions: (1) the full evolution of the planetary spin and orbit and (2) the eccentricity and obliquity fixed at their average values. We recovered previous results that higher values of fixed obliquity and eccentricity expand the habitable zone, but we also found that obliquity oscillations further expand habitable orbits in all cases. Terrestrial planets near the outer edge of the habitable zone may be more likely to support life in systems that induce rapid obliquity oscillations as opposed to fixed-spin planets. Such planets may be the easiest to directly characterize with space-borne telescopes.

habitability of exoplanets↗

Energy-Screened Many-Body Expansion for Protein–Ligand Interactions: Examining Convergence for Metalloenzymes Through Seven–Body Interactions

Fragment-based quantum chemistry is a powerful strategy for calculating protein−ligand interaction energies using quantum chemistry methods. Rigorous convergence often requires hundreds of atoms in the protein binding-site model, especially if that model is constructed using distance-based criteria to select amino acid residues, while three- and four-body calculations exhibit instability related to combinatorial proliferation in the number of subsystem calculations. Here, we report an energy-based screening protocol for the many-body expansion applied to protein−ligand interactions, implemented in the open-source FRAGME∩T code. Using a combination of aggressive screening based on semiempirical quantum chemistry, with an improved graph-theoretical algorithm to eliminate unimportant subsystems, we are able to perform n-body calculations up to n = 7 using density functional theory in triple-ζ basis sets. Distance cutoffs further reduce the cost without compromising accuracy. Rapid and stable convergence of the many-body expansion is obtained by n = 4, for a pair of metalloenzymes in which a divalent ion coordinates directly to the ligand. As compared to previous results that relied solely on distance cutoffs, oscillations in the n-body corrections are reduced or eliminated, although residual errors remain in one case. This work demonstrates that benchmark-quality protein−ligand interaction energies can be systematically converged using a method with excellent parallel efficiency and scalability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of turbulence on a kinetic auroral arc model

A plasma kinetic model of an inverted-V auroral arc structure which includes the effects of electrostatic turbulence is proposed. In the absence of turbulence, a parallel potential drop is supported by magnetic mirror forces and charge quasi neutrality, with energetic auroral ions penetrating to low altitudes; relative to the electrons, the ions' pitch angle distribution is skewed toward smaller pitch angles. The electrons energized by the potential drop form a current which excites electrostatic turbulence. In equilibrium the plasma is marginally stable. The conventional anomalous resistivity contribution to the potential drop is very small. Anomalous resistivity processes are far too dissipative to be powered by auroral particles. It is concluded that under certain circumstances equilibrium may be impossible and relaxation oscillations set in.

Cornwall, J. M.↗

An upwind differencing scheme for the time-accurate incompressible Navier-Stokes equations

The two-dimensional incompressible Navier-Stokes equations are solved in a time-accurate manner in using the method of pseudocompressibility. Using this method, subiterations in pseudotime are required to satisfy the continuity equation at each time step. An upwind differencing scheme based on flux-difference splitting is used to compute the convective terms. The upwind differencing is biased based on the sign of the local eigenvalue of the Jacobian matrix. Third-order or fifth-order spatial accuracy is maintained throughout the interior grid points. The equations are solved using an implicit line-relaxation scheme. This solution scheme is stable and is capable of running at large time steps in pseudotime, leading to fast convergence for each physical time step. A variety of computed results are presented to validate the present scheme. Results for the flow over an oscillating plate are compared with the exact analytic solution, good agreement is seen. Excellent comparison is obtained between the computed solution and the analytical results for inviscid channel flow with an oscillating back pressure. Flow solutions over a circular cylinder with vortex shedding are also presented. Finally, the flow past an airfoil at -90 deg angle-of-attack is also computed.

Rogers, Stuart E.↗

Theory of Dust Voids in Plasmas

Dusty plasmas in a gas discharge often feature a stable void, i.e., a dust-free region inside the dust cloud. This occurs under conditions relevant to both plasma processing discharges and plasma crystal experiments. The void results from a balance of the electrostatic and ion drag forces on a dust particle. The ion drag force is driven by a flow of ions outward from an ionization source and toward the surrounding dust cloud, which has a negative space charge. In equilibrium the force balance for dust particles requires that the boundary with the dust cloud be sharp, provided that the particles are cold and monodispersive. Numerical solutions of the one-dimensional nonlinear fluid equations are carried out including dust charging and dust-neutral collisions, but not ion-neutral collisions. The regions of parameter space that allow stable void equilibria are identified. There is a minimum ionization rate that can sustain a void. Spatial profiles of plasma parameters in the void are reported. In the absence of ion-neutral collisions, the ion flow enters the dust cloud's edge at Mach number M = 1. Phase diagrams for expanding or contracting voids reveal a stationary point corresponding to a single stable equilibrium void size, provided the ionization rate is constant. Large voids contract and small voids expand until they attain this stationary void size. On the other hand, if the ionization rate is not constant, the void size can oscillate. Results are compared to recent laboratory and microgravity experiments.

Goree, J.↗

Analysis and Optimization of the Characteristics of Superpower Virtode Generation

We report on the results of the 3D numerical simulation of perspective design of an oscillator based on an external-feedback virtual cathode (virtode), as well as optimization of characteristics of its generation. We have analyzed the effect of beam parameters on the virtode output characteristics. In particular, we have obtained the dependence of the efficiency on the feeding voltage and on the current being injected. We have studied the effect of beam velocity premodulation and electron velocity spread on virtode generation characteristics. It has been established that the introduction of premodulation sharply increases the generation efficiency up to 13% for optimal parameters. The possibility of stable generation in the virtode at relatively large spread in electron velocities has been demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analysis of the Spectral Stability of the Generalized Runge–Kutta Methods Applied to Initial-Boundary-Value Problems for Equations of the Parabolic Type. I. Explicit Methods

We develop a general algorithm for the analysis the spectral stability of generalized multistage Runge–Kutta methods of various orders of accuracy as applied to the numerical integration with respect to time of an initial-boundary-value problem for the second-order parabolic equation. The expression for the function of spectral stability is obtained in two alternative forms: on the basis of matrix relations and in the determinant form. We study specific realizations of various generalized explicit Runge–Kutta methods and their spectral stability. It is shown that all explicit generalized Runge–Kutta methods possess the property of conditional spectral stability and the property of conditional monotonicity of the numerical solution in time whose violation leads to the appearance of false oscillations of the approximate solution. The stability function for these methods is polynomial. It is shown that the application of two-stage generalized explicit Runge–Kutta methods leads to the appearance of predictor-corrector-type schemes. In the case of nonstationary one-dimensional heat-conduction problem, on the basis of a one-stage generalized Runge–Kutta method, we obtain the classical two-layer conditionally stable explicit finite-difference scheme on a four-point template. It is demonstrated that the five-stage generalized Runge–Kutta–Merson method is characterized by the weakest condition of spectral stability as compared with all investigated explicit generalized Runge–Kutta methods.

42 ENGINEERING↗

Investigation of dynamic noise affecting geodynamics information in a tethered subsatellite

The effects of a tethered satellite system's internal dynamics on the subsatellite were calculated including both overall motions (libration and attitude oscillations) and internal tether oscillations. The SKYHOOK tether simulation program was modified to operate with atmospheric density variations and to output quantities of interest. Techniques and software for analyzing the results were developed including noise spectral analysis. A program was begun for computing a stable configuration of a tether system subject to air drag. These configurations will be of use as initial conditions for SKYHOOK and, through linearized analysis, directly for stability and dynamical studies. A case study in which the subsatellite traverses an atmospheric density enhancement confirmed some theoretical calculations, and pointed out some aspects of the interaction with the tether system dynamics.

Gullahorn, G. E.↗

Accretion by magnetic neutron stars

The capture of material from the wind or from the atmosphere or envelope of a binary companion star is described and the resulting types of accretion flows are discussed. An explanation is given for the formation of a magnetosphere around the neutron star. Conditions for stable flow and for angular and linear momentum conservation are discussed within the framework of neutron star accretion and rough estimates of the scale of the magnetosphere are obtained. The beat-frequency magnetospheric model of quasi-periodic intensity oscillations is presented.

Lamb, F. K.↗

Direct numerical simulation of instabilities in parallel flow with spherical roughness elements

Results from a direct numerical simulation of laminar flow over a flat surface with spherical roughness elements using a spectral-element method are given. The numerical simulation approximates roughness as a cellular pattern of identical spheres protruding from a smooth wall. Periodic boundary conditions on the domain's horizontal faces simulate an infinite array of roughness elements extending in the streamwise and spanwise directions, which implies the parallel-flow assumption, and results in a closed domain. A body force, designed to yield the horizontal Blasius velocity in the absence of roughness, sustains the flow. Instabilities above a critical Reynolds number reveal negligible oscillations in the recirculation regions behind each sphere and in the free stream, high-amplitude oscillations in the layer directly above the spheres, and a mean profile with an inflection point near the sphere's crest. The inflection point yields an unstable layer above the roughness (where U''(y) is less than 0) and a stable region within the roughness (where U''(y) is greater than 0). Evidently, the instability begins when the low-momentum or wake region behind an element, being the region most affected by disturbances (purely numerical in this case), goes unstable and moves. In compressible flow with periodic boundaries, this motion sends disturbances to all regions of the domain. In the unstable layer just above the inflection point, the disturbances grow while being carried downstream with a propagation speed equal to the local mean velocity; they do not grow amid the low energy region near the roughness patch. The most amplified disturbance eventually arrives at the next roughness element downstream, perturbing its wake and inducing a global response at a frequency governed by the streamwise spacing between spheres and the mean velocity of the most amplified layer.

Deanna, R. G.↗