No abstract uploaded!

Let $K$ be a self-similar set satisfying the open set condition. Following Kaimanovich’s elegant idea, it has been proved that on the symbolic space $X$ of $K$ a natural augmented tree structure $E$ exists; it is hyperbolic, and the hyperbolic boundary $\partial_H X$ with the Gromov metric is H\"older equivalent to $K$. In this paper we consider certain reversible random walks with return ratio $0 < \lambda < 1$ on $(X,E)$. We show that the Martin boundary ${\mathcal M}$ can be identified with $\partial_H X$ and $K$. With this setup and a device of Silverstein, we obtain precise estimates of the Martin kernel and the Na\"im kernel in terms of the Gromov product. Moreover, the Naïm kernel turns out to be a jump kernel satisfying the estimate $\Theta (\xi,\eta) \asymp |\xi-\eta|^{-(\alpha+\beta)}$, where $\alpha$ is the Hausdorff dimension of $K$ and $\beta$ depends on $\lambda$. For suitable $\beta$, the kernel defines a regular non-local Dirichlet form on $K$. This extends the results of Kigami concerning random walks on certain trees with Cantor-type sets as boundaries.

The immersed boundary method is one of the most useful computational methods in studying fluid structure interaction. On the other hand, the Immersed Boundary method is also known to require small time steps to maintain stability when solved with an explicit method. Many implicit or approximately implicit methods have been proposed in the literature to remove this severe time step stability constraint, but none of them give satisfactory performance. In this paper, we propose an efficient semi-implicit scheme to remove this stiffness from the immersed boundary method for the Navier–Stokes equations. The construction of our semi-implicit scheme consists of two steps. First, we obtain a semiimplicit discretization which is proved to be unconditionally stable. This unconditionally stable semi-implicit scheme is still quite expensive to implement in practice. Next, we apply the small scale decomposition to the unconditionally stable semi-implicit scheme to construct our efficient semi-implicit scheme. Unlike other implicit or semi-implicit schemes proposed in the literature, our semi-implicit scheme can be solved explicitly in the spectral space. Thus the computational cost of our semi-implicit schemes is comparable to that of an explicit scheme. Our extensive numerical experiments show that our semiimplicit scheme has much better stability property than an explicit scheme. This offers a substantial computational saving in using the immersed boundary method.

For a bipolar hydrodynamic model of semiconductors in the form of EulerPoisson equations with Dirichlet or Neumann boundary conditions, in this paper we first heuristically analyze the most probable asymptotic profile (the so-called diffusion waves) and then prove this long-time behavior rigorously. For this, we construct correction functions to show the convergence of the original solution to the diffusion wave with optimal convergence rates by the energy method. Moreover, in the case with Dirichlet boundary condition, when the initial perturbation is in some weighted L^1 space, a faster and optimal convergence rate is also given.

In earlier work, Maire developed a class of cell-centered Lagrangian schemes for solving Euler equations of compressible gas dynamics in cylindrical coordinates. These schemes use a node-based discretization of the numerical fluxes. The control volume version has several distinguished properties, including the conservation of mass, momentum and total energy and compatibility with the geometric conservation law (GCL). However it also has a limitation in that it cannot preserve spherical symmetry for one-dimensional spherical flow. An alternative is also given to use the first order area-weighted approach which can ensure spherical symmetry, at the price of sacrificing conservation of momentum. In this paper, we apply the methodology proposed in our recent work to the first order control volume scheme of Maire to obtain the spherical symmetry property. The modified scheme can preserve one-dimensional spherical symmetry in a two-dimensional cylindrical geometry when computed on an equal-angle-zoned initial grid, and meanwhile it maintains its original good properties such as conservation and GCL. Several two-dimensional numerical examples in cylindrical coordinates are presented to demonstrate the good performance of the scheme in terms of symmetry, non-oscillation and robustness properties.