Researcher profile

Cuiyu He

Cuiyu He contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
2topics
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

5 published item(s)

preprint2026arXiv

A Unified Equilibrated Flux Recovery Framework with Robust A Posteriori Error Estimation

We introduce the Equilibrated Averaging Residual Method (EARM), a unified equilibrated flux-recovery framework for elliptic interface problems that applies to a broad class of finite element discretizations. The method is applicable in both two and three dimensions and for arbitrary polynomial orders, and it enables the construction of computationally efficient recovered fluxes. We develop EARM for both discontinuous Galerkin (DG) and conforming finite element discretizations. For DG methods, EARM can be applied directly and yields an explicit recovered flux that coincides with state-of-the-art conservative flux reconstructions. For conforming discretizations, we further propose the Orthogonal Null-space--Eliminated EARM (ON-EARM), which ensures uniqueness by restricting the correction flux to the orthogonal complement of the divergence-free null space. We prove local conservation and establish a robust a~posteriori error estimator for the recovered flux in two dimensions, with robustness measured with respect to jumps in the diffusion coefficient. Numerical results in two and three dimensions confirm the theoretical findings.

preprint2021arXiv

Towards fast weak adversarial training to solve high dimensional parabolic partial differential equations using XNODE-WAN

Due to the curse of dimensionality, solving high dimensional parabolic partial differential equations (PDEs) has been a challenging problem for decades. Recently, a weak adversarial network (WAN) proposed in (Y.Zang et al., 2020) offered a flexible and computationally efficient approach to tackle this problem defined on arbitrary domains by leveraging the weak solution. WAN reformulates the PDE problem as a generative adversarial network, where the weak solution (primal network) and the test function (adversarial network) are parameterized by the multi-layer deep neural networks (DNNs). However, it is not yet clear whether DNNs are the most effective model for the parabolic PDE solutions as they do not take into account the fundamentally different roles played by time and spatial variables in the solution. To reinforce the difference, we design a novel so-called XNODE model for the primal network, which is built on the neural ODE (NODE) model with additional spatial dependency to incorporate the a priori information of the PDEs and serve as a universal and effective approximation to the solution. The proposed hybrid method (XNODE-WAN), by integrating the XNODE model within the WAN framework, leads to significant improvement in the performance and efficiency of training. Numerical results show that our method can reduce the training time to a fraction of that of the WAN model.

preprint2020arXiv

A Mesh-free Method Using Piecewise Deep Neural Network for Elliptic Interface Problems

In this paper, we propose a novel mesh-free numerical method for solving the elliptic interface problems based on deep learning. We approximate the solution by the neural networks and, since the solution may change dramatically across the interface, we employ different neural networks in different sub-domains. By reformulating the interface problem as a least-squares problem, we discretize the objective function using mean squared error via sampling and solve the proposed deep least-squares method by standard training algorithms such as stochastic gradient descent. The discretized objective function utilizes only the point-wise information on the sampling points and thus no underlying mesh is required. Doing this circumvents the challenging meshing procedure as well as the numerical integration on the complex interface. To improve the computational efficiency for more challenging problems, we further design an adaptive sampling strategy based on the residual of the least-squares function and propose an adaptive algorithm. Finally, we present several numerical experiments in both 2D and 3D to show the flexibility, effectiveness, and accuracy of the proposed deep least-square method for solving interface problems.

preprint2020arXiv

Comparison of Shape Derivatives using CutFEM for Ill-posed Bernoulli Free Boundary Problem

In this paper we discuss a level set approach for the identification of an unknown boundary in a computational domain. The problem takes the form of a Bernoulli problem where only the Dirichlet datum is known on the boundary that is to be identified, but additional information on the Neumann condition is available on the known part of the boundary. The approach uses a classical constrained optimization problem, where a cost functional is minimized with respect to the unknown boundary, the position of which is defined implicitly by a level set function. To solve the optimization problem a steepest descent algorithm using shape derivatives is applied. In each iteration the cut finite element method is used to obtain high accuracy approximations of the pde-model constraint for a given level set configuration without re-meshing. We consider three different shape derivatives. First the classical one, derived using the continuous optimization problem (optimize then discretize). Then the functional is first discretized using the CutFEM method and the shape derivative is evaluated on the finite element functional (discretize then optimize). Finally we consider a third approach, also using a discretized functional. In this case we do not perturb the domain, but consider a so-called boundary value correction method, where a small correction to the boundary position may be included in the weak boundary condition. Using this correction the shape derivative may be obtained by perturbing a distance parameter in the discrete variational formulation. The theoretical discussion is illustrated with a series of numerical examples showing that all three approaches produce similar result on the proposed Bernoulli problem.

preprint2020arXiv

Generalized Prager-Synge Inequality and Equilibrated Error Estimators for Discontinuous Elements

The well-known Prager-Synge identity is valid in $H^1(Ω)$ and serves as a foundation for developing equilibrated a posteriori error estimators for continuous elements. In this paper, we introduce a new inequality, that may be regarded as a generalization of the Prager-Synge identity, to be valid for piecewise $H^1(Ω)$ functions for diffusion problems. The inequality is proved to be identity in two dimensions. For nonconforming finite element approximation of arbitrary odd order, we propose a fully explicit approach that recovers an equilibrated flux in $H(div; Ω)$ through a local element-wise scheme and that recovers a gradient in $H(curl;Ω)$ through a simple averaging technique over edges. The resulting error estimator is then proved to be globally reliable and locally efficient. Moreover, the reliability and efficiency constants are independent of the jump of the diffusion coefficient regardless of its distribution.