pith. sign in

arxiv: 2605.30077 · v1 · pith:UYNYFRPInew · submitted 2026-05-28 · ⚛️ physics.flu-dyn

Two-way coupling of gravity waves and wind farm wakes: a reduced-order boundary-layer model

Hossein A. Kafiabad , Majid Bastankhah This is my paper

Pith reviewed 2026-06-29 05:42 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords gravity waveswind farm wakesboundary layer modelreduced-order modelcapping inversiontwo-way couplinglarge-eddy simulationpressure gradients
0
0 comments X

The pith

A reduced-order boundary-layer model captures two-way coupling between gravity waves and wind-farm wakes through a dynamic condition at the capping inversion.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a reduced-order framework for large wind farms that includes how gravity waves in the stratified atmosphere affect turbulent wakes and how the wakes in turn affect the waves. Linearized equations are written separately for the boundary layer and the free atmosphere, then joined by a dynamic boundary condition at the inversion that transmits pressure-gradient feedback. This setup lets the model predict upstream blockage from adverse pressure and faster wake recovery from favorable pressure. Comparisons with large-eddy simulations confirm that both the internal farm flow and the large-scale wave effects are reproduced. The approach supplies a computationally efficient alternative to full atmospheric simulations while keeping vertical structure.

Core claim

The two-way coupling between gravity waves and wind-farm wakes is represented by separate governing equations for the boundary layer and the overlying stratified free atmosphere that remain coupled through a dynamic boundary condition at the capping inversion; a mixed spectral-finite-difference discretization keeps the model efficient while preserving vertical structure, and the resulting predictions match the upstream blockage and accelerated wake recovery observed in large-eddy simulations.

What carries the argument

The dynamic boundary condition at the capping inversion that transmits the feedback of gravity waves onto the boundary-layer flow.

If this is right

  • The model predicts upstream blockage induced by adverse pressure gradients from gravity waves.
  • It reproduces accelerated wake recovery within and downwind of the farm driven by favourable pressure gradients.
  • The framework matches both internal wind-farm flow and large-scale gravity-wave effects in LES comparisons.
  • The mixed discretization retains vertical boundary-layer structure at modest computational cost.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Wind-farm power estimates that omit this coupling may understate losses from blockage in very large arrays.
  • The same linear-coupling structure could be tested in time-dependent or weakly nonlinear regimes by relaxing the linearization.
  • Integration into mesoscale weather models might improve forecasts for offshore farms where stratification is pronounced.
  • Analogous reduced-order treatments could be explored for other wave-boundary-layer problems such as urban or oceanic flows.

Load-bearing premise

Linearizing the non-hydrostatic Boussinesq equations and applying simplifications appropriate to each region produces equations that remain valid when the regions are joined by the dynamic boundary condition.

What would settle it

Large-eddy simulations of a wind farm in stratified flow that show neither upstream blockage nor accelerated wake recovery under conditions where the model predicts both would falsify the coupling framework.

Figures

Figures reproduced from arXiv: 2605.30077 by Hossein A. Kafiabad, Majid Bastankhah.

Figure 1
Figure 1. Figure 1: Interaction of a wind farm with a CNBL. The red and blue curves show schematic potential-temperature and background velocity profiles, respectively. the form of an abrupt temperature jump at the capping inversion and a more gradual variation in the free atmosphere above (figure 1). This structure is called a conventionally neutral boundary layer (CNBL). In this case, the wind farm’s interaction with the ov… view at source ↗
Figure 2
Figure 2. Figure 2: Schematic of pressures before (left) and after (right) displacing the interface. mountains that has been extensively studied in mountain-wave theories (see e.g., Durran 1990; Gill 2016). Here, instead of a hill disturbing the stratified flow, the displacement of the capping inversion acts as the disturbance to the free atmosphere (Khan et al. 2026). 2.4. Coupled dynamics at inversion layer So far, the boun… view at source ↗
Figure 3
Figure 3. Figure 3: Normalised perturbations of vertical velocity (a,b), streamwise pressure gradient (c,d), and streamwise velocity (e,f) for a TNBL against a CNBL. Turbine rows are shown by black lines. 0.01 s−1 ) with the presence of internal and interfacial waves. In both scenarios, the semi￾infinite wind farm consists of 𝑁𝑟 =16 turbine rows with a constant thrust coefficient 𝐶𝑇 = 0.75, rotor diameter 𝐷 = 100 m, hub heigh… view at source ↗
Figure 4
Figure 4. Figure 4: Comparison of normalised velocity perturbation at hub height with the LES data of Zhu et al. (2025). (b) Comparison of normalised velocity perturbation vertically averaged over 𝑧 = (0, 2𝑧ℎ) with the LES data reported in Allaerts & Meyers (2019). (c) Normalised capping inversion displacement 𝜂/𝐻 comparing with the same dataset as in (b). In (b) and (c), the shaded regions indicate the streamwise extents of … view at source ↗
read the original abstract

This paper develops a reduced-order framework for modelling the two-way coupling between gravity waves and turbulent wakes in large-scale wind farms. Linearising the non-hydrostatic Boussinesq equations and introducing simplifications appropriate to the boundary layer and the overlying stratified free atmosphere yield separate governing equations for the two regions. These are coupled through a dynamic boundary condition at the capping inversion, which directly captures the feedback of gravity waves on the boundary-layer flow. A mixed spectral-finite-difference discretisation yields a computationally efficient model while retaining vertical boundary-layer structure. Comparisons with large-eddy simulations (LES) confirm the model successfully reproduces both internal wind-farm flow and large-scale gravity-wave effects. It captures the upstream blockage induced by adverse pressure gradients, as well as the accelerated wake recovery within and downwind of the farm, driven by favourable pressure gradients.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper develops a reduced-order model for two-way coupling between gravity waves and wind-farm wakes. It linearizes the non-hydrostatic Boussinesq equations, applies boundary-layer and free-atmosphere simplifications to obtain separate governing equations, couples the regions via a dynamic condition at the capping inversion, discretizes with a mixed spectral-finite-difference scheme, and reports that the resulting model reproduces LES results for upstream blockage and accelerated wake recovery.

Significance. If the linearization and coupling remain valid, the framework supplies an efficient tool for capturing gravity-wave feedback on farm-scale flows that is absent from standard wake models; this would be useful for assessing blockage, recovery, and large-scale atmospheric interactions in wind-energy applications.

major comments (3)
  1. [Abstract, §4] Abstract and §4 (results): the claim that the model 'successfully reproduces' LES results for blockage and wake recovery is stated without quantitative error metrics (e.g., L2 norms, point-wise velocity or pressure differences, or R² values) or a stated validity range for the linearization; this prevents assessment of whether the reproduction is within acceptable engineering tolerance.
  2. [§2, §3] §2 (model derivation) and §3 (coupling): the statement that the linearized equations 'remain valid when coupled through the dynamic boundary condition' is asserted without an a-posteriori evaluation of the neglected nonlinear advection terms evaluated on the computed solution; for finite farm forcing this check is load-bearing for the central claim.
  3. [§3] §3 (discretization and boundary conditions): no discussion is provided of how the dynamic inversion condition is discretized or of any numerical tests confirming that the coupling does not re-introduce O(1) perturbations that violate the small-amplitude premise inside the farm.
minor comments (2)
  1. Notation for the inversion displacement and pressure jump should be defined once at first use and used consistently thereafter.
  2. Figure captions for the LES comparison panels should state the exact farm layout, inversion height, and stratification parameters used in each case.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We address each major comment below, indicating where revisions will be made to improve the manuscript.

read point-by-point responses

  1. Referee: [Abstract, §4] Abstract and §4 (results): the claim that the model 'successfully reproduces' LES results for blockage and wake recovery is stated without quantitative error metrics (e.g., L2 norms, point-wise velocity or pressure differences, or R² values) or a stated validity range for the linearization; this prevents assessment of whether the reproduction is within acceptable engineering tolerance.

    Authors: We agree that quantitative metrics are needed for a rigorous assessment. In the revised manuscript we will add L2 norms of velocity and pressure differences between the model and LES, selected point-wise comparisons, and an explicit statement of the linearization validity range based on perturbation amplitudes relative to the mean flow. These additions will clarify performance within engineering tolerances. revision: yes

  2. Referee: [§2, §3] §2 (model derivation) and §3 (coupling): the statement that the linearized equations 'remain valid when coupled through the dynamic boundary condition' is asserted without an a-posteriori evaluation of the neglected nonlinear advection terms evaluated on the computed solution; for finite farm forcing this check is load-bearing for the central claim.

    Authors: We acknowledge that an explicit a-posteriori check strengthens the central claim. The revised manuscript will include estimates of the neglected nonlinear advection terms evaluated on the computed solution, compared against the retained linear terms, to confirm they remain small for the farm forcings examined. revision: yes

  3. Referee: [§3] §3 (discretization and boundary conditions): no discussion is provided of how the dynamic inversion condition is discretized or of any numerical tests confirming that the coupling does not re-introduce O(1) perturbations that violate the small-amplitude premise inside the farm.

    Authors: We will expand the discretization section to detail the implementation of the dynamic inversion condition within the mixed spectral-finite-difference scheme. We will also add numerical tests showing that the coupled solution keeps perturbations inside the farm within the small-amplitude regime assumed by the linearization. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation from linearized equations is self-contained

full rationale

The paper derives its reduced-order model directly from linearization of the non-hydrostatic Boussinesq equations plus boundary-layer simplifications, followed by coupling at the capping inversion. No equations or parameters are shown to be fitted to the target data and then re-presented as predictions; no self-citation chain is invoked to justify uniqueness or an ansatz; and the central coupling step is presented as following from the linearized governing equations rather than reducing to the inputs by construction. LES comparisons function as external validation, not as the source of the model itself.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

No free parameters, axioms, or invented entities are identifiable from the abstract alone; the linearization and simplifications are stated without explicit listing of fitted constants or new postulated quantities.

pith-pipeline@v0.9.1-grok · 5677 in / 1189 out tokens · 21647 ms · 2026-06-29T05:42:50.157400+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

7 extracted references · 1 canonical work pages · 1 internal anchor

  1. [1]

    & Meyers, J.2017 Boundary-layer development and gravity waves in conventionally neutral wind farms.Journal of Fluid Mechanics814, 95–130

    Allaerts, D. & Meyers, J.2017 Boundary-layer development and gravity waves in conventionally neutral wind farms.Journal of Fluid Mechanics814, 95–130

    2017

  2. [2]

    ˚A., Berkesten H ¨agglund, P

    Ebenhoch, R., Muro, B., Dahlberg, J. ˚A., Berkesten H ¨agglund, P. & Segalini, A.2017 A linearized numerical model of wind-farm flows.Wind Energy20(5), 859–875

    2017

  3. [3]

    Modelling Farm-to-Farm Interaction Using a Fast Linearised Numerical Approach

    Everley, Alexia, Kafiabad, Hossein A. & Bastankhah, Majid2026 Modelling farm-to-farm interaction using a fast linearised numerical approach, arXiv: 2605.04782

    work page internal anchor Pith review Pith/arXiv arXiv

  4. [4]

    Khan, Mehtab Ahmed, Allaerts, Dries, Watson, Simon J. & Churchfield, Matthew J.2025 Investigating the relationship between simulation parameters and flow variables in simulating atmospheric gravity waves for wind energy applications.Wind Energy Science10(6), 1167–1185

    2025

  5. [5]

    Khan, Mehtab Ahmed, Churchfield, Matthew J. & Watson, Simon J.2026 Dependence of wind- farm-induced gravity waves and wind farm performance on non-dimensional atmospheric parameters and simulation configuration.Wind Energy Science11(5), 1631–1652

    2026

  6. [6]

    & Meyers, J.2024 A parametric large-eddy simulation study of wind-farm blockage and gravity waves in conventionally neutral boundary layers.Journal of Fluid Mechanics979, A54

    Lanzilao, L. & Meyers, J.2024 A parametric large-eddy simulation study of wind-farm blockage and gravity waves in conventionally neutral boundary layers.Journal of Fluid Mechanics979, A54

    2024

  7. [7]

    & Meneveau, Charles 2025 JHTDB-wind: a web-accessible large-eddy simulation database of a wind farm with virtual sensor querying.Wind Energy Science10(12), 2821–2840

    Lemson, Gerard, Yao, Hanxun, Szalay, Alexander S., Gayme, Dennice F. & Meneveau, Charles 2025 JHTDB-wind: a web-accessible large-eddy simulation database of a wind farm with virtual sensor querying.Wind Energy Science10(12), 2821–2840. 0X0-11

    2025