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.
The convective stability of a viscous liquid in spherical layers is investigated taking into consideration rotation, the latitudinal temperature gradient, and shear flow. The results of calculating nonlinear convective motion in spherical layers are examined. A discussion is given of the applicability of the results obtained to studying convection in astrophysical objects.
The mean-flow structure of supersonic, turbulent, adiabatic-wall flow in a square duct is investigated experimentally over a development length x/D = 0-50 for a uniform flow, Mach 3.9 condition at the duct inlet. The results show that a secondary flow cell structure develops which is similar to that for the incompressible case. Development of the primary flow is influenced by the combined effects of the secondary flow and the streamwise adverse pressure gradient. Total pressure, axial mean velocity, and Mach number profiles are presented which show that the outer flow is sensitive primarily to the streamwise pressure gradient, while flow in the near-wall region is dominated by the secondary flow. Axial mean-velocity profiles plotted in terms of van Driest-scaled variables show that a well-defined log-law region exists in the near-wall layer. This region exists in the presence of a secondary flow which continuously modifies spanwise wall shear stress behavior along the length of the duct.
In the absence of identifiable physiographic features directly associated with plate tectonics, alternate mechanisms are sought for the intense tectonic deformation observed in radar images of Venus. One possible mechanism is direct coupling into an elastic lithosphere of the stresses associated with convective flow in the interior. Spectral Green's function solutions have been obtained for stresses in an elastic lithosphere overlying a Newtonian interior with an exponential depth dependence of viscosity, and a specified surface-density distribution driving the flow. At long wavelengths and for a rigid elastic/fluid boundary condition, horizontal normal stresses in the elastic lid are controlled by the vertical shear stress gradient and are directly proportional to the depth of the density disturbance in the underlying fluid. The depth and strength of density anomalies in the Venusian interior inferred by analyses of long wavelength gravity data suggest that stresses in excess of 100 MPa would be generated in a 10 km thick elastic lid unless a low viscosity channel occurring beneath the lid or a positive viscosity gradient uncouples the flow stresses. The great apparent depth of compensation of topographic features argues against this, however, thus supporting the importance of the coupling mechanism. If there is no elastic lid, stresses will also be very high near the surface, providing also that the viscosity gradient is negative.
We investigate the origin of a significant reduction in computed numerical “noise” in the spatial variation of computed wall heat flux observed when using a node-centered weighted least squares(WLSQ) gradient approach when solving the Reynolds-averaged Navier-Stokes (RANS) equations for hypersonic boundary layer flows on highly stretched tetrahedral grids. This investigation is conducted using a 2nd-order, cell-centered, finite-volume, discretization which we employ to contrast results obtained using two WLSQ gradient approaches: a conventional 3-D, cell-centered, node-neighbor, WLSQ gradient (NN-CCG) approach and a 3-D, face-averaged, node-centered gradient (F-ANG)approach. The F-ANG approach, is shown to produce significantly less “noise” in the spatial variation of wall heat flux, than the NN-CCG approach and we demonstrate the cause of this “noise” via a series of numerical experiments where we compute hypersonic laminar and turbulent boundary layer flows.We further investigate the level of grid resolution required for a tetrahedral grid to approach the fidelity of a grid resolved hexahedral grid solution obtained using an equivalent structured grid method.
The properties of conical flow without axial symmetry are analyzed. The flow around cones of circular cross section at small angles of attack is determined by correctly considering the effect of the entropy gradients in the flow.
In this work, the effects of the particle volume fraction gradient on fluid-particle interactions are studied. The phase interaction force is decomposed into three terms. For the first term, namely the symmetrized force density, we present theoretical reasoning and numerical evidence to assume that it is independent of the particle volume fraction gradient. The second term is the particle volume fraction gradient times a newly introduced diffusion stress. The third term is the divergence of the particle-fluid-particle (PFP) stress. If this assumption of independence of the particle volume fraction gradient for the first term can be verified, to the first order of the ratio of the mean distance between particles to the macroscopic lengthscale, all three terms can be studied and modeled in flows with uniform particle distributions. Models thus obtained are applicable to statistically inhomogeneous flows, with the second and third terms accounting for statistical inhomogeneity. To verify this assumption, numerical simulations of flows passing fixed arrays of particles are performed. Both uniform and nonuniform particle volume fractions are studied and compared for disperse multiphase flows with the particle Reynolds numbers ranging from 1 to 100, and particle volume fraction ranging from 1% to 26% in statistically steady states. It is found that the symmetrized force (first) term can be well approximated by the drag force obtained from studies of uniform flows. The diffusion stress is positive along the flow direction and negative in the directions perpendicular to the flow. In the case of moving particles, this stress could potentially cause particle aggregation in the flow direction and dispersion in the directions perpendicular to the flow. Finally, the diffusion stress is only important when there is a volume fraction gradient, while the PFP stress can be important in inhomogeneous flows with either nonuniform particle concentrations or nonuniform average relative velocities between the phases.
A 2-D nodal weighted least-squares gradient method and a related face-averaged nodal gradient approach that were developed for use with triangular grids are extended to 3-D for use with tetrahedral grids. In addition, a method, developed in 2-D, to stabilize the iterative convergence of these methods on quadrilateral cells is described and extended to 3-D and remedies are investigated to determine the nodal gradient averaging approach most suitable for use with grids made up of hexahedral, prismatic, pyramidal and tetrahedral cells. Moreover, due to an interest in hypersonic flow, a robust multidimensional gradient limiter procedure that is consistent with the stencil used to construct the nodal gradients is described. Finally, we demonstrate that the resulting 3-D methods are sufficiently robust for use in scramjet computations through the solution of three canonical turbulent hypersonic flow problems as well as a physically realistic 3-D scramjet inlet geometry.
Computations are presented describing the mixing and combustion of swirling jets in a coaxial stream. It is demonstrated that the boundary layer equations represent the flow reasonably well until reversed flow is imminent. For the range of parameters investigated indications are that the edge velocity has little effect on the behavior of the flow. Furthermore, confining the flow with a constant pressure wall, or impressing a favorable pressure gradient on the coaxial flow, acts to reduce the severity of the centerline adverse pressure gradient created by the swirl decay. A simple scalar eddy viscosity model, including a potential core formulation, is shown to described the behavior of weak swirling flow in the far region but is only in fair agreement with observations in the near region. The effects of swirl on a burning hydrocarbon jet exhausting into a cold coaxial stream are shown to be intensified by the reduction of the density due to combustion. The enhanced mixing properties of high swirl flow produce rapid diffusion of the burning gases into the cold edge flow causing early cessation of the NO producing reactions. Computations show that doubling the initial jet swirl could reduce the NO production by 25 percent.
Observations suggesting that the mean solar wind azimuthal field strength B(theta) near the ecliptic plane falls off more rapidly with heliocentric distance than would be expected in a classic Parker expansion is reexamined from a theoretical perspective using a three-dimensional MHD nonlinear numerical model for steady, corotating flow. For realistic solar wind parameters, it is found that a purely axisymmetric expansion can produce sizable magnetic flux deficits only when there are substantial meridional gradients in mean flow conditions localized about the ecliptic plane near the sun. Calculations on three-dimensional cororating flows are presented which demonstrate that latitudinal transport of magnetic flux by stream interactions may be an important consideration in generating the deficits in mean B(theta).
An impedance eduction technique that is based upon a numerical solution to the linearized Euler equations coupled with the minimization of an objective function is presented. The linearized Euler equations are solved numerically using a cubic finite element method and the minimum of the objective function is obtained via a gradient based optimizer. The mean flow is allowed to have a gradient in two coordinate directions and the impedance of the liner may vary arbitrarily over the liner surface. A new regularized wall impedance boundary condition presented by Rienstra which contains the effects of the mean boundary layer thickness is implemented. Impedances are educed for a perforate over honeycomb liner using test data acquired in the Langley Grazing Flow Impedance Tube. The impedance of the liner is educed both with an upstream source and a downstream source. The method is validated by comparing the educed impedance to a benchmark method. The primary conclusion of the study is that the impedances educed for upstream and downstream sources are not well matched and that the primary effect of the boundary thickness is to reduce the resistance of the liner.
Researchers at the NASA Langley Research Center are applying laser velocimetry (LV) techniques to characterize the fluid dynamics of non-isothermal flows inside fused silica chambers designed for chemical vapor deposition (CVD). Experimental issues involved in the application of LV techniques to this task include thermophoretic effects on the LV seed particles, seeding the hazardous gases, index of refraction gradients in the flow field and surrounding media, optical access, relatively low flow velocities, and analysis and presentation of sparse data. An overview of the practical difficulties these issues represent to the use of laser velocimetry instrumentation for CVD applications is given. A fundamental limitation on the application of LV techniques in non-isothermal systems is addressed which involves a measurement bias due to the presence of thermal gradients. This bias results from thermophoretic effects which cause seed particle trajectories to deviate from gas streamlines. Data from a research CVD reactor are presented which indicate that current models for the interaction of forces such as Stokes drag, inertia, gravity, and thermophoresis are not adequate to predict thermophoretic effects on particle-based velocimetry measurements in arbitrary flow configurations.
A dual laser optical system and method is disclosed for visualization of phenomena in transport substances which induce refractive index gradients such as fluid flow and pressure and temperature gradients in fluids and gases. Two images representing mutually perpendicular components of refractive index gradients may be viewed simultaneously on screen. Two lasers having wave lengths in the visible range but separated by about 1000 angstroms are utilized to provide beams which are collimated into a beam containing components of the different wave lengths. The collimated beam is passed through a test volume of the transparent substance. The collimated beam is then separated into components of the different wave lengths and focused onto a pair of knife edges arranged mutually perpendicular to produce and project images onto the screen.
Experimental evidence for the existence of shear induced migration processes is reviewed and the mechanism by Leighton and Acrivos (1987b) is described in detail. The proposed mechanism is shown to lead to the existence of an additional shear induced migration in the presence of gradients in shear stress such as would be found in Poiseuille flow, and which may be used to predict the amplitude of the observed short-term viscosity increase. The concentration and velocity profiles which result from such a migration are discussed in detail and are compared to the experimental observations of Karnis, Goldsmith and Mason (1966).
Previous studies on cross-shelf exchange, predominantly adopted an Eulerian perspective, struggled to identify water sources and pathways. Using a high-resolution regional ocean modeling system (ROMS) and Lagrangian particle tracking, this study systematically investigates the seasonal variation and dynamics of surface cross-shelf exchange in the northern South China Sea (NSCS) from a Lagrangian perspective. Based on daily released 30-day drifter trajectories we identify the key pathways, source regions for surface cross-shelf exchange, revealing pronounced seasonal variability. Results show the surface cross-shelf exchange generally following monsoon-driven Ekman transport. However, local dynamics, especially topographic modulation, can outweigh the expected Ekman-driven transport, producing surface exchange patterns opposite to that predicted from the prevailing winds. Topographic effects vary across different scales. In the coastal waters of western Guangdong during winter, despite downwelling-favorable winds, the modulation of alongshore currents by island topography induces an eastward pressure gradient. This gradient, through geostrophic balance, drives offshore flows opposite to wind-driven Ekman onshore transport. Furthermore, the eastern widened shelf exhibits a distinct seasonal variation of cross-shelf exchange, with strong offshore transport (opposite to the direction of Ekman transport) in winter and exceedingly weak exchange in summer. Analysis of the underlying mechanisms reveals that this winter offshore transport is primarily attributed to geostrophic flows driven by surface pressure gradient that is jointly modulated by the Kuroshio intrusion and local widened shelf topography, and enhanced by cumulative submesoscale processes. In summer, a persistent and strong along-isobath jet acts as a dynamic barrier, effectively suppressing the exchange. These findings highlight the important role of topography in regulating surface material transport, and have important implications for predicting the advection and dispersion of plankton or oil spills over the continental shelf influenced by monsoon.
Experimental data were interpreted using two supersonic combustion computer programs. The P1 program is based on a conventional boundary layer treatment of the mixing of concentric gas streams and complete combustion chemistry. The H1 program is based on a modified boundary layer approach which accounts for radial pressure gradients in the flow and also incorporates a finite rate chemistry calculation. The objective of the investigation was to compare the experimental data with theoretical predictions of the two programs with special emphasis on the prediction of radial pressure gradients by the H1 program. A test of the H1 program was also desired through comparison with the experimental data and with the P1 program.
Two papers are included in this progress report. In the first, the compressible Navier-Stokes equations have been used to compute leading edge receptivity of boundary layers over parabolic cylinders. Natural receptivity at the leading edge was simulated and Tollmien-Schlichting waves were observed to develop in response to an acoustic disturbance, applied through the farfield boundary conditions. To facilitate comparison with previous work, all computations were carried out at a free stream Mach number of 0.3. The spatial and temporal behavior of the flowfields are calculated through the use of finite volume algorithms and Runge-Kutta integration. The results are dominated by strong decay of the Tollmien-Schlichting wave due to the presence of the mean flow favorable pressure gradient. The effects of numerical dissipation, forcing frequency, and nose radius are studied. The Strouhal number is shown to have the greatest effect on the unsteady results. In the second paper, a transition model for low-speed flows, previously developed by Young et al., which incorporates first-mode (Tollmien-Schlichting) disturbance information from linear stability theory has been extended to high-speed flow by incorporating the effects of second mode disturbances. The transition model is incorporated into a Reynolds-averaged Navier-Stokes solver with a one-equation turbulence model. Results using a variable turbulent Prandtl number approach demonstrate that the current model accurately reproduces available experimental data for first and second-mode dominated transitional flows. The performance of the present model shows significant improvement over previous transition modeling attempts.
The study of edge flames has received increased attention in recent years. This work reports the results of a recent study into two-dimensional, planar, propagating edge flames that are remote from solid surfaces (called here, free-layer flames, as opposed to layered flames along floors or ceilings). They represent an ideal case of a flame propagating down a flammable plume, or through a flammable layer in microgravity. The results were generated using a new apparatus in which a thin stream of gaseous fuel is injected into a low-speed laminar wind tunnel thereby forming a flammable layer along the centerline. An airfoil-shaped fuel dispenser downstream of the duct inlet issues ethane from a slot in the trailing edge. The air and ethane mix due to mass diffusion while flowing up towards the duct exit, forming a flammable layer with a steep lateral fuel concentration gradient and smaller axial fuel concentration gradient. We characterized the flow and fuel concentration fields in the duct using hot wire anemometer scans, flow visualization using smoke traces, and non-reacting, numerical modeling using COSMOSFloWorks. In the experiment, a hot wire near the exit ignites the ethane air layer, with the flame propagating downwards towards the fuel source. Reported here are tests with the air inlet velocity of 25 cm/s and ethane flows of 967-1299 sccm, which gave conditions ranging from lean to rich along the centerline. In these conditions the flame spreads at a constant rate faster than the laminar burning rate for a premixed ethane air mixture. The flame spread rate increases with increasing transverse fuel gradient (obtained by increasing the fuel flow rate), but appears to reach a maximum. The flow field shows little effect due to the flame approach near the igniter, but shows significant effect, including flow reversal, well ahead of the flame as it approaches the airfoil fuel source.
The flow fields near a 45 degree swept-wing-fuselage combination at moderate angles of sideslip (plus-or-minus 8 degrees), as determined experimentally at low speed, are presented as variations with chordwise distance for various spanwise and vertical locations and angles of attack. The results indicated that for positions close to the fuselage (on and near the plane of symmetry) changes in the angle of sideslip caused large changes in the flow-field characteristics and particularly in the local angles of sideslip, which in some cases were nearly double the static angle of sideslip. In general, the effects of changing the angle of sideslip on the flow-field characteristics for all of the outboard underwing locations were qualitatively similar, although conditions at the more inboard and outboard locations were somewhat more severe for lifting conditions than at the one-half semispan location. The chordwise gradients in the flow parameters for the underwing locations were more severe than for the fuselage locations although the effect of changing the angle of sideslip was less severe, in that the incremental changes in the local angles of sideslip were approximately equal to the static angle of sideslip. Flow conditions near the wing tip were found to be critically dependent on vertical location, with the largest sideslip-induced variations occurring at the nearest vertical locations. The results also indicated that for the outboard underwing locations the wing was the predominant factor in disturbing the field of flow for the conditions investigated.