BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Xiangmin Jiao

Xiangmin Jiao appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysismath.APmath.DG

Catalog footprint

What is connected

11works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

ARPIST: Provably Accurate and Stable Numerical Integration over Spherical Triangles

Numerical integration on spheres, including the computation of the areas of spherical triangles, is a core computation in geomathematics. The commonly used techniques sometimes suffer from instabilities and significant loss of accuracy. We describe a new algorithm, called ARPIST, for accurate and stable integration of functions on spherical triangles. ARPIST is based on an easy-to-implement transformation to the spherical triangle from its corresponding linear triangle via radial projection to achieve high accuracy and efficiency. More importantly, ARPIST overcomes potential instabilities in computing the Jacobian of the transformation, even for poorly shaped triangles that may occur at poles in regular longitude-latitude meshes, by avoiding potential catastrophic rounding errors. We compare our proposed technique with L'Huilier's Theorem for computing the area of spherical triangles, and also compare it with the recently developed LSQST method (J. Beckmann, H.N. Mhaskar, and J. Prestin, GEM - Int. J. Geomath., 5:143-162, 2014) and a radial-basis-function-based technique (J. A. Reeger and B. Fornberg, Stud. Appl. Math., 137:174-188, 2015) for integration of smooth functions on spherical triangulations. Our results show that ARPIST enables superior accuracy and stability over previous methods while being orders of magnitude faster and significantly easier to implement.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Approximate Generalized Inverses with Iterative Refinement for $ε$-Accurate Preconditioning of Singular Systems

We introduce a new class of preconditioners to enable flexible GMRES to find a least-squares solution, and potentially the pseudoinverse solution, of large-scale sparse, asymmetric, singular, and potentially inconsistent systems. We develop the preconditioners based on a new observation that generalized inverses (i.e., $\boldsymbol{A}^{g}\in\{\boldsymbol{G}\mid\boldsymbol{A}\boldsymbol{G}\boldsymbol{A}=\boldsymbol{A}\}$) enable the preconditioned Krylov subspaces to converge in a single step. We then compute an approximate generalized inverse (AGI) efficiently using a hybrid incomplete factorization (HIF), which combines multilevel incomplete LU with rank-revealing QR on its final Schur complement. We define the criteria of $ε$-accuracy and stability of AGI to guarantee the convergence of preconditioned GMRES for consistent systems. For inconsistent systems, we fortify HIF with iterative refinement to obtain HIFIR, which allows accurate computations of the null-space vectors. By combining the two techniques, we then obtain a new solver, called PIPIT, for obtaining the pseudoinverse solutions for systems with low-dimensional null spaces. We demonstrate the robustness of HIF and HIFIR and show that they improve both accuracy and efficiency of the prior state of the art by orders of magnitude for systems with up to a million unknowns.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Optimal and Low-Memory Near-Optimal Preconditioning of Fully Implicit Runge-Kutta Schemes for Parabolic PDEs

Runge-Kutta (RK) schemes, especially Gauss-Legendre and some other fully implicit RK (FIRK) schemes, are desirable for the time integration of parabolic partial differential equations due to their A-stability and high-order accuracy. However, it is significantly more challenging to construct optimal preconditioners for them compared to diagonally implicit RK (or DIRK) schemes. To address this challenge, we first introduce mathematically optimal preconditioners called block complex Schur decomposition (BCSD), block real Schur decomposition (BRSD), and block Jordan form (BJF), motivated by block-circulant preconditioners and Jordan form solution techniques for IRK. We then derive an efficient, near-optimal singly-diagonal approximate BRSD (SABRSD) by approximating the quasi-triangular matrix in real Schur decomposition using an optimized upper-triangular matrix with a single diagonal value. A desirable feature of SABRSD is that it has comparable memory requirements and factorization (or setup) cost as singly DIRK (SDIRK). We approximate the diagonal blocks in these preconditioning techniques using an incomplete factorization with (near) linear complexity, such as multilevel ILU, ILU(0), or a multigrid method with an ILU-based smoother. We apply the block preconditioners in right-preconditioned GMRES to solve the advection-diffusion equation in 3D using finite element and finite difference methods. We show that BCSD, BRSD, and BJF significantly outperform other preconditioners in terms of GMRES iterations, and SABRSD is competitive with them and the prior state of the art in terms of computational cost while requiring the least amount of memory.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Hybrid Method and Unified Analysis of Generalized Finite Differences and Lagrange Finite Elements

Finite differences, finite elements, and their generalizations are widely used for solving partial differential equations, and their high-order variants have respective advantages and disadvantages. Traditionally, these methods are treated as different (strong vs. weak) formulations and are analyzed using different techniques (Fourier analysis or Green's functions vs. functional analysis), except for some special cases on regular grids. Recently, the authors introduced a hybrid method, called Adaptive Extended Stencil FEM or AES-FEM (Int. J. Num. Meth. Engrg., 2016, DOI:10.1002/nme.5246), which combines features of generalized finite differences and Lagrange finite elements to achieve second-order accuracy over unstructured meshes. However, its analysis was incomplete due to the lack of existing mathematical theory that unifies the formulations and analysis of these different methods. In this work, we introduce the framework of generalized weighted residuals to unify the formulation of finite differences, finite elements, and AES-FEM. In addition, we propose a unified analysis of the well-posedness, convergence, and mesh-quality dependency of these different methods. We also report numerical results with AES-FEM to verify our analysis. We show that AES-FEM improves the accuracy of generalized finite differences while reducing the mesh-quality dependency and simplifying the implementation of high-order finite elements.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

WLS-ENO Remap: Superconvergent and Non-Oscillatory Weighted Least Squares Data Transfer on Surfaces

Data remap between non-matching meshes is a critical step in multiphysics coupling using a partitioned approach. The data fields being transferred often have jumps in function values or derivatives. It is important but very challenging to avoid spurious oscillations (a.k.a. the Gibbs Phenomenon) near discontinuities and at the same time to achieve high-order accuracy away from discontinuities. In this work, we introduce a new approach, called WLS-ENOR, or Weighted-Least-Squares-based Essentially Non-Oscillatory Remap, to address this challenge. Based on the WLS-ENO reconstruction technique proposed by Liu and Jiao (J. Comput. Phys. vol 314, pp 749--773, 2016), WLS-ENOR differs from WLS-ENO and other WENO schemes in that it resolves not only the O(1) oscillations due to C 0 discontinuities, but also the accumulated effect of O(h) oscillations due to C 1 discontinuities. To this end, WLS-ENOR introduces a robust detector of discontinuities and a new weighting scheme for WLS-ENO near discontinuities. We also optimize the weights at smooth regions to achieve superconvergence. As a result, WLS-ENOR is more than fifth-order accurate and highly conservative in smooth regions, while being non-oscillatory and minimally diffusive near discontinuities. We also compare WLS-ENOR with some commonly used methods based on L 2 projection, moving least squares, and radial basis functions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

AP-Cloud: Adaptive Particle-in-Cloud Method for Optimal Solutions to Vlasov-Poisson Equation

We propose a new adaptive Particle-in-Cloud (AP-Cloud) method for obtaining optimal numerical solutions to the Vlasov-Poisson equation. Unlike the traditional particle-in-cell (PIC) method, which is commonly used for solving this problem, the AP-Cloud adaptively selects computational nodes or particles to deliver higher accuracy and efficiency when the particle distribution is highly non-uniform. Unlike other adaptive techniques for PIC, our method balances the errors in PDE discretization and Monte Carlo integration, and discretizes the differential operators using a generalized finite difference (GFD) method based on a weighted least square formulation. As a result, AP-Cloud is independent of the geometric shapes of computational domains and is free of artificial parameters. Efficient and robust implementation is achieved through an octree data structure with 2:1 balance. We analyze the accuracy and convergence order of AP-Cloud theoretically, and verify the method using an electrostatic problem of a particle beam with halo. Simulation results show that the AP-Cloud method is substantially more accurate and faster than the traditional PIC, and it is free of artificial forces that are typical for some adaptive PIC techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Overcoming Element Quality Dependence of Finite Elements with Adaptive Extended Stencil FEM (AES-FEM)

The finite element methods (FEM) are important techniques in engineering for solving partial differential equations, but they depend heavily on element shape quality for stability and good performance. In this paper, we introduce the Adaptive Extended Stencil Finite Element Method (AES-FEM) as a means for overcoming this dependence on element shape quality. Our method replaces the traditional basis functions with a set of generalized Lagrange polynomial (GLP) basis functions, which we construct using local weighted least-squares approximations. The method preserves the theoretical framework of FEM, and allows imposing essential boundary conditions and integrating the stiffness matrix in the same way as the classical FEM. In addition, AES-FEM can use higher-degree polynomial basis functions than the classical FEM, while virtually preserving the sparsity pattern of the stiffness matrix. We describe the formulation and implementation of AES-FEM, and analyze its consistency and stability. We present numerical experiments in both 2D and 3D for the Poison equation and a time-independent convection-diffusion equation. The numerical results demonstrate that AES-FEM is more accurate than linear FEM, is also more efficient than linear FEM in terms of error versus runtime, and enables much better stability and faster convergence of iterative solvers than linear FEM over poor-quality meshes

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

WLS-ENO: Weighted-Least-Squares Based Essentially Non-Oscillatory Schemes for Finite Volume Methods on Unstructured Meshes

ENO (Essentially Non-Oscillatory) and WENO (Weighted Essentially Non-Oscillatory) schemes are widely used high-order schemes for solving partial differential equations (PDEs), especially hyperbolic conservation laws with piecewise smooth solutions. For structured meshes, these techniques can achieve high order accuracy for smooth functions while being non-oscillatory near discontinuities. For unstructured meshes, which are needed for complex geometries, similar schemes are required but they are much more challenging. We propose a new family of non-oscillatory schemes, called WLS-ENO, in the context of solving hyperbolic conservation laws using finite-volume methods over unstructured meshes. WLS-ENO is derived based on Taylor series expansion and solved using a weighted least squares formulation. Unlike other non-oscillatory schemes, the WLS-ENO does not require constructing sub-stencils, and hence it provides more flexible framework and is less sensitive to mesh quality. We present rigorous analysis of the accuracy and stability of WLS-ENO, and present numerical results in 1-D, 2-D, and 3-D for a number of benchmark problems, and also report some comparisons against WENO.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

OPINS: An Orthogonally Projected Implicit Null-space Method for Singular and Nonsingular Saddle-point Systems

Saddle-point systems appear in many scientific and engineering applications. The systems can be sparse, symmetric or nonsymmetric, and possibly singular. In many of these applications, the number of constraints is relatively small compared to the number of unknowns. The traditional null-space method is inefficient for these systems, because it is expensive to find the null space explicitly. Some alternatives, notably constraint-preconditioned or projected Krylov methods, are relatively efficient, but they can suffer from numerical instability and even nonconvergence. In addition, most existing methods require the system to be nonsingular or be reducible to a nonsingular system. In this paper, we propose a new method, called OPINS, for singular and nonsingular saddle-point systems. OPINS is equivalent to the null-space method with an orthogonal projector, without forming the orthogonal basis of the null space explicitly. OPINS can not only solve for the unique solution for nonsingular saddle-point problems, but also find the minimum-norm solution in terms of the solution variables for singular systems. The method is efficient and easy to implement using existing Krylov solvers for singular systems. At the same time, it is more stable than the other alternatives, such as projected Krylov methods. We present some preconditioners to accelerate the convergence of OPINS for nonsingular systems, and compare OPINS against some present state-of-the-art methods for various types of singular and nonsingular systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An Embedded Boundary Method for Two Phase Incompressible Flow

We develop an embedded boundary method (EBM) to solve the two-phase incompressible flow with piecewise constant density. The front tracking method is used to track the interface. The fractional step methods are used to solve the incompressible Navier-Stokes equations while the EBM is used in the projection step to solve an elliptic interface problem for the pressure with a jump equal to the surface tension force across the interface. Several examples are used to verify the accuracy of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Consistent Computation of First- and Second-Order Differential Quantities for Surface Meshes

Differential quantities, including normals, curvatures, principal directions, and associated matrices, play a fundamental role in geometric processing and physics-based modeling. Computing these differential quantities consistently on surface meshes is important and challenging, and some existing methods often produce inconsistent results and require ad hoc fixes. In this paper, we show that the computation of the gradient and Hessian of a height function provides the foundation for consistently computing the differential quantities. We derive simple, explicit formulas for the transformations between the first- and second-order differential quantities (i.e., normal vector and principal curvature tensor) of a smooth surface and the first- and second-order derivatives (i.e., gradient and Hessian) of its corresponding height function. We then investigate a general, flexible numerical framework to estimate the derivatives of the height function based on local polynomial fittings formulated as weighted least squares approximations. We also propose an iterative fitting scheme to improve accuracy. This framework generalizes polynomial fitting and addresses some of its accuracy and stability issues, as demonstrated by our theoretical analysis as well as experimental results.

0 citations0 reviewsNo code linkedOpen access