Related papers: Convergence and round-off errors in a two-dimensio…
In this work, we determine the full expression for the global truncation error of hyperbolic partial differential equations (PDEs). In particular, we use theoretical analysis and symbolic algebra to find exact expressions for the…
The first step when solving an infinite-dimensional eigenvalue problem is often to discretize it. We show that one must be extremely careful when discretizing nonlinear eigenvalue problems. Using examples, we show that discretization can:…
In two and three dimensional Lipschitz, but not necessarily convex, polytopal domains, we propose and analyze a posteriori error estimators for an optimal control problem involving the stationary Navier--Stokes equations; control…
Due to discretization effects and truncation to finite domains, many electromagnetic simulations present non-physical modifications of Maxwell's equations in space that may generate spurious signals affecting the overall accuracy of the…
In this paper, conditional stability estimates are derived for unique continuation and Cauchy problems associated to the Poisson equation in ultra-weak variational form. Numerical approximations are obtained as minima of regularized least…
We study regularity and numerical methods for two-sided fractional diffusion equations with a lower-order term. We show that the regularity of the solution in weighted Sobolev spaces can be greatly improved compared to that in standard…
We present a new approach to compute selected eigenvalues and eigenvectors of the two-parameter eigenvalue problem. Our method requires computing generalized eigenvalue problems of the same size as the matrices of the initial two-parameter…
High order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution $\Delta x$) with only a moderate increase to computational expense. Significant…
The Nikiforov-Uvarov method is a simple, yet elegant and powerful method for solving second-order differential equations of generalized hypergeometric type. In the past, it has been used to solve many problems in quantum mechanics and…
In this paper, we investigate the AMLI-cycle method and make two contributions. First, we revisit the AMLI-cycle using the Chebyshev polynomials and establish a theory for its uniform convergence, assuming the two-grid method converges…
In this paper, we consider the initial boundary value problem of the two dimensional multi-term time fractional mixed diffusion and diffusion-wave equations. An alternating direction implicit (ADI) spectral method is developed based on…
We consider two-point non-self-adjoint boundary eigenvalue problems for linear matrix differential operators. The coefficient matrices in the differential expressions and the matrix boundary conditions are assumed to depend analytically on…
We are interested in the simulation of open quantum systems governed by the Lindblad master equation in an infinite-dimensional Hilbert space. To simulate the solution of this equation, the standard approach involves two sequential…
We consider a class of shell models of 2D-turbulence. They conserve inertially the analogues of energy and enstrophy, two quadratic forms in the shell amplitudes. Inertially conserving two quadratic integrals leads to two spectral ranges.…
We construct one and two parameter deformations of the two dimensional Chebyshev polynomials with simple recurrence coefficients, following the algorithm in [3]. Using inverse scattering techniques, we compute the corresponding…
A spectral method is developed for the direct solution of linear ordinary differential equations with variable coefficients. The method leads to matrices which are almost banded, and a numerical solver is presented that takes O(m^2n)…
The proposed two-dimensional geometrically exact beam element extends our previous work by including the effects of shear distortion, and also of distributed forces and moments acting along the beam. The general flexibility-based…
We study bifurcation behavior in periodic perturbations of two-dimensional symmetric systems exhibiting codimension-two bifurcations with a double eigenvalue when the frequencies of the perturbation terms are small. We transform the…
In many contexts the modal properties of a structure change, either due to the impact of a changing environment, fatigue, or due to the presence of structural damage. For example during flight, an aircraft's modal properties are known to…
Non-smoothness at optimal points is a common phenomenon in many eigenvalue optimization problems. We consider two recent algorithms to minimize the largest eigenvalue of a Hermitian matrix dependent on one parameter, both proven to be…