Momentum Advection Research Articles

Quantifying a common inconsistency in RANS-VoF modeling of water and oil core annular flow

Current trends in modeling CAF of oil and water combine the interface capturing methodology of VoF with the computational savings of a RANS approach. As demonstrated mathematically, this results in an inconsistency in the overall RANS-VoF treatment and incurs a number of omissions in the solution of oil fraction advection and momentum. To quantify these inconsistencies, five CAF cases with increasing Rew are considered and solved via DNS. Symptoms of the inconsistencies include qualitative errors in the prediction of the flow behavior in the entrance and fully-developed regions, as well as in the transition between these two regions. From a momentum perspective, the more troubling issue is the absence of unclosed terms resulting from fluctuating viscous and surface tension forces in the RANS equations. For momentum advection, terms associated with fluctuating density can be safely ignored due to the similar magnitude between oil and water density.

A regional NEMO 4.0 configuration of the subpolar North Atlantic

The dynamics of the North Atlantic subpolar gyre (SPG) largely control the Atlantic Meridional Overturning Circulation (AMOC), thus affecting the long-term climate variability of the global ocean. Many key processes in the North Atlantic SPG, such as Gulf Stream detachment, eddy activity, deep convection and overflows, are frequently incorrectly presented or misrepresented in modern ocean general circulation models. Here, we present a new regional eddy-resolving model of the subpolar North Atlantic based on the NEMO4 configuration and analyze its ability to adequately represent the dynamics of the subpolar North Atlantic. The model performance is assessed using in situ data (cross-section hydrography, Argo and sea gliders), remote sensing data (satellite sea surface temperature and salinity) and the ANDRO and ARMOR3D products based on optimal interpolation products. In the next step, we analyze the effects of high-resolution forcing and several recent developments in ocean modeling to improve the model solution. The implementation of the self-diffusive momentum advection scheme with additional viscosity allows for the simulation of more realistic Deep Western Boundary Current and Irminger Rings lifecycle characteristics. The surface current feedback parameterization implemented in the forcing function and the modification of the kinetic energy budget in the upper ocean, especially in intensively eddying areas, simulate the behavior predicted by coupled models. While the model is sensitive to the implementation of cool skin–warm layer parameterization, this parameterization does not significantly improve the solution compared to the observational data. Regional implementation of the local-sigma vertical coordinate system in cascading areas allows overflow waters to spread over the ocean bottom more realistically. The remote effect of the implementation of the local-sigma coordinate is traced in the vertical thermohaline structure at the western slopes of the Irminger Sea and Iceland Basin as well as in the Deep Western Boundary Current in the Labrador Sea. Compared to z-partial step setting, local-sigma allows for correct representation of continuous densification of the water in the overturning part of the AMOC and thus North Atlantic Deep-Water formation. The potential of the new model configuration for analyzing the subpolar North Atlantic, including convective activity, is discussed.

Influence of freshwater on the vertical structure of tidal currents: A case study of the Pearl River Estuary

Field observations in the Pearl River Estuary (PRE) indicate that vertical structures of tidal currents vary in different regions, which is associated with freshwater-induced stratification. A three-dimensional unstructured grid model (PRE-model) based on FVCOM is developed to study the mechanism of vertical structural changes in tidal currents. With the effect of freshwater, the amplitude of tidal currents increases in the upper layers and decreases in the bottom layers, and in regions with geometrical complexity, it is maximum in the middle depths. Momentum analysis indicates that vertical friction (VF), pressure gradient forces (PGF), and momentum advection (ADV) are the main factors leading to changes in tidal vertical structures. Variations of VF induced by stratification lead to an increase of tidal currents amplitude in the upper layers while a decrease in the bottom layers. Changes in PGF and ADV can also have significant impacts on tidal currents. Tidal frequency variations of density gradients result in pronounced tidal frequency baroclinic PGF in the bottom depths and enlarge total PGF. The amplitude of tidal currents increases at corresponding depths, which leads to subsurface maxima of tidal currents. ADV is special as its phase is variable relative to tidal currents, which can enhance or weaken tidal currents in different regions. The generation of tidal frequency baroclinic PGF is highly correlated with nonlinear processes, and the horizontal advection is considered the primary source of baroclinic PGF. The similar sources of ADV and baroclinic PGF also make them dominant in regions with geometrical complexity or large river plumes.

A Puzzling Quasi‐Periodic Variability in the Tropical Middle Atmosphere

AbstractThe Quasi‐Biennial Oscillation and the Semiannual Oscillation have been identified to be the leading modes of variability in the tropical middle atmosphere. With reanalysis data and independent rocket soundings from a low latitude site, we report the existence of yet another variability in the tropical lower mesosphere which is primarily evident as easterly bursts in zonal winds during the months of May‐July. It occurs with a variable interval of 2–5 yrs in the late 20th century and 7–9 yrs in the early 21st century. These Quasi‐Periodic Easterly Bursts are found to have remote influences on the Antarctic polar vortex as well as residual circulation in the lower mesosphere. We identify a potential causative mechanism for the easterly bursts that involve enhanced cross equatorial advection of momentum as well as gravity wave drag. A close association with Quasi Biennial Oscillation winds is observed, however, cause of the observed periodicity remains elusive.

A New WENO‐Based Momentum Advection Scheme for Simulations of Ocean Mesoscale Turbulence

AbstractCurrent eddy‐permitting and eddy‐resolving ocean models require dissipation to prevent a spurious accumulation of enstrophy at the grid scale. We introduce a new numerical scheme for momentum advection in large‐scale ocean models that involves upwinding through a weighted essentially non‐oscillatory (WENO) reconstruction. The new scheme provides implicit dissipation and thereby avoids the need for an additional explicit dissipation that may require calibration of unknown parameters. This approach uses the rotational, “vector invariant” formulation of the momentum advection operator that is widely employed by global general circulation models. A novel formulation of the WENO “smoothness indicators” is key for avoiding excessive numerical dissipation of kinetic energy and enstrophy at grid‐resolved scales. We test the new advection scheme against a standard approach that combines explicit dissipation with a dispersive discretization of the rotational advection operator in two scenarios: (a) two‐dimensional turbulence and (b) three‐dimensional baroclinic equilibration. In both cases, the solutions are stable, free from dispersive artifacts, and achieve increased “effective” resolution compared to other approaches commonly used in ocean models.

Changes in the tropical upper-tropospheric zonal momentum balance due to global warming

Abstract. We study the zonal momentum budget of the deep upper tropics in the context of present and future climates. In the zonal mean, as is known, a robust balance exists between the acceleration by the horizontal eddy momentum flux convergence and the deceleration by the mean meridional momentum advection. During summer, climatological stationary Rossby waves over the Asian monsoon longitudes converge westerly momentum into the tropics and are the primary contributors to the eddy term. During winter, anomalous westerly winds over the tropical east Pacific allow extratropical waves to propagate into the deep tropics, where they tend to break and decelerate the flow. When integrated over all longitudes, we find that eddies from these two regions sum constructively in summer and destructively in winter, always yielding a net positive momentum forcing that balances the mean flow term. The state-of-the-art CMIP6 suite qualitatively captures these features in the historical run and shows that the momentum fluxes change due to global warming. In summer, stationary eddy circulations in the Asian monsoon zone weaken in the upper troposphere (UT) but strengthen in the lower stratosphere (LS). Greater upward mass flux from the troposphere forces a stronger divergence and a more intense circulation in the LS following a Sverdrup vorticity balance. We observe this strengthening of summertime tropical and subtropical stationary waves in the LS over all longitudes and verify it in an idealized aquaplanet general circulation model experiment. In winter, we find that eddy westerlies over the east Pacific longitudes decrease in strength due to the expected weakening of the subtropical stationary waves with warming. This causes a significant decrease in the propagation of extratropical waves into this region, along with a drop in eddy potential vorticity fluxes associated with these waves. Thus, apart from the mean meridional flux, which weakens due to the projected weakening of the Hadley cells, our analysis of warming simulations clearly suggests significant and robust changes in the eddy momentum fluxes in the deep tropics. Potential implications of these changes in the context of the zonal mean flow and regional circulations are discussed.

Nonlinear processes in tsunami simulations for the Peruvian coast with focus on Lima and Callao

Abstract. This investigation addresses the tsunami inundation in Lima and Callao caused by the massive 1746 earthquake (Mw 9.0) along the Peruvian coast. Numerical modeling of the tsunami inundation processes in the nearshore includes strong nonlinear numerical terms. In a comparative analysis of the calculation of the tsunami wave effect, two numerical codes are used, Tsunami-HySEA and TsunAWI, which both solve the shallow water (SW) equations but with different spatial approximations. The comparison primarily evaluates the flow velocity fields in inundated areas. The relative importance of the various parts of the SW equations is determined, focusing on the nonlinear terms. Particular attention is paid to the contribution of momentum advection, bottom friction, and volume conservation. The influence of the nonlinearity on the degree and volume of inundation, flow velocity, and small-scale fluctuations is determined. The sensitivity of the solution concerning the bottom friction parameter is also investigated, showing the effects of nonlinearity processes in the inundated areas, wave heights, current velocity, and the spatial structure variations shown in tsunami inundation maps.

Extracting engineering-relevant wind information from research radar measurements

Historically, engineering-oriented characterization of low-level atmospheric winds has been dependent upon obtaining relevant measurements from anemometry, but recent research has shown promise in extracting relevant scales using radar measurements. Central to the success of these new analysis concepts was the proper correction of the acquired data fields to account for momentum advection and the implementation of a robust space-to-time conversion methodology. Combined, these facets allowed researchers to interpolate high resolution “tower-like” time histories from across the entire radar analysis domain, resulting in a more comprehensive evaluation of flow field evolution. However, the previous methodologies failed to properly characterize highly complex wind flows that contain wind speed and direction gradients. This study built upon these previous methodologies to better characterize complex wind flows, while also preserving engineering-relevant scales of motion. The new methodology was implemented on a wind record collected from a passing thunderstorm and included several sharp wind speed and direction changes. Validation of the resulting dual-Doppler-derived second-by-second projections was completed through comparison with available anemometry. The results highlight the significant advantages of using research radar information to bolster engineering perspectives, while also documenting the current limitations of the method.

Improvement of Global Forecast of Tropical Cyclone Intensity by Spray Heat Flux and Surface Roughness

AbstractIn global forecasting systems, tropical cyclone (TC) intensity is usually underestimated, which is one of the major challenges for TC forecasting. One possible reason is the deficiency of physical parameterizations associated with ocean surface waves. To improve the TC intensity forecast, effects of two powerful wave‐related processes, the spray‐mediated heat flux and surface roughness, are considered in a global coupled ocean‐atmosphere‐wave system (CFSv2.0‐WW3). Serial 24 hr forecasting experiments for all TCs from May to October 2019 are conducted, and comparisons are made against the International Best Track TC data. The results show that for strong TCs the underestimation of TC intensity in CFSv2.0‐WW3 is significantly improved by the modified surface roughness parameterization and the accelerated spray heat flux parameterization based on Gaussian Quadrature. For major hurricanes with initial wind speeds above 30 m/s, the error of TC maximum wind speed decreases by about 34.4%. To understand the associated dynamic and thermodynamic processes, a series of 120 hr simulations for major hurricane Hagibis is conducted. The effect of spray‐mediated heat flux accounts for 91.7% of the TC intensity improvement. The increased heat fluxes offset the negative contributions of vertical advection and ventilation, leading to enhanced entropy. The low‐level pressure is then decreased, resulting in enhanced convergence of absolute vorticity, which overpowers the negative effects of frictional dissipation and vertical momentum advection. The study indicates that the application of these two wave‐related parameterizations can potentially improve the global forecast and reanalysis of TC intensity.

Impact of Direct Radar Reflectivity Data Assimilation on the Simulation of Mesoscale Descending Inflow and Secondary Eyewall Formation in Hurricane Matthew (2016)

AbstractThe impact of assimilating ground‐based radar reflectivity on the rainband structure and secondary eyewall formation (SEF) of Hurricane Matthew (2016) is investigated within the framework of the Hurricane Weather Research and Forecasting model and its hybrid three‐dimensional ensemble‐variational data assimilation (DA) system. Compared to the control experiment (no radar reflectivity DA), the radar reflectivity DA experiment shows a clear signal of concentric eyewall and eyewall replacement cycle. Results demonstrate that radar reflectivity DA improves the stratiform rainband analysis, resulting in the mid‐level cooling associated with mesoscale descending inflow (MDI). The MDI further contributes to the low‐level acceleration maximum with boundary layer dynamics and triggers new convective updrafts in the SEF region. Momentum budget analysis also suggests that the mean vertical advection of absolute angular momentum plays an important role in the local momentum tendency in the SEF region in Hurricane Matthew (2016).

Angular momentum transport near convective-core boundaries of Gamma Doradus stars

Recent asteroseismic studies have revealed that the convective core of γ Doradus stars rotates faster than their radiative interior. We study the development of differential rotation near the convective core to test angular momentum transport processes that are typically adopted in stellar evolution models. Models that only include the advection of angular momentum by meridional circulation and shear instabilities cannot reproduce current rotational constraints, irrespective of the initial conditions. The latest formulation of internal magnetic fields based on the Tayler instability is indeed able to reproduce the internal rotation rate of post-main sequence stars; however, it appears too efficient during the main sequence and has thus been disfavoured. A less efficient version of the same transport process can simultaneously reproduce the rotation rate of the convective core, the rotation rate in radiative regions as probed by gravity-modes, and the surface rotational velocities of γ Doradus stars. Our work suggests that there are additional physical processes apart from internal magnetic fields at work in the stellar interiors of post-main sequence stars.

A Perturbative Solution for Nonlinear Stratified Upwelling over a Frictional Slope

Abstract A perturbative solution of simplified primitive equations for nonlinear weakly stratified upwelling over a frictional slope is found that resolves the vertical structure of velocity fields and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The solution uses assumptions consistent with the model proposed by Lentz and Chapman, including a steady-state, constant cross-shore density gradient, no alongshore gradients, laterally inviscid, and consideration of cross-shore advection of alongshore momentum. The solution resolves the vertical structure of velocity fields (including subsurface maxima of compensational flow, not resolved by Lentz and Chapman) and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The dynamics are similar to Lentz and Chapman; bottom stress balances alongshore wind stress in a homogeneous density ocean and is replaced by nonlinear cross-shore transport of alongshore momentum as the Burger number (S = αN/f, where α, N, and f are the bottom slope, buoyancy frequency, Coriolis frequency, respectively) increases. When the solution uses the empirical relation between cross-shore and vertical density gradients proposed by Lentz and Chapman, vorticity conservation is not satisfied and the nonlinear momentum transport estimated by the solution linearly increases with S, asymptotically matching Lentz and Chapman for S < 1. When the solution conserves interior potential vorticity, the momentum transport is proportional to S2 for S < 1 and is in better agreement with numerical simulations.

COMPARATIVE ASSESSMENT OF NON-CONSERVATIVE AND CONSERVATIVE RANS FORMULATIONS FOR COASTAL APPLICATIONS INVOLVING BREAKING WAVES

Two-phase Reynolds Averaged Navier-Stokes (RANS) simulations of breaking waves are susceptible to an unphysical thickening of the plunging crest. This is often times incorrectly attributed to deficiencies in the interface capturing scheme although, in reality, the issue stems from the nature of density treatment in momentum advection. If the density is considered face-centered in the advection term, the resulting formulation is conservative whilst if the density is modeled as a cell-centered quantity, the resulting formulation is non-conservative. Despite both approaches having been extensively applied to wave-breaking simulations in the literature, there is no study comparing both formulations for the same breaking scenario. In the present paper, we extensively compare both formulations for several depth- and steepness-induced breaking problems simulated using our in-house solver: IITM-RANS3D. Through these simulations, our work successfully addresses the following research questions: (a) how and why density treatment affects the physics of overturning and subsequently the plunging jet’s topology and (b) what are the implications of choosing a particular formulation for simulating a violent wave-structure interaction scenario?

An unstructured finite-volume level set / front tracking method for two-phase flows with large density-ratios

We extend the unstructured LEvel set / froNT tracking (LENT) method [1,2] for handling two-phase flows with strongly different densities (high-density ratios) by providing the theoretical basis for the numerical consistency between the mass and momentum conservation in the collocated Finite Volume discretization of the single-field two-phase Navier-Stokes equations. Our analysis provides the theoretical basis for the mass conservation equation introduced by Ghods and Herrmann [3] and used in [4–8]. We use a mass flux that is consistent with mass conservation in the implicit Finite Volume discretization of the two-phase momentum convection term, and solve the single-field Navier-Stokes equations with our SAAMPLE segregated solution algorithm [2]. The proposed ρLENT method recovers exact numerical stability for the two-phase momentum advection of a spherical droplet with density ratios ρ−/ρ+∈[1,104]. Numerical stability is demonstrated for in terms of the relative L∞ velocity error norm, for density-ratios in the range of [1,104], dynamic viscosity-ratios in the range of [1,104] and very strong surface tension forces, for challenging mercury/air and water/air fluid pairings. In addition, the solver performs well in cases characterized by strong interaction between two phases, i.e., oscillating droplets and rising bubbles. The proposed ρLENT method1 is applicable to any other two-phase flow simulation method that discretizes the single-field two-phase Navier-Stokes Equations using the collocated unstructured Finite Volume Method but does not solve an advection equation for the phase indicator using a flux-based approach, by adding the proposed geometrical approximation of the mass flux and the auxiliary mass conservation equation to the solution algorithm.

Variational RANS modeling of hydraulic jumps

The hydraulic jump is a phenomenon which frequently appears in nearshore, estuarine and riverine regions; it plays a major role for example in the development and migration of sand waves in mobile beds, and in controlling the propagation of infra-gravity water waves from the nearshore into the river through bar-built, shallow estuaries. Although the prediction of these natural flows becomes crucial in many circumstances, it is difficult because hydraulic jumps are complex flows, and their modeling is far from simple. The steady-state hydraulic jump may be undular or broken, and their physical properties differ significantly. Undular jump profiles are dominated by the deviation of the fluid pressure from hydrostatic conditions, with turbulence playing a secondary role. In contrast, broken hydraulic jumps are determined by the differential advection of momentum due to breaking, and turbulence stresses. Between both types, there is a complex transition involving the interaction of non-hydrostaticity and turbulence, which is difficult to mimic with depth-averaged models. Simulations of these flows are typically conducted by resorting to vertically-averaged models either using the Saint-Venant or the Serre-Green-Naghdi equations, but none of them are fully satisfactory. In this work, an alternative approach for hydraulic jumps is proposed by developing a variational Reynolds-Averaged Navier-Stokes (RANS) model which uses Kantorovich-Krylov expansions of the flow variables, and a novel depth-averaged, eddy viscosity approach suitable for non-hydrostatic flows. The model is validated through results of simulations of both undular and broken jumps compared against detailed experiments. New experiments were conducted in a simplified bar-built shallow river inlet; the variational RANS model was found to reproduce well the flow profile over the bar, the hydraulic jump position and dynamic pressures over the bar surface.

Investigating the physical mechanisms that modify wind plant blockage in stable boundary layers

Abstract. Wind plants slow down the approaching wind, a phenomenon known as blockage. Wind plant blockage undermines turbine performance for front-row turbines and potentially for turbines deeper into the array. We use large-eddy simulations to characterize blockage upstream of a finite-size wind plant in flat terrain for different atmospheric stability conditions and investigate the physical mechanisms modifying the flow upstream of the turbines. To examine the influence of atmospheric stability, we compare simulations of two stably stratified boundary layers using the Weather Research and Forecasting model in large-eddy simulation mode, representing wind turbines using the generalized actuator disk approach. For a wind plant, a faster cooling rate at the surface, which produces stronger stably stratified flow in the boundary layer, amplifies blockage. As a novelty, we investigate the physical mechanisms amplifying blockage by evaluating the different terms in the momentum conservation equation within the turbine rotor layer. The velocity deceleration upstream of a wind plant is caused by an adverse pressure gradient and momentum advection out of the turbine rotor layer. The cumulative deceleration of the flow upstream of the front-row turbines instigates vertical motions. The horizontal flow is diverted vertically, reducing momentum availability in the turbine rotor layer. Although the adverse pressure gradient upstream of the wind plant remains unchanged with atmospheric stability, vertical advection of horizontal momentum is amplified in the more strongly stable boundary layer, mainly by larger shear of the horizontal velocity, thus increasing the blockage effect.

Magnetic Braking during Direct Collapse Black Hole Formation

Magnetic fields are expected to be efficiently amplified during the formation of the first massive black holes via the small-scale dynamo and in the presence of strong accretion shocks occurring during gravitational collapse. Here, we analyze high-resolution cosmological magnetohydrodynamical simulations of gravitational collapse in atomic cooling halos, exploring the dynamical role of magnetic fields, particularly concerning the effect of magnetic braking and angular momentum transport. We find that after the initial amplification, magnetic fields contribute to the transport of angular momentum and reduce it compared to pure hydrodynamical simulations. However, the magnetic and Reynolds torques do not fully compensate for the inward advection of angular momentum, which still accumulates over timescales of ∼1 Myr. A Jeans analysis further shows that magnetic pressure strongly contributes to suppressing fragmentation on scales of 0.1–10 pc. Overall, the presence of magnetic fields thus aids in the transport of angular momentum and favors the formation of massive objects.

Evidence for the Influence of the Quasi-Biennial Oscillation on the Semiannual Oscillation in the Tropical Middle Atmosphere

Abstract The semiannual oscillation (SAO) in zonally averaged zonal winds develops just above the quasi-biennial oscillation (QBO) and dominates the seasonal variability in the tropical upper stratosphere and lower mesosphere. The magnitude, seasonality, and latitudinal structure of the SAO vary with the phase of the QBO. There is also an annual oscillation (AO) whose magnitude at the equator is smaller than those of the SAO and QBO but not negligible. This work presents the relation between the SAO, QBO, AO, and time-mean wind in the tropical upper stratosphere and lower mesosphere using winds derived from satellite geopotential height observations. The winds are generally more westerly during the easterly phase of the QBO. The SAO extends to lower altitudes during periods where the QBO is characterized by deep easterly winds. The differences in the SAO associated with the QBO are roughly confined to the latitudes where the QBO has appreciable amplitude, suggesting that the mechanism is controlled by vertical coupling. The westerly phases of the SAO and AO show downward propagation with time. This analysis suggests that forcing by dissipation of waves with westerly momentum is responsible for the westerly acceleration of both the SAO and AO. The timing and structure of the easterly phases of the SAO and AO near the stratopause are consistent with the response to meridional advection of momentum across the equator during solstices; it is not apparent that local wave processes play important roles in the easterly phases in the region of the stratopause.

On the magnitude and sensitivity of the quasi-biennial oscillation response to a tropical volcanic eruption

Abstract. Volcanic eruptions that inject sulfur dioxide into the stratosphere have the potential to alter large-scale circulation patterns, such as the quasi-biennial oscillation (QBO), which can affect weather and transport of chemical species. Here, we conduct simulations of tropical volcanic eruptions using the UM-UKCA aerosol-climate model with an explicit representation of the QBO. Eruptions emitting 60 Tg of SO2 (i.e. the magnitude of the 1815 Mt. Tambora eruption) and 15 Tg of SO2 (i.e. the magnitude of the 1991 Mt. Pinatubo eruption) were initiated at the Equator during two different QBO states. We show that tropical eruptions delay the progression of the QBO phases, with the magnitude of the delay dependent on the initial wind shear in the lower stratosphere and a much longer delay when the shear is easterly than when it is westerly. The QBO response in our model is driven by vertical advection of momentum by the stronger tropical upwelling caused by heating due to the increased volcanic sulfate aerosol loading. Direct aerosol-induced warming with subsequent thermal wind adjustment, as proposed by previous studies, is found to only play a secondary role. This interpretation of the response is supported by comparison with a simple dynamical model. The dependence of the magnitude of the response on the initial QBO state results from differences in the QBO secondary circulation. In the easterly shear zone of the QBO, the vertical component of the secondary circulation is upward and reinforces the anomalous upwelling driven by volcanic aerosol heating, whereas in the westerly shear zone the vertical component is downward and opposes the aerosol-induced upwelling. We also find a change in the latitudinal structure of the QBO, with the westerly phase of the QBO strengthening in the hemisphere with the lowest sulfate aerosol burden. Overall, our study suggests that tropical eruptions of Pinatubo magnitude or larger could force changes to the progression of the QBO, with particularly disruptive outcomes for the QBO if the eruption occurs during the easterly QBO shear.

Understanding Energy Pathways in the Gulf Stream

Abstract The Gulf Stream (GS) is one of the strongest ocean currents on the planet. Eddy-rich resolution models are needed to properly represent the dynamics of the GS; however, kinetic energy (KE) can be in excess in these models if not dissipated efficiently. The question of how and how much energy is dissipated and in particular how it flows through ocean scales thus remains an important and largely unanswered question. Using a high-resolution (∼2 km) ocean model [Coastal and Regional Ocean Community (CROCO)], we characterize the spatial and temporal distribution of turbulent cascades in the GS based on a coarse-grained method. We show that the balanced flow is associated with an inverse cascade while the forward cascade is explained by ageostrophic advection associated with frontogenesis. Downscale fluxes are dominant at scales smaller than about 20 km near the surface and most intense at the GS North Wall. There is also strong seasonal variability in KE flux, with the forward cascade intensifying in winter and the inverse cascade later in spring. The forward cascade, which represents an interior route to dissipation, is compared with both numerical and boundary dissipation processes. The contribution of interior dissipation is an order of magnitude smaller than that of the other energy sinks. We thus evaluate the sensitivity of horizontal momentum advection schemes on energy dissipation and show that the decrease in numerical dissipation in a high-order scheme leads to an increase in dissipation at the boundaries, not in the downscale flux.