pith. sign in

arxiv: 2606.07625 · v2 · pith:EBENIS6Lnew · submitted 2026-05-30 · ⚛️ physics.ao-ph · physics.flu-dyn

Building drag and shielding in a realistic urban environment

Jingzi Huang , Omduth Coceal , Marco Placidi , Zheng-Tong Xie , Maarten van Reeuwijk This is my paper

Pith reviewed 2026-06-28 17:49 UTC · model grok-4.3

classification ⚛️ physics.ao-ph physics.flu-dyn
keywords urban canopy dragwind directionbuilding shieldinglarge-eddy simulationdrag coefficientfrontal areadirectional anisotropycampus geometry
0
0 comments X

The pith

Accounting for shielding by upstream buildings reduces directional variation in urban canopy drag.

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

The study uses 24 large-eddy simulations of flow over the University of Bristol campus with 110 buildings to examine how wind direction changes total drag. Individual building drag varies strongly because many structures sit in the wakes of others, while the campus-wide drag coefficient fluctuates only moderately. Two ratios, upstream fetch Ls/Hs and height Hs/H, classify buildings into four shielding regimes; near-wake shielded buildings contribute almost no drag. A modified drag coefficient that partially or fully excludes these shielded buildings lowers directional anisotropy and produces an effective frontal area that stays roughly constant regardless of wind direction.

Core claim

Classifying buildings by the ratios Ls/Hs and Hs/H into four regimes and then computing a modified drag coefficient by excluding shielded buildings reduces directional anisotropy and yields an effective frontal area more consistent across wind directions.

What carries the argument

Shielding regimes defined by upstream fetch ratio Ls/Hs and relative height ratio Hs/H with thresholds 5 and 1, used to exclude or down-weight shielded buildings when calculating the campus drag coefficient.

If this is right

  • 20 percent of buildings account for 80 percent of total campus drag.
  • Buildings in the near-wake shielded regime experience negligible drag.
  • Buildings in the far-wake non-shielded regime experience the highest drag.
  • The campus-wide drag coefficient fluctuates only moderately with wind direction.
  • The modified drag coefficient produces a more consistent effective frontal area across directions.

Where Pith is reading between the lines

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

  • The 80/20 contribution split implies that urban models could focus computational effort on the few dominant buildings.
  • Similar shielding corrections might be tested in other dense urban areas to see if they improve wind-load predictions.
  • Climate models that currently use direction-independent drag might gain accuracy by incorporating a simple shielding filter.

Load-bearing premise

The thresholds Ls/Hs = 5 and Hs/H = 1 correctly identify shielding regimes and apply beyond this specific campus layout.

What would settle it

Repeat the simulations on a second campus geometry or with altered thresholds and check whether the modified drag coefficient still shows markedly lower directional variation than the standard one.

Figures

Figures reproduced from arXiv: 2606.07625 by Jingzi Huang, Maarten van Reeuwijk, Marco Placidi, Omduth Coceal, Zheng-Tong Xie.

Figure 1
Figure 1. Figure 1: Definition sketch (a) of solid-fluid interface 𝜕Ω𝑓 and 3-D normal vectors 𝑵 of interface. The solid domain is in grey, while the white indicates the fluid domain Ω𝑓 . The surface of an individual building 𝑞 is outlined in red; (b) of 𝑵 decomposition, the vector 𝑵 is unit and points into the fluid domain. but also allows the evaluation of the drag acting on an individual building by specifying the integrati… view at source ↗
Figure 2
Figure 2. Figure 2: (a) A satellite plane view of the campus of the University of Bristol, overlaid with a footprint of the simulation morphology. From Google Maps. (b) A 3-D view of the simulated campus morphology (Bi et al., 2025), with the colour indicating the building height level. A global Cartesian coordinate system (𝑋, 𝑌 , 𝑍) is defined with the origin at the campus centre of the plane view, where positive 𝑋, 𝑌 are al… view at source ↗
Figure 3
Figure 3. Figure 3: The vertical profiles of plane-averaged velocity (a, e), total kinematic stress (b, f) and disperisive stress (c, g) in streamwise direction (upper row) and spanwise (bottom row) direction, respectively. The streamwise distributed drag 𝑓𝐷 (d). The solid line represents the overall mean value across the wind directions, while the shaded band denotes the range between the minimum and maximum values. The dash… view at source ↗
Figure 4
Figure 4. Figure 4: (a) The frontal area index 𝜆𝑓 and the drag coefficient 𝐶𝑑 of the entire morphology varying with the wind direction 𝜃, where the solid and dashed vertical lines label the principal direction and secondary direction of the morphology, respectively. (b) A plane view of the morphology with the principal direction (thick solid) and secondary direction (thick dashed) from the centre and marked by solid line and … view at source ↗
Figure 5
Figure 5. Figure 5: (a) Cumulative distribution function (CDF) of direction-averaged building drag, where buildings are ranked in descending order of drag. (b) Highlight of the top 20 buildings with the highest drag, where colour indicates their proportion in the total drag of all buildings. The top 3 buildings are labelled with names. 4.2. Direction-averaged drag on individual buildings To investigate the drag acting on indi… view at source ↗
Figure 6
Figure 6. Figure 6: (a) Sketch of the open space (the purple patch) and shielding height in front of a target T-shaped building 42 under the wind direction 𝜃 = 180° shown as the purple arrow. The double-head arrow marks the maximum width of the fetch perpendicular to the wind direction. Variation of the shielding effect parameters (b) 𝐿𝑠∕𝐻𝑠 and (c) 𝐻𝑠∕𝐻 of the target building with the wind direction. The (d) cumulative integr… view at source ↗
Figure 7
Figure 7. Figure 7: (a) Dependence of the averaged drag coefficient of an individual building on the shielding parameters 𝐿𝑠∕𝐻𝑠 and 𝐻𝑠∕𝐻. (b) The probability of the shielding parameters 𝐿𝑠∕𝐻𝑠 and 𝐻𝑠∕𝐻 occurring on the campus. The red dashed lines classify the domain into four regimes (S1 -S4) according to the wake location and shielding situation. Near wake (NW) Far wake (FW) Percentage Average 𝐶𝑑 Percentage Average 𝐶𝑑 No Shi… view at source ↗
Figure 8
Figure 8. Figure 8: (a) Example of the modified drag procedure at 𝜃 = 180° (wind direction indicated by the arrow), showing buildings excluded from the calculation. All excluded buildings lie within the near-wake shielding regime (S1) and are greyed out in this case. A highlight of the modified frontal area in regime S3 is presented in the right corner. (b-d) Comparison of the modified drag coefficient, total drag stress, and… view at source ↗
Figure 9
Figure 9. Figure 9: (a) Facet mesh of an example building, with a representative facet highlighted together with its associated nearest Cartesian grid cells. (b) The facet pressure of building 42 on the Bristol campus. (c) The drag density field 𝜌𝐷 converted from the facet data of building 42. straightforward to show that ∑ 𝑖,𝑗,𝑘 𝜌𝜙;𝑖𝑗𝑘Δ𝑥Δ𝑦Δ𝑧𝑘 = ∑ 𝑚 𝜙𝑚𝐴𝑚, (23) The left-hand side is a discrete volume integral, which we can rew… view at source ↗
Figure 10
Figure 10. Figure 10: The height-masking footprint ℎ of building 42, overlaid with its principal (𝒗1 in black) and secondary (𝒗2 in grey) directions, originating at the building centroid (indicated by a circle). The dashed arrow indicates the wind direction 𝜃, while 𝛼 denotes the angle between the wind direction and 𝒗1 . represents the contribution of facet 𝑚 to the frontal area. Denoting the associated volumetric frontal area… view at source ↗
read the original abstract

Shielding by upstream buildings is a fundamental control on urban drag, yet its influence remains poorly quantified in realistic urban environments. Here, we investigate shielding effects using building-resolved large-eddy simulations of the University of Bristol campus, comprising 110 buildings of varying height, shape and orientation. Twenty-four wind directions are considered, allowing each building to experience a wide range of upstream shielding conditions. While the total drag of the campus exhibits only moderate directional variability, the drag acting on individual buildings varies substantially. In the present case, approximately $20\%$ of buildings account for $80\%$ of the total drag, which is primarily attributed to a small number of large buildings that contribute disproportionately high drag forces. To quantify shielding, we introduce two dimensionless parameters: the upstream fetch ratio, $L_s/H_s$, and the relative height ratio, $H_s/H$, where $L_s$ is the distance to the nearest upstream obstacle, $H_s$ is the height of the upstream obstacle, and $H$ is the height of the target building. These parameters distinguish between near- and far-wake conditions and between sheltered and exposed buildings, providing a simple method to characterise shielding effects in realistic urban environments. The study provides valuable quantitative insight into drag and shielding in the Bristol campus morphology; more importantly, it establishes a general framework for analysing drag and shielding that can be applied in other complex urban environments. The results identify shielding as a primary control on building drag and motivate shielding-aware measures of effective frontal area and drag coefficient

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

1 major / 1 minor

Summary. The manuscript reports results from 24 building-resolved large-eddy simulations of the University of Bristol campus (110 buildings) under constant pressure gradient forcing. It finds that the campus-scale drag coefficient shows moderate directional variation while individual-building drag varies strongly due to shielding; 20% of buildings account for ~80% of total drag. Two dimensionless ratios (upstream fetch L_s/H_s and relative height H_s/H) are introduced with fixed thresholds of 5 and 1 to classify buildings into four shielding regimes. A modified drag coefficient obtained by partially or fully excluding buildings in the near-wake shielded regime is shown to reduce directional anisotropy and produce a more directionally consistent effective frontal area.

Significance. If the reduction in anisotropy holds under the proposed modification, the work supplies a concrete, simulation-derived route to improve urban-canopy drag parametrizations that currently ignore directional shielding. The use of 24 independent, building-resolved LES runs under identical forcing supplies a reproducible dataset that directly supports the central claim without circular fitting.

major comments (1)
  1. [Abstract and shielding-regime definition] Abstract and the paragraph defining the shielding regimes: the thresholds L_s/H_s = 5 and H_s/H = 1 are stated without derivation, calibration against the LES drag values, or any sensitivity test. Because the modified drag coefficient is obtained precisely by excluding buildings classified as near-wake shielded under these exact cut-offs, the reported reduction in directional anisotropy and the improved consistency of effective frontal area rest on an untested choice; modest shifts in either threshold could reclassify enough of the high-drag buildings to alter the quantitative improvement.
minor comments (1)
  1. [Methods] The notation for L_s, H_s and the four regime labels should be accompanied by an explicit equation or schematic in the methods section to avoid ambiguity when readers attempt to reproduce the classification.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their constructive assessment and for highlighting the importance of justifying the regime thresholds. We respond to the single major comment below and will revise the manuscript to address the concern.

read point-by-point responses

  1. Referee: [Abstract and shielding-regime definition] Abstract and the paragraph defining the shielding regimes: the thresholds L_s/H_s = 5 and H_s/H = 1 are stated without derivation, calibration against the LES drag values, or any sensitivity test. Because the modified drag coefficient is obtained precisely by excluding buildings classified as near-wake shielded under these exact cut-offs, the reported reduction in directional anisotropy and the improved consistency of effective frontal area rest on an untested choice; modest shifts in either threshold could reclassify enough of the high-drag buildings to alter the quantitative improvement.

    Authors: We agree that the thresholds require explicit justification and testing. The values were chosen to reflect typical near-wake recovery distances (L_s/H_s ≈ 5) and height ratios (H_s/H = 1) drawn from prior urban canopy studies on wake extent, but the original manuscript does not derive them from the present LES data or test sensitivity. In revision we will add a dedicated paragraph in the methods section explaining the physical basis and include a new sensitivity analysis (varying each threshold by ±20 % and recomputing the modified drag coefficient and anisotropy metrics). This will demonstrate robustness of the reported reduction in directional variability. The revised abstract will briefly note the physical motivation for the cut-offs. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results computed directly from LES outputs

full rationale

The paper performs 24 building-resolved LES runs under fixed pressure gradient and computes drag coefficients and frontal areas directly from the velocity and pressure fields. The two ratios L_s/H_s and H_s/H are introduced as classification criteria with fixed numerical thresholds; buildings are then binned and a modified C_d is obtained by excluding or down-weighting the near-wake-shielded subset. This is a post-processing choice applied to the simulation data, not a fitted parameter whose value is adjusted to force the reported reduction in directional anisotropy, nor a self-citation chain, nor a renaming of a known result. No equation in the supplied text equates the final modified drag to any input quantity by algebraic identity. The derivation chain therefore remains independent of its own outputs.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 1 invented entities

The two thresholds are chosen values that define the regimes; the simulations rest on standard LES assumptions and the constant-pressure-gradient driving condition.

free parameters (2)
  • Ls/Hs threshold = 5
    Chosen cutoff to separate near-wake shielded regime from far-wake non-shielded regime
  • Hs/H threshold = 1
    Chosen cutoff to separate buildings taller or shorter than upstream structures
axioms (2)
  • domain assumption Large-eddy simulations with building-resolved grids accurately capture the drag forces exerted by urban structures
    All quantitative results rest on the 24 LES runs
  • domain assumption A constant imposed pressure gradient produces a representative mean flow equivalent to real atmospheric forcing
    Used to drive all simulations
invented entities (1)
  • Four shielding regimes (near-wake shielded, far-wake non-shielded, etc.) no independent evidence
    purpose: Categorize buildings according to expected drag levels under different wind directions
    Defined by applying the two ratios and chosen thresholds to the campus geometry

pith-pipeline@v0.9.1-grok · 5750 in / 1525 out tokens · 33241 ms · 2026-06-28T17:49:04.599024+00:00 · methodology

discussion (0)

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