No abstract uploaded!

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.

Adapting the convex integration technique introduced in \cite{dLSz5} and subsequently developed in \cite{IsettVicol,IsettOh16,Buckmaster2013transporting}, we construct H\"{o}lder continuous weak solutions to the three dimensional Prandtl system and some other models with vertical viscosity.

In this work, we consider the numerical solution of the nonlinear Schrödinger equation with a highly oscillatory potential (NLSE-OP). The NLSE-OP is a model problem which frequently occurs in recent studies of some multiscale dynamical systems, where the potential introduces wide temporal oscillations to the solution and causes numerical difficulties. We aim to analyze rigorously the error bounds of the splitting schemes for solving the NLSE-OP to a fixed time. Our theoretical results show that the Lie–Trotter splitting scheme is uniformly and optimally accurate at the first order provided that the oscillatory potential is integrated exactly, while the Strang splitting scheme is not. Our results apply to general dispersive or wave equations with an oscillatory potential. The error estimates are confirmed by numerical results.

The Boltzmann equation which describes the time evolution of a large number of particles through the binary collision in statistics physics has close relation to the systems of fluid dynamics, that is, Euler equations and NavierStokes equations. As for a basic wave pattern to Euler equations, we consider the nonlinear stability of contact discontinuities to the Boltzmann equation. Even though the stability of the other two nonlinear waves, i.e., shocks and rarefaction waves has been extensively studied, there are few stability results on the contact discontinuity because unlike shock waves and rarefaction waves, its derivative has no definite sign, and decays slower than a rarefaction wave. Moreover, it behaves like a linear wave in a nonlinear setting so that its coupling with other nonlinear waves reveals a complicated interaction mechanism. Based on the new definition of contact waves to the Boltzmann equation