On the use of asymptotically motivated gauge functions to obtain convergent series solutions to nonlinear ODEs

• Published in: IMA journal of applied mathematics Vol. 88; no. 1; pp. 43 - 66
• Main Authors: Naghshineh, Nastaran, Reinberger, W Cade, Barlow, Nathaniel S, Samaha, Mohamed A, Weinstein, Steven J
• Format: Journal Article
• Language: English
• Published: 01.04.2023
• ISSN: 0272-4960; 1464-3634
• DOI: 10.1093/imamat/hxad006

We examine the power series solutions of two classical nonlinear ordinary differential equations of fluid mechanics that are mathematically related by their large-distance asymptotic behaviours in semi-infinite domains. The first problem is that of the ‘Sakiadis’ boundary layer over a moving flat wall, for which no exact analytic solution has been put forward. The second problem is that of a static air–liquid meniscus with surface tension that intersects a flat wall at a given contact angle and limits to a flat pool away from the wall. For the latter problem, the exact analytic solution—given as distance from the wall as a function of meniscus height—has long been known (Batchelor, G. K. (1967) An Introduction to Fluid Dynamics, chapter 1: The physical properties of fluids. Cambridge). Here, we provide an explicit solution as meniscus height versus distance from the wall to elucidate structural similarities to the Sakiadis boundary layer. Although power series solutions are readily obtainable to the governing nonlinear ordinary differential equations, we show that—in both problems—the series diverge due to non-physical complex or negative real-valued singularities. In both cases, these singularities can be moved by expanding in exponential gauge functions motivated by their respective large distance asymptotic behaviours to enable series convergence over their full semi-infinite domains. For the Sakiadis problem, this not only provides a convergent Taylor series (and conjectured exact) solution to the ODE, but also a means to evaluate the wall shear parameter (and other properties) to within any desired precision. Although the nature of nonlinear ODEs precludes general conclusions, our results indicate that asymptotic behaviours can be useful when proposing variable transformations to overcome power series divergence.