In this paper, we study the problem of computing the effective diffusivity for a particle moving in chaotic and stochastic flows. In addition, we numerically investigate the residual diffusion phenomenon in chaotic advection. The residual diffusion refers to the nonzero effective (homogenized) diffusion in the limit of zero molecular diffusion as a result of chaotic mixing of the streamlines. In this limit, traditional numerical methods typically fail since the solutions of the advection-diffusion equations develop sharp gradients. Instead of solving the Fokker--Planck equation in the Eulerian formulation, we compute the motion of particles in the Lagrangian formulation, which is modeled by stochastic differential equations (SDEs). We propose an effective numerical integrator based on a splitting method to solve the corresponding SDEs in which the deterministic subproblem is symplectic preserving while the random subproblem

A semi-explicit formula of solution to the boundary layer system for thermal layer derived from the compressible Navier-Stokes equations with the non-slip boundary condition when the viscosity coefficients vanish is given, in particular in three space dimension. In contrast to the inviscid Prandtl system studied by [7] in two space dimension, the main difficulty comes from the coupling of the velocity field and the temperature field through a degenerate parabolic equation. The convergence of these boundary layer equations to the inviscid Prandtl system is justified when the initial temperature goes to a constant. Moreover, the time asymptotic stability of the linearized system around a shear flow is given, and in particular, it shows that in three space dimension, the asymptotic stability depends on whether the direction of tangential velocity field of the shear flow is invariant in the normal direction respective to the boundary.

No abstract uploaded!

We introduce a new technique for isolating components on the spectral side of the trace formula. By applying it to the Jacquet--Rallis relative trace formula, we complete the proof of the global Gan--Gross--Prasad conjecture and its refinement Ichino--Ikeda conjecture for U(n)×U(n+1) in the stablecase.

For any vertex algebra V and any subalgebra A of V, there is a new subalgebra of V known as the commutant of A in V. This construction was introduced by Frenkel-Zhu, and is a generalization of an earlier construction due to Kac-Peterson and Goddard-Kent-Olive known as the coset construction. In this paper, we interpret the commutant as a vertex algebra notion of invariant theory. We present an approach to describing commutant algebras in an appropriate category of vertex algebras by reducing the problem to a question in commutative algebra. We give an interesting example of a Howe pair (ie, a pair of mutual commutants) in the vertex algebra setting.

The purpose of this paper is to study the dynamics of the interaction among a special class of solutions of the one-dimensional Camassa–Holm equation. The equation yields soliton solutions whose identity is preserved through nonlinear interactions. These solutions are characterized by a discontinuity at the peak in the wave shape and are thus called peakon solutions. We apply a particle method to the Camassa–Holm equation and show that the nonlinear interaction among the peakon solutions resembles an elastic collision, i.e., the total energy and momentum of the system before the peakon interaction is equal to the total energy and momentum of the system after the collision. From this result, we provide several numerical illustrations which support the analytical study, as well as showcase the merits of using a particle method to simulate solutions to the Camassa–Holm equation under a wide class of initial data.

The convergence of a fourth order finite difference method for the 2-D unsteady, viscous incompressible Boussinesq equations, based on the vorticity-stream function formulation, is established in this article. A compact fourth order scheme is used to discretize the momentum equation, and long-stencil fourth order operators are applied to discretize the temperature transport equation. A local vorticity boundary condition is used to enforce the no-slip boundary condition for the velocity. One-sided extrapolation is used near the boundary, dependent on the type of boundary condition for the temperature, to prescribe the temperature at ‘‘ghost’’ points lying outside of the computational domain. Theoretical results of the stability and accuracy of the method are also provided. In numerical experiments the method has been shown to be capable of producing highly resolved solutions at a reasonable computational cost.

We study the collapsing behavior of the Kähler–Ricci flow on a compact Kähler manifold X admitting a holomorphic submersion inherited from its canonical bundle. We show that the flow metric degenerates at exactly the rate of e^{-t} as predicted by the cohomology information, and so the fibres collapse at the optimal rate. Consequently, it leads to some analytic and geometric extensions to the regular case of works by J. Song and G. Tian. Its applicability to general Calabi–Yau fibrations will also be discussed in local settings.

We study the Kahler-Ricci flow on a class of projective bundles over a compact Kahler-Einstein manifold. Assuming the initial Kahler metric admits a U(1)-invariant momentum profile, we give a criterion, characterized by the triple (Σ,L,[ω0]), under which the P1-fiber collapses along the Kahler-Ricci flow and the projective bundle converges to Σ in the Gromov-Hausdorff sense. Furthermore, the Kahler-Ricci flow must have Type I singularity and is of (Cn ×P1)-type. This generalizes and extends part of Song-Weinkove’s work on Hirzebruch surfaces.

Simulation of electromagnetic wave propagation in metamaterials leads to more complicated time domain Maxwell's equations than the standard Maxwell's equations in free space. In this paper, we develop and analyze a non-dissipative discontinuous Galerkin (DG) method for solving the Maxwell's equations in Drude metamaterials. Previous discontinuous Galerkin methods in the literature for electromagnetic wave propagation in metamaterials were either non-dissipative but sub-optimal, or dissipative and optimal. Our method uses a different and simple choice of numerical fluxes, achieving provable non-dissipative stability and optimal error estimates simultaneously. We prove the stability and optimal error estimates for both semi- and fully discrete DG schemes, with the leap-frog time discretization for the fully discrete case. Numerical results are given to demonstrate that the DG method can solve metamaterial Maxwell's equations effectively.