The Diffuse Interface Method and Applications to Coupled Systems in Fluid Dynamics

涉及运动域的系统，其中流体与相邻区域相互作用，通常使用偏微分方程建模，其中耦合条件在公共界面上成立。这种系统出现在医学（肝脏灌注，淋巴循环，心脏瓣膜的关闭和打开），地质力学（裂缝传播，地表和地下水流之间的耦合）和其他应用中。这种系统的数值模拟是基于离散近似的控制方程。为了准确地描述动力学，界面跟踪方法，其中计算网格的节点与界面的参数表示对齐，经常使用。然而，界面跟踪方法迅速变得难以应用时域变形很大。为了防止数值故障，计算昂贵的网格再生或类似的技术是必要的，当网格元素变得高度歪斜。漫射界面方法是基于固定网格方法的替代策略。对于这种方法，使用相场函数重新制定模型，该相场函数从一个区域中的零平滑地过渡到另一个区域中的一个。计算网格节点不必与界面对齐，其位置现在使用相场函数捕获。这种方法是有用的，即使当域不随时间变化，或在两个区域之间的界面很难准确确定的情况下，或当界面的几何形状是复杂的。然而，扩散界面方法在界面处引入了额外的误差，这需要仔细控制。本计画旨在为扩散界面法在流体力学中的应用建立数学基础。在这项工作中开发的技术预计将适用于其他耦合系统，涉及流体和多孔弹性和/或弹性结构以及。该项目包括通过参与研究对研究生进行培训。该项目的重点是发展适用于流体动力学耦合系统的扩散界面方法的数学理论和数值方法。这将通过研究一系列耦合流动模型来实现，包括流体-多孔介质相互作用、流体-多孔弹性结构相互作用和流体-弹性结构相互作用。对于每个模型，工作需要证明的适定性的基本扩散界面问题，显示的扩散界面模型收敛到相应的尖锐的界面模型的界面层的宽度为零，并计算收敛速度，包括建模误差和基于有限元方法的离散解的近似误差。将使用加权Sobolev空间进行分析。该奖项反映了NSF的法定使命，并被认为是值得通过使用基金会的知识价值和更广泛的影响审查标准进行评估的支持。

Systems involving moving domains in which a fluid interacts with a neighboring region are often modeled using partial differential equations with coupling conditions that hold across a common interface. Such systems occur in medicine (liver perfusion, lymphatic circulation, the closing and opening of heart valves), geomechanics (fracture propagation, coupling between the surface and groundwater flows), and other applications. Numerical simulations of such systems are based on discrete approximations to the governing equations. To accurately describe the dynamics, interface tracking methods, in which the nodes of the computational mesh are aligned with the parametric representation of the interface, are often used. However, interface-tracking methods rapidly become difficult to apply when domain deformations are large. To prevent numerical failures, computationally expensive mesh regenerations or similar techniques are necessary when mesh elements become highly skewed. The diffuse interface method is an alternative strategy based on a fixed mesh approach. For this method, the model is reformulated using a phase-field function that smoothly transitions from zero in one region to one in the other region. The computational mesh nodes do not have to be aligned with the interface, whose location is now captured using the phase-field function. This approach is useful even when the domain does not change in time, or in cases where the interface between the two regions is difficult to determine exactly, or when the geometry of the interface is complex. However, the diffuse interface method introduces an additional error at the interface, which needs to be carefully controlled. This project aims to establish mathematical foundations for application of the diffuse interface method in fluid dynamics. The techniques developed in this work are expected to be applicable to other coupled systems involving fluids and poroelastic and/or elastic structures as well. The project includes training of graduate students through involvement in the research.This project focuses on the development of mathematical theory and numerical methods for the diffuse interface method applied to coupled systems in fluid dynamics. This will be achieved by studying a hierarchy of coupled flow models, including the fluid-porous medium interaction, fluid-poroelastic structure interaction, and fluid-elastic structure interaction. For each model, the work entails proving the well-posedness of the underlying diffuse interface problem, showing the convergence of the diffuse interface model to the corresponding sharp interface model as the width of the interfacial layer goes to zero, and calculating the rate of convergence including the modeling error and the approximation error of the discrete solution based on the finite element method. The analysis will be performed using weighted Sobolev spaces. Numerical methods will be developed and implemented for each model.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.