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

Basic Information

• Paper ID: 2510.08729
• Title: An investigation into the combined effect of steady and pulsatile flow within a wavy channel filled with an anisotropic porous medium
• Authors: Sanchita Pramanik, Timir Karmakar (National Institute of Technology Meghalaya)
• Classification: physics.flu-dyn math-ph math.MP
• Publication Date: October 13, 2025
• Paper Link: https://arxiv.org/abs/2510.08729

Abstract

This study provides an analytical solution to investigate the effects of steady and pulsatile flow in a symmetric wavy channel filled with an anisotropic porous medium. The flow within the channel is governed by the Darcy-Brinkman equation. The study employs the perturbation method to solve the problem, assuming that the ratio of channel width to wavelength is small (δ² ≪ 1). Under the assumption λ²δ² ≪ 1, the problem is solved to O(δ²), where λ is the anisotropic permeability ratio. The primary objective is to investigate how the anisotropic permeability ratio affects flow reversal near the crests and troughs of wavy walls, including both steady flow conditions and steady flow superimposed with pulsatile components.

Research Background and Motivation

Problem Background

1. Scientific Significance: Flow in wavy channels is a key problem in fluid mechanics with broad applications in mass transfer, heat transfer, marine engineering, and biomedical fields
2. Engineering Applications:
   • Blood oxygenation processes in membrane oxygenators
   • Blood flow simulation in the cardiovascular system (arterial stenosis and aneurysms)
   • Geothermal reservoirs and carbon capture storage systems
   • Corrugated surfaces in heat transfer enhancement devices

Research Motivation

1. Existing Limitations: Previous studies primarily focused on isotropic porous media, whereas porous media in practical applications often exhibit anisotropic characteristics
2. Technical Requirements: Fiber bundles in membrane oxygenators form anisotropic networks, requiring shape optimization to achieve efficient gas exchange
3. Theoretical Gap: The combined effects of wavy walls and pressure pulsations in anisotropic porous media leading to flow separation phenomena lack in-depth investigation

Core Contributions

1. Theoretical Innovation: First systematic investigation of the combined effects of steady and pulsatile flow in anisotropic porous media
2. Analytical Method: Provides O(δ²)-order analytical solutions based on perturbation theory, applicable to channels with small aspect ratios
3. Flow Mechanism: Reveals the mechanism by which anisotropic permeability ratio affects flow reversal near crests and troughs
4. Application Value: Provides theoretical foundation for membrane oxygenator design and cardiovascular disease diagnosis
5. Boundary Layer Analysis: Establishes Darcy approximation theory in regions away from wavy walls

Methodology Details

Problem Definition

Investigation of incompressible fluid flow characteristics in a two-dimensional symmetric wavy channel filled with anisotropic porous medium, with channel wall equation: $H^*(x^*) = b(1 + a\sin(2πx^*))$ where b is the mean channel width, a is the amplitude parameter, and L is the wavelength.

Governing Equations

Flow is governed by the Darcy-Brinkman equation: $\rho_f \frac{\partial \mathbf{V}^*}{\partial t^*} = -\nabla P^* + \mu_{eff}\nabla^2\mathbf{V}^* - \mu\mathbf{K}^{-1}\mathbf{V}^*$

Anisotropic permeability tensor: $\mathbf{K} = \begin{pmatrix} K_x & 0 \\ 0 & K_y \end{pmatrix}$

Dimensionless Parameters

• Aspect ratio: $δ = b/L$
• Darcy number: $Da = K_x/b^2$
• Anisotropy ratio: $λ^2 = K_x/K_y$
• Viscosity ratio: $M = μ_{eff}/μ$
• Womersley number: $Wo = ρ_f Ω_p b^2/μ$

Perturbation Expansion

Velocity and pressure are expanded in powers of δ²: $u_{(s,p)} = u_{(s0,p0)} + δ^2u_{(s1,p1)} + O(δ^4)$$v_{(s,p)} = v_{(s0,p0)} + δ^2v_{(s1,p1)} + O(δ^4)$$P_{(s,p)} = P_{(s0,p0)} + δ^2P_{(s1,p1)} + O(δ^4)$

Lubrication Approximation

Under the assumption δ² ≪ 1, lubrication theory is employed to solve boundary value problems. For pulsatile flow, the parameter ξ² = η² + α² is introduced, where η = √(iWo) and α² = 1/Da.

Experimental Setup

Parameter Range

• Aspect ratio: δ = 0.1-0.3
• Amplitude parameter: a = 0.3
• Darcy number: Da = 0.001-1
• Anisotropy ratio: λ = 1-2
• Womersley number: Wo = 1-100
• Pulsation amplitude: χ = 1

Boundary Conditions

1. No-slip condition at walls: $u = 0, v = 0$ at $y = H(x)$
2. Symmetry condition: $\partial u/\partial y = 0, v = 0$ at $y = 0$
3. Volume flow rate conservation: $\int_0^{H(x)} u(x,y)dy = 1$

Analysis Methods

1. Steady Flow: Solve steady Darcy-Brinkman equation by setting η = 0
2. Pulsatile Flow: Consider pure pulsatile pressure-driven Womersley flow
3. Combined Effects: Unsteady flow with steady component superimposed by pulsatile component

Experimental Results

Steady Flow Results

1. Streamline Characteristics:
   • For λ ≤ 1.5, streamlines are essentially distributed along the wall shape
   • For λ ≥ 1.75, separation and recirculation zones appear near the crest
   • Recirculation zone size expands with increasing λ
2. Pressure Gradient:
   • Negative pressure gradient (favorable gradient) at small λ values
   • Positive pressure gradient (adverse gradient) appears at the crest as λ increases
   • Adverse pressure gradient induces flow separation
3. Wall Shear Stress:
   • Shear stress approaches zero at λ = 1.5
   • For λ ≥ 1.75, shear stress becomes negative near x = 0.25, indicating flow separation

Pulsatile Flow Results

1. Vortex Evolution:
   • t = 0-Tp/4: Fluid flows in the main flow direction
   • t = 2Tp/4: Vortex formation near the crest
   • t = 3Tp/4: Vortex extends into the main flow, complete flow reversal
2. Anisotropy Effects:
   • λ = 1 (isotropic): No vortex formation at t = 0 and Tp/4
   • λ = 2 (strong anisotropy): Obvious circulation observed at all time instances

Combined Effects Results

1. Time Dependence: Onset time and location of flow reversal depend on the interaction between λ and time t
2. Pressure Gradient Variation: Adverse pressure gradient appears at t = 3Tp/4 and 4Tp/4
3. Pulsatile Boundary Layer Thickness: $δ_{puls} = 1/S$, where $S = \frac{1}{\sqrt{M}}(α^4 + Wo^2)^{1/4}\cos\left(\frac{1}{2}\tan^{-1}\left(\frac{Wo}{α^2}\right)\right)$

Pressure Drop Analysis

Steady pressure drop increases rapidly with α² (inverse of Darcy number), indicating that low permeability significantly increases flow resistance. Increase in amplitude a also leads to higher pressure drop, reflecting the effect of surface roughness.

Related Work

Wavy Channel Flow Research

• Tsangaris & Leiter (1984): Perturbation method for viscous flow
• Pozrikidis (1987): Numerical study of Stokes flow
• Wang & Chen (2002): Forced convection heat transfer

Porous Media Flow

• Ng & Wang (2010): Darcy-Brinkman flow
• Wei et al. (2003): Wavy channel with porous elastic layer
• Karmakar & Raja Sekhar (2017): Flow reversal in anisotropic porous media

Pulsatile Flow Applications

• Womersley (1955): Theoretical foundation of arterial pulsatile flow
• Sobey (1980-1983): Flow in membrane oxygenators
• Griffith et al. (2009): Pulsatile flow in stenotic geometries

Conclusions and Discussion

Main Conclusions

1. Anisotropy Effects: Increasing λ significantly promotes flow separation, forming recirculation zones near crests and troughs
2. Temporal Evolution: Vortex formation and dissipation in pulsatile flow exhibit obvious periodicity
3. Combined Mechanism: The combined effects of steady and pulsatile components alter the timing and location of flow separation
4. Boundary Layer Structure: Dual-layer structure of near-wall Brinkman layer and far-field Darcy region

Engineering Applications

1. Membrane Oxygenators: Anisotropy and pulsation can optimize oxygen transfer efficiency
2. Cardiovascular Diagnosis: Provides theoretical foundation for blood flow analysis in arterial stenosis and aneurysms
3. Heat Transfer Enhancement: Flow separation at corrugated surfaces enhances convective mixing

Limitations

1. Small Parameter Assumption: Only applicable when δ² ≪ 1 and λ²δ² ≪ 1
2. Two-Dimensional Simplification: Three-dimensional effects in practical applications are not considered
3. Linearization: Perturbation method limits applicability to large amplitude cases

Future Directions

1. Three-Dimensional Extension: Consider three-dimensional wavy channels with spanwise variations
2. Nonlinear Analysis: Nonlinear effects in large amplitude cases
3. Experimental Validation: Comparison with experimental data
4. Optimization Design: Geometric parameter optimization based on flow control

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Complete mathematical derivation based on perturbation theory
2. Physical Insight: Reveals the underlying mechanism of anisotropic permeability on flow separation
3. Application Value: Provides practical theoretical tools for membrane oxygenators and cardiovascular flow
4. Systematicity: Comprehensive analysis covering steady, pulsatile, and combined effects
5. Boundary Layer Theory: Establishes Brinkman-Darcy dual-layer structure theory

Weaknesses

1. Parameter Constraints: Small parameter assumptions limit the generality of the method
2. Experimental Absence: Lacks experimental data to validate theoretical predictions
3. Numerical Comparison: No comparison with CFD numerical simulation results
4. Material Characteristics: Insufficient discussion of anisotropic characteristics of actual porous materials

Impact

1. Academic Contribution: Provides new analytical framework for porous media flow theory
2. Engineering Value: Direct application prospects in biomedical engineering and heat transfer enhancement
3. Methodological Innovation: Perturbation-lubrication combined method can be generalized to other complex geometric flows

Applicable Scenarios

1. Microfluidic Devices: Flow control in small-scale corrugated channels
2. Biomedical Engineering: Vascular flow and artificial organ design
3. Industrial Applications: Optimization of heat exchangers and separation devices
4. Geological Engineering: Fluid transport analysis in porous rock layers

References

The paper cites 62 relevant references spanning multiple fields including fluid mechanics, porous media theory, and biomedical engineering, providing a solid theoretical foundation for this research.