The focus of this work is on rigorous mathematical analysis of the topological derivative based detection algorithms for the localization of an elastic inclusion of vanishing characteristic size. A filtered quadratic misfit is considered, and the performance of the topological derivative imaging functional resulting therefrom is analyzed. Our analysis reveals that the imaging functional may not attain its maximum at the location of the inclusion. Moreover, the resolution of the image is below the diffraction limit. Both phenomena are due to the coupling of pressure and shear waves propagating with different wave speeds and polarization directions. A novel imaging functional based on the weighted Helmholtz decomposition of the topological derivative is, therefore, introduced. It is thereby substantiated that the maximum of the imaging functional is attained at the location of the inclusion and the resolution is enhanced and proves to be the diffraction limit. Finally, we investigate the stability of the proposed imaging functionals with respect to measurement and medium noises.

We consider reflector imaging in a weakly random waveguide. We address the situation in which the source is farther from the reflector to be imaged than the energy equipartition distance, but the receiver array is closer to the reflector to be imaged than the energy equipartition distance. As a consequence, the reflector is illuminated by a partially coherent field and the signals recorded by the receiver array are noisy. This paper shows that migration of the recorded signals cannot give a good image, but an appropriate migration of the cross correlations of the recorded signals can give a very good image. The resolution and stability analysis of this original functional shows that the reflector can be localized with an accuracy of the order of the wavelength even when the receiver array has small aperture, and that broadband sources are necessary to ensure statistical stability, whatever the aperture of the array.

This paper analyzes the random fluctuations obtained by a heterogeneous multi-scale first-order finite element method applied to solve elliptic equations with a random potential. Several multiscale numerical algorithms have been shown to correctly capture the homogenized limit of solutions of elliptic equations with coefficients modeled as stationary and ergodic random fields. Because theoretical results are available in the continuum setting for such equations, we consider here the case of a secondorder elliptic equations with random potential in two dimensions of space. We show that the random fluctuations of such solutions are correctly estimated by the heterogeneous multi-scale algorithm when appropriate fine-scale problems are solved on subsets that cover the whole computational domain. However, when the fine-scale problems are solved over patches that do not cover the entire domain, the random fluctuations may or may not be estimated accurately. In the case of random potentials with short-range interactions, the variance of the random fluctuations is amplified as the inverse of the fraction of the medium covered by the patches. In the case of random potentials with long-range interactions, however, such an amplification does not occur and random fluctuations are correctly captured independent of the (macroscopic) size of the patches. These results are consistent with those obtained in [9] for more general equations in the one-dimensional setting and provide indications on the loss in accuracy that results from using coarser, and hence computationally less intensive, algorithms.

The aim of this paper is to develop a mathematical framework for opto-elastography. In opto-elastography, a mechanical perturbation of the medium produces a decorrelation of optical speckle patterns due to the displacements of optical scatterers. To model this, we consider two optically random media, with the second medium obtained by shifting the first medium in some local region. We derive the radiative transfer equation for the cross-correlation of the wave fields in the media. Then we derive its diffusion approximation. In both the radiative transfer and the diffusion regimes, we relate the correlation of speckle patterns to the solutions of the radiative transfer and the diffusion equations. We present numerical simulations based on our model which are in agreement with recent experimental measurements.

The aim of this paper is to provide a fast and efficient procedure for (real-time) target identification in imaging based on matching on a dictionary of precomputed generalized polarization tensors (GPTs). The approach is based on some important properties of the GPTs and new invariants. A new shape representation is given and numerically tested in the presence of measurement noise. The stability and resolution of the proposed identification algorithm is numerically quantified. We compare the proposed GPT-based shape representation with a moment-based one.

This paper aims to mathematically advance the field of quantitative thermo-acoustic imaging. Given several electromagnetic data sets, we establish for the first time an analytical formula for reconstructing the absorption coefficient from thermal energy measurements. Since the formula involves derivatives of the given data up to the third order, it is unstable in the sense that small measurement noises may cause large errors. However, in the presence of measurement noise, the obtained formula, together with a noise regularization technique, provides a good initial guess for the true absorption coefficient. We finally correct the errors by deriving a reconstruction formula based on the least square solution of an optimal control problem and prove that this optimization step reduces the errors occurring and enhances the resolution.

We consider the optimization approach to the acousto-electric imaging problem. Assuming that the electric conductivity distribution is a small perturbation of a constant, we investigate the least-squares estimate analytically using (multiple) Fourier series, and confirm the widely observed fact that acousto-electric imaging has high resolution and is statistically stable. We also analyze the case of partial data and the case of limited-view data, in which some singularities of the conductivity can still be imaged.

We consider an elliptic pseudo-differential equation with a highly oscillating linear potential modeled as a stationary ergodic random field. The random field is a function composed with a centered long-range correlated Gaussian process. In the limiting of vanishing correlation length, the heterogeneous solution converges to a deterministic solution obtained by averaging the random potential. We characterize the deterministic and stochastic correctors. With proper rescaling, the mean-zero stochastic corrector converges to a Gaussian random process in probability and weakly in the spatial variables. In addition, for two prototype equations involving the Laplacian and the fractional Laplacian operators, we prove that the limit holds in distribution in some Hilbert spaces. We also determine the size of the deterministic corrector when it is larger than the stochastic corrector. Depending on the correlation structure of the random field and on the singularities of the Green’s function, we show that either the deterministic or the random part of the corrector dominates.

We consider the problem of the random fluctuations in the solutions to elliptic PDEs with highly oscillatory random coefficients. In our setting, as the correlation length of the fluctuations tends to zero, the heterogeneous solution converges to a deterministic solution obtained by averaging. When the Green’s function to the unperturbed operator is sufficiently singular (i.e., not square integrable locally), the leading corrector to the averaged solution may be either deterministic or random, or both in a sense we shall explain. Our main application is the solution of an elliptic problem with random Robin boundary condition that may be used to model diffusion of signaling molecules through a layer of cells into a bulk of extracellular medium. The problem is then described by an elliptic pseudo-differential operator (a Dirichlet-to-Neumann operator) on the boundary of the domain with random potential. In the physical setting of a three dimensional extracellular medium on top of a two-dimensional surface of cells forming a layer of epithelium, we show that the approximate corrector to averaging consists of a deterministic correction plus a Gaussian field of amplitude proportional to the correlation length of the random medium. The result is obtained under some assumptions on the four-point correlation function in the medium. We provide examples of such random media based on Gaussian and Poisson statistics.

We consider the theory of correctors to homogenization in stationary transport equations with rapidly oscillating, random coefficients. Let ε << 1 be the ratio of the correlation length in the random medium to the overall distance of propagation. As ε↓0, we show that the heterogeneous transport solution is well-approximated by a homogeneous transport solution. We then show that the rescaled corrector converges in (probability) distribution and weakly in the space and velocity variables, to a Gaussian process as an application of a central limit result. The latter result requires strong assumptions on the statistical structure of randomness and is proved for random processes constructed by means of a Poisson point process.

We consider the stabilization of linear kinetic equations with a random relaxation term. The well-known framework of hypocoercivity by J. Dolbeault, C. Mouhot and C. Schmeiser (2015) ensures the stability in the deterministic case. This framework, however, cannot be applied directly for arbitrarily small random relaxation parameters. Therefore, we introduce a Galerkin formulation, which reformulates the stochastic system as a sequence of deterministic ones. We prove for the γ-distribution that the hypocoercivity framework ensures the stability of this series and hence the stochastic stability of the underlying random kinetic equation. The presented approach also yields a convergent numerical approximation.

In this work, we consider the numerical solution of the nonlinear Schrödinger equation with a highly oscillatory potential (NLSE-OP). The NLSE-OP is a model problem which frequently occurs in recent studies of some multiscale dynamical systems, where the potential introduces wide temporal oscillations to the solution and causes numerical difficulties. We aim to analyze rigorously the error bounds of the splitting schemes for solving the NLSE-OP to a fixed time. Our theoretical results show that the Lie–Trotter splitting scheme is uniformly and optimally accurate at the first order provided that the oscillatory potential is integrated exactly, while the Strang splitting scheme is not. Our results apply to general dispersive or wave equations with an oscillatory potential. The error estimates are confirmed by numerical results.

We introduce two exponential-type integrators for the “good” Boussinesq equation. They are of orders one and two, respectively, and they require lower spatial regularity of the solution compared to classical exponential integrators. For the first integrator, we prove first-order convergence in H^r for solutions in H^{r+1} with r>1/2. This new integrator even converges (with lower order) in H^r for solutions in H^r, i.e., without any additional smoothness assumptions. For the second integrator, we prove second-order convergence in H^r for solutions in H^{r+3} with r>1/2 and convergence in L^2 for solutions in H^3. Numerical results are reported to illustrate the established error estimates. The experiments clearly demonstrate that the new exponential-type integrators are favorable over classical exponential integrators for initial data with low regularity.

We present and analyze two numerical methods for the logarithmic Schrödinger equation (LogSE) consisting of a regularized splitting method and a regularized conservative Crank–Nicolson finite difference method (CNFD). In order to avoid numerical blow-up and/or to suppress round-off error due to the logarithmic nonlinearity in the LogSE, a regularized logarithmic Schrödinger equation (RLogSE) with a small regularized parameter 0<ε≪1 is adopted to approximate the LogSE with linear convergence rate O(ε). Then we use the Lie–Trotter splitting integrator to solve the RLogSE and establish its error bound O(τ^{1/2}ln(ε^{−1})) with τ>0 the time step, which implies an error bound at O(ε+τ^{1/2}ln(ε^{−1})) for the LogSE by the Lie–Trotter splitting method. In addition, the CNFD is also applied to discretize the RLogSE, which conserves the mass and energy in the discretized level. Numerical results are reported to confirm our error bounds and to demonstrate rich and complicated dynamics of the LogSE.

We consider the nonlinear Schr{ö}dinger equation with a defocusing nonlinearity which is mass-(super)critical and energy-subcritical. We prove uniform in time error estimates for the Lie-Trotter time splitting discretization. This uniformity in time is obtained thanks to a vectorfield which provides time decay estimates for the exact and numerical solutions. This vectorfield is classical in scattering theory, and requires several technical modifications compared to previous error estimates for splitting methods.

We consider the Schrödinger equation with a nondispersive logarithmic nonlinearity and a repulsive harmonic potential. For a suitable range of the coefficients, there exist two positive stationary solutions, each one generating a continuous family of solitary waves. These solutions are Gaussian, and turn out to be orbitally unstable. We also discuss the notion of ground state in this setting: for any natural definition, the set of ground states is empty.

We consider the kinetic Fokker-Planck equation with a class of general force. We prove the existence and uniqueness of a positive normalized equilibrium (in the case of a gen- eral force) and establish some exponential rate of convergence to the equilibrium (and the rate can be explicitly computed). Our results improve results about classical force to general force case. Our result also improve the rate of convergence for the Fitzhugh-Nagumo equation from non-quantitative to quantitative explicit rate.

No abstract uploaded!

No abstract uploaded!

No abstract uploaded!

We propose a conservation-dissipation formalism (CDF) for coarse-grained descriptions of irreversible processes. This formalism is based on a stability criterion for non-equilibrium thermodynamics. The criterion ensures that non-equilibrium states tend to equilibrium in long time. As a systematic methodology, CDF provides a feasible procedure in choosing non-equilibrium state variables and determining their evolution equations. The equations derived in CDF have a unified elegant form. They are globally hyperbolic, allow a convenient definition of weak solutions, and are amenable to existing numerics. More importantly, CDF is a genuinely nonlinear formalism and works for systems far away from equilibrium. With this formalism, we formulate novel thermodynamics theories for heat conduction in rigid bodies and non-isothermal compressible Maxwell fluid flows as two typical examples. In these examples, the non-equilibrium variables are exactly the conjugate variables of the heat fluxes or stress tensors. The new theory generalizes Cattaneo’s law or Maxwell’s law in a regularized and nonlinear fashion.

In this article, we propose a novel approach to construct macroscopic balance equations and constitutive equations describing various irreversible phenomena. It is based on the general principles of non-equilibrium thermodynamics and consists of four basic steps: picking suitable state variables, choosing a strictly concave entropy function, properly separating entropy fuxes and production rates, and determining a dissipation matrix. Our approach takes advantage of both extended irreversible thermodynamics and GENERIC formalisms and shows a direct correspondence with Levermore’s moment-closure hierarchies for the Boltzmann equation. As a direct application, a new ten-moment model beyond the classical hierarchies is constructed and is shown to recover the Euler equations in the equilibrium state. These interesting results may put various macroscopic modeling approaches, starting from the general principles of non-equilibrium thermodynamics, on a solid microscopic foundation based on the Boltzmann equation.

We investigate the fractional Schrödinger equation with a periodic PT-symmetric potential. In the inverse space, the problem transfers into a first-order nonlocal frequency-delay partial differential equation. We show that at a critical point, the band structure becomes linear and symmetric in the one-dimensional case, which results in a nondiffracting propagation and conical diffraction of input beams. If only one channel in the periodic potential is excited, adjacent channels become uniformly excited along the propagation direction, which can be used to generate laser beams of high power and narrow width. In the two-dimensional case, there appears conical diffraction that depends on the competition between the fractional Laplacian operator and the PT-symmetric potential. This investigation may find applications in novel on-chip optical devices.

In this paper we study the zero-viscosity limit of 2-D Boussinesq equations with vertical viscosity and zero diffusivity, which is a nonlinear system with partial dissipation arising in atmospheric sciences and oceanic circulation. The domain is taken as R2+ with Navier-type boundary. We prove the nonlinear stability of the approximate solution constructed by boundary layer expansion in conormal Sobolev space. The optimal expansion order and convergence rates for the inviscid limit are also identified in this paper. Our paper extends the partial zero-dissipation limit results of Boussinesq system with full dissipation by Chae D. [Adv.Math.203,no.2,2006] in the whole space to the case with partial dissipation and Navier boundary in the half plane.

We prove the nonlinear asymptotic stability of the gravitational hydrostatic equilibrium for the general equation of state of pressure-density relation in the framework of vacuum free boundary problem of spherically symmetric compressible Navier-Stokes-Poisson equations in three dimensions.The results apply to white dwarfs and polytropes with γ≥2, for which even the linearized asymptotic stability results are not available in literature, to the best of the authors' knowledge. Detailed decay rates of perturbations are given.