No abstract uploaded!

We derive bounds for eigenvalues of the Laplacian of graphs using discrete versions of the Sobolev inequalities and heat kernel estimates.

In this paper, a structure-preserving algorithm is developed for the computation of a semi-stabilizing solution of a Generalized Algebraic Riccati Equation (GARE). The semi-stabilizing solution of GAREs has been used to characterize the solvability of the (J, J ′ )-spectral factorization problem in control theory for general rational matrices which may have poles and zeros on the extended imaginary axis. The main difficulty in solving such a GARE lies in the fact that its associated Hamiltonian/skew-Hamiltonian pencil has eigenvalues on the extended imaginary axis. Consequently, it is not clear which eigenspace of the associated Hamiltonian/skew-Hamiltonian pencil can characterize the desired semi-stabilizing solution. That is, it is not clear which eigenvectors and principal vectors corresponding to the eigenvalues on the extended imaginary axis should be contained in the eigenspace that we wish to compute. Hence, the well-known generalized eigenspace approach for the classical algebraic Riccati equations cannot be employed directly. The proposed algorithm consists of a structure-preserving doubling algorithm (SDA) and a postprocessing procedure to determine the desired eigenvectors and principal vectors corresponding to the purely imaginary and infinite eigenvalues. Under mild assumptions, linear convergence of rate 1/2 for the SDA is proved. Numerical experiments illustrate that the proposed algorithm performs efficiently and reliably.

For decades, the widely used finite difference method on staggered grids, also known as the marker and cell (MAC) method, has been one of the simplest and most effective numerical schemes for solving the Stokes equations and Navier–Stokes equations. Its superconvergence on uniform meshes has been observed by Nicolaides (SIAM J Numer Anal 29(6):1579–1591, 1992), but the rigorous proof is never given. Its behavior on non-uniform grids is not well studied, since most publications only consider uniform grids. In this work, we develop the MAC scheme on non-uniform rectangular meshes, and for the first time we theoretically prove that the superconvergence phenomenon (i.e., second order convergence in the L2 norm for both velocity and pressure) holds true for the MAC method on non-uniform rectangular meshes. With a careful and accurate analysis of various sources of errors, we observe that even though the local truncation errors are only first order in terms of mesh size, the global errors after summation are second order due to the amazing cancellation of local errors. This observation leads to the elegant superconvergence analysis even with non-uniform meshes. Numerical results are given to verify our theoretical analysis.

No abstract uploaded!