Shear-induced force and dispersion due to buoyancy in a horizontal Hele-Shaw cell

Prabakaran Rajamanickam

2024-11-07

Abstract:This paper investigates shear flow in a Hele-Shaw cell, driven by varying horizontal buoyancy forces resulting from a horizontal density gradient induced by a scalar field. By employing asymptotic methods and taking the dependence of density and transport coefficients on the scalar field into account, effective two-dimensional hydrodynamic equations coupled with the scalar conservation equation are derived. These equations determine an effective diffusion coefficient for the scalar field accounting for shear-induced diffusion, and an effective shear-induced buoyancy force that modifies the classical Darcy's law. The derived equations provide a foundation for future research into various problems involving scalar transport in horizontal Hele-Shaw cells.

Fluid Dynamics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the shear flow caused by the horizontal density gradient in the horizontal Hele - Shaw cell and its influence on mass transfer. Specifically, the paper focuses on how the buoyancy change due to the horizontal density gradient induces shear - induced diffusion (Taylor or shear - induced diffusion) in a narrow horizontal geometry. This phenomenon is caused by the interaction between the strong shear flow on the local scale and the slowly varying scalar field, and macroscopically appears as a diffusion process.

### Main research contents
1. **Shear - induced buoyancy force**:
- The paper derives effective two - dimensional hydrodynamic equations, which take into account the dependence of density and transport coefficients on the scalar field.
- An effective shear - induced buoyancy force is introduced, which modifies the classical Darcy's law.

2. **Shear - induced diffusion coefficient**:
- The paper derives an effective diffusion coefficient, which takes into account the shear - induced diffusion effect.
- This diffusion coefficient appears as a diffusion process on a large scale, and its value is proportional to the square of the Péclet number of the flow.

### Research methods
- **Asymptotic method**: The paper uses the asymptotic analysis method to derive effective two - dimensional depth - averaged control equations from the complete three - dimensional Navier - Stokes equations.
- **Non - linear diffusion equation**: Under the Boussinesq approximation, a non - linear diffusion equation, the Erdogan–Chatwin equation, is derived.

### Research background
- **Natural convection**: When a horizontal density gradient exists, mechanical equilibrium is impossible because fluid elements at different positions but the same height are subject to different buoyancy forces.
- **Temperature or solute concentration gradient**: Such a gradient can be caused by a spatially varying scalar field (such as the fluid temperature or the concentration of a dispersed solute), which affects the density.
- **Taylor's shear - induced diffusion**: Taylor first expounded the shear - induced diffusion phenomenon in 1953, in which the shear flow is driven by an externally applied pressure gradient.

### Research significance
- **Theoretical basis**: The paper provides a theoretical basis for future research on various problems involving scalar transfer in horizontal Hele - Shaw cells.
- **Practical applications**: The research results can be applied to a variety of practical problems, such as the evolution of non - axisymmetric hot spots and the stability analysis of horizontal convection.

### Key formulas
- **Effective diffusion coefficient**:
\[
D_{\text{eff},x}^* = D^* \left[ 1 + \frac{\gamma g^2 h^8 \alpha^2}{\nu^{*2} D^{*2}} \left( \frac{\partial \theta}{\partial x^*} \right)^2 \right]
\]
where \( D^* \) is the molecular diffusion coefficient of the scalar, \(\gamma=\frac{2}{2835}\) is a numerical factor, \( g \) is the acceleration due to gravity, \(\nu^*\) is the dynamic viscosity of the fluid, \( h \) is the half - width of the channel, and \(\alpha\) is the expansion coefficient.

- **Effective buoyancy force**:
\[
F_{\text{eff}} = \gamma R_a^2 \left[ -\rho \nabla \left( \frac{(\nabla \rho)^2}{\mu \lambda} \right) - \frac{3 \rho \nabla \rho}{2 \mu^2 \lambda} \nabla \rho \cdot \nabla \mu + \rho (\nabla \rho \cdot \nabla) \left( \frac{\nabla \rho}{\mu} \right) + \frac{3 \nabla \rho}{2 \mu \Pr} \nabla \cdot \left( \frac{\rho}{\mu} \nabla \rho \right) \right]
\]