Showing 1–10 of 10 results for author: Kafiabad, H

The exponential Lagrangian mean

Authors: Abhijeet Minz, Lois E. Baker, Hossein A. Kafiabad, Jacques Vanneste

Abstract: Lagrangian averaging is a valuable tool for the analysis and modelling of multiscale processes in fluid dynamics. The numerical computation of Lagrangian averages from simulation data is considered challenging, however. In response, we develop a straightforward form of Lagrangian (time) averaging -- exponential averaging -- and derive simple partial differential equations that govern the evolution… ▽ More Lagrangian averaging is a valuable tool for the analysis and modelling of multiscale processes in fluid dynamics. The numerical computation of Lagrangian averages from simulation data is considered challenging, however. In response, we develop a straightforward form of Lagrangian (time) averaging -- exponential averaging -- and derive simple partial differential equations that govern the evolution of the Lagrangian mean fields. These equations can be solved at minimum cost as part of the numerical simulation of fluid models, using the same time and space discretisation as for the dynamical equations. We implement exponential averaging in the rotating shallow-water model and demonstrate its effectiveness at filtering out large-amplitude Poincaré waves while retaining the salient features of an underlying turbulent flow. We generalise the computation of exponential Lagrangian means from scalar to tensor fields and apply the generalisation to the computation of the Lagrangian mean momentum and of the related pseudomomentum. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Inhomogeneity-induced wavenumber diffusion

Authors: Michael R. Cox, Hossein A. Kafiabad, Jacques Vanneste

Abstract: Inertia-gravity waves are scattered by background flows as a result of Doppler shift by a non-uniform velocity. In the WKB regime, the scattering process reduces to a diffusion in spectral space. Other inhomogeneities the waves encounter, such as density variations, also cause scattering and spectral diffusion. We generalise the spectral diffusion equation to account for these inhomogeneities. We… ▽ More Inertia-gravity waves are scattered by background flows as a result of Doppler shift by a non-uniform velocity. In the WKB regime, the scattering process reduces to a diffusion in spectral space. Other inhomogeneities the waves encounter, such as density variations, also cause scattering and spectral diffusion. We generalise the spectral diffusion equation to account for these inhomogeneities. We apply the result to the rotating shallow water system, for which height inhomogeneities arise from velocity inhomogeneities through geostrophy, and to the Boussinesq system for which buoyancy inhomogeneities arise similarly. We compare the contributions that height and buoyancy variations make to the spectral diffusion with the contribution of the Doppler shift. In both systems, we find regimes where all contributions are significant. We support our findings with exact solutions of the diffusion equation and with ray tracing simulations in the shallow water case. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Lagrangian filtering for wave-mean flow decomposition

Authors: Lois E. Baker, Hossein A. Kafiabad, Jacques Vanneste

Abstract: Geophysical flows are typically composed of wave and mean motions with a wide range of overlapping temporal scales, making separation between the two types of motion in wave-resolving numerical simulations challenging. Lagrangian filtering - whereby a temporal filter is applied in the frame of the flow - is an effective way to overcome this challenge, allowing clean separation of waves from mean f… ▽ More Geophysical flows are typically composed of wave and mean motions with a wide range of overlapping temporal scales, making separation between the two types of motion in wave-resolving numerical simulations challenging. Lagrangian filtering - whereby a temporal filter is applied in the frame of the flow - is an effective way to overcome this challenge, allowing clean separation of waves from mean flow based on frequency separation in a Lagrangian frame. Previous implementations of Lagrangian filtering have used particle tracking approaches, which are subject to large memory requirements or difficulties with particle clustering. Kafiabad and Vanneste (2023, KV23) recently proposed a novel method for finding Lagrangian means without particle tracking by solving a set of partial differential equations alongside the governing equations of the flow. In this work, we adapt the approach of KV23 to develop a flexible, on-the-fly, PDE-based method for Lagrangian filtering using arbitrary convolutional filters. We present several different wave-mean decompositions, demonstrating that our Lagrangian methods are capable of recovering a clean wave-field from a nonlinear simulation of geostrophic turbulence interacting with Poincaré waves. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Computing Lagrangian means

Authors: Hossein A. Kafiabad, Jacques Vanneste

Abstract: Lagrangian averaging plays an important role in the analysis of wave--mean-flow interactions and other multiscale fluid phenomena. The numerical computation of Lagrangian means, e.g. from simulation data, is however challenging. Typical implementations require tracking a large number of particles to construct Lagrangian time series which are then averaged using a low-pass filter. This has drawback… ▽ More Lagrangian averaging plays an important role in the analysis of wave--mean-flow interactions and other multiscale fluid phenomena. The numerical computation of Lagrangian means, e.g. from simulation data, is however challenging. Typical implementations require tracking a large number of particles to construct Lagrangian time series which are then averaged using a low-pass filter. This has drawbacks that include large memory demands, particle clustering and complications of parallelisation. We develop a novel approach in which the Lagrangian means of various fields (including particle positions) are computed by solving partial differential equations (PDEs) that are integrated over successive averaging time intervals. We propose two strategies, distinguished by their spatial independent variables. The first, which generalises the algorithm of Kafiabad (2022, J. Fluid Mech. 940, A2), uses end-of-interval particle positions; the second directly uses the Lagrangian mean positions. The PDEs can be discretised in a variety of ways, e.g. using the same discretisation as that employed for the governing dynamical equations, and solved on-the-fly to minimise the memory footprint. We illustrate the new approach with a pseudospectral implementation for the rotating shallow-water model. Two applications to flows that combine vortical turbulence and Poincare waves demonstrate the superiority of Lagrangian averaging over Eulerian averaging for wave--vortex separation. △ Less

Submitted 11 March, 2023; v1 submitted 4 August, 2022; originally announced August 2022.

Inertia-gravity-wave diffusion by geostrophic turbulence: the impact of flow time dependence

Authors: Michael R. Cox, Hossein A. Kafiabad, Jacques Vanneste

Abstract: The scattering of three-dimensional inertia-gravity waves by a turbulent geostrophic flow leads to the redistribution of their action through what is approximately a diffusion process in wavevector space. The corresponding diffusivity tensor was obtained by Kafiabad, Savva & Vanneste (2019, J. Fluid Mech., 869, R7) under the assumption of a time-independent geostrophic flow. We relax this assumpti… ▽ More The scattering of three-dimensional inertia-gravity waves by a turbulent geostrophic flow leads to the redistribution of their action through what is approximately a diffusion process in wavevector space. The corresponding diffusivity tensor was obtained by Kafiabad, Savva & Vanneste (2019, J. Fluid Mech., 869, R7) under the assumption of a time-independent geostrophic flow. We relax this assumption to examine how the weak diffusion of wave action across constant-frequency cones that results from the slow time dependence of the geostrophic flow affects the distribution of wave energy. We find that the stationary wave-energy spectrum that arises from a single-frequency wave forcing is localised within a thin boundary layer around the constant-frequency cone, with a thickness controlled by the acceleration spectrum of the geostrophic flow. We obtain an explicit analytic formula for the wave-energy spectrum which shows good agreement with the results of a high-resolution simulation of the Boussinesq equations. △ Less

Submitted 1 November, 2022; v1 submitted 19 July, 2022; originally announced July 2022.

Grid-based calculation of Lagrangian mean

Authors: Hossein Kafiabad

Abstract: Lagrangian averaging has been shown to be more effective than Eulerian mean in separating waves from slow dynamics in two-timescale flows. It also appears in many reduced models that capture the wave feedback on the slow flow. Its calculation, however, requires tracking particles in time, which imposes several difficulties in grid-based numerical simulations or estimation from fixed-point measurem… ▽ More Lagrangian averaging has been shown to be more effective than Eulerian mean in separating waves from slow dynamics in two-timescale flows. It also appears in many reduced models that capture the wave feedback on the slow flow. Its calculation, however, requires tracking particles in time, which imposes several difficulties in grid-based numerical simulations or estimation from fixed-point measurements. To circumvent these difficulties, we propose a grid-based iterative method to calculate Lagrangian mean without tracking particles in time. We explain how this method reduces computation, memory footprint and communication between processors in parallelised numerical models. To assess the accuracy of this method several examples are examined and discussed. We also explore an application of this method in the context of shallow-water equations by quantifying the validity of wave-averaged geostrophic balance -- a modified form of geostrophic balance accounting for the effect of strong waves on slow dynamics. △ Less

Submitted 8 March, 2022; v1 submitted 12 October, 2021; originally announced October 2021.

Interaction of near-inertial waves with an anticyclonic vortex

Authors: Hossein A. Kafiabad, Jacques Vanneste, William R. Young

Abstract: Anticyclonic vortices focus and trap near-inertial waves so that near-inertial energy levels are elevated within the vortex core. Some aspects of this process, including the nonlinear modification of the vortex by the wave, are explained by the existence of trapped near-inertial eigenmodes. These vortex eigenmodes are easily excited by an initial wave with horizontal scale much larger than that of… ▽ More Anticyclonic vortices focus and trap near-inertial waves so that near-inertial energy levels are elevated within the vortex core. Some aspects of this process, including the nonlinear modification of the vortex by the wave, are explained by the existence of trapped near-inertial eigenmodes. These vortex eigenmodes are easily excited by an initial wave with horizontal scale much larger than that of the vortex radius. We study this process using a wave-averaged model of near-inertial dynamics and compare its theoretical predictions with numerical solutions of the three-dimensional Boussinesq equations. In the linear approximation, the model predicts the eigenmode frequencies and spatial structures, and a near-inertial wave energy signature that is characterized by an approximately time-periodic, azimuthally invariant pattern. The wave-averaged model represents the nonlinear feedback of the waves on the vortex via a wave-induced contribution to the potential vorticity that is proportional to the Laplacian of the kinetic energy density of the waves. When this is taken into account, the modal frequency is predicted to increase linearly with the energy of the initial excitation. Both linear and nonlinear predictions agree convincingly with the Boussinesq results. △ Less

Submitted 16 October, 2020; originally announced October 2020.

Inertia-gravity-wave scattering by geostrophic turbulence

Authors: Miles A. C. Savva, Hossein A. Kafiabad, Jacques Vanneste

Abstract: In rotating stratified flows including in the atmosphere and ocean, inertia-gravity waves (IGWs) often coexist with a geostrophically balanced turbulent flow. Advection and refraction by this flow lead to wave scattering, redistributing IGW energy in the position--wavenumber phase space. We give a detailed description of this process by deriving a kinetic equation governing the evolution of the IG… ▽ More In rotating stratified flows including in the atmosphere and ocean, inertia-gravity waves (IGWs) often coexist with a geostrophically balanced turbulent flow. Advection and refraction by this flow lead to wave scattering, redistributing IGW energy in the position--wavenumber phase space. We give a detailed description of this process by deriving a kinetic equation governing the evolution of the IGW phase-space energy density. The derivation relies on the smallness of the Rossby number characterising the geostrophic flow, which is treated as a random field with known statistics, and makes no assumption of spatial scale separation. The kinetic equation describes energy transfers that are restricted to IGWs with the same frequency, as a result of the timescale separation between waves and flow. We formulate the kinetic equation on the constant-frequency surface -- a double cone in wavenumber space -- using polar spherical coordinates, and we examine the form of the two scattering cross sections involved, which quantify energy transfers between IGWs with, respectively, the same and opposite directions of vertical propagation. The kinetic equation captures both the horizontal isotropisation and the cascade of energy across scales that result from scattering. We focus our attention on the latter to assess the predictions of the kinetic equation against direct simulations of the three-dimensional Boussinesq equations, finding good agreement. △ Less

Submitted 5 August, 2020; originally announced August 2020.

Wave-averaged balance: a simple example

Authors: Hossein A. Kafiabad, Jacques Vanneste, William R. Young

Abstract: In the presence of inertia-gravity waves, the geostrophic and hydrostatic balance that characterises the slow dynamics of rapidly rotating, strongly stratified flows holds in a time-averaged sense and applies to the Lagrangian-mean velocity and buoyancy. We give an elementary derivation of this wave-averaged balance and illustrate its accuracy in numerical solutions of the three-dimensional Boussi… ▽ More In the presence of inertia-gravity waves, the geostrophic and hydrostatic balance that characterises the slow dynamics of rapidly rotating, strongly stratified flows holds in a time-averaged sense and applies to the Lagrangian-mean velocity and buoyancy. We give an elementary derivation of this wave-averaged balance and illustrate its accuracy in numerical solutions of the three-dimensional Boussinesq equations, using a simple configuration in which vertically planar near-inertial waves interact with a barotropic anticylonic vortex. We further use the conservation of the wave-averaged potential vorticity to predict the change in the barotropic vortex induced by the waves. △ Less

Submitted 11 November, 2020; v1 submitted 6 March, 2020; originally announced March 2020.

Journal ref: J. Fluid Mech. 911 (2021) R1

Diffusion of inertia-gravity waves by geostrophic turbulence

Authors: Hossein Kafiabad, Miles A. C. Savva, Jacques Vanneste

Abstract: The scattering of inertia-gravity waves by large-scale geostrophic turbulence in a rapidly rotating, strongly stratified fluid leads to the diffusion of wave energy on the constant-frequency cone in wavenumber space. We derive the corresponding diffusion equation and relate its diffusivity to the wave characteristics and the energy spectrum of the turbulent flow. We check the predictions of this e… ▽ More The scattering of inertia-gravity waves by large-scale geostrophic turbulence in a rapidly rotating, strongly stratified fluid leads to the diffusion of wave energy on the constant-frequency cone in wavenumber space. We derive the corresponding diffusion equation and relate its diffusivity to the wave characteristics and the energy spectrum of the turbulent flow. We check the predictions of this equation against numerical simulations of the three-dimensional Boussinesq equations in initial-value and forced scenarios with horizontally isotropic wave and flow fields. In the forced case, wavenumber diffusion results in a $k^{-2}$ wave energy spectrum consistent with as-yet-unexplained features of observed atmospheric and oceanic spectra. △ Less

Submitted 9 April, 2019; v1 submitted 11 February, 2019; originally announced February 2019.