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

Emergence of a nocturnal low-level jet from a broad baroclinic zone

An analytical model is presented for the generation of a Blackadar-like nocturnal low-level jet in a broad baroclinic zone. The flow is forced from below (flat ground) by a surface buoyancy gradient and from above (free atmosphere) by a constant pressure gradient force. Diurnally-varying mixing coefficients are specified to increase abruptly at sunrise and decrease abruptly at sunset. With attention restricted to a surface buoyancy that varies linearly with a horizontal coordinate, the Boussinesq-approximated equations of motion, thermal energy, and mass conservation reduce to a system of one-dimensional equations that can be solved analytically. Sensitivity tests with southerly jets suggest that (i) stronger jets are associated with larger decreases of the eddy viscosity at sunset (as in Blackadar theory), (ii) the nighttime surface buoyancy gradient has little impact on jet strength, and (iii) for pure baroclinic forcing (no free-atmosphere geostrophic wind), the nighttime eddy diffusivity has little impact on jet strength, but the daytime eddy diffusivity is very important and has a larger impact than the daytime eddy viscosity. The model was applied to a jet that developed in fair weather conditions over the Great Plains from southern Texas to northern South Dakota on 1 May 2020. The ECMWF Reanalysis v5 (ERA5) for the afternoon prior to jet formation showed that a broad north-south-oriented baroclinic zone covered much of the region. The peak model-predicted winds were in good agreement with ERA5 winds and lidar data from the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) central facility in north-central Oklahoma.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Entrainment, Detrainment, and Dilution of Dry and Moist Atmospheric Thermals

Here this study examines the entrainment, detrainment, and dilution of dry and moist (cloud) atmospheric thermals in large-eddy simulations. In a neutrally stable environment (with respect to dry dynamics), moist thermals have an increase in radius R with thermal height z t (α ≡ dR/dz t ) about 4 times smaller compared to dry thermals when density stratification is considered and ~2.4 times smaller without density stratification (i.e., applying the Boussinesq approximation). An analytic expression relating α to several dimensionless parameters is derived from the thermal impulse–circulation relation to clarify the factors impacting α. This expression shows that the difference in buoyancy structure between moist and dry thermals, with buoyancy concentrated in the central cores of moist thermals owing to latent heating, explains their smaller spreading rates. Individual contributions of entrainment and detrainment are analyzed using a direct parcel-based approach in the simulations. Moist thermals have similar fractional detrainment but much smaller fractional entrainment rates compared to dry thermals, consistent with the differences in α. Despite having smaller α, moist thermals are similarly dilute (quantified by a passive tracer) as dry thermals because of their greater mixing efficiency with the environment. Thus, moist thermals are substantially dilute but expand much less in size/volume as they rise compared to dry thermals. The α values for moist thermals in (dry) neutral and statically stable environments are similar, but fractional entrainment and especially detrainment rates are greater in the stable environment. Large detrainment rates are associated with a breakdown of the broader thermal vortex ring structure, especially with low environmental relative humidity, attributed in part to evaporation and buoyancy reversal.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The extended auxiliary equation mapping method to determine novel exact solitary wave solutions of the nonlinear fractional PDEs

Abstract In this paper, some new nonlinear fractional partial differential equations (PDEs) have been considered.Three models are including the space-time fractional-order Boussinesq equation, space-time (2 + 1)-dimensional breaking soliton equations, and space-time fractional-order SRLW equation describe the behavior of these equations in the diverse applications. Meanwhile, the fractional derivatives in the sense of β -derivative are defined. Some fractional PDEs will convert to the considered ordinary differential equations by the help of transformation of β -derivative. These equations are analyzed utilizing an integration scheme, namely, the extended auxiliary equation mapping method. The different kinds of traveling wave solutions, solitary, topological, dark soliton, periodic, kink, and rational, fall out as a by-product of this scheme. Finally, the existence of the solutions for the constraint conditions is also shown. The outcome indicates that some fractional PDEs are used as a growing finding in the engineering sciences, mathematical physics, and so forth.

Engineering↗
DOE OSTI2021

Deployment of the Finite Volume Method in Pronghorn for Gas and Salt cooled Pebble Bed Reactors

This report summarizes the activities related to ”Complete FHR and HTGR pebble bed simulator, including initial validation” funded by the NEAMS thermal-hydraulics focus area. The activity revolves around the coarse-mesh thermal-hydraulics code Pronghorn and its application to gas and salt cooled Pebble bed reactor (PBR). The main difference between gas and salt cooled PBR from a thermal-hydraulics perspective is the fluid. To address the difference in fluid behavior, two separate approaches are implemented in the MOOSE Navier Stokes module: 1) a Boussinesq approximation and 2) a fully compressible formulation. The developed finite volume method capabilities are used for improving pre-existing gas-cooled and salt-cooled pebble-bed reactor models. A steady-state, multiphysics gas-cooled pebble-bed reactor model is created that couples the equilibrium core depletion capability developed in previous work, and the finite volume method capability developed for this report. The salt-cooled pebble-bed reactor model is upgraded to use the incompressible finite volume method capability and then extended to three spatial dimensions. Finally, several verification-and-validation exercises performed with Pronghorn are documented using the verification-and-validation report of the MooseDoc system. The goal of this effort to document the verification-and-validation level of Pronghorn and improve stakeholder confidence in the results obtained with Pronghorn.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

On the approximations underpinning fast nuclear cloud models

Fast nuclear cloud models have been used for many decades to predict cloud rise and growth following a nuclear detonation for both emergency response and debris collections mission planning. We present derivations and discussions of the approximations made by such fast nuclear cloud models in common use. In particular, the models we discuss all fundamentally rely on an empirical entrainment parameter to inexpensively model the cloud rise and growth of a buoyant bubble due to entrainment of an ambient air mass on its surface. We show the particular circumstances under which these models are equivalent as well as apply a virtual mass correction throughout them. Aside from some minor differences in representation and the use of the Boussinesq approximation, these models are all very similar and necessitate the specification of an entrainment parameter to accurately determine the cloud rise.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1973

An analytical and numerical study of the Martian planetary boundary layer over slopes.

A one-dimensional model of the Martian planetary boundary layer over sloping terrain is analyzed under a variety of conditions. Analytical results for the steady and diurnal components of the temperature and wind fields are found when a Boussinesq model with a Newtonian cooling law is considered. These results form a basis for understanding the numerical results which include more realistic representations for the heating and parametrizations for the eddy transfer of momentum and heat. The diurnal boundary layer thickness is determined primarily by radiative processes, and the amplitudes of the wind and temperature oscillations are found to depend in an important way on the latitude and slope magnitude. Typically, oscillations in the temperature of plus or minus 15 K and in the upslope wind of plus or minus 25 m/sec are found 1 km above a Martian slope of 0.005.

Blumsack, S. L.↗
NASA NTRS1973

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

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

Bohachevsky, I. O.↗
NASA NTRS1973

Radiative damping of trapped gravity waves in the solar atmosphere.

The radiative damping of trapped gravity waves in an optically thin atmosphere is studied for a stratified Boussinesq fluid. The character of the atmospheric eigenmodes depends on the distribution of the Brunt-Vaisala frequency N and the radiative relaxation time tau. The calculations for simple layer models show that if N x tau is large over some finite fraction of the trapping region, then modes of long lifetime can exist. In order to suppress gravity waves entirely, it is necessary that N x tau be less than or equal to 1 over the entire trapping region. Qualitative application of the results to the solar atmosphere leads to the conclusion that gravity wave eigenmodes of the solar atmosphere, although damped, are by no means eliminated by radiative effects.

Clark, P. A.↗
NASA NTRS1973

Numerical simulation of small-scale thermal convection in the atmosphere

A Boussinesq system is integrated numerically in three dimensions and time in a study of nonhydrostatic convection in the atmosphere. Simulation of cloud convection is achieved by the inclusion of parametrized effects of latent heat and small-scale turbulence. The results are compared with the cell structure observed in Rayleigh-Benard laboratory conversion experiments in air. At a Rayleigh number of 4000, the numerical model adequately simulates the experimentally observed evolution, including some prominent transients of a flow from a randomly perturbed initial conductive state into the final state of steady large-amplitude two-dimensional rolls. At Rayleigh number 9000, the model reproduces the experimentally observed unsteady equilibrium of vertically coherent oscillatory waves superimposed on rolls.

Somerville, R. C. J.↗
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 NTRS1974

Laser Doppler velocimeter system simulation for sensing aircraft wake vortices

A hydrodynamic model of aircraft vortex wakes in an irregular wind shear field near the ground is developed and used as a basis for modeling the characteristics of a laser Doppler detection and vortex location system. The trailing vortex sheet and the wind shear are represented by discrete free vortices distributed over a two-dimensional grid. The time dependent hydrodynamic equations are solved by direct numerical integration in the Boussinesq approximation. The ground boundary is simulated by images, and fast Fourier Transform techniques are used to evaluate the vorticity stream function. The atmospheric turbulence was simulated by constructing specific realizations at time equal to zero, assuming that Kolmogoroff's law applies, and that the dissipation rate is constant throughout the flow field. The response of a simulated laser Doppler velocimeter is analyzed by simulating the signal return from the flow field as sensed by a simulation of the optical/electronic system.

Thomson, J. A. L.↗
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 NTRS1975

Pressure- and buoyancy-driven thermal convection in a rectangular enclosure

Results are presented for unsteady laminar thermal convection in compressible fluids at various reduced levels of gravity in a rectangular enclosure which is heated on one side and cooled on the opposite side. The results were obtained by solving numerically the equations of conservation for a viscous, compressible, heat-conducting, ideal gas in the presence of a gravitational body force. The formulation differs from the Boussinesq simplification in that the effects of variable density are completely retained. A conservative, explicit, time-dependent, finite-difference technique was used and good agreement was found for the limited cases where direct comparison with previous investigations was possible. The solutions show that the thermally induced motion is acoustic in nature at low levels of gravity and that the unsteady-state rate of heat transfer is thereby greatly enhanced relative to pure conduction. The nonlinear variable density profile skews the streamlines towards the cooler walls but is shown to have little effect on the steady-state isotherms.

Spradley, L. W.↗
NASA NTRS1976

Studies of earth simulation experiments

The low gravity environment of earth orbit offers the potential for performing experiments involving baroclinic Geophysical Fluid Dynamics (GFD) on spherical surfaces. These experiments in turn have the potential for providing deeper understanding of large scale planetary and solar circulations. However, to perform these experiments, one requires an experimental technique whereby a radially directed body force can be generated to simulate a radial gravitational force field. One viable technique is the use of dielectric fluids with temperature dependent dielectric permittivity in a radially directed electric field. Application of the Boussinesq approximation to the equations of motion for this system and restrictions on the size of certain electrodynamic terms in the energy equations yields a set of equations which are analogous to the equations of motions of geophysical systems like the earth's atmosphere on term by term basis. The theoretical design of GFD experiments for performance in earth orbit are described along with results of preliminary tests of a prototype.

Hart, J. E.↗
NASA NTRS1977

Convection of viscous fluids: Energetics, self-similarity, experiments, geophysical applications and analogies

The main results were the formulas for the mean convection velocities, of a viscous fluid and for the mean temperature difference in the bulk of the convecting fluid. These were obtained: by scaling analysis of the Boussinesq equations, by analysis of the energetics of the process, and by using similarity and dimensional arguments. The last approach defines the criteria of similarity and allows the proposition of some self-similarity hypotheses. By several simple new ways, an expression for the efficiency coefficient gamma of the thermal convection was also obtained. An analogy is pointed out between non-turbulent convection of a viscous fluid and the structure of turbulence for scales less than Kolmogorov's internal viscous microscale of turbulence.

Golitsyn, G. S.↗
NASA NTRS1977

Numerical solutions of single-mode convection equations

In the Boussinesq approximation, single-mode equations describing thermal convection are constructed by expanding the fluctuating velocity and temperature fields in a complete set of functions (or planforms) of the horizontal coordinates and retaining just one term. Numerical solutions of the single-mode equations are investigated, chief consideration being given to hexagonal planforms. Extensive surveys of steady solutions are presented for various Rayleigh numbers, Prandtl numbers, and horizontal wavenumbers. The dependences on Rayleigh number and Prandtl number at very large Rayleigh number are in satisfactory agreement with the results of asymptotic expansions.

Toomre, J.↗
NASA NTRS1979

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

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

Pfister, L.↗