2025-11-20T17:19:15.817110

An investigation into the combined effect of steady and pulsatile flow within a wavy channel filled with an anisotropic porous medium

Pramanik, Karmakar

We provide an analytical solution to examine the impact of steady and pulsatile flow in a symmetric wavy channel filled with an anisotropic porous medium. The flow inside the wavy channel with a porous substrate is governed by the Darcy-Brinkman equation. We used the perturbation approach to ascertain the solution to the problem, assuming that the ratio between the channel width and the wavelength is very small (i.e., $Î´^2 \ll 1$). We solved the problem up to $O(Î´^2)$ by assuming $Î»^2 Î´^2 \ll 1$, where $Î»$ is the ratio of anisotropic permeability. Our primary objective is to investigate how the anisotropic permeability ratio influences flow reversal near the crest and trough of a wavy wall under both steady flow conditions and when a pulsatile component is superimposed on the steady flow. In the case of steady flow, flow separation occurs near the crest region of the wavy wall, leading to the formation of recirculating zones where fluid becomes trapped. In contrast, for pulsatile flow as well as combined effects of steady and pulsatile flow, vortices form transiently at certain phases of the pulsation cycle. Subsequently, these vortices bulge into the mainstream, resulting in a complete backflow. We explored the approximate flow reversal zone using Darcy's approximation theory far from the wavy wall. The flow recirculation associated with steady flow is critical for mass transport and enhances convective mixing. On the other hand, flow reversal at different time instances can optimize oxygen transfer in devices such as membrane oxygenators.

academic

An investigation into the combined effect of steady and pulsatile flow within a wavy channel filled with an anisotropic porous medium

基本信息

论文ID: 2510.08729
标题: An investigation into the combined effect of steady and pulsatile flow within a wavy channel filled with an anisotropic porous medium
作者: Sanchita Pramanik, Timir Karmakar (National Institute of Technology Meghalaya)
分类: physics.flu-dyn math-ph math.MP
发表时间: 2025年10月13日
论文链接: https://arxiv.org/abs/2510.08729

摘要

本研究提供了一个解析解来研究各向异性多孔介质填充的对称波浪通道中稳态和脉动流的影响。通道内的流动由Darcy-Brinkman方程控制。研究采用摄动方法求解问题，假设通道宽度与波长之比很小(δ² ≪ 1)。在λ²δ² ≪ 1的假设下，将问题求解至O(δ²)阶，其中λ为各向异性渗透率比。主要目标是研究各向异性渗透率比如何影响波浪壁面波峰和波谷附近的流动反转，包括稳态流动条件和稳态流叠加脉动分量的情况。

研究背景与动机

问题背景

科学意义: 波浪通道中的流动是流体力学中的关键问题，在传质传热、海洋工程、生物医学等领域具有广泛应用
工程应用:
- 膜式氧合器中的血液氧合过程
- 心血管系统中的血流模拟(动脉狭窄和动脉瘤)
- 地热储层和碳捕获存储系统
- 传热强化设备中的波纹表面

研究动机

现有局限: 以往研究主要关注各向同性多孔介质，而实际应用中多孔介质往往呈现各向异性特征
技术需求: 膜式氧合器中纤维束形成各向异性网络，需要优化形状以实现高效气体交换
理论缺口: 各向异性多孔介质中波浪壁面和压力脉动的联合效应导致的流动分离现象缺乏深入研究

核心贡献

理论创新: 首次系统研究了各向异性多孔介质中稳态和脉动流的联合效应
解析方法: 基于摄动理论提供了O(δ²)阶的解析解，适用于小长宽比通道
流动机理: 揭示了各向异性渗透率比对波峰和波谷附近流动反转的影响机制
应用价值: 为膜式氧合器设计和心血管疾病诊断提供了理论基础
边界层分析: 建立了远离波浪壁面区域的Darcy近似理论

长宽比: $δ = b/L$
Darcy数: $Da = K_x/b^2$
各向异性比: $λ^2 = K_x/K_y$
粘度比: $M = μ_{eff}/μ$
Womersley数: $Wo = ρ_f Ω_p b^2/μ$

长宽比: δ = 0.1-0.3
振幅参数: a = 0.3
Darcy数: Da = 0.001-1
各向异性比: λ = 1-2
Womersley数: Wo = 1-100
脉动振幅: χ = 1

边界条件

壁面无滑移条件: $u = 0, v = 0$ at $y = H(x)$
对称条件: $\partial u/\partial y = 0, v = 0$ at $y = 0$
体积流量守恒: $\int_0^{H(x)} u(x,y)dy = 1$

分析方法

稳态流动: 设置η = 0求解稳态Darcy-Brinkman方程
脉动流动: 考虑纯脉动压力驱动的Womersley流
联合效应: 稳态流叠加脉动分量的非定常流动

实验结果

稳态流动结果

流线特征:
- λ ≤ 1.5时，流线基本沿壁面形状分布
- λ ≥ 1.75时，波峰附近出现分离和回流区
- 回流区大小随λ增加而扩大
压力梯度:
- 小λ值时压力梯度为负(有利梯度)
- λ增大时波峰处出现正压力梯度(逆压梯度)
- 逆压梯度导致流动分离
壁面剪应力:
- λ = 1.5时剪应力接近零
- λ ≥ 1.75时剪应力在x = 0.25附近变负，表明流动分离

脉动流动结果

涡旋演化:
- t = 0-Tp/4: 流体沿主流方向流动
- t = 2Tp/4: 波峰附近形成涡旋
- t = 3Tp/4: 涡旋扩展至主流，完全反向流动
各向异性影响:
- λ = 1(各向同性): t = 0和Tp/4时无涡旋形成
- λ = 2(强各向异性): 各时刻均可观察到明显环流

联合效应结果

时间依赖性: 流动反转的起始时间和位置取决于λ和时间t的相互作用
压力梯度变化: t = 3Tp/4和4Tp/4时出现逆压梯度
脉动边界层厚度: $δ_{puls} = 1/S$ ，其中 $S = \frac{1}{\sqrt{M}}(α^4 + Wo^2)^{1/4}\cos\left(\frac{1}{2}\tan^{-1}\left(\frac{Wo}{α^2}\right)\right)$

各向异性效应: λ增大显著促进流动分离，在波峰和波谷附近形成回流区
时间演化: 脉动流中涡旋的形成和消失呈现明显的时间周期性
联合机制: 稳态和脉动的联合效应改变了流动分离的时间和位置
边界层结构: 近壁面Brinkman层和远场Darcy区的双层结构