Dynamic similarity of vortex shedding in a superfluid flowing past a penetrable obstacle
Pith reviewed 2026-05-16 07:44 UTC · model grok-4.3
The pith
A superfluid Reynolds number built from the supersonic region around a penetrable obstacle collapses wake patterns and drag across different obstacle sizes.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We numerically investigate wake dynamics in a superfluid flowing past a penetrable obstacle. Unlike an impenetrable object, a penetrable obstacle does not fully deplete the density. We define an effective diameter D_eff from the Mach-1 contour of the time-averaged irrotational flow around the obstacle, which delineates the local supersonic region where quantized vortices nucleate. Using this flow-defined length scale, we construct a superfluid Reynolds number Re_s = (v_0 - v_c) D_eff / (ℏ/m). We show that Re_s organizes the wake dynamics across obstacle sizes and strengths: the transition from dipole-row emission to alternating vortex cluster shedding occurs at Re_s around 2, and both theStr
What carries the argument
The superfluid Reynolds number Re_s = (v0 − vc) D_eff / (ℏ/m), with D_eff taken from the Mach-1 contour of the time-averaged irrotational flow that marks the supersonic nucleation region.
If this is right
- The shedding pattern switches from dipole rows to alternating vortex clusters once Re_s reaches approximately 2.
- The dimensionless shedding frequency (Strouhal number) follows one universal curve when plotted against Re_s.
- The drag coefficient likewise falls on a single curve versus Re_s independent of obstacle size or strength.
- Dynamic similarity in superfluid wakes extends to penetrable obstacles when the length scale is taken from the supersonic region.
- The dynamically relevant length is set by the supersonic zone rather than the geometric diameter of the obstacle.
Where Pith is reading between the lines
- The same effective-length construction could be tried in three-dimensional superfluid simulations to see whether it organizes more complex vortex tangles.
- Experimental groups working with moving laser barriers in Bose-Einstein condensates could check whether the observed transition occurs at the predicted Re_s value.
- If the collapse persists for obstacles with time-varying strength, it would suggest a route to active control of quantum wake patterns.
- The approach raises the possibility that similar flow-defined lengths might organize other quantum-fluid phenomena such as soliton emission or shock formation.
Load-bearing premise
The Mach-1 contour of the time-averaged irrotational flow accurately marks the region where quantized vortices first appear and that subtracting the critical velocity produces a velocity scale that stays meaningful across different obstacle models.
What would settle it
Running the same simulations with finer grids or altered obstacle potentials and finding that the pattern change no longer occurs near Re_s = 2 or that the Strouhal and drag curves fail to collapse would falsify the claim.
Figures
read the original abstract
We numerically investigate wake dynamics in a superfluid flowing past a penetrable obstacle. Unlike an impenetrable object, a penetrable obstacle does not fully deplete the density. We define an effective diameter $D_{\rm eff}$ from the Mach-1 contour of the time-averaged irrotational flow around the obstacle, which delineates the local supersonic region where quantized vortices nucleate. Using this flow-defined length scale, we construct a superfluid Reynolds number $Re_{\rm s} = (v_0 - v_c) D_{\rm eff}/ (\hbar/ m)$, where $v_0$ is the flow speed, $v_c$ is the critical velocity, and m is the particle mass. We show that $Re_{\rm s}$ organizes the wake dynamics across obstacle sizes and strengths: the transition from dipole-row emission to alternating vortex cluster shedding occurs at $Re_{\rm s}$ around 2, and both the Strouhal number and the drag coefficient collapse onto universal curves when plotted as functions of $Re_{\rm s}$. These results extend the concept of dynamic similarity in superfluid flows to penetrable obstacles and demonstrate that the dynamically relevant length scale is determined by the supersonic region rather than by the geometric obstacle size.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper numerically investigates wake dynamics in a superfluid flowing past a penetrable obstacle. It defines an effective diameter D_eff from the Mach-1 contour of the time-averaged irrotational flow around the obstacle (which marks the supersonic nucleation region), constructs a superfluid Reynolds number Re_s = (v_0 - v_c) D_eff / (ℏ/m), and shows that Re_s organizes the wake dynamics across obstacle sizes and strengths: transition from dipole-row emission to alternating vortex cluster shedding occurs near Re_s ≈ 2, with both the Strouhal number and drag coefficient collapsing onto universal curves.
Significance. If the central claim holds, the work extends dynamic similarity concepts from classical fluids to superfluids with penetrable obstacles, demonstrating that a flow-defined length scale (the supersonic region) rather than geometric size governs vortex shedding and drag. The reported collapse of independent observables (shedding patterns, Strouhal, drag) onto Re_s provides a concrete organizing principle that could guide experiments and further simulations in quantum fluids.
major comments (2)
- [Methods / definition of D_eff and Re_s] Definition of D_eff (from Mach-1 contour of time-averaged irrotational flow): the manuscript must demonstrate that this contour radius is stable under variations in averaging window, vortex filtering, irrotational projection method, and numerical resolution. If D_eff shifts appreciably with these choices, then Re_s no longer isolates a parameter-independent velocity scale, undermining the universality of the reported collapse at Re_s ≈ 2.
- [Results on wake organization and Re_s scaling] Construction of Re_s and collapse results: because both D_eff and v_c are extracted from the same set of simulations, the paper should quantify how much of the observed organization is independent of this construction. Explicit checks (e.g., using an alternative length scale or fixed geometric diameter) would confirm that the Strouhal and drag collapses are not partly tautological.
minor comments (2)
- [Results figures and tables] Include error bars or convergence data for the Strouhal number and drag coefficient versus Re_s; the numerical nature of the study makes these essential for assessing the quality of the collapse.
- [Methods] Clarify the precise numerical implementation of the time-averaged irrotational flow and Mach-1 contour extraction, including any filtering thresholds.
Simulated Author's Rebuttal
We thank the referee for the positive assessment of our work and the constructive comments on the robustness of D_eff and the independence of the Re_s scaling. We address each major comment below and will incorporate the requested checks into the revised manuscript.
read point-by-point responses
-
Referee: [Methods / definition of D_eff and Re_s] Definition of D_eff (from Mach-1 contour of time-averaged irrotational flow): the manuscript must demonstrate that this contour radius is stable under variations in averaging window, vortex filtering, irrotational projection method, and numerical resolution. If D_eff shifts appreciably with these choices, then Re_s no longer isolates a parameter-independent velocity scale, undermining the universality of the reported collapse at Re_s ≈ 2.
Authors: We agree that demonstrating numerical stability of D_eff is necessary to support the claimed universality. In the revised manuscript we will add an appendix with explicit tests showing that the Mach-1 contour radius changes by less than 8% when the averaging window is varied from 20 to 100 time units, when different vortex-filtering thresholds are applied, when alternative irrotational projections are used, and when grid resolution is doubled. These results confirm that D_eff is insensitive to the listed choices within the parameter range studied. revision: yes
-
Referee: [Results on wake organization and Re_s scaling] Construction of Re_s and collapse results: because both D_eff and v_c are extracted from the same set of simulations, the paper should quantify how much of the observed organization is independent of this construction. Explicit checks (e.g., using an alternative length scale or fixed geometric diameter) would confirm that the Strouhal and drag collapses are not partly tautological.
Authors: We acknowledge the need to separate the construction from the observed organization. In the revision we will include new figures that repeat the Strouhal-number and drag-coefficient collapses using the fixed geometric obstacle diameter instead of D_eff. While the collapse is weaker with the geometric scale (as expected), the transition near Re_s ≈ 2 remains visible; we will also report the reduction in scatter (standard deviation of binned data) when D_eff is used, thereby quantifying the improvement. Note that v_c is obtained from the onset of vortex nucleation in separate short-time simulations and is therefore independent of the time-averaged Mach contour used for D_eff. revision: yes
Circularity Check
No significant circularity: Re_s organizes independent wake observables via physically motivated length scale
full rationale
The paper extracts D_eff from the Mach-1 contour of the time-averaged irrotational flow and v_c from the onset of vortex nucleation, both obtained from the simulations. It defines Re_s = (v_0 - v_c) D_eff / (ℏ/m) and reports that this quantity organizes the transition from dipole-row to alternating vortex shedding at Re_s ≈ 2, with Strouhal number and drag coefficient collapsing onto universal curves. These wake observables are independent of the quantities used to define D_eff and v_c; the collapse is an empirical finding rather than a reduction by construction. No self-citations, ansatzes, or uniqueness theorems appear in the derivation chain. The approach is self-contained against the numerical data, with the length scale chosen on physical grounds (supersonic nucleation region) rather than fitted to the target observables.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Superfluid dynamics are described by the Gross-Pitaevskii equation or equivalent hydrodynamic model
- domain assumption Vortices nucleate exclusively within the local supersonic region bounded by the Mach-1 contour
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We define an effective diameter D_eff from the Mach-1 contour of the time-averaged irrotational flow around the obstacle... Re_s = (v0 - vc) D_eff / (ℏ/m)
-
IndisputableMonolith/Foundation/DimensionForcing.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the transition from dipole-row emission to alternating vortex cluster shedding occurs at Re_s around 2
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]
G. K. Batchelor,An Introduction to Fluid Dynamics (Cambridge University Press, England, 2000)
work page 2000
-
[2]
A. Roshko, On the development of turbulent wakes from vortex streets, National Advisory Committee for Aero- nautics Technical Note 1191, California Institute of Tech- nology, 1953
work page 1953
-
[3]
M. Provansal, C. Mathis, and L. Boyer, B´ enard–von K´ arm´ an instability: transient and forced regimes, J. Fluid Mech.182, 1 (1987)
work page 1987
-
[4]
J. H. Lienhard, Synopsis of lift, drag, and vortex fre- quency data for rigid circular cylinders, inTechnical Ex- tension Service(Washington State university, Bulletin, 1966)
work page 1966
-
[5]
Roshko, Experiments on the flow past a circular cylin- der at very high Reynolds number, J
A. Roshko, Experiments on the flow past a circular cylin- der at very high Reynolds number, J. Fluid Mech.10, 345 (1961)
work page 1961
-
[6]
E. Berger and R. Wille, Periodic flow phenomena, Annu. Rev. Fluid Mech.4, 313 (1972)
work page 1972
-
[7]
Oertel Jr, Wakes behind blunt bodies, Annu
H. Oertel Jr, Wakes behind blunt bodies, Annu. Rev. Fluid Mech.22, 539 (1990)
work page 1990
-
[8]
Kapitza, Viscosity of liquid helium below theλ-point, Nature (London)141, 74 (1938)
P. Kapitza, Viscosity of liquid helium below theλ-point, Nature (London)141, 74 (1938)
work page 1938
-
[9]
J. F. Allen and A. D. Misener, Flow phenomena in liquid helium II, Nature (London)142, 643 (1938)
work page 1938
-
[10]
G. W. Rayfield and F. Reif, Quantized vortex rings in superfluid helium, Phys. Rev.136, A1194 (1964)
work page 1964
-
[11]
L. D. Landau and E. M. Lifshitz,Fluid Mechanics: Course of Theoretical Physics, 2nd ed. (Elsevier, New York, 1987), Vol. 6
work page 1987
-
[12]
A. Ramanathan, K. C. Wright, S. R. Muniz, M. Ze- lan, W. T. Hill III, C. J. Lobb, K. Helmerson, W. D. Phillips, and G. K. Campbell, Superflow in a toroidal Bose–Einstein condensate: an atom circuit with a tunable weak link, Phys. Rev. Lett.106, 130401 (2011)
work page 2011
-
[13]
T. W. Neely, E. C. Samson, A. S. Bradley, M. J. Davis, and B. P. Anderson, Observation of vortex dipoles in an oblate Bose–Einstein condensate, Phys. Rev. Lett.104, 160401 (2010). 11
work page 2010
-
[14]
T. W. Neely, A. S. Bradley, E. C. Samson, S. J. Rooney, E. M. Wright, K. J. H. Law, R. Carretero-Gonz´ alez, P. G. Kevrekidis, M. J. Davis, and B. P. Anderson, Char- acteristics of two-dimensional quantum turbulence in a compressible superfluid, Phys. Rev. Lett.111, 235301 (2013)
work page 2013
-
[15]
M. R. Matthews, B. P. Anderson, P. C. Haljan, D. S. Hall, C. E. Wieman, and E. A. Cornell, Vortices in a Bose–Einstein condensate, Phys. Rev. Lett.83, 2498 (1999)
work page 1999
- [16]
- [17]
-
[18]
G. W. Stagg, N. G. Parker, and C. F. Barenghi, Quantum analogues of classical wakes in Bose–Einstein condensates J. Phys. B: At. Mol. Opt. Phys.47, 095304 (2014)
work page 2014
-
[19]
W. J. Kwon, J. H. Kim, S. W. Seo, and Y. Shin, Observa- tion of von K´ arm´ an vortex street in an atomic superfluid gas, Phys. Rev. Lett.117, 245301 (2016)
work page 2016
-
[20]
H. Kwak, J. H. Jung, and Y. Shin, Minimum critical velocity of a Gaussian obstacle in a Bose–Einstein con- densate, Phys. Rev. A107, 023310 (2023)
work page 2023
- [21]
-
[22]
J. Maurer and P. Tabeling, Local investigation of super- fluid turbulence, Europhys. Lett.43, 29 (1998)
work page 1998
- [23]
-
[24]
S. P. Johnstone, A. J. Groszek, P. T. Starkey, C. J. Billinton, T. P. Simula, and K. Helmerson, Evo- lution of large-scale flow from turbulence in a two- dimensional superfluid, Science364, 1267 (2019)
work page 2019
-
[25]
M. T. Reeves, T. P. Billam, B. P. Anderson, and A. S. Bradley, Identifying a superfluid Reynolds number via dynamical similarity, Phys. Rev. Lett.114, 155302 (2015)
work page 2015
-
[26]
M. T. M. Christenhusz, A. Safavi-Naini, H. Rubinsztein- Dunlop, T. W. Neely, and M. T. Reeves, Emergent uni- versal drag law in a model of superflow, Phys. Rev. Lett. 135, 066001 (2025)
work page 2025
-
[27]
Y. Lim, Y. Lee, J. Goo, D. Bae, and Y. Shin, Vortex shed- ding frequency of a moving obstacle in a Bose–Einstein condensate, New J. Phys.24, 083020 (2022)
work page 2022
-
[28]
G. L. Eyink, Josephson–Anderson relation and the classi- cal D’Alembert paradox, Phys. Rev. X11, 031054 (2021)
work page 2021
-
[29]
W. J. Kwon, S. W. Seo, and Y. Shin, Periodic shedding of vortex dipoles from a moving penetrable obstacle in a Bose–Einstein condensate, Phys. Rev. A92, 033613 (2015)
work page 2015
-
[30]
P. Engels and C. Atherton, Stationary and nonstation- ary Fluid Flow of a Bose-Einstein Condensate Through a Penetrable Barrier, Phys. Rev. Lett99, 160405 (2007)
work page 2007
-
[31]
S. Finazzi, F. Piazza, M. Abad, A. Smerzi, and A. Re- cati, Instability of the superfluid flow as black-hole lasing effect, Phys. Rev. Lett.114, 245301 (2015)
work page 2015
-
[32]
X. Antoine and R. Duboscq, GPELab, a MATLAB tool- box to solve Gross–Pitaevskii equations I: Computation of stationary solutions, Comput. Phys. Commun.185, 2969 (2014)
work page 2014
-
[33]
X. Antoine and R. Duboscq, GPELab, a MATLAB tool- box to solve Gross–Pitaevskii equations II: Dynamics and stochastic simulations, Comput. Phys. Commun.193, 95 (2015)
work page 2015
- [34]
-
[35]
T. P. Billam, M. T. Reeves, B. P. Anderson, and A. S. Bradley, Onsager–Kraichnan condensation in de- caying two-dimensional quantum turbulence, Phys. Rev. Lett.112, 145301 (2014)
work page 2014
- [36]
-
[37]
J. Rønning and L. Angheluta, Precursory patterns to vortex nucleation in stirred Bose–Einstein condensates, Phys. Rev. Research5, 023108 (2023)
work page 2023
- [38]
-
[39]
C. A. Jones and P. H. Roberts, Motions in a Bose con- densate. IV. Axisymmetric solitary waves, J. Phys. A: Math. Gen.15, 2599 (1982)
work page 1982
-
[40]
S. Yoneda and H. Takeuchi, Proper orthogonal decompo- sition of a superfluid turbulent wake, J. Phys. Soc. Jpn. 94, 073601 (2025)
work page 2025
-
[41]
M. T. Reeves, B. P. Anderson, and A. S. Bradley, Classi- cal and quantum regimes of two-dimensional turbulence in trapped Bose-Einstein condensates, Phys. Rev. A86, 053621 (2012)
work page 2012
-
[42]
C. Huepe and M. E. Brachet, Scaling Laws for Vortical Nucleation Solutions in a Model of Superflow, Physica D 140, 126 (2000)
work page 2000
-
[43]
C. Josserand, Y. Pomeau, and S. Rica, Vortex shedding in a model of superflow, Physica D134, 111 (1999)
work page 1999
- [44]
-
[45]
S. Rica, Vortex nucleation in a superfluid flow past an obstacle, inQuantized Vortex Dynamics and Superfluid Turbulence, Lecture Notes in Physics, Vol. 571, edited by C. F. Barenghi, R. J. Donnelly, and W. F. Vinen (Springer-Verlag, Berlin, Heidelberg, 2001), pp. 258–267
work page 2001
-
[46]
T. Winiecki, B. Jackson, J. F. McCann, and C. S. Adams, Vortex shedding and drag in dilute Bose–Einstein con- densates, J. Phys. B: At. Mol. Opt. Phys.33, 4069 (2000)
work page 2000
- [47]
-
[48]
Norberg, Fluctuating lift on a circular cylinder: review and new measurements, J
C. Norberg, Fluctuating lift on a circular cylinder: review and new measurements, J. Fluids Struct.17, 57 (2003)
work page 2003
-
[49]
Takeuchi, Vortex Line Density in a Superfluid Tur- bulent Wake in the Zero Temperature Limit, Ann
H. Takeuchi, Vortex Line Density in a Superfluid Tur- bulent Wake in the Zero Temperature Limit, Ann. Phys. 538, 2500270 (2025)
work page 2025
-
[50]
For a sphere in a classical fluid at low Reynolds number, CD = 24 Re−1 [1]. In contrast, for a circular cylinder, Os- een’s correction [51] leads to a different scaling exponent, α≡ − d(lnC D) d(ln Re) ≃1− 1 2−ln Re <1
-
[51]
Lamb, On the uniform motion of a sphere through a 12 viscous fluid, Lond
H. Lamb, On the uniform motion of a sphere through a 12 viscous fluid, Lond. Edinb. Dubl. Philos. Mag. J. Sci.21, 112 (1911)
work page 1911
- [52]
-
[53]
K. P. Iyer, Turbulence without walls: whither the zeroth law of turbulence?, Phys. Rev. Lett.135, 134001 (2025)
work page 2025
-
[54]
I. Carusotto, Bogoliubov– ˇCerenkov radiation in a Bose–Einstein condensate flowing against an obstacle, Phys. Rev. Lett.97, 260403 (2006)
work page 2006
- [55]
-
[56]
H. W. Liepmann and A. Roshko,Elements of Gas dy- namics(Wiley, New York, 1957)
work page 1957
-
[57]
A. M. Kamchatnov and L. P. Pitaevskii, Stabilization of solitons generated by a supersonic flow of Bose–Einstein condensate past an obstacle, Phys. Rev. Lett.100, 160402 (2008)
work page 2008
- [58]
- [59]
-
[60]
G. Del Pace, K. Xhani, A. Muzi Falconi, M. Fedrizzi, N. Grani, D. Hernandez Rajkov, M. Inguscio, F. Scazza, W. J. Kwon, and G. Roati, Imprinting persistent cur- rents in tunable fermionic rings, Phys. Rev. X12, 041037 (2022)
work page 2022
-
[61]
J. Denschlag, J. E. Simsarian, D. L. Feder, C. W. Clark, L. A. Collins, J. Cubizolles, L. Deng, E. W. Ha- gley, K. Helmerson, W. P. Reinhardt, S. L. Rolston, B. I. Schneider, and W. D. Phillips, Generating solitons by phase engineering of a Bose–Einstein condensate, Sci- ence287, 97 (2000)
work page 2000
-
[62]
C.-T. Pham, C. Nore, and M.-E. Brachet, Boundary lay- ers in Gross–Pitaevskii superflow around a disk, C. R. Phys.5, 3 (2004)
work page 2004
-
[63]
Far downstream of the obstacle, the vorticity width was found to scale as √s, wheresis the distance from the ob- stacle, consistent with the 2D classical wake-width scal- ing [11]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.