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.

A two species interacting system motivated by the density functional theory for triblock copolymers contains long range interaction that affects the two species differently. In a two species periodic assembly of discs, the two species appear alternately on a lattice. The minimal two species periodic assembly is the one with the least energy per lattice cell area. There is a parameter $b$ in $[0,1]$ and the type of the lattice associated with the minimal assembly varies depending on $b$. There are several threshold defined by a number $B=0.1867...$ If $b \in [0, B)$, the minimal assembly is associated with a rectangular lattice whose ratio of the longer side and the shorter side is in $[\sqrt{3}, 1)$; if $b \in [B, 1-B]$, the minimal assembly is associated with a square lattice; if $b \in (1-B, 1]$, the minimal assembly is associated with a rhombic lattice with an acute angle in $[\frac{\pi}{3}, \frac{\pi}{2})$. Only when $b=1$, this rhombic lattice is the hexagonal lattice. None of the other values of $b$ yields the hexagonal lattice, a sharp contrast to the situation for one species interacting systems, where the hexagonal lattice is ubiquitously observed.

In this article, we continue the work in Guan-Li and study a normalized hypersurface flow in the more general ambient setting of warped product spaces. This flow preserves the volume of the bounded domain enclosed by a graphical hypersurface, and monotonically decreases the hypersurface area. As an application, the isoperimetric problem in warped product spaces is solved for such domains.

Given a compact Riemannian manifold $M$, we consider a warped product manifold $\bar M = I \times_h M$, where $I$ is an open interval in $\Bbb R$. For a positive function $\psi$ defined on $\bar M$, we generalize the arguments in \cite{GRW2015} and \cite{RW16}, to obtain the curvature estimates for Hessian equations $\sigma_k(\kappa)=\psi(V,\nu(V))$. We also obtain some existence results for the starshaped compact hypersurface $\Sigma$ satisfying the above equation with various assumptions.

An analogue of the well-known $$ \frac{3}{{16}} $$ lower bound for the first eigenvalue of the Laplacian for a congruence hyperbolic surface is proven for a congruence tower associated with any non-elementary subgroup$L$of SL(2,$Z$). The proof in the case that the Hausdorff of the limit set of$L$is bigger than $$ \frac{1}{2} $$ is based on a general result which allows one to transfer such bounds from a combinatorial version to this archimedian setting. In the case that delta is less than $$ \frac{1}{2} $$ we formulate and prove a somewhat weaker version of this phenomenon in terms of poles of the corresponding dynamical zeta function. These “spectral gaps” are then applied to sieving problems on orbits of such groups.