Determining the limiting solutions of nonstationary axisymmetric Hele–Shaw problems

Numerical solution of the Hele–Shaw problem reduces to solution of three boundary-value problems of determining analytic functions of a complex variable in each time step: conformal mapping of the range of the parametric variable to the physical plane, the Dirichlet problems for determining the electric-field strength, and the Riemann–Hilbert problem for calculating partial time derivatives of the coordinates of points of the interelectrode space (the images of the points on the boundary of the parametric plane are fixed). Unlike in the two-dimensional problem, the electric-field strength is determined using integral transformations of an analytic function. Approximation by spline function is performed, and more accurate and steady (than the well-known ones) general solution algorithms for the nonstationary axisymmetric problems are described. Results of a numerical study of the formation of stationary and self-similar configurations are presented.