Numerical Investigation of Discontinuous Ice Effects on Swept Wings
Pith reviewed 2026-05-18 04:59 UTC · model grok-4.3
The pith
Discontinuous ice on swept wings reduces lift more than continuous ice by disrupting vortices with gap jets.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Discontinuous ice causes a more severe reduction in lift than continuous ice. Continuous ice forms a large separation bubble that helps maintain lift, while discontinuous ice disrupts leading-edge vortex formation through gap jets, resulting in greater lift loss but a smaller drag penalty. The discontinuous-ice wing does not exhibit a sudden stall-induced lift drop. Its flow is characterized by a separating shear layer and Kármán vortex shedding that becomes irregular due to gap jets, with three chord-based Strouhal numbers identified at 11.3, 22.6, and 33.9.
What carries the argument
Gap jets formed between segments of discontinuous ice that interfere with leading-edge vortex formation and create irregular shear layers.
If this is right
- Discontinuous ice produces greater lift loss but smaller drag penalty than continuous ice at the same conditions.
- The discontinuous-ice configuration avoids the sudden lift drop associated with stall on continuous-ice wings.
- Lift and drag fluctuations occur mainly at twice the vortex-shedding frequency because of the gap jets.
- The separating shear layer becomes irregular specifically due to interference from the gap jets.
Where Pith is reading between the lines
- Partial ice removal or de-icing that leaves gaps could sometimes worsen lift performance compared with leaving the ice in place.
- The higher Strouhal number when scaled by ice width suggests gap jets create a more energetic wake than a simple cylinder.
- Design of ice-protection systems might need to account for gap-induced unsteady loads in addition to mean aerodynamic penalties.
Load-bearing premise
The artificially simulated discontinuous ice shapes and the enhanced delayed detached-eddy simulation accurately capture real ice accretion and the unsteady three-dimensional flows on swept wings.
What would settle it
Wind-tunnel measurements of lift, drag, and surface pressures on a physical swept-wing model with naturally accreted or precisely replicated discontinuous ice at the same angles of attack would directly test the reported lift reductions and Strouhal numbers.
read the original abstract
This study investigates the aerodynamic performance and flow structures of infinite swept wings with artificially simulated discontinuous ice using an enhanced delayed detached-eddy simulation. Comparisons are made among clean, continuous-ice, and discontinuous-ice configurations. Results show that discontinuous ice causes a more severe reduction in lift than continuous ice. While continuous ice forms a large separation bubble that helps maintain lift, discontinuous ice disrupts leading-edge vortex formation through gap jets, resulting in greater lift loss but a smaller drag penalty. Unlike the continuous-ice wing, the discontinuous-ice case does not exhibit a sudden stall-induced lift drop. The flow over the discontinuous-ice wing can be characterized by two canonical patterns: a separating shear layer and K\'arm\'an vortex shedding. However, the separating shear layer becomes irregular due to the interference of gap jets. Three characteristic chord-based Strouhal numbers (St)-11.3, 22.6, and 33.9-are identified. The lowest (St=11.3) corresponds to the shedding of vortex pairs; when nondimensionalized by the ice width, it yields St = 0.58, which is higher than that of a canonical cylinder wake. Furthermore, lift and drag fluctuations occur predominantly at St = 22.6, twice the shedding frequency, primarily induced by the gap jets-a phenomenon absent in the continuous-ice case.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript numerically investigates the aerodynamic performance and flow structures of infinite swept wings with artificially simulated discontinuous ice using enhanced delayed detached-eddy simulation (ED-DES). It compares clean, continuous-ice, and discontinuous-ice configurations and claims that discontinuous ice causes a more severe reduction in lift than continuous ice. Continuous ice forms a large separation bubble that helps maintain lift, while discontinuous ice disrupts leading-edge vortex formation through gap jets, yielding greater lift loss but a smaller drag penalty and no sudden stall-induced lift drop. The flow over the discontinuous-ice wing exhibits separating shear layers and Kármán vortex shedding, with three characteristic chord-based Strouhal numbers (11.3, 22.6, 33.9) identified; lift and drag fluctuations occur predominantly at St = 22.6 due to gap jets.
Significance. If the ED-DES predictions prove accurate, the results would advance understanding of three-dimensional unsteady interactions on iced swept wings by contrasting the effects of continuous versus discontinuous ice accretion. The detailed characterization of gap-jet interference with vortex formation, the absence of sudden stall in the discontinuous case, and the specific Strouhal numbers (including the ice-width-normalized value of 0.58) provide concrete, falsifiable observations on force fluctuations and flow patterns that could inform icing-related aerodynamic modeling.
major comments (3)
- [Abstract] Abstract: The central comparative claim that discontinuous ice produces more severe lift reduction than continuous ice (via gap-jet disruption of leading-edge vortices versus a stabilizing separation bubble) is not anchored by any reported validation against icing-tunnel or flight-test data at the relevant sweep and Reynolds number; without such anchoring the lift-ordering result remains conditional on the specific numerical idealization.
- [Methods/Results] Methods/Results: No sensitivity study is presented on the artificial discontinuous-ice parameters (gap width, protrusion height, or spanwise periodicity) that control the gap-jet momentum and vortex coherence; because these parameters are imposed rather than grown from an icing model, their variation could alter the reported lift penalty and the identified flow patterns.
- [Results] Results: Mesh-convergence data and uncertainty quantification for the lift, drag, and Strouhal-number predictions are absent, which directly affects confidence in the quantitative differences in stall behavior and the claim that fluctuations occur at twice the shedding frequency due to gap jets.
minor comments (2)
- [Abstract] Abstract: The definition of the Strouhal number when nondimensionalized by ice width (yielding St = 0.58) should be stated explicitly and compared to literature values for cylinder wakes to strengthen the interpretation.
- Ensure consistent first-use definition of acronyms (ED-DES, DDES) and clarify whether the reported Kármán vortex shedding is spanwise or streamwise in the discontinuous-ice case.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed comments on our numerical study of discontinuous ice on swept wings. We address each major comment point by point below, with honest indications of where revisions will be incorporated.
read point-by-point responses
-
Referee: [Abstract] Abstract: The central comparative claim that discontinuous ice produces more severe lift reduction than continuous ice (via gap-jet disruption of leading-edge vortices versus a stabilizing separation bubble) is not anchored by any reported validation against icing-tunnel or flight-test data at the relevant sweep and Reynolds number; without such anchoring the lift-ordering result remains conditional on the specific numerical idealization.
Authors: Our investigation is a high-fidelity numerical study employing ED-DES on infinite swept wings. The reported lift reduction ordering and associated flow mechanisms (gap jets disrupting leading-edge vortices in the discontinuous case versus a stabilizing separation bubble in the continuous case) are direct outcomes of the simulated flow fields under identical numerical conditions. We agree that the quantitative results are conditional on the numerical idealizations and lack direct experimental anchoring at the exact sweep and Reynolds number. In the revised manuscript we will explicitly qualify the abstract and discussion sections to emphasize the computational nature of the findings and the physical mechanisms identified, while noting the desirability of future experimental validation. revision: partial
-
Referee: [Methods/Results] Methods/Results: No sensitivity study is presented on the artificial discontinuous-ice parameters (gap width, protrusion height, or spanwise periodicity) that control the gap-jet momentum and vortex coherence; because these parameters are imposed rather than grown from an icing model, their variation could alter the reported lift penalty and the identified flow patterns.
Authors: The discontinuous-ice geometry parameters were chosen to represent typical values drawn from icing literature for swept-wing accretion, enabling the formation of gap jets while remaining computationally tractable. A comprehensive sensitivity study on gap width, protrusion height, and periodicity was not performed owing to the substantial cost of additional ED-DES runs. We will add a paragraph in the methods section of the revised manuscript that explains the parameter selection rationale and qualitatively discusses expected sensitivities of the gap-jet momentum and vortex coherence based on our preliminary observations. revision: partial
-
Referee: [Results] Results: Mesh-convergence data and uncertainty quantification for the lift, drag, and Strouhal-number predictions are absent, which directly affects confidence in the quantitative differences in stall behavior and the claim that fluctuations occur at twice the shedding frequency due to gap jets.
Authors: Grid-independence checks were conducted during the simulation campaign to confirm adequate resolution of the separating shear layers, gap jets, and vortex shedding. Detailed convergence data were not reported in the original manuscript. We will include a new subsection (or appendix) in the revised version that presents mesh-convergence results for lift and drag coefficients together with the extracted Strouhal numbers, along with a brief statement on numerical uncertainty derived from the observed variations across grids. revision: yes
- Direct experimental validation of the lift-reduction ordering and Strouhal-number values at the specific sweep and Reynolds number, which would require new icing-tunnel or flight-test data beyond the scope of this numerical investigation.
Circularity Check
No significant circularity in direct CFD simulation outputs
full rationale
The paper reports aerodynamic coefficients, flow structures, and Strouhal numbers obtained by solving the unsteady Navier-Stokes equations with an enhanced delayed detached-eddy simulation turbulence closure on fixed, artificially imposed ice geometries. These quantities are direct numerical outputs, not parameters fitted to the simulation data and then re-labeled as predictions, nor self-defined quantities. No load-bearing self-citations, uniqueness theorems, or ansatzes imported from prior author work appear in the derivation chain. The central comparative claims (discontinuous ice producing greater lift loss via gap-jet disruption) follow from the computed solutions without reduction to the inputs by construction.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Enhanced delayed detached-eddy simulation accurately captures the unsteady separated flow and vortex shedding induced by ice geometries on swept wings.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Results show that discontinuous ice causes a more severe reduction in lift than continuous ice. While continuous ice forms a large separation bubble that helps maintain lift, discontinuous ice disrupts leading-edge vortex formation through gap jets...
-
IndisputableMonolith/Foundation/DimensionForcing.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Three characteristic chord-based Strouhal numbers (St)—11.3, 22.6, and 33.9—are identified.
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]
M. B. Bragg, A. P. Broeren, L. A. Blumenthal, "Iced-airfoil aerodynamics," Progress in Aerospace Sciences 41(5) (2005) 323–362. https://doi.org/10.1016/j.paerosci.2005.07.001
-
[2]
Physical mechanisms of glaze ice scallop formations on swept wings,
M. Vargas, E. Reshotko, "Physical mechanisms of glaze ice scallop formations on swept wings," 36th AIAA Aerospace Sciences Meeting and Exhibit (1998)
work page 1998
-
[3]
Current experimental basis for modeling ice accretions on swept wings,
M. Vargas, "Current experimental basis for modeling ice accretions on swept wings," Journal of Aircraft 44(1) (2007) 274–290
work page 2007
-
[4]
Scallop ice shape characteristics of swept wing based on large-scale icing wind tunnel experiment,
W. A. Qiang, C. H. Ningli, W. A. Yuanbo, W. H. Li, Y. Li, X. Yi, "Scallop ice shape characteristics of swept wing based on large-scale icing wind tunnel experiment," Chinese Journal of Aeronautics 36(12) (2023) 214–230
work page 2023
-
[5]
Aerodynamic measurements of an airfoil with simulated glaze ice,
M. B. Bragg, W. J. Coirier, "Aerodynamic measurements of an airfoil with simulated glaze ice," AIAA Paper 86-0484, 24th Aerospace Sciences Meeting (1986)
work page 1986
-
[6]
An experimental study of the flow field about an airfoil with glaze ice,
M. B. Bragg, S. A. Spring, "An experimental study of the flow field about an airfoil with glaze ice," AIAA Paper 87-0100, 25th Aerospace Sciences Meeting (1987)
work page 1987
-
[7]
Experimental measurements in a large separation bubble due to a simulated glaze ice accretion,
M. B. Bragg, A. Khodadoust, "Experimental measurements in a large separation bubble due to a simulated glaze ice accretion," AIAA Paper 88-0116, 26th Aerospace Sciences Meeting (1988)
work page 1988
-
[8]
Flowfield measurements about an airfoil with leading-edge ice shapes,
A. P. Broeren, M. B. Bragg, H. E. Addy Jr., "Flowfield measurements about an airfoil with leading-edge ice shapes," Journal of Aircraft 43(4) (2006) 1226– 1234
work page 2006
-
[9]
Numerical simulation of iced wing using separating shear layer fixed turbulence models,
H. Li, Y. Zhang, H. Chen, "Numerical simulation of iced wing using separating shear layer fixed turbulence models," AIAA Journal 59(6) (2021) 3667–3681. https://doi.org/10.2514/1.J060143. 45
-
[10]
Aerodynamic prediction and roughness implementation toward ice accretion simulation,
J. Chen, P. Yang, Y. Zhang, S. Fu, "Aerodynamic prediction and roughness implementation toward ice accretion simulation," Physics of Fluids 36(1) (2024)
work page 2024
-
[11]
J. Chen, Y. Zhang, S. Fu, "Ice accretion simulation incorporating roughness effects and separation correction through a transition model," Physics of Fluids 36(4) (2024)
work page 2024
-
[12]
Measurement of unsteady flow reattachment on an airfoil with an ice shape,
P. J. Ansell, M. B. Bragg, "Measurement of unsteady flow reattachment on an airfoil with an ice shape," AIAA Journal 52(3) (2014) 656–659. https://doi.org/10.2514/1.J052519
-
[13]
Numerical study of iced swept-wing performance degradation using RANS,
I. Ozcer, A. Pueyo, F. Menter, S. Hafid, H. Yang, "Numerical study of iced swept-wing performance degradation using RANS," in: International Conference on Icing of Aircraft, Engines, and Structures, Vienna, Austria, 2023, pp. 2023–01–1402. https://doi.org/10.4271/2023-01-1402
-
[14]
Improved prediction of flow around airfoil accreted with horn or ridge ice,
M. Xiao, Y. Zhang, "Improved prediction of flow around airfoil accreted with horn or ridge ice," AIAA Journal 59(6) (2021) 2318–2327
work page 2021
-
[15]
Numerical study on the aerodynamics of an iced airfoil with scale-resolving simulations,
M. L. Wong, A. S. Ghate, G. D. Stich, G. K. Kenway, C. C. Kiris, "Numerical study on the aerodynamics of an iced airfoil with scale-resolving simulations," AIAA Scitech 2023 Forum (2023)
work page 2023
-
[16]
Large-eddy simulations of complex aerodynamic flows over multi-element iced airfoils,
Y. M. Lee, J. H. Lee, L. P. Raj, J. H. Jo, R. S. Myong, "Large-eddy simulations of complex aerodynamic flows over multi-element iced airfoils," Aerospace Science and Technology 109 (2021) 106417
work page 2021
-
[17]
Aerodynamics of a swept wing with ice accretion at low Reynolds number,
J. Diebold, M. Monastero, M. Bragg, "Aerodynamics of a swept wing with ice accretion at low Reynolds number," AIAA Paper No. 2012-2795 (2012)
work page 2012
-
[18]
Effect of simulated scalloped ice on the aerodynamics of a swept-wing at low-Reynolds number,
N. Sandhu, M. R. Soltani, M. B. Bragg, C. W. Lum, B. Woodard, A. P. Broeren, et al., "Effect of simulated scalloped ice on the aerodynamics of a swept-wing at low-Reynolds number," in: 2018 Atmospheric and Space Environments Conference (NASA, 2018) p. 3495. 46
work page 2018
-
[19]
Effect of spanwise discontinuities on the aerodynamics of a swept wing with scalloped ice accretion,
B. S. Woodard, M. B. Bragg, "Effect of spanwise discontinuities on the aerodynamics of a swept wing with scalloped ice accretion," AIAA Paper No. 2020-2848 (2020)
work page 2020
-
[20]
Aerodynamic simulation of artificial scallop ice shapes on NACA 23012 airfoil,
R. Silva Reghin, T. Borges, R. C. de Sousa Sorbilli, T. Barbosa de Araújo, A. F. de Castro da Silva, "Aerodynamic simulation of artificial scallop ice shapes on NACA 23012 airfoil," AIAA SCITECH 2022 Forum (2022)
work page 2022
-
[21]
Numerical study of the effects of discontinuous ice on three-dimensional wings,
J. Chen, Y. Zhang, S. Fu, "Numerical study of the effects of discontinuous ice on three-dimensional wings," Physics of Fluids 36(7) (2024)
work page 2024
-
[22]
On the use of artificial ice shapes for large-eddy simulations in aircraft icing,
B. Bornhoft, P. Moin, S. S. Jain, S. T. Bose, "On the use of artificial ice shapes for large-eddy simulations in aircraft icing," Journal of Aircraft (2025) 1–18
work page 2025
-
[23]
Minimum-dissipation models for large-eddy simulation,
W. Rozema, H. J. Bae, P. Moin, R. Verstappen, "Minimum-dissipation models for large-eddy simulation," Physics of Fluids 27(1) (2015). https://doi.org/10.1063/1.4928700
-
[24]
Enhanced delayed detached-eddy simulation with anisotropic minimum dissipation subgrid length scale,
Z. Zhou, M. Xiao, D. Li, Y. Zhang, "Enhanced delayed detached-eddy simulation with anisotropic minimum dissipation subgrid length scale," Physics of Fluids 37(2) (2025)
work page 2025
-
[25]
A numerical investigation of sweep effects on turbulent flow over iced wings,
Z. Zhou, M. Xiao, J. Chen, L. Li, Y. Zhang, "A numerical investigation of sweep effects on turbulent flow over iced wings," arXiv preprint arXiv:2507.15182 (2025)
-
[26]
A hybrid RANS-LES approach with delayed-DES and wall-modelled LES capabilities,
M. L. Shur, P. R. Spalart, M. K. Strelets, A. K. Travin, "A hybrid RANS-LES approach with delayed-DES and wall-modelled LES capabilities," Int. J. Heat Fluid Flow 29(6) (2008) 1638–1649
work page 2008
-
[27]
Development of DDES and IDDES formulations for the K-x shear stress transport model,
M. S. Gritskevich, A. V. Garbaruk, J. Schütze, F. R. Menter, "Development of DDES and IDDES formulations for the K-x shear stress transport model," Flow. Turbul. Combust. 88(3) (2012) 431–449. 47
work page 2012
-
[28]
NASA Langley Research Center, "CFL3D version 6.7," https://nasa.github.io/CFL3D/ (2017)
work page 2017
-
[29]
M. Kloker, M. Borsboom, "A fully implicit scheme for unsteady flow calculations solved by the approximate factorization technique - implicit dual time stepping," Technical Report VKI Project Report 1986-16, VKI, Rhode St. Genese, Belgium (1986)
work page 1986
-
[30]
Enhanced prediction of three-dimensional finite iced wing separated flow near stall,
M. Xiao, Y. Zhang, F. Zhou, "Enhanced prediction of three-dimensional finite iced wing separated flow near stall," International Journal of Heat and Fluid Flow 98 (2022) 109067. https://doi.org/10.1016/j.ijheatfluidflow.2022.109067
-
[31]
DDES with adaptive coefficient for stalled flows past a wind turbine airfoil,
J. Liu, W. Zhu, Z. Xiao, H. Sun, Y. Huang, Z. Liu, "DDES with adaptive coefficient for stalled flows past a wind turbine airfoil," Energy 161 (2018). https://doi.org/10.1016/j.energy.2018.07.176
-
[32]
Aerodynamics of a finite wing with simulated ice,
A. Khodadoust, M. B. Bragg, "Aerodynamics of a finite wing with simulated ice," Journal of Aircraft 32(1) (1995) 137–144. https://doi.org/10.2514/3.46694
-
[33]
Effect of a simulated ice accretion on the aerodynamics of a swept wing,
M. Bragg, A. Khodadoust, R. Soltani, S. Wells, M. Kerho, "Effect of a simulated ice accretion on the aerodynamics of a swept wing," 29th Aerospace Sciences Meeting, AIAA Paper 1991-0442 (1991)
work page 1991
-
[34]
A. P. Broeren, S. Lee, C. M. LaMarre, M. B. Bragg, "Effect of airfoil geometry on performance with simulated ice accretions (volume 1: experimental investigation)," Report No. DOT/FAA/AR-03/64 (2003)
work page 2003
-
[35]
Velocity gradient partitioning in turbulent flows,
R. Arun, T. Colonius, "Velocity gradient partitioning in turbulent flows," Journal of Fluid Mechanics 1000 (2024) R5
work page 2024
-
[36]
R. Arun, M. Kamal, T. Colonius, P. L. Johnson "Normality-based analysis of multiscale velocity gradients and energy transfer in direct and large-eddy simulations of isotropic turbulence," arXiv preprint arXiv:2504.19356 (2025). 48
-
[37]
The triple decomposition of the velocity gradient tensor as a standardized real Schur form,
J. Kronborg, J. Hoffman, "The triple decomposition of the velocity gradient tensor as a standardized real Schur form," Phys. Fluids 35(3) (2023) 031703
work page 2023
-
[38]
A series in 1/√Re to represent the Strouhal– Reynolds number relationship of the cylinder wake,
C. H. K. Williamson, G. L. Brown, "A series in 1/√Re to represent the Strouhal– Reynolds number relationship of the cylinder wake," Journal of Fluids and Structures 1 (8) (1998) 1073–1085
work page 1998
-
[39]
Numerical simulation of aerodynamic features with ice shapes via high-fidelity CFD method,
H. Liu, C. Zhang, "Numerical simulation of aerodynamic features with ice shapes via high-fidelity CFD method," in: W. G. Habashi (Ed.), Handbook of Numerical Simulation of In-Flight Icing, Springer International Publishing, Cham (2023) 1–40. https://doi.org/10.1007/978-3-030-64725-4_45-1
-
[40]
C. Zhang, X. Tan, W. Xu, G. Li, W. Fuxin, H. Liu, "High-fidelity modeling of turbulent shear flow downstream of a 3-D airfoil with spanwise ice contamination leading stall," Computers & Fluids 240 (2022) 105423. https://doi.org/10.1016/j.compfluid.2022.105423
-
[41]
Numerical investigation of the unsteady flow past an iced multi-element airfoil,
M. Xiao, Y. Zhang, F. Zhou, "Numerical investigation of the unsteady flow past an iced multi-element airfoil," AIAA Journal 58(10) (2020) 3848–3862. https://doi.org/10.2514/1.J059114
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.