Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FREE 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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Unsteady Extinction of Opposed Jet Ethylene/Methane HIFiRE Surrogate Fuel Mixtures vs Air

A unique idealized study of the subject fuel vs. air systems was conducted using an Oscillatory-input Opposed Jet Burner (OOJB) system and a newly refined analysis. Extensive dynamic-extinction measurements were obtained on unanchored (free-floating) laminar Counter Flow Diffusion Flames (CFDFs) at 1-atm, stabilized by steady input velocities (e.g., U(sub air)) and perturbed by superimposed in-phase sinusoidal velocity inputs at fuel and air nozzle exits. Ethylene (C2H4) and methane (CH4), and intermediate 64/36 and 15/85 molar percent mixtures were studied. The latter gaseous surrogates were chosen earlier to mimic ignition and respective steady Flame Strengths (FS = U(sub air)) of vaporized and cracked, and un-cracked, JP-7 "like" kerosene for a Hypersonic International Flight Research Experimentation (HIFiRE) scramjet. For steady idealized flameholding, the 100% C2H4 flame is respectively approx. 1.3 and approx.2.7 times stronger than a 64/36 mix and CH4; but is still 12.0 times weaker than a 100% H2-air flame. Limited Hot-Wire (HW) measurements of velocity oscillations at convergent-nozzle exits, and more extensive Probe Microphone (PM) measurements of acoustic pressures, were used to normalize Dynamic FSs, which decayed linearly with pk/pk U(sub air) (velocity magnitude, HW), and also pk/pk P (pressure magnitude, PM). Thus Dynamic Flame Weakening (DFW) is defined as % decrease in FS per Pascal of pk/pk P oscillation, namely, DFW = -100 d(U(sub air)/U(sub air),0Hz)/d(pkpk P). Key findings are: (1) Ethylene flames are uniquely strong and resilient to extinction by oscillating inflows below 150 Hz; (2) Methane flames are uniquely weak; (3) Ethylene / methane surrogate flames are disproportionately strong with respect to ethylene content; and (4) Flame weakening is consistent with limited published results on forced unsteady CFDFs. Thus from 0 to approx. 10 Hz and slightly higher, lagging diffusive responses of key species led to progressive phase lags (relative to inputs) in the oscillating flames, and caused maximum weakening. At 20 to 150 Hz, diffusion-rate-limited effects diminished, causing flames to "regain strengnth," and eventually become completely insensitive beyond 300 Hz. Detailed mechanistic understanding is needed. Overall, ethylene flames are remarkably resilient to dynamic extinction by oscillating inflows. They are the strongest, with the notable exception of H2. For HIFiRE tests, the 64%/36% surrogate disproportionally retains the high dynamic FS of ethylene, so the potential for loss of scramjet flameholding (flameout) due to low frequency oscillations is significantly mitigated.

Vaden, Sarah N.↗
NASA NTRS1985

Force-free equilibria of magnetized jets

Force-free equilibrium configurations of magnetic-pressure-dominated magnetized supersonic jets confined by slowly varying external pressure are investigated analytically. For the case where internal dissipation mechanisms are active, the lowest-energy field configuration is found to be the superposition of an axisymmetric mode and a helical mode with a wavelength equal to 5 times the jet radius, and the pressure below which the nonaxisymmetric mode becomes energetically favorable is given as 2700 times the product of the 4th power of the magnetic helicity per unit length and the -6th power of the magnetic flux. A model of the total and polarized emission of such a configuration is developed and applied to the extended well-collimated astronomically resolved jet NGC 6251. The model is shown to reproduce significant features such as transverse oscillations of the ridge line, width oscillations and emission knots, the projected magnetic-field configuration, oscillations of the degree of polarization, and the distribution of the Faraday rotation measure.

Koenigl, A.↗
NASA NTRS1961

Stability and Control Characteristics of a Small-Scale Model of an Aerial Vehicle Supported by Two Ducted Fans

An investigation has been made to determine the stability and control characteristics in hovering and in forward flight of a free-flight model representing a type of vertical-take-off-and-landing aircraft which utilizes two fixed ducted fans as its sole source of lift and propulsion. The model, having fans 28 inches in diameter, was considered to be approximately one-third the size of a full-scale aircraft. Control moments for most of the hovering tests and all the forward-flight tests were provided by remotely controlled compressed-air jets at the sides and ends of the model. For one brief phase of the hovering investigation a system of vanes in the duct slipstreams was substituted for the jets as a means of roll control. During the forward-flight tests, the model was flown with both the tandem and side-by-side duct arrangements. In hovering the model exhibited strongly divergent oscillations about the pitch and roll axes. The pitching oscillation of the tandem configuration was of a fairly long period and was not particularly difficult to control; the rolling oscillation, however, was of a relatively short period and was extremely difficult to control. Both oscillations could be completely eliminated by the addition of a sufficient amount of artificial damping. The control moments produced by the vane-type roll control system were weak and were accompanied by a side force of appreciable magnitude and undesirable direction. In forward flight the model required an undesirably large nose-down tilt angle for equilibrium at any appreciable speed. A vane was placed transversely in the slipstream of the forward duct of the tandem configuration in an attempt to reduce this tilt angle. The vane was effective in reducing the tilt angle but apparently caused an increase in the power requirements and in the angle-of-attack instability. Without the vane, a forward speed of 30 knots (full scale) required a nose-down tilt angle of about 300. A powerful pitch control moment was required not only to maintain the trim attitude but also to 2 overcome the effects of instability with angle of attack. Less pitch control moment was required for the tandem configuration than for the side-by-side configuration at any given forward speed. The instability in roll increased with forward speed. No forward speeds in excess of about 20 knots (full scale) were achieved until the artificial damping in roll and the yaw control moment were increased appreciably above values which had proved satisfactory for hovering flight.

Parlett, Lysle P.↗
NASA NTRS2011

Excitation of a Parallel Plate Waveguide by an Array of Rectangular Waveguides

This work addresses the problem of excitation of a parallel plate waveguide by an array of rectangular waveguides that arises in applications such as the continuous transverse stub (CTS) antenna and dual-polarized parabolic cylindrical reflector antennas excited by a scanning line source. In order to design the junction region between the parallel plate waveguide and the linear array of rectangular waveguides, waveguide sizes have to be chosen so that the input match is adequate for the range of scan angles for both polarizations. Electromagnetic wave scattered by the junction of a parallel plate waveguide by an array of rectangular waveguides is analyzed by formulating coupled integral equations for the aperture electric field at the junction. The integral equations are solved by the method of moments. In order to make the computational process efficient and accurate, the method of weighted averaging was used to evaluate rapidly oscillating integrals encountered in the moment matrix. In addition, the real axis spectral integral is evaluated in a deformed contour for speed and accuracy. The MoM results for a large finite array have been validated by comparing its reflection coefficients with corresponding results for an infinite array generated by the commercial finite element code, HFSS. Once the aperture electric field is determined by MoM, the input reflection coefficients at each waveguide port, and coupling for each polarization over the range of useful scan angles, are easily obtained. Results for the input impedance and coupling characteristics for both the vertical and horizontal polarizations are presented over a range of scan angles. It is shown that the scan range is limited to about 35 for both polarizations and therefore the optimum waveguide is a square of size equal to about 0.62 free space wavelength.

Rengarajan, Sembiam↗
NASA NTRS2008

Supersonic Disk Gap Band Parachute Performance in the Wake of a Viking-Type Aeroshell from Mach 2 to 2.5

Supersonic wind tunnel testing of 0.813 m diameter Disk-Gap-Band parachutes is being conducted in the NASA Glenn Research Center (GRC) 10' x 10' wind-tunnel. The tests are conducted in support of the Mars Science Laboratory Parachute Decelerator System development and qualification. Four percent of full-scale parachutes were constructed similarly to the flight-article in material and construction techniques. The parachutes are attached to a 4% scale MSL entry-vehicle to simulate the free-flight configuration. The parachutes are tested from Mach 2 to 2.5 over a Reynolds number (Re) range of 1 to 3 x 10(exp 6), representative of the MSL deployment envelope. Constrained and unconstrained test configurations are investigated to quantify the effects of parachute trim, suspension line interaction, and alignment with the capsule wake. The parachute is constrained horizontally through the vent region, to measure canopy breathing and wake interaction for fixed trim angles of 0 and 10 degrees from the velocity vector. In the unconstrained configuration the parachute is permitted to trim and cone, similar to the free-flight varying its alignment relative to the entry-vehicle wake. Test diagnostics were chosen to quantify parachute performance and to provide insight into the flow field structure. An in-line load cell provided measurement of unsteady and mean drag as a function of Mach and Re. High-speed shadowgraph video of the upstream parachute flow field was used to capture bow-shock motion and stand of distance. Particle image velocimetry of the upstream parachute flow field provides spatially and temporally resolved measurement velocity and turbulent statistics. Multiple high speed video views of targets placed in the interior of the canopy enable photo-grammetric measurement of the fabric motion in time and space from reflective. High speed video is also used to document the supersonic inflation and measure trim angle, projected area, and frequency of area oscillations.

EDL↗
NASA NTRS2011

Intraseasonal Variations in Tropical Energy Balance: Relevance to Climate Sensitivity?

Intraseasonal variability of deep convection represents a fundamental mode of organization for tropical convection. While most studies of intraseasonal oscillations (ISOs) have focused on the spatial propagation and dynamics of convectively coupled circulations, here we examine the projection of ISOs on the tropically-averaged heat and moisture budget. One unresolved question concerns the degree to which observable variations in the "fast" processes (e.g. convection, radiative / turbulent fluxes) can inform our understanding of feedback mechanisms operable in the context of climate change. Our analysis use daily data from satellite observations, the Modern Era analysis for Research and Applications (MERRA), and other model integrations to address these questions: (i) How are tropospheric temperature variations related to that tropical deep convection and the associated ice cloud fractional amount (ICF), ice water path (IWP), and properties of warmer liquid clouds? (ii) What role does moisture transport play vis-a-vis ocean latent heat flux in enabling the evolution of deep convection to mediate PBL - free atmospheric temperature equilibration? (iii) What affect do convectively generated upper-tropospheric clouds have on the TOA radiation budget? Our methodology is similar to that of Spencer et al., (2007 GRL ) whereby a composite time series of various quantities over 60+ ISO events is built using tropical mean tropospheric temperature signal as a reference to which the variables are related at various lag times (from -30 to +30 days). The area of interest encompasses the global oceans between 20oN/S. The increase of convective precipitation cannot be sustained by evaporation within the domain, implying strong moisture transports into the tropical ocean area. The decrease in net TOA radiation that develops after the peak in deep convective rainfall, is part of the response that constitutes a "discharge" / "recharge" mechanism that facilitates tropical heat balance maintenance on these time scales. However, water vapor and hydrologic scaling relationships for this mode of variability cast doubt on the utility of ISO variations as proxies for climate sensitivity response to external radiatively forced (e.g. greenhouse gas-induced) climate change.

Robertson, Franklin R.↗
NASA NTRS2020

Detection and Mitigation of Transient Instabilities in Multi-Agent Systems and Swarms

We first introduce the novel concept of transient instabilities in multi-agent systems and swarms, i.e., a small disturbance leads to increasing-amplitude oscillations throughout the swarm, which results in a large number of inter-agent collisions. This instability is dominant in the transient phase of the system and it does not appear in the steady-state behavior of the system, as each agent uses Lyapunov-stable feedback control laws. We present a rigorous definition of transient instability in swarms, and discuss its key properties. We also present a sufficient condition to check if a swarm will be transient stable. We study the behaviour of different control laws under this condition. We also show how transient instability phenomenon could impact the dynamics of deployable structures. Next, we present a novel control architecture that augments the baseline formation maintenance controller to mitigate transient instabilities. At its heart, the proposed architecture consists of a projection operator based estimator disguised as a reference model that generates collision-free trajectories for the agents to follow. We present numerical simulation results to demonstrate the effectiveness of our proposed approaches.

Quadrelli, Marco B.↗
NASA NTRS2004

Transitional Bubble in Periodic Flow Phase Shift

One particular characteristic observed in unsteady shear layers is the phase shift relative to the main flow. In attached boundary layers this will have an effect both on the instantaneous skin friction and heat transfer. In separation bubbles the contribution to the drag is dominated by the pressure distribution. However, the most significant effect appears to be the phase shift on the transition process. Unsteady transition behaviour may determine the bursting of the bubble resulting in an un-recoverable full separation. An early analysis of the phase shift was performed by Stokes for the incompressible boundary layer of an oscillating wall and an oscillating main flow. An amplitude overshoot within the shear layer as well as a phase shift were observed that can be attributed to the relatively slow diffusion of viscous stresses compared to the fast change of pressure. Experiments in a low speed facility with the boundary layer of a flat plate were evaluated in respect to phase shift. A pressure distribution similar to that on the suction surface of a turbomachinery aerofoil was superimposed generating a typical transitional separation bubble. A periodically unsteady main flow in the suction type wind tunnel was introduced via a rotating flap downstream of the test section. The experiments covered a range of the three similarity parameters of momentum-loss-thickness Reynolds-number of 92 to 226 and Strouhal-number (reduced frequency) of 0.0001 to 0.0004 at the separation point, and an amplitude range up to 19 %. The free stream turbulence level was less than 1% .Upstream of the separation point the phase shift in the laminar boundary layer does not appear to be affected significantly bay either of the three parameters. The trend perpendicular to the wall is similar to the Stokes analysis. The problem scales well with the wave velocity introduced by Stokes, however, the lag of the main flow near the wall is less than indicated analytically. The separation point comes closest to the Stokes analysis but the phase is still 20 degrees lower at the wall.

Talan, M.↗
NASA NTRS1998

Inner Core Rotation from Geomagnetic Westward Drift and a Stationary Spherical Vortex in Earth's Core

The idea that geomagnetic westward drift indicates convective leveling of the planetary momentum gradient within Earth's core is pursued in search of a differentially rotating mean state, upon which various oscillations and secular effects might be superimposed. The desired state conforms to roughly spherical boundary conditions, minimizes dissipative interference with convective cooling in the bulk of the core, yet may aid core cooling by depositing heat in the uppermost core and lower mantle. The variational calculus of stationary dissipation applied to a spherical vortex within the core yields an interesting differential rotation profile, akin to spherical Couette flow bounded by thin Hartmann layers. Four boundary conditions are required. To concentrate shear induced dissipation near the core-mantle boundary, these are taken to be: (i) no-slip at the core-mantle interface; (ii) geomagnetically estimated bulk westward flow at the base of the core-mantle boundary layer; (iii) no-slip at the inner-outer core interface; and, to describe magnetic locking of the inner core to the deep outer core; (iv) hydrodynamically stress-free at the inner-outer core boundary. By boldly assuming the axial core angular momentum anomaly to be zero, the super-rotation of the inner core relative to the mantle is calculated to be at most 1.5 deg./yr.

Voorhies, Coerte V.↗
NASA NTRS1999

Inner Core Rotation from Geomagnetic Westward Drift and a Stationary Spherical Vortex in Earth's Core

The idea that geomagnetic westward drift indicates convective leveling of the planetary momentum gradient within Earth's core is pursued in search of a differentially rotating mean state, upon which various oscillations and secular effects might be superimposed. The desired state conforms to roughly spherical boundary conditions, minimizes dissipative interference with convective cooling in the bulk of the core, yet may aide core cooling by depositing heat in the uppermost core and lower mantle. The variational calculus of stationary dissipation applied to a spherical vortex within the core yields an interesting differential rotation profile akin to spherical Couette flow bounded by thin Hartmann layers. Four boundary conditions are required. To concentrate shear induced dissipation near the core-mantle boundary, these are taken to be: (i) no-slip at the core-mantle interface; (ii) geomagnetically estimated bulk westward flow at the base of the core-mantle boundary layer; (iii) no-slip at the inner-outer core interface; and, to describe magnetic locking of the inner core to the deep outer core, (iv) hydrodynamically stress-free at the inner-outer core boundary. By boldly assuming the axial core angular momentum anomaly to be zero, the super-rotation of the inner core is calculated to be at most 1.5 degrees per year.

Voorhies, C. V.↗
NASA NTRS2014

Variable Acceleration Force Calibration System (VACS)

Conventionally, force balances have been calibrated manually, using a complex system of free hanging precision weights, bell cranks, and/or other mechanical components. Conventional methods may provide sufficient accuracy in some instances, but are often quite complex and labor-intensive, requiring three to four man-weeks to complete each full calibration. To ensure accuracy, gravity-based loading is typically utilized. However, this often causes difficulty when applying loads in three simultaneous, orthogonal axes. A complex system of levers, cranks, and cables must be used, introducing increased sources of systematic error, and significantly increasing the time and labor intensity required to complete the calibration. One aspect of the VACS is a method wherein the mass utilized for calibration is held constant, and the acceleration is changed to thereby generate relatively large forces with relatively small test masses. Multiple forces can be applied to a force balance without changing the test mass, and dynamic forces can be applied by rotation or oscillating acceleration. If rotational motion is utilized, a mass is rigidly attached to a force balance, and the mass is exposed to a rotational field. A large force can be applied by utilizing a large rotational velocity. A centrifuge or rotating table can be used to create the rotational field, and fixtures can be utilized to position the force balance. The acceleration may also be linear. For example, a table that moves linearly and accelerates in a sinusoidal manner may also be utilized. The test mass does not have to move in a path that is parallel to the ground, and no re-leveling is therefore required. Balance deflection corrections may be applied passively by monitoring the orientation of the force balance with a three-axis accelerometer package. Deflections are measured during each test run, and adjustments with respect to the true applied load can be made during the post-processing stage. This paper will present the development and testing of the VASC concept.

Rhew, Ray D.↗
NASA NTRS1985

Magentically actuated compressor

A vibration free fluid compressor particularly adapted for Stirling cycle cryogenic refrigeration apparatus comprises a pair of identical opposing ferromagnetic pistons located in a housing and between a gas spring including a sealed volume of a working fluid such as gas under pressure. The gas compresses and expands in accordance with movement of the pistons to generate a compression wave which can be vented to other apparatus, for example, a displacer unit in a Stirling cycle engine. The pistons are urged outwardly due to the pressure of the gas; however, a fixed electromagnetic coil assembly located in the housing adjacent the pistons, is periodically energized to produce a magnetic field which interlinks the pistons in such a fashion that the pistons are mutually attracted to one another. The mass of the pistons, in conjunction with the compressed gas between them, form a naturally resonant system which, when the pistons are electromagnetically energized, produces an oscillating compression wave in the entrapped fluid medium.

Evans, J.↗
NASA NTRS1992

Time simulation of flutter with large stiffness changes

Time simulation of flutter, involving large local structural changes, is formulated with a state-space model that is based on a relatively small number of generalized coordinates. Free-free vibration modes are first calculated for a nominal finite-element model with relatively large fictitious masses located at the area of structural changes. A low-frequency subset of these modes is then transformed into a set of structural modal coordinates with which the entire simulation is performed. These generalized coordinates and the associated oscillatory aerodynamic force coefficient matrices are used to construct an efficient time-domain, state-space model for a basic aeroelastic case. The time simulation can then be performed by simply changing the mass, stiffness, and damping coupling terms when structural changes occur. It is shown that the size of the aeroelastic model required for time simulation with large structural changes at a few apriori known locations is similar to that required for direct analysis of a single structural case. The method is applied to the simulation of an aeroelastic wind-tunnel model. The diverging oscillations are followed by the activation of a tip-ballast decoupling mechanism that stabilizes the system but may cause significant transient overshoots.

Karpel, Mordechay↗
NASA NTRS1992

Time simulation of flutter with large stiffness changes

Time simulation of flutter, involving large local structural changes, is formulated with a state-space model that is based on a relatively small number of generalized coordinates. Free-free vibration modes are first calculated for a nominal finite-element model with relatively large fictitious masses located at the area of structural changes. A low-frequency subset of these modes is then transformed into a set of structural modal coordinates with which the entire simulation is performed. These generalized coordinates and the associated oscillatory aerodynamic force coefficient matrices are used to construct an efficient time-domain, state-space model for basic aeroelastic case. The time simulation can then be performed by simply changing the mass, stiffness and damping coupling terms when structural changes occur. It is shown that the size of the aeroelastic model required for time simulation with large structural changes at a few a priori known locations is similar to that required for direct analysis of a single structural case. The method is applied to the simulation of an aeroelastic wind-tunnel model. The diverging oscillations are followed by the activation of a tip-ballast decoupling mechanism that stabilizes the system but may cause significant transient overshoots.

Karpel, M.↗
NASA NTRS1994

Effects of turbulence compressibility and unsteadiness in compression corner flow

The structure of the separated flow region over a 20 degree compression corner at a free-stream Mach number of 2.84 is investigated computationally using a Reynolds averaged Navier Stokes (R.A.N.S.) solver and kappa-epsilon model. At this Mach number and ramp angle, a steady-state recirculation region of order delta(sub o) is observed, with onset of a 'plateau' in the wall pressure distribution near the corner. At lower ramp angles, separation is negligible, while at an angle of 24 degrees, separation regions of length 2 delta(sub o) are expected. Of interest here is the response of the mathematical model to inclusion of the pressure dilatation term for turbulent kinetic energy. Compared with the experimental data of Smits and Muck (1987), steady-state computations show improvement when the pressure dilatation term is included. Unsteady computations, using both unforced and then forced inlet conditions, did not predict the oscillatory motion of the separation bubble as observed in laboratory experiments. An analysis of the separation bubble oscillation and the turbulent boundary layer (T.B.L.) frequencies for this flow suggests that the bubble oscillations are of nearly the same order as the turbulent frequencies, and therefore difficult for the model to separate and resolve.

Brankovic, A.↗
NASA NTRS2014

Tropospheric Ozone Increases over the Southern Africa Region: Bellwether for Rapid Growth in Southern Hemisphere Pollution?

Increases in free-tropospheric (FT) ozone based on ozonesonde records from the early 1990s through 2008 over two subtropical stations, Irene (near Pretoria, South Africa) and Réunion (21 deg. S, 55 deg. E; approx. 2800 km NE of Irene in the Indian Ocean), have been reported. Over Irene a large increase in the urban-influenced boundary layer (BL, 1.5-4 km) was also observed during the 18-year period, equivalent to 30%decade−1. Here we show that the Irene BL trend is at least partly due to a gradual change in the sonde launch times from early morning to the midday period. The FT ozone profiles over Irene in 1990-2007 are re-examined, filling in a 1995-1999 gap with ozone profiles taken during the Measurements of Ozone by Airbus In-service Aircraft (MOZAIC) project over nearby Johannesburg. A multivariate regression model that accounts for the annual ozone cycle, El Niño-Southern Oscillation (ENSO) and possible tropopause changes was applied to monthly averaged Irene data from 4 to 11 km and to 1992-2011 Réunion sonde data from 4 to 15 km. Statistically significant trends appear predominantly in the middle and upper troposphere (UT; 4-11 km over Irene, 4-15 km over Réunion) in winter (June-August), with increases 1 ppbv yr(exp. −1) over Irene and approx. 2 ppbv yr(exp. −1) over Réunion. These changes are equivalent to approx. 25 and 35-45%decade( exp. −1), respectively. Both stations also display smaller positive trends in summer, with a 45%decade(exp. −1) ozone increase near the tropopause over Réunion in December. To explain the ozone increases, we investigated a time series of dynamical markers, e.g., potential vorticity (PV) at 330-350 K. PV affects UT ozone over Irene in November-December but displays little relationship with ozone over Réunion. A more likely reason for wintertime FT ozone increases over Irene and Réunion appears to be long-range transport of growing pollution in the Southern Hemisphere. The ozone increases are consistent with trajectory origins of air parcels sampled by the sondes and with recent NOx emissions trends estimated for Africa, South America and Madagascar. For Réunion trajectories also point to sources from the eastern Indian Ocean and Asia.

Trends↗
NASA NTRS1994

Screech tones from free and ducted supersonic jets

It is well known that screech tones from supersonic jets are generated by a feedback loop. The loop consists of three main components. They are the downstream propagating instability wave, the shock cell structure in the jet plume, and the feedback acoustic waves immediately outside the jet. Evidence will be presented to show that the screech frequency is largely controlled by the characteristics of the feedback acoustic waves. The feedback loop is driven by the instability wave of the jet. Thus the tone intensity and its occurrence are dictated by the characteristics of the instability wave. In this paper the dependence of the instability wave spectrum on the azimuthal mode number (axisymmetric or helical/flapping mode, etc.), the jet-to-ambient gas temperature ratio, and the jet Mach number are studied. The results of this study provide an explanation for the observed screech tone mode switch phenomenon (changing from axisymmetric to helical mode as Mach number increases) and the often-cited experimental observation that tone intensity reduces with increase in jet temperature. For ducted supersonic jets screech tones can also be generated by feedback loops formed by the coupling of normal duct modes to instability waves of the jet. The screech frequencies are dictated by the frequencies of the duct modes. Super resonance, resonance involving very large pressure oscillations, can occur when the feedback loop is powered by the most amplified instability wave. It is proposed that the observed large amplitude pressure fluctuations and tone in the test cells of Arnold Engineering Development Center were generated by super resonance. Estimated super-resonance frequency for a Mach 1.3 axisymmetric jet tested in the facility agrees well with measurement.

Tam, C. K. W.↗
NASA NTRS2002

Entropy Splitting for High Order Numerical Simulation of Compressible Turbulence

A stable high-order numerical scheme for direct numerical simulation (DNS) of shock-free compressible turbulence is presented. The method is applicable to general geometries. It contains no upwinding, artificial dissipation, or filtering. Instead the method relies on the stabilizing mechanisms of an appropriate conditioning of the governing equations and the use of compatible spatial difference operators for the interior points (interior scheme) as well as the boundary points (boundary scheme). An entropy-splitting approach splits the inviscid flux derivatives into conservative and nonconservative portions. The spatial difference operators satisfy a summation-by-parts condition, leading to a stable scheme (combined interior and boundary schemes) for the initial boundary value problem using a generalized energy estimate. A Laplacian formulation of the viscous and heat conduction terms on the right hand side of the Navier–Stokes equations is used to ensure that any tendency to odd–even decoupling associated with central schemes can be countered by the fluid viscosity. The resulting methods are able to minimize the spurious high-frequency oscillations associated with pure central schemes, especially for long time integration applications such as DNS. For validation purposes, the methods are tested in a DNS of compressible turbulent plane channel flow at low values of friction Mach number, where reference turbulence data bases exist. It is demonstrated that the methods are robust in terms of grid resolution, and in good agreement with published channel data. Accurate turbulence statistics can be obtained with moderate grid sizes. Stability limits on the range of the splitting parameter are determined from numerical tests.

N.D. Sandham↗