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 flow over wave shaped surfaces is studied using a Navier Stokes solver. Detailed comparisons with theoretical results are presented, including the stability of a laminar flow over wavy surfaces. Drag characteristics of nonplanar surfaces are predicted using the Navier-Stokes solver. The secondary instabilities of wall bounded and free shear flows are also discussed.
One-dimensional oscillating flow models are very useful for designing pulse tubes. They are simple to use, not computationally intensive, and the physical relationship between temperature, pressure and mass flow are easy to understand when used in conjunction with phasor diagrams. They do not possess, however, the ability to directly calculate thermal and momentum diffusion in the direction transverse to the oscillating flow. To account for transverse effects, lumped parameter corrections, which are obtained though experiment, must be used. Or two-dimensional solutions of the differential fluid equations must be obtained. A linear two-dimensional solution to the fluid equations has been obtained. The solution provides lumped parameter corrections for one-dimensional models. The model accounts for heat transfer and shear flow between the gas and the tube. The complex Nusselt number and complex shear wall are useful in describing these corrections, with phase relations and amplitudes scaled with the Prandtl and Valensi numbers. The calculated ratio, a, between a two-dimensional solution of the oscillating temperature and velocity and a one-dimensional solution for the same shows a scales linearly with Va for Va less than 30. In this region alpha less than 0.5, that is, the enthalpy flow calculated with a two-dimensional model is 50% of a calculation using a one-dimensional model. For Va greater than 250, alpha = 0.8, showing that diffusion is still important even when it is confined to a thing layer near the tube wall.
This work will enable improved understanding of tokamak edge flows and impurity-species dynamics and their impact on MHD stability as it relates to ELMs, RMPs and QH-mode by leveraging computational developments of the NIMROD code. Experimentally, it is established that the tokamak flows and impurities species can have a substantial effect on MHD stability. This is particularly true for the edge plasma where, for example, large flow shear is correlated with the occurrence of the Quiescent H-mode (QH- mode) state as opposed to operation with edge localized modes (ELMs). Understanding plasma flows is challenging because of the many physical effects that come into play: particle orbits loss, neutrals dynamics, and interaction between multiple ion species with decoupled motion. We will develop a model that incorporates multiple ion species into the MHD framework that includes self-consistent magnetic-field evolution. With this model we will understand how the presence of impurities impacts the tokamak edge-pedestal flows. Finally, we will quantify the impact of the new multiple species modeling on 3D QH-mode simulations. This work will have broader impacts on all MHD studies when multiple species are present (as is the case in modern tokamaks), as well as contributing to the edge and transport communities who also are interested in the studies of multiple-ion-species flows in the edge pedestal region.
A method of identification of coherent structures in a random turbulent flow is outlined with emphasis on the structure best correlated with a random energy-integrable field and application of the proper orthogonal decomposition to the shear flow of the wall region. The wall region of a pipe flow with almost pure glycerine is investigated experimentally, dynamic equations are obtained, and an energy-transfer model is employed. Numerical solutions are obtained, implications for the flow in the wall region are discussed with focus placed on the sets of solutions exhibiting intermittency. The most significant finding of the study is considered to be the suggestion of the etiology of the bursting phenomenon; presuming that abrupt transitions from one fixed point to the other can be identified with a burst, these bursts appear to be produced autonomously by the wall region, but be triggered by pressure signals from the outer layer.
A new sub-grid mixing model for use in large eddy simulations of turbulent combustion is presented and applied to a hydrogen-air diffusion flame. The sub-grid model is based on Kerstein's Linear Eddy Model (1988). The model is first used to predict the mixing of a conserved scalar in a turbulent shear flow. The model correctly predicts the behavior of the pdf of the scalar field and displays a nonmarching peak at the preferred mixture fraction as the shear layer is traversed. It is then illustrated how a reduced chemical mechanism can be implemented within the linear eddy subgrid model formulation. The model is used to predict NO formation in hydrogen-air diffusion flame using a reduced chemical mechanism involving nine reactive scalars.
In investigating the potential of a new actuator for use in an active flow control system, several objectives had to be accomplished, the largest of which was the experimental setup. The work was conducted at the NASA Langley 20x28 Shear Flow Control Tunnel. The actuator named Thunder, is a high deflection piezo device recently developed at Langley Research Center. This research involved setting up the instrumentation, the lighting, the smoke, and the recording devices. The instrumentation was automated by means of a Power Macintosh running LabVIEW, a graphical instrumentation package developed by National Instruments. Routines were written to allow the tunnel conditions to be determined at a given instant at the push of a button. This included determination of tunnel pressures, speed, density, temperature, and viscosity. Other aspects of the experimental equipment included the set up of a CCD video camera with a video frame grabber, monitor, and VCR to capture the motion. A strobe light was used to highlight the smoke that was used to visualize the flow. Additional effort was put into creating a scale drawing of another tunnel on site and a limited literature search in the area of active flow control.
A three-dimensional numerical simulation of plasma density irregularities in the postsunset equatorial F region ionosphere leading to equatorial spread F (ESF) is described. The simulation evolves under realistic background conditions including bottomside plasma shear flow and vertical current. It also incorporates C/NOFS satellite data which partially specify the forcing. A combination of generalized Rayleigh-Taylor instability (GRT) and collisional shear instability (CSI) produces growing waveforms with key features that agree with C/NOFS satellite and ALTAIR radar observations in the Pacific sector, including features such as gross morphology and rates of development. The transient response of CSI is consistent with the observation of bottomside waves with wavelengths close to 30 km, whereas the steady state behavior of the combined instability can account for the 100+ km wavelength waves that predominate in the F region.
The Ring-Sheared Drop module was conceived as a centimeter-scale, containerless biochemical reactor to study shearing flow at gas-liquid interfaces without complications associated with solid walls. The drop is constrained by a thin, stationary contact ring in one hemisphere and sheared by the steady rotation of another ring in the other hemisphere. The original experiment considers a 2.5 cm diameter drop of human insulin solution in order to examine the formation of amyloid fibrils without wall nucleation. Shear-induced amyloidogenesis is of wide interest due to the pathogenic role amyloid proteins play in diseases such as Alzheimer's, Parkinson's, prion diseases, and type-2 diabetes. Pinning full-scale drops of native protein solutions - prior to any fibrillization - was challenging and required a collaborative effort. In this talk we report on laboratory experiments, multiphase CFD simulations, and several parabolic flights to find strategies that are effective for constraining large drops of protein in microgravity.
The basic objective of this research is to extend the capabilities of Large Eddy Simulations (LES) and Direct Numerical Simulations (DNS) for the computational analyses of high speed reacting flows. In the efforts related to LES, we were primarily involved with assessing the performance of the various modern methods based on the Probability Density Function (PDF) methods for providing closures for treating the subgrid fluctuation correlations of scalar quantities in reacting turbulent flows. In the work on DNS, we concentrated on understanding some of the relevant physics of compressible reacting flows by means of statistical analysis of the data generated by DNS of such flows. In the research conducted in the second year of this program, our efforts focused on the modeling of homogeneous compressible turbulent flows by PDF methods, and on DNS of non-equilibrium reacting high speed mixing layers. Some preliminary work is also in progress on PDF modeling of shear flows, and also on LES of such flows.
A wave-1 minor warming is simulated in a mechanistic, global, primitive equation model. The modification of the zonal flow by the rearrangement of potential vorticity on a middle stratopheric isentrope is compared in a fully nonlinear model and in a model with only one wave and the zonal flow (a quasi-linear model). The permanent rearrangement of potential vorticity during the wave episode is more intense and more localized meridionally in the fully nonlinear model, which is able to capture the process of planetary wave breaking in some detail. Additional experiments reveal that the differences between the quasi-linear and nonlinear models persist for a broad range of wave amplitudes, and that the quasi-linear model can qualitatively reproduce the modification of the zonal flow by the wave when the diffusive dissipation of the wave is enhanced. These results are discussed in the context of the theory of barotropic Rossby waves in shear flows, and in comparison with recent numerical simulations of the middle atmosphere.
One of the long-standing uncertainties in the wave-resonance theory of coronal heating is the stability of the resonance layer. The wave motions in the resonance layer produce highly localized shear flows which vary sinusoidally in time with the resonance period. This configuration is potentially susceptible to the Kelvin-Helmholtz instability (KHI), which can enhance small-scale structure and turbulent broadening of shear layers on relatively rapid ideal timescales. We have investigated numerically the response of a characteristic velocity profile, derived from resonance absorption models, to finite fluid perturbations comparable to photospheric fluctuations. We find that the KHI primarily should affect long (approximately greater than 6 x 10(exp 4) km) loops where higher velocity flows (M approximately greater than 0.2) exist in resonance layers of order 100 km wide. There, the Kelvin-Helmholtz growth time is comparable to or less than the resonance quarter-period, and the potentially stabilizing magnetic effects are not felt until the instability is well past the linear growth stage. Not only is the resonance layer broadened by the KHI, but also the convective energy transport out of the resonance layer is increased, thus adding to the efficiency of the wave-resonance heating process. In shorter loops, e.g., those in bright points and compact flares, the stabilization due to the magnetic field and the high resonance frequency inhibit the growth of the Kelvin-Helmholtz instability beyond a minimal level.
In Quiescent H-mode (QH-mode), edge-localized modes (ELMs) are naturally replaced by a low-n edge harmonic oscillation (EHO), yet the self-regulating transport mechanism driven by the EHO remains insufficiently understood. Using high-spatiotemporal-resolution imaging diagnostics on KSTAR—electron cyclotron emission imaging and broadband electron cyclotron emission—we resolve the eigenmode structure of the EHO and elucidate its regulatory role in edge transport. The EHO is localized within the pedestal near the maximum pressure gradient, propagates in the ion-diamagnetic direction, and its radial envelope expands with increasing shear, suggesting that rotational shear is closely associated with the structural evolution of the EHO. Information-theoretic Transfer Entropy analysis identifies a distinct ‘dual-stabilization’ regulatory interaction pattern: the EHO is associated with enhanced outward energy transport to limit the pressure gradient, while the background shear flow is associated with preferential regulatory influence on the EHO saturation amplitude and energy flux. These results demonstrate that the QH-mode pedestal is sustained by a shear-associated regulation mechanism, where the coupling between rotational shear and mode structure appears closely linked to the sustained ELM-free state.
Self-generating magnetic field research is performed through the operation of the newly built Madison Plasma Dynamo Experiment (MPDX). MPDX is designed to explore a hitherto unexplored part of parameter space where dynamos operate in nature. Dynamos are systems which continuously transform kinetic energy from plasma flow into magnetic energy. Before this project, MPDX had demonstrated the ability to create steady-state, unmagnetized high conductivity, low viscosity plasmas in the lab and shown that sheared flows can be driven from the plasma boundary using an electrostatic stirring technique; the resulting Reynolds numbers and Alfvén Mach numbers easily match those required for our dynamo scenarios. The research activities included: (1) systematic undertaking of a set of plasma hydrodynamics experiments, exploring several different flow geometries, measuring plasma viscosity, how toroidal and poloidal flows self-consistently interact, observing the transition to turbulence and experimentally understanding the boundary conditions (on flow and magnetic field) that make this system unique, (2) searching for slow, laminar dynamos, and (3) searching for fast, turbulent dynamos. Discovering the conditions under which dynamos self-generate magnetic fields and then understanding how this field changes plasma dynamics is one of the most compelling questions in all of plasma astrophysics.
A newly developed non-isotropic multiple scale turbulence model (MS/ASM) is described for complex flow calculations. This model focuses on the direct modeling of Reynolds stresses and utilizes split-spectrum concepts for modeling multiple scale effects in turbulence. Validation studies on free shear flows, rotating flows and recirculating flows show that the current model perform significantly better than the single scale k-epsilon model. The present model is relatively inexpensive in terms of CPU time which makes it suitable for broad engineering flow applications.
The principles of two flow control strategies, those of Huebler (Luescher and Huebler, 1989) and of Ott et al. (1990) are discussed, and the two schemes are compared for their ability to control shear flow, using fully developed and transitional solutions of the Ginzburg-Landau equation as models for such flows. It was found that the effectiveness of both methods in obtaining control of fully developed flows depended strongly on the 'distance' in state space between the uncontrolled flow and goal dynamics. There were conceptual difficulties in applying the Ott et al. method to transitional convectively unstable flows. On the other hand, the Huebler method worked well, within certain limitations, although at a large cost in energy terms.
The maintenance of variance and attendant heat flux in linear, forced dissipative baroclinic shear flows subject to stochastic excitation is examined. The baroclinic problem, is intrinsically nonnormal and its stochastic dynamics is found to differ significantly from the more familiar stochastic dynamics of normal systems. When the shear is sufficiently great in comparison to dissipative effects, stochastic excitation supports highly enhanced variance levels in these nonnormal systems compared to variance levels supproted by the same forcing and dissipation in related normal systems. The eddy variance and associated heat flux are found to arise in response to transient amplification of a subset of forcing functions that obtain energy from the mean flow and project this energy on a distinct subset of response functions (EOFs) that are in turn distinct from the set of normal modes of the system. A method for obtaining the dominant forcing and response functions as well as the distribution of heat flux for a given flow is described.
Electrical transport in semiconducting and metallic particle suspensions is an enabling feature of emerging grid-scale battery technologies. Although the physics of the transport process plays a key role in these technologies, no universal framework has yet emerged. Here, we examine the important contribution of shear flow to the electrical transport of non-Brownian suspensions. We find that these suspensions exhibit a strong dependence of the transport rate on the particle volume fraction and applied shear rate, which enables the conductivity to be dynamically changed by over 10 7 decades based on the applied shear rate. We combine experiments and simulations to conclude that the transport process relies on a combination of charge and particle diffusion with a rate that can be predicted using a quantitative physical model that incorporates the self-diffusion of the particles.
The laser-Doppler ambiguities predicted by George and Lumley (1973) have been verified experimentally for turbulent pipe flows. Experiments were performed at Reynolds numbers from 5000 to 15,000 at the center line and near the wall. Ambiguity levels were measured from power spectral densities of FM demodulated laser signals and were compared with calculations based on the theory. The turbulent spectra for these water flows after accounting for the ambiguity were equivalent to hot-film measurements at similar Reynolds numbers. The feasibility of laser-Doppler measurements very close to the wall in shear flows is demonstrated.