Researcher profile

Chen-Song Zhang

Chen-Song Zhang contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 published item(s)

preprint2026arXiv

R-PINN: Recovery-type a-posteriori estimator enhanced adaptive PINN

In recent years, with the advancements in machine learning and neural networks, algorithms using physics-informed neural networks (PINNs) to solve PDEs have gained widespread applications. While these algorithms are well-suited for a wide range of equations, they often exhibit suboptimal performance when applied to equations with large local gradients, resulting in substantial localized errors. To address this issue, this paper proposes an adaptive PINN algorithm designed to improve accuracy in such cases. The core idea of the algorithm is to adaptively adjust the distribution of collocation points based on the recovery-type a-posterior error of the current numerical solution, enabling a better approximation of the true solution. This approach is inspired by the adaptive finite element method. By combining the recovery-type a-posteriori estimator, a gradient-recovery estimator commonly used in the adaptive finite element method (FEM) with PINNs, we introduce the Recovery-type a-posteriori estimator enhanced adaptive PINN (R-PINN) and compare its performance with a typical adaptive PINN algorithm, FI-PINN. Our results demonstrate that R-PINN achieves faster convergence with fewer adaptive points and significantly outperforms in the cases with multiple regions of large errors than FI-PINN. Notably, our method is a hybrid numerical approach for solving partial differential equations, integrating adaptive FEM with PINNs.

preprint2022arXiv

A new analytical framework for the convergence of inexact two-grid methods

Two-grid methods with exact solution of the Galerkin coarse-grid system have been well studied by the multigrid community: an elegant identity has been established to characterize the convergence factor of exact two-grid methods. In practice, however, it is often too costly to solve the Galerkin coarse-grid system exactly, especially when its size is large. Instead, without essential loss of convergence speed, one may solve the coarse-grid system approximately. In this paper, we develop a new framework for analyzing the convergence of inexact two-grid methods: two-sided bounds for the energy norm of the error propagation matrix of inexact two-grid methods are presented. In the framework, a restricted smoother involved in the identity for exact two-grid convergence is used to measure how far the actual coarse-grid matrix deviates from the Galerkin one. As an application, we establish a unified convergence theory for multigrid methods.

preprint2020arXiv

Lowest-degree robust finite element scheme for a fourth-order elliptic singular perturbation problem on rectangular grids

In this paper, a piecewise quadratic nonconforming finite element method on rectangular grids for a fourth-order elliptic singular perturbation problem is presented. This proposed method is robustly convergent with respect to the perturbation parameter. Numerical results are presented to verify the theoretical findings. The new method uses piecewise quadratic polynomials, and is of the lowest degree possible. Optimal order approximation property of the finite element space is proved by means of a locally-averaged interpolation operator newly constructed. This interpolator, however, is not a projection. Indeed, we establish a general theory and show that no locally defined interpolation associated with the locally supported basis functions can be projective for the finite element space in use. Particularly, the general theory gives an answer to a long-standing open problem presented in [Demko, J. Approx. Theory, $\bf{43}$(2):151--156, 1985].

preprint2014arXiv

An adaptive hybrid stress transition quadrilateral finite element method for linear elasticity

In this paper, we discuss an adaptive hybrid stress finite element method on quadrilateral meshes for linear elasticity problems. To deal with hanging nodes arising in the adaptive mesh refinement, we propose new transition types of hybrid stress quadrilateral elements with 5 to 7 nodes. In particular, we derive a priori error estimation for the 5-node transition hybrid stress element to show that it is free from Poisson-locking, in the sense that the error bound in the a priori estimate is independent of the Lame constant $λ$. We introduce, for quadrilateral meshes, refinement/coarsening algorithms, which do not require storing the refinement tree explicitly, and give an adaptive algorithm. Finally we provide some numerical results.

preprint2013arXiv

An Error-Resilient Redundant Subspace Correction Method

As we stride toward the exascale era, due to increasing complexity of supercomputers, hard and soft errors are causing more and more problems in high-performance scientific and engineering computation. In order to improve reliability (increase the mean time to failure) of computing systems, a lot of efforts have been devoted to developing techniques to forecast, prevent, and recover from errors at different levels, including architecture, application, and algorithm. In this paper, we focus on algorithmic error resilient iterative linear solvers and introduce a redundant subspace correction method. Using a general framework of redundant subspace corrections, we construct iterative methods, which have the following properties: (1) Maintain convergence when error occurs assuming it is detectable; (2) Introduce low computational overhead when no error occurs; (3) Require only small amount of local (point-to-point) communication compared to traditional methods and maintain good load balance; (4) Improve the mean time to failure. With the proposed method, we can improve reliability of many scientific and engineering applications. Preliminary numerical experiments demonstrate the efficiency and effectiveness of the new subspace correction method.

preprint2013arXiv

Combined Preconditioning with Applications in Reservoir Simulation

We develop a simple algorithmic framework to solve large-scale symmetric positive definite linear systems. At its core, the framework relies on two components: (1) a norm-convergent iterative method (i.e. smoother) and (2) a preconditioner. The resulting preconditioner, which we refer to as a combined preconditioner, is much more robust and efficient than the iterative method and preconditioner when used in Krylov subspace methods. We prove that the combined preconditioner is positive definite and show estimates on the condition number of the preconditioned system. We combine an algebraic multigrid method and an incomplete factorization preconditioner to test the proposed framework on problems in petroleum reservoir simulation. Our numerical experiments demonstrate noticeable speed-up when we compare our combined method with the standalone algebraic multigrid method or the incomplete factorization preconditioner.

preprint2013arXiv

Numerical Study of Geometric Multigrid Methods on CPU--GPU Heterogeneous Computers

The geometric multigrid method (GMG) is one of the most efficient solving techniques for discrete algebraic systems arising from elliptic partial differential equations. GMG utilizes a hierarchy of grids or discretizations and reduces the error at a number of frequencies simultaneously. Graphics processing units (GPUs) have recently burst onto the scientific computing scene as a technology that has yielded substantial performance and energy-efficiency improvements. A central challenge in implementing GMG on GPUs, though, is that computational work on coarse levels cannot fully utilize the capacity of a GPU. In this work, we perform numerical studies of GMG on CPU--GPU heterogeneous computers. Furthermore, we compare our implementation with an efficient CPU implementation of GMG and with the most popular fast Poisson solver, Fast Fourier Transform, in the cuFFT library developed by NVIDIA.

preprint2012arXiv

A Scalable Auxiliary Space Preconditioner for High-Order Finite Element Methods

In this paper, we revisit an auxiliary space preconditioning method proposed by Xu [Computing 56, 1996], in which low-order finite element spaces are employed as auxiliary spaces for solving linear algebraic systems arising from high-order finite element discretizations. We provide a new convergence rate estimate and parallel implementation of the proposed algorithm. We show that this method is user-friendly and can play an important role in a variety of Poisson-based solvers for more challenging problems such as the Navier--Stokes equation. We investigate the performance of the proposed algorithm using the Poisson equation and the Stokes equation on 3D unstructured grids. Numerical results demonstrate the advantages of the proposed algorithm in terms of efficiency, robustness, and parallel scalability.

preprint2012arXiv

Self-sustaining oscillations of a falling sphere through Johnson-Segalman fluids

We confirm numerically that the Johnson-Segalman model is able to reproduce the continual oscillations of the falling sphere observed in some viscoelastic models. The empirical choice of parameters used in the Johnson-Segalman model is from the ones that show the non-monotone stress-strain relation of the steady shear flows of the model. The carefully chosen parameters yield continual, self-sustaining, (ir)regular and periodic oscillations of the speed for the falling sphere through the Johnson-Segalman fluids. In particular, our simulations reproduce the phenomena: the falling sphere settles slower and slower until a certain point at which the sphere suddenly accelerates and this pattern is repeated continually.