SEARCH · Search NASA
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.
A plasma instability which may lead to the thermal coupling of ions and electrons in two-temperature accretion flows on a time scale shorter than the Coulomb coupling time has been found. The processes involved with the instability are described. It is shown that the instability has time to grow only in regions of large-amplitude, small-scale MHD turbulence.
Flow induced instability force due to a labyrinth seal is analyzed. An approximate solution is given for the partial differential equation representing the flow in labyrinth seal and it is compared with the finite difference method in order to verify the accuracy of both methods. The effects of difference of inlet and outlet pressures of the seal, deflection of pressure and mass flow from the steady state, rotor diameter, seal clearance, seal interval and seal number on the flow induced force of the seal are investigated and it is known that some of these factors are very influential on the flow induced force.
The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.
We present observations of the formation, propagation, and decay of vortex-shaped features in coronal images from the Solar Dynamics Observatory associated with an eruption starting at about 2:30 UT on 2010 April 8. The series of vortices were formed along the interface between an erupting (dimming) region and the surrounding corona. They ranged in size from several to 10 arcsec and traveled along the interface at 6-14 km s{sup -1}. The features were clearly visible in six out of the seven different EUV wave bands of the Atmospheric Imaging Assembly. Based on the structure, formation, propagation, and decay of these features, we identified the event as the first observation of the Kelvin-Helmholtz instability (KHI) in the corona in EUV. The interpretation is supported by linear analysis and by a nonlinear 2.5-dimensional magnetohydrodynamic model of KHI. We conclude that the instability is driven by the velocity shear between the erupting and closed magnetic field of the coronal mass ejection. The shear-flow-driven instability can play an important role in energy transfer processes in coronal plasma.
Computational analyses of fluid flow through packed pebble bed domains using the Reynolds-averaged NavierStokes framework have had limited success in the past. Because of a lack of high-fidelity experimental or computational data, optimization of Reynolds-averaged closure models for these geometries has not been extensively developed. In the present study, direct numerical simulation was employed to develop a high-fidelity database that can be used for optimizing Reynolds-averaged closure models for pebble bed flows. A face-centered cubic domain with periodic boundaries was used. Flow was simulated at a Reynolds number of 9308 and cross-verified by using available quasi-DNS data. During the simulations, low-frequency instability modes were observed that affected the stationary solution. Furthermore, these instabilities were investigated by using the method of proper orthogonal decomposition, and a correlation was found between the time-dependent asymmetry of the averaged velocity profile data and the behavior of the highest energy eigenmodes.
We conduct a global, weakly nonlinear analysis of the magnetorotational instability (MRI) in a Taylor–Couette flow. This is a multiscale, perturbative treatment of the nonideal, axisymmetric MRI near threshold, subject to realistic radial boundary conditions and cylindrical geometry. We analyze both the standard MRI, initialized by a constant vertical background magnetic field, and the helical MRI, with an azimuthal background field component. This is the first weakly nonlinear analysis of the MRI in a global Taylor–Couette geometry, as well as the first weakly nonlinear analysis of the helical MRI. We find that the evolution of the amplitude of the standard MRI is described by a real Ginzburg–Landau equation (GLE), whereas the amplitude of the helical MRI takes the form of a complex GLE. This suggests that the saturated state of the helical MRI may itself be unstable on long spatial and temporal scales.
As demonstrated by IMAGE-EUV observations, plasmaspheric drainage plumes are a common feature emanating from the plasmasphere during periods of enhanced convection. Concurrent IMAGE-EUV and LANL-MPA observations reveal that plasmaspheric plumes are not a uniform density structure that extends out to geosynchronous orbit and beyond. Instead, in–situ observations reveal a large spatial-scale structure populated by small-scale density structures suggestive of the presence of a turbulent process. Here, we investigate the occurrence and nature of plasmaspheric drainage plumes as observed by concurrent EUV-MPA observations, as well as explore the possibility that the in–situ observed small-scale density structures are the signature of an instability produced by the sheared velocity found within the plasmaspheric layers separating the corotation of the main plasmasphere from the convection-driven flow generating the plume. Particle-in-cell simulations indicate that the generation of shear-flow driven instabilities is possible under plasmaspheric plasma conditions. The initial results of these simulations indicate that the plasmaspheric instabilities generated are electrostatic in nature, consist with in–situ plume observations.
We present a new technique for the investigation of shock-driven hydrodynamic phenomena in gases, liquids, and solids in arbitrary geometries. The technique consists of a pulsed power-driven resistive wire array explosion in combination with multi-MHz synchrotron radiography. Compared to commonly used techniques, it offers multiple advantages: (1) the shockwave geometry can be shaped to the requirements of the experiment, (2) the pressure (P > 300 MPa) generated by the exploding wires enables the use of liquid and solid hydrodynamic targets with well-characterized initial conditions (ICs), (3) the multi-MHz radiography enables data acquisition to occur within a single experiment, eliminating uncertainties regarding repeatability of the ICs and subsequent dynamics, and (4) the radiographic measurements enable estimation of compression ratios from the x-ray attenuation. In addition, the use of a synchrotron x-ray source allows the hydrodynamic samples to be volumetrically characterized at a high spatial resolution with synchrotron-based microtomography. This experimental technique is demonstrated by performing a planar Richtmyer–Meshkov instability (RMI) experiment on an aerogel–water interface characterized by Atwood number A0∼−0.8 and Mach number M∼1.5. The qualitative and quantitative features of the experiment are discussed, including the energy deposition into the exploding wires, shockwave generation, compression of the interface, startup phase of the instability, and asymptotic growth consistent with Richtmyer's impulsive theory. Additional effects unique to liquids and solids—such as cavitation bubbles caused by rarefaction flows or initial jetting due to small perturbations—are observed. It is also demonstrated that the technique is not shape dependent by driving a cylindrically convergent RMI experiment.
We analyze the local wave-number (LWN) model, a two-point spectral closure model for turbulence, as applied to the Rayleigh-Taylor (RT) instability, the flow induced by the relaxation of a statically-unstable density stratification. Model outcomes are validated against data from 3D simulations of the RT instability. In the first part of the study we consider the minimal model terms required to capture inhomogeneous mixing and show that this version, with suitable model coefficients, is sufficient to capture the evolution of important mean global quantities including mix-width, turbulent mass flux velocity, and Reynolds stress, if the start time is chosen such that the earliest transitions are avoided. However, this simple model does not permit the expected finite asymptote of the density-specific-volume covariance b. In the second part of the study, we investigate two forms for a source term for the evolution of the spectrum of density-specific-volume covariance for the LWN model. The first includes an empirically motivated calibration of the source to achieve the final asymptotic state of constant b. Here, the second form does not require calibration but, in conjunction with enhanced diffusion and drag captures the full evolution of all the dynamical quantities, namely, the mix-layer growth, turbulent mass-flux velocity, Reynolds stress, as well as the desired behavior of b.
Previous work has indicated that inelastic grains undergoing homogeneous cooling may be unstable, giving rise to the formation of velocity vortices and particle clusters for sufficiently large systems. Such instabilities are observed in industrial coal and biomass gasifiers and are known to influence gas-solid contact area, mixing dynamics, and heat/mass transfer rates. However, the driving mechanisms that lead to vortices and clusters are not well understood. Discrete-particle simulations provide a well-established method for understanding such mechanisms but are not a feasible technique for predicting the behavior of large-scale systems. Kinetic-theory-based continuum models offer an effective means of describing such flows, and instabilities present a stringent test of such models due to the transient, three-dimensional nature of instabilities and the large range of time and length scales over which these mechanisms occur.This work begins with the study, via a combination of continuum models and discrete- particle simulations, of a relatively simple flow and includes additional complexities in a stepwise manner to assess various driving mechanisms. Comparisons with discrete-particle simulations, which offer detailed, well-established (but computationally limited) descriptions of particle flows, indicate the ability of continuum models to accurately incorporate each mechanism. Specifically, the critical length scale for velocity vortices and/or particle clusters are studied via direct numerical simulation, molecular dynamics simulations, linear stability analyses of the continuum model, and transient simulations of the continuum model in a range of flow complexities, including moderate dissipation and particle concentration, frictional particles collisions, high gradients, and gas-solid flows. Strong agreement between kinetic-theory-based continuum models and discrete-particle simulations is found for a range for conditions. Furthermore, discrete-particle simulations offer insight into the relative importance of various mechanisms for instability in both granular and gas-solid flows. Normal inelasticity and friction are both found to be important in granular flows. Particle collisions, thermal drag, and mean drag are found to be important in gas-solid flows.
We propose a mechanism of state transition in wind-fed black hole (BH) binaries (high-mass X-ray binaries) such as Cyg X-1 and LMC X-1. Modeling a line-driven stellar wind from the companion by two-dimensional hydrodynamical calculations, we investigate the processes of wind capture by, and accretion onto, the BH. We assume that the wind acceleration is terminated at the He ii ionization front because ions responsible for line-driven acceleration are ionized within the front, i.e., the He iii region. It is found that the mass accretion rate inferred from the luminosity is remarkably smaller than the capture rate. Considering the difference, we construct a model for the state transition based on the accretion flow being controlled by magnetorotational instability. The outer flow is torus-like, and plays an important role to trigger the transition. The model can explain why state transition does occur in Cyg X-1, while not in LMC X-1. Cyg X-1 exhibits a relatively low luminosity, and then the He ii ionization front is located and can move between the companion and BH, depending on its ionizing photon flux. On the other hand, LMC X-1 exhibits too high luminosity for the front to move considerably; the front is too close to the companion atmosphere. The model also predicts that each state of high-soft or low-hard would last fairly long because the luminosity depends weakly on the wind velocity. In the context of the model, the state transition is triggered by a fluctuation of the magnetic field when its amplitude becomes comparable to the field strength in the torus-like outer flow.
The dust-acoustic waves and the stability theory for the permeating dusty plasma with power-law distributions are studied by using nonextensive q-statistics. In two limiting physical cases, when the thermal velocity of the flowing dusty plasma is much larger than, and much smaller than the phase velocity of the waves, we derived the dust-acoustic wave frequency, the instability growth rate, and the instability critical flowing velocity. As compared with the formulae obtained in part I [Gong et al., Phys. Plasmas 19, 043704 (2012)], all formulae of the present cases and the resulting plasma characteristics are q-dependent, and the power-law distribution of each plasma component of the permeating dusty plasma has a different q-parameter and thus has a different nonextensive effect. Further, we make numerical analyses of an example that a cometary plasma tail is passing through the interplanetary space dusty plasma and we show that these power-law distributions have significant effects on the plasma characteristics of this kind of plasma environment.
Tertiary modes in electrostatic drift-wave turbulence are localized near extrema of the zonal velocity U(x) with respect to the radial coordinate x. We argue that these modes can be described as quantum harmonic oscillators with complex frequencies, so their spectrum can be readily calculated. The corresponding growth rate $γ$ TI is derived within the modified Hasegawa-Wakatani model. We show that $γ$TI equals the primary-instability growth rate plus a term that depends on the local U''; hence, the instability threshold is shifted compared to that in homogeneous turbulence. This provides a generic explanation of the well-known yet elusive Dimits shift, which we find explicitly in the Terry-Horton limit. Linearly unstable tertiary modes either saturate due to the evolution of the zonal density or generate radially propagating structures when the shear |U'| is sufficiently weakened by viscosity. The Dimits regime ends when such structures are generated continuously.
With the advancing miniaturization of technological applications, processes on the mesoscale become increasingly important. This is the scale where the individual movement of particles transforms into cooperative behavior-behavior that cannot be explained by investigating the motion of individual particles alone.Complex plasmas are ideally suited to study the limits of cooperative behavior. The time scales of the dynamics of the microparticles embedded in the plasma are such that their movement can be fully resolved, and an investigation on the atomistic (kinetic) level is possible. In addition, complex plasmas can be considered a model system for ordinary fluids: The internal microparticle dynamics is basically undamped and is characterized by the similarity parameters matching those of other fluids. This similarity does not break down even at small scales: For instance, in [2], microparticle droplets comprised of only a few 1000-10000 particles were examined. In these experiments, the Weber number (the ratio of inertia to surface tension forces) matches that of falling water drops. As another example, the onset of a Rayleigh-Taylor instability in a complex plasma can be described by the ordinary dispersion relation, even at scales of only few particle layers. This allows investigating the 'nanoscale' of fluid flows, and, hence, the limits of cooperative behavior.
Various secondary and tertiary instabilities in plane channael flow are explored via time-dependent numerical simulations using the incompressible Navier-Stokes equations. Comparisons are made between transitional flows at Reynolds numbers 1500, 5000, and 8000. The lambda vortex, detached shear layer, and inverted vortex regions are identified, and the origin of the latter is explained. The laminar breakdown of the Re = 1500 flow is computed with high resolution, and the nature of its ensuing hairpin eddies is clarified by numerical particle paths. The potential of center-mode, rather than wall-mode transitions is proposed, and the resulting flow structure is described.
Topics covered include: Flow Physics and control for Internal and External Aerodynamics (not in TOC...starts on pg13); Breaking CFD Bottlenecks in Gas-Turbine Flow-Path Design; Streamwise Vortices on the Convex Surfaces of Circular Cylinders and Turbomachinery Blading; DNS and Embedded DNS as Tools for Investigating Unsteady Heat Transfer Phenomena in Turbines; Cavitation, Flow Structure and Turbulence in the Tip Region of a Rotor Blade; Development and Application of Plasma Actuators for Active Control of High-Speed and High Reynolds Number Flows; Active Flow Control of Lifting Surface With Flap-Current Activities and Future Directions; Closed-Loop Control of Vortex Formation in Separated Flows; Global Instability on Laminar Separation Bubbles-Revisited; Very Large-Scale Motions in Smooth and Rough Wall Boundary Layers; Instability of a Supersonic Boundary-Layer With Localized Roughness; Active Control of Open Cavities; Amplitude Scaling of Active Separation Control; U.S. Air Force Research Laboratory's Need for Flow Physics and Control With Applications Involving Aero-Optics and Weapon Bay Cavities; Some Issues Related to Integrating Active Flow Control With Flight Control; Active Flow Control Strategies Using Surface Pressure Measurements; Reduction of Unsteady Forcing in a Vaned, Contra-Rotating Transonic Turbine Configuration; Active Flow Control Stator With Coanda Surface; Controlling Separation in Turbomachines; Flow Control on Low-Pressure Turbine Airfoils Using Vortex Generator Jets; Reduced Order Modeling Incompressible Flows; Study and Control of Flow Past Disk, and Circular and Rectangular Cylinders Aligned in the Flow; Periodic Forcing of a Turbulent Axisymmetric Wake; Control of Vortex Breakdown in Critical Swirl Regime Using Azimuthal Forcing; External and Turbomachinery Flow Control Working Group; Boundary Layers, Transitions and Separation; Efficiency Considerations in Low Pressure Turbines; Summary of Conference; and Final Plenary Session Transcript.
The interaction of the drift wave (DW) turbulence and zonal flow (ZF) is investigated with the modified Hasegawa–Mima equation taking into account the backreaction of ZF velocity on DW turbulence. It is shown that the y-averaged enstrophy of DW turbulence and the velocity of ZF are intrinsically related. By utilizing this feature, a nonlinear stage of DW modulational instability is considered within the framework of the wave kinetic equation. It is shown that in this approximation, the nonlinear stage of the modulational instability results in the collapsing solutions, accompanied by the “wave breaking” phenomenon. Numerical simulations based on the Hasegawa–Mima equation show that for a weak DW turbulence, Φ ̃ = ( e φ ̃ / T e ) ( L n / ρ s ) ⪝ 1 , the collapsing-like features on both ZF and y-averaged enstrophy of DW turbulence decay in time and then re-emerge again at different locations. For the case of a strong DW turbulence, Φ ̃ > 1 , where nonlinear interactions of DW harmonics dominate, stable spatial structures of ZF and y-averaged enstrophy of DW turbulence emerge.
A mechanism for explosive magnetic reconnection is investigated by analyzing the nonlinear evolution of a collisionless tearing mode in a two-fluid model that includes the effects of electron inertia and temperature. These effects cooperatively enable a fast reconnection by forming an X-shaped current-vortex layer centered at the reconnection point. A high-resolution simulation of this model for an unprecedentedly small electron skin depth d{sub e} and ion-sound gyroradius ρ{sub s}, satisfying d{sub e}=ρ{sub s}, shows an explosive tendency for nonlinear growth of the tearing mode, where it is newly found that the explosive widening of the X-shaped layer occurs locally around the reconnection point with the length of the X shape being shorter than the domain length and the wavelength of the linear tearing mode. The reason for the onset of this locally enhanced reconnection is explained theoretically by developing a novel nonlinear and nonequilibrium inner solution that models the local X-shaped layer, and then matching it to an outer solution that is approximated by a linear tearing eigenmode with a shorter wavelength than the domain length. This theoretical model proves that the local reconnection can release the magnetic energy more efficiently than the global one and the estimated scaling of the explosive growth rate agrees well with the simulation results.