In this article, we give a unified theory for constructing boundary layer expansions for discretized transport equations with homogeneous Dirichlet boundary conditions. We exhibit a natural assumption on the discretization under which the numerical solution can be written approximately as a two-scale boundary layer expansion. In particular, this expansion yields discrete semigroup estimates that are compatible with the continuous semigroup estimates in the limit where the space and time steps tend to zero. The novelty of our approach is to cover numerical schemes with arbitrarily many time levels.

In this paper, a diffuse-interface immersed boundary method (IBM) is proposed to treat three different thermal boundary conditions (Dirichlet, Neumann, Robin) in thermal flow problems. The novel IBM is implemented combining with the lattice Boltzmann method (LBM). The present algorithm enforces the three types of thermal boundary conditions at the boundary points. Concretely speaking, the IBM for the Dirichlet boundary condition is implemented using an iterative method, and its main feature is to accurately satisfy the given temperature on the boundary. The Neumann and Robin boundary conditions are implemented in IBM by distributing the jump of the heat flux on the boundary to surrounding Eulerian points, and the jump is obtained by applying the jump interface conditions in the normal and tangential directions. A simple analysis of the computational accuracy of IBM is developed. The analysis indicates that the Taylor-Green vortices problem which was used in many previous studies is not an appropriate accuracy test example. The capacity of the present thermal immersed boundary method is validated using four numerical experiments: (1) Natural convection in a cavity with a circular cylinder in the center; (2) Flows over a heated cylinder; (3) Natural convection in a concentric horizontal cylindrical annulus; (4) Sedimentation of a single isothermal cold particle in a vertical channel. The numerical results show good agreements with the data in the previous literatures.

We consider the obstacle problem for the Schrödinger evolution in the exterior of the unit ball with Dirichlet boundary condition. Under radial symmetry we compute explicitly the fundamental solution for the linear Dirichlet Schrödinger propagator ${{e}^{it{{\Delta }_{D}}}}$ and give a robust algorithm to prove sharp ${{L}^{1}}\,\to \,{{L}^{\infty }}$ dispersive estimates. We showcase the analysis in dimensions $n\,=\,5,\,7$. As an application, we obtain global well-posedness and scattering for defocusing energy-critical $\text{NLS}$ on $\Omega \,=\,{{\mathbb{R}}^{n}}\backslash \overline{B\left( 0,\,1 \right)}$ with Dirichlet boundary condition and radial data in these dimensions.

In this paper, we study numerically quantized vortex dynamics and their interaction in the two-dimensional (2D) Ginzburg-Landau equation (GLE) with a dimensionless parameter ε > 0 on bounded domains under either Dirichlet or homogeneous Neumann boundary condition. We begin with a review of the reduced dynamical laws for time evolution of quantized vortex centers in GLE and show how to solve these nonlinear ordinary differential equations numerically. Then we present efficient and accurate numerical methods for discretizing the GLE on either a rectangular or a disk domain under either Dirichlet or homogeneous Neumann boundary condition. Based on these efficient and accurate numerical methods for GLE and the reduced dynamical laws, we simulate quantized vortex interaction of GLE with different ε and under different initial setups including single vortex, vortex pair, vortex dipole and vortex lattice, compare them with those obtained from the corresponding reduced dynamical laws, and identify the cases where the reduced dynamical laws agree qualitatively and/or quantitatively as well as fail to agree with those from GLE on vortex interaction. Finally, we also obtain numerically different patterns of the steady states for quantized vortex lattices under the GLE dynamics on bounded domains.