SEARCH · Search NASA

Results for “Boussinesq equation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Solitons induced by boundary conditions from the Boussinesq equation

The behavior of solitons induced by boundary excitation is investigated at various time-dependent conditions and different unperturbed water depths, using the Korteweg-de Vries (KdV) equation. Then, solitons induced from Boussinesq equations under similar conditions were studied, making it possible to remove the restriction in the KdV equation and to treat soliton head-on collisions (as well as overtaking collisions) and reflections. It is found that the results obtained from the KdV and the Boussinesq equations are in good agreement.

Chou, Ru Ling↗

NASA NTRS1997

On the Asymptotic Regimes and the Strongly Stratified Limit of Rotating Boussinesq Equations

Asymptotic regimes of geophysical dynamics are described for different Burger number limits. Rotating Boussinesq equations are analyzed in the asymptotic limit, of strong stratification in the Burger number of order one situation as well as in the asymptotic regime of strong stratification and weak rotation. It is shown that in both regimes horizontally averaged buoyancy variable is an adiabatic invariant for the full Boussinesq system. Spectral phase shift corrections to the buoyancy time scale associated with vertical shearing of this invariant are deduced. Statistical dephasing effects induced by turbulent processes on inertial-gravity waves are evidenced. The 'split' of the energy transfer of the vortical and the wave components is established in the Craya-Herring cyclic basis. As the Burger number increases from zero to infinity, we demonstrate gradual unfreezing of energy cascades for ageostrophic dynamics. The energy spectrum and the anisotropic spectral eddy viscosity are deduced with an explicit dependence on the anisotropic rotation/stratification time scale which depends on the vertical aspect ratio parameter. Intermediate asymptotic regime corresponding to strong stratification and weak rotation is analyzed where the effects of weak rotation are accounted for by an asymptotic expansion with full control (saturation) of vertical shearing. The regularizing effect of weak rotation differs from regularizations based on vertical viscosity. Two scalar prognostic equations for ageostrophic components (divergent velocity potential and geostrophic departure ) are obtained.

Babin, A.↗

NASA NTRS1966

Nonlinear stability of the Boussinesq equations by the method of energy.

Energy method applied to stability of flows governed by nonlinear Boussinesq equations, establishing universal stability region in Reynolds-Rayleigh number plane

REYNOLDS NUMBER↗

DOE OSTI2023

Mean Flow from Phase Averages in the 2D Boussinesq Equations

The atmosphere and ocean are described by highly oscillatory PDEs that challenge both our understanding of their dynamics and their numerical approximation. This paper presents a preliminary numerical study of one type of phase averaging applied to mean flows in the 2D Boussinesq equations that also has application to numerical methods. The phase averaging technique, well-known in dynamical systems theory, relies on a mapping using the exponential operator, and then an averaging over the phase. The exponential operator has connections to the Craya–Herring basis pioneered by Jack Herring to study the fluid dynamics of oscillatory, nonlinear fluid dynamics. In this paper, we perform numerical experiments to study the effect of this averaging technique on the time evolution of the solution. We explore its potential as a definition for mean flows. We also show that, as expected from theory, the phase-averaging method can reduce the magnitude of the time rate of change in the PDEs, making them potentially suitable for time stepping methods.

Wingate, Beth A.↗

DOE OSTI2025

Notes on the non-dimensionalization of the MHD Boussinesq convection equations

The objective of these notes is to demonstrate that the formulation of the MHD Boussinesq convection equations that follows the conventions of the turbulence community used in our previous work is consistent and compatible with the formulation from the convection community used in Busse and Pesch.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Statistical Analysis of the Limiting Dynamics of Two dimensional Boussinesq Turbulent Systems

In the current study, we investigate rotational and stratified turbulent systems within the framework of the Boussinesq equations. The effects of rotation and stratification give rise to a variety of dynamical regimes in atmospheric and oceanic turbulence. To develop a better understanding of the fundamental flow properties through a more straightforward mathematical framework, we commence our research with a simplified two-dimensional Boussinesq model that incorporates rotation and stratification effects. We explore two distinct limiting dynamics, considering pure rotation and pure stratification separately. We examine these intriguing dynamics through numerical investigations, applying the exact solution theory of Boussinesq equations and examining small-scale perturbations to the exact solution under rotation and stratification. We conduct an analysis of statistical quantities associated with the turbulent state variables. The results indicate that, under weak rotation at the equilibrium statistical state, the flow field exhibits vortex flow. However, strong rotation transforms the flow into vertical shear flow. We observe both downscale and upscale energy transfers with a decay rate of k −3 in the absence of external forcing and dissipation. Furthermore, we delve into the study of internal gravity wave mode interactions in the simplified two-dimensional Boussinesq system, which is a crucial aspect of geophysical turbulence due to gravity. Our findings reveal that the time series of the mode coefficients exhibit wave-like interactions between wave modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1990

Recurrence in truncated Boussinesq models for nonlinear waves in shallow water

The rapid spatial recurrence of weakly nonlinear and weakly dispersive progressive shallow-water waves is examined using a numerical integration technique on the discretized and truncated form of the Boussinesq equations. This study primarily examines recurrence in wave fields with Ursell number O(1) and characterizes the sensitivity of recurrence to initial spectral shape and number of allowed frequency modes. It is shown that the rapid spatial recurrence is not an inherent property of the considered Boussinesq systems for evolution distances of 10-50 wavelengths. The main result of the study is that highly truncated Boussinesq models of resonant shallow-water ocean surface gravity waves predict rapid multiple recurrence cycles, but that this is an artifact dependent on the number of allowed modes. For initial conditions consisting of essentially all energy concentrated in a single mode, damping of the recurrence cycles increases as the number of low-power background modes increases. When more than 32 modes are allowed, the recurrence behavior is relatively insensitive to the number of allowed modes.

Elgar, Steve↗

DOE OSTI2025

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗

DOE OSTI2025

A nonhydrostatic formulation for MPAS-Ocean

3D Poisson equation↗

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

Multiple buoyancy driven flows in a vertical cylinder heated from below

The structure of axisymmetric buoyancy-driven convection in a vertical cylinder heated from below is probed by finite element solution of the Boussinesq equations coupled with computed-implemented perturbation techniques for detecting and tracking multiple flows and for determining flow stability. Results are reported for fluids with Prandtl number of one and for cylinders with aspect ratio (Lambda) (defined as the height to radius of the cylinder) between 0.5 and 2.25. Extensive calculations of the neutral stability curve for the static solution and of the nonlinear motions along the bifurcating flow families show a continuous evolution of the primary cellular motion from a single toroidal cell to two and three cells nested radially in the cylinder, instead of the sharp transitions found for a cylinder with shear-free sidewalls. The smooth transitions in flow structure with Rayleigh number and lambda are explained by nonlinear connectivity between the first two bifurcating flow families formed either by a secondary bifurcation point for Lambda or = Lambda * approximately 0.80 or by a limit point for Lambda Lambda *. The transition between these two modes may be described by the theory of multiple limit point bifurcation.

Yamaguchi, Y.↗

NASA NTRS1985

Numerical simulation in three space dimensions of time-dependent thermal convection in a rotating fluid

Three-dimensional time-dependent convection in a plane layer of fluid, uniformly heated from below and subject to vertical shear and to rotation about an axis tilted from the vertical, was simulated by the numerical solution of the Boussinesq equations, including all Coriolis terms. Rotation about a vertical axis produces smaller convection cells with diminished heat fluxes and considerable vorticity. When the rotation axis is tilted from the vertical to represent tropical latitudes, the convection cells become elongated in a N-S direction. Imposed flows with constant vertical shear produce convective rolls aligned with the mean flow. When the rotation vector is tilted from the vertical, the competing effects due to rotation and shear can stabilize the convective motions.

Hathaway, D. H.↗

NASA NTRS1992

Surface tension and buoyancy-driven flow in a non-isothermal liquid bridge

The Navier-Stokes-Boussinesq equations governing the transport of momentum, mass and heat in a nonisothermal liquid bridge with a temperature-dependent surface tension are solved using a vorticity-stream-function formulation together with a nonorthogonal coordinate transformation. The equations are discretized using a pseudo-unsteady semi-implicit finite difference scheme and are solved by the ADI method. A Picard-type iteration is adopted which consists of inner and outer iterative processes. The outer iteration is used to update the shape of the free surface. Two schemes have been used for the outer iteration; both use the force balance normal to the free surface as the distinguished boundary condition. The first scheme involves successive approximation by the direct solution of the distinguished boundary condition. The second scheme uses the artificial force imbalance between the fluid pressure, viscous and capillary forces at the free surface which arises when the boundary condition for force balance normal to the surface is not satisfied. This artificial imbalance is then used to change the surface shape until the distinguished boundary condition is satisfied. These schemes have been used to examine a variety of model liquid bridge situations including purely thermocapillary-driven flow situations and mixed thermocapillary- and bouyancy-driven flow.

Zhang, Yiqiang↗

NASA NTRS2001

Experimental and Numerical Investigation of Buoyancy Driven Convection During PDAMNA Thin Film Growth

This paper presents results from numerical simulations as well as laboratory experiments of buoyancy driven convection in an ampoule under varying heating and gravitational acceleration loadings. The modeling effort in this work resolves the large scale natural convective motion that occurs in the fluid during photodeposition of polydiacetelene films which is due to energy absorbed by the growth solution from a UV source. Consequently, the growth kinetics of the film are ignored in the model discussed here, and also a much simplified ampoule geometry is considered. The objective of this work is to validate the numerical prediction on the strength and structure of buoyancy driven convection that could occur under terrestrial conditions during nonlinear optical film growth. The validation is used to enable a reliable predictive capability on the nature and strength of the convective motion under low gravity conditions. The ampoule geometry is in the form of a parallelepiped with rectangular faces. The numerical results obtained from the solution to the Boussinesq equations show that natural convection will occur regardless of the orientation of the UV source with respect to the gravity vector. The least strong convective motion occurred with the UV beam directed at the top face of the parallelepiped. The strength of the convective motion was found to be almost linearly proportional to the total power of the UV source. Also, it was found that the strength of the convective motion decreased linearly with the gravity due to acceleration. The pattern of the convective flow on the other hand, depended on the source location.

Antar, Basil N.↗

NASA NTRS1987

The thermal-vortex equations

The Boussinesq approximation is extended so as to explicitly account for the transfer of fluid energy through viscous action into thermal energy. Ideal and dissipative integral invariants are discussed, in addition to the general equations for thermal-fluid motion.

Shebalin, John V.↗

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.↗

DOE OSTI2022

Stability in integrable nonlocal nonlinear equations

Recently a variety of nonlocal integrable systems has been introduced that besides fields located at particular space-time points simultaneously also contain fields that are located at different, but symmetrically related, points. Here we investigate different types of soliton solutions with regard to their stability against linear perturbations obtained for the nonlocal version of the Hirota/nonlinear Schrödinger equation and the so-called Alice and Bob versions of the Korteweg-de Vries and Bousinesq equations. In this work, we encounter different types of scenarios: Solition solutions that are linearly stable or unstable and also solutions that change their stability properties depending on the parameter regime they are in.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗