Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boussinesq”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
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

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 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 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 NTRS1993

Steady, three-dimensional, internally heated convection

Numerical calculations have been carried out of steady, symmetric, 3D modes of convection in internally heated, infinite Prandtl number, Boussinesq fluids at a Rayleigh number of 1.4 x 10 exp 4 in a spherical shell with inner/outer radius of 0.55 and in a 3 x 3 x 1 rectangular box. Multiple patterns of convection occur in both geometries. In the Cartesian geometry, the patterns are dominated by cylindrical cold downflows and a broad hot upwelling. In the spherical geometry the patterns consist of cylindrical cold downwellings centered either at the vertices of a tetrahedron or the centers of the faces of a cube. The cold downflow cylinders are immersed in a background of upwelling within which there are cylindrical hot concentrations (plumes) and hot halos around the downflows. The forced hot upflow return plumes of internally heated spherical convection are fundamentally different from the buoyancy-driven plumes of heated from below convection.

Schubert, G.↗
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 NTRS1955

The Theories of Turbulence

The theory of turbulence reached its full growth at the end of the 19th century as a result of the work by Boussinesq and Reynolds. It then underwent a long period of stagnation which ended under the impulse given to it by the development of wind tunnels caused by the needs of aviation. Numerous researchers, attempted to put Reynolds' elementary statistical theory into a more precise form. During the war, some isolated scientists - von Weizsacker and Heisenberg in Germany, Kolmogoroff in Russia, Onsager in the U.S.A. - started a program of research. By a system of assumptions which make it possible to approach the structure of turbulence in well-defined limiting conditions quantitatively, they obtained a certain number of laws on the correlations and the spectrum. Since the late reports have improved the mathematical language of turbulence, it was deemed advisable to start with a detailed account of the mathematical methods applicable to turbulence, inspired at first by the work of the French school, above all for the basic principles, then the work of the foreigners, above all for the theory of the spectrum.

MECHANICS, STATISTICAL↗
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 NTRS1993

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

The Control-B simulation experiment was designed to improve on the Control-A experiment performed with the GMASS model. This experiment addressed several inadequacies with the first smooth terrain numerical simulation by including: (1) increased nested-grid resolution to better define the simulated gravity waves, (2) increased horizontal diffusion to remove outflow boundary condition noise, and (3) the use of reanalyzed rawinsonde data and surface observations in the initial state to increase the definition of the observed jet streak as well as other low-level features. A smoothed-terrain dry simulation with the nested-grid GMASS model has revealed many important aspects of the processes which resulted in the generation of gravity waves in the region and time when and where they were observed. However, the vertical structure, number, and characteristics of the waves are still quite different from observed waves as diagnosed thus necessitating future improved simulations. However, this control simulation has produced substantial insight into processes which occur on many spatial scales over a 30 hour time period thus allowing one to draw promising inferences as to the mechanisms for the complex process which occurred in nature during the CCOPE case study. The theoretical aspects of the project have focus on understanding the nature of the ageostrophic circulations which are produced in idealized models of the atmosphere in which the troposphere is modeled in one of two ways. The first model assumes that the lower atmosphere can be represented as a single layer of homogeneous fluid whose upper surface is free to exhibit vertical displacement. Two-dimensional internal convergence (divergence) occurring during the adjustment to an asymptotic equilibrium state from an ageostrophic initial state whose momentum structure is representative of a midlatitude localized zonal wind anomaly will cause the free upper surface of the homogeneous atmosphere to rise (fall), and therefore the response can be viewed as being physically three-dimensional. The second model assumes that the troposphere can be represented by an unbounded continuously stratified Boussinesq fluid of constant Brunt-Vaisala frequency N, where the vertical gradient of the basic state potential temperature profile allows for the existence of vertically propagating internal inertia-gravity waves.

Lin, Yuh-Lang↗
NASA NTRS1994

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

An in-depth analysis of observed gravity waves and their relationship to precipitation bands over the Montana mesonetwork during the 11-12 July 1981 CCOPE case study indicated two episodes of coherent waves. While geostrophic adjustment, shearing instability, and terrain were all implicated separately or in combination as possible wave generation mechanisms, the lack of upper-air data within the wave genesis region made it difficult to define the genesis processes from observations alone. The first part of this paper, 3D Numerical Modeling Studies of Terrain-Induced Mass/Momentum Perturbations, employs a mesoscale numerical model to help diagnose the intricate early wave generation mechanisms during the first observed gravity wave episode. The meso-beta scale numerical model is used to study various simulations of the role of multiple geostrophic adjustment processes in focusing a region for gravity wave genesis. The second part of this paper, Linear Theory and Theoretical Modeling, investigates the response of non-resting rotating homogeneous and continuously stratified Boussinesq models of the terrestrial atmosphere to temporally impulsive and uniformly propagating three-dimensional localized zonal momentum sources representative of midlatitude jet streaks. The methods of linear perturbation theory applied to the potential vorticity (PV) and wave field equations are used to study the geostrophic adjustment dynamics. The total zonal and meridional wind perturbations are separated into geostrophic and ageostrophic components in order to define and follow the evolution of both the primary and secondary mesocirculations accompanying midlatitude jetogenesis forced by geostrophic adjustment processes. This problem is addressed to help fill the gap in understanding the dynamics and structure of mesoscale inertia-gravity waves forced by geostrophic adjustment processes in simple two-dimensional quiescent current systems and those produced by mesoscale numerical models simulating the orographic and diabatic perturbation of three-dimensional quasi-geostrophically balanced synoptic scale jet streaks associated with complex baroclinic severe storm producing environments.

Lin, Yuh-Lang↗
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↗
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 NTRS1995

Visualization of 2-D and 3-D Tensor Fields

In previous work we have developed a novel approach to visualizing second order symmetric 2-D tensor fields based on degenerate point analysis. At degenerate points the eigenvalues are either zero or equal to each other, and the hyperstreamlines about these points give rise to trisector or wedge points. These singularities and their connecting hyperstreamlines determine the topology of the tensor field. In this study we are developing new methods for analyzing and displaying 3-D tensor fields. This problem is considerably more difficult than the 2-D one, as the richness of the data set is much larger. Here we report on our progress and a novel method to find, analyze and display 3-D degenerate points. First we discuss the theory, then an application involving a 3-D tensor field, the Boussinesq problem with two forces.

Hesselink, Lambertus↗
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 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 NTRS1997

Visualization of 2-D and 3-D Tensor Fields

In previous work we have developed a novel approach to visualizing second order symmetric 2-D tensor fields based on degenerate point analysis. At degenerate points the eigenvalues are either zero or equal to each other, and the hyper-streamlines about these points give rise to tri-sector or wedge points. These singularities and their connecting hyper-streamlines determine the topology of the tensor field. In this study we are developing new methods for analyzing and displaying 3-D tensor fields. This problem is considerably more difficult than the 2-D one, as the richness of the data set is much larger. Here we report on our progress and a novel method to find , analyze and display 3-D degenerate points. First we discuss the theory, then an application involving a 3-D tensor field, the Boussinesq problem with two forces.

Hesselink, Lambertus↗
NASA NTRS1999

Modeling Studies of PVT Growth of ZnSe: Current Status and Future Course

Bulk growth of wide band gap II-VI semiconductors by physical vapor transport (PVT) has been developed and refined over the past several years at NASA Marshall Space Flight Center. Results from a modeling study of PVT crystal growth of ZnSe are reported in this paper. The PVT process is numerically investigated using a two-dimensional formulation of the governing equations and associated boundary conditions. Both the incompressible Boussinesq approximation and a compressible model are tested to determine the influence of gravity on the process and to discern the differences between the two approaches. The influence of a residual gas is included in the models. The results show that both the incompressible and compressible approximations provide comparable results and the presence of a residual gas tends to measurably reduce the mass flux in the system. Detailed flow, thermal and concentration profiles are provided. The simulations show that the Stefan flux dominates the system flow field and the subtle gravitational effects can be gauged by subtracting this flux from the calculated profiles. Shear flows, due to solutal buoyancy, of the order of 50 microns/s for the liorizont,-d growth orientation and 10 microns/s for the vertical orientation are predicted. Whether these flows can fully account for the observed gravity related growth morphological effects and inhomogeneous solute and dopant distributions is a matter of conjecture. A template for future modeling efforts in this area is suggested which incorporates a mathematical approach to the tracking of the growth front based on energy of formation concepts.

Ramachandran, N.↗
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↗