Researcher profile

Haijun Wu

Haijun Wu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

preprint2022arXiv

A combined multiscale finite element method based on the LOD technique for the multiscale elliptic problems with singularities

In this paper, we construct a combined multiscale finite element method (MsFEM) using the Local Orthogonal Decomposition (LOD) technique to solve the multiscale problems which may have singularities in some special portions of the computational domain. For example, in the simulation of steady flow transporting through highly heterogeneous porous media driven by extraction wells, the singularities lie in the near-well regions. The basic idea of the combined method is to utilize the traditional finite element method (FEM) directly on a fine mesh of the problematic part of the domain and using the LOD-based MsFEM on a coarse mesh of the other part. The key point is how to define local correctors for the basis functions of the elements near the coarse and fine mesh interface, which require meticulous treatment. The proposed method takes advantages of the traditional FEM and the LOD-based MsFEM, which uses much less DOFs than the standard FEM and may be more accurate than the LOD-based MsFEM for problems with singularities. The error analysis is carried out for highly varying coefficients, without any assumptions on scale separation or periodicity. {Numerical examples with periodic and random highly oscillating coefficients}, as well as the multiscale problems on the L-shaped domain, and multiscale problems with high-contrast channels or well-singularities are presented to demonstrate the efficiency and accuracy of the proposed method.

preprint2022arXiv

Dispersion Analysis of CIP-FEM for Helmholtz Equation

When solving the Helmholtz equation numerically, the accuracy of numerical solution deteriorates as the wave number $k$ increases, known as `pollution effect' which is directly related to the phase difference between the exact and numerical solutions, caused by the numerical dispersion. In this paper, we propose a dispersion analysis for the continuous interior penalty finite element method (CIP-FEM) and derive an explicit formula of the penalty parameter for the $p^{\rm th}$ order CIP-FEM on tensor product (Cartesian) meshes, with which the phase difference is reduced from $\mathcal{O}\big(k(kh)^{2p}\big)$ to $\mathcal{O}\big(k(kh)^{2p+2}\big)$. Extensive numerical tests show that the pollution error of the CIP-FE solution is also reduced by two orders in $kh$ with the same penalty parameter.

preprint2022arXiv

Finite Element Method for a Nonlinear PML Helmholtz Equation with High Wave Number

A nonlinear Helmholtz equation (NLH) with high wave number and Sommerfeld radiation condition is approximated by the perfectly matched layer (PML) technique and then discretized by the linear finite element method (FEM). Wave-number-explicit stability and regularity estimates and the exponential convergence are proved for the nonlinear truncated PML problem. Preasymptotic error estimates are obtained for the FEM, where the logarithmic factors in h required by the previous results for the NLH with impedance boundary condition are removed in the case of two dimensions. Moreover, local quadratic convergences of the Newton's methods are derived for both the NLH with PML and its FEM. Numerical examples are presented to verify the accuracy of the FEM, which demonstrate that the pollution errors may be greatly reduced by applying the interior penalty technique with proper penalty parameters to the FEM. The nonlinear phenomenon of optical bistability can be successfully simulated.

preprint2022arXiv

Topological Transformation and Free-Space Transport of Photonic Hopfions

Structured light fields embody strong spatial variations of polarisation, phase and amplitude. Understanding, characterization and exploitation of such fields can be achieved through their topological properties. Three-dimensional (3D) topological solitons, such as hopfions, are 3D localized continuous field configurations with nontrivial particle-like structures, that exhibit a host of important topologically protected properties. Here, we propose and demonstrate photonic counterparts of hopfions with exact characteristics of Hopf fibration, Hopf index, and Hopf mapping from real-space vector beams to homotopic hyperspheres representing polarisation states. We experimentally generate photonic hopfions with on-demand high-order Hopf indices and independently controlled topological textures, including Néel-, Bloch-, and anti-skyrmionic types. We also demonstrate a robust free-space transport of photonic hopfions, thus, showing potential of hopfions for developing optical topological informatics and communications.

preprint2020arXiv

Hybridizable Discontinuous Galerkin Methods for Helmholtz Equation with High Wave Number. Part I: Linear case

This paper addresses several aspects of the linear Hybridizable Discontinuous Galerkin Method (HDG) for the Helmholtz equation with impedance boundary condition at high frequency. First, error estimates with explicit dependence on the wave number $k$ for the HDG approximations to the exact solution $u$ and its negative gradient $\mathbf{q}=-\nabla u$ are derived. It is shown that $k\Vert u - u_h \Vert_{L^2(Ω)} + \Vert \mathbf{q} -\mathbf{q}_h \Vert_{L^2(Ω)} = O(k^2h^2+k^4h^3)$ under the conditions that $k^3h^2$ is sufficiently small and that the penalty parameter $τ\eqsim k$, where $h$ is the mesh size. Note that the convergence order in $\mathbf{q}_h$ is full and the pollution error is $O(k^4h^3)$, which improve the existent results. Secondly, by using a standard postprocessing procedure from the HDG method for elliptic problems, a piecewise quadratic function $u_h^*$ is obtained so that $k\Vert u-u_h^*\Vert_{L^2(Ω)}=O(k^3h^3+k^4h^3)$. Note that the postprocessing procedure improves only the interpolation error (from $O(k^2h^2)$ to $O(k^3h^3)$) but leaves the pollution error $O(k^4h^3)$ unchanged. Thirdly, dispersion analyses and extensive numerical tests show that the pollution effect can be eliminated completely in 1D case and reduced greatly in 2D case by selecting appropriate penalty parameters. The preasymptotic error analysis of the higher order HDG method for the Helmholtz equation with high wave number is studied in Part II.