Numerical Analysis (math.NA)

  • PDF
    We introduce a cousin of the DPG method - the DPG* method - discuss their relationship and compare the two methods through numerical experiments.
  • PDF
    In this work, we propose a new Galerkin-Petrov method for the numerical solution of the classical spatially homogeneous Boltzmann equation. This method is based on an approximation of the distribution function by associated Laguerre polynomials and spherical harmonics and test an a variational manner with globally defined three-dimensional polynomials. A numerical realization of the algorithm is presented. The algorithmic developments are illustrated with the help of several numerical tests.
  • PDF
    We study an iterative low-rank approximation method for the solution of the steady-state stochastic Navier--Stokes equations with uncertain viscosity. The method is based on linearization schemes using Picard and Newton iterations and stochastic finite element discretizations of the linearized problems. For computing the low-rank approximate solution, we adapt the nonlinear iterations to an inexact and low-rank variant, where the solution of the linear system at each nonlinear step is approximated by a quantity of low rank. This is achieved by using a tensor variant of the GMRES method as a solver for the linear systems. We explore the inexact low-rank nonlinear iteration with a set of benchmark problems, using a model of flow over an obstacle, under various configurations characterizing the statistical features of the uncertain viscosity, and we demonstrate its effectiveness by extensive numerical experiments.
  • PDF
    In this paper, we develop a stochastic Asymptotic-Preserving (sAP) scheme for the kinetic chemotaxis system with random inputs, which will converge to the modified Keller-Segel model with random inputs in the diffusive regime. Based on the generalized Polynomial Chaos (gPC) approach, we design a high order stochastic Galerkin method using implicit-explicit (IMEX) Runge-Kutta (RK) time discretization with a macroscopic penalty term. The new schemes improve the parabolic CFL condition to a hyperbolic type when the mean free path is small, which shows significant efficiency especially in uncertainty quantification (UQ) with multi-scale problems. The stochastic Asymptotic-Preserving property will be shown asymptotically and verified numerically in several tests. Many other numerical tests are conducted to explore the effect of the randomness in the kinetic system, in the aim of providing more intuitions for the theoretic study of the chemotaxis models.
  • PDF
    Finite element meshing of geometric domains by affine simplices results in polytopial approximation of these domains. The attendant error in the representation of the domain limits the accuracy of finite element methods based on such meshes. Yet, simplicity of affine meshes makes them a desirable modeling tool in many applications. In this paper we develop and analyze finite element methods that remain stable and optimally accurate on polytopial approximations of domains with smooth boundaries. This is achieved by constraining the extension of finite element solution on the polytopial domain to weakly match the prescribed boundary condition on the true geometric boundary. We provide numerical examples illustrating key properties of the new method, including optimal convergence rates in $H^1$ and $L^2$-norm.
  • PDF
    The Bayesian perspective on inverse problems has attracted much mathematical attention in recent years. Particular attention has been paid to Bayesian inverse problems (BIPs) in which the parameter to be inferred lies in an infinite-dimensional space, a typical example being a scalar or tensor field coupled to some observed data via an ODE or PDE. This article gives an introduction to the framework of well-posed BIPs in infinite-dimensional parameter spaces, as advocated by Stuart (Acta Numer. 19:451--559, 2010) and others. This framework has the advantage of ensuring uniformly well-posed inference problems independently of the finite-dimensional discretisation used for numerical solution. Recently, this framework has been extended to the case of a heavy-tailed prior measure in the family of stable distributions, such as an infinite-dimensional Cauchy distribution, for which polynomial moments are infinite or undefined. It is shown that analogues of the Karhunen--Loève expansion for square-integrable random variables can be used to sample such measures on quasi-Banach spaces. Furthermore, under weaker regularity assumptions than those used to date, the Bayesian posterior measure is shown to depend Lipschitz continuously in the Hellinger and total variation metrics upon perturbations of the misfit function and observed data.
  • PDF
    We present an adaptation of the MA-LBR scheme to the Monge-Ampère equation with second boundary value condition, provided the target is a convex set. This yields a fast adaptive method to numerically solve the Optimal Transport problem between two absolutely continuous measures, the second of which has convex support. The proposed numerical method actually captures a specific Brenier solution which is minimal in some sense. We prove the convergence of the method as the grid stepsize vanishes and we show with numerical experiments that it is able to reproduce subtle properties of the Optimal Transport problem.
  • PDF
    In this work, we prove rigorous convergence properties for a semi-discrete, moment-based approximation of a model kinetic equation in one dimension. This approximation is equivalent to a standard spectral method in the velocity variable of the kinetic distribution and, as such, is accompanied by standard algebraic estimates of the form $N^{-q}$, where $N$ is the number of modes and $q$ depends on the regularity of the solution. However, in the multiscale setting, we show that the error estimate can be expressed in terms of the scaling parameter $\epsilon$, which measures the ratio of the mean-free-path to the characteristic domain length. In particular, we show that the error in the spectral approximation is $\mathcal{O}(\epsilon^{N+1})$. More surprisingly, for isotropic initial conditions, the coefficients of the expansion satisfy super convergence properties. In particular, the error of the $\ell^{th}$ coefficient of the expansion scales like $\mathcal{O}(\epsilon^{2N})$ when $\ell =0$ and $\mathcal{O}(\epsilon^{2N+2-\ell})$ for all $1\leq \ell \leq N$. This result is significant, because the low-order coefficients correspond to physically relevant quantities of the underlying system. All the above estimates involve constants depending on $N$, the time $t$, and the initial condition. We investigate specifically the dependence on $N$, in order to assess whether increasing $N$ actually yields an additional factor of $\epsilon$ in the error. Numerical tests will also be presented to support the theoretical results.
  • PDF
    The modified biharmonic equation is encountered in a variety of application areas, including streamfunction formulations of the Navier-Stokes equations. We develop a separation of variables representation for this equation in polar coordinates, for either the interior or exterior of a disk, and derive a new class of special functions which makes the approach stable. We discuss how these functions can be used in conjunction with fast algorithms to accelerate the solution of the modified biharmonic equation or the "bi-Helmholtz" equation in more complex geometries.
  • PDF
    Acceleration schemes can dramatically improve existing optimization procedures. In most of the work on these schemes, such as nonlinear GMRES, acceleration is based on minimizing the $\ell_2$ norm of some target on subspaces of $\mathbb{R}^n$. There are many numerical examples that show how accelerating general purpose and domain-specific optimizers with N-GMRES results in large improvements. We propose a natural modification to N-GMRES, which significantly improves the performance in a testing environment originally used to advocate N-GMRES. Our proposed approach, which we refer to as O-ACCEL, is novel in that it minimizes an approximation to the \emphobjective function on subspaces of $\mathbb{R}^n$. A convergence theorem is proved for our proposed approach. Comparisons with L-BFGS and N-CG indicate the competitiveness of O-ACCEL. As it can be combined with domain-specific optimizers, it may also be beneficial in areas where L-BFGS or N-CG are not suitable.
  • PDF
    In this article, we solve a deterministically generalized interpolation problem by a stochastic approach. We introduce a kernel-based probability measure on a Banach space by a covariance kernel which is defined on the dual space of the Banach space. The kernel-based probability measure provides a numerical tool to construct and analyze the kernel-based estimators conditioned on non-noise data or noisy data including algorithms and error analysis. Same as meshfree methods, we can also obtain the kernel-based approximate solutions of elliptic partial differential equations by the kernel-based probability measure.
  • PDF
    We consider computing the $k$-th eigenvalue and its corresponding eigenvector of large sparse generalized Hermitian eigenvalue problems. In electronic structure calculations, several properties of materials, such as those of optoelectronic device materials, are governed by the eigenpair with a material-specific index $k.$ Building on spectral bisection, we present a framework for computing the $k$-th eigenpair with validation of its index. Numerical results with problem sizes up to 1.5 million are reported, and the results demonstrate the accuracy and efficiency of the proposed framework.
  • PDF
    Nonlocal gradient operators are basic elements of nonlocal vector calculus that play important roles in nonlocal modeling and analysis. In this work, we extend earlier analysis on nonlocal gradient operators. In particular, we study a nonlocal Dirichlet integral that is given by a quadratic energy functional based on nonlocal gradients. Our main finding, which differs from claims made in previous studies, is that the coercivity and stability of this nonlocal continuum energy functional may hold for some properly chosen nonlocal interaction kernels but may fail for some other ones. This can be significant for possible applications of nonlocal gradient operators in various nonlocal models. In particular, we discuss some important implications for the peridynamic correspondence material models.
  • PDF
    A triangular membrane element with a new multi-resolution analysis (MRA) is proposed and a new finite element method is hence presented. The MRA framework is formulated out of a displacement subspace sequence whose basis functions are built out of scaling and shifting on the element domain of basic full node shape function. The basic full node shape function is constructed by means of extending the shape function triangle domain for a split node at the zero coordinates to the hexagon area enclosing the zero coordinates. As a result, a new split-full node notion is presented and a novel rational MRA concept together with the resolution level (RL) is constituted for the element. Via practical examples, it is found that the traditional triangular membrane element and method is a mono-resolution one and also a special case of the proposed element and method. The meshing for the monoresolution membrane element model is based on the empiricism while the RL adjusting for the multiresolution is laid on the rigorous mathematical basis. The analysis clarity of a plane-stress structure is actually determined by the RL, not by the mesh. Thus, the accuracy of a structural analysis is replaced by the clarity, the irrational MRA by the rational and the mesh model by the RL that is the discretized model by the integrated. The continuity of the full shape function unveils secrets behind assembling artificially of node-related items in global matrix formation by the conventional FEM.
  • PDF
    This paper presents a nonnegative polynomial that cannot be represented with nonnegative coefficients in the simplicial Bernstein basis by subdividing the standard simplex. The example shows that Bernstein Theorem cannot be extended to certificates of nonnegativity for polynomials with zeros at isolated points.
  • PDF
    Hyperspectral imaging is a cutting-edge type of remote sensing used for mapping vegetation properties, rock minerals and other materials. A major drawback of hyperspectral imaging devices is their intrinsic low spatial resolution. In this paper, we propose a method for increasing the spatial resolution of a hyperspectral image by fusing it with an image of higher spatial resolution that was obtained with a different imaging modality. This is accomplished by solving a variational problem in which the regularization functional is the directional total variation. To accommodate for possible mis-registrations between the two images, we consider a non-convex blind super-resolution problem where both a fused image and the corresponding convolution kernel are estimated. Using this approach, our model can realign the given images if needed. Our experimental results indicate that the non-convexity is negligible in practice and that reliable solutions can be computed using a variety of different optimization algorithms. Numerical results on real remote sensing data from plant sciences and urban monitoring show the potential of the proposed method and suggests that it is robust with respect to the regularization parameters, mis-registration and the shape of the kernel.