Steady advection–diffusion around finite absorbers in two-dimensional potential flows

We consider perhaps the simplest non-trivial problem in advection–diffusion – a finite absorber of arbitrary cross-section in a steady two-dimensional potential flow of concentrated fluid. This problem has been studied extensively in the theory of solidification from a flowing melt, and it also arises in advection–diffusion-limited aggregation. In both cases, the fundamental object is the flux to a circular disk, obtained by conformal mapping from more complicated shapes. Here, we construct an accurate numerical solution by an efficient method that involves mapping to the interior of the disk and using a spectral method in polar coordinates. The method combines exact asymptotics and an adaptive mesh to handle boundary layers. Starting from a well-known integral equation in streamline coordinates, we also derive high-order asymptotic expansions for high and low Péclet numbers ($Pe$). Remarkably, the ‘high’$Pe$ expansion remains accurate even for such low $Pe$ as $10^{-3}$. The two expansions overlap well near $Pe\,{=}\,0.1$, allowing the construction of an analytical connection formula that is uniformly accurate for all $Pe$ and angles on the disk with a maximum relative error of 1.75%. We also obtain an analytical formula for the Nusselt number ($Nu$) as a function of $Pe$ with a maximum relative error of 0.53% for all possible geometries after conformal mapping. Considering the concentration disturbance around a disk, we find that the crossover from a diffusive cloud (at low $Pe$) to an advective wake (at high $Pe$) occurs at $Pe\,{\approx}\,60$.

[1]  T. Carleman Über die Abelsche Integralgleichung mit konstanten Integrationsgrenzen , 1922 .

[2]  N. Goldenfeld Lectures On Phase Transitions And The Renormalization Group , 1972 .

[3]  D. F. Hays,et al.  Table of Integrals, Series, and Products , 1966 .

[4]  L. Trefethen Numerical conformal mapping , 1986 .

[5]  E. C. Titchmarsh Introduction to the Theory of Fourier Integrals , 1938 .

[6]  V. M. Aleksandrov,et al.  Asymptotic solution of a class of integral equations and its application to contact problems for cylindrical elastic bodies , 1967 .

[7]  Benny Davidovitch,et al.  Dynamics of conformal maps for a class of non-Laplacian growth phenomena. , 2003, Physical review letters.

[8]  Andreas Acrivos,et al.  Heat and Mass Transfer from Single Spheres in Stokes Flow , 1962 .

[9]  A. Erdélyi,et al.  Higher Transcendental Functions , 1954 .

[10]  T. MacRobert Higher Transcendental Functions , 1955, Nature.

[11]  J. Ahner On the solution of a class of integral equations*1 , 1984 .

[12]  Christopher G. Phillips,et al.  HEAT AND MASS TRANSFER FROM A FILM INTO STEADY SHEAR FLOW , 1990 .

[13]  G. Batchelor,et al.  An Introduction to Fluid Dynamics , 1968 .

[14]  J. Conway,et al.  Functions of a Complex Variable , 1964 .

[15]  V. M. Aleksandrov,et al.  Asymptotic solution of a class of integral equations encountered in the investigation of mixed problems of the mathematical physics for regions with cylindrical boundaries , 1968 .

[16]  L. Trefethen Spectral Methods in MATLAB , 2000 .

[17]  W. Hackbusch Singular Integral Equations , 1995 .

[18]  Vimal Singh,et al.  Perturbation methods , 1991 .

[19]  J. Burgers A mathematical model illustrating the theory of turbulence , 1948 .

[20]  I. S. Gradshteyn,et al.  Table of Integrals, Series, and Products , 1976 .

[21]  Sam Howison,et al.  Complex variable methods in Hele–Shaw moving boundary problems , 1992, European Journal of Applied Mathematics.

[22]  Oono,et al.  Renormalization group and singular perturbations: Multiple scales, boundary layers, and reductive perturbation theory. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[23]  J. Myers Wave Scattering and the Geometry of a Strip , 1965 .

[24]  T. Pedley,et al.  Nutrient Uptake by a Self‐Propelled Steady Squirmer , 2003 .

[25]  Redner,et al.  Tracer dispersion in planar multipole flows. , 1994, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[26]  The equilibrium shape of an ice-soil body formed by liquid flow past a pair of freezing columns☆ , 1994 .

[27]  I. Stakgold Green's Functions and Boundary Value Problems , 1979 .

[28]  S. Howison,et al.  Two-dimensional solidification and melting in potential flows , 1999, Journal of Fluid Mechanics.

[29]  V. L. Rvachev,et al.  Plate in the form of an infinite strip on an elastic half-space: PMM vol. 40, n≗ 2, 1976, pp. 298–305 , 1976 .

[30]  J. H. PEARCE,et al.  Functions of a Complex Variable , 1947, Nature.

[31]  C. Pearson On the finite strip problem , 1957 .

[32]  Martin Z. Bazant Conformal mapping of some non-harmonic functions in transport theory , 2004, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[33]  L. Gary Leal,et al.  Laminar Flow and Convective Transport Processes , 1992 .

[34]  Determination of the equilibrium shape of the bodies formed during the solidification of filtration flow , 1988 .

[35]  Ian Eames,et al.  Longitudinal dispersion by bodies fixed in a potential flow , 1999, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[36]  George Weiss,et al.  Methods Based on the Wiener-Hopf Technique for the Solution of Partial Differential Equations , 1958 .

[37]  G. Mukhamadullina,et al.  Hysteretic Effects in the Problems of Artificial Freezing , 1998, SIAM J. Appl. Math..

[38]  H. Brenner Forced convection heat and mass transfer at small Péclet numbers from a particle of arbitrary shape , 1963 .

[39]  K. Kornev,et al.  Mathematical theory of freezing for flow in porous media , 1994, Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences.

[40]  Adrian Bejan,et al.  Heatline visualization of forced convection in porous media , 1994 .

[41]  J. Myers Derivation of a matrix Painlevé equation germane to wave scattering by a broken corner , 1984 .

[42]  L. Cummings,et al.  Evolution of flat crystallisation front in forced hydrodynamic flow—some explicit solutions , 1999 .

[43]  M. G. Krein,et al.  Integral equations on a half-line with kernels de-pending upon the difference of the arguments , 1962 .