SEARCH · Engineering Papers
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.
The simple procedure presented for spatial stability computations can substantially reduce the computational requirements of such analyses, as illustrated for the cases of both internal and external cases of compressible and incompressible flows, and both viscous and inviscid instability modes. Excellent estimates of spatial eigenvalues are obtained.
The fluid dynamics of a high-speed jet are governed by the instability waves that form in the free-shear boundary layer of the jet. Jet excitation manipulates the growth and saturation of particular instability waves to control the unsteady flow structures that characterize the energy cascade in the jet.The results may include jet noise mitigation or a reduction in the infrared signature of the jet. The Localized Arc Filament Plasma Actuators (LAFPA) have demonstrated the ability to excite a high-speed jets in laboratory experiments. Extending and optimizing this excitation technology, however, is a complex process that will require many tests and trials. Computational simulations can play an important role in understanding and optimizing this actuator technology for real-world applications. Previous research has focused on developing a suitable actuator model and coupling it with the appropriate computational fluid dynamics (CFD) methods using two-dimensional spatial flow approximations. This work is now extended to three-dimensions (3-D) in space. The actuator model is adapted to a series of discrete actuators and a 3-D LES simulation of an excited jet is run. The results are used to study the fluid dynamics near the actuator and in the jet plume.
Two dual-stream, heated jet, Compact Jet Engine Simulator (CJES) units are designed for wind tunnel acoustic experiments involving a Hybrid Wing Body (HWB) vehicle. The newly fabricated CJES units are characterized with a series of acoustic and flowfield investigations to ensure successful operation with minimal rig noise. To limit simulator size, consistent with a 5.8% HWB model, the CJES units adapt Ultra Compact Combustor (UCC) technology developed at the Air Force Research Laboratory. Stable and controllable operation of the combustor is demonstrated using passive swirl air injection and backpressuring of the combustion chamber. Combustion instability tones are eliminated using nonuniform flow conditioners in conjunction with upstream screens. Through proper flow conditioning, rig noise is reduced by more than 20 dB over a broad spectral range, but it is not completely eliminated at high frequencies. The low-noise chevron nozzle concept designed for the HWB test shows expected acoustic benefits when installed on the CJES unit, and consistency between CJES units is shown to be within 0.5 dB OASPL.
An automated cloud tracking algorithm is applied to Cassini Imaging Science Subsystem high-resolution apoapsis images of Saturn from 2005 and 2007 and moderate resolution images from 2011 and 2012 to define the near-global distribution of zonal winds and eddy momentum fluxes at the middle troposphere cloud level and in the upper troposphere haze. Improvements in the tracking algorithm combined with the greater feature contrast in the northern hemisphere during the approach to spring equinox allow for better rejection of erroneous wind vectors, a more objective assessment at any latitude of the quality of the mean zonal wind, and a population of winds comparable in size to that available for the much higher contrast atmosphere of Jupiter. Zonal winds at cloud level changed little between 2005 and 2007 at all latitudes sampled. Upper troposphere zonal winds derived from methane band images are approx. 10 m/s weaker than cloud level winds in the cores of eastward jets and approx. 5 m/s stronger on either side of the jet core, i.e., eastward jets appear to broaden with increasing altitude. In westward jet regions winds are approximately the same at both altitudes. Lateral eddy momentum fluxes are directed into eastward jet cores, including the strong equatorial jet, and away from westward jet cores and weaken with increasing altitude on the flanks of the eastward jets, consistent with the upward broadening of these jets. The conversion rate of eddy to mean zonal kinetic energy at the visible cloud level is larger in eastward jet regions (5.2x10(exp -5) sq m/s) and smaller in westward jet regions (1.6x10(exp -5) sqm/s) than the global mean value (4.1x10(ep -5) sq m/s). Overall the results are consistent with theories that suggest that the jets and the overturning meridional circulation at cloud level on Saturn are maintained at least in part by eddies due to instabilities of the large-scale flow near and/or below the cloud level.
Surface waves at the magnetopause flanks typically feature steeper, i.e., more inclined leading (antisunward facing) than trailing (sunward facing) edges. This is expected for Kelvin-Helmholtz instability (KHI) amplified waves. Very rarely, during northward interplanetary magnetic field (IMF) conditions, anomalous inverse steepening has been observed. The small-scale tetrahedral configuration of the Magnetospheric Multiscale spacecraft and their high time resolution measurements enable us to routinely ascertain magnetopause boundary inclinations during surface wave passage with high accuracy by four-spacecraft timing analysis. At the dusk flank magnetopause, 77%/23% of the analyzed wave intervals exhibit regular inverse steepening. Inverse steepening happens during northward IMF conditions, as previously reported and, in addition, during intervals of dominant equatorial IMF. Inverse steepening observed under the latter conditions may be due to the absence of KHI or due to instabilities arising from the alignment of flow and magnetic fields in the magnetosheath.
The laminar boundary layer which forms on a sharp cone in a supersonic freestream, where lateral curvature plays a key role in the physics of the problem is considered. This flow is then analyzed from the point of view of linear, temporal, inviscid stability. The basic, non-axisymmetric disturbance equations are derived for general flows of this class, and a so called triply generalized inflexion condition is found for the existence of subsonic neutral modes of instability. This condition is analogous to the well-known generalized inflexion condition found in planar flows, although in the present case the condition depends on both axial and aximuthal wavenumbers. Extensive numerical results are presented for the stability problem at a freestream Mach number of 3.8, for a range of streamwise locations. These results reveal that a new mode of instability may occur, peculiar to flows of this type involving curvature. Additionally, asymptotic analyses valid close to the tip of the cone, far downstream of the cone are presented, and these give a partial (asymptotic) description of this additional mode of instability.
In nature and engineering applications, the Rayleigh–Taylor instability (RTI) occurs over a wide range of Atwood, Reynolds, Mach, and Knudsen numbers. At low Atwood, Mach, and Knudsen numbers, the classic advective instability causes quasi-symmetric bubble and spike growth on the two sides of the interface. However, recent findings suggest that at high degrees of rarefaction, advective effects are suppressed and molecular diffusion leads to planar growth of the density fronts on either side of the interface. This study aims to investigate the flow physics of the transition from advective to diffusive behavior, focusing on the onset of the kinetic effects. Using the gas kinetic methodology, RTI is simulated over a range of Knudsen and Mach numbers in the transition regime. The simulation results reveal the various stages of transformation from advective instability to diffusive transport. For the first time, the study demonstrates the existence of a Knudsen–Mach parameter regime where the bubble side exhibits advective instability, while the other side shows a planar density front due to molecular diffusion, rather than the canonical advective spike shape. The dominance of different mechanisms on the two sides of the interface leads to the advective–diffusive asymmetry. In conclusion, the findings of this study can lead to a more comprehensive understanding of RTI over a wide range of Mach and Knudsen numbers.
We examine the excitation of inviscid stationary crossflow instabilities near an isolated surface hump (or indentation) underneath a three-dimensional boundary layer. As the hump height (or indentation depth) is increased from zero, the receptivity process becomes nonlinear even before the stability characteristics of the boundary layer are modified to a significant extent. This behavior contrasts sharply with earlier findings on the excitation of the lower branch Tollmien-Schlichting modes and is attributed to the inviscid nature of the crossflow modes, which leads to a decoupling between the regions of receptivity and stability. As a result of this decoupling, similarity transformations exist that allow the nonlinear receptivity of a general three-dimensional boundary layer to be studied with a set of canonical solutions to the viscous sublayer equations. The parametric study suggests that the receptivity is likely to become nonlinear even before the hump height becomes large enough for flow reversal to occur in the canonical solution. We also find that the receptivity to surface humps increases more rapidly as the hump height increases than is predicted by linear theory. On the other hand, receptivity near surface indentations is generally smaller in comparison with the linear approximation. Extension of the work to crossflow receptivity in compressible boundary layers and to Gortler vortex excitation is also discussed.
Computations are performed to study the flow past an isolated, spanwise symmetric roughness element in zero pressure gradient boundary layers at Mach 3.5 and 5.9, with an emphasis on roughness heights of less than 55 percent of the local boundary layer thickness. The Mach 5.9 cases include flow conditions that are relevant to both ground facility experiments and high altitude flight ("cold wall" case). Regardless of the Mach number, the mean flow distortion due to the roughness element is characterized by long-lived streamwise streaks in the roughness wake, which can support instability modes that did not exist in the absence of the roughness element. The higher Mach number cases reveal a variety of instability mode shapes with velocity fluctuations concentrated in different localized regions of high base flow shear. The high shear regions vary from the top of a mushroom shaped structure characterizing the centerline streak to regions that are concentrated on the sides of the mushroom. Unlike the Mach 3.5 case with nearly same values of scaled roughness height k/delta and roughness height Reynolds number Re(sub kk), the odd wake modes in both Mach 5.9 cases are significantly more unstable than the even modes of instability. Additional computations for a Mach 3.5 boundary layer indicate that the presence of a roughness element can also enhance the amplification of first mode instabilities incident from upstream. Interactions between multiple roughness elements aligned along the flow direction are also explored.
The stability of the parallel flow between a vertical crystal-melt interface and a vertical wall held at a temperature above the melting point of the crystal is analyzed for Prandtl numbers, P, ranging from 0.01 to 100. Three modes of instability occur: (1) a buoyant mode, (2) a shear mode, and (3) a coupled crystal-melt mode. The buoyant and shear modes are similar to those that occur for flow between two vertical rigid walls held at different temperatures. For Prandtl numbers greater than approximately two, the coupled crystal-melt mode occurs at a lower Grashof number than the other two modes. Specific results are given for succinonitrile (P = 22.8) and lead (P = 0.0225). These calculations and similar calculations for a cylindrical geometry were motivated by and are in general agreement with recent experiments on succinonitrile.
An adaptive control algorithm has been developed for the suppression of combustion thermo-acoustic instabilities. This technique involves modulating the fuel flow in the combustor with a control phase that continuously slides within the stable phase region, in a back and forth motion. The control method is referred to as Adaptive Sliding Phasor Averaged Control (ASPAC). The control method is evaluated against a simplified simulation of the combustion instability. Plans are to validate the control approach against a more physics-based model and an actual experimental combustor rig.
Transient behavior of flare-associated solar wind in the nonradial open field region is numerically investigated, taking into account the thermal and dynamical coupling between the chromosphere and the corona. A realistic steady solar wind is constructed which passes through the inner X-type critical point in the rapidly diverging region. The wind speed shows a local maximum at the middle, O-type, critical point. The wind's density and pressure distributions decrease abruptly in the rapidly diverging region of the flow tube. The transient behavior of the wind following flare energy deposition includes ascending and descending conduction fronts. Thermal instability occurs in the lower corona, and ascending material flows out through the throat after the flare energy input ceases. A local density distribution peak is generated at the shock front due to the pressure deficit just behind the shock front.
Analyses of the x ray surface brightness profiles of cluster cooling flows suggest that the mass flow rate decreases towards the center of the cluster. It is often suggested that this decrease results from thermal instabilities, in which denser blobs of gas cool rapidly and drop below x ray emitting temperatures. If the seeds for the thermal instabilities are entropy perturbations, these perturbations must enter the flow already in the nonlinear regime. Otherwise, the blobs would take too long to cool. Here, researchers suggest that such nonlinear perturbations might start as blobs of interstellar gas which are stripped out of cluster galaxies. Assuming that most of the gas produced by stellar mass loss in cluster galaxies is stripped from the galaxies, the total rate of such stripping is roughly M sub Interstellar Matter (ISM) approx. 100 solar mass yr(-1). It is interesting that the typical rates of cooling in cluster cooling flows are M sub cool approx. 100 solar mass yr(-1). Thus, it is possible that a substantial portion of the cooling gas originates as blobs of interstellar gas stripped from galaxies. The magnetic fields within and outside of the low entropy perturbations can help to maintain their identities, both by suppressing thermal conduction and through the dynamical effects of magnetic tension. One significant question concerning this scenario is: Why are cooling flows seen only in a fraction of clusters, although one would expect gas stripping to be very common. It may be that the density perturbations only survive and cool efficiently in clusters with a very high intracluster gas density and with the focusing effect of a central dominant galaxy. Inhomogeneities in the intracluster medium caused by the stripping of interstellar gas from galaxies can have a number of other effects on clusters. For example, these density fluctuations may disrupt the propagation of radio jets through the intracluster gas, and this may be one mechanism for producing Wide-Angle-Tail radio galaxies.
Understanding gas migration in compacted clay materials, e.g., bentonite and claystone, is important for the design and performance assessment of an engineered barrier system of a radioactive waste repository system, as well as many practical applications. Existing field and laboratory data on gas migration processes in low-permeability clay materials demonstrate the complexity of flow and transport processes, including various types of instabilities, caused by nonlinear dynamics of coupled processes of liquid–gas exchange, dilation, fracturing, fracture healing, etc., which cannot be described by classical models of fluid dynamics in porous media. We here show that the complexity of gas migration processes can be explained using a phenomenological concept of nonlinear dynamics and deterministic chaos theory. To do so, we analyzed gas pressure and gas influx (i.e., input) and outflux (i.e., output), recorded during the gas injection experiment in the compact Mx80-D bentonite sample, and calculated a set of the diagnostic parameters of nonlinear dynamics and chaos, such a global embedding dimension, a correlation dimension, an information dimension, and a spectrum of Lyapunov exponents, as well as plotted 2D and 3D pseudo-phase-space strange attractors, based on the univariate influx and outflux time series data. These results indicate the presence of phenomena of low-dimensional deterministic chaotic behavior of gas migration in bentonite. In particular, during the onset of gas influx in the bentonite core, before the breakthrough, the development of gas flow pathways is characterized by the process of chaotic gas diffusion. After the breakthrough, with inlet-to-outlet movement of gas, the prevailing process is chaotic advection. During the final phase of the experiment, with no influx to the sample, the relaxation pattern of gas outflux is resumed back to a process of chaotic diffusion. The types of data analysis and a proposed phenomenological model can be used to establish the basic principles of experimental data-gathering, modeling predictions, and a research design.
Even though slow convective flow is difficult to detect in solidifying metals, it can readily be observed in transparent materials by observing the motion of small neutrally buoyant particles. Succinonitrile, which solidifies with an unfaceted solid/liquid interface and has well characterized physical properties, is considered an excellent material for such studies. For studies of solute-induced convection, ethanol is a useful addition to succinonitrile since it has a lower density and a somewhat similar molecular structure. Samples of high purity and ethanol-doped succinonitrile are unidirectionally solidified in a vertical temperature gradient. Latex mimcrospheres 2 microns in diameter are suspended in the liquid to reveal the convective flow. Convective and morphological stability is observed as a function of solute concentration and growth velocity. These measurements are compared with theoretical calculations that predict the transition from stability to instability as a function of solidification conditions. The predicted transitions occur at low concentrations and solidification velocities; for this reason, extreme care must be taken in order to eliminate the effects of impurities or thermally induced convection.
The Dimits shift is the shift between the threshold of the drift-wave primary instability and the actual onset of turbulent transport in magnetized plasma. It is generally attributed to the suppression of turbulence by zonal flows, but developing a more detailed understanding calls for consideration of specific reduced models. The modified Terry--Horton system has been proposed by St-Onge [J. Plasma Phys. {\bf 83}, 905830504 (2017)] as a minimal model capturing the Dimits shift. Here, we use this model to develop an analytic theory of the Dimits shift and a related theory of the tertiary instability of zonal flows. We show that tertiary modes are localized near extrema of the zonal velocity $U(x)$, where $x$ is the radial coordinate. By approximating $U(x)$ with a parabola, we derive the tertiary-instability growth rate using two different methods and show that the tertiary instability is essentially the primary drift-wave instability modified by the local $U''$. Then, depending on $U''$, the tertiary instability can be suppressed or unleashed. The former corresponds to the case when zonal flows are strong enough to suppress turbulence (Dimits regime), while the latter corresponds to the case when zonal flows are unstable and turbulence develops. This understanding is different from the traditional paradigm that turbulence is controlled by the flow shear $U'$. Our analytic predictions are in agreement with direct numerical simulations of the modified Terry--Horton system.
The stability of the flow between two vertical, infinite, rigid, coaxial cylinders held at different temperatures is analyzed by linear stability theory. For a Prandtl number of 22.8 and a radius ratio of 0.02, the flow is unstable to an axisymmetric perturbation at a critical Grashof number of 2150; the wave speed of the instability is comparable to the maximum velocity of the unperturbed flow. When the rigid outer cylinder-fluid interface is replaced by a crystal-melt interface which can change shape, two new modes of instability occur at lower Grashof numbers.
The two-dimensional steady and time dependent properties of plane Poiseuille and plane Couette flows are analyzed using iterative techniques and full numerical simulation of the Navier-Stokes equations. It is shown that the finite-amplitude two-dimensional states investigated are strongly unstable to very small three-dimensional perturbations. It is also shown, through full numerical simulation, that this explosive secondary instability can explain the subcritical transitions that occur in real flows. Finally, it is shown that the three-dimensional instability can be analyzed by a linear stability analysis of a two-dimensional flow consisting of the basic parallel flow and a steady (or quasi-steady) finite-amplitude two-dimensional cellular motion.