{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,1,13]],"date-time":"2026-01-13T04:05:15Z","timestamp":1768277115885,"version":"3.49.0"},"reference-count":40,"publisher":"MDPI AG","issue":"6","license":[{"start":{"date-parts":[[2023,6,2]],"date-time":"2023-06-02T00:00:00Z","timestamp":1685664000000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Algorithms"],"abstract":"<jats:p>To enforce the conservation of mass principle, a pressure Poisson equation arises in the numerical solution of incompressible fluid flow using the pressure-based segregated algorithms such as projection methods. For unsteady flows, the pressure Poisson equation is solved at each time step usually in physical space using iterative solvers, and the resulting pressure gradient is then applied to make the velocity field divergence-free. It is generally accepted that this pressure-correction stage is the most time-consuming part of the flow solver and any meaningful acceleration would contribute significantly to the overall computational efficiency. The objective of the present work was to develop a fast hybrid pressure-correction algorithm for numerical simulation of incompressible flows around obstacles in the context of projection methods. The key idea is to adopt different numerical methods\/discretisations in the sub-steps of projection methods. Here, a classical second-order time-marching projection method, which consists of two sub-steps, was chosen for the purposes of demonstration. In the first sub-step, the momentum equations were discretised on unstructured grids and solved by conventional numerical methods, here a meshless method. In the second sub-step (pressure-correction), the proposed algorithm adopts a double-discretisation system and combines the weighted least-squares approximation with the essence of immersed boundary methods. Such a design allowed us to develop an FFT-based solver to speed up the solution of the pressure Poisson equation for flow cases with obstacles, while keeping the implementation of the boundary conditions for the momentum equations as easy as conventional numerical methods do with unstructured grids. The numerical experiments of five test cases were performed to verify and validate the proposed hybrid algorithm and evaluate its computational performance. The results showed that the new FFT-based hybrid algorithm works and is robust, and it was significantly faster than the multigrid-based reference method. The hybrid algorithm opens an avenue for the development of next-generation high-performance parallel computational fluid dynamics solvers for incompressible flows.<\/jats:p>","DOI":"10.3390\/a16060287","type":"journal-article","created":{"date-parts":[[2023,6,2]],"date-time":"2023-06-02T01:33:54Z","timestamp":1685669634000},"page":"287","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":3,"title":["A Fast Hybrid Pressure-Correction Algorithm for Simulating Incompressible Flows by Projection Methods"],"prefix":"10.3390","volume":"16","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-9205-6444","authenticated-orcid":false,"given":"Jiannong","family":"Fang","sequence":"first","affiliation":[{"name":"Wind Engineering and Renewable Energy Laboratory (WiRE), School of Architecture, Civil and Environmental Engineering (ENAC), \u00c9cole Polytechnique F\u00e9d\u00e9rale de Lausanne (EPFL), 1015 Lausanne, Switzerland"},{"name":"Scientific IT and Application Support (SCITAS), \u00c9cole Polytechnique F\u00e9d\u00e9rale de Lausanne (EPFL), 1015 Lausanne, Switzerland"}]}],"member":"1968","published-online":{"date-parts":[[2023,6,2]]},"reference":[{"key":"ref_1","doi-asserted-by":"crossref","first-page":"745","DOI":"10.1090\/S0025-5718-1968-0242392-2","article-title":"Numerical solution of the navier-stokes equations","volume":"22","author":"Chorin","year":"1968","journal-title":"Math. Comput."},{"key":"ref_2","doi-asserted-by":"crossref","first-page":"115","DOI":"10.24033\/bsmf.1662","article-title":"Une m\u00e9thode d\u2019approximation de la solution des \u00e9quations de Navier\u2013Stokes","volume":"98","author":"Temam","year":"1968","journal-title":"Bull. Soc. Math. France"},{"key":"ref_3","doi-asserted-by":"crossref","first-page":"6011","DOI":"10.1016\/j.cma.2005.10.010","article-title":"An overview of projection methods for incompressible flows","volume":"195","author":"Guermond","year":"2006","journal-title":"Comput. Methods Appl. Mech. Eng."},{"key":"ref_4","doi-asserted-by":"crossref","first-page":"257","DOI":"10.1016\/0021-9991(89)90151-4","article-title":"A second-order projection method for the incompressible navier-stokes equations","volume":"85","author":"Bell","year":"1989","journal-title":"J. Comput. Phys."},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"464","DOI":"10.1006\/jcph.2001.6715","article-title":"Accurate Projection Methods for the Incompressible Navier\u2013Stokes Equations","volume":"168","author":"Brown","year":"2001","journal-title":"J. Comput. Phys."},{"key":"ref_6","doi-asserted-by":"crossref","first-page":"308","DOI":"10.1016\/0021-9991(85)90148-2","article-title":"Application of a fractional-step method to incompressible Navier\u2013Stokes equations","volume":"59","author":"Kim","year":"1985","journal-title":"J. Comput. Phys."},{"key":"ref_7","doi-asserted-by":"crossref","first-page":"2121","DOI":"10.1029\/1999WR900094","article-title":"Surface length scales and shear stress: Implications for land-atmosphere interaction over complex terrain","volume":"35","author":"Albertson","year":"1999","journal-title":"Water Resour. Res."},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"2178","DOI":"10.2514\/1.J054569","article-title":"Approximate projection method for the incompressible navier-stokes equations","volume":"54","author":"Capuano","year":"2016","journal-title":"AIAA J."},{"key":"ref_9","doi-asserted-by":"crossref","first-page":"325","DOI":"10.1016\/j.jcp.2004.04.006","article-title":"A class of fully second order accurate projection methods for solving the incompressible Navier\u2013Stokes equations","volume":"200","author":"Liu","year":"2004","journal-title":"J. Comput. Phys."},{"key":"ref_10","doi-asserted-by":"crossref","first-page":"5989","DOI":"10.1016\/j.jcp.2009.05.010","article-title":"High-order compact schemes for incompressible flows: A simple and efficient method with quasi-spectral accuracy","volume":"228","author":"Laizet","year":"2009","journal-title":"J. Comput. Phys."},{"key":"ref_11","doi-asserted-by":"crossref","first-page":"1735","DOI":"10.1002\/fld.2480","article-title":"Incompact3d: A powerful tool to tackle turbulence problems with up to O(105) computational cores","volume":"67","author":"Laizet","year":"2011","journal-title":"Int. J. Numer. Methods Fluids"},{"key":"ref_12","doi-asserted-by":"crossref","first-page":"A20","DOI":"10.1017\/jfm.2019.869","article-title":"Non-equilibrium three-dimensional boundary layers at moderate Reynolds numbers","volume":"883","author":"Duran","year":"2019","journal-title":"J. Fluid Mech."},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"273","DOI":"10.1016\/j.nucengdes.2016.12.018","article-title":"Direct numerical simulation of a low momentum round jet in channel crossflow","volume":"313","author":"Wu","year":"2017","journal-title":"Nucl. Eng. Des."},{"key":"ref_14","doi-asserted-by":"crossref","first-page":"2052","DOI":"10.1175\/1520-0469(1984)041<2052:ALESMF>2.0.CO;2","article-title":"A large- eddy-simulation model for the study of planetary boundary-layer turbulence","volume":"41","author":"Moeng","year":"1984","journal-title":"J. Atmos. Sci."},{"key":"ref_15","doi-asserted-by":"crossref","first-page":"125","DOI":"10.1007\/s10546-022-00749-y","article-title":"Scaling of flows over realistic urban geometries: A large-eddy simulation study","volume":"186","author":"Cheng","year":"2023","journal-title":"Bound.-Layer Meteorol."},{"key":"ref_16","doi-asserted-by":"crossref","first-page":"584","DOI":"10.1006\/jcph.1999.6246","article-title":"An SPH Projection Method","volume":"152","author":"Cummins","year":"1999","journal-title":"J. Comput. Phys."},{"key":"ref_17","doi-asserted-by":"crossref","first-page":"275","DOI":"10.1016\/S0141-1187(03)00002-6","article-title":"Simulation of near-shore solitary wave mechanics by an incompressible SPH method","volume":"24","author":"Lo","year":"2002","journal-title":"Appl. Ocean Res."},{"key":"ref_18","doi-asserted-by":"crossref","first-page":"15","DOI":"10.1007\/s11804-022-00260-3","article-title":"A Coupled ISPH-TLSPH Method for Simulating Fluid-Elastic Structure Interaction Problems","volume":"21","author":"Salehizadeh","year":"2022","journal-title":"J. Mar. Sci. Appl."},{"key":"ref_19","doi-asserted-by":"crossref","first-page":"369","DOI":"10.1016\/0021-9991(91)90215-7","article-title":"An improvement of fractional step methods for the incompressible Navier\u2013Stokes equations","volume":"92","author":"Le","year":"1991","journal-title":"J. Comput. Phys."},{"key":"ref_20","doi-asserted-by":"crossref","unstructured":"de Michele, C., Capuano, F., and Coppola, G. (2020). Fast-Projection Methods for the Incompressible Navier\u2013Stokes Equations. Fluids, 5.","DOI":"10.3390\/fluids5040222"},{"key":"ref_21","doi-asserted-by":"crossref","first-page":"109693","DOI":"10.1016\/j.jcp.2020.109693","article-title":"A fast pressure-correction method for incompressible flows over curved walls","volume":"421","author":"Aithal","year":"2020","journal-title":"J. Comput. Phys."},{"key":"ref_22","doi-asserted-by":"crossref","first-page":"110518","DOI":"10.1016\/j.jcp.2021.110518","article-title":"Low-cost Runge-Kutta integrators for incompressible flow simulations","volume":"443","author":"Karam","year":"2021","journal-title":"J. Comput. Phys."},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"416","DOI":"10.1016\/j.jcp.2014.05.024","article-title":"A fast pressure-correction method for incompressible two-fluid flows","volume":"273","author":"Dodd","year":"2014","journal-title":"J. Comput. Phys."},{"key":"ref_24","doi-asserted-by":"crossref","first-page":"108028","DOI":"10.1016\/j.cpc.2021.108028","article-title":"Towards pseudo-spectral incompressible smoothed particle hydrodynamics (ISPH)","volume":"266","author":"Fourtakas","year":"2021","journal-title":"Comput. Phys. Commun."},{"key":"ref_25","doi-asserted-by":"crossref","first-page":"239","DOI":"10.1146\/annurev.fluid.37.061903.175743","article-title":"Immersed boundary methods","volume":"37","author":"Mittal","year":"2005","journal-title":"Annu. Rev. Fluid Mech."},{"key":"ref_26","doi-asserted-by":"crossref","first-page":"427","DOI":"10.1016\/j.jcp.2006.12.009","article-title":"Modeling turbulent flow over fractal trees with renormalized numerical simulation","volume":"225","author":"Chester","year":"2007","journal-title":"J. Comput. Phys."},{"key":"ref_27","doi-asserted-by":"crossref","first-page":"109971","DOI":"10.1016\/j.buildenv.2022.109971","article-title":"Effects of unstable thermal stratification on the flow characteristics in an idealized rural-to-urban transition region: A large-eddy simulation study","volume":"230","author":"Ma","year":"2023","journal-title":"Build. Environ."},{"key":"ref_28","doi-asserted-by":"crossref","first-page":"421","DOI":"10.1007\/s10546-017-0283-9","article-title":"Large-Eddy Simulations of Atmospheric Flows Over Complex Terrain Using the Immersed-Boundary Method in the Weather Research and Forecasting Model","volume":"165","author":"Ma","year":"2017","journal-title":"Bound.-Layer Meteorol."},{"key":"ref_29","doi-asserted-by":"crossref","first-page":"104604","DOI":"10.1016\/j.compfluid.2020.104604","article-title":"Wall modeled immersed boundary method for high Reynolds number flow over complex terrain","volume":"208","author":"Liu","year":"2020","journal-title":"Comput. Fluids"},{"key":"ref_30","doi-asserted-by":"crossref","first-page":"082008","DOI":"10.1088\/1742-6596\/753\/8\/082008","article-title":"Intercomparison of terrain-following coordinate transformation and immersed boundary methods in large-eddy simulation of wind fields over complex terrain","volume":"753","author":"Fang","year":"2016","journal-title":"J. Phys. Conf. Ser."},{"key":"ref_31","doi-asserted-by":"crossref","first-page":"8894","DOI":"10.1016\/j.jcp.2008.06.031","article-title":"A regularized Lagrangian finite point method for the simulation of incompressible viscous flows","volume":"227","author":"Fang","year":"2008","journal-title":"J. Comput. Phys."},{"key":"ref_32","unstructured":"Hockney, R.W., and Eastwood, J.W. (1981). Computer Simulation Using Particles, McGraw-Hill."},{"key":"ref_33","unstructured":"Press, W.H., Teukolsky, S.A., Vetterling, W.T., and Flannery, B.P. (1992). Numerical Recipes in FORTRAN: The Art of Scientific Computing, Cambridge University Press. [2nd ed.]."},{"key":"ref_34","first-page":"123","article-title":"An aggregation-based algebraic multigrid method","volume":"37","author":"Notay","year":"2010","journal-title":"Electron. Trans. Numer. Anal."},{"key":"ref_35","doi-asserted-by":"crossref","first-page":"A1079","DOI":"10.1137\/100818509","article-title":"An algebraic multigrid method with guaranteed convergence rate","volume":"34","author":"Napov","year":"2012","journal-title":"SIAM J. Sci. Comput."},{"key":"ref_36","doi-asserted-by":"crossref","first-page":"A2288","DOI":"10.1137\/110835347","article-title":"Aggregation-based algebraic multigrid for convection-diffusion equations","volume":"34","author":"Notay","year":"2012","journal-title":"SIAM J. Sci. Comput."},{"key":"ref_37","doi-asserted-by":"crossref","first-page":"8179","DOI":"10.1016\/j.jcp.2011.07.017","article-title":"Towards oscillation-free implementation of the immersed boundary method with spectral-like methods","volume":"230","author":"Fang","year":"2011","journal-title":"J. Comput. Phys."},{"key":"ref_38","doi-asserted-by":"crossref","first-page":"356","DOI":"10.1016\/j.amc.2015.03.011","article-title":"PoisFFT\u2014A free parallel fast Poisson solver","volume":"267","author":"Fuka","year":"2015","journal-title":"Appl. Math. Comput."},{"key":"ref_39","doi-asserted-by":"crossref","first-page":"C280","DOI":"10.1137\/15M1010798","article-title":"FFT, FMM, or Multigrid? A comparative Study of State-Of-the-Art Poisson Solvers for Uniform and Nonuniform Grids in the Unit Cube","volume":"38","author":"Gholami","year":"2016","journal-title":"SIAM J. Sci. Comput."},{"key":"ref_40","doi-asserted-by":"crossref","first-page":"311","DOI":"10.1002\/fld.1650050403","article-title":"Numerical and experimental study of driven flow in a polar cavity","volume":"5","author":"Fuchs","year":"1985","journal-title":"Int. J. Numer. Methods Fluids"}],"container-title":["Algorithms"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/1999-4893\/16\/6\/287\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,10,10]],"date-time":"2025-10-10T19:47:33Z","timestamp":1760125653000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/1999-4893\/16\/6\/287"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2023,6,2]]},"references-count":40,"journal-issue":{"issue":"6","published-online":{"date-parts":[[2023,6]]}},"alternative-id":["a16060287"],"URL":"https:\/\/doi.org\/10.3390\/a16060287","relation":{},"ISSN":["1999-4893"],"issn-type":[{"value":"1999-4893","type":"electronic"}],"subject":[],"published":{"date-parts":[[2023,6,2]]}}}