BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Haijun Wu

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

ResearcherUnclaimed source record
math.NANumerical Analysisphysics.optics

Catalog footprint

What is connected

13works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

13 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Preasymptotic error analysis of higher order FEM and CIP-FEM for Helmholtz equation with high wave number

A preasymptotic error analysis of the finite element method (FEM) and some continuous interior penalty finite element method (CIP-FEM) for Helmholtz equation in two and three dimensions is proposed. $H^1$- and $L^2$- error estimates with explicit dependence on the wave number $k$ are derived. In particular, it is shown that if $k^{2p+1}h^{2p}$ is sufficiently small, then the pollution errors of both methods in $H^1$-norm are bounded by $O(k^{2p+1}h^{2p})$, which coincides with the phase error of the FEM obtained by existent dispersion analyses on Cartesian grids, where $h$ is the mesh size, $p$ is the order of the approximation space and is fixed. The CIP-FEM extends the classical one by adding more penalty terms on jumps of higher (up to $p$-th order) normal derivatives in order to reduce efficiently the pollution errors of higher order methods. Numerical tests are provided to verify the theoretical findings and to illustrate great capability of the CIP-FEM in reducing the pollution effect.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Multilevel Preconditioner with Stable Coarse Grid Corrections for the Helmholtz Equation

In this paper we consider a class of robust multilevel precontioners for the Helmholtz equation with high wave number. The key idea in this work is to use the continuous interior penalty finite element methods (CIP-FEM) studied in \cite{Wu12,Wu12-hp} to construct the stable coarse grid correction problems. The multilevel methods, based on GMRES smoothing on coarse grids, are then served as a preconditioner in the outer GMRES iteration. In the one dimensional case, convergence property of the modified multilevel methods is analyzed by the local Fourier analysis. From our numerical results, we find that the proposed methods are efficient for a reasonable range of frequencies. The performance of the algorithms depends relatively mildly on wave number. In particular, only one GMRES smoothing step may guarantee the optimal convergence of our multilevel algorithm, which remedies the shortcoming of the multilevel algorithm in \cite{EEO01}.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A combined finite element and multiscale finite element method for the multiscale elliptic problems

The oversampling multiscale finite element method (MsFEM) is one of the most popular methods for simulating composite materials and flows in porous media which may have many scales. But the method may be inapplicable or inefficient in some portions of the computational domain, e.g., near the domain boundary or near long narrow channels inside the domain due to the lack of permeability information outside of the domain or the fact that the high-conductivity features cannot be localized within a coarse-grid block. In this paper we develop a combined finite element and multiscale finite element method (FE-MsFEM), which deals with such portions by using the standard finite element method on a fine mesh and the other portions by the oversampling MsFEM. The transmission conditions across the FE-MSFE interface is treated by the penalty technique. A rigorous convergence analysis for this special FE-MsFEM is given under the assumption that the diffusion coefficient is periodic. Numerical experiments are carried out for the elliptic equations with periodic and random highly oscillating coefficients, as well as multiscale problems with high contrast channels, to demonstrate the accuracy and efficiency of the proposed method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An absolutely stable discontinuous Galerkin method for the indefinite time-harmonic Maxwell equations with large wave number

This paper develops and analyzes an interior penalty discontinuous Galerkin (IPDG) method using piecewise linear polynomials for the indefinite time harmonic Maxwell equations with the impedance boundary condition in the three dimensional space. The main novelties of the proposed IPDG method include the following: first, the method penalizes not only the jumps of the tangential component of the electric field across the element faces but also the jumps of the tangential component of its vorticity field; second, the penalty parameters are taken as complex numbers of negative imaginary parts. For the differential problem, we prove that the sesquilinear form associated with the Maxwell problem satisfies a generalized weak stability (i.e., inf-sup condition) for star-shaped domains.Such a generalized weak stability readily infers wave-number explicit a priori estimates for the solution of the Maxwell problem, which plays an important role in the error analysis for the IPDG method. For the proposed IPDG method, we show that the discrete sesquilinear form satisfies a coercivity for all positive mesh size $h$ and wave number $k$ and for general domains including non-star-shaped ones. In turn, the coercivity easily yields the well-posedness and stability estimates (i.e., a priori estimates) for the discrete problem without imposing any mesh constraint. Based on these discrete stability estimates, by adapting a nonstandard error estimate technique of Fung and Wu (2009), we derive both the energy-norm and the $L^2$-norm error estimates for the IPDG method in all mesh parameter regimes including pre-asymptotic regime (i.e., $k^2 h\gtrsim 1$). Numerical experiments are also presented to gauge the theoretical results and to numerically examine the pollution effect (with respect to $k$) in the error bounds.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Continuous Interior Penalty Finite Element Method for Helmholtz Equation with High Wave Number: One Dimensional Analysis

This paper addresses the properties of Continuous Interior Penalty (CIP) finite element solutions for the Helmholtz equation. The $h$-version of the CIP finite element method with piecewise linear approximation is applied to a one-dimensional model problem. We first show discrete well posedness and convergence results, using the imaginary part of the stabilization operator, for the complex Helmholtz equation. Then we consider a method with real valued penalty parameter and prove an error estimate of the discrete solution in the $H^1$-norm, as the sum of best approximation plus a pollution term that is the order of the phase difference. It is proved that the pollution can be eliminated by selecting the penalty parameter appropriately. As a result of this analysis, thorough and rigorous understanding of the error behavior throughout the range of convergence is gained. Numerical results are presented that show sharpness of the error estimates and highlight some phenomena of the discrete solution behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Pre-asymptotic Error Analysis of CIP-FEM and FEM for Helmholtz Equation with High Wave Number. Part II: $hp$ version

In this paper, which is part II in a series of two, the pre-asymptotic error analysis of the continuous interior penalty finite element method (CIP-FEM) and the FEM for the Helmholtz equation in two and three dimensions is continued. While part I contained results on the linear CIP-FEM and FEM, the present part deals with approximation spaces of order $p \ge 1$. By using a modified duality argument, pre-asymptotic error estimates are derived for both methods under the condition of $\frac{kh}{p}\le C_0\big(\frac{p}{k}\big)^{\frac{1}{p+1}}$, where $k$ is the wave number, $h$ is the mesh size, and $C_0$ is a constant independent of $k, h, p$, and the penalty parameters. It is shown that the pollution errors of both methods in $H^1$-norm are $O(k^{2p+1}h^{2p})$ if $p=O(1)$ and are $O\Big(\frac{k}{p^2}\big(\frac{kh}{σp}\big)^{2p}\Big)$ if the exact solution $u\in H^2(\Om)$ which coincide with existent dispersion analyses for the FEM on Cartesian grids. Here $\si$ is a constant independent of $k, h, p$, and the penalty parameters. Moreover, it is proved that the CIP-FEM is stable for any $k, h, p>0$ and penalty parameters with positive imaginary parts. Besides the advantage of the absolute stability of the CIP-FEM compared to the FEM, the penalty parameters may be tuned to reduce the pollution effects.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Continuous Interior Penalty Finite Element Methods for the Helmholtz Equation with Large Wave Number

This paper develops and analyzes some continuous interior penalty finite element methods (CIP-FEMs) using piecewise linear polynomials for the Helmholtz equation with the first order absorbing boundary condition in two and three dimensions. The novelty of the proposed methods is to use complex penalty parameters with positive imaginary parts. It is proved that, if the penalty parameter is a pure imaginary number $ı\ga$ with $0<\ga\le C$, then the proposed CIP-FEM is stable (hence well-posed) without any mesh constraint. Moreover the method satisfies the error estimates $C_1kh+C_2k^3h^2$ in the $H^1$-norm when $k^3h^2\le C_0$ and $C_1kh+\frac{C_2}{\ga}$ when $k^3h^2> C_0$ and $kh$ is bounded, where $k$ is the wave number, $h$ is the mesh size, and the $C$'s are positive constants independent of $k$, $h$, and $\ga$. Optimal order $L^2$ error estimates are also derived. The analysis is also applied if the penalty parameter is a complex number with positive imaginary part. By taking $\ga\to 0+$, the above estimates are extended to the linear finite element method under the condition $k^3h^2\le C_0$. Numerical results are provided to verify the theoretical findings. It is shown that the penalty parameters may be tuned to greatly reduce the pollution errors.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An unfitted $hp$-interface penalty finite element method for elliptic interface problems

An $hp$ version of interface penalty finite element method ($hp$-IPFEM) is proposed for elliptic interface problems in two and three dimensions on unfitted meshes. Error estimates in broken $H^1$ norm, which are optimal with respect to $h$ and suboptimal with respect to $p$ by half an order of $p$, are derived. Both symmetric and non-symmetric IPFEM are considered. Error estimates in $L^2$ norm are proved by the duality argument.

0 citations0 reviewsNo code linkedOpen access