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 523 records · Page 29

Hydrodynamic Stability of Liquid-Propellant Combustion: Landau's Problem Revisited

Hydrodynamic, or Landau, instability in combustion is typically associated with the onset of wrinkling of a flame surface, corresponding to the formation of steady cellular structures as the stability threshold is crossed. As its name suggests, it stems from hydrodynamic effects connected with thermal expansion across the reaction region. In the context of liquid-propellant combustion, the classical models that originally predicted this phenomenon have been extended to include the important effects that arise from a dynamic dependence of the burning rate on the local pressure and temperature fields. Thus, the onset of Landau instability has now been shown to occur for sufficiently small negative values of the pressure sensitivity of the burning rate, significantly generalizing previous classical results for this problem that assumed a constant normal burning rate. It has also been shown that the onset of instability occurs for decreasing values of the disturbance wave number as the gravitational-acceleration parameter decreases. Consequently, in an appropriate weak-gravity limit, Landau instability becomes a long-wave phenomena associated with the formation of large cells on the liquid-propellant surface. Additionally, a pulsating form of hydrodynamic instability has been shown to occur as well, corresponding to the onset of temporal oscillations in the location of the liquid/gas interface. This instability occurs for sufficiently large negative values of the pressure sensitivity, and is enhanced by increasing values of the burning-rate temperature sensitivity. It is further shown that for sufficiently small values of this parameter, there exists a stable range of pressure sensitivities for steady, planar burning such that the classical cellular form of hydrodynamic instability and the more recent pulsating form of hydrodynamic instability can each occur as the corresponding stability threshold is crossed. For larger thermal sensitivities, however, the pulsating stability boundary evolves into a C-shaped curve in the (disturbance-wave number, pressure-sensitivity) plane, indicating loss of stability to pulsating perturbations for all sufficiently large disturbance wavelengths. It is thus concluded, based on characteristic parameter values that an equally likely form of hydrodynamic instability in liquid-propellant combustion is of a nonsteady, long-wave nature, distinct from the steady, cellular form originally predicted by Landau.

Margolis, S. B.↗

Quantification of Variation in Combustion Instability Amplitude in a Multi-Nozzle Can Combustor

In this work, we quantify the level of variation in unstable combustion oscillation amplitudes and identify the source of this variation in a multi-nozzle can combustor. At conditions where pollutant emissions are reduced, lean-premixed combustors can undergo thermoacoustic instability when pressure and heat release rate oscillations couple. A commonly used method for suppressing instability is fuel staging, which is a method where fuel is unevenly distributed between nozzles in a multi-nozzle combustor. Our work follows others who have characterized the effect of fuel staging on combustion instability and the mechanisms by which it works during both steady-state and transient operation. One of the outcomes from our previous work was that certain instability operating points display a high level of pressure oscillation amplitude variation. Instead of oscillating at a constant limit-cycle amplitude, the pressure oscillation amplitude varies significantly in time. In this work, we use the concept of permutation entropy to quantify the level of variation in the pressure fluctuation amplitude. We correlate the level of variation with a number of state variables over a range of operating conditions, 291 test cases in all. These state variables include mixture equivalence ratio; transient timescale, amplitude, and direction; hardware temperatures; gas temperatures; and thermoacoustic damping and growth rates. Significant pressure oscillation amplitude variation occurs when the thermoacoustic damping rate is low. The damping and growth rates can be low for a number of reasons, but they are highly correlated with the metal temperature of the centerbody, where the flames are anchored; lower temperatures result in lower damping rates during stable operation and lower growth rates during unstable operation. These results show the importance of the thermal boundary condition on the time-dependent behavior of the thermoacoustic instability.

Westfall, Seth↗

Evolution of the Heliotail Lobes over a Solar Cycle as Measured by IBEX

Energetic neutral atom (ENA) measurements by IBEX reveal that the heliotail comprises an energy-dependent multilobe structure. We examine the heliotail evolution over 11 yr of IBEX observations covering a full solar cycle (SC). We find the following: (1) The heliotail structure persists over the entire SC, comprising three ENA-enhanced and two ENA-suppressed lobes. (2) Lobe sizes and locations are generally stable but exhibit variations in ENA fluxes driven by the SC. (3) Lobe centers follow a cyclic behavior over multiple SC phases, indicating direct signatures of slow and fast solar wind (SW) interactions in the inner heliosheath (IHS). (4) The tilted plane passing through the port–starboard lobes' centers oscillates in latitude but maintains its tilt from the ecliptic plane, likely a consequence of the interstellar magnetic field draping around the heliosphere. (5) The transition of the central heliotail from a single lobe at ~1.1 keV to two lobes above ~2 keV is SC-dependent and directly reflects the IHS plasma properties, i.e., when ENA fluxes from fast SW from the polar coronal holes change over time. (6) The central lobe exhibits a substructure that is enhanced and offset from the downwind direction, possibly indicating an asymmetric ENA emission or an asymmetry in the parent plasma distribution. These results reveal the general stability of the heliotail structure over time and distinct variations in individual lobes' properties in relation to the SC phases. Furthermore, results show the effects of multiple SC phases in the tail, reflecting different ENA travel times and source histories.

79 ASTRONOMY AND ASTROPHYSICS↗

Pulsating Hydrodynamic Instability and Thermal Coupling in an Extended Landau/Levich Model of Liquid-Propellant Combustion: Inviscid Analysis - 1

Hydrodynamic (Landau) instability in combustion is typically associated with the onset of wrinkling of a flame surface, corresponding to the formation of steady cellular structures as the stability threshold is crossed. In the context of liquid-propellant combustion, such instability has recently been shown to occur for critical values of the pressure sensitivity of the burning rate and the disturbance wavenumber, significantly generalizing previous classical results for this problem that assumed a constant normal burning rate. Additionally, however, a pulsating form of hydrodynamic instability has been shown to occur as well, corresponding to the onset of temporal oscillations in the location of the liquid/gas interface. In the present work, we consider the realistic influence of a non-zero temperature sensitivity in the local burning rate on both types of stability thresholds. It is found that for sufficiently small values of this parameter, there exists a stable range of pressure sensitivities for steady, planar burning such that the classical cellular form of hydrodynamic instability and the more recent pulsating form of hydrodynamic instability can each occur as the corresponding stability threshold is crossed. For larger thermal sensitivities, however, the pulsating stability boundary evolves into a C-shaped curve in the (disturbance-wavenumber, pressure-sensitivity) plane, indicating loss of stability to pulsating perturbations for all sufficiently large disturbance wavelengths. It is thus concluded, based on characteristic parameter values, that an equally likely form of hydrodynamic instability in liquid-propellant combustion is of a non-steady, long-wave nature, distinct from the steady, cellular form originally predicted by Landau.

Margolis, Stephen B.↗

Sonoluminescence in Space: The Critical Role of Buoyancy in Stability and Emission Mechanisms

Sonoluminescence is the term used to describe the emission of light from a violently collapsing bubble. Sonoluminescence ("light from sound") is the result of extremely nonlinear pulsations of gas/vapor bubbles in liquids when subject to sufficiently high amplitude acoustic pressures. In a single collapse, a bubble's volume can be compressed more than a thousand-fold in the span of less than a microsecond. Even the simplest consideration of the thermodynamics yields pressures on the order of 10,000 ATM, and temperatures of at least 10,000K. On the face of things, it is not surprising that light should be emitted from such an extreme process. Since 1990 (the year that Gaitan discovered light from a single bubble) there has been a tremendous amount of experimental and theoretical research in stable, single-bubble sonoluminescence (SBSL), yet there remain at least four unexplained phenomena associated with SBSL in 1g: the light emission mechanism itself, the existence of anisotropies in the emitted light, the disappearance of the bubble at some critical acoustic pressure, and the appearance of quasiperiodic and chaotic oscillations in the flash timing. Gravity, in the context of the buoyant force, is implicated in all four of these. We are developing KC-135 experiments probing the effect of gravity on single bubble sonoluminescence. By determining the stability boundaries experimentally in microgravity, and measuring not only light emission but mechanical bubble response, we will be able to directly test the predictions of existing theories.

Thomas, Charles R.↗

The prediction of nonlinear longitudinal combustion instability in liquid propellant rockets

An analytical technique was developed to solve nonlinear longitudinal combustion instability problems. The analysis yields the transient and limit cycle behavior of unstable motors and the threshold amplitude required to trigger a linearly stable motor into unstable operation. The limit cycle waveforms were found to exhibit shock wave characteristics for most unstable engine operating conditions. A method of correlating the analytical solutions with experimental data was developed. Calculated results indicate that a second-order solution adequately describes the behavior of combustion instability oscillations over a broad range of engine operating conditions, but that higher order effects must be accounted for in order to investigate engine triggering.

Lores, M. E.↗

Implicit Total Variation Diminishing (TVD) schemes for steady-state calculations

The application of a new implicit unconditionally stable high resolution total variation diminishing (TVD) scheme to steady state calculations. It is a member of a one parameter family of explicit and implicit second order accurate schemes developed by Harten for the computation of weak solutions of hyperbolic conservation laws. This scheme is guaranteed not to generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments show that this scheme not only has a rapid convergence rate, but also generates a highly resolved approximation to the steady state solution. A detailed implementation of the implicit scheme for the one and two dimensional compressible inviscid equations of gas dynamics is presented. Some numerical computations of one and two dimensional fluid flows containing shocks demonstrate the efficiency and accuracy of this new scheme.

Yee, H. C.↗

Implicit total variation diminishing (TVD) schemes for steady-state calculations

The application of a new implicit unconditionally stable high resolution total variation diminishing (TVD) scheme to steady state calculations. It is a member of a one parameter family of explicit and implicit second order accurate schemes developed by Harten for the computation of weak solutions of hyperbolic conservation laws. This scheme is guaranteed not to generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments show that this scheme not only has a rapid convergence rate, but also generates a highly resolved approximation to the steady state solution. A detailed implementation of the implicit scheme for the one and two dimensional compressible inviscid equations of gas dynamics is presented. Some numerical computations of one and two dimensional fluid flows containing shocks demonstrate the efficiency and accuracy of this new scheme. Previously announced in STAR as N83-23085

Yee, H. C.↗

Three-dimensional instability of rotating flows with oscillating axial strain

The equations of motion for perturbed uniformly rotating flows with uniform axial-time periodic strain, are derived from the Navier-Stokes equations in the low Mach number limit. The perturbation equations admit exponentially growing three-dimensional solutions for which the amplification factors per period are computed for a range of compression and swirl ratios. It is found that for a given compression ratio, the flow is stable for swirl ratios, but at high swirl ratios the flow is unstable with the amplification factor dependent on wave angle but independent of wavelength. For an unstable swirl ratio, higher compression ratios yield larger amplification factors.

Mansour, Nagi N.↗

The effect of artificial bulk viscosity in simulations of forced compressible turbulence

The use of an artificial bulk viscosity for shock stabilization is a common approach employed in turbulence simulations with high-order numerics. The effect of the artificial bulk viscosity is analyzed in the context of large eddy simulations by using as a test case simulations of linearly-forced compressible homogeneous turbulence (Petersen and Livescu, 2010 [12]). This case is unique in that it allows for the specification of a priori target values for total dissipation and ratio of solenoidal to dilatational dissipation. A comparison between these target values and the true predicted levels of dissipation is thus used to investigate the performance of the artificial bulk viscosity. Results show that the artificial bulk viscosity is effective at achieving stable solutions, but also leads to large values of artificial dissipation that outweigh the physical dissipation caused by fluid viscosity. An alternate approach, which employs the artificial thermal conductivity only, shows that the dissipation of dilatational modes is entirely due to the fluid viscosity. However, this method leads to unwanted Gibbs oscillations around the shocklets. The use of shock sensors that further localize the artificial bulk viscosity did not reduce the amount of artificial dissipation introduced by the artificial bulk viscosity. Finally, an improved forcing function that explicitly accounts for the role of the artificial bulk viscosity in the budget of turbulent kinetic energy was explored.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Pygmy resonances: what’s in a name?

The centroid, width and percentage of energy weighted sum rule of dipole resonances can be strongly affected by dynamical fluctuations and static deformations of the nuclear surface, deformations and fluctuations which, in turn, depend on pairing, and thus on Cooper pairs. Because of angular momentum conservation, such insight is restricted, to lowest order, to fluctuations/deformations of quadrupole and monopole type. The latter being closely connected with the neutron (excess) skin and thus with soft dipole modes. From the values ( N - Z)/ A ≈ 0.18, 0.21, and 0.45 for the nuclei 122 Sn, 208 Pb, and 11 Li, it is expected that the latter system, which is weakly bound by pairing effects (spatially extended single Cooper pair and odd proton acting as spectator), constitutes an attractive laboratory to study the properties of soft E 1-modes and thus of isospin nuclear deformation. From the calculation of the full dipole response function in QRPA, discretizing the continuum in a spherical box of radius of 40 fm, one finds a GDR with centroid E x ≈ 24 MeV, width Γ ≈11 MeV and carrying 90% of the EWSR, and a low-lying collective resonance characterized by E X = 0.75 MeV, Γ = 0.5 MeV and 6.2% EWSR. The wave function of the latter resonance is built out of about fifteen components (both protons and neutrons), typical of a collective mode. The transition densities indicate this soft E 1 - mode to be generated by surface density oscillation of the neutron (halo) skin (Δ r np ≈ 1.71 fm) relative to an approximately isospin-saturated core. Through a detailed study of the full dipole response of 11 Li we will draw a comparison between the soft E 1 - mode of this halo nucleus and the PDR of heavy stable nuclei, pointing to the physical similarities and also to the basic differences.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analytical and experimental studies of impinging liquid jets

Impinging injectors are a common type of injector used in liquid propellant rocket engines and are typically used in engines where both propellants are injected as a liquid, e.g., engines using LOX/hydrocarbon and storable propellant combinations. The present research program is focused on providing the requisite fundamental understanding associated with impinging jet injectors for the development of an advanced a priori combustion stability design analysis capability. To date, a systematic study of the atomization characteristics of impinging liquid jets under cold-flow conditions have been completed. Effects of orifice diameter, impingement angle, pre-impingement length, orifice length-to-diameter ratio, fabrication procedure, jet flow condition and jet velocity under steady and oscillating, and atmospheric- and high-pressure environments have been investigated. Results of these experimental studies have been compared to current models of sheet breakup and drop formation. In addition, the research findings have been scrutinized to provide a fundamental explanation for a proven empirical correlation used in the design of stable impinging injector-based rocket engines.

Ryan, H. M.↗

Analysis of High Order Difference Methods for Multiscale Complex Compressible Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes with incremental studies was initiated. Here we further refine the analysis on, and improve the understanding of the adaptive numerical dissipation control strategy. Basically, the development of these schemes focuses on high order nondissipative schemes and takes advantage of the progress that has been made for the last 30 years in numerical methods for conservation laws, such as techniques for imposing boundary conditions, techniques for stability at shock waves, and techniques for stable and accurate long-time integration. We concentrate on high order centered spatial discretizations and a fourth-order Runge-Kutta temporal discretizations as the base scheme. Near the bound-aries, the base scheme has stable boundary difference operators. To further enhance stability, the split form of the inviscid flux derivatives is frequently used for smooth flow problems. To enhance nonlinear stability, linear high order numerical dissipations are employed away from discontinuities, and nonlinear filters are employed after each time step in order to suppress spurious oscillations near discontinuities to minimize the smearing of turbulent fluctuations. Although these schemes are built from many components, each of which is well-known, it is not entirely obvious how the different components be best connected. For example, the nonlinear filter could instead have been built into the spatial discretization, so that it would have been activated at each stage in the Runge-Kutta time stepping. We could think of a mechanism that activates the split form of the equations only at some parts of the domain. Another issue is how to define good sensors for determining in which parts of the computational domain a certain feature should be filtered by the appropriate numerical dissipation. For the present study we employ a wavelet technique introduced in as sensors. Here, the method is briefly described with selected numerical experiments.

Sjoegreen, Bjoern↗

Emergence of 4 H $J^π$ = 1 - resonance in contact theories

We obtain the s- and p-wave low-energy scattering parameters for n 3 H elastic scattering and the position of the 4 H J π = 1 - resonance using the pionless effective field theory at leading order. Results are extracted with three numerical techniques: confining the system in a harmonic oscillator trap, solving the Faddeev-Yakubovsky equations in configuration space, and using an effective two-body cluster approach. The renormalization of the theory for the relevant amplitudes is assessed in a cutoff-regulator range between 1 fm -1 and 10 fm -1 . Most remarkably, we find a cutoff-stable/RG-invariant resonance in the 4 H J π = 1 - system. This p-wave resonance is a universal consequence of a shallow two-body state and the introduction of a three-body s-wave scale set by the triton binding energy. The stabilization of a resonant state in a few-fermion system through pure contact interactions has a significant consequence for the powercounting of the pionless theory. Specifically, it suggests the appearance of similar resonant states also in larger nuclei, like 16-oxygen, in which the theory’s leading order does not predict stable states. Those resonances would provide a starting state to be moved to the correct physical position by the perturbative insertion of sub-leading orders, possibly resolving the discrepancy between data and contact EFT.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of Transformer Magnetizing Reactance Using Describing Function in Direct Power Control of Back-to-Back Modular Multilevel Converter with Advanced Grid Support

This paper presents a Port Controlled Hamiltonian based direct power control architecture for a back-to-back modular multilevel converter system connecting two ac sources at different frequencies. The system features advanced grid support functionalities based on IEEE 1547-2018, implemented on the inverter side. The rectifier side controller ensures reference following for active and/or reactive powers and maintains the commanded dc bus voltage. The inverter side controller ensures power command following for active and reactive powers.The proposed control architecture is designed to suppress second harmonic oscillations in powers during unbalanced grid voltage sags by dynamically adjusting the currents on each ac side. This also ensures effective elimination of any second harmonic oscillations in the equivalent dc bus voltage. Validation is performed on an OPAL-RT real-time platform with case studies on unbalanced and balanced sags, demonstrating the controller's effectiveness during real time implementation. A reduced-scale laboratory prototype further verifies these case studies, with experimental results for balanced sags due to grid simulator limitations. The results confirm the robustness and efficiency of the proposed control strategy in ensuring stable and reliable operation under various grid conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Dynamic analysis of an open-cycle gas-core nuclear rocket

Reactivity and flow disturbances were used to investigate the transient response of a conceptual open cycle gas core nuclear rocket engine. The disturbances were made with the system initially operating at its steady state design point. Results of the study show that the feedbacks associated with the propellant density and propellant temperature have a dominant effect on the response of the system. Furthermore, there appears to be a rather limited range of values of these propellant feedback coefficients for which the gas core nuclear rocket has a stable response. The system was rather insensitive to a fuel flow rate disturbance, whereas a similar disturbance in the propellant flow rate caused large changes in reactor power. For a similar disturbance in the propellant flow rate caused large changes in reactor power. For most reactivity and flow rate disturbances, the response showed oscillations of various intensity.

Turney, G. E.↗

Nonlinear calculation of a radial mode of oscillation of a hot white dwarf

A linearly unstable first harmonic radial mode of oscillation (period 0.86 s) of a hot white dwarf with a hydrogen envelope is analyzed in the nonlinear regime by application of the Stellingwerf technique to a nonlinear hydrocode. The mass and equilibrium luminosity of the hot white dwarf are 1.3 solar masses, and 1.6 x 10 to the 37th erg/s, respectively. A stable nonlinear mode is found to exist, in which appreciable shock damping occurs at and near the surface - the shock damping being the dominant stabilizing effect. This enhances the possibility that such an object can produce a pulsating X-ray source.

Degregoria, A. J.↗

Robust Optimal Adaptive Control Method with Large Adaptive Gain

In the presence of large uncertainties, a control system needs to be able to adapt rapidly to regain performance. Fast adaptation is referred to the implementation of adaptive control with a large adaptive gain to reduce the tracking error rapidly. However, a large adaptive gain can lead to high-frequency oscillations which can adversely affect robustness of an adaptive control law. A new adaptive control modification is presented that can achieve robust adaptation with a large adaptive gain without incurring high-frequency oscillations as with the standard model-reference adaptive control. The modification is based on the minimization of the Y2 norm of the tracking error, which is formulated as an optimal control problem. The optimality condition is used to derive the modification using the gradient method. The optimal control modification results in a stable adaptation and allows a large adaptive gain to be used for better tracking while providing sufficient stability robustness. Simulations were conducted for a damaged generic transport aircraft with both standard adaptive control and the adaptive optimal control modification technique. The results demonstrate the effectiveness of the proposed modification in tracking a reference model while maintaining a sufficient time delay margin.

Nguyen, Nhan T.↗