pith. sign in

arxiv: 1907.09912 · v1 · pith:6OMHUA5Enew · submitted 2019-07-17 · ⚛️ physics.comp-ph

Asymptotic Approximant for the Falkner-Skan Boundary-Layer equation

Elizabeth R. Belden , Zachary A. Dickman , Steven J. Weinstein , Alex D. Archibee , Ethan Burroughs , Nathaniel S. Barlow This is my paper

Pith reviewed 2026-05-24 19:47 UTC · model grok-4.3

classification ⚛️ physics.comp-ph
keywords Falkner-Skan equationasymptotic approximantboundary layer flowwedge flowBlasius equationanalytic solutionsradius of convergencepressure gradient parameter
0
0 comments X

The pith

The asymptotic approximant developed for Blasius flow extends to deliver accurate closed-form analytic solutions for the Falkner-Skan equation across wedge angles with beta from -0.1988 to 1.

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

The paper shows that the same recursive asymptotic approximant previously constructed for the Blasius flat-plate boundary layer produces accurate analytic closed-form solutions when applied to the Falkner-Skan equation that governs flow over a wedge. The resulting expressions cover the full practical range of the pressure-gradient parameter beta and correctly capture both the unique solutions for favorable gradients and the non-unique positive- and negative-shear solutions that appear for adverse gradients. A reader would care because the construction is performed by simple recursion, avoids the matrix inversion required by Padé approximants, directly locates the complex-plane singularities that limit the radius of convergence, and supplies benchmark far-field constants.

Core claim

The asymptotic approximant applied to the Blasius boundary layer flow over a flat plate yields accurate analytic closed-form solutions to the Falkner-Skan boundary layer equation for flow over a wedge having angle βπ/2 to the horizontal, for β∈[-0.198837735,1], accurately representing the previously established non-unique solutions for β<0 having positive and negative shear rates along the wedge.

What carries the argument

The asymptotic approximant: a recursively constructed analytic function that matches the near-wedge power series to the far-field asymptotic behavior of the Falkner-Skan solution.

If this is right

  • Accurate closed-form representations exist for both unique and non-unique solutions across the full practical range of wedge angles.
  • Complex-plane singularities that set the radius of convergence of the power series can be read directly from the approximant.
  • The approximant is constructed by recursion without matrix inversion, in contrast to traditional Padé methods.
  • Benchmark constants characterizing the far-field asymptotic decay are obtained for each beta.

Where Pith is reading between the lines

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

  • The same recursive matching procedure may transfer without modification to other nonlinear boundary-layer ODEs that possess similar far-field decay.
  • Rapid analytic parameter studies over beta become feasible once the approximant coefficients are tabulated.
  • The extracted complex singularities could guide the construction of improved series resummations for related fluid-mechanics problems.

Load-bearing premise

The recursive construction rules and matching conditions developed for the Blasius case remain valid when the pressure-gradient parameter beta is introduced and when multiple solutions exist for beta less than zero.

What would settle it

High-precision numerical integration of the Falkner-Skan equation for any beta inside the claimed interval that deviates from the approximant by more than the expected truncation error would falsify the accuracy claim.

Figures

Figures reproduced from arXiv: 1907.09912 by Alex D. Archibee, Elizabeth R. Belden, Ethan Burroughs, Nathaniel S. Barlow, Steven J. Weinstein, Zachary A. Dickman.

Figure 1
Figure 1. Figure 1: Geometry of flow over a wedge showing potential flow and viscous boundary layer regions. a) Domain for the potential flow solution for b>0, from which the nonzero velocity U(x) and associated Bernouli pressure P(x) along y=0 is extracted. b) Domain for the thin viscous boundary layer (Falkner-Skan problem) that is driven by the potential flow solutions for U(x) and P(x). There are no pressure variations in… view at source ↗
Figure 2
Figure 2. Figure 2: Solution properties as a function of the wedge angle parameter, [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: (a) The N-term series (4) labeled SN and approximant (6) labeled AN compared with numerical solution (•). (b) Derivatives of approximant (6) for N=20. Data shown here corresponds to conditions at point I in [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: (a) The [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: (a) The N-term series (4) labeled SN and approximant (6) labeled AN compared with numerical solution (•). (b) Derivatives of approximant (6) for N=20. Data shown here corresponds to conditions at point III in [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: (a) The N-term series (4) labeled SN and approximant (6) labeled AN compared with numerical solution (•). (b) Derivatives of approximant (6) for N=20. Data shown here corresponds to conditions at point IV in [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: (a) The [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Error associated with approximant (6), measured against numerical solutions with [PITH_FULL_IMAGE:figures/full_fig_p011_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: (a) The N-term series (4) labeled SN and approximant (6) labeled AN compared with numerical solution (•). (b) Derivatives of approximant (6) for N=30. Data shown here corresponds to conditions at point VI in [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: (a) The N-term series (4) labeled SN and approximant (7) labeled AN compared with numerical solution (•). (b) Derivatives of approximant (7) for N=25. Data shown here corresponds to conditions at point VI in [PITH_FULL_IMAGE:figures/full_fig_p012_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: (a) Location of singularities that limit convergence of the Falkner-Skan series (4), shown [PITH_FULL_IMAGE:figures/full_fig_p012_11.png] view at source ↗
read the original abstract

We demonstrate that the asymptotic approximant applied to the Blasius boundary layer flow over a flat plat (Barlow et al., 2017 Q. J. Mech. Appl. Math., 70(1): 21-48) yields accurate analytic closed-form solutions to the Falkner-Skan boundary layer equation for flow over a wedge having angle $\beta\pi/2$ to the horizontal. A wide range of wedge angles satisfying $\beta\in[-0.198837735, 1]$ are considered, and the previously established non-unique solutions for $\beta<0$ having positive and negative shear rates along the wedge are accurately represented. The approximant is used to determine the singularities in the complex plane that prescribe the radius of convergence of the power series solution to the Falkner-Skan equation. An attractive feature of the approximant is that it may be constructed quickly by recursion compared with traditional Pad\'e approximants that require a matrix inversion. The accuracy of the approximant is verified by numerical solutions, and benchmark numerical values are obtained that characterize the asymptotic behavior of the Falkner-Skan solution at large distances from the wedge.

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 manuscript claims that the asymptotic approximant previously developed for the Blasius boundary-layer equation (β=0) directly produces accurate closed-form analytic solutions to the Falkner-Skan equation f''' + f f'' + β(1 - (f')²) = 0 for wedge flows over the range β ∈ [-0.198837735, 1]. It asserts that the same recursive construction accurately represents the non-unique solutions for β < 0, locates singularities in the complex plane that determine the power-series radius of convergence, and supplies benchmark values for the far-field asymptotic behavior, all verified against independent numerical solutions.

Significance. If the extension of the approximant holds, the work supplies a computationally inexpensive analytic representation for a classic family of boundary-layer problems, with the practical advantage of recursive construction that avoids the matrix inversion required by Padé approximants. The provision of benchmark asymptotic constants and complex-plane singularity locations would also be of direct use to the computational fluid-dynamics community.

major comments (3)
  1. [§3] §3 (construction of the approximant): the recursion and matching conditions are presented as a direct carry-over from the Blasius case (Barlow et al., 2017). Because the β term modifies both the recurrence relation for the power-series coefficients a_n and the far-field decay rate of f' → 1, the manuscript must explicitly derive or adjust the recursion for nonzero β; without this step the claimed applicability to the full Falkner-Skan family is not established.
  2. [§4] §4 (multiple solutions for β < 0): the paper states that both the positive- and negative-shear branches are accurately represented, yet no separate error tables or convergence plots are supplied for the negative-shear branch at, e.g., β = −0.1. The central claim that the same approximant captures both branches therefore lacks the quantitative support required to be load-bearing.
  3. [§5] §5 (complex singularities): the radius-of-convergence results are obtained by locating poles of the approximant, but the manuscript does not demonstrate that these poles remain accurate when β is varied away from zero; a direct comparison of the predicted radius versus the numerically determined radius of convergence for at least two additional β values is needed.
minor comments (2)
  1. [Abstract] The abstract cites the range β ∈ [-0.198837735, 1] but does not state the precise numerical criterion used to truncate the interval at the lower end; a brief sentence in §2 would clarify this choice.
  2. [Figures] Figure captions should explicitly label which curve corresponds to the positive-shear and which to the negative-shear branch when β < 0.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We address each major comment below.

read point-by-point responses

  1. Referee: [§3] §3 (construction of the approximant): the recursion and matching conditions are presented as a direct carry-over from the Blasius case (Barlow et al., 2017). Because the β term modifies both the recurrence relation for the power-series coefficients a_n and the far-field decay rate of f' → 1, the manuscript must explicitly derive or adjust the recursion for nonzero β; without this step the claimed applicability to the full Falkner-Skan family is not established.

    Authors: We agree that an explicit derivation of the recursion for general β is required. In the revised manuscript we will substitute the power series into the Falkner-Skan equation, collect terms, and display the resulting recurrence relation that incorporates the β-dependent contributions to the coefficients a_n together with the adjusted far-field matching conditions. revision: yes

  2. Referee: [§4] §4 (multiple solutions for β < 0): the paper states that both the positive- and negative-shear branches are accurately represented, yet no separate error tables or convergence plots are supplied for the negative-shear branch at, e.g., β = −0.1. The central claim that the same approximant captures both branches therefore lacks the quantitative support required to be load-bearing.

    Authors: We acknowledge the need for dedicated quantitative support. The revised manuscript will include separate error tables and convergence plots for the negative-shear branch, with explicit results shown at β = −0.1 and at least one additional negative value. revision: yes

  3. Referee: [§5] §5 (complex singularities): the radius-of-convergence results are obtained by locating poles of the approximant, but the manuscript does not demonstrate that these poles remain accurate when β is varied away from zero; a direct comparison of the predicted radius versus the numerically determined radius of convergence for at least two additional β values is needed.

    Authors: We will add a direct comparison of approximant-predicted radii against numerically computed radii of convergence for at least two additional β values (β = 0.5 and β = −0.1) to confirm accuracy of the singularity locations away from the Blasius case. revision: yes

Circularity Check

0 steps flagged

Minor self-citation of approximant method; central extension to Falkner-Skan is independently verified

full rationale

The paper cites its own prior work (Barlow et al. 2017) for the asymptotic approximant construction developed on the Blasius equation (β=0). However, the present manuscript applies the method to the Falkner-Skan equation for a range of β values, determines singularities, and explicitly verifies accuracy against independent numerical solutions. The load-bearing claim does not reduce solely to the self-citation; the numerical benchmarks and explicit demonstration for nonzero β supply external grounding. This qualifies as normal self-citation of a method with independent validation rather than circularity.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the transferability of the recursive approximant construction from the Blasius case; no free parameters or invented entities are mentioned in the abstract.

axioms (1)
  • domain assumption The asymptotic approximant form and recursion rules developed for the Blasius equation remain valid when the Falkner-Skan pressure-gradient parameter β is introduced.
    The paper states that the approximant 'yields accurate analytic closed-form solutions' for the Falkner-Skan equation, implying this transferability.

pith-pipeline@v0.9.0 · 5754 in / 1300 out tokens · 20217 ms · 2026-05-24T19:47:05.687409+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

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

  1. [1]

    Introduction The Falkner-Skan equation describes boundary layer flow over a wedge of angle βπ/2 to the horizontal that is driven by an external pressure gradient predicted from potential flow (see Fig. 1). The equation also applies to regimes where the pressure gradient opposes the flow when β < 0 (Fig. 1c,d) for which boundary layer separation may occur. Th...

    work page internal anchor Pith review Pith/arXiv arXiv 1907

  2. [2]

    reversed flow

    Asymptotic properties and series expansion of the F alkner-Skan equation Solutions to the Falkner-Skan equation system (1) are found for a given flow by fixing the parameter β, which is related to the wedge angle shown in Fig. 1. Note that this parameter incorporates the effect of wedge angle on the potential flow that drives the fluid motion in the boundary l...

    work page 1988

  3. [3]

    Asymptotic Approximant The divergence of the Falkner-Skan series (4) demonstrated in Figs. 3a through 9a is overcome using the method of asymptotic approximants, which constrains the analytic continuation of the series via an asymptotic behavior away from the point of expansion [15]. An approximant that satisfies the η→∞ behavior (5) is given as fA = η + B...

    work page

  4. [4]

    This adds to the increasing number of problems in disparate areas of mathematical physics to which asymptotic approximants have been applied successfully [12, 13, 14, 15, 18, 17]

    Conclusions Asymptotic approximants provide nearly exact closed-form solutions to the Falkner-Skan boundary layer equation for varying wedge angle. This adds to the increasing number of problems in disparate areas of mathematical physics to which asymptotic approximants have been applied successfully [12, 13, 14, 15, 18, 17]. Advantages of asymptotic appr...

    work page

  5. [5]

    Schlichting

    H. Schlichting. Boundary Layer Theory . McGraw-Hill, seventh edition, 1979

    work page 1979

  6. [6]

    Cebeci and H

    T. Cebeci and H. B. Keller. Shooting and parallel shooting methods for solving the Falkner-Skan boundary-layer equation. J. Comp. Phys. , 7:289–300, 1971

    work page 1971

  7. [7]

    Laine and L

    C. Laine and L. Reinhart. Further numerical methods for the Falkner-Skan equaitons: shooting and continuation techniques. I. J. Num. Meth. Fluids , 4:833–852, 1984

    work page 1984

  8. [8]

    R. Fazio. A novel approach to the numerical solution of boundary value problems on infinite intervals. SIAM J. Num. Anal. , 33(4):1473–1483, 1996

    work page 1996

  9. [9]

    Asaithambi

    A. Asaithambi. A finite-difference method for the Falkner-Skan equation. Appl. Math. Comp. , 92:135–141, 1998

    work page 1998

  10. [10]

    S. S. Motsa and P. Sibanda. An efficient numerical method for solving Falkner-Skan boundary layer flows. I. J. Num. Meth. Fluids , 69:499–508, 2012

    work page 2012

  11. [11]

    R. Fazio. Blasius problem and Falkner-Skan model: T¨ opfer’s algorithm and its extension. Comp. Fluids , 73(15):202–209, 2013

    work page 2013

  12. [12]

    C-S. Liu. An iterative method based-on eigenfunctions and adjoint eigenfunctions for solving the Falkner-Skan equation. Appl. Math. Lett. , 67:33–39, 2017

    work page 2017

  13. [13]

    Bararnia, N

    H. Bararnia, N. Haghparast, M. Miansari, and A. Barari. Flow analysis for the Falkner-Skan wedge flow. Curr. Sci., 103(2):169–177, 2012

    work page 2012

  14. [14]

    B. I. Yun. New approximate analytical solutions of the Falkner-Skan equation. J. Appl. Math. , 2012:170802: 1–12, 2012

    work page 2012

  15. [15]

    A. Khidir. A note on the solution of general Falkner-Skan problem by two novel semi-analytical techniques. Propul. Power Res. , 4(4):212–220, 2015

    work page 2015

  16. [16]

    N. S. Barlow, A. J. Schultz, S. J. Weinstein, and D. A. Kofke. An asymptotically consistent approximant method with application to soft- and hard-sphere fluids. J. Chem. Phys. , 137:204102, 2012

    work page 2012

  17. [17]

    N. S. Barlow, A. J. Schultz, S. J. Weinstein, and D. A. Kofke. Critical isotherms from virial series using asymptotically consistent approximants. AIChE J. , 60(9):3336–3349, 2014

    work page 2014

  18. [18]

    N. S. Barlow, A. J. Schultz, S. J. Weinstein, and D. A. Kofke. Communication: Analytic continutation of the virial series through the critical point using parametric approximants. J. Chem. Phys. , 143:071103:1–5, 2015

    work page 2015

  19. [19]

    N. S. Barlow, C. R. Stanton, N Hill, S. J. Weinstein, and A. G. Cio. On the summation of divergent, truncated, and underspecified power series via asymptotic approximants. Q. J. Mech. Appl. Math. , 70(1):21–48, 2017

    work page 2017

  20. [20]

    C. M. Bender and S. A. Orszag. Advanced Mathematical Methods for Scientists and Engineers I: Asymptotic Methods and Perturbation Theory . McGraw-Hill, 1978

    work page 1978

  21. [21]

    R. J. Beachley, M. Mistysyn, J. A. Faber, S. J. Weinstein, and N. S. Barlow. Accurate closed-form trajectories of light around a kerr black hole using asymptotic approximants. Class. Quant. Grav. , 35(20):1–28, 2018

    work page 2018

  22. [22]

    N. S. Barlow, S. J. Weinstein, and J. A. Faber. An asymptotically consistent approximant for the equatorial bending angle of light due to kerr black holes. Class. Quant. Grav. , 34(135017):1–16, 2017

    work page 2017

  23. [23]

    J. P. Boyd. The Blasius function: Computations before computers, the value of tricks, undergraduate projects, and open research problems. SIAM Rev., 50(4):791–804, 2008

    work page 2008

  24. [24]

    Anil Lal and Neeraj Paul M

    S. Anil Lal and Neeraj Paul M. An accurate Taylor series solution with high radius of convergence for the Blasius function and parameters of asymptotic variation. JAFM, 7(4):557–564, 2014

    work page 2014

  25. [25]

    Stewartson

    K. Stewartson. Further solutions of the Falkner-Skan equation. Proc. Roy. Phil. Soc. , 50:454–465, 1954

    work page 1954

  26. [26]

    See (insert url provided by journal) for a spreadsheet of values that were computed in this work

    work page

  27. [27]

    P. L. Sachdev, R. B. Kudenatti, and N. M. Bujurke. Exact analytic solution of a boundary value problem for the Falkner-Skan equation. Stud. Appl. Math , 120:1–16, 2008

    work page 2008

  28. [28]

    W. A. Coppel. On a differential equation of boundary layer theory. Phil. Trans. R. Soc. Lond. A , 253(1023):101– 136, 1960

    work page 1960

  29. [29]

    R. V. Churchill. Complex Variables. McGraw-Hill, 1948

    work page 1948

  30. [30]

    Gonnet, S

    P. Gonnet, S. G¨ uttel, and L. N. Trefethen. Robust Pad´ e approximation via svd. SIAM Rev. , 55(1):101–117, asymptotic approximants 14 2013

    work page 2013

  31. [31]

    G. A. Baker Jr. and J. L. Gammel. The Pad´ e approximant. J. Math. Anal. Appl. , 2:21–30, 1961

    work page 1961

  32. [32]

    E. C. Titchmarsh. The Theory of Functions . Oxford University Press, 2 edition, 1968

    work page 1968

  33. [33]

    Van Dyke

    M. Van Dyke. Perturbation Methods in Fluid Mechanics . Parabolic, 1975

    work page 1975

  34. [34]

    A. J. Guttmann and I. Jensen. Series analysis. In Polygons, Polyominoes, and Polycubes , pages 181–202. Springer, 2009

    work page 2009

  35. [35]

    J. P. Boyd. The Blasius function in the complex plane. Exper. Math., 8(4):381–394, 1999

    work page 1999