A second-order projection method for the incompressible navier-stokes equations

Abstract

In this paper we describe a second-order projection method for the time-dependent, incompressible Navier-Stokes equations. As in the original projection method developed by Chorin, we first solve diffusion-convection equations to predict intermediate velocities which are then projected onto the space of divergence-free vector fields. By introducing more coupling between the diffusion--{;onvection step and the projection step we obtain a temporal discretization that is second-order accurate. Our treatment of the diffusion-convection step uses a specialized higher order Godunov method for differencing the nonlinear convective terms that provides a robust treatment of these terms at high Reynolds number. The Godunov procedure is second-order accurate for smooth flow and remains stable for discontinuous initial data, even in the zero-viscosity limit. We approximate the projection directly using a Galerkin procedure that uses a local basis for discretely divergence-free vector fields. Numerical results are presented validating the convergence properties of the method. We also apply the method to doubly periodic shear-layers to assess the performance of the method on more difficult applications

Computing multi-mode shock-induced compressible turbulent mixing at late times

“…High-resolution methods were originally designed to address issues of accuracy, and physically-correct behaviour in the proximity of discontinuities such as shock waves, as well as contact discontinuities. ILES is the established numerical approach for compressible turbulent mixing but is also widely used in many other Bell and Colella (1989) and Drikakis et al (1994) have shown that non-oscillatory methods can also be used for incompressible flows. In this study, the incompressible and compressible Euler equations are solved by nonoscilatory methods (Drikakis and Rider 2004).…”

Numerical Study on Mobilization of Oil Slugs in Capillary Model with Level set Approach

“…Typically, the projection method only requires the boundary condition for velocity, and no boundary conditions for pressure are required. Actually, specifying the pressure boundary condition will overestimate the system (Bell et al, 1989). The pressure intrinsically satisfies the Neumann condition, which reads…”

Hemodynamic Assessment of Pre- and Post-operative Aortic Coarctation from MRI

“…The Inlet and Wall cells are all interior to the domain, while the Outlet cells are situated on the domain boundaries, by extending the lumen of each vessel in its centerline direction [13] until it reaches the cartesian boundary. Our embedded boundary Navier-Stokes solver uses a fractional step method [14] that computes in a first step an intermediate velocity field, using the nonlinear advection-diffusion equation for velocity, and then projects the intermediate velocity onto the field of divergence free and tangent to the vessel boundary vector fields. For the velocity advection we use second-order upwind, Van-Leer slope limiting methods, while for the diffusion force components we use a semi-implicit approach as in [15] which is first order accurate and unconditionally stable in 3D.…”