Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boussinesq flow”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Differential rotation and turbulent convection: A new Reynolds stress model and comparison with solar data

In most hydrodynamic cases, the existence of a turbulent flow superimposed on a mean flow is caused by a shear instability in the latter. Boussinesq suggested the first model for the turbulent Reynolds stresses bar-(u(sub i)u(sub j)) in which the mean shear S(sub ij) is the cause (or source) of turbulence represented by the stress bar-(u(sub i)u(sub j)). In the case of solar differential rotation, exactly the reverse physical process occurs: turbulence (which must pre-exist) generates a mean flow which manifests itself in the form of differential rotation. Thus, the Boussinesq model is wholly inadequate because in the solar case, cause and effect are reversed. Since the Boussinesq model is inadequate, one needs an alternative model for the Reynolds stresses. We present a new dynamical model for the Reynolds stresses, convective fluxes, turbulent kinetic energy, and temperature fluctuations. The complete model requires the solution of 11 differential equations. We then introduce a set of simplifying assumptions which reduce the full dynamical model to a set of algebraic Reynolds stress models. We explicitly solve one of these models that entails only one differential equation. The overall agreement with the data is obtained with a model that is neither phenomenological nor one that requires a full numerical simulation, since it is algebraic in nature. The new model can play an important role in understanding the complex physics underlying the interplay between solar differential rotation and convection, as many physical processes can naturally be incorporated into the model.

Canuto, V. M.↗
NASA NTRS1992

A unified approach for numerical simulation of viscous compressible and incompressible flows over adiabatic and isothermal walls

A new formulation (including the choice of variables, their non-dimensionalization, and the form of the artificial viscosity) is proposed for the numerical solution of the full Navier-Stokes equations for compressible and incompressible flows with heat transfer. With the present approach, the same code can be used for constant as well as variable density flows. The changes of the density due to pressure and temperature variations are identified and it is shown that the low Mach number approximation is a special case. At zero Mach number, the density changes due to the temperature variation are accounted for, mainly through a body force term in the momentum equation. It is also shown that the Boussinesq approximation of the buoyancy effects in an incompressible flow is a special case. To demonstrate the new capability, three examples are tested. Flows in driven cavities with adiabatic and isothermal walls are simulated with the same code as well as incompressible and supersonic flows over a wall with and without a groove. Finally, viscous flow simulations of an oblique shock reflection from a flat plate are shown to be in good agreement with the solutions available in literature.

Hafez, M.↗
NASA NTRS1992

Surface tension and buoyancy-driven flow in a non-isothermal liquid bridge

The Navier-Stokes-Boussinesq equations governing the transport of momentum, mass and heat in a nonisothermal liquid bridge with a temperature-dependent surface tension are solved using a vorticity-stream-function formulation together with a nonorthogonal coordinate transformation. The equations are discretized using a pseudo-unsteady semi-implicit finite difference scheme and are solved by the ADI method. A Picard-type iteration is adopted which consists of inner and outer iterative processes. The outer iteration is used to update the shape of the free surface. Two schemes have been used for the outer iteration; both use the force balance normal to the free surface as the distinguished boundary condition. The first scheme involves successive approximation by the direct solution of the distinguished boundary condition. The second scheme uses the artificial force imbalance between the fluid pressure, viscous and capillary forces at the free surface which arises when the boundary condition for force balance normal to the surface is not satisfied. This artificial imbalance is then used to change the surface shape until the distinguished boundary condition is satisfied. These schemes have been used to examine a variety of model liquid bridge situations including purely thermocapillary-driven flow situations and mixed thermocapillary- and bouyancy-driven flow.

Zhang, Yiqiang↗
NASA NTRS1983

Multiple buoyancy driven flows in a vertical cylinder heated from below

The structure of axisymmetric buoyancy-driven convection in a vertical cylinder heated from below is probed by finite element solution of the Boussinesq equations coupled with computed-implemented perturbation techniques for detecting and tracking multiple flows and for determining flow stability. Results are reported for fluids with Prandtl number of one and for cylinders with aspect ratio (Lambda) (defined as the height to radius of the cylinder) between 0.5 and 2.25. Extensive calculations of the neutral stability curve for the static solution and of the nonlinear motions along the bifurcating flow families show a continuous evolution of the primary cellular motion from a single toroidal cell to two and three cells nested radially in the cylinder, instead of the sharp transitions found for a cylinder with shear-free sidewalls. The smooth transitions in flow structure with Rayleigh number and lambda are explained by nonlinear connectivity between the first two bifurcating flow families formed either by a secondary bifurcation point for Lambda or = Lambda * approximately 0.80 or by a limit point for Lambda Lambda *. The transition between these two modes may be described by the theory of multiple limit point bifurcation.

Yamaguchi, Y.↗
NASA NTRS1994

Incompressible Navier-Stokes Computations with Heat Transfer

The existing pseudocompressibility method for the system of incompressible Navier-Stokes equations is extended to heat transfer problems by including the energy equation. The solution method is based on the pseudo compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Current computations use one-equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard k-epsilon model equations. Both forced and natural convection problems are examined. Numerical results from turbulent reattaching flow behind a backward-facing step will be compared against experimental measurements for the forced convection case. The validity of Boussinesq approximation to simplify the buoyancy force term will be investigated. The natural convective flow structure generated by heat transfer in a vertical rectangular cavity will be studied. The numerical results will be compared by experimental measurements by Morrison and Tran.

Kiris, Cetin↗
NASA NTRS1982

Bifurcation and stability of low-order steady flows in horizontally and vertically forced convection

A nonlinear spectral model of two-dimensional, shallow Boussinesq convection which responds to heating in both the horizontal and vertical directions is examined. The governing partial differential system is converted to an infinite set of ordinary differential equations and truncated to a small set to permit detailed study of the number and types of transitions from one flow configuration to another. The Hadley number and the Rayleigh number are defined as the horizontal and vertical thermal forcing mechanisms, respectively, for inclusion in the nonlinear spectral model, which is composed of three equations. The model is then used to describe steady states, linearly stable solutions, and balancing factors in unstable stratification. The number and the distribution of the steady states are found to be qualitatively independent of the aspect ratio and the Prandtl number.

Yost, D. A.↗
NASA NTRS1993

Banded surface flow maintained by convection in a model of the rapidly rotating giant planets

In three-dimensional numerical simulations of a rapidly rotating Boussinesq fluid shell, thermally driven convection in the form of columns parallel to the rotation axis generates an alternately directed mean zonal flow with a cylindrical structure. The mean structure at the outer spherical surface consists of a broad eastward flow at the equator and alternating bands of westward and eastward flows at higher latitudes in both hemispheres. The banded structure persists even though the underlying convective motions are time-dependent. These results, although still far from the actual motions seen on Jupiter and Saturn, provide support for theoretical suggestions that thermal convection can account for the remarkable banded flow structures on these planets.

Sun, Zi-Ping↗
NASA NTRS1991

Nonadiabatic and three-dimensional effects in compressible turbulent boundary layers

A defect stream function formulation for nonadiabatic flow with small crossflow is developed. The first-integral property of this formulation provides for two removal of the streamline curvature term in the governing equation so that the form of the reduced equation for small crossflow is the same as that for two dimensional flow. The combined law of the wall and wake is used in place of the no-slip boundary condition. The tangential velocity equation for law-of-the-wall flow is shown to be the same for three-dimensions as for two when the Boussinesq approximation applies, and a closed form solution for the crossflow angle in the inner region is obtained. Analytic solutions for nonadiabatic, compressible, equilibrium flow with a Clauser outer-region eddy-viscosity model are obtained, and excellent agreement with experimental skin friction and velocity profile data for nonadiabatic, compressible flat-plate flow is achieved. An analytic solution for a linear inner-region eddy-viscosity model is also obtained; the wake function part of this solution is found to be inconsistent with the empirically established law of the wake.

Barnwell, Richard W.↗
NASA NTRS2009

Enhancement of USM3D Unstructured Flow Solver for High-Speed High-Temperature Shear Flows

Large temperature and pressure fluctuations have a profound effect on turbulence development in transonic and supersonic jets. For high-speed, high-temperature jet flows, standard turbulence models lack the ability to predict the observed mixing rate of a shear layer. Several proposals to address this deficiency have been advanced in the literature to modify the turbulence transport equations in a variety of ways. In the present study, some of the most proven and simple modifications to two-equation turbulence models have been selected and implemented in NASA's USM3D tetrahedral Navier-Stokes flow solver. The modifications include the addition of compressibility correction and pressure dilatation terms in the turbulence transport equations for high-speed flows, and the addition of a simple modification to the Boussinesq's closure model coefficient for high-temperature jets. The efficacy of the extended models is demonstrated by comparison with experimental data for two supersonic axisymmetric jet test cases at design pressure ratio.

Pandya, Mohagna J.↗
NASA NTRS1989

The stability of a compressible stratified shear layer

The stability of a shear layer under the effect of gravity is investigated using the compressible magnetohydrodynamic (MHD) equations, including an effective gravity term to represent the curvature effects of the flow and magnetic field line geometry. A general eigenmode equation is derived for a two-dimensional MHD fluid, and an energy-principle analysis to explain the effect of compressibility on the critical Richardson number is presented. For the case of a hyperbolic tangent shear flow and exponential density profile, it was found that, in the Boussinesq approximation, the compressibility raises the critical Richardson number from 1/4 to as much as 1/2, with the exact value depending on the value of the magnetic field at infinity. Under approximation of a strong asymptotic magnetic field, without invoking the Boussinesq approximation, it is shown both analytically and numerically that the density gradient terms cause the shear instability to be dispersive. The long-wavelength stability boundary for the Richardson number J = 0 is characterized by a normalized phase velocity c =

Wang, Z.↗
NASA NTRS1992

Direct simulation of the stably stratified turbulent Ekman layer

The Navier-Stokes equations and the Boussinesq approximation were used to compute a 3D time-dependent turbulent flow in the stably stratified Ekman layer over a smooth surface. The simulation data are found to be in very good agreement with atmospheric measurements when nondimensionalized according to Nieuwstadt's local scaling scheme. Results suggest that, when Reynolds number effects are taken into account, the 'constant Froud number' stable layer model (Brost and Wyngaard, 1978) and the 'shearing length' stable layer model (Hunt, 1985) for the dissipitation rate of turbulent kinetic energy are both valid. It is concluded that there is good agreement between the direct numerical simulation results and large-eddy simulation results obtained by Mason and Derbyshire (1990).

Coleman, G. N.↗
NASA NTRS1992

The program FANS-3D (finite analytic numerical simulation 3-dimensional) and its applications

In this study, the program named FANS-3D (Finite Analytic Numerical Simulation-3 Dimensional) is presented. FANS-3D was designed to solve problems of incompressible fluid flow and combined modes of heat transfer. It solves problems with conduction and convection modes of heat transfer in laminar flow, with provisions for radiation and turbulent flows. It can solve singular or conjugate modes of heat transfer. It also solves problems in natural convection, using the Boussinesq approximation. FANS-3D was designed to solve heat transfer problems inside one, two and three dimensional geometries that can be represented by orthogonal planes in a Cartesian coordinate system. It can solve internal and external flows using appropriate boundary conditions such as symmetric, periodic and user specified.

Bravo, Ramiro H.↗
NASA NTRS1980

A simple approximate result for the maximum growth rate of baroclinic instabilities

The Charney problem for baroclinic instability involves the quasi-geostrophic instability of a zonal flow on a beta plane where the zonal flow is characterized by a constant vertical shear. The atmosphere is non-Boussinesq and continuous. The solution of this problem involves confluent hypergeometric functions, and the mathematical difficulty of the problem has precluded extracting simple results of generality. It is shown that there exists very simple, powerful approximate result for the growth rate of the most rapidly growing instability, viz., that this growth rate is linearly proportional to the surface meridional temperature gradient. The coefficient of proportionality is also easily determined. The result extends to substantially more general profiles than those in the Charney problem.

Lindzen, R. S.↗
NASA NTRS1979

A theoretical study of three-dimensional barotropic instability with applications to the upper stratosphere

A numerical method for tropospheric baroclinic instability is used to establish the structures, phase speeds and growth rates of the normal modes of a number of idealized barotropically unstable mean zonal wind fields in a quasi-geostrophic Boussinesq framework. The objective was to gain insight into the effects of vertical mean flow variation on barotropically unstable stratosphere. Both horizontally symmetric and asymmetric flows are considered, but in all cases the flows are vertically symmetric. Major findings show that (1) the wave growth rates associated with the strongest expected horizontal shears decrease with increasing vertical shear and are at least 28% smaller for the weakest expected vertical shear than with no shear; (2) differences in the form of either the vertical or the horizontal shear flow can affect this growth rate reduction by about 25%; and (3) the vertical scale of the barotropic unstable waves is of the order of the geometric mean of the vertical mean flow scale and the vertical penetrat

Pfister, L.↗
NASA NTRS1997

A Numerical Method for Incompressible Flow with Heat Transfer

A numerical method for the convective heat transfer problem is developed for low speed flow at mild temperatures. A simplified energy equation is added to the incompressible Navier-Stokes formulation by using Boussinesq approximation to account for the buoyancy force. A pseudocompressibility method is used to solve the resulting set of equations for steady-state solutions in conjunction with an approximate factorization scheme. A Neumann-type pressure boundary condition is devised to account for the interaction between pressure and temperature terms, especially near a heated or cooled solid boundary. It is shown that the present method is capable of predicting the temperature field in an incompressible flow.

Sa, Jong-Youb↗
NASA NTRS1996

Non-linear k-epsilon-v(sup 2)(bar) modeling with application to high-lift

The k-epsilon-v(sup 2)(bar) model has been investigated to quantify its predictive performance on two high-lift configurations: 2D flow over a single-element aerofoil, involving closed-type separation; 3D flow over a prolate spheroid, involving open-type separation. A 'code-friendly' modification has been proposed which enhances the numerical stability, in particular, for explicit and uncoupled flow solvers. As a result of introducing Reynolds-number dependence into a coefficient of the s-equation, the skin-friction distribution for the by-pass transitional flow over a flat plate is better predicted. In order to improve deficiencies arising from the Boussinesq approximation, a nonlinear stress-strain constitutive relation was adopted, in which the only one free constant is calibrated on the basis of DNS data, and the Reynolds-stress anisotropy near the wall is fairly well represented.

Lien, F. S.↗
NASA NTRS2000

Measurement of the Surface Dilatational Viscosity of an Insoluble Surfactant Monolayer at the Air/Water Interface Using a Pendant Drop Apparatus

When a fluid interface with surfactants is at rest, the interfacial stress is isotropic (as given by the equilibrium interfacial tension), and is described by the equation of state which relates the surface tension to the surfactant surface concentration. When surfactants are subjected to shear and dilatational flows, flow induced interaction of the surfactants; can create interfacial stresses apart from the equilibrium surface tension. The simplest relationship between surface strain rate and surface stress is the Boussinesq-Scriven constitutive equation completely characterized by three coefficients: equilibrium interfacial tension, surface shear viscosity, and surface dilatational viscosity Equilibrium interfacial tension and surface shear viscosity measurements are very well established. On the other hand, surface dilatational viscosity measurements are difficult because a flow which change the surface area also changes the surfactant surface concentration creating changes in the equilibrium interfacial tension that must be also taken into account. Surface dilatational viscosity measurements of existing techniques differ by five orders of magnitude and use spatially damped surface waves and rapidly expanding bubbles. In this presentation we introduce a new technique for measuring the surface dilatational viscosity by contracting an aqueous pendant drop attached to a needle tip and having and insoluble surfactant monolayer at the air-water interface. The isotropic total tension on the surface consists of the equilibrium surface tension and the tension due to the dilation. Compression rates are undertaken slow enough so that bulk hydrodynamic stresses are small compared to the surface tension force. Under these conditions we show that the total tension is uniform along the surface and that the Young-Laplace equation governs the drop shape with the equilibrium surface tension replaced by the constant surface isotropic stress. We illustrate this technique using DPPC as the insoluble surfacant monolayer and measured for it a surface dilatational viscosity in the LE phase that is 20 surface poise.

Lorenzo, Jose↗
NASA NTRS1973

General atmospheric circulation driven by polar and diurnal surface temperature variations.

Described is a global circulation model for the Venus atmosphere that includes the effects of both polar cooling and diurnal temperature variation. It is based on a linearized Boussinesq approximation and boundary conditions derived from theoretical and empirical considerations. The time-dependent, three-dimensional flow field is deduced without any a priori assumptions about its configuration. Results show that the mean atmospheric motions are essentially zonal in a narrow belt near the equator and change to become meridional over most of the globe. The circulation pattern is not symmetrical and rotates about the polar axis of the planet with the period of the solar day.

Bohachevsky, I. O.↗