On the linear stability study of zonal incompressible flows on a sphere

The normal mode (linear) stability of zonal flows of a nondivergent fluid on a rotating sphere is considered. The spherical harmonics are used as the basic functions on the sphere. The stability matrix representing in this basis the vorticity equation operator linearized about a zonal flow is analyzed in detail using the recurrent formula derived for the nonlinear triad interaction coefficients. It is shown that the zonal flow having the form of a Legendre polynomial P_n(μ) of degree n is stable to infinitesimal perturbations of every invariant set I_m with |m| ≥ n. For each zonal number m, I_m is here the span of all the spherical harmonics $Y^{m}_{k}(x)$, whose degree k is greater than or equal to m. It is also shown that such small-scale perturbations are stable not only exponentially, but also algebraically. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 649–665, 1998     

Hydromagnetic instability of a power-law liquid film flowing down a vertical cylinder using numerical approximation approach techniques

The long-wave perturbation method is employed to investigate the hydromagnetic stability of a thin electrically-conductive power-law liquid film flowing down the external surface of a vertical cylinder in a magnetic field. The validity of the numerical results is improved through the introduction of the flow index and the magnetic force into the governing equation. In contrast to most previous studies presented in the literature, the solution scheme employed in this study is based on a numerical approximation approach rather than an analytical method. The normal mode approach is used to analyze the stability of the film flow. The modeling results reveal that the stability of the film flow system is weakened as the radius of the cylinder is reduced. However, the flow stability can be enhanced by increasing the intensity of the magnetic field and the flow index, respectively. In general, the optimum conditions can be found through the use of a system to alter stability of the film flow by controlling the applied magnetic field.     

On the spectral problem in the linear stability study of flows on a sphere

A normal mode instability study of a steady nondivergent flow on a rotating sphere is considered. A real-order derivative and family of the Hilbert spaces of smooth functions on the unit sphere are introduced, and some embedding theorems are given. It is shown that in a viscous fluid on a sphere, the operator linearized about a steady flow has a compact resolvent, that is, a discrete spectrum with the only possible accumulation point at infinity, and hence, the dimension of the unstable manifold of a steady flow is finite. Peculiarities of the operator spectrum in the case of an ideal flow on a rotating sphere are also considered. Finally, as examples, we consider the normal mode stability of polynomial (zonal) basic flows and discuss the role of the linear drag, turbulent diffusion and sphere rotation in the normal mode stability study.     

A boundary integral method for the numerical computation of the forces exerted on a sphere in viscous incompressible flows near a plane wall

Slow uniform flows of a viscous, incompressible fluid past a rigid sphere near a plane wall are considered. The drag and lateral forces exerted on the sphere by the fluid are computed. The numerical results are compared with existing theoretical and experimental data for these and related fluid flows. They are based on a boundary element method for various linearized boundary value problems.

---

Zusammenfassung Es werden langsame, gleichförmige Strömungen einer zähen, inkompressiblen Flüssigkeit um eine starre Kugel in der Nähe einer ebenen Wand betrachtet. Die von der Flüssigkeit auf die Kugel ausgeübten Widerstands- und Querkräfte werden berechnet. Die numerischen Ergebnisse werden mit bestehenden, theoretischen und experimentellen Daten für diese und verwandte Strömungen verglichen. Sie basieren auf einer Randelementmethode für verschiedene, linearisierte Randwertprobleme.

     

Stable Lagrangian numerical differentiation with the highest order of approximation

Some asymptotic representations for the truncation error for the Lagrangian numerical differentiation are presented, when the ratio of the distance between each interpolation node and the differentiated point to step-parameter h is known. Furthermore, if the sampled values of the function at these interpolation nodes have perturbations which are bounded by ε, a method for determining step-parameter h by means of perturbation bound ε and order n of interpolation is provided to saturate the order of approximation. And all the investigations in this paper can be generalized to the set of quasi-uniform nodes.     

Polynomial approximation on the sphere using scattered data

We consider the problem of approximately reconstructing a function f defined on the surface of the unit sphere in the Euclidean space ℝ^{q +1} by using samples of f at scattered sites. A central role is played by the construction of a new operator for polynomial approximation, which is a uniformly bounded quasi‐projection in the de la Vallée Poussin style, i.e. it reproduces spherical polynomials up to a certain degree and has uniformly bounded L^p operator norm for 1 ≤ p ≤ ∞. Using certain positive quadrature rules for scattered sites due to Mhaskar, Narcowich and Ward, we discretize this operator obtaining a polynomial approximation of the target function which can be computed from scattered data and provides the same approximation degree of the best polynomial approximation. To establish the error estimates we use Marcinkiewicz–Zygmund inequalities, which we derive from our continuous approximating operator. We give concrete bounds for all constants in the Marcinkiewicz–Zygmund inequalities as well as in the error estimates. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Spectral approximation for drainage of an Elastico‐viscous liquid and error analysis

A Fourier‐Galerkin spectral method is proposed and used to analyze a system of quasilinear partial differential equations governing the drainage of liquids of the Oldroyd four‐constant type. It is shown that, Fourier‐Galerkin approximations are convergent with spectral accuracy. An efficient and accurate algorithm based on the Fourier‐Galerkin approximations to the system of quasilinear partial differential equations are developed and implemented. Numerical results indicating the high accuracy and effectiveness of this algorithm are presented. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 492–505, 2012     

The capability of approximation for neural networks interpolant on the sphere

Compared with planar hyperplane, fitting data on the sphere has been an important and active issue in geoscience, metrology, brain imaging, and so on. In this paper, using a functional approach, we rigorously prove that for given distinct samples on the unit sphere there exists a feed‐forward neural network with single hidden layer which can interpolate the samples, and simultaneously near best approximate the target function in continuous function space. Also, by using the relation between spherical positive definite radial basis functions and the basis function on the Euclidean space ?^{d + 1}, a similar result in a spherical Sobolev space is established. Copyright © 2010 John Wiley & Sons, Ltd.     

Spectral Galerkin approximation and rigorous error analysis for the Steklov eigenvalue problem in circular domain

In this paper, we present spectral Galerkin approximation and rigorous error analysis for the Steklov eigenvalue problem in a circular domain. First of all, we use the polar coordinate transformation and technique of separation of variables to reduce the problem to a sequence of equivalent 1‐dimensional eigenvalue problems that can be solved individually in parallel. Then, we derive the pole conditions and introduce weighted Sobolev space according to pole conditions. Together with the approximate properties of orthogonal polynomials, we prove the error estimates of approximate eigenvalues for each 1‐dimensional eigenvalue problem. Finally, we provide some numerical experiments to validate the theoretical results and algorithms.     

Investigation of the stability of periodic flows of a viscous fluid

An exact expression is obtained for the critical Reynolds number (R_*) for loss of stability in a wide class of one-dimensional periodic flows. An evolutionary equation is derived in the case of a small subcritically (R − R_* 1) which describes the dynamics of the secondary vortex structure.     

Legendre rational approximation on the whole line

The Legendre rational approximation is investigated. Some approximation results are established, which form the mathematical foundation of a new spectral method on the whole line. A model problem is considered. Numerical results show the efficiency of this new approach.     

Spectral approximation of Wiener-Hopf operators with almost periodic generating function

The paper deals with spectral approximation of Wiener-Hopf operators acting on L^p -spaces by their

finite sections. The generating functions of the Wiener-Hopf operators are supposed to be continuous plus almost

periodic.While the usual spectra of the finite sections drastically fail to converge to the spectrum of the Wiener-Hopf

operator,it turns out that other spectral approximants, viz. the pseudospectra and the numerical ranges, do converge

perfectly.The proof requires a modified approach to the finite section method for Wiener-Hopf operators. This note

generalizes results obtained by Böttcher, Grudsky and Silbermann for the case of continuous generating

functions.     

Numerical approximation for a Baer-Nunziato model of two-phase flows

We present a well-balanced numerical scheme for approximating the solution of the Baer-Nunziato model of two-phase flows by balancing the source terms and discretizing the compaction dynamics equation. First, the system is transformed into a new one of three subsystems: the first subsystem consists of the balance laws in the gas phase, the second subsystem consists of the conservation law of the mass in the solid phase and the conservation law of the momentum of the mixture, and the compaction dynamic equation is considered as the third subsystem. In the first subsystem, stationary waves are used to build up a well-balanced scheme which can capture equilibrium states. The second subsystem is of conservative form and thus can be numerically treated in a standard way. For the third subsystem, the fact that the solid velocity is constant across the solid contact suggests us to compose the technique of the Engquist-Osher scheme. We show that our scheme is capable of capturing exactly equilibrium states. Moreover, numerical tests show the convergence of approximate solutions to the exact solution.     

Finite element approximation of spectral problems with Neumann boundary conditions on curved domains

This paper deals with the finite element approximation of the spectral problem for the Laplace equation with Neumann boundary conditions on a curved nonconvex domain . Convergence and optimal order error estimates are proved for standard piecewise linear continuous elements on a discrete polygonal domain in the framework of the abstract spectral approximation theory.

     

On the steady flow of a viscous fluid past a sphere at small reynolds numbers using Oseen's approximation

Sommaire On a obtenu une expression générale pour la fonction de courant deStokes pour l'écoulement d'un liquide visqueux autour d'une sphère. Les lignes de courant ont été tracés pour les nombres deReynolds=1, 4 et 10. Le décollement du courant fluide n'a pas lieu pour les N. R.=1 et 4, tandis que pour le N. R.=10 les tourbillons attachés se forment dans le sillage.     

On the approximation of Laplacian eigenvalues in graph disaggregation

Graph disaggregation is a technique used to address the high cost of computation for power law graphs on parallel processors. The few high-degree vertices are broken into multiple small-degree vertices, in order to allow for more efficient computation in parallel. In particular, we consider computations involving the graph Laplacian, which has significant applications, including diffusion mapping and graph partitioning, among others. We prove results regarding the spectral approximation of the Laplacian of the original graph by the Laplacian of the disaggregated graph. In addition, we construct an alternate disaggregation operator whose eigenvalues interlace those of the original Laplacian. Using this alternate operator, we construct a uniform preconditioner for the original graph Laplacian.     

On a two-fluid model of two-phase compressible flows and its numerical approximation

We consider a two-fluid model of two-phase compressible flows. First, we derive several forms of the model and of the equations of state. The governing equations in all the forms contain source terms representing the exchanges of momentum and energy between the two phases. These source terms cause unstability for standard numerical schemes. Using the above forms of equations of state, we construct a stable numerical approximation for this two-fluid model. That only the source terms cause the oscillations suggests us to minimize the effects of source terms by reducing their amount. By an algebraic operator, we transform the system to a new one which contains only one source term. Then, we discretize the source term by making use of stationary solutions. We also present many numerical tests to show that while standard numerical schemes give oscillations, our scheme is stable and numerically convergent.     

Non‐Hermitian perturbations of Hermitian matrix‐sequences and applications to the spectral analysis of the numerical approximation of partial differential equations

This article concerns the spectral analysis of matrix‐sequences which can be written as a non‐Hermitian perturbation of a given Hermitian matrix‐sequence. The main result reads as follows. Suppose that for every n there is a Hermitian matrix X_n of size n and that {X_n}_n～_λf, that is, the matrix‐sequence {X_n}_n enjoys an asymptotic spectral distribution, in the Weyl sense, described by a Lebesgue measurable function f; if $\Vert Y_n{\Vert }_2=o(\sqrt{n})$ with ‖·‖₂ being the Schatten 2 norm, then {X_n+Y_n}_n～_λf. In a previous article by Leonid Golinskii and the second author, a similar result was proved, but under the technical restrictive assumption that the involved matrix‐sequences {X_n}_n and {Y_n}_n are uniformly bounded in spectral norm. Nevertheless, the result had a remarkable impact in the analysis of both spectral distribution and clustering of matrix‐sequences arising from various applications, including the numerical approximation of partial differential equations (PDEs) and the preconditioning of PDE discretization matrices. The new result considerably extends the spectral analysis tools provided by the former one, and in fact we are now allowed to analyze linear PDEs with (unbounded) variable coefficients, preconditioned matrix‐sequences, and so forth. A few selected applications are considered, extensive numerical experiments are discussed, and a further conjecture is illustrated at the end of the article.     

On existence of global solutions to the Navier-Stokes equations for compressible and viscous flows on the surface of a sphere

Following tecniques proposed by A. V. Khazikhov and V. A. Weigant in 1995, we prove the global, with respect to time, existence and uniqueness of the solution to the Navier-Stokes equations for a compressible, viscous and barotropic fluid which moves on the surface of a sphere. In obtaining the main estimates we make use of the Hodge decomposition and the generalized potential theory due to K. Kodaira.

---

Sunto Seguendo le tecniche proposte da A. V. Khazikhov and V. A. Weigant nel 1995, si prova l'esistenza ed unicità della soluzione per le equazioni di Navier-Stokes per un fluido comprimibile, viscoso e barotropico che si muova sulla superficie di una sfera. Nell'ottenere le principali stime si utilizzano la decomposizione di Hodge e la teoria del potenziale generalizzato, dovuta a K. Kodaira.