Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boussinesq approximation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Suppression of Buoyancy in Gaseous Media at High Temperatures

Consider a rectangular box filled with a fluid having a heated bottom and a cold top surface, and insulated side-walls (Benard problem). As the temperature difference between the horizontal top and bottom surfaces increases, a critical condition, defined quantitatively by the Rayleigh number, is reached beyond which density stratification can no longer be sustained by conduction and the fluid disrupts from its stable, quiescent state into an unstable, convective mode in which lighter and heavier gas mix. This paper suggests that such a statement is not necessarily true for gaseous media under normalized temperature differences that are much larger than justifiable for the Boussinesq approximation! In fact, there may be situations where a system cannot ever be made unstable with respect to the onset on buoyant convection no matter how large the temperature (density) difference becomes at a given pressure even under normal gravity! This unexpected behavior is primarily attributed to highly temperature-sensitive kinematic viscosity which counteracts the tendency toward instability and dampens convection by making the gas more viscous at higher temperatures. This compensation of the buoyant force by the viscous force exhibits itself by the formation of a peak hot-surface temperature beyond which a system will tend to be more stable as the hot-surface temperature increases.

Gokoglu, Suleyman A.↗
NASA NTRS2005

Partially-Averaged Navier Stokes Model for Turbulence: Implementation and Validation

Partially-averaged Navier Stokes (PANS) is a suite of turbulence closure models of various modeled-to-resolved scale ratios ranging from Reynolds-averaged Navier Stokes (RANS) to Navier-Stokes (direct numerical simulations). The objective of PANS, like hybrid models, is to resolve large scale structures at reasonable computational expense. The modeled-to-resolved scale ratio or the level of physical resolution in PANS is quantified by two parameters: the unresolved-to-total ratios of kinetic energy (f(sub k)) and dissipation (f(sub epsilon)). The unresolved-scale stress is modeled with the Boussinesq approximation and modeled transport equations are solved for the unresolved kinetic energy and dissipation. In this paper, we first present a brief discussion of the PANS philosophy followed by a description of the implementation procedure and finally perform preliminary evaluation in benchmark problems.

Girimaji, Sharath S.↗
NASA NTRS2012

Magnetic Helicity and Planetary Dynamos

A model planetary dynamo based on the Boussinesq approximation along with homogeneous boundary conditions is considered. A statistical theory describing a large-scale MHD dynamo is found, in which magnetic helicity is the critical parameter

Shebalin, John V.↗
NASA NTRS1996

Inertial Range Dynamics in Boussinesq Turbulence

L'vov and Falkovich have shown that the dimensionally possible inertial range scaling laws for Boussinesq turbulence, Kolmogorov and Bolgiano scaling, describe steady states with constant flux of kinetic energy and of entropy respectively. These scaling laws are treated as similarity solutions of the direct interaction approximation for Boussinesq turbulence. The Kolmogorov scaling solution corresponds to a weak perturbation by gravity of a state in which the temperature is a passive scalar but in which a source of temperature fluctuations exists. Using standard inertial range balances, the renormalized viscosity and conductivity, turbulent Prandtl number, and spectral scaling law constants are computed for Bolgiano scaling.

Rubinstein, Robert↗
NASA NTRS1985

A numerical study of global axisymmetric states and their stability to quasi-geostrophic disturbances

A steady, axisymmetric model of the general circulation is developed as a basis for climate stability studies. The model includes the effects of heating, rotation, and internal friction, but neglects topography. It is assumed that the axisymmetric flow may be modeled by making the Boussinesq and deep convection approximations. The hydrostatic assumption is not made, thus permitting the advective terms to be included in the vertical equation of motion. The initial set of five primitive equations is reduced to three equations in terms of the zonal velocity, meridional streamfunction, and the potential temperature perturbation. The application of the Boussinesq, deep convection, and quasi-geostrophic assumptions limits the ranges of the heating and rotation rates. For values not too far from typical atmospheric values, the model produces a stability boundary separating Hadley from Rossby flow. The boundary is characterized by a particular value of vertical wind shear, which suggests that baroclinic instability is the primary mechanism for the loss of stability. The initial growth rates are largest for longitudinal waves 4-7, also in agreement with studies of baroclinic instability.

Henderson, H. W.↗
NASA NTRS2000

Modeling of a Non-Dilute Alloy Solidification Under Terrestrial and Microgravity Conditions

A general finite element model of non-dilute alloy solidification was used to simulate growth of HgCdTe from the melt. 2-D and 3-D models with and without Oberbeck-Boussinesq (O-B) approximation were used. Verification of the 3-D model was based on experimental data. The location and the shape of the interface were determined using melting temperatures obtained from the HgCdTe liquidus curve. Double-diffusive convection in the melt is caused by radial temperature gradients and by rejection at the solid-melt interface of the heavier HgTe-rich solute. Because of these complicating factors, dimensional rather than non-dimensional modeling was performed. Models based on complete set of Navier-Stokes equations for compressible media are superior to the O-B approximation in microgravity, because they provide more accurate estimation of convection.

Bune, Andris V.↗
NASA NTRS1980

Nonlinear axially symmetric circulations in a nearly inviscid atmosphere

The structure of certain axially symmetric circulations in a stably stratified, differentially heated, rotating Boussinesq fluid on a sphere is analyzed. A simple approximate theory (similar to that introduced by Schneider (1977)) is developed for the case in which the fluid is sufficiently inviscid that the poleward flow in the Hadley cell is nearly angular momentum conserving. The theory predicts the width of the Hadley cell, the total poleward heat flux, the latitude of the upper level jet in the zonal wind, and the distribution of surface easterlies and westerlies. Fundamental differences between such nearly inviscid circulations and the more commonly studied viscous axisymmetric flows are emphasized. The theory is checked against numerical solutions to the model equations.

Held, I. M.↗
NASA NTRS1992

Convection treatment for high Rayleigh number, laminar, natural convection calculation

The problem of natural convection in a square cavity of a Boussinesq incompressible fluid is computed using a consistent central-difference approximation, a consistent second-order upwind scheme, and the hybrid scheme for the energy equation. It was found that, for this problem, the second-order upwind scheme yields satisfactory solutions and that an apparently low-order combination of schemes, with the second-order upwind scheme for the momentum equations and the hybrid scheme for the energy equation, can produce satisfactory predictions for this problem.

Shyy, Wei↗
NASA NTRS1991

Numerical simulations of turbulent thermal convection with differential rotation

The solar nebula, from which the planets in our solar system formed, featured a disk of gas and dust grains in rapid, differential rotation, and at some stage was likely to have been unstable to thermal convection. This situation is suspected by many to lead to significant turbulent Reynolds stress production and angular momentum transport in such systems, and estimates of transport rates have been attempted from unsubstantiated phenomenological models. In order to determine the circumstances and physical conditions under which our own planetary system formed and to explain recent observations of young stellar systems, it is necessary to develop realistic models of heat and angular momentum transport for such flows. Developing an understanding of complicated flows featuring thermal convection, rotation, and shear is also of wide interest in stellar astrophysics and in planetary and terrestrial atmospheric studies. The ultimate objective is to develop workable models based on the numerical simulations for constructing global solar nebula models; viz., relatively simple prescriptions for heat and angular momentum fluxes from given system parameters (e.g., ratios of rotation, shear, and convective lapse rates) are characterized, quantified, and developed. Toward this end, our program has been to attempt to understand the behavior of the direct numerical simulations of Boussinesq convection, which, despite the complexity of the results, is still an overly simplified approximation to the real system and should be more amenable to analysis. These results are also intended to be tested against turbulence models, especially those designed for atmospheric boundary layers, and may provide a basis for subgrid-scale models. In order to make the numerical simulations more realistic with regard to the solar nebula problem, a fully compressible code that will allow incorporation of large density stratifications and realistic thermodynamic and radiative properties is developed. In order to explore the properties of these flows at the very high values of Re found in natural systems and the very low values of Pr found in most astrophysical contexts, we will need to employ large-eddy simulations for which we want to determine the most appropriate subgrid-scale model to incorporate.

Cabot, William 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 NTRS1975

Modal equations for cellular convection

We expand the fluctuating flow variables of Boussinesq convection in the planform functions of linear theory. Our proposal is to consider a drastic truncation of this expansion as a possible useful approximation scheme for studying cellular convection. With just one term included, we obtain a fairly simple set of equations which reproduces some of the qualitative properties of cellular convection and whose steady-state form has already been derived by Roberts (1966). This set of 'modal equations' is analyzed at slightly supercritical and at very high Rayleigh numbers. In the latter regime the Nusselt number varies with Rayleigh number just as in the mean-field approximation with one horizontal scale when the boundaries are rigid. However, the Nusselt number now depends also on the Prandtl number in a way that seems compatible with experiment. The chief difficulty with the approach is the absence of a deductive scheme for deciding which planforms should be retained in the truncated expansion.

Gough, D. O.↗
NASA NTRS1989

Blow-up of unsteady two-dimensional Euler and Navier-Stokes solutions having stagnation-point form

The time-dependent form of the classic, two-dimensional stagnation-point solution of the Navier-Stokes equations is considered. If the viscosity is zero, a class of solutions of the initial-value problem can be found in closed form using Lagrangian coordinates. These solutions exhibit singular behavior in finite time, because of the infinite domain and unbounded initial vorticity. Thus, the blow-up found by Stuart in three dimensions using the stagnation-point form, also occurs in two. The singularity vanishes under a discrete, finite-dimensional 'point vortex' approximation, but is recovered as the number of vortices tends to infinity. We find that a small positive viscosity does not arrest the breakdown, but does strongly alter its form. Similar results are summarized for certain Boussinesq stratified flows.

Childress, S.↗
NASA NTRS1991

GEOSIM: A numerical model for geophysical fluid flow simulation

A numerical model which simulates geophysical fluid flow in a wide range of problems is described in detail, and comparisons of some of the model's results are made with previous experimental and numerical studies. The model is based upon the Boussinesq Navier-Stokes equations in spherical coordinates, which can be reduced to a cylindrical system when latitudinal walls are used near the pole and the ratio of latitudinal length to the radius of the sphere is small. The equations are approximated by finite differences in the meridional plane and spectral decomposition in the azimuthal direction. The user can specify a variety of boundary and initial conditions, and there are five different spectral truncation options. The results of five validation cases are presented: (1) the transition between axisymmetric flow and baroclinic wave flow in the side heated annulus; (2) the steady baroclinic wave of the side heated annulus; (3) the wave amplitude vacillation of the side heated annulus; (4) transition to baroclinic wave flow in a bottom heated annulus; and (5) the Spacelab Geophysical Fluid Flow Cell (spherical) experiment.

Butler, Karen A.↗
NASA NTRS1994

Renormalization Group Theory of Bolgiano Scaling in Boussinesq Turbulence

Bolgiano scaling in Boussinesq turbulence is analyzed using the Yakhot-Orszag renormalization group. For this purpose, an isotropic model is introduced. Scaling exponents are calculated by forcing the temperature equation so that the temperature variance flux is constant in the inertial range. Universal amplitudes associated with the scaling laws are computed by expanding about a logarithmic theory. Connections between this formalism and the direct interaction approximation are discussed. It is suggested that the Yakhot-Orszag theory yields a lowest order approximate solution of a regularized direct interaction approximation which can be corrected by a simple iterative procedure.

Rubinstein, Robert↗
NASA NTRS1990

Model-data comparisons of moments of nonbreaking shoaling surface gravity waves

The predictions of linear and nonlinear (Boussinesq) shoaling wave models for nonbreaking unidirectional surface gravity waves are compared to field observations, with particular emphasis on quantities that may be important for cross-shore sediment transport. The extensive data sets were obtained on two natural beaches, span water depths between 1 and 10 m, and include incident wave power spectra with narrow, broad, and bimodal shapes. Significant wave heights varied between approximately 30 and 100 cm, and peak periods between approximately 8 and 18 s. The evolution of total variances of sea surface elevation, cross-shore velocity, and horizontal acceleration is modeled at least qualitatively well by both linear and nonlinear theories. Only the nonlinear theory predicts the increasingly asymmetric sea surface elevations and horizontal velocities (pitched-forward wave shapes) and the weaker variation of skewness (difference between crest and trough profiles) which are observed to occur during shoaling. The nonlinear theory also models qualitatively well the large skewed accelerations which occur during the passage of asymmetric waves.

Elgar, Steve↗
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 NTRS2007

Antarctic Circumpolar Current Transport Variability during 2003-05 from GRACE

Gravity Recovery and Climate Experiment (GRACE) gravity data spanning January 2003-November 2005 are used as proxies for ocean bottom pressure (BP) averaged over 1 month, spherical Gaussian caps 500 km in radius, and along paths bracketing the Antarctic Circumpolar Current's various fronts. The GRACE BP signals are compared with those derived from the Estimating the Circulation and Climate of the Ocean (ECCO) ocean modeling-assimilation system, and to a non-Boussinesq version of the Regional Ocean Model System (ROMS). The discrepancy found between GRACE and the models is 1.7 cm H2O (1 cm H2O approx. 1 hPa), slightly lower than the 1.9 cmH2O estimated by the authors independently from propagation of GRACE errors. The northern signals are weak and uncorrelated among basins. The southern signals are strong, with a common seasonality. The seasonal cycle GRACE data observed in the Pacific and Indian Ocean sectors of the ACC are consistent, with annual and semiannual amplitudes of 3.6 and 0.6 cmH2O (1.1 and 0.6 cm H2O with ECCO), the average over the full southern path peaks (stronger ACC) in the southern winter, on days of year 197 and 97 for the annual and semiannual components, respectively; the Atlantic Ocean annual peak is 20 days earlier. An approximate conversion factor of 3.1 Sv (Sv equiv 10(exp 6)cu m/s) of barotropic transport variability per cm H2O of BP change is estimated. Wind stress data time series from the Quick Scatterometer (QuikSCAT), averaged monthly, zonally, and over the latitude band 40(deg)- 65(deg)S, are also constructed and subsampled at the same months as with the GRACE data. The annual and semiannual harmonics of the wind stress peak on days 198 and 82, respectively. A decreasing trend over the 3 yr is observed in the three data types.

Gravity Recovery and Climate Experiment (GRACE)↗