MathSciDoc: An Archive for Mathematician ∫

In this paper, we establish the optimal error estimate of the particle method for a family of nonlinear evolutionary partial differential equations, or the so-called b-equation. The b-equation, including the Camassa-Holm equation and the Degasperis-Procesi equation, has many applications in diverse scientific fields. The particle method is an approximation of the b-equation in Lagrangian representation. We also prove short-time existence, uniqueness and regularity of the Lagrangian representation of the b-equation.

In biological transport mechanisms such as insect respiration and renal filtration, fluid travels along a leaky channel allowing material exchange with systems exterior to the channel. The channels in these systems may undergo peristaltic pumping which is thought to enhance the material exchange. To date, little analytic work has been done to study the effect of pumping on material extraction across the channel walls. In this paper, we examine a fluid extraction model in which fluid flowing through a leaky channel is exchanged with fluid in a reservoir. The channel walls are allowed to contract and expand uniformly, simulating a pumping mechanism. In order to efficiently determine solutions of the model, we derive a formal power series solution for the Stokes equations in a finite channel with uniformly contracting/expanding permeable walls. This flow has been well studied in the case in which the normal velocity at the channel walls is proportional to the wall velocity. In contrast we do not assume flow that is proportional to the wall velocity, but flow that is driven by hydrostatic pressure, and we use Darcy’s law to close our system for normal wall velocity. We incorporate our flow solution into a model that tracks the material pressure exterior to the channel. We use this model to examine flux across the channel-reservoir barrier and demonstrate that pumping can either enhance or impede fluid extraction across channel walls. We find that associated with each set of physical flow and pumping parameters, there are optimal reservoir conditions that maximize the amount of material flowing from the channel into the reservoir.

Abstract In this note we establish the uniform (Formula presented.)-bound for the weak solutions to a degenerate Keller-Segel equation with the diffusion exponent (Formula presented.) under a sharp condition on the initial data for the global existence. As a consequence, the uniqueness of the weak solutions is also proved.

This note considers the hyper-contractivity and the uniqueness of the weak solutions to the two dimensional Keller–Segel equations. We prove a refined hyper-contractive property and consequently obtain the uniqueness of global weak solutions provided that the initial data satisfy , and . We also extend the result to higher dimensions.

We derive an exact solution for Stokes flow in an in a channel with permeable walls. We assume that at the channel walls, the normal component of the fluid velocity is described by Darcy's law and the tangential component of the fluid velocity is described by the no slip condition. The pressure exterior to the channel is assumed to be constant. Although this problem has been well studied, typical studies assume that the permeability of the wall is small relative to other non-dimensional parameters; this work relaxes this assumption and explores a regime in parameter space that has not yet been well studied. A consequence of this relaxation is that transverse velocity is no longer necessarily small when compared with the axial velocity. We use our result to explore how existing asymptotic theories break down in the limit of large permeability.

In this paper, we apply a simple finite element numerical scheme, proposed in an earlier work (Liu in Math Comput 70(234):579-593, 2000), to perform a high resolution numerical simulation of incompressible flow over an irregular domain and analyze its boundary layer separation. Compared with many classical finite element fluid solvers, this numerical method avoids a Stokes solver, and only two Poisson-like equations need to be solved at each time step/stage. In addition, its combination with the fully explicit fourth order Runge-Kutta (RK4) time discretization enables us to compute high Reynolds number flow in a very efficient way. As an application of this robust numerical solver, the dynamical mechanism of the boundary layer separation for a triangular cavity flow with Reynolds numbers \(Re=10^4\) and \(Re=10^5\) , including the precise values of bifurcation location and critical time, are reported in this paper. In addition, we provide a super-convergence analysis for the simple finite element numerical scheme, using linear elements over a uniform triangulation with right triangles.

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 Burton-Cabrera-Frank (BCF) model for the flow of line defects (steps) on crystal surfaces has offered useful insights into nanostructure evolution. This model has rested on phenomenological grounds. Our goal is to show via scaling arguments the emergence of the BCF theory for noninteracting steps from a stochastic atomistic scheme of a kinetic restricted solid-on-solid model in one spatial dimension. Our main assumptions are: adsorbed atoms (adatoms) form a dilute system, and elastic effects of the crystal lattice are absent. The step edge is treated as a front that propagates via probabilistic rules for atom attachment and detachment at the step. We formally derive a quasistatic step flow description by averaging out the stochastic scheme when terrace diffusion, adatom desorption, and deposition from above are present.

We provide a complete and rigorous description of phase transitions for kinetic models of self-propelled particles interacting through alignment. These models exhibit a competition between alignment and noise. Both the alignment frequency and noise intensity depend on a measure of the local alignment. We show that, in the spatially homogeneous case, the phase transition features (number and nature of equilibria, stability, convergence rate, phase diagram, hysteresis) are totally encoded in how the ratio between the alignment and noise intensities depend on the local alignment. In the spatially inhomogeneous case, we derive the macroscopic models associated to the stable equilibria and classify their hyperbolicity according to the same function.

Nonlinear hyperbolic systems with relaxations may encounter different scales of relaxation time, which is a prototype multiscale phenomenon that arises in many applications. In such a problem the relaxation time is of O(1) in part of the domain and very small in the remaining domain in which the solution can be approximated by the zero relaxation limit which can be solved numerically much more efficiently. For the Jin–Xin relaxation system in such a two-scale setting, we establish its wellposedness and singular limit as the (smaller) relaxation time goes to zero. The limit is a multiscale coupling problem which couples the original Jin–Xin system on the domain when the relaxation time is O(1) with its relaxation limit in the other domain through interface conditions which can be derived by matched interface layer analysis.As a result, we also establish the well-posedness and regularity (such as boundedness in sup norm with bounded total variation and L 1-contraction) of the coupling problem, thus providing a rigorous mathematical foundation, in the general nonlinear setting, to the multiscale domain decomposition method for this two-scale problem originally proposed in Jin et al. in Math. Comp. 82, 749–779, 2013.