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!

No abstract uploaded!

We consider the Schrodinger equation with a logarithmic nonlinearity and a repulsive harmonic potential. Depending on the parameters of the equation, the solution may or may not be dispersive. When dispersion occurs, it does with an exponential rate in time. To control this, we change the unknown function through a generalized lens transform. This approach neutralizes the possible boundary effects, and could be used in the case of the Schrodinger equation without potential. We then employ standard splitting methods on the new equation via a nonuniform grid, after the logarithmic nonlinearity has been regularized. We also discuss the case of a power nonlinearity and give some results concerning the error estimates of the first-order Lie-Trotter splitting method for both cases of nonlinearities. Finally extensive numerical experiments are reported to investigate the dynamics of the equations.

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.