No abstract uploaded!

Two-dimensional interaction between two Taylor vortices is simulated systematically through solving the two-dimensional, unsteady compressible Navier-Stokes equations using a fifth order weighted essentially nonoscillatory finite difference scheme. The main purpose of this study is to reveal the mechanism of sound generation in two-dimensional interaction of two Taylor vortices. Based on an extensive parameter study on the evolution of the vorticity field, we classify the interaction of two Taylor vortices into four types. The first type is the interaction of two counter-rotating vortices with similar strengths. The second type is the interaction of two co-rotating vortices without merging. The third type is the merging of two co-rotating vortices. The fourth type is the interaction of two vortices with a large difference in their strengths or scales. The mechanism of sound generation is analyzed.

No abstract uploaded!

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.

We present a new formulaiton of the incompressible Navier-Stokes equation in terms of an auxiliary field that differs from the velocity by a gauge transformation. The gauge freedom allows us to assign simple and specific boundary conditions for both the auxiliary field and the gauge field, thus eliminating the issue of pressure boundary condition in the usual primitive variable formulation. The resulting dynamic and kinematic equations can then be solved by standard methods for heat and Poisson equations. A normal mode analysis suggests that in contrast to the classical projection method, the gauge method does not suffer from the problem of numerical boundary layers. Thus the subtleties on the spatial discretization for the projection methods are removed. Consequently, the projection step in the gauge method can be accomplished by standard Poissin solves. We demonstrate the efficiency and accuracy of the gauge method by several numerical examples, including the flow past cylinder.