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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Modeling and new equipment definition for the vibration isolation box equipment system

Our MSAD-funded research project is to provide numerical modeling support for the VIBES (Vibration Isolation Box Experiment System) which is an IML2 flight experiment being built by the Japanese research team of Dr. H. Azuma of the Japanese National Aerospace Laboratory. During this reporting period, the following have been accomplished: A semi-consistent mass finite element projection algorithm for 2D and 3D Boussinesq flows has been implemented on Sun, HP And Cray Platforms. The algorithm has better phase speed accuracy than similar finite difference or lumped mass finite element algorithms, an attribute which is essential for addressing realistic g-jitter effects as well as convectively-dominated transient systems. The projection algorithm has been benchmarked against solutions generated via the commercial code FIDAP. The algorithm appears to be accurate as well as computationally efficient. Optimization and potential parallelization studies are underway. Our implementation to date has focused on execution of the basic algorithm with at most a concern for vectorization. The initial time-varying gravity Boussinesq flow simulation is being set up. The mesh is being designed and the input file is being generated. Some preliminary 'small mesh' cases will be attempted on our HP9000/735 while our request to MSAD for supercomputing resources is being addressed. The Japanese research team for VIBES was visited, the current set up and status of the physical experiment was obtained and ongoing E-Mail communication link was established.

Sani, Robert L.↗
NASA NTRS1992

Flow effects in a vertical CVD reactor

A model is presented to simulate the non-Boussinesq flow in a vertical, two-dimensional, chemical vapor deposition reactor under atmospheric pressure. Temperature-dependent conductivity, mass diffusivity, viscosity models, and reactive species mass transfer to the substrate are incorporated. In the limits of small Mach number and small aspect ratio, asymptotic expressions for the flow, temperature, and species fields are developed. Soret diffusion effects are also investigated. Analytical solutions predict an inverse relationship between temperature field and concentration field due to Soret effects. This finding is consistent with numerical simulations, assisting in the understanding of the complex interactions amongst the flow, thermal, and species fields in a chemically reacting system.

Young, G. W.↗
NASA NTRS1991

Frontogenesis driven by horizontally quadratic distributions of density

Attention is given to the quadratic density distribution in a channel, which has been established by Simpson and Linden to be the simplest case of the horizontally nonlinear distribution of fluid density required for the production of frontogenesis. The porous-media and Boussinesq flow models are examined, and their evolution equations are reduced to one-dimensional systems. While both the porous-media and the inviscid/nondiffusive Boussinesq systems exhibit classic frontogenesis behavior, the viscous Boussinesq system exhibits a more complex behavior: boundary-layer effects force frontogenesis away from the lower boundary, and at late times the steepest density gradients are close to mid-channel.

Jacqmin, David↗
NASA NTRS1992

Mixed convective-dynamic roll vortices and their effects on initial wind and temperature profiles

A new nonlinear 14-coefficient spectral model of two-dimensional shallow Boussinesq flow is developed and used to investigate the onset and development of both dynamically and convectively forced boundary-layer rolls. The model is developed to accept arbitrary basic-state wind profiles as dynamic forcing, using an Ekman profile to provide a means for easy comparison with other studies. The results are qualitatively compared with those of previous theoretical and observational investigations. The rolls are shown to significantly alter the initial wind profile in the sense found by Faller and Kaylor (1967) and Brown (1970), but via a mechanism independent of the Coriolis force.

Haack, Tracy↗
NASA NTRS1991

Mixed convective/dynamic roll vortices and their effects on initial wind and temperature profiles

The onset and development of both dynamically and convectively forced boundary layer rolls are studied with linear and nonlinear analyses of a truncated spectral model of shallow Boussinesq flow. Emphasis is given here on the energetics of the dominant roll modes, on the magnitudes of the roll-induced modifications of the initial basic state wind and temperature profiles, and on the sensitivity of the linear stability results to the use of modified profiles as basic states. It is demonstrated that the roll circulations can produce substantial changes to the cross-roll component of the initial wind profile and that significant changes in orientation angle estimates can result from use of a roll-modified profile in the stability analysis. These results demonstrate that roll contributions must be removed from observed background wind profiles before using them to investigate the mechanisms underlying actual secondary flows in the boundary layer. The model is developed quite generally to accept arbitrary basic state wind profiles as dynamic forcing. An Ekman profile is chosen here merely to provide a means for easy comparison with other theoretical boundary layer studies; the ultimate application of the model is to study observed boundary layer profiles. Results of the analytic stability analysis are validated by comparing them with results from a larger linear model. For an appropriate Ekman depth, a complete set of transition curves is given in forcing parameter space for roll modes driven both thermally and dynamically. Preferred orientation angles, horizontal wavelengths and propagation frequencies, as well as energetics and wind profile modifications, are all shown to agree rather well with results from studies on Ekman layers as well as with studies on near-neutral and convective atmospheric boundary layers.

Haack, Tracy↗
NASA NTRS1986

Boundary-layer computational model for predicting the flow and heat transfer in sudden expansions

Fully developed turbulent and laminar flows through symmetric planar and axisymmetric expansions with heat transfer were modeled using a finite-difference discretization of the boundary-layer equations. By using the boundary-layer equations to model separated flow in place of the Navier-Stokes equations, computational effort was reduced permitting turbulence modelling studies to be economically carried out. For laminar flow, the reattachment length was well predicted for Reynolds numbers as low as 20 and the details of the trapped eddy were well predicted for Reynolds numbers above 200. For turbulent flows, the Boussinesq assumption was used to express the Reynolds stresses in terms of a turbulent viscosity. Near-wall algebraic turbulence models based on Prandtl's-mixing-length model and the maximum Reynolds shear stress were compared.

Lewis, J. P.↗
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 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 NTRS1979

Sunspots and the physics of magnetic flux tubes. VI - Convective propulsion. VII - Heat flow in a convective downdraft

The effect of negative aerodynamic drag in an ideal fluid subject to convective instability is considered. It is shown that a cylinder moving in such a fluid is propelled forward in its motion by the convective forces and that the characteristic acceleration time is comparable to the onset time of convective motions in the fluid. It is suggested that convective propulsion plays an important role in the dynamics of flux tubes extending through the surface of the sun. The suppression of the upward heat flow in a Boussinesq convective cell with free upper and lower boundaries by a downdraft is then analyzed. Application to the solar convection zone indicates that downdrafts of 1 to 2 km/s at depths of 1000 to 4000 km beneath the visible surface of the sun are sufficient to reduce the upward heat flux to a small fraction of the ambient value.

Parker, E. N.↗
NASA NTRS2009

Lessons Learned from Numerical Simulations of the F-16XL Aircraft at Flight Conditions

Nine groups participating in the Cranked Arrow Wing Aerodynamics Project International (CAWAPI) project have contributed steady and unsteady viscous simulations of a full-scale, semi-span model of the F-16XL aircraft. Three different categories of flight Reynolds/Mach number combinations were computed and compared with flight-test measurements for the purpose of code validation and improved understanding of the flight physics. Steady-state simulations are done with several turbulence models of different complexity with no topology information required and which overcome Boussinesq-assumption problems in vortical flows. Detached-eddy simulation (DES) and its successor delayed detached-eddy simulation (DDES) have been used to compute the time accurate flow development. Common structured and unstructured grids as well as individually-adapted unstructured grids were used. Although discrepancies are observed in the comparisons, overall reasonable agreement is demonstrated for surface pressure distribution, local skin friction and boundary velocity profiles at subsonic speeds. The physical modeling, steady or unsteady, and the grid resolution both contribute to the discrepancies observed in the comparisons with flight data, but at this time it cannot be determined how much each part contributes to the whole. Overall it can be said that the technology readiness of CFD-simulation technology for the study of vehicle performance has matured since 2001 such that it can be used today with a reasonable level of confidence for complex configurations.

Rizzi, Arthur↗
NASA NTRS2009

What was Learned from Numerical Simulations of F-16XL (CAWAPI) at Flight Conditions

Nine groups participating in the CAWAPI project have contributed steady and unsteady viscous simulations of a full-scale, semi-span model of the F-16XL aircraft. Three different categories of flight Reynolds/Mach number combinations were computed and compared with flight-test measurements for the purpose of code validation and improved understanding of the flight physics. Steady-state simulations are done with several turbulence models of different complexity with no topology information required and which overcome Boussinesq-assumption problems in vortical flows. Detached-eddy simulation (DES) and its successor delayed detached-eddy simulation (DDES) have been used to compute the time accurate flow development. Common structured and unstructured grids as well as individually-adapted unstructured grids were used. Although discrepancies are observed in the comparisons, overall reasonable agreement is demonstrated for surface pressure distribution, local skin friction and boundary velocity profiles at subsonic speeds. The physical modeling, be it steady or unsteady flow, and the grid resolution both contribute to the discrepancies observed in the comparisons with flight data, but at this time it cannot be determined how much each part contributes to the whole. Overall it can be said that the technology readiness of CFD-simulation technology for the study of vehicle performance has matured since 2001 such that it can be used today with a reasonable level of confidence for complex configurations.

Rizzi, Arthur↗
NASA NTRS2014

NASA Trapezoidal Wing Simulation Using Stress-w and One- and Two-Equation Turbulence Models

The Wilcox 2006 stress-omega model (also referred to as WilcoxRSM-w2006) has been implemented in the NASA Langley code CFL3D and used to study a variety of 2-D and 3-D configurations. It predicted a variety of basic cases reasonably well, including secondary flow in a supersonic rectangular duct. One- and two-equation turbulence models that employ the Boussinesq constitutive relation were unable to predict this secondary flow accurately because it is driven by normal turbulent stress differences. For the NASA trapezoidal wing at high angles of attack, the WilcoxRSM-w2006 model predicted lower maximum lift than experiment, similar to results of a two-equation model.

Rodio, J. J.↗
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 NTRS1979

Stellar convection 2: A multi-mode numerical solution for convection in spheres

The convective flow of a self gravitating sphere of Boussinesq fluid for small Reynolds and Peclet numbers is numerically determined. The decomposition of the equations of motion into modes is reviewed and a relaxation method is developed and presented to compute the solutions to these equations. The stable equilibrium flow for a Rayleigh number of 10 to the 4th power and a Prandtl number of 10 is determined. The 2 and 3 dimensional spectra of the kinetic and thermal energies and the convective flux as a function of wavelengths are calculated in terms of modes. The anisotropy of the flow as a function of wavelength is defined.

Marcus, P. S.↗
NASA NTRS1978

Turbulent convection within rapidly rotating superadiabatic fluids with horizontal temperature gradients

A model is developed for small scale thermally driven turbulent convection in rapidly rotating systems. The Boussinesq equations of adiabatic and inviscid flow are perturbed about a basic state which is superadiabatic and which contains a horizontal temperature gradient balanced by a geostrophic wind. The growth rates and length scales of the most unstable modes are determined. It is found that these modes are either zonally or latitudinally symmetric, depending on the latitude and the horizontal temperature gradient. The transports of heat and momentum by the modes are computed. The linear amplitudes are assumed to be limited by shear instability. The model predicts a latitude dependent eddy diffusivity. The applicability of the model to astrophysical and planetary bodies is briefly discussed.

Flasar, F. M.↗
NASA NTRS1980

Normal and abnormal turbulent transport in a stratified medium

The influences of memory effects, coupling between velocity and temperature fluctuations and tensorial transport properties on momentum and heat transfers in turbulent flow which do not follow the Boussinesq relation are examined. It is shown that the memory effect, represented by the Lagrangian of the velocity gradients, can account for the decoupling between the flux and the gradient, while the tensorial properties of the transport coefficients allow a normal Boussinesq-type transfer with memory and anomalous counter-gradient or gradient-less transport.

Tchen, C.-M.↗
NASA NTRS1995

Meso-beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

Mesoscale model simulations provide insight into the complex jet streak adjustments on 11-12 July 1981 that preceded the first of two significant gravity wave events to have been generated over the Rocky Mountains in Montana. Simulations employing a variety of terrain treatments indicate that prior to wave formation, geostrophic adjustment processes modified the structure of the mid-upper tropospheric jet streak by creating secondary jetlets to the southeast of the polar jet streak in proximity to the gravity wave generation region. This simulated restructuring of the mid-upper tropospheric jet streak is the result of a four stage process. During stage 1, the wind adjusts to the mass field as the jet streak exit region propagates into the inflection point between the upstream trough and downstream ridge in the height field. Stage 2 is initiated as the mass field is forced to adjust to the new ageostrophic wind field created during stage 1. Stage 3 is defined by a second geostrophic adjustment process occurring in a similar manner but to the south and east of the adjustment which occurs during stage 1. A low-level mesoscale jetlet is formed during stage 4 in response to the low-level pressure falls that are established during stage 3. The perturbation of this jetlet, caused by orographically-induced adiabatic and diabatic physical processes, is the likely mechanism responsible for the generation of the first and second episode of observed gravity waves. The dynamics responsible for this wave episode are discussed as differential surface sensible heating inducing an orographically-forced mountain-plains solenoid, resulting in the formation of additional mesoscale jetlets and internal gravity waves. Also discussed is how convective latent heating modifies the numerically simulated terrain-induced internal gravity waves, especially their amplitude and phase velocities, which provide better agreement with those wave characteristics observed in nature. Finally, the three-dimensional linear response of a zonally uniform barotropic flow in a vertically unbounded, continuously stratified, Boussinesq atmosphere which is perturbed from geostrophic equilibrium is investigated.

Lin, Yuh-Lang↗