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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Evaluation of a distributed catchment scale water balance model

The validity of some of the simplifying assumptions in a conceptual water balance model is investigated by comparing simulation results from the conceptual model with simulation results from a three-dimensional physically based numerical model and with field observations. We examine, in particular, assumptions and simplifications related to water table dynamics, vertical soil moisture and pressure head distributions, and subsurface flow contributions to stream discharge. The conceptual model relies on a topographic index to predict saturation excess runoff and on Philip's infiltration equation to predict infiltration excess runoff. The numerical model solves the three-dimensional Richards equation describing flow in variably saturated porous media, and handles seepage face boundaries, infiltration excess and saturation excess runoff production, and soil driven and atmosphere driven surface fluxes. The study catchments (a 7.2 sq km catchment and a 0.64 sq km subcatchment) are located in the North Appalachian ridge and valley region of eastern Pennsylvania. Hydrologic data collected during the MACHYDRO 90 field experiment are used to calibrate the models and to evaluate simulation results. It is found that water table dynamics as predicted by the conceptual model are close to the observations in a shallow water well and therefore, that a linear relationship between a topographic index and the local water table depth is found to be a reasonable assumption for catchment scale modeling. However, the hydraulic equilibrium assumption is not valid for the upper 100 cm layer of the unsaturated zone and a conceptual model that incorporates a root zone is suggested. Furthermore, theoretical subsurface flow characteristics from the conceptual model are found to be different from field observations, numerical simulation results, and theoretical baseflow recession characteristics based on Boussinesq's groundwater equation.

Troch, Peter A.↗
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 NTRS1987

The effects of adiabatic and viscous heatings on plumes

The influences of viscous and adiabatic heatings on the transient development of plumes have been assessed by employing the time-dependent convection equations in the extended-Boussinesq limit. Both temperature-dependent and temperature- and pressure-dependent viscosities are employed. For constant thermal expansion of coefficient, the vigor of plume flows is diminished with larger magnitudes of the thermal expansivity. This effect is caused by the local domination of adiabatic decompression over viscous heating inside the developing plume for both the Cartesian and axisymmetric domains. For deep mantle flows a thermal expansivity decreasing with depth is required for thin plume structures to exist. Values of the thermal expansivity contrast, which can promote plume growth, are in accord with estimates from recent studies in seismic equation of state.

Zhao, Wuling↗
NASA NTRS1999

Double Diffusive Convection in Materials Processing

A great number of crystals grown in space are plagued by convective motions which contribute to structural flaws. The character of these instabilities is not well understood but is associated with density variations in the presence of residual gravity (g-jitter). As a specific example, past HgCdTe crystal growth space experiments by Lehoczky and co-workers indicate radial compositional asymmetry in the grown crystals. In the case of HgCdTe the rejected component into the melt upon solidification is HgTe which is denser than the melt. The space grown crystals indicate the presence of three dimensional flow with the heavier HgTe-rich material clearly aligned with the residual gravity (0.55-1.55 micro g) vector. This flow stems from double-diffusive convection, namely, thermal and solutal buoyancy driven flow in the melt. The study of double-diffusive convection is multi-faceted and rather vast. In our investigation, we seek to focus on one specific aspect of this discipline that is of direct relevance to materials processing especially crystal growth, namely, the side ways heating regime. This problem has been widely studied, both experimentally and numerically, in the context of solar ponds wherein the system is characterized by a linear salt (solutal) gradient with an imposed lateral temperature gradient. The induced flow instabilities arise from the wide disparity between the fluid thermal diffusivity and the solute diffusivity. The extension of the analysis to practical crystal growth applications has however not been rigorously made and understood. One subtle but important difference in crystal growth systems is the fact that die system solute gradient is non-linear (typically exponential). Besides, the crystal growth problem has the added complexities of solidification, both lateral and longitudinal thermal gradients and segregation phenomena in systems where binary and ternary compounds are being grown. This paper treats the side ways heating problem alone in a model fluid system. Results from detailed numerical calculations, mainly two dimensional are provided. The interactions between a non-linear solute gradient and an imposed transverse thermal gradient are investigated. The buoyancy effects are treated in the traditional Boussinesq approximation and also in a more complete density formulation to address recent concerns of the first approach especially in simulations of the system response in a reduced gravity environment. Detailed flow, temperature and solute field plots along with heat and mass transfer results are presented in the paper. Implications to practical crystal growth systems as discerned from the modeling results are also explored and reported.

Ramachandra, Narayanan↗
NASA NTRS1985

Instability of time-periodic flows

The instabilities of some spatially and/or time-periodic flows are discussed, in particular, flows with curved streamlines which can support Taylor-Gortler vortices are described in detail. The simplest flow where this type of instability can occur is that due to the torsional oscillations of an infinitely long circular cylinder. For more complicated spatially varying time-periodic flows, a similar type of instability can occur and is spatially localized near the most unstable positions. When nonlinear effects are considered it is found that the instability modifies the steady streaming boundary layer induced by the oscillatory motion. It is shown that a rapidly rotating cylinder in a uniform flow is susceptible to a related type of instability; the appropriate stability equations are shown to be identical to those which govern the instability of a boussinesq fluid of Prandtl number unity heated time periodically from below.

Hall, P.↗
NASA NTRS1987

Instability of time-periodic flows

The instabilities of some spatially and/or time-periodic flows are discussed, in particular, flows with curved streamlines which can support Taylor-Gortler vortices are described in detail. The simplest flow where this type of instability can occur is that due to the torsional oscillations of an infinitely long circular cylinder. For more complicated spatially varying time-periodic flows, a similar type of instability can occur and is spatially localized near the most unstable positions. When nonlinear effects are considered it is found that the instability modifies the steady streaming boundary layer induced by the oscillatory motion. It is shown that a rapidly rotating cylinder in a uniform flow is susceptible to a related type of instability; the appropriate stability equations are shown to be identical to those which govern the instability of a Boussinesq fluid of Prandtl number unity heated time periodically from below.

Hall, Philip↗
NASA NTRS1991

Prediction of jet mean flow structure in support of HSCT noise suppression concepts

The paper describes the application of techniques based on computational fluid dynamics to the simulation of jet flowfields. A solution code for the Reynolds-averaged Navier-Stokes equations is supplemented by conventional two-equation turbulence models based on the Boussinesq approximation. The axisymmetric SCIPVIS code is enhanced with the PARCH and CRAFT codes to examine plug-jet flowfields and imperfectly expanded axisymmetric free round jets. The sensitivity of shock/boundary layer interactions is observed in simulations of the plug case, and the adaptive gridding in the disk region and turbulence levels generated at the triple point are identified as areas in the Mach case that require improvement. Jet-wave structure in the region beyond the first several shock cells can be predicted, and turbulence modeling can be undertaken with respect to improving compressibility, length scale, vorticity, and energy budget. The mean flow structure of imperfectly expanded jets can be studied to develop related noise suppression concepts for the High-Speed Civilian Transport (HSCT).

Sinha, N.↗
NASA NTRS1981

The Lagrangian-mean motions forced by steady, dissipating equatorial waves. I

Waves are treated with a normal mode structure in order to determine the steady mean motion of the atmosphere that can be induced by dissipating equatorial waves. A model is developed which comprises a continuously stratified atmosphere at rest on the equatorial beta-plane. It is assumed that waves are excited by the corrugated bottom and are in a steady state, that dissipation is due to Newtonian cooling and Rayleigh friction, steadiness in wave magnitude is up to the second order, the waves have a long wave length, wave induced mean flows do not affect the waves, mean flows are steady, and dissipation mechanisms for the mean flows are the same as for the waves. Disturbance equations are formulated, along with Eulerian- and Lagrangian-mean flows, and the nonexistence of cross equatorial mean flows is demonstrated. Kelvin waves are shown to possess a Lagrangian-mean meridional circulation which is the same as the Eulerian-mean circulation. In the Boussinesq limit, however, neither the Eulerian- nor the Lagrangian-mean meridional circulations are caused by Kelvin waves. Further examination is made of Rossby-gravity waves and n = 1 westward propagating inertio-gravity waves.

Takahashi, M.↗
NASA NTRS1990

Buoyancy-driven instability in a vertical cylinder: Binary fluids with Soret effect. I - General theory and stationary stability results

The buoyancy-driven instability of a monocomponent or binary fluid completely contained in a vertical circular cylinder is investigated, including the influence of the Soret effect for the binary mixture. The Boussinesq approximation is used, and the resulting linear stability problem is solved using a Galerkin technique. The analysis considers fluid mixtures ranging from gases to liquid metals. The flow structure is found to depend strongly on both the cylinder aspect ratio and the magnitude of the Soret effect. The predicted stability limits are shown to agree closely with experimental observations.

Hardin, G. R.↗
NASA NTRS1992

Numerical study of the small scale structures in Boussinesq convection

Two-dimensional Boussinesq convection is studied numerically using two different methods: a filtered pseudospectral method and a high order accurate Essentially Nonoscillatory (ENO) scheme. The issue whether finite time singularity occurs for initially smooth flows is investigated. The numerical results suggest that the collapse of the bubble cap is unlikely to occur in resolved calculations. The strain rate corresponding to the intensification of the density gradient across the front saturates at the bubble cap. We also found that the cascade of energy to small scales is dominated by the formulation of thin and sharp fronts across which density jumps.

Weinan, E.↗
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 NTRS1994

Local shear instabilities in weakly ionized, weakly magnetized disks

We extend the analysis of axisymmetric magnetic shear instabilities from ideal magnetohydrodynamic (MHD) flows to weakly ionized plasmas with coupling between ions and neutrals caused by collisions, ionization, and recombination. As part of the analysis, we derive the single-fluid MHD dispersion relation without invoking the Boussinesq approximation. This work expands the range of applications of these instabilities from fully ionized accretion disks to molecular disks in galaxies and, with somewhat more uncertainty, to protostellar disks. Instability generally requires the angular velocity to decrease outward, the magnetic field strengths to be subthermal, and the ions and neutrals to be sufficiently well coupled. If ionization and recombination processes can be neglected on an orbital timescale, adequate coupling is achieved when the collision frequency of a given neutral with the ions exceeds the local epicyclic freqency. When ionization equilibrium is maintained on an orbital timescale, a new feature is present in the disk dynamics: in contrast to a single-fluid system, subthermal azimuthal fields can affect the axisymmetric stability of weakly ionized two-fluid systems. We discuss the underlying causes for this behavior. Azimuthal fields tend to be stabilizing under these circumstances, and good coupling between the neutrals and ions requires the collision frequency to exceed the epicyclic frequency by a potentially large secant factor related to the magnetic field geometry. When the instability is present, subthermal azimuthal fields may also reduce the growth rate unless the collision frequency is high, but this is important only if the field strengths are very subthermal and/or the azimuthal field is the dominant field component. We briefly discuss our results in the context of the Galactic center circumnuclear disk, and suggest that the shear instability might be present there, and be responsible for the observed turbulent motions.

Blaes, Omer M.↗
NASA NTRS1991

Time-accurate finite difference solutions to the incompressible Navier-Stokes/energy equations

Two new algorithms for solving coupled Navier-Stokes and energy equations are presented. The algorithms are compared with available Navier-Stokes solutions for the model problem of a driven-cavity-flow over a range of Reynolds numbers. It is noted that the algorithms represent two different implementations of the fractional step approach to the solution of the coupled Navier-Stokes and energy equations using the Boussinesq approximation.

Mirfakhraee, Ali↗
NASA NTRSDate not supplied

Simulation and Modeling of Hypersonic Turbulent Boundary Layers Subject to Adverse Pressure Gradients due to Concave Streamline Curvature

Direct numerical simulations (DNS) of adverse-pressure-gradient turbulent boundary layers over a planar concave wall are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models. The wall geometry and flow conditions of the DNS are representative of the experimental data for a Mach 4.9 turbulent boundary layer that was tested on a two-dimensional planar concave wall model in the high-speed blow-down wind tunnel located at the National Aerothermochemistry Laboratory at Texas A&M University (TAMU). The DNS was validated against the experimental results of TAMU for the same flow conditions and wall geometry. An analysis of the DNS datasets was also conducted to provide an assessment of the validity of Morkovin’s hypothesis and the strong Reynolds analog for turbulence subject to mechanical nonequilibrium. In addition to the DNS results, RANS predictions are obtained by using the Baldwin-Lomax (BL), Spalart-Allmaras (SA), and the k - w SST turbulence models. The comparisons between RANS and DNS showed little impact of an adverse pressure gradient on the accuracy of these models, at least up to an incompressible Clauser pressure gradient parameter of beta(sub inc) 1.22. While the Boussinesq assumption provided reasonable predictions for the Reynolds shear stress, it failed to adequately predict the normal components of the Reynolds stress.

turbulent boundary layers↗
DOE OSTI2020

Quantitative assessment of eddy viscosity rans models for turbulent mixed convection in a differentially heated plane channel

Turbulent mixed convection between two vertical, infinite parallel plates at different temperatures is studied using various two-equation turbulence models. The numerical simulations are performed at a turbulent Reynolds number of Re τ = 150 and a Grashof number of Gr = 9.6 × 10 5 . Comparisons are made against the highly trusted DNS results. Consistent with the DNS approach, the current simulations are performed using constant properties and the Boussinesq approximation to predict the influence of buoyancy. Previous studies have provided assessments of two-equation turbulence models for various scenarios, but often rely on a qualitative “eye” test in order to determine the most appropriate model to predict a given flow. This study aims to provide a new form of quantitative assessment that accounts for both the physics captured by the turbulence model as well as the magnitude of the system response quantities (SRQ) using a modified symmetric mean absolute percent error (SMAPE) method. This method is designed to be approachable to researchers at any level and can be applied to system response quantities from multiple research fields. Uncertainty quantification is also performed to determine the discretization error for each turbulence model. Recommendations are made as to which turbulence models best capture the physics – hydrodynamically and thermally – using both local and global validation metrics. Lastly, a sensitivity analysis is performed on the damping functions used in the most accurate models. This underpins the potential of model developments and adjustments most worth pursuing for buoyant flows. Finally, this framework provides a more physics-based comparative analysis of the selected turbulence models.

42 ENGINEERING↗
DOE OSTI2022

Non-Boussinesq subgrid-scale model with dynamic tensorial coefficients

A major drawback of Boussinesq-type subgrid-scale stress models used in large-eddy simulations is the inherent assumption of alignment between large-scale strain rates and filtered subgrid-stresses. A priori analyses using direct numerical simulation (DNS) data have shown that this assumption is invalid locally as subgrid-scale stresses are poorly correlated with the large-scale strain rates [J. Bardina, J. Ferziger, and W. Reynolds, Improved subgrid-scale models for large-eddy simulation, in Proceedings of the 13th Fluid and Plasmadynamics Conference, AIAA (1980); C. Meneveau and K. Katz, Scale-invariance and turbulence models for large-eddy simulation, Ann. Rev. Fluid Mech. 32, 1 (2000)]. In the present work, a new, non-Boussinesq subgrid-scale model is presented where the model coefficients are computed dynamically. Some previous non-Boussinesq models have observed issues in providing adequate dissipation of turbulent kinetic energy [e.g., Bardina et al., Proceedings of the 13th Fluid and Plasmadynamics Conference (1980); R. A. Clark, J. Ferziger, and W.C. Reynolds. Evaluation of subgrid-scale models using an accurately simulated turbulent flow, J. Fluid Mech. 91, 1 (1979); S. Stolz and N. A. Adams, An approximate deconvolution procedure for large-eddy simulation, Phys. Fluids 11, 1699 (1999)]; however, the present model is shown to provide sufficient dissipation using dynamic coefficients. Modeled subgrid-scale Reynolds stresses satisfy the consistency requirements of the governing equations for large-eddy simulation (LES), vanish in laminar flow and at solid boundaries, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The new model, referred to as the dynamic tensor-coefficient Smagorinsky model (DTCSM), has been tested in simulations of canonical flows: decaying and forced homogeneous isotropic turbulence, and wall-modeled turbulent channel flow at high Reynolds numbers. The results show favorable agreement with DNS data. It has been shown that DTCSM offers similar predictive capabilities as the dynamic Smagorinsky model for canonical flows. In order to assess the performance of DTCSM in more complex flows, wall-modeled simulations of high Reynolds number flow over a Gaussian bump (Boeing speed bump) exhibiting smooth-body flow separation are performed. Predictions of surface pressure and skin friction, compared against DNS and experimental data, show improved accuracy from DTCSM in comparison to existing static coefficient (Vreman) and dynamic Smagorinsky model. The computational cost of performing LES with this model is up to 15% higher than the dynamic Smagorinsky model.

42 ENGINEERING↗
NASA NTRS1974

Finite-amplitude thermal convection in a spherical shell

The properties of finite-amplitude thermal convection for a Boussinesq fluid contained in a spherical shell are investigated. All nonlinear terms are retained in the equations, and both axisymmetric and nonaxisymmetric solutions are studied. The velocity is expanded in terms of poloidal and toroidal vectors. Spherical surface harmonics resolve the horizontal structure of the flow, but finite differences are used in the vertical. With a few modifications, the transform method developed by Orszag (1970) is used to calculate the nonlinear terms, while Green's function techniques are applied to the poloidal equation and diffusion terms.

Young, R. E.↗
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.↗