In this paper, we develop an efficient and energy stable fully-discrete local discontinuous Galerkin (LDG) method for the Cahn–Hilliard–Hele–Shaw (CHHS) system. The semi-discrete energy stability of the LDG method is proved firstly. Due to the strict time step restriction ( t = O( x4)) of explicit time discretization methods for stability, we introduce a semi- implicit time integration scheme which is based on a convex splitting of the discrete Cahn– Hilliard energy. The unconditional energy stability has been proved for this fully-discrete LDG scheme. The fully-discrete equations at the implicit time level are nonlinear. Thus, the nonlinear Full Approximation Scheme (FAS) multigrid method has been applied to solve this system of algebraic equations, which has been shown the nearly optimal complexity numerically. Numerical results are also given to illustrate the accuracy and capability of the LDG method coupled with the multigrid solver.

We propose an efficient threshold dynamics method for topology optimization for fluids modeled with the Stokes equation. The proposed algorithm is based on minimization of an objective energy function that consists of the dissipation power in the fluid and the perimeter approximated by nonlocal energy, subject to a fluid volume constraint and the incompressibility condition. We show that the minimization problem can be solved with an iterative scheme in which the Stokes equation is approximated by a Brinkman equation. The indicator functions of the fluid-solid regions are then updated according to simple convolutions followed by a thresholding step. We demonstrate mathematically that the iterative algorithm has the total energy decaying property. The proposed algorithm is simple and easy to implement. A simple adaptive time strategy is also used to accelerate the convergence of the iteration. Extensive numerical experiments in both two and three dimensions show that the proposed iteration algorithm is quite robust and converges in much fewer iterations and is more efficient than many existing methods.

In this paper we prove a new matrix Li-Yau-Hamilton (LYH) estimate for K¨ahler-Ricci flow on manifolds with nonnegative bi- sectional curvature. The form of this new LYH estimate is obtained by the interpolation consideration originated in [Ch] by Chow. This new inequality is shown to be connected with Perelman’s entropy formula through a family of differential equalities. In the rest of the paper, we show several applications of this new estimate and its corresponding estimate for linear heat equation. These include a sharp heat kernel comparison theorem, generalizing the earlier result of Li and Tian [LT], a manifold version of Stoll’s theorem [St] on the characterization of ‘algebraic divisors’, and a localized monotonicity formula for analytic subvari- eties, which sharpens the Bishop volume comparison theorem. Motivated by the connection between the heat kernel estimate and the reduced volume monotonicity of Perelman [P], we prove a sharp lower bound of the fundamental solution to the forward conjugate heat equation, which in a certain sense dual to Perelman’s monotonicity of the reduced volume. As an application of this new monotonicity formula, we show that the blow-down limit of a certain type of long-time solution is a gradient expanding soliton, generalizing an earlier result of Cao. We also illustrate the connection between the new LYH estimate and the Hessian comparison theorem of [FIN] on the forward reduced distance. Localized monotonicity formulae on entropy and forward reduced volume are also derived.

Solitary waves bifurcated from edges of Bloch bands in two-dimensional periodic media are determined both analytically and numerically in the context of a two-dimensional nonlinear Schrödinger equation with a periodic potential. Using multiscale perturbation methods, the envelope equations of solitary waves near Bloch bands are analytically derived. These envelope equations reveal that solitary waves can bifurcate from edges of Bloch bands under either focusing or defocusing nonlinearity, depending on the signs of the second-order dispersion coefficients at the edge points. Interestingly, at edge points with two linearly independent Bloch modes, the envelope equations lead to a host of solitary wave structures, including reduced-symmetry solitons, dipole-array solitons, vortex-cell solitons, and so on—many of which have not been reported before to our knowledge. It is also shown analytically that the centers of envelope solutions can be positioned at only four possible locations at or between potential peaks. Numerically, families of these solitary waves are directly computed both near and far away from the band edges. Near the band edges, the numerical solutions spread over many lattice sites, and they fully agree with the analytical solutions obtained from the envelope equations. Far away from the band edges, solitary waves are strongly localized, with intensity and phase profiles characteristic of individual families.

In the present paper, we prove the 2-part of Birch and Swinnerton-Dyer conjecture for an explicit infinite family of rank 0 quadratic twists of the modular elliptic curve χ_0(14), using an explicit form of the Waldspurger formula. We also give an explicit infinite family of rank 1 quadratic twists of χ_0(14) whose Tate–Shafarevich group is of odd cardinality.