Source author record

Shuonan Wu

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

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

7 published item(s)

preprint2022arXiv

A $C^{0}$ finite element approximation of planar oblique derivative problems in non-divergence form

This paper proposes a $C^{0}$ (non-Lagrange) primal finite element approximation of the linear elliptic equations in non-divergence form with oblique boundary conditions in planar, curved domains. As an extension of [Calcolo, 58 (2022), No. 9], the Miranda-Talenti estimate for oblique boundary conditions at a discrete level is established by enhancing the regularity on the vertices. Consequently, the coercivity constant for the proposed scheme is exactly the same as that from PDE theory. The quasi-optimal order error estimates are established by carefully studying the approximation property of the finite element spaces. Numerical experiments are provided to verify the convergence theory and to demonstrate the accuracy and efficiency of the proposed methods.

preprint2022arXiv

A monotone discretization for integral fractional Laplacian on bounded Lipschitz domains: Pointwise error estimates under Hölder regularity

We propose a monotone discretization for the integral fractional Laplace equation on bounded Lipschitz domains with the homogeneous Dirichlet boundary condition. The method is inspired by a quadrature-based finite difference method of Huang and Oberman, but is defined on unstructured grids in arbitrary dimensions with a more flexible domain for approximating singular integral. The scale of the singular integral domain not only depends on the local grid size, but also on the distance to the boundary, since the Hölder coefficient of the solution deteriorates as it approaches the boundary. By using a discrete barrier function that also reflects the distance to the boundary, we show optimal pointwise convergence rates in terms of the Hölder regularity of the data on both quasi-uniform and graded grids. Several numerical examples are provided to illustrate the sharpness of the theoretical results.

preprint2022arXiv

Discontinuous Galerkin methods for magnetic advection-diffusion problems

We devise and analyze a class of the primal discontinuous Galerkin methods for the magnetic advection-diffusion problems based on the weighted-residual approach. In addition to the upwind stabilization, we find a new mechanism under the vector case that provides more flexibility in constructing the schemes. For the more general Friedrichs system, we show the stability and optimal error estimate, which boil down to two core ingredients -- the weight function and the special projection -- that contain information of advection. Numerical experiments are provided to verify the theoretical results.

preprint2020arXiv

$C^0$ finite element approximations of linear elliptic equations in non-divergence form and Hamilton-Jacobi-Bellman equations with Cordes coefficients

This paper is concerned with $C^0$ (non-Lagrange) finite element approximations of the linear elliptic equations in non-divergence form and the Hamilton-Jacobi-Bellman (HJB) equations with Cordes coefficients. Motivated by the Miranda-Talenti estimate, a discrete analog is proved once the finite element space is $C^0$ on the $(n-1)$-dimensional subsimplex (face) and $C^1$ on $(n-2)$-dimensional subsimplex. The main novelty of the non-standard finite element methods is to introduce an interior penalty term to argument the PDE-induced variational form of the linear elliptic equations in non-divergence form or the HJB equations. As a distinctive feature of the proposed methods, no penalization or stabilization parameter is involved in the variational forms. As a consequence, the coercivity constant (resp. monotonicity constant) for the linear elliptic equations in non-divergence form (resp. the HJB equations) at discrete level is exactly the same as that from PDE theory. The quasi-optimal order error estimates as well as the convergence of the semismooth Newton method are established. Numerical experiments are provided to validate the convergence theory and to illustrate the accuracy and computational efficiency of the proposed methods.

preprint2020arXiv

On the unisolvence for the quasi-polynomial spaces of differential forms

We consider quasi-polynomial spaces of differential forms defined as weighted (with a positive weight) spaces of differential forms with polynomial coefficients. We show that the unisolvent set of functionals for such spaces on a simplex in any spatial dimension is the same as the set of such functionals used for the polynomial spaces. The analysis in the quasi-polynomial spaces, however, is not standard and requires a novel approach. We are able to prove our results without the use of Stokes' Theorem, which is the standard tool in showing the unisolvence of functionals in polynomial spaces of differential forms. These new results provide tools for studying exponentially-fitted discretizations stable for general convection-diffusion problems in Hilbert differential complexes.

preprint2019arXiv

A mixed discontinuous Galerkin method with symmetric stress for Brinkman problem based on the velocity-pseudostress formulation

The Brinkman equations can be regarded as a combination of the Stokes and Darcy equations which model transitions between the fast flow in channels (governed by Stokes equations) and the slow flow in porous media (governed by Darcy's law). The numerical challenge for this model is the designing of a numerical scheme which is stable for both the Stokes-dominated (high permeability) and the Darcy-dominated (low permeability) equations. In this paper, we solve the Brinkman model in $n$ dimensions ($n = 2, 3$) by using the mixed discontinuous Galerkin (MDG) method, which meets this challenge. This MDG method is based on the pseudostress-velocity formulation and uses a discontinuous piecewise polynomial pair $\underline{\bm{\mathcal{P}}}_{k+1}^{\mathbb{S}}$-$\bm{\mathcal{P}}_k$ $(k \geq 0)$, where the stress field is symmetric. The main unknowns are the pseudostress and the velocity, whereas the pressure is easily recovered through a simple postprocessing. A key step in the analysis is to establish the parameter-robust inf-sup stability through specific parameter-dependent norms at both continuous and discrete levels. Therefore, the stability results presented here are uniform with respect to the permeability. Thanks to the parameter-robust stability analysis, we obtain optimal error estimates for the stress in broken $\underline{\bm{H}}(\bm{\rm div})$-norm and velocity in $\bm{L}^2$-norm. Furthermore, the optimal $\underline{\bm{L}}^2$ error estimate for pseudostress is derived under certain conditions. Finally, numerical experiments are provided to support the theoretical results and to show the robustness, accuracy, and flexibility of the MDG method.

preprint2016arXiv

A Second Order Time Homogenized Model for Sediment Transport

A multi-scale method for the hyperbolic systems governing sediment transport in subcritical case is developed. The scale separation of this problem is due to the fact that the sediment transport is much slower than flow velocity. We first derive a zeroth order homogenized model, and then propose a first order correction. It is revealed that the first order correction for hyperbolic systems has to be applied on the characteristic speed of slow variables in one dimensional case. In two dimensional case, besides the characteristic speed, the source term is also corrected. We develop a second order numerical scheme following the framework of heterogeneous multi-scale method. The numerical results in both one and two dimensional cases demonstrate the effectiveness and efficiency of our method.