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

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↗
DOE OSTI2020

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↗
DOE OSTI2020

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

Positive Configuration Space

Abstract We define and study the totally nonnegative part of the Chow quotient of the Grassmannian, or more simply the nonnegative configuration space . This space has a natural stratification by positive Chow cells , and we show that nonnegative configuration space is homeomorphic to a polytope as a stratified space. We establish bijections between positive Chow cells and the following sets: (a) regular subdivisions of the hypersimplex into positroid polytopes, (b) the set of cones in the positive tropical Grassmannian, and (c) the set of cones in the positive Dressian. Our work is motivated by connections to super Yang–Mills scattering amplitudes, which will be discussed in a sequel.

Physics↗
DOE OSTI2022

Degradation Reduces Microbial Richness and Alters Microbial Functions in an Australian Peatland

Abstract Peatland ecosystems cover only 3% of the world’s land area; however, they store one-third of the global soil carbon (C). Microbial communities are the main drivers of C decomposition in peatlands, yet we have limited knowledge of their structure and function. While the microbial communities in the Northern Hemisphere peatlands are well documented, we have limited understanding of microbial community composition and function in the Southern Hemisphere peatlands, especially in Australia. We investigated the vertical stratification of prokaryote and fungal communities from Wellington Plains peatland in the Australian Alps. Within the peatland complex, bog peat was sampled from the intact peatland and dried peat from the degraded peatland along a vertical soil depth gradient (i.e., acrotelm, mesotelm, and catotelm). We analyzed the prokaryote and fungal community structure, predicted functional profiles of prokaryotes using PICRUSt, and assigned soil fungal guilds using FUNGuild. We found that the structure and function of prokaryotes were vertically stratified in the intact bog. Soil carbon, manganese, nitrogen, lead, and sodium content best explained the prokaryote composition. Prokaryote richness was significantly higher in the intact bog acrotelm compared to degraded bog acrotelm. Fungal composition remained similar across the soil depth gradient; however, there was a considerable increase in saprotroph abundance and decrease in endophyte abundance along the vertical soil depth gradient. The abundance of saprotrophs and plant pathogens was two-fold higher in the degraded bog acrotelm. Soil manganese and nitrogen content, electrical conductivity, and water table level (cm) best explained the fungal composition. Our results demonstrate that both fungal and prokaryote communities are shaped by soil abiotic factors and that peatland degradation reduces microbial richness and alters microbial functions. Thus, current and future changes to the environmental conditions in these peatlands may lead to altered microbial community structures and associated functions which may have implications for broader ecosystem function changes in peatlands.

Environmental Sciences & Ecology↗
DOE OSTI2023

Another Year of Record Heat for the Oceans

Abstract Changes in ocean heat content (OHC), salinity, and stratification provide critical indicators for changes in Earth’s energy and water cycles. These cycles have been profoundly altered due to the emission of greenhouse gasses and other anthropogenic substances by human activities, driving pervasive changes in Earth’s climate system. In 2022, the world’s oceans, as given by OHC, were again the hottest in the historical record and exceeded the previous 2021 record maximum. According to IAP/CAS data, the 0–2000 m OHC in 2022 exceeded that of 2021 by 10.9 ± 8.3 ZJ (1 Zetta Joules = 10 21 Joules); and according to NCEI/NOAA data, by 9.1 ± 8.7 ZJ. Among seven regions, four basins (the North Pacific, North Atlantic, the Mediterranean Sea, and southern oceans) recorded their highest OHC since the 1950s. The salinity-contrast index, a quantification of the “salty gets saltier—fresh gets fresher” pattern, also reached its highest level on record in 2022, implying continued amplification of the global hydrological cycle. Regional OHC and salinity changes in 2022 were dominated by a strong La Niña event. Global upper-ocean stratification continued its increasing trend and was among the top seven in 2022.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Drivers of coupled climate model biases in representing Labrador Sea convection

Abstract This study investigates the representation of ocean convection in the Labrador Sea in seven Earth System Models (ESMs) from the Coupled Model Intercomparison Project Phase 5 and 6 datasets. The relative role of the oceanic and atmospheric biases in the subpolar North Atlantic gyre are explored using regional ocean simulations where the atmospheric forcing or the ocean initial and boundary conditions are replaced by reanalysis data in the absence of interactive air-sea coupling. Commonalities and differences among model behaviors are discussed with the objective of finding a pathway forward to improve the representation of the ocean mean state and variability in a region of fundamental importance for climate variability and change. Results highlight that an improved representation of ocean stratification in the North Atlantic subpolar gyre is urgently needed to constrain future climate change projections. While improving the ocean model resolution in the North Atlantic alone may contribute a better representation of both boundary currents and propagation of heat and freshwater anomalies into the Labrador Sea, it may not be sufficient. Addressing the atmospheric heat flux bias with better resolution in the atmosphere and land topography may allow for deep convection to occur in the Labrador Sea in some of the models that miss it entirely, but the greatest priority remains improving the representation of ocean stratification.

Liu, Guangpeng↗
DOE OSTI2026

Climate-eutrophication-anoxia interactions in Late Glacial Soppensee, Switzerland: Forcings, non-linear responses and recovery

Combined effects of climate warming and anthropogenic nutrient loadings lead to lake eutrophication and anoxia globally. Because of chemical feedbacks, lakes under multiple stressors often respond in non-linear ways. However, it remains unclear whether climate change alone can lead to non-linear lake responses in the absence of anthropogenic nutrient disturbances. Here, we investigate the interactions between climate variability, nutrient cycling and trophic state changes, mixing regimes, anoxia and related chemical feedback in a small kettle-hole lake in Switzerland during Late Glacial times (15.2–12.6 cal ka BP), a period known for high-amplitude climate change in pre-anthropogenic times. After its formation during Heinrich Stadial 1 (>15 cal ka BP), Soppensee was oligotrophic and well-mixed. Soppensee became eutrophic and developed anoxia at 14.25 cal ka BP. Phosphorus (P) was released from sediments through the reductive dissolution of Fe-oxyhydroxides, fuelling eutrophication. Eutrophication lagged the Bølling warming (14.65 cal ka BP) by 400 years, suggesting that rising temperatures were not the trigger for eutrophication. Instead, eutrophication responded non-linearly to forest closure (threshold at 76 % arboreal pollen AP), which shielded Soppensee from wind mixing, enhancing lake stratification, anoxia and P release, intensifying eutrophication. These conditions ended during the 200-years cold period of the Aegelsee Oscillation (GI-1d, ca. 14.0 cal ka BP) when the landscape regionally opened (AP<76 %); the lake became well-mixed, oxygenated and P was efficiently sequestered. Throughout the Allerød (13.9–12.8 cal ka BP), enhanced Fe input prompted diagenetic vivianite formation, sequestering P in sediments, naturally remediating lake eutrophication despite closed forests, warm temperatures, lake stratification and anoxia.

Environmental sciences↗
DOE OSTI2022

Aboveground biomass density models for NASA’s Global Ecosystem Dynamics Investigation (GEDI) lidar mission

NASA's Global Ecosystem Dynamics Investigation (GEDI) is collecting spaceborne full waveform lidar data with a primary science goal of producing accurate estimates of forest aboveground biomass density (AGBD). This paper presents the development of the models used to create GEDI's footprint-level (~25 m) AGBD (GEDI04_A) product, including a description of the datasets used and the procedure for final model selection. The data used to fit our models are from a compilation of globally distributed spatially and temporally coincident field and airborne lidar datasets, whereby we simulated GEDI-like waveforms from airborne lidar to build a calibration database. We used this database to expand the geographic extent of past waveform lidar studies, and divided the globe into four broad strata by Plant Functional Type (PFT) and six geographic regions. GEDI's waveform-to-biomass models take the form of parametric Ordinary Least Squares (OLS) models with simulated Relative Height (RH) metrics as predictor variables. From an exhaustive set of candidate models, we selected the best input predictor variables, and data transformations for each geographic stratum in the GEDI domain to produce a set of comprehensive predictive footprint-level models. We found that model selection frequently favored combinations of RH metrics at the 98th, 90th, 50th, and 10th height above ground-level percentiles (RH98, RH90, RH50, and RH10, respectively), but that inclusion of lower RH metrics (e.g. RH10) did not markedly improve model performance. Second, forced inclusion of RH98 in all models was important and did not degrade model performance, and the best performing models were parsimonious, typically having only 1-3 predictors. Third, stratification by geographic domain (PFT, geographic region) improved model performance in comparison to global models without stratification. Fourth, for the vast majority of strata, the best performing models were fit using square root transformation of field AGBD and/or height metrics. There was considerable variability in model performance across geographic strata, and areas with sparse training data and/or high AGBD values had the poorest performance. These models are used to produce global predictions of AGBD, but will be improved in the future as more and better training data become available.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Nondimensional convection numbers modeling thermally stratified storage tanks: Richardson's number and hot-water tanks

Thermally stratified storage tank studies have spanned over 50 years to increase the thermal storage efficiency and accurate prediction of the outlet temperature particularly for solar applications. The studies have reviewed and modeled the jet and plume flow phenomena inside the tank due to the inlet mixing and stratification level. Kelvin–Helmholtz and Rayleigh–Taylor instabilities are the major drivers of the mixing in these tanks. Momentum jets deflecting off walls at the bottom of the tank also create significant mixing. Reviewing Richardson models shows that the categorization was based on the range of Reynolds numbers at the inlet. Unfortunately, the use of superficial velocity in calculating the Richardson number results in critical values in the literature ranging from below 0.25 to 100. The most used length scale associated with these flows is an inertial scale based on the tank height or diameter although the mixing can occur at a relatively smaller scale. The various inlet devices and a large span of flow rates experienced in thermally stratified storage tanks requisite the use of the Reynolds number in combination with a convection number for accurate one dimensional models that predict performance over the long-term. The evaluation of peak shifting of electric loads leveraging renewable sources for applications, including residential hot-water tanks, commercial water tanks, and large-scale chilled water storage tanks, require these models. Finally, this paper is focused on establishing the significance of the convection numbers in conjunction with the Reynolds number for modeling the thermal stratification in storage tanks.

14 SOLAR ENERGY↗
DOE OSTI2023

Uniform momentum and temperature zones in unstably stratified turbulent flows

Wall-bounded turbulent flows exhibit a zonal arrangement, in which streamwise velocity organizes into uniform momentum zones (UMZs), separated by thin layers of elevated interfacial shear. While significant research efforts have focused on these structural features in neutrally stratified flows, the effects of unstable thermal stratification on UMZs and on analogous uniform temperature zones (UTZs) have not been considered previously. In this article, statistical properties of UMZs and UTZs are investigated using a suite of large eddy simulations of unstably stratified turbulent channel flow spanning weakly to highly convective conditions. When normalized by the friction velocity and stability-dependent mixing length, the mean velocity gradient based on UMZ interfacial velocity jumps and the vorticity thickness exhibits good collapse for all stabilities, establishing a link between UMZ properties and scaling predictions from Monin–Obukhov similarity theory. A similar relationship is found between UTZ properties and surface-layer scaling of the mean temperature gradient. In the mixed layer, mean UMZ depth is quasi-constant with wall-normal distance, while the deepest UTZs are found in the centre of the boundary layer. These instantaneous structures are found to be linked to the well-mixed velocity and temperature profiles in the convective mixed layer. Conditional averaging indicates that both UMZ and UTZ interfaces are associated with ejections of momentum and warm updrafts below the interface and sweeps of momentum and cool downdrafts above the interface. These results demonstrate a tangible connection between instantaneous structural features, mean properties and scaling laws in unstably stratified flows.

Mechanics↗
DOE OSTI2025

Oceanic mixing and waves in the presence of a suspended canopy

Large-eddy simulations are analysed to determine the influence of suspended canopies, such as those formed in macroalgal farms, on ocean mixed layer (OML) deepening and internal wave generation. In the absence of a canopy, we show that Langmuir turbulence, when compared with wind-driven shear turbulence, results in a deeper OML and more pronounced internal waves beneath the OML. Subsequently, we examine simulations with suspended canopies of varying densities located in the OML, in the presence of a background geostrophic current. Intensified turbulence occurs in the shear layer at the canopy’s bottom edge, arising from the interaction between the geostrophic current and canopy drag. Structures resembling Kelvin–Helmholtz (KH) instability emerge as the canopy shear layer interacts with the underlying stratification, radiating internal waves beneath the OML. Both intensified turbulence and lower-frequency motions associated with KH-type structures are critical factors in enhancing mixing. Consequently, the OML depth increases by up to a factor of two compared with cases without a canopy. Denser canopies and stronger geostrophic currents lead to more pronounced KH-type structures and internal waves, stronger turbulence and greater OML deepening. Additionally, vertical nutrient transport is enhanced as the OML deepens due to the presence of the canopy. Considering that the canopy density investigated in this study closely represents offshore macroalgal farms, these findings suggest a mechanism for passive nutrient entrainment conducive to sustainable farming. Overall, this study reveals the intricate interactions between the suspended canopy, turbulent mixing and stratification, underscoring their importance in reshaping OML characteristics.

Bo, Tong (ORCID:0000000260300561)↗
DOE OSTI2020

Baroclinic Residual Circulation and Mass Transport Due to Internal Tides

Baroclinic (BC) tidal residual circulation due to internal tides is investigated around islands over a shallow ridge using a numerical ocean model. Internal tides enhance vertical mixing over shallow slopes, leading to horizontal density gradients that drive BC residual circulation along the main thermocline. For a strongly stratified summer case, the vertical diffusivity estimated by the Mellor and Yamada turbulence closure model exceeds 1 × 10 -2 m 2 s -1 , and the velocity of BC residual circulations reaches 0.2 m s -1 . The magnitude of BC residual circulation is larger than that of barotropic residual tidal circulation, implying that BC residual circulation due to internal tides plays an important role in coastal ocean circulation. Furthermore, BC residual flow accounts for an equal percentage (5%) of the total tidal kinetic energy as the barotropic residual flow under summer stratification conditions. Results from a coupled sediment resuspension and transport model show the growth of intermediate nepheloid layers formed by strong bottom shear stress and BC residual flow. The residual component contributes largely (23% of the total in summer) to the total suspended sediment flux. Seasonal variability is explored, with weaker winter stratification leading to reduced mixing, and thus weaker BC residual circulation and sediment flux. Lastly, the magnitude of the BC residual circulation is also shown to be proportional to the square of the tidal amplitude.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Towards a self-consistent model of the convective core boundary in upper main sequence stars: I. 2.5D and 3D simulations

There is strong observational evidence that the convective cores of intermediate-mass and massive main sequence stars are substantially larger than those predicted by standard stellar-evolution models. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-timescale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number S EVEN -L EAGUE H YDRO code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a 15 M ⊙ zero-age main sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor, ranging from 10 3 to 10 6 . After reaching thermal equilibrium, all of our models show a clear penetration layer; its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. The penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification in simulations at the lower end of our luminosity range. This structure corresponds to the adiabatic ‘step overshoot’ model often employed in stellar-evolution calculations. The simulations with the highest and lowest luminosity differ by less than a factor of two in the penetration distance. The high computational cost of 3D simulations makes our current 3D data set rather sparse. Depending on how we extrapolate the 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from 0.09 to 0.44 pressure scale heights, which is broadly compatible with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Observed declines in upper ocean phosphate-to-nitrate availability

Climate warming is increasing ocean stratification, which in turn should decrease the nutrient flux to the upper ocean. This may slow marine primary productivity, causing cascading effects throughout food webs. However, observing changes in upper ocean nutrients is challenging because surface concentrations are often below detection limits. We show that the nutricline depth, where nutrient concentrations reach well-detected levels, is tied to productivity and upper ocean nutrient availability. Next, we quantify nutricline depths from a global database of observed vertical nitrate and phosphate profiles to assess contemporary trends in global nutrient availability (1972–2022). We find strong evidence that the P-nutricline (phosphacline) is mostly deepening, especially throughout the southern hemisphere, but the N-nutricline (nitracline) remains mostly stable. Earth System Model (ESM) simulations support the hypothesis that reduced iron stress and increased nitrogen fixation buffer the nitracline, but not phosphacline, against increasing stratification. These contemporary trends are expected to continue in the coming decades, leading to increasing phosphorus but not nitrogen stress for marine phytoplankton, with important ramifications for ocean biogeochemistry and food web dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Striations in moderate pressure dc driven nitrogen glow discharge

Plasma stratification has been studied for more than a century. Despite the many experimental studies reported on this topic, theoretical analyses and numerical modeling of this phenomenon have been mostly limited to rare gases. In this work, a one-dimensional fluid model with detailed kinetics of electrons and vibrationally excited molecules is employed to simulate moderate-pressure (i.e. a few Torrs) dc discharge in nitrogen in a 15.5 cm long tube of radius 0.55 cm. The model also considers ambipolar diffusion to account for the radial loss of ions and electrons to the wall. The proposed model predicts self-excited standing striations in nitrogen for a range of discharge currents. The impact of electron transport parameters and reaction rates obtained from a solution of local two-term and a multi-term Boltzmann equation on the predictions are assessed. In-depth kinetic analysis indicates that the striations result from the undulations in electron temperature caused due to the interaction between ionization and vibrational reactions. Furthermore, the vibrationally excited molecules associated with the lower energy levels are found to influence nitrogen plasma stratification and the striation pattern strongly. A balance between ionization processes and electron energy transport allows the formation of the observed standing striations. Simulations were conducted for a range of discharge current densities from ~0.018 to 0.080 mA cm –2 , for an operating pressure of 0.7 Torr. Parametric studies show that the striation length decreases with increasing discharge current. Furthermore, the predictions from the model are compared against experimental measurements and are found to agree favorably.

42 ENGINEERING↗
DOE OSTI2022

The principle of minimal power

This article is devoted to the memory of Yuri P Raizer, who passed away in 2021. He left a noticeable trace in gas discharge physics. The principle of minimal power (the state that requires minimal power is most probable) is thoroughly used in his books. Although the fundamental laws of physics do not imply this ad hoc principle, a detailed analysis of underlying phenomena can often reveal why nature prefers this path. Raizer illustrated this principle for plasma stratification, formation of electrode spots, discharge constriction, the shape of an arc channel, etc. We argue that the nonlinearity of equations describing gas discharges can often justify the realization of a plasma state maintained at minimal electric power. This nonlinearity appears because small groups of energetic electrons often control the ionization processes. The number of these electrons depends strongly on the ratio of the electric field to gas density, E/N. Under certain conditions, the ionization rate can also depend nonlinearly on electron density due to stepwise ionization and Coulomb collisions. Here, we use the principle of minimal power to illustrate some of Raizer's contributions to gas discharge physics from a single point of view. We demonstrate that nonlinearity of ionization processes in gas discharges can substantiate this principle for plasma stratification. However, striations of s, p, and r types in neon could exist with minimal or no ionization enhancement. This reminds us of Raizer's warning that applying the minimal power principle could lead to erroneous predictions, and a proper theory is required in each case to justify its use. 'The phenomenon of striations satisfies the principle of minimal power' – Yuri Raizer

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Large-eddy simulations of the Northeastern US coastal marine boundary layer

In this study, large eddy simulations (LES) of offshore boundary layers near the Nantucket coast are performed using Nalu-Wind. The marine boundary layer conditions are chosen to match the predominantly unstable and neutral conditions measured by the Cape Wind platform. The appropriate domain, resolution, and boundary condition settings required for the LES are established through this work. Differences between stable and unstable cases are found in the wind speed profiles, averaged statistics, and wind spectra, and explained in terms of stratification effects.

54 ENVIRONMENTAL SCIENCES↗