This paper introduces a unified concept and algorithm for the fractional-step (FS), artificial compressibility (AC) and pressure-projection (PP) methods for solving the incompressible Navier-Stokes equations. The proposed FSAC-PP approach falls into the group of pseudo-time splitting high-resolution methods incorporating the characteristics-based (CB) Godunov-type treatment of convective terms with PP methods. Due to the fact that the CB Godunov-type methods are applicable directly to the hyperbolic AC formulation and not to the elliptical FS-PP (split) methods, thus the straightforward coupling of CB Godunov-type schemes with PP methods is not possible. Therefore, the proposed FSAC-PP approach unifies the fully-explicit AC and semi-implicit FS-PP methods of Chorin including a PP step in the dual-time stepping procedure to a) overcome the numerical stiffness of the classical AC approach at (very) low and moderate Reynolds numbers, b) incorporate the accuracy and convergence properties of CB Godunov-type schemes with PP methods, and c) further improve the stability and efficiency of the AC method for steady and unsteady flow problems. The FSAC-PPmethod has also been coupled with a non-linear, full-multigrid and fullapproximation storage (FMG-FAS) technique to further increase the efficiency of the solution. For validating the proposed FSAC-PP method, computational examples are presented for benchmark problems. The overall results show that the unified FSAC-PP approach is an efficient algorithm for solving incompressible flow problems.

We use the heat equation as an illustrative example to show that the unified method introduced by one of the authors can be employed for constructing analytical solutions for linear evolution partial differential equations in one spatial dimension involving non-separable boundary conditions as well as non-local constraints. Furthermore, we show that for the particular case in which the boundary conditions become separable, the unified method provides an easier way for constructing the relevant classical spectral representations avoiding the classical spectral analysis approach. We note that the unified method always yields integral expressions which, in contrast to the series or integral expressions obtained by the standard transform methods, are uniformly convergent at the boundary. Thus, even for the cases that the standard transform methods can be implemented, the unified method provides alternative solution expressions which have advantages for both numerical and asymptotic considerations. The former advantage is illustrated by providing the numerical evaluation of typical boundary value problems.