The paper presents results on piecewise polynomial approximations of tensor product type in Sobolev-Slobodecki spaces by various interpolation and projection techniques, on error estimates for quadrature rules and projection operators based on hierarchical bases, and on inverse inequalities.

In this paper we present a new computationally efficient numerical scheme for the minimizing flow for the computation of the optimal L2 mass transport mapping using the fluid approach. We review the method and discuss its numerical properties. We then derive a new scaleable, efficient discretization and a solution technique for the problem and show that the problem is equivalent to a mixed form formulation of a nonlinear fluid flow in porous media. We demonstrate the effectiveness of our approach using a number of numerical experiments.

An efficient multigrid solver for the Oseen problems discretized by Marker and Cell (MAC) scheme on staggered grid is developed in this paper. Least squares commutator distributive Gauss-Seidel (LSC-DGS) relaxation is generalized and developed for Oseen problems. Residual overweighting technique is applied to further improve the performance of the solver and a defect correction method is suggested to improve the accuracy of the discretization. Some numerical results are presented to demonstrate the efficiency and robustness of the proposed solver.

We analyze here, a two-grid finite element method for the two dimensional time-dependent incompressible Navier-Stokes equations with non-smooth initial data. It involves solving the non-linear Navier-Stokes problem on a coarse grid of size H and solving a Stokes problem on a fine grid of size h, h « H. This method gives optimal convergence for velocity in H1-norm and for pressure in L2-norm. The analysis mainly focuses on the loss of regularity of the solution at t = 0 of the Navier-Stokes equations.

This paper provides an overview of the main ideas driving the bootstrap algebraic multigrid methodology, including compatible relaxation and algebraic distances for defining effective coarsening strategies, the least squares method for computing accurate prolongation operators and the bootstrap cycles for computing the test vectors that are used in the least squares process. We review some recent research in the development, analysis and application of bootstrap algebraic multigrid and point to open problems in these areas. Results from our previous research as well as some new results for some model diffusion problems with highly oscillatory diffusion coefficient are presented to illustrate the basic components of the BAMG algorithm.

This paper proposes and analyzes semi-discrete and fully discrete hybrid stress finite element methods for elastodynamic problems. A hybrid stress quadrilateral finite element approximation is used in the space directions. A second-order center difference is adopted in the time direction for the fully discrete scheme. Error estimates of the two schemes, as well as a stability result for the fully discrete scheme, are derived. Numerical experiments are done to verify the theoretical results.

A finite-element multigrid scheme for elliptic Nash-equilibrium multiobjective optimal control problems with control constraints is investigated. The multigrid computational framework implements a nonlinear multigrid strategy with collective smoothing for solving the multiobjective optimality system discretized with finite elements. Error estimates for the optimal solution and two-grid local Fourier analysis of the multigrid scheme are presented. Results of numerical experiments are presented to demonstrate the effectiveness of the proposed framework.

Algebraic Multigrid (AMG) methods were developed originally for numerically solving Partial Differential Equations (PDE), not necessarily on structured grids. In the last two decades solvers inspired by the AMG approach, were developed for non PDE problems, including data and image analysis problems, such as clustering, segmentation, quantization and others. These solvers share a common principle in that there is a crosstalk between fine and coarse representations of the problems, with flow of information in both directions, fine-to-coarse and coarse-to-fine. This paper surveys some of these problems and the AMG-inspired algorithms for their solution.

In this paper we present a fully discrete A-ø finite element method to solve Maxwell’s equations with a nonlinear degenerate boundary condition, which represents a generalization of the classical Silver-Müller condition for a non-perfect conductor. The relationship between the normal components of the electric field E and the magnetic field H obeys a power-law nonlinearity of the type H x n = n x (|E x n|α-1E x n) with α ∈ (0,1]. We prove the existence and uniqueness of the solutions of the proposed A-ø scheme and derive the error estimates. Finally, we present some numerical experiments to verify the theoretical result.

A shifted Laplacian operator is obtained from the Helmholtz operator by adding a complex damping. It serves as a basic tool in the most successful multigrid approach for solving highly indefinite Helmholtz equations — a Shifted Laplacian preconditioner for Krylov-type methods. Such preconditioning significantly accelerates Krylov iterations, much more so than the multigrid based on original Helmholtz equations. In this paper, we compare approximation and relaxation properties of the Helmholtz operator with and without the complex shift, and, based on our observations, propose a new hybrid approach that combines the two. Our analytical conclusions are supported by two-dimensional numerical results.

We present a conservative formulation and a numerical algorithm for the reduced-gravity shallow-water equations on a beta plane, subjected to a constant wind forcing that leads to the formation of double-gyre circulation in a closed ocean basin. The novelty of the paper is that we reformulate the governing equations into a nonlinear hyperbolic conservation law plus source terms. A second-order fractional-step algorithm is used to solve the reformulated equations. In the first step of the fractional-step algorithm, we solve the homogeneous hyperbolic shallow-water equations by the wave-propagation finite volume method. The resulting intermediate solution is then used as the initial condition for the initial-boundary value problem in the second step. As a result, the proposed method is not sensitive to the choice of viscosity and gives high-resolution results for coarse grids, as long as the Rossby deformation radius is resolved. We discuss the boundary conditions in each step, when no-slip boundary conditions are imposed to the problem. We validate the algorithm by a periodic flow on an f-plane with exact solutions. The order-of-accuracy for the proposed algorithm is tested numerically. We illustrate a quasi-steady-state solution of the double-gyre model via the height anomaly and the contour of stream function for the formation of double-gyre circulation in a closed basin. Our calculations are highly consistent with the results reported in the literature. Finally, we present an application, in which the double-gyre model is coupled with the advection equation for modeling transport of a pollutant in a closed ocean basin.