Search NASA⌕ Search

SEARCH · Search NASA

Results for “eddying model”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Three-step LES-C models for flows at high Reynolds numbers

Here, we investigate the need for the second correction step in the recently proposed LES-C models for fluid flows at high Reynolds numbers. These models use a predictor-corrector idea to enhance the efficiency of the existing Large Eddy Simulation models. Different three-step(one defect step, two corrections) LES-C models, based on the Leray-α, ADM and NS-ω LES models, are tested in three different situations. The new Leray-α-C2 model (C2 stands for two correction steps) is applied to the Navier–Stokes equations; the ADC2 is applied to the MagnetoHydroDynamic flow; and the NS-ω-C2 is used in the fluid-fluid interaction problem. We evaluate the effectiveness of the second correction step in all these settings, using qualitative and quantitative numerical tests.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Idealized Large-Eddy Simulations of Stratocumulus Advecting over Cold Water. Part I: Boundary Layer Decoupling

Abstract We explore the decoupling physics of a stratocumulus-topped boundary layer (STBL) moving over cooler water, a situation mimicking warm-air advection (WADV). We simulate an initially well-mixed STBL over a doubly periodic domain with the sea surface temperature decreasing linearly over time using the System for Atmospheric Modeling large-eddy model. Due to the surface cooling, the STBL becomes increasingly stably stratified, manifested as a near-surface temperature inversion topped by a well-mixed cloud-containing layer. Unlike the stably stratified STBL in cold-air advection (CADV) that is characterized by cumulus coupling, the stratocumulus deck in the WADV is unambiguously decoupled from the sea surface, manifested as weakly negative buoyancy flux throughout the subcloud layer. Without the influxes of buoyancy from the surface, the convective circulation in the well-mixed cloud-containing layer is driven by cloud-top radiative cooling. In such a regime, the downdrafts propel the circulation, in contrast to that in CADV regime for which the cumulus updrafts play a more determinant role. Such a contrast in convection regime explains the difference in many aspects of the STBLs including the entrainment rate, cloud homogeneity, vertical exchanges of heat and moisture, and lifetime of the stratocumulus deck, with the last being subject to a more thorough investigation in Part II. Finally, we investigate under what conditions a secondary stratus near the surface (or fog) can form in the WADV. We found that weaker subsidence favors the formation of fog whereas a more rapid surface cooling rate does not. Significant Statement The low-lying blanketlike clouds, called stratocumulus (Sc), reflect much incoming sunlight, substantially modulating Earth’s temperature. While much is known about how the Sc evolves when it moves over warmer water, few studies examine the opposite situation of Sc moving over colder water. We used a high-resolution numerical model to simulate such a case. When moving over cold water, the Sc becomes unambiguously decoupled from the water surface, distinctive from its warm counterpart in which the Sc interacts with the water surface via intermittent cauliflower-like clouds called cumulus clouds. Such decoupling influences many aspects of the Sc–sea surface system, which combine to alter the ability of the Sc to reflect sunlight, thereby influencing the climate. This work laid the foundation for future work that quantifies the contribution of such a decoupled Sc regime to Earth’s radiative budget and climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Flow Rate Dependent 1D Model for Thermally Stratified Hot-Water Energy Storage

Stratified tank models are used to simulate thermal storage in applications such as residential or commercial hot-water storage tanks, chilled-water storage tanks, and solar thermal systems. The energy efficiency of these applications relates to the system components and the level of stratification maintained during various flow events in the tank. One-dimensional (1D) models are used in building energy simulations because of the short computation time but often do not include flow-rate dependent mixing. The accuracy of 1D models for plug flow, plug flow with axial conduction, and two convection eddy-diffusivity models were compared with experimental data sets for discharging a 50-gal residential tank and recharging the tank with hot water from an external hot-water source. A minimum and maximum relationship for the eddy diffusivity factor were found at Re <2100 and >10,000 for recirculation of hot water to the top of the tank and vertical tubes inletting cold water at the bottom. The root mean square error decreased from >4 °C to near 2 °C when considering flow-based mixing models during heating, while the exponential decay of the eddy diffusion results in a root mean square error reduction of 1 °C for cone-shaped diffusers that begin to relaminarize flow at the inlet.

1D flow model↗
DOE OSTI2025

A New Framework for the Attribution of Air‐Sea CO 2 Exchange

Abstract The air‐sea transfer of carbon dioxide can be viewed as a dynamic system through which atmospheric and oceanic processes push surface waters away from thermodynamic equilibrium, while diffusive gas transfer pulls them back toward local equilibrium. These push/pull processes drive significant sub‐seasonal, seasonal, and interannual variability in air‐sea carbon fluxes, the quantification of which is critical both for diagnosing the ocean response to fossil fuel emissions and for attempts to mitigate anthropogenic climate disruption through intentional modification of surface ocean biogeochemistry. In this study, we present a new approach for attributing air‐sea carbon fluxes to specific mechanisms. The new framework is first applied to a two‐box ocean nutrient and carbon cycle model as an illustrative example. Next, outputs from a regional eddy‐resolving model of the Southern Ocean are analyzed. The roles of multiple physical and biogeochemical processes are identified. The decomposition of the seasonal air‐sea carbon flux shows the dominant role of biological carbon pumps that are partially compensated by the transport convergence. Finally, the framework is used to diagnose the response to mesoscale iron and alkalinity release, explicitly quantifying transport feedback and eventual impacts on net air‐sea carbon flux. Ocean carbon transport has divergent influences between iron and alkalinity release, due to opposing near‐surface gradients of dissolved inorganic carbon. We suggest that our attribution framework may be a useful analytical technique for monitoring natural ocean carbon fluxes and quantifying the impacts of human intervention on the ocean carbon cycle.

Ito, Takamitsu [School of Earth and Atmospheric Sc↗
DOE OSTI2023

A baseline for ensemble-based, time-resolved inflow reconstruction for a single turbine using large-eddy simulations and latent diffusion models

We are interested in reconstructing winds flowing through a turbine on a second-by-second basis over a 10 min window. Previously, we developed a machine learning algorithm that takes in a snapshot of wind speed measurements and generates ensembles of three-dimensional wind field estimates. Here, we use these estimates as initial conditions in large-eddy simulations and reconstruct atmospheric and turbine response dynamic quantities in a synthetic field campaign. In doing so, we establish a baseline for model validation that future time-aware data assimilation techniques will be compared to. In turbine-free case studies, ground truth wind speeds consistently fall within our estimated wind speed distribution for the first 100 s after the simulation start. In simulations with turbines, the wind estimates show a small bias of 0.10 m s 1 and good correlation of 0.80 during the first 100 s. During this window, our estimates of the Blade 1 bending moment and generator power typically span the ground truth, with the estimate of the former performing better overall. In summary, this approach shows promise as a stand-alone technique for reconstructing real-world inflow and turbine dynamics in 1-2 min windows and as a foundation for future time-aware data assimilation techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Upscaling Wetland Methane Emissions From the FLUXNET–CH4 Eddy Covariance Network (UpCH4 v1.0): Model Development, Network Assessment, and Budget Comparison

Wetlands are responsible for 20%–31% of global methane (CH 4 ) emissions and account for a large source of uncertainty in the global CH 4 budget. Data-driven upscaling of CH 4 fluxes from eddy covariance measurements can provide new and independent bottom-up estimates of wetland CH 4 emissions. Here, we develop a six-predictor random forest upscaling model (UpCH4), trained on 119 site-years of eddy covariance CH 4 flux data from 43 freshwater wetland sites in the FLUXNET-CH4 Community Product. Network patterns in site-level annual means and mean seasonal cycles of CH 4 fluxes were reproduced accurately in tundra, boreal, and temperate regions (Nash-Sutcliffe Efficiency ~0.52–0.63 and 0.53). UpCH4 estimated annual global wetland CH 4 emissions of 146 ± 43 TgCH 4 y –1 for 2001–2018 which agrees closely with current bottom-up land surface models (102–181 TgCH 4 y –1 ) and overlaps with top-down atmospheric inversion models (155–200 TgCH 4 y –1 ). However, UpCH4 diverged from both types of models in the spatial pattern and seasonal dynamics of tropical wetland emissions. We conclude that upscaling of eddy covariance CH 4 fluxes has the potential to produce realistic extra-tropical wetland CH 4 emissions estimates which will improve with more flux data. To reduce uncertainty in upscaled estimates, researchers could prioritize new wetland flux sites along humid-to-arid tropical climate gradients, from major rainforest basins (Congo, Amazon, and SE Asia), into monsoon (Bangladesh and India) and savannah regions (African Sahel) and be paired with improved knowledge of wetland extent seasonal dynamics in these regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modeling and simulation of transitional Rayleigh–Taylor flow with partially averaged Navier–Stokes equations

In this work, the partially averaged Navier–Stokes (PANS) equations are used to predict the variable-density Rayleigh–Taylor (RT) flow at Atwood number 0.5 and maximum Reynolds number 500. This is a prototypical problem of material mixing, featuring laminar, transitional, and turbulent flow, instabilities and coherent structures, density fluctuations, and production of turbulence kinetic energy by both shear and buoyancy mechanisms. These features pose numerous challenges to modeling and simulation, making the RT flow ideal to develop the validation space of the recently proposed PANS Besnard–Harlow–Rauenzahn-linear eddy viscosity model closure. The numerical simulations are conducted at different levels of physical resolution and test three approaches to set the parameters $f_\phi$ defining the range of physically resolved scales. The computations demonstrate the efficiency (accuracy vs cost) of the PANS model predicting the spatiotemporal development of the RT flow. Results comparable to large-eddy simulations and direct numerical simulations are obtained at significantly lower physical resolution without the limitations of the Reynolds-averaged Navier–Stokes equations in these transitional flows. The data also illustrate the importance of appropriate selection of the physical resolution and the resolved fraction of each dependent quantity $\phi$ of the turbulent closure, $f_\phi$. These two aspects determine the ability of the model to resolve the flow phenomena not amenable to modeling by the closure and, as such, the computations’ fidelity.

Navier Stokes equations↗
DOE OSTI2021

Simulation of flame acceleration and deflagration-to-detonation transition with heat transfer in HE product and fractured HE

A simple heat transfer model has been developed for high explosives. For the gaseous high explosive (HE) product, the Smargorinsky eddy viscosity model has been adopted for estimating the turbulent thermal conductivity. Convective heat transfer in the fractured solid HE is modeled as the effective thermal conduction in the porous media, which is equivalent to the enthalpy flux carried by the velocity determined by Darcy’s law with Forchheimer term.Other features of the model include a multi-reaction chain, equation-of-state for each species involved, and a two-phase hotspot model. Test calculations are performed with various initial disturbances. Onset and progress of flame acceleration (FA) and chemical reactions are strongly affected by the heat transfer and initial disturbance. Staged development of FA also affirms the necessity of modeling based on a multi-reaction chain. When conditions warrant,e.g., stronger initial disturbance, FA evolves to deflagration-to-detonation transition (DDT),demonstrating automatic capture of the DDT.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Pressure-Gradient-Based RANS Model for Predicting Separation in Transitional and Turbulent Flows

Predicting flow separation poses a significant challenge for RANS models, particularly in transitional flows over airfoils. We propose a novel improvement to RANS models to predict incipient separation in both transitional and fully turbulent flows. Our approach modifies the eddy viscosity model in regions indicated by a pressure-gradient criterion that accounts for intermittency - determining whether the boundary layer is laminar or turbulent. This model demonstrates robust generalization across various airfoil shapes and Reynolds numbers. Applied to the NREL Phase VI wind turbine rotor, our model shows improved aerodynamic performance predictions compared to the baseline RANS model.

k-omega SST↗
DOE OSTI2022

Sensitivity Analysis of Wind and Turbulence Predictions With Mesoscale‐Coupled Large Eddy Simulations Using Ensemble Machine Learning

Abstract Coupling between mesoscale models and large‐eddy simulation (LES) models is increasingly used to more realistically represent the wide range of scales of atmospheric motions affecting boundary layer winds and turbulence that need to be simulated accurately for applications such as wind energy. However, such mesoscale‐to‐microscale coupled modeling frameworks are potentially affected by a large number of uncertain closure parameters. Here, we investigate the sensitivity associated with six closure parameters related to a 1.5‐order subgrid‐scale turbulence closure for an ensemble of mesoscale‐coupled LES. The simulations are performed using the Weather Research and Forecasting model nested from horizontal resolutions of greater than a kilometer down to tens of meters. Closure parameters are varied to generate perturbed parameter ensembles for two case studies of highly sheared, convective boundary layers observed in the Columbia Basin of Oregon and Washington during the Second Wind Forecast Improvement Project. Machine learning algorithms are used to explore the sensitivity of LES predictions, considering the effects of the perturbed physical parameters alongside categorical factors such as the case study identity, measurement location, and LES resolution. For the conditions we examine, a single parameter, the eddy viscosity coefficient, is the dominant source of parametric sensitivity and its importance is comparable to the categorical factors for several of the simulation response variables we examine.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Water level changes in Lake Erie drive 21st century CO 2 and CH 4 fluxes from a coastal temperate wetland

Wetland water depth influences microbial and plant communities, which can alter the above- and below-ground carbon cycling of a wetland. Wetland water depths are likely to change due to shifting precipitation patterns, which will affect projections of greenhouse gas emissions; however, these effects are rarely incorporated into wetland greenhouse gas models. Seeking to address this gap, we used a mechanistic model, ecosys, to simulate a range of water depth scenarios in a temperate wetland, and analyzed simulated predictions of carbon dioxide (CO 2 ) andmethane (CH 4 ) fluxes over the 21st century. We tested our model using eddy covariance measurements of CO 2 and CH 4 fluxes collected at the Old Woman Creek National Estuarine Research Reserve (OWC) during 2015 and 2016. OWC is a lacustrine, estuarine, freshwater, mineral-soil marsh. An empirical model found that the wetland water depth is highly dependent on the water depth of the nearby Lake Erie. Future wetland surface water depths were modeled based on projection of Lake Erie's water depth using four separate NOAA projections, resulting in four wetland water-depth scenarios. Two of the four 21st century projections for Lake Erie water depths used in this study indicated that the water depth of the wetland would remain nearly steady; however, the other two indicated decreases in the wetland water depth. In our scenario where the wetland dries out, we project the wetland's climatological warming effect will decrease due to smaller CH 4 fluxes to the atmosphere and larger CO 2 uptake by the wetland. We also found that increased water level can lower emissions by shifting the site towards more open water areas, which have lower CH 4 emissions. We found that decreased water depths would cause more widespread colonization of the wetland by macrophyte vegetation. Using an empirical relationship, we also found that further drying could result in other, non-wetland vegetation to emerge, dramatically altering soil carbon cycling. In three of our four projections, we found that in general the magnitude of CO 2 and CH 4 fluxes steadily increase over the next 100 years in response to higher temperatures. However, in our driest simulations, we projected a different response due to increased oxidation of soil carbon, with CH 4 emissions decreasing substantially from an annual cumulative peak of 224.6 to a minimum of 104.7 gC m –2 year –1 . In that same simulation, net cumulative flux of CO 2 changed from being a sink of 56.5 gC m –2 year –1 to a source of 369.6 gCm –2 year –1 over the same period, despite a temperature increase from 13.7 °C to 14.2 °C. This temperature shift in our other three cases with greater water depths increased the source strength of CH 4 and the sink strength of CO 2 . Here, we conclude that the magnitude of wetland greenhouse-gas fluxes depended on the water depth primarily as it affected the areal percentage of the wetland available for plant colonization, but dramatic decreases in water depths could cause significant reductions in the wetland CH 4 fluxes, while simultaneously altering the wetland vegetation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Direct Numerical Simulation of Flow Over a Wall-Mounted Cube with the Nek5000 Spectral Element Code: DNS at Re h = 3900

The direct numerical simulation (DNS) of the canonical wall-mounted unit cube subjected to two distinct incident velocity profiles is performed at a Reynolds number where the bulk flow characteristics are known to become relatively Reynolds number insensitive. The aim of this work is to highlight the sensitivity of such bluff-body flows to mean shear and to provide a representative set of data for such flow scenarios where common turbulence modeling techniques often fail. Simple boundary conditions and a small domain are selected to allow for cost-effective and easy comparisons for model development purposes. In addition to mean velocity profiles, select turbulence statistics are presented in detail. Further, the basic efficacy of eddy viscosity models is evaluated and found to be adequate only for the shear stress components for the bluff-body flows examined here. Other failure mechanisms for RANS are proposed.

42 ENGINEERING↗
DOE OSTI2020

Theoretical Understanding of the Linear Relationship between Convective Updrafts and Cloud-Base Height for Shallow Cumulus Clouds. Part II: Continental Conditions

This is the Part II of a two-part study that seeks a theoretical understanding of an empirical relationship for shallow cumulus clouds: subcloud updraft velocity covaries linearly with the cloud-base height. This work focuses on continental cumulus clouds that are more strongly forced by surface fluxes and more deviated from equilibrium than those over oceans (Part I). We use a simple analytical model for shallow cumulus that is well tested against a high-resolution (25 m in the horizontal) large-eddy simulation model. Consistent with a conventional idea, we find that surface Bowen ratio is the key variable that regulates the covariability of both parameters: under the same solar insolation, a drier surface allows for stronger buoyancy flux, triggering stronger convection that deepens the subcloud layer. We find that the slope of the Bowen-ratio-regulated relationship between the two parameters (defined as λ) is dependent on both the local time and the stability of the lower free atmosphere. The value of λ decreases with time exponentially from sunrise to early afternoon and linearly from early afternoon to sunset. The value of λ is larger in a more stable atmosphere. In addition, continental λ in the early afternoon more than doubles the oceanic λ. Validation of the theoretical results against ground observations over the Southern Great Plains shows a reasonable agreement. Physical mechanisms underlying the findings are explained from the perspective of different time scales at which updrafts and cloud-base height respond to a surface flux forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Revisiting AMOC Transport Estimates From Observations and Models

Abstract Reference level assumptions used to calculate the Atlantic meridional overturning circulation transports at the RAPID and Meridional Overturning Variability Experiment (MOVE) observing arrays are revisited in an eddying ocean model. Observational transport calculation methods are complemented by several alternative approaches. At RAPID, the model transports from the observational method and the model truth (based on the actual model velocities) agree well in their mean and variability. There are substantial differences among the transport estimates obtained with various methods at the MOVE site. These differences result from relatively large and time‐varying reference velocities at depth in the model, not supporting a level‐of‐no‐motion. The methods that account for these reference velocities properly at MOVE produce transports that are in good agreement with the model truth. In contrast with the observational estimates, the model transport trends at MOVE and RAPID largely agree with each other on pentadal to multi‐decadal time scales.

Danabasoglu, Gokhan↗
DOE OSTI2021

Simulated Dust Transport in the Convective Boundary Layer

Soil dusts are an important source of aerosol in agricultural regions and can affect the Earth's radiation budget through the modification of cloud properties, and in particular, through their ability to act as ice nucleating particles. In order to impact cloud properties, agricultural soil dusts need to be transported from the point of emission to cloud–relevant altitudes. Vertical transport within the planetary boundary layer is strongly controlled by turbulence and is challenging to represent accurately in regional and global models. Large–eddy simulations (LES) are run at resolutions capable of resolving most of the turbulent energy directly and can thus better simulate vertical transport. In this study, we leverage the LES ARM Symbiotic Simulation and Observation (LASSO) large–eddy simulations to simulate vertical transport of agricultural dust within the turbulent boundary layer using a modified version of the stochastic dispersion model FLEXPART–WRF. We find that the modified model is better capable of simulating particle transport due to turbulence, and that particle size was the greatest factor in determining particle lifetime. Individual meteorology and particle density had intermediate effects upon particle transport, while release height had little effect upon simulation results. Finally, we utilize a quasi–single column model (QSCM) approach to determine how our results compare to a parameterized treatment of turbulence. The QSCM simulations led to greater tracer transport out of the boundary layer, with ramifications for any studies utilizing a Lagrangian stochastic model to understand tracer dispersion. Furthermore, these results highlight the importance of accurately simulating turbulence for understanding particle transport.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Increased power gains from wake steering control using preview wind direction information

Abstract. Yaw controllers typically rely on measurements taken at the wind turbine, resulting in a slow reaction to wind direction changes and subsequent power losses due to misalignments. Delayed yaw action is especially problematic in wake steering operation because it can result in power losses when the yaw misalignment angle deviates from the intended one due to a changing wind direction. This study explores the use of preview wind direction information for wake steering control in a two-turbine setup with a wind speed in the partial load range. For these conditions and a simple yaw controller, results from an engineering model identify an optimum preview time of 90 s. These results are validated by forcing wind direction changes in a large-eddy simulation model. For a set of six simulations with large wind direction changes, the average power gain from wake steering increases from only 0.44 % to 1.32 %. For a second set of six simulations with smaller wind direction changes, the average power gain from wake steering increases from 1.24 % to 1.85 %. Low-frequency fluctuations are shown to have a larger impact on the performance of wake steering and the effectiveness of preview control, in particular, than high-frequency fluctuations. From these results, it is concluded that the benefit of preview wind direction control for wake steering is substantial, making it a topic worth pursuing in future work.

17 WIND ENERGY↗
DOE OSTI2022

Evaluation of the Fitch Wind-Farm Wake Parameterization with Large-Eddy Simulations of Wakes Using the Weather Research and Forecasting Model

Abstract Wind-farm parameterizations in weather models can be used to predict both the power output and farm effects on the flow; however, their correctness has not been thoroughly assessed. We evaluate the wind-farm parameterization of the Weather Research and Forecasting Model with large-eddy simulations (LES) of the wake performed with the same model. We study the impact on the velocity and turbulence kinetic energy (TKE) of inflow velocity, roughness, resolution, number of turbines (one or two), and inversion height and strength. We compare the mesoscale with the LES by spatially averaging the LES within areas correspondent to the mesoscale horizontal spacing: one covering the turbine area and two downwind. We find an excellent agreement of the velocity within the turbine area between the two types of simulations. However, within the same area, we find the largest TKE discrepancies because in mesoscale simulations, the turbine-added TKE has to be highest at the turbine position to be advected downwind. Within the downwind areas, differences between velocities increase as the wake recovers faster in the LES, whereas for the TKE both types of simulations show similar levels. From the various configurations, the impact of inversion height and strength is small for these heights and inversion levels. The highest impact for the one-turbine simulations appears under the low-speed case due to the higher thrust, whereas the impact of resolution is low for the large-eddy simulations but high for the mesoscale simulations. Our findings demonstrate that higher-fidelity simulations are needed to validate wind-farm parameterizations.

17 WIND ENERGY↗
DOE OSTI2021

Generation of mean flows in rotating anisotropic turbulence: The case of solar near-surface shear layer

Results from helioseismology indicate that the radial gradient of the rotation rate in the near-surface shear layer (NSSL) of the Sun is independent of latitude and radius. Theoretical models using the mean-field approach have been successful in explaining this property of the NSSL, while global direct or large-eddy magnetoconvection models have so far been unable to reproduce this. We investigate the reason for this discrepancy by measuring the mean flows, Reynolds stress, and turbulent transport coefficients under conditions mimicking those in the solar NSSL. Simulations with as few ingredients as possible to generate mean flows were studied. These ingredients are inhomogeneity due to boundaries, anisotropic turbulence, and rotation. The parameters of the simulations were chosen such that they matched the weakly rotationally constrained NSSL. The simulations probe locally Cartesian patches of the star at a given depth and latitude. The depth of the patch was varied by changing the rotation rate such that the resulting Coriolis numbers covered the same range as in the NSSL. We measured the turbulent transport coefficient relevant for the nondiffusive (Λ-effect) and diffusive (turbulent viscosity) parts of the Reynolds stress and compared them with predictions of current mean-field theories. A negative radial gradient of the mean flow is generated only at the equator where meridional flows are absent. At other latitudes, the meridional flow is comparable to the mean flow corresponding to differential rotation. We also find that the meridional components of the Reynolds stress cannot be ignored. Additionally, we find that the turbulent viscosity is quenched by rotation by about 50% from the surface to the bottom of the NSSL. Our local simulations do not validate the explanation for the generation of the NSSL from mean-field theory where meridional flows and stresses are neglected. However, the rotational dependence of the turbulent viscosity in our simulations agrees well with theoretical predictions. Moreover, our results agree qualitatively with global convection simulations in that an NSSL can only be obtained near the equator.

79 ASTRONOMY AND ASTROPHYSICS↗