Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NATURAL CONVECTION”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

A Revised Method of Presenting Wavenumber-Frequency Power Spectrum Diagrams That Reveals the Asymmetric Nature of Tropical Large-scale Waves

The popular method of presenting wavenumber-frequency power spectrum diagrams for studying tropical large-scale waves in the literature is shown to give an incomplete presentation of these waves. The so-called "convectively-coupled Kelvin (mixed Rossby-gravity) waves" are presented as existing only in the symmetric (antisymmetric) component of the diagrams. This is obviously not consistent with the published composite/regression studies of "convectively-coupled Kelvin waves," which illustrate the asymmetric nature of these waves. The cause of this inconsistency is revealed in this note and a revised method of presenting the power spectrum diagrams is proposed. When this revised method is used, "convectively-coupled Kelvin waves" do show anti-symmetric components, and "convectively-coupled mixed Rossby-gravity waves (also known as Yanai waves)" do show a hint of symmetric components. These results bolster a published proposal that these waves be called "chimeric Kelvin waves," "chimeric mixed Rossby-gravity waves," etc. This revised method of presenting power spectrum diagrams offers a more rigorous means of comparing the General Circulation Models (GCM) output with observations by calling attention to the capability of GCMs in correctly simulating the asymmetric characteristics of the equatorial waves.

Chao, Winston C.↗
NASA NTRS1994

Principal modes of variation of rain-rate probability distributions

Radar or satellite observations of an area generate sequences of rain-rate maps. From a gridded map a histogram of rain rates can be obtained representing the relative areas occupied by rain rates of various strengths. The histograms vary with time as precipitating systems in the area evolve and decay and amounts of convective and stratiform rain in the area change. A method of decomposing the histograms into linear combinations of a few empirical distributions with time-dependent coefficients is developed, using principal component analysis as a starting point. When applied to a tropical Atlantic dataset (GATE), two distributions emerge naturally from the analysis, resembling stratiform and convective rain-rate distributions in that they peak at low and high rain rates, respectively. The two 'modes' have different timescales and only the high-rain-rate mode has a statistically significant diurnal cycle. The ability of just two modes to describe rain variabiltiy over an area can explain why methods of estimating area-averaged rain rate from the area covered by rain rates above a certain threshold are so successful.

Bell, Thomas L.↗
NASA NTRS1987

The dynamo dilemma

The recent determination that the angular velocity Omega of the sun declines downward through the convective zone raises serious questions about the nature of the solar dynamo. The principal qualitative features of the sun are the azimuthal fields that migrate toward the equator in association with an oscillating poloidal field which reverses at about the time of maximum appearance of bipolar magnetic regions. If Omega decreases downward, or is negligible, the horizontal gradient in Omega produces a dynamo with some of these essential characteristics. There is reason to think that the dynamo is confined to the lower half of the convective zone, where alpha has the opposite sign from the usual (alpha of greater than 0 in the northern hemisphere) producing equatorward migration but reversing the sign of the associated poloidal field. Meridional circulation may play an essential role in shaping the dynamo. At the present time it is essential to measure Omega accurately and determine the nature of the meridional circulation.

Parker, E. N.↗
NASA NTRS1994

Empirical convection models for northward IMF

It is clear that polar cap convection during times of northward Interplanetary Magnetic Field (IMF) is more structured and of lower mean speed than at times of southward IMF. This, coupled with the fact that the polar cap is smaller, means that empirical models are more difficult to construct with certainty. It is also clear that sunward flow deep in the polar cap is often observed, but its connection with the rest of the flow pattern is controversial. At present, empirical models are of three types: 'statistical' models wherein data from different days but with similar IMF conditions are averaged together; 'pattern recognition' models, which are built up by examining individually hundreds of passes to derive a 'typical' pattern which embodies features frequently observed; and 'assimilative' models, which use data of different types and from as many locations as possible, but all taken at the same time, in order to derive a snapshot (or series of snapshots) of the entire pattern. Each type of model has its own difficulties. Statistical models, by their very nature, smooth out flow features (e.g. the convection reversal, and the locus of sunward flow deep in the polar cap) which are not found at precisely the same invariant latitudes and magnetic local times on different days. Pattern recognition models are better at reproducing small-scale features, but the large-scale pattern can be a matter of interpretation. Assimilative models (such as AMIE) hold out the best hope for creating instantaneous, global convection patterns; however, the analysis technique tends to be most irregular (and least reliable) in the regions which are not well covered by in situ data. It appears that, at least at times, a four cell model with sunward flow at the highest and lowest latitudes, and antisunward flow in between, is consistent with the observations. At other times, the observations may be consistent with a two-cell convection pattern, but which includes significant meanders within the polar cap.

Moses, Julie J.↗
NASA NTRS2011

Tropical and Subtropical Cloud Transitions in Weather and Climate Prediction Models: The GCSS/WGNE Pacific Cross-Section Intercomparison (GPCI)

A model evaluation approach is proposed in which weather and climate prediction models are analyzed along a Pacific Ocean cross section, from the stratocumulus regions off the coast of California, across the shallow convection dominated trade winds, to the deep convection regions of the ITCZ the Global Energy and Water Cycle Experiment Cloud System Study/Working Group on Numerical Experimentation (GCSS/ WGNE) Pacific Cross-Section Intercomparison (GPCI). The main goal of GPCI is to evaluate and help understand and improve the representation of tropical and subtropical cloud processes in weather and climate prediction models. In this paper, a detailed analysis of cloud regime transitions along the cross section from the subtropics to the tropics for the season June July August of 1998 is presented. This GPCI study confirms many of the typical weather and climate prediction model problems in the representation of clouds: underestimation of clouds in the stratocumulus regime by most models with the corresponding consequences in terms of shortwave radiation biases; overestimation of clouds by the 40-yr ECMWF Re-Analysis (ERA-40) in the deep tropics (in particular) with the corresponding impact in the outgoing longwave radiation; large spread between the different models in terms of cloud cover, liquid water path and shortwave radiation; significant differences between the models in terms of vertical cross sections of cloud properties (in particular), vertical velocity, and relative humidity. An alternative analysis of cloud cover mean statistics is proposed where sharp gradients in cloud cover along the GPCI transect are taken into account. This analysis shows that the negative cloud bias of some models and ERA-40 in the stratocumulus regions [as compared to the first International Satellite Cloud Climatology Project (ISCCP)] is associated not only with lower values of cloud cover in these regimes, but also with a stratocumulus-to-cumulus transition that occurs too early along the trade wind Lagrangian trajectory. Histograms of cloud cover along the cross section differ significantly between models. Some models exhibit a quasi-bimodal structure with cloud cover being either very large (close to 100%) or very small, while other models show a more continuous transition. The ISCCP observations suggest that reality is in-between these two extreme examples. These different patterns reflect the diverse nature of the cloud, boundary layer, and convection parameterizations in the participating weather and climate prediction models.

Teixeira, J.↗
NASA NTRS1993

Order of accuracy of QUICK and related convection-diffusion schemes

This report attempts to correct some misunderstandings that have appeared in the literature concerning the order of accuracy of the QUICK scheme for steady-state convective modeling. Other related convection-diffusion schemes are also considered. The original one-dimensional QUICK scheme written in terms of nodal-point values of the convected variable (with a 1/8-factor multiplying the 'curvature' term) is indeed a third-order representation of the finite volume formulation of the convection operator average across the control volume, written naturally in flux-difference form. An alternative single-point upwind difference scheme (SPUDS) using node values (with a 1/6-factor) is a third-order representation of the finite difference single-point formulation; this can be written in a pseudo-flux difference form. These are both third-order convection schemes; however, the QUICK finite volume convection operator is 33 percent more accurate than the single-point implementation of SPUDS. Another finite volume scheme, writing convective fluxes in terms of cell-average values, requires a 1/6-factor for third-order accuracy. For completeness, one can also write a single-point formulation of the convective derivative in terms of cell averages, and then express this in pseudo-flux difference form; for third-order accuracy, this requires a curvature factor of 5/24. Diffusion operators are also considered in both single-point and finite volume formulations. Finite volume formulations are found to be significantly more accurate. For example, classical second-order central differencing for the second derivative is exactly twice as accurate in a finite volume formulation as it is in single-point.

Leonard, B. P.↗
DOE OSTI2021

Pronghorn: A Multidimensional Coarse Mesh Application for Advanced Reactor Thermal-Hydraulics

This paper presents an overview of Pronghorn, a multiscale thermal-hydraulic (T/H) application developed by Idaho National Laboratory and the University of California, Berkeley. Pronghorn, built on the open-source finite element Multiphysics Object-Oriented Simulation Environment (MOOSE), leverages state-of-the-art physical models, numerical methods, and nonlinear solvers to deliver fast-running advanced reactor T/H simulation capabilities within a modern software engineering environment. This work summarizes the physical models, multiphysics and multiscale coupling, and numerical discretization in Pronghorn with emphasis on our initial target application to pebble bed reactors (PBRs). A diverse set of applications are shown to depressurized natural circulation in the SANA experiments, forced convection in the Pebble Bed Modular Reactor, three-dimensional (3-D)/one-dimensional coupling of Pronghorn and RELAP-7 systems T/H for loop analysis in the High Temperature Reactor Power Module, and forced convection in the Mark-1 Pebble Bed Fluoride-Salt-Cooled High-Temperature Reactor. A multiphysics coupling of Pronghorn, RELAP-7, and Griffin deterministic neutronics for a gas-cooled PBR demonstrates the capability of the MOOSE framework for reactor design calculations. These applications highlight the verification and validation underlying Pronghorn’s software development while emphasizing features that improve upon capabilities offered by legacy tools in areas such as 3-D unstructured meshing, physics modeling, and multiphysics coupling.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Computational Fluid Dynamics and Heat Transfer Modeling of a Dimpled Heat Exchanger

Multiphysics and analytical calculations were conducted for a heat exchanger with passive, natural circulation flow. A glycol/water working fluid convects the heat to a dimpled heat exchanger shell, which subsequently transfers the heat to the soil, which acts as the ultimate heat sink. Because the system is fully-passive, it is not subject to the expenses, maintenance, and mechanical breakdowns associated with moving parts. Density, heat capacity, and thermal conductivity material properties were measured for various soil samples, and subsequently included as input for the soil heat conduction model. The soil model was coupled to a computational fluid dynamics (CFD) heat exchanger model that included the dynamic Smagorinsky large eddy simulation and k- omega turbulence models. The analysis showed that the fluid dynamics and heat transfer models worked properly, albeit at a slow pace. Nevertheless, the coupled CFD/heat conduction simulation ran long enough to determine a key parameter—the amount of heat conducted from the heat exchanger to the ground. This unique performance value, along with experimental data, was used as input for stand-alone, fast-running CFD models, as well as boundaries to obtain solutions to partial differential equations for soil heat conduction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1980

Sound attenuation by liners in a blown flap environment

Sound propagation through a hot wall-jet flow over an absorbing wall is studied. The radiated sound field subject to the influence of flow convection and refraction is evaluated, and the nature of acoustic attenuation attributable to a sound absorbing liner is determined. Using a two-dimensional model, the noise field under the aircraft is also determined, and a slug-flow model is used to describe the influence of flow, density, and temperature on acoustic sources in jets. Results show significant changes in the radiated source due to the interference phenomenon, and a good absorber has the potential of changing the sound pressure range of variation to unity. A liner is also found to increase or decrease sound pressure, depending on the frequency.

Parthasarathy, S. P.↗
NASA NTRS1995

Analysis of helioseismic wave fields to examine horizontal structures

A preliminary inversion procedure that carries out a local area analysis on simulated oscillation data is presented and discussed. The procedure is carried out in order to deduce two dimensional subsurface structures in the horizontal plane, representative of thermal variations, as a function of depth. The aim is the evaluation of an inversion procedure that utilizes information gained from the phase distortion occurring in artificially generated acoustic waves in order to determine the subsurface thermal structure. These distortions would naturally occur as a direct consequence of convective motions in the solar interior.

Julien, K. A.↗
NASA NTRS2001

Large Lewis No. Edge-Flame Instabilities

Edge-flames play an important role in a number of microgravity investigations, and in the general study of flames. Examples include the candle-flame experiments carried out on board both the Space Shuttle and the Mir Space Station; the flame-spread-over-liquid work carried out by H. Ross and W. Sirignano amongst others and lifted turbulent diffusion flames. In all of these configurations a local two-dimensional flame structure can be identified which looks like a flame-sheet with an edge, and these structures exhibit dynamical behavior which characterizes them and distinguishes them from ad hoc 2D flame structures. Edge-flames can exist in both a non-premixed context (edges of diffusion flames) and in a premixed context (edges of deflagrations), but the work reported here deals with the edges of diffusion flames. It is particularly relevant, we believe, to oscillations that have been seen in both the candle-flame context, and the flame-spread-over-liquid context. These oscillations are periodic edge-oscillations (in an appropriate reference frame), sans oscillation of the trailing diffusion flame. It is shown that if the Lewis number of the fuel is sufficiently large (the Lewis number of the oxidizer is taken to be 1), and the Damkohler number is sufficiently small, oscillating-edge solutions can be found. Oscillations are encouraged by an on-edge convective flow and the insertion of a cold probe, discouraged by an off-edge convective flow. In the present work, the nature of these oscillations is examined in more depth, using a variety of numerical strategies.

Buckmaster, J.↗
NASA NTRSDate not supplied

Maintaining the Continuity of the Global Climate Data Record of Lightning from Space

Lightning is a natural hazard, an indicator of changes in convective storm behavior across multiple spatial and temporal scales, and a major influence on climate-affecting chemicals such as ozone and methane. Lightning also has substantial impacts on the near-space environment. As such, the World Meteorological Organization (WMO) has designated lightning as an Essential Climate Variable (ECV); that is, a measurement that is really important for understanding Earth’s climate. Spaceborne lightning observations are the gold standard for measuring global lightning, as the detection efficiency is more spatially uniform compared to ground-based lightning measurements. In addition, spaceborne is the only practical method for measuring total lightning (i.e., both intracloud and cloud-to-ground) at the global scale, as ground-based networks have limited ability to detect intracloud lightning at long distances (e.g., over oceans, far from land-based sensors). Since 1995, with only a two-year gap (2015-2017), lightning has been observed globally using a combination of Optical Transient Detector (OTD) and Lightning Imaging Sensor (LIS) instruments. However, the year 2023 is the planned end for the International Space Station (ISS) LIS mission, and there is no designated low-Earth orbit (LEO) successor yet identified. Since 2017, geostationary (GEO) observations of lightning have become increasingly available. These include the U.S. Geostationary Lightning Mappers (GLMs), China’s Lightning Mapping Imager (LMI), and Europe’s Lightning Imager (LI). While GEO provides continuous, quasi-hemispheric coverage of lightning and thunderstorms, there are current gaps preventing global coverage, particularly in portions of Asia and the western Pacific, as well as high-latitude regions. Moreover, integration of the GEO constellation into a Climate Data Record (CDR) of lightning is a major challenge only recently begun to be addressed by the WMO. Furthermore, recent work has demonstrated some deficiencies in the global record of lightning from space. For example, small lightning flashes are very common in severe weather, but are difficult to observe using the present set of spaceborne lightning imagers. Moreover, some types of lightning are best observed using frequencies other than the commonly used near-infrared (777-nm), including ultraviolet and the radio spectrum. Finally, it is not possible to assess the vertical distribution of lightning globally, even though this measurement is highly relevant to thunderstorm behavior as well as lightning’s effects on greenhouse gases and the near-space environment. This presentation will provide a current assessment of lightning measurements from space, will identify significant current and future gaps in these observations, and will evaluate potential solutions for improving the overall CDR of global lightning.

Timothy Lang↗
DOE OSTI2026

Local convection characteristics of inline arrangement of Kagome-shaped unit cells in a square duct

Transient liquid crystal thermography experiments have been conducted to determine detailed convective heat transfer coefficients at the endwalls of lattice-frame configurations based on Kagome-shaped unit cells. Kagome unit cells with porosity of 0.88 were arranged in an inline manner, where two such arrangements were studied. In the first arrangement, a total of ten unit cells were placed next to each other along the streamwise direction resulting in a continuous configuration. In the 2nd arrangement, alternate unit cells from the continuous configuration were dropped, resulting in a discrete configuration which featured a total of 5 Kagome unit cells. Due to the asymmetric nature of the strut connections within the unit cell, the convective heat transfer coefficient maps were determined for the two opposite walls where the struts meet the endwalls. Transient liquid crystal experiments were conducted for Reynolds number ranging between 10,000 and 30,000. Here, the study was focused on the developing nature of flow along the streamwise direction and the local convection characteristics for continuous and discrete arrangement of unit cells. For the continuous configuration, the Nusselt number ratios (N u /Nu 0 ) varied between 3.15–3.46 and 2.75–2.89 for the two walls A and B, respectively. For the discrete configuration, convective heat transfer coefficients varied between 2.56–2.71 and 2.09–2.40 for the two walls. Kagome unit cells have the potential to be fabricated through inexpensive manufacturing routes such as wire-woven method and these unit cells find their applications in the areas which require different heat transfer levels on opposite walls.

Kagome↗
NASA NTRS2006

Can High-resolution WRF Simulations Be Used for Short-term Forecasting of Lightning?

A number of research teams have begun to make quasi-operational forecast simulations at high resolution with models such as the Weather Research and Forecast (WRF) model. These model runs have used horizontal meshes of 2-4 km grid spacing, and thus resolved convective storms explicitly. In the light of recent global satellite-based observational studies that reveal robust relationships between total lightning flash rates and integrated amounts of precipitation-size ice hydrometeors in storms, it is natural to inquire about the capabilities of these convection-resolving models in representing the ice hydrometeor fields faithfully. If they do, this might make operational short-term forecasts of lightning activity feasible. We examine high-resolution WRF simulations from several Southeastern cases for which either NLDN or LMA lightning data were available. All the WRF runs use a standard microphysics package that depicts only three ice species, cloud ice, snow and graupel. The realism of the WRF simulations is examined by comparisons with both lightning and radar observations and with additional even higher-resolution cloud-resolving model runs. Preliminary findings are encouraging in that they suggest that WRF often makes convective storms of the proper size in approximately the right location, but they also indicate that higher resolution and better hydrometeor microphysics would be helpful in improving the realism of the updraft strengths, reflectivity and ice hydrometeor fields.

Goodman, S. J.↗
NASA NTRS2002

A Method to Retrieve Rainfall Rate over Land from TRMM Observations

Tropical Rainfall Measuring Mission (TRMM) Precipitation Radar (PR) observations over mesoscale convective systems (MCSs) reveal that there are localized maxima in the rain rate with a scale of about 10 to 20 km that represent thunderstorms (Cbs). Some of these Cbs are developing or intense, while others are decaying or weak. These Cbs constitute only about 20 % of the rain area of a given MCS. Outside of Cbs, the average rain rate is much weaker than that within Cbs. From an analysis of the PR data, we find that the spatial distribution of rain and its character, convective or stratiform, is highly inhomogeneous. This complex nature of rain exists on a scale comparable to that of a Cb. The 85 GHz brightness temperature, T85, observations of the TRMM Microwave Imager (TMI) radiometer taken over an MCS reflect closely the PR rain rate pattern over land. Local maxima in rain rate shown by PR are observed as local minima in T85. Where there are no minima in T85, PR observations indicate there is light rain. However, the TMI brightness temperature measurements (Tbs) have poor ability to discriminate convective rain from stratiform rain. For this reason, a TMI rain retrieval procedure that depends primarily on the magnitude of Tbs performs poorly. In order to retrieve rain rate from TMI data on land one has to include the spatial distribution information deduced from the T85 data in the retrieval method. Then, quantitative estimation of rain rate can be accomplished. A TMI rain retrieval method developed along these lines can yield estimates of rain rate and its frequency distribution which agree closely with that given by PR. We find the current TRMM project TMI (Version 5) rain retrieval algorithm on land could be improved with the retrieval scheme developed here. To support the conceptual frame work of the rain retrieval method developed here, a theoretical analysis of the TMI brightness temperatures in convective and stratiform regions is presented.

Prabhakara, C.↗
NASA NTRS1980

High-latitude ion convection in the nighttime F region

A large data base of Atmosphere Explorer C measurements of the bulk ion velocity vector in the high-latitude F region has been examined. The concept of shear and rotational flow reversals is expanded upon to determine the extent to which the boundary between antisunward and sunward convection (i.e., a polar cap boundary) is an electrostatic equipotential. A most important factor in characterizing the nature of the large-scale geometrical configuration of the nighttime F region convection is the location and extent of quasi-equipotential regions of the convection reversal boundary.

Heelis, R. A.↗
NASA NTRS2000

Onset of Convection Due to Surface Tension Variations in Multicomponent and Binary Fluid Layers

Under certain conditions, such as in thin liquid films or microgravity, surface tension variations along a free surface can induce convection. Convection onset due to surface tension variation is important to many terrestrial technological processes in addition to microgravity materials processing applications. Examples include coating, drying crystallization, solidification, liquid surface contamination, and containerless processing. In double-diffusive and multicomponent systems, the spatial variations of surface tension are associated with two or more stratifying agencies, respectively. For example, both temperature and species (concentration) gradients are associated with convection in the solidification of binary alloys or salt ponds. The direction of the two (or more) gradients has a profound effect on the nature of the flow at or slightly beyond the onset of convection. Our recent work at the NASA Lewis Research Center focused on characterizing surface-tension-induced onset of convection, often referred to as Marangoni-Benard convection. Exact solutions for the stationary neutral stability of multicomponent fluid layers with interfacial deformation were derived. These solutions also permit the computation of a boundary curve that separates the long and finite wavelength instabilities. Computing points along this boundary using the exact solution (when possible) is more efficient than the typical numerical approaches, such as finite difference or spectral methods. Above the curve, a long wavelength instability was predicted, suggesting that convection would occur principally through one large flow cell in the layer, whereas below the curve, finite wavelength instabilities occur which suggest multiple finite-sized circulation cells. For many common liquids with layer depths greater than 100 mm, finite wave instability is predicted under terrestrial conditions; however, with little exception, long wavelength instability is predicted in microgravity for the identical fluid systems.

Skarda, J. Raymond Lee↗
NASA NTRS1974

Thermal instability and energy relations in the convective envelopes of slowly rotating stars

The nonstationary convection that arises due to thermal instability in the envelopes of lower Main Sequence stars (and the sun in particular) is investigated. The convection that arises in the envelopes as a result of thermal instability is of nonstationary turbulent nature. It is shown that construction of a consistent model of the zones of turbulent convection in the stars requires the use of the complete heat flux equation with consideration of the terms that appear because of turbulence. Expressions for these additional terms in terms of averaged characteristics of the motion and eddy viscosity coefficients are derived on the basis of the Prandtl-Wasiutinsky theory.

Yavorskaya, I. M.↗