BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Shi Shu

Shi Shu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath-phmath.MPmath.OC

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
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

7 published item(s)

preprint2022arXiv

Error Analysis of Virtual Element Methods for the Time-dependent Poisson-Nernst-Planck Equations

We discuss and analyze the virtual element method on general polygonal meshes for the time-dependent Poisson-Nernst-Planck equations, which are a nonlinear coupled system widely used in semiconductors and ion channels. The spatial discretization is based on the elliptic projection and the $L^2$ projection operator, and for the temporal discretization, the backward Euler scheme is employed. After presenting the semi and fully discrete schemes, we derive the a priori error estimates in the $L^2$ and $H^1$ norms. Finally, a numerical experiment verifies the theoretical convergence results.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Parallel Multi-Stage Preconditioners with Adaptive Setup for the Black Oil Model

The black oil model is widely used to describe multiphase porous media flow in the petroleum industry. The fully implicit method features strong stability and weak constraints on time step-sizes; hence, commonly used in the current mainstream commercial reservoir simulators. In this paper, a CPR-type preconditioner with an adaptive "setup phase" is developed to improve parallel efficiency of petroleum reservoir simulation. Furthermore, we propose a multi-color Gauss-Seidel (GS) algorithm for algebraic multigrid method based on the coefficient matrix of strong connections. Numerical experiments show that the proposed preconditioner can improve the parallel performance for both OpenMP and CUDA implements. Moreover, the proposed algorithm yields good parallel speedup as well as same convergence behavior as the corresponding single-threaded algorithm. In particular, for a three-phase benchmark problem, the parallel speedup of the OpenMP version is over 6.5 with 16 threads and the CUDA version reaches more than 9.5.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A multigrid-reduction-in-time solver with a new two-level convergence for unsteady fractional Laplacian problems

The multigrid-reduction-in-time (MGRIT) technique has proven to be successful in achieving higher run-time speedup by exploiting parallelism in time. The goal of this article is to develop and analyze a MGRIT algorithm, using FCF-relaxation with time-dependent time-grid propagators, to seek the finite element approximations of unsteady fractional Laplacian problems. The multigrid with line smoother proposed in [L. Chen, R. H. Nochetto, E. Ot{á}rola, A. J. Salgado, Math. Comp. 85 (2016) 2583--2607] is chosen to be the spatial solver. Motivated by [B. S. Southworth, SIAM J. Matrix Anal. Appl. 40 (2019) 564--608], we provide a new temporal eigenvalue approximation property and then deduce a generalized two-level convergence theory which removes the previous unitary diagonalization assumption on the fine and coarse time-grid propagators required in [X. Q. Yue, S. Shu, X. W. Xu, W. P. Bu, K. J. Pan, Comput. Math. Appl. 78 (2019) 3471--3484]. Numerical computations are included to confirm the theoretical predictions and demonstrate the sharpness of the derived convergence upper bound.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Dynamic Subspace Based BFGS Method for Large Scale Optimization Problem

Large-scale unconstrained optimization is a fundamental and important class of, yet not well-solved problems in numerical optimization. The main challenge in designing an algorithm is to require a few storage locations or very inexpensive computations while preserving global convergence. In this work, we propose a novel approach solving large-scale unconstrained optimization problem by combining the dynamic subspace technique and the BFGS update algorithm. It is clearly demonstrated that our approach has the same rate of convergence in the dynamic subspace as the BFGS and less memory than L-BFGS. Further, we give the convergence analysis by constructing the mapping of low-dimensional Euclidean space to the adaptive subspace. We compare our hybrid algorithm with the BFGS and L-BFGS approaches. Experimental results show that our hybrid algorithm offers several significant advantages such as parallel computing, convergence efficiency, and robustness.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Adaptive-Multilevel BDDC algorithm for three-dimensional plane wave Helmholtz systems

In this paper, we are concerned with the weighted plane wave least-squares (PWLS) method for three-dimensional Helmholtz equations, and develop the multi-level adaptive BDDC algorithms for solving the resulting discrete system. In order to form the adaptive coarse components, the local generalized eigenvalue problems for each common face and each common edge are carefully designed. The condition number of the two-level adaptive BDDC preconditioned system is proved to be bounded above by a user-defined tolerance and a constant which is dependent on the maximum number of faces and edges per subdomain and the number of subdomains sharing a common edge. The efficiency of these algorithms is illustrated on a benchmark problem. The numerical results show the robustness of our two-level adaptive BDDC algorithms with respect to the wave number, the number of subdomains and the mesh size, and illustrate that our multi-level adaptive BDDC algorithm can reduce the scale of the coarse problem and can be used to solve large wave number problems efficiently.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Algebraic multigrid block preconditioning for multi-group radiation diffusion equations

The paper focuses on developing and studying efficient block preconditioners based on classical algebraic multigrid for the large-scale sparse linear systems arising from the fully coupled and implicitly cell-centered finite volume discretization of multi-group radiation diffusion equations, whose coefficient matrices can be rearranged into the $(G+2)\times(G+2)$ block form, where $G$ is the number of energy groups. The preconditioning techniques are based on the monolithic classical algebraic multigrid method, physical-variable based coarsening two-level algorithm and two types of block Schur complement preconditioners. The classical algebraic multigrid is applied to solve the subsystems that arise in the last three block preconditioners. The coupling strength and diagonal dominance are further explored to improve performance. We use representative one-group and twenty-group linear systems from capsule implosion simulations to test the robustness, efficiency, strong and weak parallel scaling properties of the proposed methods. Numerical results demonstrate that block preconditioners lead to mesh- and problem-independent convergence, and scale well both algorithmically and in parallel.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Local Averaging Type a Posteriori Error Estimates for the Nonlinear Steady-state Poisson-Nernst-Planck Equations

The a posteriori error estimates are studied for a class of nonlinear stead-state Poisson-Nernst-Planck equations, which are a coupled system consisting of the Nernst-Planck equation and the Poisson equation. Both the global upper bounds and the local lower bounds of the error estimators are obtained by using a local averaging operator. Numerical experiments are given to confirm the reliability and efficiency of the error estimators.

0 citations0 reviewsNo code linkedOpen access