RUI: effective multi-method and domain decomposition approaches to solving multiphysics problems of fluid-porous flow with transport

The National Science Foundation Division of Mathematical Sciences awarded $166,327 to The Research Foundation For The State University Of New York, operating as Farmingdale State College, on September 1, 2026, under the Mathematical and Physical Sciences program (CFDA 47.049). The award supports development of numerical algorithms and computational models for fluid-porous flow systems in which fluid partly flows freely and partly filters through a porous medium. The research addresses fundamental mathematical and algorithmic challenges in coupled free-flow and porous-media modeling, including incompatibilities between mathematical models in the two regions, accurate treatment of free-porous boundaries, solute transport, and porous-structure elasticity. The project will develop multi-method computational approaches combining finite element and boundary integral methods to effectively model flow in heterogeneous porous media in three dimensions. Applications span biotechnology (groundwater remediation, tumor treatment through blood filtration modeling, bacterial biofilm behavior in wastewater treatment), advanced manufacturing (oil and air filter design), and infectious disease contexts. The award includes support for undergraduate research experiences. Performance is located in Farmingdale, New York, with a period of performance from September 1, 2026, through August 31, 2029. The Research Foundation for the State University of New York serves as the fiscal and administrative agent for SUNY institutions.

Abstract

A fluid partly flowing freely and partly filtrating through a porous medium is fundamental in many important environmental and biological phenomena that occur in our daily lives. The development of state of the art computer models for such fluid systems is relevant in multiple areas of national interest. In biotechnology, such a model can be used in studies of groundwater remediation, filtration of blood through arterial walls in the context of optimal treatment of tumors, and the role of bacterial biofilms in wastewater treatment and infectious diseases. In advanced manufacturing, the models can help design improved oil and air filters. To study and predict such processes, IT is important to develop accurate and efficient numerical algorithms to solve the systems of underlying mathematical equations. This project will address some of the fundamental difficulties in the development of such algorithms, for example the incompatibilities of the mathematical models in the free and porous regions, an accurate treatment of transitions at the free/porous boundary, and accounting for solute transport and elasticity of the porous structure. This project will support education through undergraduate research experiences. The free flow is often modeled by the incompressible stokes equations and the porous medium by darcy?S law, coupled across the interface by the balance of forces, normal velocities, and a slip condition for the tangential velocity. This project aims to develop, analyze, and implement effective multi-method computational models in three dimensions. The main goals of the project are to (I) develop the ability to combine different methods of discretization to effectively model flow in heterogeneous porous media (e.g., using finite element methods) with boundary integral methods that have many advantages in the free flow region, (II) incorporate the transport of material by the fluid between the two domains using an advection-diffusion equation, and (III) investigate extensions to poroelastic media using the quasi-static biot system (which combines the darcy flow with elasticity of the porous structure). Robust parallel domain decomposition algorithms will be developed and analyzed to provide a theoretical basis for optimal convergence. Novel boundary integral-based methods will be developed for the advection-diffusion equation. To extend to poroelastic media, this project will investigate loosely coupled iterative schemes with time lagging and domain decomposition to solve independently the three problems: the stokes, darcy, and elasticity. The stability and accuracy will be studied in order to design algorithms with optimal convergence properties. This award reflects NSF's statutory mission and has been deemed worthy of support through evaluation using the foundation's intellectual merit and broader impacts review criteria.- subawards are not planned for this award.