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

異方性多孔質媒体で満たされた波状チャネル内の定常流と脈動流の複合効果に関する研究

基本情報

論文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近似理論を確立

方法論の詳細

問題設定

異方性多孔質媒体で満たされた二次元対称波状チャネル内の非圧縮性流体の流動特性を研究する。チャネル壁面方程式は以下の通りである： $H^*(x^*) = b(1 + a\sin(2πx^*))$ ここでbは平均チャネル幅、aは振幅パラメータ、Lは波長である。

支配方程式

流動はDarcy-Brinkman方程式によって支配される： $\rho_f \frac{\partial \mathbf{V}^*}{\partial t^*} = -\nabla P^* + \mu_{eff}\nabla^2\mathbf{V}^* - \mu\mathbf{K}^{-1}\mathbf{V}^*$

異方性透過率テンソル： $\mathbf{K} = \begin{pmatrix} K_x & 0 \\ 0 & K_y \end{pmatrix}$

無次元化パラメータ

長幅比: $δ = 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領域の二層構造