Search

A new efficient analytical-numerical method is developed for solving a problem for the potential vorticity equation in the quasi-geostrophic approximation with allowance for vertical diffusion of mass and momentum. The method is used to analyze small perturbations of ocean currents of finite transverse scale with a general parabolic vertical profile of velocity. For the arising spectral non-self-adjoint problem, asymptotic expansions of the eigenfunctions and eigenvalues are constructed for small wave numbers and the existence of a countable set of complex eigenvalues with an unboundedly decreasing imaginary part is shown. On the integration interval , a system of three neighborhoods is introduced and a solution in each of them is constructed in the form of power series expansions, which are matched smoothly, so that the eigenfunctions and eigenvalues are efficiently calculated with high accuracy. For a varying wave number , the trajectories of complex eigenvalues are computed for various parameters of the problem and the existence of double eigenvalues is shown. The complex picture of instability developing in the simulated flow depending on physical parameters of the problem is briefly described. [ABSTRACT FROM AUTHOR]

For the problem of conformal mapping of a half-plane onto a -shaped domain with arbitrary geometry, an efficient method is developed for finding parameters of the Schwarz–Christoffel integral, i.e., the preimages of the vertices (prevertices) and the pre-integral multiplier. Special attention is given to the situation of crowding prevertices, in which case conventional integration methods face significant difficulties. For this purpose, the concept of a cluster is introduced, its center is determined, and all integrand binomials with prevertices from this cluster are expanded into a fast-convergent series by applying a unified scheme. Next, the arising integrals are reduced to single or double series in terms of Gauss hypergeometric functions . The fast convergence of the resulting expansions is ensured by applying formulas for analytic continuation of to a neighborhood of the point q = 1 and using numerically stable recurrence relations. The constructed expansions are also fairly efficient for choosing initial approximations for prevertices in Newton's iteration method. By using the leading terms of these expansions, the approximations for the prevertices are expressed in explicit form in terms of elementary functions, and the subsequent iterations ensure the fast convergence of the algorithm. After finding the parameters in the integral, the desired mapping is constructed as a combination of power series expansions at prevertices, regular expansions at the preimage of the center of symmetry, a Laurent series in a semi-annulus, and special series near the preimages of the vertical segments. Numerical results demonstrate the high efficiency of the developed method, especially in the case of strong crowding of prevertices. [ABSTRACT FROM AUTHOR]