Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boussinesq approximation”

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
DOE OSTI2021

Parallelized POD-based suboptimal economic model predictive control of a state-constrained Boussinesq approximation

Motivated by an energy efficient building application, we want to optimize a quadratic cost functional subject to the Boussinesq approximation of the Navier-Stokes equations and to bilateral state and control constraints. Since the computation of such an optimal solution is numerically costly, we design an efficient strategy to compute a sub-optimal (but applicationally acceptable) solution with significantly reduced computational effort. We employ an economic Model Predictive Control (MPC) strategy to obtain a feedback control. The MPC sub-problems are based on a linear-quadratic optimal control problem subjected to mixed control and state constraints and a convection-diffusion equation, reduced with proper orthogonal decomposition. Finally, to solve each sub-problem, we apply a primal-dual active set strategy. The method can be fully parallelized, which enables the solution of large problems with real-world parameters.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

The stability of a compressible stratified shear layer

The stability of a shear layer under the effect of gravity is investigated using the compressible magnetohydrodynamic (MHD) equations, including an effective gravity term to represent the curvature effects of the flow and magnetic field line geometry. A general eigenmode equation is derived for a two-dimensional MHD fluid, and an energy-principle analysis to explain the effect of compressibility on the critical Richardson number is presented. For the case of a hyperbolic tangent shear flow and exponential density profile, it was found that, in the Boussinesq approximation, the compressibility raises the critical Richardson number from 1/4 to as much as 1/2, with the exact value depending on the value of the magnetic field at infinity. Under approximation of a strong asymptotic magnetic field, without invoking the Boussinesq approximation, it is shown both analytically and numerically that the density gradient terms cause the shear instability to be dispersive. The long-wavelength stability boundary for the Richardson number J = 0 is characterized by a normalized phase velocity c =

Wang, Z.↗
DOE OSTI2025

Assessment of integral models for non-Boussinesq lazy plumes using numerical simulations

Integral modelling of turbulent buoyant plumes is crucial for rapid predictions of plume characteristics. While the governing equations are typically derived using self-similarity and a Boussinesq approximation, these assumptions may not hold for plumes originating from finite-area sources with large density ratios. Here, this work evaluates the accuracy of integral-scale models for non-Boussinesq lazy plumes using high-fidelity numerical simulations of turbulent helium plumes. We analyse the plume kinematics by computing vertical fluxes, plume radius and radial profiles, establishing some disparities between common practice and physical accuracy. We identify how the definition of the plume radius changes the perception of the plume structure when the flow is not self-similar and derive a relationship between the flux-based and threshold-based definitions without requiring self-similarity. We then examine the plume dynamics by evaluating the source terms from the governing plume equations. Our results support neglecting diffusive and viscous effects but emphasise the importance of the mean pressure gradient, even in the self-similar regime. Two coefficients need to be modelled: the well-known entrainment coefficient and the lesser-known momentum correction coefficient, which is a correction required for the momentum equation to account for self-similar and slender approximations. The momentum correction coefficient is found to be approximately constant and slightly greater than the assumed value of 1. The standard entrainment coefficient models perform well up to a local Richardson number three times the asymptotic value but overpredict entrainment for larger Richardson numbers. We propose a correction using the known finite limit of entrainment at infinite Richardson number.

Meehan, Michael Alexander [Sandia National Laborat↗
NASA NTRS1992

A unified approach for numerical simulation of viscous compressible and incompressible flows over adiabatic and isothermal walls

A new formulation (including the choice of variables, their non-dimensionalization, and the form of the artificial viscosity) is proposed for the numerical solution of the full Navier-Stokes equations for compressible and incompressible flows with heat transfer. With the present approach, the same code can be used for constant as well as variable density flows. The changes of the density due to pressure and temperature variations are identified and it is shown that the low Mach number approximation is a special case. At zero Mach number, the density changes due to the temperature variation are accounted for, mainly through a body force term in the momentum equation. It is also shown that the Boussinesq approximation of the buoyancy effects in an incompressible flow is a special case. To demonstrate the new capability, three examples are tested. Flows in driven cavities with adiabatic and isothermal walls are simulated with the same code as well as incompressible and supersonic flows over a wall with and without a groove. Finally, viscous flow simulations of an oblique shock reflection from a flat plate are shown to be in good agreement with the solutions available in literature.

Hafez, M.↗
NASA NTRS1992

The effects of Venusian mantle convection with multiple phase transitions

Recently there was a flurry of activities in studying the effects of phase transitions in the Earth's mantle. From petrological and geophysical considerations, phase-transitions would also play an important role in venusian dynamics. The basic differences between the two planets are the surface boundary conditions, both thermally and mechanically. In this vein we have studied time-dependent mantle convection with multiple phase transitions and depth-dependent thermal expansivity (alpha is approximately rho(exp -6)), based on high-pressure and temperature measurements. Both the olivine-spinel and spinel-perovskite transitions were simulated by introducing an effective thermal expansivity, as described. Used together with the extended Boussinesq Approximation this method serves as a powerful tool to examine the effects of phase transitions on convection at relatively low computational costs.

Steinbach, V.↗
DOE OSTI2022

Isogeometric large-eddy simulations of turbulent particle-laden flows

In recent years, isogeometric analysis (IGA) has attracted significant attention from the computational mechanics community due to its ability to integrate design and analysis. Besides, IGA is also a higher-order discretization technique for solving partial differential equations, showing high approximation capability per degree of freedom. In this paper, we extend the application realm of IGA to particle-laden flows based on Eulerian–Eulerian description that couples Navier–Stokes equations with a density transport equation through a Boussinesq approximation. The coupled systems are solved by using quadratic non-uniform rational B-spline (NURBS) functions and a recently developed residual-based variational multiscale (VMS) formulation, which introduces coupling between the fine velocity scales and density equation residuals. We deploy the proposed approach to perform large-eddy simulations (LES) of dilute particle-laden flows over a flat surface at Reynolds number = 10,000. We compare the simulation results against direct numerical simulation (DNS) results from the literature. We find that combining VMS and IGA, the proposed approach enables accurate prediction of a wide range of flow/particle statistics with a relatively lower mesh resolution.

Mathematics↗
DOE OSTI2022

Eddy diffusivity operator in homogeneous isotropic turbulence

Here, we use the recently developed macroscopic forcing method to compute the scale-dependent eddy diffusivity characterizing ensemble-averaged scalar and momentum transport in incompressible homogeneous isotropic turbulence. For scales larger than the energy containing eddies, eddy diffusivity is found to be constant and consistent with the Boussinesq approximation. However, for small scales eddy diffusivity is found to vanish inversely proportional to the wave number. Behavior at all scales is reasonably captured by a nonlocal eddy diffusivity operator modeled as D/$\sqrt{I - l^2∇^2}$, where D is the eddy diffusivity in the Boussinesq limit, and l is a constant on the order of the large-eddy length. These results present a direct measurement of eddy diffusivity in turbulence with implications in turbulence modeling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Impacts of Permeability Heterogeneity and Background Flow on Supercritical CO 2 Dissolution in the Deep Subsurface

Motivated by CO 2 capture and sequestration (CCS) design considerations, we consider the coupled effects of permeability heterogeneity and background flow on the dissolution of a supercritical CO 2 lens into an underlying deep, confined aquifer. We present the results of a large-scale Monte Carlo simulation study examining the interaction of background flow rate and three parameters describing multi-Gaussian log-permeability fields: mean, variance, and correlation length. Hundreds of high-resolution simulations were performed using the PFLOTRAN finite volume software to model CO 2 dissolution in a kilometer-scale aquifer over 1,000 years. Predictive dimensionless scaling relationships relating CO 2 dissolution rate to heterogeneity statistics, Rayleigh (Ra) and Péclet (Pe) numbers were developed for both gravitationally dominated free convection to background flow-dominated forced convection regimes. An empirical criterion, Pe = Ra 3/4 , was discovered for regime transition. All simulations converged quickly to a quasi-steady, approximately linear dissolution rate. However, this rate displayed profound variability between permeability field realizations sharing the same heterogeneity statistics, even under mild permeability heterogeneity. In general, increased heterogeneity was associated with a lower mean and higher variance of dissolution rate, undesirable from a CCS design perspective. The relationship between dissolution rate and background flow was found to be complex and nonlinear. Dimensionless scaling relationships were uncovered for a number of special cases. Results call into question the validity of the Boussinesq approximation in the context of modest-to-high background flow rates and the general applicability of numerical simulations without background flow.

42 ENGINEERING↗
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 NTRS2002

Application of the generalized vertical coordinate ocean model for better representing satellite data

It is found that two adaptive parametric functions can be introduced into the basic ocean equations for utilizing the optimal or hybrid features of commonly used z-level, terrain- following, isopycnal, and pressure coordinates in numerical ocean models. The two parametric functions are formulated by combining three techniques: the arbitrary vertical coordinate system of Kasahara (1 974), the Jacobian pressure gradient formulation of Song (1 998), and a newly developed metric factor that permits both compressible (non-Boussinesq) and incompressible (Boussinesq) approximations. Based on the new formulation, an adaptive modeling strategy is proposed and a staggered finite volume method is designed to ensure conservation of important physical properties and numerical accuracy. Implementation of the combined techniques to SCRUM (Song and Haidvogel1994) shows that the adaptive modeling strategy can be applied to any existing ocean model without incurring computational expense or altering the original numerical schemes. Such a generalized coordinate model is expected to benefit diverse ocean modelers for easily choosing optimal vertical structures and sharing modeling resources based on a common model platform. Several representing oceanographic problems with different scales and characteristics, such as coastal canyons, basin-scale circulation, and global ocean circulation, are used to demonstrate the model's capability for multiple applications. New results show that the model is capable of simultaneously resolving both Boussinesq and non-Boussinesq, and both small- and large-scale processes well. This talk will focus on its applications of multiple satellite sensing data in eddy-resolving simulations of Asian Marginal Sea and Kurosio. Attention will be given to how Topex/Poseidon SSH, TRMM SST; and GRACE ocean bottom pressure can be correctly represented in a non- Boussinesq model.

numerical ocean model generalized coordinate syste↗
DOE OSTI2024

The inviscid incompressible limit of Kelvin–Helmholtz instability for plasmas

The Kelvin–Helmholtz Instability (KHI) is an interface instability that develops between two fluids or plasmas flowing with a common shear layer. KHI occurs in astrophysical jets, solar atmosphere, solar flows, cometary tails, planetary magnetospheres. Two applications of interest, encompassing both space and fusion applications, drive this study: KHI formation at the outer flanks of the Earth’s magnetosphere and KHI growth from non-uniform laser heating in magnetized direct-drive implosion experiments. Here, we study 2D KHI with or without a magnetic field parallel to the flow. We use both the GAMERA code, which solves the compressible Euler equations, and the STRATOSPEC code, which solves the Navier-Stokes equations under the Boussinesq approximation, coupled with the magnetic field dynamics. GAMERA is a global three-dimensional MHD code with high-order reconstruction in arbitrary nonorthogonal curvilinear coordinates, which is developed for a large range of astrophysical applications. STRATOSPEC is a three-dimensional pseudo-spectral code with an accuracy of infinite order (no numerical diffusion). Magnetized KHI is a canonical case for benchmarking hydrocode simulations with extended MHD options. An objective is to assess whether or not, and under which conditions, the incompressibility hypothesis allows to describe a dynamic compressible system. For comparing both codes, we reach the inviscid incompressible regime, by decreasing the Mach number in GAMERA, and viscosity and diffusion in STRATOSPEC. Here, we specifically investigate both single-mode and multi-mode initial perturbations, either with or without magnetic field parallel to the flow. The method relies on comparisons of the density fields, 1D profiles of physical quantities averaged along the flow direction, and scale-by-scale spectral densities. We also address the triggering, formation and damping of filamentary structures under varying Mach number or Atwood number, with or without a parallel magnetic field. Comparisons show very satisfactory results between the two codes. The vortices dynamics is well reproduced, along with the breaking or damping of small-scale structures. We end with the extraction of growth rates of magnetized KHI from the compressible regime to the incompressible limit in the linear regime assessing the effects of compressibility under increasing magnetic field. The observed differences between the two codes are explained either from diffusion or non-Boussinesq effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2026

Numerical Simulation of a Natural Convection–Driven Air-Cooled Reactor Cavity Cooling System Experiment

Ensuring the efficient removal of decay heat from the reactor vessel is essential for the safety of advanced reactor technologies. Several Generation-IV concepts incorporate variations in the reactor vessel cooling systems to achieve this objective. High-temperature gas-cooled reactors utilize a reactor cavity cooling system (RCCS), a passive ex-vessel system designed to operate without active components or external power during accident conditions. The RCCS removes decay heat primarily through radiative and convective heat transfer mechanisms. Here, this study presents a comprehensive validation of a computational fluid dynamics Reynolds-averaged Navier-Stokes model for the University of Wisconsin-Madison air-cooled RCCS facility. Validation was conducted for both high- and low-power natural convection cases under a uniform heating profile. Near-wall resolution was found to be critical for accurately modeling natural convection in the RCCS; employing an all-𝑦 + wall treatment resulted in wall temperature discrepancies exceeding 50 °⁢𝐶 compared to a wall-resolved mesh. Thermal-hydraulic behaviors under natural and forced convection conditions were compared within the heated cavity and RCCS. A turbulence model sensitivity analysis indicated that low-Reynolds number k-ɛ, k-ω shear stress transport (SST), and Reynolds stress transport models produce similar wall temperature predictions. A buoyancy modeling sensitivity study revealed that the Boussinesq approximation significantly underpredicted thermal-hydraulic behavior in the RCCS. Based on these findings, modeling recommendations are provided. The validated data set along with identified sensitivities refine the modeling of natural convection in the RCCS. The information produced by this study supports RCCS design, optimization, and safety evaluations, enabling the calibration and verification of reduced-order thermal-hydraulic models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Numerical investigation of two-dimensional buoyancy-driven eddies in liquid metal magnetohydrodynamic flows in breeding blankets *

Abstract Lead–lithium flows are key features in the design of tokamak breeding blanket concepts such as the dual-coolant lead–lithium (DCLL). Since they flow under magnetic fields, they are affected by magnetohydrodynamic (MHD) effects. The neutron flux originating in the tokamak plasma heats the breeding blanket channels in a non-uniform manner, inducing buoyancy forces in the liquid metal. Buoyancy may become a source of quasi-two-dimensional (Q2D) turbulence, and the appearance of eddies may affect the transport of heat and tritium across the blanket. Blankets characterized by high-speed liquid metal flows (such as DCLL) will need ceramic insulating walls to reduce the MHD-related pressure drop in the channels. In our simulations, we have used the Q2D model proposed by Sommeria and Moreau (SM82) which is especially suitable for modeling electrically insulating channel flows. On top of that, we have modelled buoyancy forces in the momentum equation using the Oberbeck–Boussinesq approximation. In this work, we include a validation of the implemented Q2D model in buoyancy-driven cases and identify a thermohydraulic configuration that promotes the generation of eddies and the accumulation of tritium. We also show the results provided by our post-processing tool based on the bi-dimensional fast Fourier transform for eddy detection and characterization. We complete our investigation by performing an initial assessment of how eddies can accumulate tritium in the breeding blanket. We conclude by discussing the relevance of preventing accumulation of tritium within the blanket and suggest a possible solution.

Suarez, Daniel (ORCID:0000000251740100)↗
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 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↗
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.↗