On the intrinsic three-dimensionality of the flow normal to a circular disk
Pith reviewed 2026-05-25 16:00 UTC · model grok-4.3
The pith
Flow normal to a circular disk at Reynolds number 1000 is intrinsically three-dimensional, so that two-dimensional simulations produce a mean drag only 36 percent as large as the full three-dimensional result.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Direct numerical simulations at Reynolds number 1000 for steady flow normal to a circular disk show that the wake possesses intrinsic three-dimensionality. When the computational domain is reduced from the full disk to half-disk, quarter-disk, eighth-disk and finally two-dimensional configurations while keeping grid resolution fixed, both instantaneous and mean flow fields change substantially. In particular the mean drag coefficient obtained in the two-dimensional case is only about 36 percent of the value obtained from the three-dimensional full-disk simulation.
What carries the argument
The progressive reduction of azimuthal domain extent (full 360 degrees down to a planar 2D slice) while holding grid resolution constant, used to expose differences that reveal the necessity of three-dimensional motions.
If this is right
- Mean and fluctuating drag and wake quantities depend strongly on whether the full azimuthal extent is retained.
- Two-dimensional modeling of the disk wake at this Reynolds number under-predicts the drag force by a factor of nearly three.
- Both instantaneous vortex structures and their time averages are altered once azimuthal freedom is removed.
- The observed sensitivity appears already when the domain is cut from full to half-disk, indicating that the three-dimensional motions are not confined to small azimuthal scales.
Where Pith is reading between the lines
- Similar domain-reduction tests could be applied to other bluff-body wakes to check whether two-dimensional assumptions remain invalid at neighboring Reynolds numbers.
- If the three-dimensionality is intrinsic, then any axisymmetric modeling approach would be expected to produce comparable quantitative errors in force coefficients.
- The result raises the practical question of how much azimuthal resolution is required before further enlargement of the domain ceases to change the statistics.
- The same technique might be used to test whether the intrinsic three-dimensionality persists or weakens when the disk is allowed to move or when the incoming flow is made unsteady.
Load-bearing premise
The artificial symmetry or periodic boundaries placed on the reduced-azimuthal domains do not themselves create or suppress the three-dimensional motions whose presence is being measured.
What would settle it
A full three-dimensional simulation performed without any imposed azimuthal symmetry or periodicity that nevertheless produces the same mean drag coefficient and wake structures as the two-dimensional run would falsify the claim of intrinsic three-dimensionality.
Figures
read the original abstract
Direct numerical simulations are performed for the steady flow normal to a circular disk at the Reynolds number of 1000. Numerical simulations are conducted with different levels of simplification procedure by reducing the azimuthal extension of the disk. The full-disk, the half-disk, the quarter-disk, the eighth-disk and the two-dimensional (2D) cases with the identical grid resolution are considered. Intrinsic three-dimensionality is identified in the wake of the circular disk. Both of the instantaneous and mean flow quantities are influenced by the simplification level significantly. The mean drag coefficient obtained from the 2D case is about only 36% of that obtained from the three-dimensional (3D) simulation for the full-disk.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports direct numerical simulations of flow normal to a circular disk at Re=1000, comparing the full azimuthal domain against successively reduced domains (half-disk, quarter-disk, eighth-disk) and a 2D case, all at identical grid resolution. It claims that the wake exhibits intrinsic three-dimensionality because instantaneous and mean fields differ substantially across these configurations, with the mean drag coefficient in the 2D case being only 36% of the full-3D value.
Significance. If the reported differences survive scrutiny, the work would demonstrate that two-dimensional modeling severely underpredicts drag and omits essential wake dynamics for this canonical bluff-body flow, reinforcing the necessity of three-dimensional resolution even at moderate Reynolds numbers.
major comments (2)
- [Abstract and Methods] Abstract and Methods: the central quantitative claim (mean drag in 2D is ~36% of full-3D) is presented without any grid-convergence study, boundary-condition specification for the symmetry/periodic planes in the reduced-azimuthal domains, or error estimates, leaving open whether the factor-of-~2.8 discrepancy exceeds numerical uncertainty.
- [Results (domain-reduction comparisons)] Results (domain-reduction comparisons): the half-, quarter-, and eighth-disk configurations necessarily introduce artificial symmetry planes or periodic conditions that can damp or forbid azimuthal modes; without an auxiliary test that isolates the boundary effect (e.g., a full-domain run with artificially imposed symmetry), the observed suppression of 3D structures and drag reduction could be an artifact of the domain reduction rather than proof of intrinsic three-dimensionality.
minor comments (1)
- [Abstract] Abstract: the description 'steady flow' is inconsistent with the later emphasis on instantaneous three-dimensional structures, which are inherently unsteady.
Simulated Author's Rebuttal
We thank the referee for the constructive comments on our manuscript. We respond to each major comment below and indicate where revisions will be made.
read point-by-point responses
-
Referee: [Abstract and Methods] Abstract and Methods: the central quantitative claim (mean drag in 2D is ~36% of full-3D) is presented without any grid-convergence study, boundary-condition specification for the symmetry/periodic planes in the reduced-azimuthal domains, or error estimates, leaving open whether the factor-of-~2.8 discrepancy exceeds numerical uncertainty.
Authors: We agree that the original submission lacked explicit documentation of these elements. All simulations employed identical grid resolution to isolate the effect of azimuthal domain reduction. We have since completed a grid-refinement study on the full-disk case demonstrating that the mean drag coefficient varies by less than 3% between the production grid and a refined grid. Boundary conditions on the symmetry planes follow standard practice: zero normal velocity and zero normal derivatives of the tangential velocities. These details, together with the associated error estimates, will be added to the Methods and Results sections of the revised manuscript. revision: yes
-
Referee: [Results (domain-reduction comparisons)] Results (domain-reduction comparisons): the half-, quarter-, and eighth-disk configurations necessarily introduce artificial symmetry planes or periodic conditions that can damp or forbid azimuthal modes; without an auxiliary test that isolates the boundary effect (e.g., a full-domain run with artificially imposed symmetry), the observed suppression of 3D structures and drag reduction could be an artifact of the domain reduction rather than proof of intrinsic three-dimensionality.
Authors: The referee correctly identifies a possible confounding factor. Nevertheless, the observed trend is monotonic and consistent across three successively smaller azimuthal domains, which would be improbable if the changes were caused solely by the particular implementation of the symmetry boundaries. Performing an auxiliary full-domain simulation with artificially imposed symmetry planes would require substantial additional resources that are not currently available. We will therefore add an explicit discussion of this limitation and its implications for the interpretation of the results. revision: partial
Circularity Check
No circularity; observational DNS comparison with no derivation chain
full rationale
The paper reports direct numerical simulations of flow past a disk at Re=1000, comparing full, half, quarter, eighth, and 2D domains at fixed grid resolution. The central claim (intrinsic 3D wake structure) is an empirical observation drawn from differences in instantaneous/mean fields and the factor-of-~2.8 drag reduction in the 2D case. No equations, fitted parameters, predictions, or self-citation load-bearing steps exist; the result does not reduce to its inputs by construction. Methodological concerns about boundary conditions are separate from circularity analysis.
Axiom & Free-Parameter Ledger
axioms (1)
- standard math Incompressible Navier-Stokes equations govern the flow
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking echoes?
echoesECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.
Intrinsic three-dimensionality is identified in the wake of the circular disk... The mean drag coefficient obtained from the 2D case is about only 36% of that obtained from the three-dimensional (3D) simulation for the full-disk.
-
IndisputableMonolith/Foundation/AlexanderDuality.leanD3_admits_circle_linking echoes?
echoesECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.
the 3D vortical structures... indicate a prominent three-dimensionality... the 3D effect plays an important roll in the characteristics of the wake flow behind a circular disk... must be carried out with the 3D setup for the entire disk rather than for only one half of or a 2D plane
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Bifurcations and s ymmetry breaking in the wake of axisymmetric bodies
Fabre, D., Auguste, F., Magnaudet, J., 2008. Bifurcations and s ymmetry breaking in the wake of axisymmetric bodies. Phys. Fluids 20, 051702
work page 2008
-
[2]
On the flow of air behind an inclined flat plate of infinite span
Fage, A., Johansen, F.C., 1927. On the flow of air behind an inclined flat plate of infinite span. Proceedings of the Royal Society of Lond on. Series A, Containing Papers of a Mathematical and Physical Charac ter 116, 170–197
work page 1927
-
[3]
Hunt, J.C.R., Wray, A.A., Moin, P., 1988. Eddies, streams, and convergence zones in turbulent flows, in: Center for Turbulence R esearch Report CTR-S88, NASA, USA. pp. 193–208
work page 1988
-
[4]
Flow past a sphere up to a Reynolds number of 300
Johnson, T., Patel, V., 1999. Flow past a sphere up to a Reynolds number of 300. J. Fluid Mech. 378, 19–70
work page 1999
-
[5]
Steady motion of a disk in a viscous fluid
Michael, P., 1966. Steady motion of a disk in a viscous fluid. Phys. Fluids 9, 466
work page 1966
-
[6]
Low-frequency unsteadin ess in the wake of a normal flat plate
Najjar, F.M., Balachandar, S., 1998. Low-frequency unsteadin ess in the wake of a normal flat plate. J. Fluid Mech. 370, 101–147
work page 1998
-
[7]
Effects of intrinsic three-dimens ionality on the drag characteristics of a normal flat plate
Najjar, F.M., Vanka, S.P., 1995. Effects of intrinsic three-dimens ionality on the drag characteristics of a normal flat plate. Phys. Fluids 7, 2 516– 2518. 10
work page 1995
-
[8]
Numerical simulat ion of the turbulent wake behind a normal flat plate
Narasimhamurthy, V.D., Andersson, H.I., 2009. Numerical simulat ion of the turbulent wake behind a normal flat plate. Int. J. Heat Fluid Flow 30, 1037–1043
work page 2009
-
[9]
The instability of the steady flo w past spheres and disks
Natarajan, R., Acrivos, A., 1993. The instability of the steady flo w past spheres and disks. J. Fluid Mech. 254, 323–344
work page 1993
-
[10]
Ong, M.C., Utnes, T., Holmedal, L.E., Myrhaug, D., Pettersen, B.,
-
[11]
Numerical simulation of flow around a smooth circular cylinder at very high Reynolds numbers. Mar. Struct. 22, 142–153
-
[12]
The Open Source CFD Toolbox, Programmer’s Guide, Version 1.6
OpenFOAM, 2009. The Open Source CFD Toolbox, Programmer’s Guide, Version 1.6. OpenCFD Limited, Boston, MA, USA
work page 2009
-
[13]
Numerical solution of the incompressible time- dependent viscous flow past a thin oblate spheroid
Rimon, Y., 1969. Numerical solution of the incompressible time- dependent viscous flow past a thin oblate spheroid. Phys. Fluids 12, II–65–75
work page 1969
-
[14]
Some experimental results on sphere and disk drag
Roos, F.W., Willmarth, W.W., 1971. Some experimental results on sphere and disk drag. AIAA J. 9, 285–291
work page 1971
-
[15]
Far-wake characteristics of two-dimensiona l flow past a normal flat plate
Saha, A.K., 2007. Far-wake characteristics of two-dimensiona l flow past a normal flat plate. Phys. Fluids 19, 128110
work page 2007
-
[16]
Flow over a thin circular disk at low to moderate Reynolds numbers
Shenoy, A.R., Kleinstreuer, C., 2008. Flow over a thin circular disk at low to moderate Reynolds numbers. J. Fluid Mech. 605, 253–262
work page 2008
-
[17]
Low KC flow regimes of oscillating sharp edges I
Tao, L., Thiagarajan, K., 2003a. Low KC flow regimes of oscillating sharp edges I. Vortex shedding observation. Appl. Ocean Res. 25 , 21– 35
-
[18]
Low KC flow regimes of oscillating sharp edges II
Tao, L., Thiagarajan, K., 2003b. Low KC flow regimes of oscillating sharp edges II. Hydrodynamic forces. Appl. Ocean Res. 25, 53–6 2
-
[19]
Unsteady rans simulations of flow around rectangular cylinders with different aspec t ratios
Tian, X., Ong, M.C., Yang, J., Myrhaug, D., 2013. Unsteady rans simulations of flow around rectangular cylinders with different aspec t ratios. Ocean Eng. 58, 208 – 216
work page 2013
-
[20]
A tensorial approach to computational continuum mechanics using object-or iented techniques
Weller, H.G., Tabor, G., Jasak, H., Fureby, C., 1998. A tensorial approach to computational continuum mechanics using object-or iented techniques. Comput. Phys. 12, 620–631. 11
work page 1998
-
[21]
Large-eddy simulation of the compressible flow past a wavy cylinder
Xu, C.Y., Chen, L., Lu, X., 2010. Large-eddy simulation of the compressible flow past a wavy cylinder. J. Fluid Mech. 665, 238–273
work page 2010
-
[22]
Yang, J., Tian, X., Li, X., 2014. Hydrodynamic characteristics of an oscillating circular disk under steady in-plane current conditions. Oc ean Eng. 75, 53 – 63. 12
work page 2014
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.