{"title":"A pressure-projection pre-conditioning multi-fractional-step method for Navier–Stokes Flow in Porous Media","authors":"Lateef T. Akanji","doi":"10.1016/j.jocs.2024.102313","DOIUrl":null,"url":null,"abstract":"<div><p>A new pressure-projection pre-conditioning multi-fractional-step (PPP-MFS) method for incompressible Navier–Stokes flow in porous media is presented. This fractional step method is applied to decouple the pressure and the velocity; thereby, overcoming the computational costs and difficulty associated with the resolution of the nonlinear term in the Navier–Stokes equation for fine geometric models. Specifically, time evolution is decomposed into a sequence of multi-fractional solution steps. In the first step, an elliptic problem is solved for pressure (<span><math><mi>p</mi></math></span>) with a no-slip boundary condition. This gives the Stokes pressure and velocity fields. In the second step, the obtained pressure <span><math><mi>p</mi></math></span> is then projected onto the field <span><math><mrow><mi>p</mi><mo>∗</mo></mrow></math></span> and used to solve for velocity field (<span><math><mi>u</mi></math></span>) required for the pre-conditioning of the solution to the Navier–Stokes equation. The pressure and velocity fields are obtained from the solution of the Navier–Stokes equation in the third step. Numerical and geometric discretisation of porous samples were carried out using finite-element method. For flow in simple channel models represented by two- and three-dimensions and in systems with high conductivity, the Stokes and Navier–Stokes numerical solutions produced close pressure and velocity field approximations. For flow around a cylinder, computation time is consistently higher in the Navier–Stokes equation by a factor of 2 with a pronounced non-symmetric pressure field at high mesh refinements. This computation time is desirable given that Navier–Stokes computation without preconditioning can be orders of magnitude more expensive.</p></div>","PeriodicalId":48907,"journal":{"name":"Journal of Computational Science","volume":null,"pages":null},"PeriodicalIF":3.1000,"publicationDate":"2024-05-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S1877750324001066/pdfft?md5=5f19a37cb40cd73380263125dfe6d958&pid=1-s2.0-S1877750324001066-main.pdf","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational Science","FirstCategoryId":"94","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S1877750324001066","RegionNum":3,"RegionCategory":"计算机科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
A new pressure-projection pre-conditioning multi-fractional-step (PPP-MFS) method for incompressible Navier–Stokes flow in porous media is presented. This fractional step method is applied to decouple the pressure and the velocity; thereby, overcoming the computational costs and difficulty associated with the resolution of the nonlinear term in the Navier–Stokes equation for fine geometric models. Specifically, time evolution is decomposed into a sequence of multi-fractional solution steps. In the first step, an elliptic problem is solved for pressure () with a no-slip boundary condition. This gives the Stokes pressure and velocity fields. In the second step, the obtained pressure is then projected onto the field and used to solve for velocity field () required for the pre-conditioning of the solution to the Navier–Stokes equation. The pressure and velocity fields are obtained from the solution of the Navier–Stokes equation in the third step. Numerical and geometric discretisation of porous samples were carried out using finite-element method. For flow in simple channel models represented by two- and three-dimensions and in systems with high conductivity, the Stokes and Navier–Stokes numerical solutions produced close pressure and velocity field approximations. For flow around a cylinder, computation time is consistently higher in the Navier–Stokes equation by a factor of 2 with a pronounced non-symmetric pressure field at high mesh refinements. This computation time is desirable given that Navier–Stokes computation without preconditioning can be orders of magnitude more expensive.
期刊介绍:
Computational Science is a rapidly growing multi- and interdisciplinary field that uses advanced computing and data analysis to understand and solve complex problems. It has reached a level of predictive capability that now firmly complements the traditional pillars of experimentation and theory.
The recent advances in experimental techniques such as detectors, on-line sensor networks and high-resolution imaging techniques, have opened up new windows into physical and biological processes at many levels of detail. The resulting data explosion allows for detailed data driven modeling and simulation.
This new discipline in science combines computational thinking, modern computational methods, devices and collateral technologies to address problems far beyond the scope of traditional numerical methods.
Computational science typically unifies three distinct elements:
• Modeling, Algorithms and Simulations (e.g. numerical and non-numerical, discrete and continuous);
• Software developed to solve science (e.g., biological, physical, and social), engineering, medicine, and humanities problems;
• Computer and information science that develops and optimizes the advanced system hardware, software, networking, and data management components (e.g. problem solving environments).