Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRATIFICATION”

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 325 records · Page 18

Thermal stability effects on the structure of the velocity field above an air-water interface

Mean velocity and turbulence measurements are described for turbulent flows above laboratory water waves, under various wind and thermal stratification conditions. Experimental results, when presented in the framework of Monin-Obukhov (1954) similarity theory, support local scaling based on evaluation of stratification effects at the same nondimensional distance from the mean water surface. Such scaling allows an extension of application of the above theory to the outer region of the boundary layer. Throughout the fully turbulent region, ratios of mean velocity gradients, eddy viscosities, and turbulence intensities under nonneutral and neutral conditions correlate well with the parameter z/Lambda (Lambda being a local Obukhov length and z the vertical coordinate of the mean air flow) and show good agreement with established field correlations. The influence of stratification on the wind-stress coefficient can be estimated from an empirical relationship in terms of its value under neutral conditions and a bulk Richardson number.

Papadimitrakis, Y. A.↗

The correspondence of surface climate parameters with satellite and terrain data

One of the goals of the research was to develop a ground sampling stragegy for calibrating remotely sensed measurements of surface climate parameters. The initial sampling strategy involved the stratification of the terrain based on important ancillary surface variables such as slope, exposure, insolation, geology, drainage, fire history, etc. For a spatially heterogeneous population, sampling error is reduced and efficiency increased by stratification of the landscape into more homogeneous sub-areas and by employing periodic random spacing of samples. These concepts were applied in the initial stratification of the study site for the purpose of locating and allocating instrumentation.

Dozier, Jeff↗

Covariance of biophysical data with digital topograpic and land use maps over the FIFE site

This paper discusses the biophysical stratification of the FIFE site, implementation of the stratification utilizing geographic information system methods, and validation of the stratification with respect to field measurements of biomass, Bowen ratio, soil moisture, and the greenness vegetation index (GVI) derived from TM satellite data. Maps of burning and topographic position were significantly associated with variation in GVI, biomass, and Bowen ratio. The stratified design did not significantly alter the estimated site-wide means for surface climate parameters but accounted for between 25 and 45 percent of the sample variance depending on the variable.

Davis, F. W.↗

Stratified turbulence in the atmosphere and oceans - A new subgrid model

Turbulence in stably and unstably stratified media is studied, and the results are used to construct a subgrid simulation (SGS) model for use in large eddy simulation (LES). It was found that, although the assumption of inertiality of the subgrid scales is reasonable for the case of unstable stratification, it ceases to be so in the case of stable stratification, where gravity removes kinetic energy from the eddies: the generation of gravity waves becomes the dominant physical process, with dissipation relegated to higher wavenumbers. Preliminary results show that the total kinetic energy dissipation length scale increases with stability, in accordance with LES results but in disagreement with Deardorff's model that suggests a decrease of all dissipation scales in presence of stratification.

Canuto, V. M.↗

Transient development of perturbations in stratified shear flow

Transient development of perturbations in inviscid stratified shear flow is investigated. Use is made of closed form analytic solutions that allow concise identification of optimally growing plane-wave solutions for the case of an unbounded flow with constant shear and stratification. For the case of channel flow, variational techniques are employed to determine the optimally growing disturbances. The maximum energy growth attained over a specific time interval decreases continuously with increasing stratification, and no special significance attaches to Ri = 0.25. Indeed, transient growth can be substantial even for Ri = O(1). A general lower bound on the energy growth attained by an optimal perturbation in a stratified flow over a given time interval is the square root of the growth attained by the corresponding perturbation in unstratified flow. Enhanced perturbation persistence is found for mean-flow stratification lying in the range Ri between 0.1 and 0.3. Small but finite perturbations in mean flow with Ri less than 0.4 produce regions with locally negative total density gradient, which are expected to overturn. Although the perturbations are of wave form, buoyancy fluxes mediate transfer between perturbation kinetic and potential energy during transient development, thus implying that buoyancy flux is not a determinative diagnostic for distinguishing between waves and turbulence in stratified flows.

Farrell, Brian F.↗

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

Numerical Modeling Studies of Wake Vortex Transport and Evolution Within the Planetary Boundary Layer

In support of the wake vortex effect of the Terminal Area Productivity program, we have put forward four tasks to be accomplished in our proposal. The first task is validation of two-dimensional wake vortex-turbulence interaction. The second task is investigation of three-dimensional interaction between wake vortices and atmospheric boundary layer (ABL) turbulence. The third task is ABL studies. The, fourth task is addition of a Klemp-Durran condition at the top boundary for TASS model. The accomplishment of these tasks will increase our understanding of the dynamics of wake vortex and improve forecasting systems responsible for air safety and efficiency. The first two tasks include following three parts: (a) Determine significant length scale for vortex decay and transport, especially the length scales associated with the onset of Crow instability (Crow, 1970); (b) Study the effects of atmospheric turbulence on the decay of the wake vortices; and (c) Determine the relationships between decay rate, transport properties and atmospheric parameters based on large eddy simulation (LES) results and the observational data. These parameters may include turbulence kinetic energy, dissipation rate, wind shear and atmospheric stratification. The ABL studies cover LES modeling of turbulence structure within planetary boundary layer under transition and stable stratification conditions. Evidences have shown that the turbulence in the stable boundary layer can be highly intermittent and the length scales of eddies are very small compared to those in convective case. We proposed to develop a nesting grid mesh scheme and a modified Klemp-Durran conditions (Klemp and Wilhelmson, 1978) at the top boundary for TASS model to simulate planetary boundary layer under stable stratification conditions. During the past year, our group has made great efforts to carry out the above mentioned four tasks simultaneously. The work accomplished in the last year will be described in the next section.

Lin, Yuh-Lang↗

Energy Spectra of Strongly Stratified and Rotating Turbulence

Turbulence under strong stratification and rotation is usually characterized as quasi-two dimensional turbulence. We develop a "quasi-two dimensional" energy spectrum which changes smoothly between the Kolmogorov -5/3 law (no stratification), the -2 scalings of Zhou for the case of strong rotation, as well as the -2 scalings for the case of strong rotation and stratification. For strongly stratified turbulence, the model may give the -2 scaling predicted by Herring; and the -5/3 scaling indicated by some mesoscale observations.

Mahalov, Alex↗

Storm-driven Mixing and Potential Impact on the Arctic Ocean

Observations of the ocean, atmosphere, and ice made by Ice-Ocean Environmental Buoys (IOEBs) indicate that mixing events reaching the depth of the halocline have occurred in various regions in the Arctic Ocean. Our analysis suggests that these mixing events were mechanically forced by intense storms moving across the buoy sites. In this study, we analyzed these mixing events in the context of storm developments that occurred in the Beaufort Sea and in the general area just north of Fram Strait, two areas with quite different hydrographic structures. The Beaufort Sea is strongly influenced by inflow of Pacific water through Bering Strait, while the area north of Fram Strait is directly affected by the inflow of warm and salty North Atlantic water. Our analyses of the basin-wide evolution of the surface pressure and geostrophic wind fields indicate that the characteristics of the storms could be very different. The buoy-observed mixing occurred only in the spring and winter seasons when the stratification was relatively weak. This indicates the importance of stratification, although the mixing itself was mechanically driven. We also analyze the distribution of storms, both the long-term climatology as well as the patterns for each year in the last two decades. The frequency of storms is also shown to be correlated- (but not strongly) to Arctic Oscillation indices. This study indicates that the formation of new ice that leads to brine rejection is unlikely the mechanism that results in the type of mixing that could overturn the halocline. On the other hand, synoptic-scale storms can force mixing deep enough to the halocline and thermocline layer. Despite a very stable stratification associated with the Arctic halocline, the warm subsurface thermocline water is not always insulated from the mixed layer.

Yang, Jiayan↗

Large Eddy Simulations of Upper-Ocean Turbulence Across the Equatorial Pacific Cold Tongue

Turbulence in the equatorial Pacific cold tongue plays multiple roles in the sea-surface temperature budget. However, observations of turbulence are limited to a handful of field experiments lasting from a few days to one month at 0 N , 140 W, and the key turbulent fluxes can only be inferred. Although explicit simulations of turbulence resolve these fluxes, the required large eddy simulations (LES) are costly and thus also limited to a few weeks. A two-fold hypothesis motivates this contribution to the PUMP pre-field modelling experiments in support of TPOS: 1) that there exist useful relationships between the surface wind, wave, and buoyancy forcing, oceanic shear and stratification, and turbulent mixing, but 2) these relationships differ between regions of equatorial dynamics (f~0) and higher latitudes, such that the state space is insufficiently sampled to understand these relationships and represent them in ocean general circulation models (OGCM). To address the lack of observations, we simulated 35 continuous days of upper ocean turbulence using LES during October-November 1985 at sites in two regions (70 days in total) of the cold tongue: 0 N, 140 W and 3 N, 140 W. The model reproduces many features of the observed turbulence at 0,140W, including the “deep-cycle” observed by Gregg et al. (1985) Nature and the turbulent pulses every few hours during the night and in the early morning. Furthermore, the results at 3 N support the hypothesis advanced by Cherian et al. (2021) JPO that strong, but intermittent deep-cycle turbulence is observed to the north of the equator in the presence of strong shear associated with Tropical Instability Waves. Comparisons with the K-profile parameterization (KPP) scheme show that additional tuning of the KPP parameters can reduce differences between the turbulent fluxes produced by KPP and LES at both sites. The results will suggest modifications to OGCMs and already emphasize the need for a more extensive measurements of vertical profiles of shear, stratification and turbulent mixing to reveal the spatio-temporal structure of turbulence and how it relates to forcing, shear and stratification across the equatorial Pacific Cold tongue.

Large Eddy Simulations↗

Comparison of plug flow and multi-node stratified tank modeling approaches regarding computational efficiency and accuracy

Residential water heaters contain water stratified by temperature-driven density differences. This implies that a water tank can reach a state in which the top and bottom sections have different temperatures, unless mixing happens. A high degree of thermal stratification can improve the efficiency of some water heaters, by saving the amount of energy required for the heat-up process. Studies of stratification became popular in the 1970s and it remains an active research topic today. The research has led to the development of different models and techniques to better predict and define a stratified tanks behavior. By comparing these models and techniques used previously to describe thermal stratification, the phenomenon could be better understood, exploited, and used to increase efficiency and thermal energy capacity in modern water tanks. From the existing models, we found the one-dimensional standard plug-flow and a multi node model to be appropriate for analyzing the processes of the heat up and cool-down in a water tank. These two models are based on energy balances. This work involved comparing the accuracy and computational effort needed to implement these models. To assess accuracy, we compared both types of existing models to experimental data (also collected in this work) which included a heat up process using an external heat pump. This external process included a layering process that has an eddy diffusivity at five times the rate of thermal diffusion. For this project, we implemented the models in MATLAB, the multi-paradigm numerical computing environment. We quantified model accuracy using the root mean squared error between modeled data and experimental data for six measured tank temperatures. Comparing the accuracy and the computational time taken to run the simulation provides a method to contrast the performance of each model and a way to rate it. The multi node model was run using from 6 to 96 spatial nodes; the plug flow model was run using 1 to 0.001 º C temperature bin sizes. Additionally, timesteps were varied from 4 to 236 s. The results quantify the tradeoff between accuracy and computational time, providing guidance for simulations to intelligently select the best model type and simulation parameters. This research can be used to validate the pre-existing models and possibly improve the modern water tank.

Bulnes, Fernando Karg↗

Two-zone stratified wetwell model development and implementation for RELAP-7

The Reactor Core Isolation Cooling (RCIC) system consists of a turbine, pump and wetwell and is used in many US Boiling Water Reactors (BWRs) as an important heat removal system. The wetwell is a major heat sink within containment. Experimental investigations have shown that thermal stratification may occur in the wetwell during RCIC operation. Current systems-level analysis codes lack stratified wetwell models. Here, this paper describes the development of a stratified wetwell model for the RELAP-7 code. This paper uses a simple two-zone model for the wetwell water that can capture thermal stratification. This model assumes that a buoyant plume transports heat and mass from the steam injection site to the upper water zone. The buoyant plume assumption requires that steam is injected at a low enough rate so it adds negligible momentum. Model validation against experimental data shows similar trends of temperature field development in the suppression pool.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A hierarchical stochastic modeling approach for representing point bar geometries and petrophysical property variations

The flow of fluids in point bars is affected by the existence of heterogeneities like shale drapes that are found on the surfaces of inclined heterolithic stratifications. In fact, these shale drapes can act as fluid flow baffles; therefore, developing a framework for modeling point bars and their associated heterogeneities is vital. In this study, a stochastic process-based modeling approach is presented for capturing the main point bar heterogeneities: accretion surfaces (i.e., the aerial heterogeneity) and inclined heterolithic stratifications (i.e., the vertical heterogeneity). The former was modeled using a sine-generation function and the latter, with a sigmoidal function after which they were combined into a 3D point bar model. To ensure proper modeling of petrophysical properties, we developed a more representative gridding scheme which generates curvilinear grids representative of the point bar geometry. This grid was then transformed into a rectilinear grid to allow for geostatistical simulation after which all petrophysical properties were mapped back into the original curvilinear grid. An essential element of this modeling approach is the stochastic representation of shale drapes at the interface between successive accretion surfaces. The workflow was tested using a real field dataset for the Cranfield field, Mississippi. The constructed model was then subjected to a flow simulation study mimicking a CO 2 storage scenario. Various sensitivities were simulated to evaluate the effect of heterogeneities on CO 2 flow within the point-bar. Results demonstrate the importance of representing point-bar related heterogeneity and the spatial distribution of shale drapes on CO 2 plume migration and storage.

58 GEOSCIENCES↗

Computationally efficient evaluation of optimum homogeneous charge compression ignition operating range with accelerated multizone engine cycle simulation

The very intensive calculations necessary to define a performance map requiring evaluation of over a hundred individual operating points can be efficiently conducted with accelerated multizone for engine cycle simulation, leading to a definition of regions of acceptable and optimum homogeneous charge compression ignition operation. Accelerated multizone for engine cycle simulation has the virtue of enabling accurate evaluation of many operating conditions based on thermal stratification data from a single fluid mechanics run at motored conditions. This is possible because thermal stratification is more sensitive to engine geometry than to operating conditions. In this article, accuracy of accelerated multizone for engine cycle simulation is demonstrated by comparison with experimental data for iso-octane homogeneous charge compression ignition operation over a broad range of lean equivalence ratios (0.14–0.28). The validated accelerated multizone for engine cycle simulation model is then applied to generating a performance map for an engine controlled by appropriately adjusting equivalence ratio and internal exhaust gas recirculation. Regions of acceptable and optimum combustion are identified. It is finally demonstrated that while indicated mean effective pressure remains low for optimum homogeneous charge compression ignition operation (1–4 bar), this is sufficient for a large fraction of typical driving in light-duty vehicles. Much driving including idle can therefore be done in homogeneous charge compression ignition mode at high efficiency and low (essentially zero) NOx and particulate matter emissions.

30 DIRECT ENERGY CONVERSION↗

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗

A NEW MODEL FOR MIXING BY DOUBLE-DIFFUSIVE CONVECTION (SEMI-CONVECTION). III. THERMAL AND COMPOSITIONAL TRANSPORT THROUGH NON-LAYERED ODDC

Oscillatory double-diffusive convection (ODDC; also known as semi-convection) refers to a type of double-diffusive instability that occurs in regions of planetary and stellar interiors that have a destabilizing thermal stratification and a stabilizing mean molecular weight stratification. In this series of papers, we use an extensive suite of three-dimensional (3D) numerical simulations to quantify the transport of heat and chemical species by ODDC. Rosenblum et al. first showed that ODDC can either spontaneously form layers that significantly enhance the transport of heat and chemical species compared to microscopic transport or remain in a state dominated by large-scale gravity waves, in which there is a more modest enhancement of the turbulent transport rates. Subsequent studies in this series focused on identifying under what conditions layers form and quantifying transport through layered systems. Here we proceed to characterize transport through systems that are unstable to ODDC, but do not undergo spontaneous layer formation. We measure the thermal and compositional fluxes in non-layered ODDC from both two-dimensional (2D) and 3D numerical simulations, and show that 3D simulations are well approximated by similar simulations in a 2D domain. We find that the turbulent mixing rate in this regime is weak and can, to a first-level approximation, be neglected. We conclude by summarizing the findings of papers I through III into a single prescription for transport systems unstable to ODDC.

79 ASTRONOMY AND ASTROPHYSICS↗

TOWARD HIGH-PRECISION SEISMIC STUDIES OF WHITE DWARF STARS: PARAMETRIZATION OF THE CORE AND TESTS OF ACCURACY

We present a prescription for parametrizing the chemical profile in the core of white dwarfs in light of the recent discovery that pulsation modes may sometimes be deeply confined in some cool pulsating white dwarfs. Such modes may be used as unique probes of the complicated chemical stratification that results from several processes that occurred in previous evolutionary phases of intermediate-mass stars. This effort is part of our ongoing quest for more credible and realistic seismic models of white dwarfs using static, parametrized equilibrium structures. Inspired by successful techniques developed in design optimization fields (such as aerodynamics), we exploit Akima splines for the tracing of the chemical profile of oxygen (carbon) in the core of a white dwarf model. A series of tests are then presented to better seize the precision and significance of the results that can be obtained in an asteroseismological context. We also show that the new parametrization passes an essential basic test, as it successfully reproduces the chemical stratification of a full evolutionary model.

79 ASTRONOMY AND ASTROPHYSICS↗

Reverberation Mapping of Optical Emission Lines in Five Active Galaxies

We present the first results from an optical reverberation mapping campaign executed in 2014 targeting the active galactic nuclei (AGNs) MCG+08-11-011, NGC 2617, NGC 4051, 3C 382, and Mrk 374. Our targets have diverse and interesting observational properties, including a “changing look” AGN and a broad-line radio galaxy. Based on continuum-H β lags, we measure black hole masses for all five targets. We also obtain H γ and He ii λ 4686 lags for all objects except 3C 382. The He ii λ 4686 lags indicate radial stratification of the BLR, and the masses derived from different emission lines are in general agreement. The relative responsivities of these lines are also in qualitative agreement with photoionization models. These spectra have extremely high signal-to-noise ratios (100–300 per pixel) and there are excellent prospects for obtaining velocity-resolved reverberation signatures.

79 ASTRONOMY AND ASTROPHYSICS↗