MathSciDoc: An Archive for Mathematician ∫

We investigate, by numerical simulation, the shear layer instability associated with the outer layer of a spiral vortex formed behind an impulsively started thin ellipse. The unstable free shear layer undergoes a secondary instability. We connect this instability with the dynamics of corner vortices adjacent to the tip of the ellipse by observing that the typical turnover time of the corner vortex matches the period of the unstable mode in the shear layer. We suggest that the corner vortex acts as a signal generator, and produces periodic perturbation which triggers the instability.

Solutions to conservation laws satisfy the monotonicity property: the number of local extrema is a non-increasing function of time, and local maximum/minimum values decrease/increase monotonically in time. This paper investigates this property from a numerical standpoint. We introduce a class of fully discrete in space and time, high order accurate, difference schemes, called generalized monotone schemes. Convergence toward the entropy solution is proven via a new technique of proof, assuming that the initial data has a finite number of extremum values only, and the flux-function is strictly convex. We define discrete paths of extrema by tracking local extremum values in the approximate solution. In the course of the analysis we establish the pointwise convergence of the trace of the solution along a path of extremum. As a corollary, we obtain a proof of convergence for a MUSCL-type scheme being second order accurate away from sonic points and extrema.

The Euler equation of compressible flows is solved by the finite volume method, where high order accuracy is achieved by the reconstruction of each component of upwind fluxes of a flux splitting using the biased averaging procedure. Compared to the solution reconstruction in Godunov-type methods, its implementation is simple and easy, and the computational complexity is relatively low. This approach is parameter-free and requires neither a Riemann solver nor field-by-field decomposition. The numerical results from both dynamic and steady state calculations demonstrate the accuracy and robustness of this approach. Some techniques for the acceleration of the convergence to the steady state are discussed, including multigrid and multistage Runge–Kutta time methods.

In this ar ti cle, the authors de vel oped an al ter na tive meth od ol ogy to cal cu late theoption- adjusted spread for the mortgage- backed se cu ri ties us ing par tial dif fer ential equa tion tech nique. The nu meri cal im ple men ta tion is dis cussed in de tail, includ ing the con ver gence and er ror analy sis. This ap proach pro vides us a fast andac cu rate way to pric ing MBS.

We consider finite difference schemes based on the impulse density variable. We show that the original velocity—impulse density formulation of Oseledets is marginally ill-posed for the inviscid flow, and this has the consequence that some ordinarily stable numerical methods in other formulations become unstable in the velocity—impulse density formulation. We present numerical evidence of this instability. We then discuss the construction of stable finite difference schemes by requiring that at the numerical level the nonlinear terms be convertible to similar terms in the primitive variable formulation. Finally we give a simplified velocity—impulse density formulation which is free of these complications and yet retains the nice features of the original velocity—impulse density formulation with regard to the treatment of boundary. We present numerical results on this simplified formulation for the driven cavity flow on both the staggered and non-staggered grids.

We are concerned with the convergence of Lax-Wendroff type schemes with high resolution to the entropy solutions for conservation laws. These schemes include the original Lax-Wendroff scheme proposed by Lax and Wendroff in 1960 and its two step versions–the Richtmyer scheme and the MacCormack scheme. For the convex scalar conservation laws with algebraic growth flux functions, we prove the convergence of these schemes to the weak solutions satisfying appropriate entropy inequalities. The proof is based on detailed L p estimates of the approximate solutions, H −1 compactness estimates of the corresponding entropy dissipation measures, and some compensated compactness frameworks. Then these techniques are generalized to study the convergence problem for the nonconvex scalar case and the hyperbolic systems of conservation

Simple, efficient, and accurate finite difference methods are introduced for 3D unsteady viscous incompressible flows in the vorticity–vector potential formulation on nonstaggered grids. Two different types of methods are discussed. They differ in the implementation of the normal component of the vorticity boundary condition and consequently the enforcement of the divergence free condition for vorticity. Both second-order and fourth-order accurate schemes are developed. A detailed accuracy test is performed, revealing the structure of the error and the effect of how the convective terms are discretized near the boundary. The influence of the divergence free condition for vorticity to the overall accuracy is studied. Results on the cubic driven cavity flow at Reynolds number 500 and 3200 are shown and compared with that of the MAC scheme.

In this paper, we study the fluid-dynamic limit for the one-dimensional Broadwell model of the nonlinear Boltzmann equation in the presence of boundaries. We consider an analogue of Maxwell's diffusive and reflective boundary conditions. The boundary layers can be classified as either compressive or expansive in terms of the associated characteristic fields. We show that both expansive and compressive boundary layers (before detachment) are nonlinearly stable and that the layer effects are localized so that the fluid dynamic approximation is valid away from the boundary. We also show that the same conclusion holds for short time without the structural conditions on the boundary layers. A rigorous estimate for the distance between the kinetic solution and the fluid-dynamic solution in terms of the mean-free path in theL∞-norm is obtained provided that the interior fluid flow is smooth. The rate of convergence is optimal.

In this paper, we investigate the boundary layer behavior of solutions to the one dimensional Broadwell model of the nonlinear Boltzmann equation for small mean free path. We consider the analogue of Maxwell's diffusive and the reflexive boundary conditions. It is found that even for such a simple model, there are boundary layers due to purely kinetic effects which cannot be detected by the corresponding Navier-Stokes system. It is also found numerically that a compressive boundary layer is not always stable in the sense that it may detach from the boundary and move into the interior of the gas as a shock layer.

We study the oscillatory behavior that arises in solutions of a dispersive numerical scheme for the Hopf equation whenever the classical solution of the equation develops a singularity. Modulation equations are derived that describe period-two oscillations so long as the solution of those equations takes values for which the equations are hyperbolic. However, those equations have an elliptic region that may be entered by its solutions in a finite time, after which the corresponding period-two oscillations are seen to break down. This kind of phenomenon has not been observed for integrable schemes. The generation and propagation of period-two oscillations are asymptotically analyzed and a matching formula is found for the transition between oscillatory and nonoscillatory regions. Modulation equations are also presented for period-three oscillations. Numerical experiments are carried out that illustrate our analysis.