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.

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.

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.

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.

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.