A key difficulty in the analysis and numerical approximation of the shallow water equations is the nonconservative product of measures due to the gravitational force acting on a sloped bottom. Solutions may be nonunique, and numerical schemes are not only consistent discretizations of the shallow water equations, but they also determine how to model the physics. Our derivation is based on a subcell reconstruction using infinitesimal singular layers at the cell boundaries, as inspired by S. Noelle, Y. Xing, and C.-W. Shu [J. Comput. Phys., 226 (2007), pp. 29-58]. One key step is to separate the singular measures. Another aspect is the reconstruction of the solution variables in the singular layers. We study three reconstructions. The first leads to the well-known scheme of Audusse et al., [SIAM J. Sci. Comput., 25 (2004), pp. 2050-2065], which introduces the hydrostatic reconstruction. The second is a modification proposed in [T. Morales de Luna, M. J. Castro D´ıaz, and C. Par´es, Appl. Math. Comput., 219 (2013), pp. 9012-9032], which analyzes whether a wave has enough energy to overcome a step. The third is our new scheme, which borrows its structure from the wet-dry front. For a number of cases discussed in recent years, where water runs down a hill, Audusse’s scheme converges slowly or fails. Morales’ scheme gives a visible improvement. Both schemes are clearly outperformed by our new scheme.

No abstract uploaded!

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.

In the paper, we study the Prandtl system with initial data admitting non-degenerate critical points. For any index \sigma\in [3/2, 2], we obtain the local in time well-posedness in the space of Gevrey class \sigma\in [3/2, 2], in the tangential variable and Sobolev class in the normal variable so that the monotonicity condition on the tangential velocity is not needed to overcome the loss of tangential derivative. This answers the open question raised in the paper of D. Grard-Varet and N. Masmoudi [{\it Ann. Sci. c. Norm. Supr}.(4) 48 (2015), no. 6, 1273-1325], in which the case \sigma\in [3/2, 2], is solved.

We establish several gradient estimates for second-order divergence type parabolic and elliptic systems. The coefficients and data are assumed to be Hlder or Dini continuous in the time variable and all but one spatial variable. This type of system arises from the problems of linearly elastic laminates and composite materials. For the proof, we use Campanatos approach in a novel way. Non-divergence type equations under similar conditions are also discussed.