We present a brief survey on (Weakly) Admissible Meshes and corresponding Discrete Extremal Sets, namely Approximate Fekete Points and Discrete Leja Points. These provide new computational tools for polynomial least squares and interpolation on multidimensional compact sets, with different applications such as numerical cubature, digital filtering, spectral and high-order methods for PDEs.

We develop and analyze an adaptive hybridized Interior Penalty Discontinuous Galerkin (IPDG-H) method for H(curl)-elliptic boundary value problems in 2D or 3D arising from a semi-discretization of the eddy currents equations. The method can be derived from a mixed formulation of the given boundary value problem and involves a Lagrange multiplier that is an approximation of the tangential traces of the primal variable on the interfaces of the underlying triangulation of the computational domain. It is shown that the IPDG-H technique can be equivalently formulated and thus implemented as a mortar method. The mesh adaptation is based on a residual-type a posteriori error estimator consisting of element and face residuals. Within a unified framework for adaptive finite element methods, we prove the reliability of the estimator up to a consistency error. The performance of the adaptive symmetric IPDG-H method is documented by numerical results for representative test examples in 2D.

This paper is concerned with the numerical approximations of semi-linear stochastic partial differential equations of elliptic type in multi-dimensions. Convergence analysis and error estimates are presented for the numerical solutions based on the spectral method. Numerical results demonstrate the good performance of the spectral method.

We provided in [14] and [15] a semilocal convergence analysis for Newton’s method on a Banach space setting, by splitting the given operator. In this study, we improve the error bounds, order of convergence, and simplify the sufficient convergence conditions. Our results compare favorably with the Newton-Kantorovich theorem for solving equations.

In this paper we discuss the extension to exponential splitting methods with respect to time-dependent operators. We concentrate on the Suzuki’s method, which incorporates ideas to the time-ordered exponential of [3,11,12,34]. We formulate the methods with respect to higher order by using kernels for an extrapolation scheme. The advantages include more accurate and less computational intensive schemes to special time-dependent harmonic oscillator problems. The benefits of the higher order kernels are given on different numerical examples.

An efficient finite difference framework based on moving meshes methods is developed for the three-dimensional free surface viscoelastic flows. The basic model equations are based on the incompressible Navier-Stokes equations and the Oldroyd-B constitutive model for viscoelastic flows is adopted. A logical domain semi-Lagrangian scheme is designed for moving-mesh solution interpolation and convection. Numerical results show that harmonic map based moving mesh methods can achieve better accuracy for viscoelastic flows with free boundaries while using much less memory and computational time compared to the uniform mesh simulations.