We consider an incompressible viscous flow without surface tension in a finite-depth domain of two dimensions, with free top boundary and fixed bottom boundary. This system is governed by the Navier-Stokes equations in this moving domain and the transport equation on the moving boundary. In this paper, we construct a stable numerical scheme to simulate the evolution of this system by discontinuous Galerkin method, and discuss the error analysis of the fluid under certain assumptions. Our formulation is mainly based on the geometric structure introduced earlier in the literature, and the natural energy estimate, which is rarely used in the numerical study of this system before.

In this paper, we introduce a general approach for proving optimal $L^2$ error estimates for the semi-discrete local discontinuous Galerkin (LDG) methods solving linear high order wave equations. The optimal order of error estimates hold not only for the solution itself but also for the auxiliary variables in the LDG method approximating the various order derivatives of the solution. Several examples including the one-dimensional third order wave equation, one-dimensional fifth order wave equation, and multi-dimensional Schr\"{o}dinger equation are explored to demonstrate this approach. The main idea is to derive energy stability for the various auxiliary variables in the LDG discretization, via using the scheme and its time derivatives with different test functions. Special projections are utilized to eliminate the jump terms at the cell boundaries in the error estimate in order to achieve the optimal order of accuracy.

Peskin's Immersed Boundary (IB) method is one of the most popular numerical methods for many years and has been applied to problems in mathematical biology, fluid mechanics, material sciences, and many other areas. Peskiness IB method is associated with discrete delta functions. It is believed that the IB method is first order accurate in the $L^{\infty}$ norm. But almost no rigorous proof could be found in the literature until recently [Mori, Comm. Pure. Appl. Math: 61:2008] in which the author showed that the velocity is indeed first order accurate for the Stokes equations with a periodic boundary condition. In this paper, we show first order convergence with a $\log h$ factor of the IB method for elliptic interface problems with Dirichlet boundary conditions.

METHODS AND APPLICATIONS OF ANALYSIS. 2000 International Press Vol. 7, No. 3, pp. 591-604, September 2000 012

The purpose of this work is to incorporate numerically, in a discontinuous Galerkin (DG) solver of a Boltzmann-Poisson model for hot electron transport, an electronic conduction band whose values are obtained by the spherical averaging of the full band structure given by a local empirical pseudopotential method (EPM) around a local minimum of the conduction band for silicon, as a midpoint between a radial band model and an anisotropic full band, in order to provide a more accurate physical description of the electron group velocity and conduction energy band structure in a semiconductor. This gives a better quantitative description of the transport and collision phenomena that fundamentally define the behaviour of the Boltzmann - Poisson model for electron transport used in this work. The numerical values of the derivatives of this conduction energy band, needed for the description of the electron group velocity, are obtained by means of a cubic spline interpolation. The EPM-Boltzmann-Poisson transport with this spherically averaged EPM calculated energy surface is numerically simulated and compared to the output of traditional analytic band models such as the parabolic and Kane bands, numerically implemented too, for the case of 1D $n^+-n-n^+$ silicon diodes with $400nm$ and $50nm$ channels. Quantitative differences are observed in the kinetic moments related to the conduction energy band used, such as mean velocity, average energy, and electric current (momentum).