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.

Given a strictly hyperbolic, genuinely nonlinear system of conservation laws, we prove the a priori bound u (t,) u (t,)= O (1)(1+ t) | ln | on the distance between an exact BV solution u and a viscous approximation u, letting the viscosity coefficient 0. In the proof, starting from u we construct an approximation of the viscous solution u by taking a mollification u* and inserting viscous shock profiles at the locations of finitely many large shocks for each fixed . Error estimates are then obtained by introducing new Lyapunov functionals that control interactions of shock waves in the same family and also interactions of waves in different families. 2004 Wiley Periodicals, Inc.

The purpose of this research monograph is to survey some recent developments in the analysis of shock reflection-diffraction, to present our original mathematical proofs of von Neumann’s conjectures for potential flow, to collect most of the related results and new techniques in the analysis of partial differential equations achieved in the last decades, and to discuss a set of fundamental open problems relevant to the directions of future research in this and related areas.

For systems of hyperbolic conservation laws, a new Glimm functional was recently constructed when the linearly degenerate manifold in each characteristic field is either the whole space or it consists of a finitely many smooth and transversal manifolds of codimension one. This new functional leads to the neat consistency and convergence rate estimation of the Glimm scheme. In this paper, by the motivation of the result in Ref. 2, we show that the corresponding new Glimm functional can be constructed for general systems only under the strict hyperbolicity assumption.

We characterize probability measure with finite moment of any order in terms of the symmetric difference operators of their Fourier transforms. By using our new characterization, we prove the continuity f (t, v)\in C ((0,\infty), L^ 1_ {2k-2+ }) , where f (t, v)\in C ((0,\infty), L^ 1_ {2k-2+ }) stands for the density of unique measure-valued solution f (t, v)\in C ((0,\infty), L^ 1_ {2k-2+ }) of the Cauchy problem for the homogeneous non-cutoff Boltzmann equation, with Maxwellian molecules, corresponding to a probability measure initial datum f (t, v)\in C ((0,\infty), L^ 1_ {2k-2+ }) satisfying\[\int| v|^{2k-2+ } dF_0 (v)<\infty, 0\leq < 2, k= 2, 3, 4,\cdots\] provided that f (t, v)\in C ((0,\infty), L^ 1_ {2k-2+ }) is not a single Dirac mass.

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.

Speed ensemble of bistable (combustion) fronts in mean zero stationary Gaussian shear flows inside two and three dimensional channels is studied with a min-max variational principle. In the small root mean square regime of shear flows, a new class of multi-scale test functions are found to yield speed asymptotics. The quadratic speed enhancement law holds with probability arbitrarily close to one under the almost sure continuity (dimension two) and mean square Hlder regularity (dimension three) of the shear flows. Remarks are made on the quadratic and linear laws of front speed expectation in the small and large root mean square regimes.

In this paper, we study the asymptotic stability of rarefaction waves for the compressible isentropic NavierStokes equations with density-dependent viscosity. First, a weak solution around a rarefaction wave to the Cauchy problem is constructed by approximating the system and regularizing the initial values which may contain vacuum states. Then some global in time estimates on the weak solution are obtained. Based on these uniform estimates, the vacuum states are shown to vanish in finite time and the weak solution we constructed becomes a unique strong one. Consequently, the stability of the rarefaction wave is proved in a weak sense. The theory holds for large-amplitudes rarefaction waves and arbitrary initial perturbations.

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.

In this article, we study wave dynamics in the fractional nonlinear Schrödinger equation with a modulated honeycomb potential. This problem arises from recent research interests in the interplay between topological materials and nonlocal governing equations. Both are current focuses in scientific research fields. We first develop the Floquet-Bloch spectral theory of the linear fractional Schrödinger operator with a honeycomb potential. Especially, we prove the existence of conical degenerate points, i.e., Dirac points, at which two dispersion band functions intersect. We then investigate the dynamics of wave packets spectrally localized at a Dirac point and derive the leading effective envelope equation. It turns out the envelope can be described by a nonlinear Dirac equation with a varying mass. With rigorous error estimates, we demonstrate that the asymptotic solution based on the effective envelope equation approximates the true solution well in the weighted-H^s space.

Under the diffusion scaling and a scaling assumption on the microscopic component, a non-classical fluid dynamic system was derived in\cite {BGLY} that is related to the system of ghost effect derived in\cite {Sone-2} in different settings. This paper aims to justify this limit system for a non-trivial background profile with slab symmetry. The result reveals not only the diffusion phenomena in the temperature and density, but also the flow of higher order in Knudsen number due to the gradient of the temperature. Precisely, we show that the solution to the Boltzmann equation converges to a diffusion wave with decay rates in both Knudsen number and time.

As a continuation of our series works on the Boltzmann equation without angular cutoff assumption, in this part, the global existence of solution to the Cauchy problem in the whole space is proved in some suitable weighted Sobolev spaces for hard potential when the solution is a small perturbation of a global equilibrium.

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.

No abstract uploaded!

In [9] Anderson’s conjecture was proven by comparing values of Bloch functions with the variation of the function. We extend that result on Bloch functions from two to arbitrary dimension and prove that ∫_{[0,x]}|∇b(ζ)|e^{b(ζ)}d|ζ|<∞.. In the second part of the paper, we show that the area or volume integral ∫_{B^d}|∇u(w)|p(w,θ)dA(w) for positive harmonic functions u is bounded by the value cu(0) for at least one θ. The integral is also transferred to simply connected domains and interpreted from the point of view of stochastics. Several emerging open problems are presented.

For the class of graph-directed self-similar measures on R, which could have overlaps but are essentially of finite type, we set up a framework for deriving a closed formula for the spectral dimension of the Laplacians defined by these measures. For the class of finitely ramified graph-directed self-similar sets, the spectral dimension of the associated Laplace operators has been obtained by Hambly and Nyberg [6]. The main novelty of our results is that the graph-directed self-similar measures we consider do not need to satisfy the graph open set condition.

We study linear and nonlinear Schrodinger equations defined by fractal measures. Under the assumption that the Laplacian has compact resolvent, we prove that there exists a uniqueness weak solution for a linear Schrodinger equation, and then use it to obtain the existence and uniqueness of weak solution of a nonlinear Schrodinger equation. We prove that for a class of self-similar measures on R with overlaps, the Schrodinger equations can be discretized so that the finite element method can be applied to obtain approximate solutions. We also prove that the numerical solutions converge to the actual solution and obtain the rate of convergence.

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.

Edge states attract more and more research interests owing to the novel topologically protected properties. In this work, we studied edge modes and traveling edge states via the linear Dirac equation with so-called edge-admissible masses. The unidirectional edge state provides a heuristic approach to more general traveling edge states through the localized behavior along slowly varying edges. We show the dominated asymptotic solutions of two typical edge states that follow circular and curved edges with small curvature by the analytic and quantitative arguments.

Mathematical analysis on electromagnetic waves in photonic graphene, a photonic topological material which has a honeycomb structure, is one of the most important current research topics. By modulating the honeycomb structure, numerous topological phenomena have been observed recently. The electromagnetic waves in such a media are generally described by the 2-dimensional wave equation. It has been shown that the corresponding elliptic operator with a honeycomb material weight has Dirac points in its dispersion surfaces. In this paper, we study the time evolution of the wave packets spectrally concentrated at such Dirac points in a modulated honeycomb material weight. We prove that such wave packet dynamics is governed by the Dirac equation with a varying mass in a large but finite time. Our analysis provides mathematical insights to those topological phenomena in photonic graphene.

This article is a sequel to our previous work [13] concerned with the derivation of high-order homogenized models for the Stokes equation in a periodic porous medium. We provide an improved asymptotic analysis of the coefficients of the higher order models in the low-volume fraction regime whereby the periodic obstacles are rescaled by a factor η which converges to zero. By introducing a new family of order k corrector tensors with a controlled growth as η→0 uniform in k∈\N, we are able to show that both the infinite order and the finite order models converge in a coefficient-wise sense to the three classical asymptotic regimes. Namely, we retrieve the Darcy model, the Brinkman equation or the Stokes equation in the homogeneous cubic domain depending on whether η is respectively larger, proportional to, or smaller than the critical size η_{crit}∼ϵ2/(d−2). For completeness, the paper first establishes the analogous results for the perforated Poisson equation, considered as a simplified scalar version of the Stokes system.

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 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!

In this paper we consider the well-posedness of the compressible Prandtl boundary layer for the Navier–Stokes–Poisson equations with nonslip boundary condition and obtain the linear stability of the approximate solution constructed by boundary layer expansion.