Source author record

Jun Zou

Jun Zou 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

37works
12topics
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

37 published item(s)

preprint2022arXiv

A family of immersed finite element spaces and applications to three dimensional $\mathbf{H}(\text{curl})$ interface problems

Maxwell interface problems are of great importance in many electromagnetic applications. Unfitted mesh methods are especially attractive in 3D computation as they can circumvent generating complex 3D interface-fitted meshes. However, many unfitted mesh methods rely on non-conforming approximation spaces, which may cause a loss of accuracy for solving Maxwell equations, and the widely-used penalty techniques in the literature may not help in recovering the optimal convergence. In this article, we provide a remedy by developing Nédélec-type immersed finite element spaces with a Petrov-Galerkin scheme that is able to produce optimal-convergent solutions. To establish a systematic framework, we analyze all the $H^1$, $\mathbf{H}(\text{curl})$ and $\mathbf{H}(\text{div})$ IFE spaces and form a discrete de Rham complex. Based on these fundamental results, we further develop a fast solver using a modified Hiptmair-Xu preconditioner which works for both the GMRES and CG methods.

preprint2022arXiv

Asymptotic Analysis for a Nonlinear Reaction-Diffusion System Modeling an Infectious Disease

In this paper we study a nonlinear reaction-diffusion system which models an infectious disease caused by bacteria such as those for cholera. One of the significant features in this model is that a certain portion of the recovered human hosts may lose a lifetime immunity and could be infected again. Another important feature in the model is that the mobility for each species is allowed to be dependent upon both the location and time. With the whole population assumed to be susceptible with the bacteria, the model is a strongly coupled nonlinear reaction-diffusion system. We prove that the nonlinear system has a unique solution globally in any space dimension under some natural conditions on the model parameters and the given data. Moreover, the long-time behavior and stability analysis for the solutions are carried out rigorously. In particular, we characterize the precise conditions on variable parameters about the stability or instability of all steady-state solutions. These new results provide the answers to several open questions raised in the literature.

preprint2022arXiv

Elastodynamical resonances and cloaking of negative material structures beyond quasistatic approximation

Given the flexibility of choosing negative elastic parameters, we construct material structures that can induce two resonance phenomena, referred to as the elastodynamical resonances. They mimic the emerging plasmon/polariton resonance and anomalous localized resonance in optics for subwavelength particles. However, we study the peculiar resonance phenomena for linear elasticity beyond the subwavelength regime. It is shown that the resonance behaviours possess distinct characters, with some similar to the subwavelength resonances, but some sharply different due to the frequency effect. It is particularly noted that we construct a core-shell material structure that can induce anomalous localized resonance as well as cloaking phenomena beyond the quasi-static limit. The study is boiled down to analyzing the so-called elastic Neumann-Poincáre (N-P) operator in the frequency regime. We provide an in-depth analysis of the spectral properties of the N-P operator beyond the quasi-static approximation,and these results are of independent interest to the spectral theory of layer potential operators.

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

Least-Squares Method for Inverse Medium Problems

We present a two-stage least-squares method to inverse medium problems of reconstructing multiple unknown coefficients simultaneously from noisy data. A direct sampling method is applied to detect the location of the inhomogeneity in the first stage, while a total least-squares method with mixed regularization is used to recover the medium profile in the second stage. The total least-squares method is designed to minimize the residual of the model equation and the data fitting, along with an appropriate regularization, in an attempt to significantly improve the accuracy of the approximation obtained from the first stage. We shall also present an analysis on the well-posedness and convergence of this algorithm. Numerical experiments are carried out to verify the accuracies and robustness of this novel two-stage least-squares algorithm, with great tolerance of noise.

preprint2022arXiv

Mittag--Leffler stability of numerical solutions to time fractional ODEs

The asymptotic stable region and long-time decay rate of solutions to linear homogeneous Caputo time fractional ordinary differential equations (F-ODEs) are known to be completely determined by the eigenvalues of the coefficient matrix. Very different from the exponential decay of solutions to classical ODEs, solutions of F-ODEs decay only polynomially, leading to the so-called Mittag-Leffler stability, which was already extended to semi-linear F-ODEs with small perturbations. This work is mainly devoted to the qualitative analysis of the long-time behavior of numerical solutions. By applying the singularity analysis of generating functions developed by Flajolet and Odlyzko (SIAM J. Disc. Math. 3 (1990), 216-240), we are able to prove that both $\mathcal{L}$1 scheme and strong $A$-stable fractional linear multistep methods (F-LMMs) can preserve the numerical Mittag-Leffler stability for linear homogeneous F-ODEs exactly as in the continuous case. Through an improved estimate of the discrete fractional resolvent operator, we show that strong $A$-stable F-LMMs are also Mittag-Leffler stable for semi-linear F-ODEs under small perturbations. For the numerical schemes based on $α$-difference approximation to Caputo derivative, we establish the Mittag-Leffler stability for semi-linear problems by making use of properties of the Poisson transformation and the decay rate of the continuous fractional resolvent operator. Numerical experiments are presented for several typical time fractional evolutional equations, including time fractional sub-diffusion equations, fractional linear system and semi-linear F-ODEs. All the numerical results exhibit the typical long-time polynomial decay rate, which is fully consistent with our theoretical predictions.

preprint2022arXiv

Performance Analysis of Wireless Network Aided by Discrete-Phase-Shifter IRS

Discrete phase shifters of intelligent reflecting surface (IRS) generates phase quantization error (QE) and degrades the receive performance at the receiver. To make an analysis of the performance loss caused by IRS with phase QE, based on the law of large numbers, the closed-form expressions of signal-to-noise ratio (SNR) performance loss (PL), achievable rate (AR), and bit error rate (BER) are successively derived under line-of-sight (LoS) channels and Rayleigh channels. Moreover, based on the Taylor series expansion, the approximate simple closed form of PL of IRS with approximate QE is also given. The simulation results show that the performance losses of SNR and AR decrease as the number of quantization bits increase, while they gradually increase with the number of IRS phase shifter elements increase. Regardless of LoS channels or Rayleigh channels, when the number of quantization bits is larger than or equal to 3, the performance losses of SNR and AR are less than 0.23dB and 0.08bits/s/Hz, respectively, and the BER performance degradation is trivial. In particular, the performance loss difference between IRS with QE and IRS with approximate QE is negligible when the number of quantization bits is not less than 2.

preprint2021arXiv

An Analysis of Stochastic Variance Reduced Gradient for Linear Inverse Problems

Stochastic variance reduced gradient (SVRG) is a popular variance reduction technique for accelerating stochastic gradient descent (SGD). We provide a first analysis of the method for solving a class of linear inverse problems in the lens of the classical regularization theory. We prove that for a suitable constant step size schedule, the method can achieve an optimal convergence rate in terms of the noise level (under suitable regularity condition) and the variance of the SVRG iterate error is smaller than that by SGD. These theoretical findings are corroborated by a set of numerical experiments.

preprint2021arXiv

Minnaert resonances for bubbles in soft elastic materials

Minnaert resonance is a well known physical phenomenon and it has many important applications, in particular in the effective realisation of acoustic metamaterials using bubbly mediums in recent years. In this paper, motived by the Minnaert resonance, we consider the low-frequency resonance for acoustic bubbles embedded in soft elastic materials. This is a hybrid physical process that couples the acoustic and elastic wave propagations. By delicately and subtly balancing the acoustic and elastic parameters as well as the geometry of the bubble, we show that Minnaert resonance can (approximately) occur for rather general constructions. Our study poses a great potential for the effective realisation of negative elastic materials by using bubbly elastic mediums.

preprint2020arXiv

A Direct Sampling Method for Simultaneously Recovering Inhomogeneous Inclusions of Different Nature

In this work, we investigate a class of elliptic inverse problems and aim to simultaneously recover multiple inhomogeneous inclusions arising from two different physical parameters, using very limited boundary Cauchy data collected only at one or two measurement events. We propose a new fast, stable and highly parallelable direct sampling method (DSM) for the simultaneous reconstruction process. Two groups of probing and index functions are constructed, and their desired properties are analyzed. In order to identify and decouple the multiple inhomogeneous inclusions of different physical nature, we introduce a new concept of mutually almost orthogonality property that generalizes the important concept of almost orthogonality property in classical DSMs for inhomogeneous inclusions of same physical nature. With the help of this new concept, we develop a reliable strategy to distinguish two different types of inhomogeneous inclusions with noisy data collected at one or two measurement events. We further improve the decoupling effect by choosing an appropriate boundary influx. Numerical experiments are presented to illustrate the robustness and efficiency of the proposed method.

preprint2020arXiv

Receive Antenna Selection for Secure Pre-coding Aided Spatial Modulation

In this paper, we make an investigation of receive antenna selection (RAS) strategies in the secure pre-coding aided spatial modulation (PSM) system with the aid of artificial noise. Due to a lack of the closed-form expression for secrecy rate (SR) in secure PSM systems, it is hard to optimize the RAS. To address this issue, the cut-off rate is used as an approximation of the SR. Further, two low-complexity RAS schemes for maximizing SR, called Max-SR-L and Max-SR-H, are derived in the low and high signal-to-noise ratio (SNR) regions, respectively. Due to the fact that the former works well in the low SNR region but becomes worse in the medium and high SNR regions while the latter also has the similar problem, a novel RAS strategy Max-SR-A is proposed to cover all SNR regions. Simulation results show that the proposed Max-SR-H and Max-SR-L schemes approach the optimal SR performances of the exhaustive search (ES) in the high and low SNR regions, respectively. In particular, the SR performance of the proposed Max-SR-A is close to that of the optimal ES and better than that of the random method in almost all SNR regions.

preprint2020arXiv

Sequential and Incremental Precoder Design for Joint Transmission Network MIMO Systems with Imperfect Backhaul

In this paper, we propose a sequential and incremental precoder design for downlink joint transmission (JT) network MIMO systems with imperfect backhaul links. The objective of our design is to minimize the maximum of the sub-stream mean square errors (MSE), which dominates the average bit error rate (BER) performance of the system. In the proposed scheme,we first optimize the precoder at the serving base station (BS), and then sequentially optimize the precoders of non-serving BSs in the JT set according to the descending order of their probabilities of participating in JT. The BS-wise sequential optimization process can improve the system performance when some BSs have to temporarily quit the JT operations because of poor instant backhaul conditions. Besides, the precoder of an additional BS is derived in an incremental way, i.e., the sequentially optimized precoders of previous BSs are fixed, thus the additional precoder plays an incremental part in the multi-BS JT operations. An iterative algorithm is designed to jointly optimize the sub-stream precoder and sub-stream power allocation for each additional BS in the proposed sequential and incremental optimization scheme. Simulations show that, under the practical backhaul link conditions, our scheme significantly outperforms the autonomous global precoding (AGP) scheme in terms of BER performance.

preprint2020arXiv

Super-resolution in recovering embedded electromagnetic sources in high contrast media

The purpose of this work is to provide a rigorous mathematical analysis of the expected super-resolution phenomenon in the time-reversal imaging of electromagnetic (EM) radiating sources embedded in a high contrast medium. It is known that the resolution limit is essentially determined by the sharpness of the imaginary part of the EM Green's tensor for the associated background. We first establish the close connection between the resolution and the material parameters and the resolvent of the electric integral operator, via the Lippmann-Schwinger representation formula. We then present an insightful characterization of the spectral structure of the integral operator for a general bounded domain and derive the pole-pencil decomposition of its resolvent in the high contrast regime. For the special case of a spherical domain, we provide some quantitative asymptotic behavior of the eigenvalues and eigenfunctions. These mathematical findings shall enable us to provide a concise and rigorous illustration of the super-resolution in the EM source reconstruction in high contrast media. Some numerical examples are also presented to verify our main theoretical results.

preprint2020arXiv

Unique continuation from a generalized impedance edge-corner for Maxwell's system and applications to inverse problems

We consider the time-harmonic Maxwell system in a domain with a generalized impedance edge-corner, namely the presence of two generalized impedance planes that intersect at an edge. The impedance parameter can be $0, \infty$ or a finite non-identically vanishing variable function. We establish an accurate relationship between the vanishing order of the solutions to the Maxwell system and the dihedral angle of the edge-corner. In particular, if the angle is irrational, the vanishing order is infinity, i.e. strong unique continuation holds from the edge-corner. The establishment of those new quantitative results involve a highly intricate and subtle algebraic argument. The unique continuation study is strongly motivated by our study of a longstanding inverse electromagnetic scattering problem. As a significant application, we derive several novel unique identifiability results in determining a polyhedral obstacle as well as it surface impedance by a single far-field measurement. We also discuss another potential and interesting application of our result in the inverse scattering theory related to the information encoding.

preprint2019arXiv

Convergence Rates of Tikhonov Regularizations for Elliptic and Parabolic Inverse Radiativity Problems

We shall study in this paper the Lipschitz type stabilities and convergence rates of Tikhonov regularization for the recovery of the radiativities in elliptic and parabolic systems with Dirichlet boundary conditions. The Lipschitz type stability estimates are derived. Due to the difficulty of the verification of the existing source conditions or nonlinearity conditions for the considered inverse radiativity problems in high dimensional spaces, some new variational source conditions are proposed. The conditions are rigorously verified in general dimensional spaces under the Lipschitz type stability estimates and the reasonable convergence rates are achieved.

preprint2016arXiv

A Multilevel Sampling Method for Detecting Sources in a Stratified Ocean Waveguide

In the reconstruction process of sound waves in a 3D stratified waveguide, a key technique is to effectively reduce the huge computational demand. In this work, we propose an efficient and simple multilevel reconstruction method to help locate the accurate position of a point source in a stratified ocean. The proposed method can be viewed as a direct sampling method since no solutions of optimizations or linear systems are involved. The novel method exhibits several strengths: fast convergence, robustness against noise, advantages in computational complexity and applicability for a very small number of receivers.

preprint2016arXiv

Cellular Communications on License-Exempt Spectrum: A Tutorial

A traditional cellular system (e.g., LTE) operates only on the licensed spectrum. This tutorial explains the concept of cellular communications on both licensed and license-exempt spectrum under a unified architecture. The purpose to extend a cellular system into the bandwidth-rich license-exempt spectrum is to form a larger cellular network for all spectrum types. This would result in an ultimate mobile converged cellular network. This tutorial examines the benefits of this concept, the technical challenges, and provides a conceptual LTE-based design example that helps to show how a traditional cellular system like the LTE can adapt itself to a different spectrum type, conform to the regulatory requirements, and harmoniously co-exist with the incumbent systems such as Wi-Fi. In order to cope with the interference and regulation rules on license-exempt spectrum, a special medium access mechanism is introduced into the existing LTE transmission frame structure to exploit the full benefits of coordinated and managed cellular architecture.

preprint2016arXiv

Elastic-net regularization versus $\ell^1$-regularization for linear inverse problems with quasi-sparse solutions

We consider the ill-posed operator equation $Ax=y$ with an injective and bounded linear operator $A$ mapping between $\ell^2$ and a Hilbert space $Y$, possessing the unique solution \linebreak $x^†=\{x^†_k\}_{k=1}^\infty$. For the cases that sparsity $x^†\in \ell^0$ is expected but often slightly violated in practice, we investigate in comparison with the $\ell^1$-regularization the elastic-net regularization, where the penalty is a weighted superposition of the $\ell^1$-norm and the $\ell^2$-norm square, under the assumption that $x^†\in \ell^1$. There occur two positive parameters in this approach, the weight parameter $η$ and the regularization parameter as the multiplier of the whole penalty in the Tikhonov functional, whereas only one regularization parameter arises in $\ell^1$-regularization. Based on the variational inequality approach for the description of the solution smoothness with respect to the forward operator $A$ and exploiting the method of approximate source conditions, we present some results to estimate the rate of convergence for the elastic-net regularization. The occurring rate function contains the rate of the decay $x^†_k \to 0$ for $k \to \infty$ and the classical smoothness properties of $x^†$ as an element in $\ell^2$.

preprint2016arXiv

Preconditioners and Their Analyses for Edge Element Saddle-point Systems Arising from Time-harmonic Maxwell Equations

We shall propose and analyze some new preconditioners for the saddle-point systems arising from the edge element discretization of the time-harmonic Maxwell equations in three dimensions. We will first consider the saddle-point systems with vanishing wave number, for which we present an important relation between the solutions of the singular curl-curl system and the non-singular saddle-point system, then demonstrate that the saddle-point system can be efficiently solved by the Hiptmair-Xu solver. For the saddle-point systems with non-vanishing wave numbers, we will show that the PCG with a new preconditioner can apply for the non-singular system when wave numbers are small, while the methods like preconditioned MINRES may apply for some existing and new preconditioners when wave numbers are large. The spectral behaviors of the resulting preconditioned systems for the existing and new preconditioners are analyzed and compared, and numerical experiments are presented to demonstrate and compare the efficiencies of these preconditioners.

preprint2015arXiv

A parallel space-time domain decomposition method for unsteady source inversion problems

In this paper, we propose a parallel space-time domain decomposition method for solving an unsteady source identification problem governed by the linear convection-diffusion equation. Traditional approaches require to solve repeatedly a forward parabolic system, an adjoint system and a system with respect to the unknowns. The three systems have to be solved one after another. These sequential steps are not desirable for large scale parallel computing. A space-time restrictive additive Schwarz method is proposed for a fully implicit space-time coupled discretization scheme to recover the time-dependent pollutant source intensity functions. We show with numerical experiments that the scheme works well with noise in the observation data. More importantly it is demonstrated that the parallel space-time Schwarz preconditioner is scalable on a supercomputer with over $10^3$ processors, thus promising for large scale applications.

preprint2015arXiv

Analysis on Non-negative Factorizations and Applications

In this work we perform some mathematical analysis on non-negative matrix factorizations (NMF) and apply NMF to some imaging and inverse problems. We will propose a sparse low-rank approximation of big positive data and images in terms of tensor products of positive vectors, and investigate its effectiveness in terms of the number of tensor products to be used in the approximation. A new concept of multi-level analysis (MLA) framework is also suggested to extract major components in the matrix representing structures of different resolutions, but still preserving the positivity of the basis and sparsity of the approximation. We will also propose a semi-smooth Newton method based on primal-dual active sets for the non-negative factorization. Numerical results are given to demonstrate the effectiveness of the proposed method to capture features in images and structures of inverse problems under no a-priori assumption on the data structure, as well as to provide a sparse low-rank representation of the data.

preprint2015arXiv

Phased and phaseless domain reconstruction in inverse scattering problem via scattering coefficients

In this work we shall review the (phased) inverse scattering problem and then pursue the phaseless reconstruction from far-field data with the help of the concept of scattering coefficients. We perform sensitivity, resolution and stability analysis of both phased and phaseless problems and compare the degree of ill-posedness of the phased and phaseless reconstructions. The phaseless reconstruction is highly nonlinear and much more severely ill-posed. Algorithms are provided to solve both the phased and phaseless reconstructions in the linearized case. Stability is studied by estimating the condition number of the inversion process for both the phased and phaseless cases. An optimal strategy is suggested to attain the infimum of the condition numbers of the phaseless reconstruction, which may provide an important guidance for efficient phaseless measurements in practical applications. To the best of our knowledge, the stability analysis in terms of condition numbers are new for the phased and phaseless inverse scattering problems, and are very important to help us understand the degree of ill-posedness of these inverse problems. Numerical experiments are provided to illustrate the theoretical asymptotic behavior, as well as the effectiveness and robustness of the phaseless reconstruction algorithm.

preprint2015arXiv

Two-level space-time domain decomposition methods for unsteady inverse problems

As the number of processor cores on supercomputers becomes larger and larger, algorithms with high degree of parallelism attract more attention. In this work, we propose a novel space-time coupled algorithm for solving an inverse problem associated with the time-dependent convection-diffusion equation in three dimensions. We introduce a mixed finite element/finite difference method and a one-level and a two-level space-time parallel domain decomposition preconditioner for the Karush-Kuhn-Tucker (KKT) system induced from reformulating the inverse problem as an output least-squares optimization problem in the space-time domain. The new full space approach eliminates the sequential steps of the optimization outer loop and the inner forward and backward time marching processes, thus achieves high degree of parallelism. Numerical experiments validate that this approach is effective and robust for recovering unsteady moving sources. We report strong scalability results obtained on a supercomputer with more than 1,000 processors.

preprint2015arXiv

Well-posedness of a Pulsed Electric Field Model in Biological Media and its Finite Element Approximation

This work aims at providing a mathematical and numerical framework for the analysis on the effects of pulsed electric fields on biological media. Biological tissues and cell suspensions are described as having a heteregeneous permittivity and a heteregeneous conductivity. Well-posedness of the model problem and the regularity of its solution are established. A fully discrete finite element scheme is proposed for the numerical approximation of the potential distribution as a function of time and space simultaneously for an arbitrary shaped pulse, and it is demonstrated to enjoy the optimal convergence order in both space and time. The proposed numerical scheme has potential applications in the fields of medicine, food sciences, and biotechnology.

preprint2014arXiv

Direct Sampling Method for Diffusive Optical Tomography

In this work, we are concerned with the diffusive optical tomography (DOT) problem in the case when only one or two pairs of Cauchy data is available. We propose a simple and efficient direct sampling method (DSM) to locate inhomogeneities inside a homogeneous background and solve the DOT problem in both full and limited aperture cases. This new method is easy to implement and less expensive computationally. Numerical experiments demonstrate its effectiveness and robustness against noise in the data. This provides a new promising numerical strategy for the DOT problem.

preprint2014arXiv

Randomized Algorithms for Large-scale Inverse Problems with General Regularizations

We shall investigate randomized algorithms for solving large-scale linear inverse problems with general regularizations. We first present some techniques to transform inverse problems of general form into the ones of standard form, then apply randomized algorithms to reduce large-scale systems of standard form to much smaller-scale systems and seek their regularized solutions in combination with some popular choice rules for regularization parameters. Then we will propose a second approach to solve large-scale ill-posed systems with general regularizations. This involves a new randomized generalized SVD algorithm that can essentially reduce the size of the original large-scale ill-posed systems. The reduced systems can provide approximate regularized solutions with about the same accuracy as the ones by the classical generalized SVD, and more importantly, the new approach gains obvious robustness, stability and computational time as it needs only to work on problems of much smaller size. Numerical results are given to demonstrated the efficiency of the algorithms.

preprint2014arXiv

Super-resolution in Imaging High Contrast Targets from the Perspective of Scattering Coefficients

In this paper we consider the inverse scattering problem for high-contrast targets. We mathematically analyze the experimentally-observed phenomenon of super-resolution in imaging the target shape. This is the first time that a mathematical theory of super-resolution has been established in the context of imaging high contrast inclusions. We illustrate our main findings with a variety of numerical examples. Our analysis is based on the novel concept of scattering coefficients. These findings may help in developing resonant structures for resolution enhancement.

preprint2013arXiv

A New Splitting Method for Time-dependent Convection-dominated Diffusion Problems

We present a new splitting method for time-dependent convection-dominated diffusion problems. The original convection diffusion system is split into two sub-systems: a pure convection system and a diffusion system. At each time step, a convection problem and a diffusion problem are solved successively. The scheme has the following nice features: the convection subproblem is solved explicitly and a multistep technique is introduced to essentially enlarge the stability region so that the resulting scheme behaves like an unconditionally stable scheme; the diffusion subproblem is always self-adjoint and coercive so that it can be solved efficiently using many existing optimal preconditioned iterative solvers. The scheme is then extended for Navier-Stokes equations, where the nonlinear convection is resolved by a linear explicit multistep scheme at the convection step, and only a generalized Stokes problem is needed to solve at the diffusion step with the resulting stiffness matrix being invariant in the time marching process. The new schemes are all free from tuning some stabilization parameters for the convection-dominated diffusion problems. Numerical simulations are presented to demonstrate the stability, convergence and performance of the single-step and multistep variants of the new scheme.

preprint2013arXiv

Convergence of an Adaptive Finite Element Method for Distributed Flux Reconstruction

We shall establish the convergence of an adaptive conforming finite element method for the reconstruction of the distributed flux in a diffusion system. The adaptive method is based on a posteriori error estimators for the distributed flux, state and costate variables. The sequence of discrete solutions produced by the adaptive algorithm is proved to converge to the true triplet satisfying the optimality conditions in the energy norm and the corresponding error estimator converges to zero asymptotically.

preprint2013arXiv

Nearly Cloaking the Full Maxwell Equations: Cloaking Active Contents with General Conducting Layers

The regularized near-cloak via the transformation optics approach in the time-harmonic electromagnetic scattering is considered. This work extends the existing studies mainly in two aspects. First, it presents a near-cloak construction by incorporating a much more general conducting layer between the cloaked and cloaking regions. This might be of significant practical importance when production fluctuations occur. Second, it allows the cloaked contents to be both passive and active with an applied current inside. In assessing the near-cloaking performance, comprehensive and sharp estimates are derived for the scattering amplitude in terms of the asymptotic regularization parameter and the material tensors of the conducting layer. The scattering estimates are independent of the passive/active contents being cloaked, which implies that one could nearly cloak arbitrary contents by using the proposed near-cloak construction.

preprint2013arXiv

Optimal Shape Design by Partial Spectral Data

In this paper, we are concerned with a shape design problem, in which our target is to design, up to rigid transformations and scaling, the shape of an object given either its polarization tensor at multiple contrasts or the partial eigenvalues of its Neumann-Poincaré operator, which are known as the Fredholm eigenvalues. We begin by proposing to recover the eigenvalues of the Neumann-Poincaré operator from the polarization tensor by means of the holomorphic functional calculus. Then we develop a regularized Gauss-Newton optimization method for the shape reconstruction process. We present numerical results to demonstrate the effectiveness of the proposed methods and to illustrate important properties of the Fredholm eigenvalues and their associated eigenfunctions. Our results are expected to have important applications in the design of plasmon resonances in nanoparticles as well as in the multifrequency or pulsed imaging of small anomalies.

preprint2013arXiv

The Concept of Heterogeneous Scattering Coefficients and Its Application in Inverse Medium Scattering

This work investigates the scattering coefficients for inverse medium scattering problems. It shows some fundamental properties of the coefficients such as symmetry and tensorial properties. The relationship between the scattering coefficients and the far-field pattern is also derived. Furthermore, the sensitivity of the scattering coefficients with respect to changes in the permittivity and permeability distributions is investigated. In the linearized case, explicit formulas for reconstructing permittivity and permeability distributions from the scattering coefficients is proposed. They relate the exponentially ill-posed character of the inverse medium scattering problem at a fixed frequency to the exponential decay of the scattering coefficients. Moreover, they show the stability of the reconstruction from multifrequency measurements. This provides a new direction for solving inverse medium scattering problems.

preprint2012arXiv

A Direct Sampling Method for Inverse Electromagnetic Medium Scattering

In this paper, we study the inverse electromagnetic medium scattering problem of estimating the support and shape of medium scatterers from scattered electric or magnetic near-field data. We shall develop a novel direct sampling method based on an analysis of electromagnetic scattering and the behavior of the fundamental solution. The method is applicable even with one incident field and needs only to compute inner products of the measured scattered field with the fundamental solutions located at sampling points. Hence it is strictly direct, computationally very efficient, and highly tolerant to the presence of noise in the data. Two- and three-dimensional numerical experiments indicate that it can provide reliable support estimates of one single and multiple scatterers in case of both exact and highly noisy data.

preprint2012arXiv

A Direct Sampling Method for Inverse Scattering Using Far-Field Data

This work is concerned with a direct sampling method (DSM) for inverse acoustic scattering problems using far-field data. The method characterizes some unknown obstacles, inhomogeneous media or cracks, directly through an indicator function computed from the measured data. Using one or very few incident waves, the DSM provides quite reasonable profiles of scatterers in time-harmonic inverse acoustic scattering without a priori knowledge of either the physical properties or the number of disconnected components of the scatterer. We shall first derive the DSM using far-field data, then carry out a systematic evaluation of the performances and distinctions of the DSM using both near-field and far-field data. The numerical simulations are shown to demonstrate interesting and promising potentials of the DSM: a) ability to identify not only medium scatterers, but also obstacles, and even cracks, using measurement data from one or few incident directions, b) robustness with respect to large noise, and c) computational efficiency with only inner products involved.

preprint2012arXiv

A Multilevel Sampling Algorithm for Locating Inhomogeneous Media

In the reconstruction process of unknown multiple scattering objects in inverse medium scattering problems, the first important step is to effectively locate some approximate domains that contain all inhomogeneous media. Without such an effective step, one may have to take a much larger computational domain than actually needed in the reconstruction of all scattering objects, thus resulting in a huge additional computational efforts. In this work we propose a simple and efficient multilevel reconstruction algorithm to help locate an accurate position and shape of each inhomogeneous medium. Then other existing effective but computationally more demanding reconstruction algorithms may be applied in these initially located computational domains to achieve more accurate shapes of the scatter and the contrast values over each medium domain. The new algorithm exhibits several strengths: robustness against noise, requiring less incidences, fast convergence, flexibility to deal with scatterers of special shapes, and advantages in computational complexity.

preprint2012arXiv

A Two-stage Method for Inverse Medium Scattering

We present a novel numerical method to the time-harmonic inverse medium scattering problem of recovering the refractive index from near-field scattered data. The approach consists of two stages, one pruning step of detecting the scatterer support, and one resolution enhancing step with mixed regularization. The first step is strictly direct and of sampling type, and faithfully detects the scatterer support. The second step is an innovative application of nonsmooth mixed regularization, and it accurately resolves the scatterer sizes as well as intensities. The model is efficiently solved by a semi-smooth Newton-type method. Numerical results for two- and three-dimensional examples indicate that the approach is accurate, computationally efficient, and robust with respect to data noise.

preprint2012arXiv

Singular perturbation of reduced wave equation and scattering from an embedded obstacle

We consider time-harmonic wave scattering from an inhomogeneous isotropic medium supported in a bounded domain $Ω\subset\mathbb{R}^N$ ($N\geq 2$). {In a subregion $D\SubsetΩ$, the medium is supposed to be lossy and have a large mass density. We study the asymptotic development of the wave field as the mass density $ρ\rightarrow +\infty$} and show that the wave field inside $D$ will decay exponentially while the wave filed outside the medium will converge to the one corresponding to a sound-hard obstacle $D\SubsetΩ$ buried in the medium supported in $Ω\backslash\bar{D}$. Moreover, the normal velocity of the wave field on $\partial D$ from outside $D$ is shown to be vanishing as $ρ\rightarrow +\infty$. {We derive very accurate estimates for the wave field inside and outside $D$ and on $\partial D$ in terms of $ρ$, and show that the asymptotic estimates are sharp. The implication of the obtained results is given for an inverse scattering problem of reconstructing a complex scatterer.}