🔢 mathematics

Influence of the Reynolds number on non-Newtonian flow in thin porous media

This paper uses homogenization theory to investigate how the Reynolds number affects the flow of a Carreau-type non-Newtonian fluid through a thin porous medium, identifying a critical scaling ( $\text{Re} \sim 1/\varepsilon$ ) that determines whether inertial effects influence the resulting Darcy law.

Original authors: Maria Anguiano, Matthieu Bonnivard, Francisco J. Suarez-Grau

Published 2026-02-10

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Maria Anguiano, Matthieu Bonnivard, Francisco J. Suarez-Grau

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to push honey through a very thin, sponge-like material. If you push it slowly, the honey flows smoothly. If you try to blast it through with extreme pressure, the "swirls" and "momentum" of the liquid start to matter, and the physics changes completely.

This scientific paper explores exactly that: How much "oomph" (speed and momentum) can a thick, non-Newtonian liquid have before the standard rules of flow break down?

Here is a breakdown of the paper using everyday analogies.

1. The "Non-Newtonian" Character: The Shape-Shifting Liquid

Most liquids we know, like water, are "Newtonian." They have a constant thickness (viscosity). No matter how hard you stir them, they stay just as "watery."

The paper studies Carreau fluids. These are "non-Newtonian," meaning they change their personality based on how much stress you put on them.

Pseudoplastic (Shear-thinning): Think of ketchup. It’s thick and stubborn in the bottle, but once you shake it or hit it, it becomes runny and flows easily.
Dilatant (Shear-thickening): Think of Oobleck (cornstarch and water). If you touch it gently, it’s a liquid. If you punch it, it turns into a solid brick.

2. The "Thin Porous Medium": The Microscopic Obstacle Course

The researchers aren't looking at an open pipe; they are looking at a thin porous medium.

Imagine a very thin sandwich made of layers of paper, but instead of solid paper, it’s a forest of tiny, microscopic pillars. The liquid has to navigate this dense, microscopic obstacle course. Because the "sandwich" is so thin, the liquid is essentially forced to move in a 2D plane, even though it's technically in a 3D world.

3. The "Reynolds Number": The Speed Limit of Physics

The core of the paper is about the Reynolds Number. In physics, this is a number that tells you a tug-of-war is happening between two forces:

Viscosity: The "stickiness" trying to keep things calm and orderly.
Inertia: The "momentum" trying to make things swirl and tumble.

The researchers wanted to find the "Critical Reynolds Number."

The Analogy: Imagine a crowd of people walking through a narrow, crowded hallway filled with pillars.

Low Reynolds Number: Everyone is walking slowly and carefully. They bump into pillars, but they don't crash into each other. You can predict exactly where everyone will be using a simple "flow rule" (called Darcy’s Law).
High Reynolds Number: Everyone starts sprinting. People are slamming into pillars and colliding with each other. The simple "flow rule" fails because the chaos (inertia) is now the boss.

4. The Big Discovery: The "Magic Threshold"

The scientists used complex math (called homogenization theory) to prove that there is a specific "speed limit" for these fluids.

They discovered that for this specific type of thin, pillar-filled environment, the magic threshold is of the order $1/\epsilon$ .

If the speed is below this limit, you can ignore the "chaos" and use a simple math formula (Darcy's Law) to predict the flow.
If the speed is above this limit, the "chaos" (inertia) takes over, and the simple formulas become useless.

Crucially, they found that this "speed limit" changes depending on whether the liquid is a "ketchup-type" (shear-thinning) or an "Oobleck-type" (shear-thickening).

5. Why does this matter? (The "So What?")

This isn't just math for the sake of math. This research is vital for industrial processes like:

Injection Molding: When factories melt plastic (which is a non-Newtonian fluid) and blast it into thin molds to make phone cases or car parts.
Oil and Mud Flow: Understanding how thick sludge moves through tiny cracks in the earth.

By knowing exactly when the "chaos" starts, engineers can design better machines and predict exactly how much pressure they need to use to get the job done without causing a mess.

This paper, titled "Influence of the Reynolds number on non-Newtonian flow in thin porous media," investigates the asymptotic behavior of generalized Newtonian fluids—specifically those following the Carreau law—as they flow through a thin, perforated domain.

1. The Problem

The authors address the modeling of non-Newtonian fluid flow through a thin porous medium $\Omega_\epsilon \subset \mathbb{R}^3$ . This medium is characterized by a thickness $\epsilon \ll 1$ and is perforated by a periodic array of solid cylinders of size $\epsilon$ .

The governing equations are the nonlinear stationary Navier–Stokes equations with a viscosity $\eta_r$ defined by the Carreau law:
$\eta_r(D[u]) = (\eta_0 - \eta_\infty)(1 + \lambda|D[u]|^2)^{\frac{r-2}{2}} + \eta_\infty$
The central scientific question is: At what Reynolds number ($Re$) do the inertial terms in the Navier–Stokes equations become negligible? Identifying this threshold is crucial for determining whether the complex Navier–Stokes system can be simplified into a more manageable Darcy’s law (either linear or nonlinear) for large-scale modeling.

2. Methodology

The study employs a combination of advanced mathematical techniques:

Homogenization Theory: To derive effective macroscopic equations from the microscopic description of the porous medium.
Dimension Reduction & Scaling: The authors use a change of variables (stretching the $x_3$ -direction by $1/\epsilon$ ) to handle the thinness of the domain.
Asymptotic Analysis: They assume the Reynolds number scales as $\epsilon^{-\gamma}$ and analyze the limit as $\epsilon \to 0$ .
Advanced Interpolation & Unfolding Method: To handle the technical difficulties of the inertial term and the pressure extension in perforated domains, they use an adapted "unfolding method" to establish compactness and convergence properties.
Numerical Simulation: A Newton-type method is proposed to solve the resulting nonlinear Darcy’s law, implemented via finite element methods.

3. Key Contributions and Results

The paper’s primary contribution is the identification of a critical Reynolds number that depends on the geometry of the domain and the type of fluid.

A. The Critical Reynolds Number

The authors prove that the inertial term has no influence in the limit if the Reynolds number is of order $1/\epsilon$ (i.e., $\gamma \leq 1$ ). If the Reynolds number exceeds this threshold, the inertial terms must be included in the model.

Table 1 Summary (Critical $Re$):

Domain Type	Pseudoplastic/Newtonian ( $1 < r \leq 2$ )	Dilatant ( $r > 2$ )
Thin Slab	$1/\epsilon^{3/2}$	$1/\epsilon^{1+(r/2)}$
Porous Media	$1/\epsilon^{3/2}$	$1/\epsilon^{1+(r/4)}$
Thin Porous Media	$1/\epsilon$	$1/\epsilon$

B. Derived Effective Laws (for $\gamma \leq 1$ )

Pseudoplastic/Newtonian Fluids:
- If $\gamma < 1$ : A linear 2D Darcy’s law is obtained (viscosity $\eta_0$ ).
- If $\gamma = 1$ : A nonlinear 2D Darcy’s law of Carreau type is obtained.
Dilatant Fluids:
- If $\gamma < 1$ : A linear 2D Darcy’s law is obtained.
- If $\gamma = 1$ : A nonlinear 2D Darcy’s law of Carreau type is obtained.

C. Numerical Validation

The authors demonstrate that their proposed Newton-type method effectively solves the nonlinear Darcy law. They show that for circular inclusions, the filtration velocity is aligned with the pressure gradient, whereas for elliptic inclusions, the anisotropy of the geometry causes a misalignment.

4. Significance

This research provides a rigorous mathematical foundation for engineers and physicists working with polymer injection molding, oil flow, and mud transport. By establishing the exact scaling of the Reynolds number required to use Darcy's law in thin geometries, the paper prevents the misuse of simplified models in high-velocity regimes and provides a validated numerical framework for simulating complex non-Newtonian flows in thin, structured environments.