A Quasi-Analytical Solution for "Carreau-Yasuda-like" Shear-thinning Fluids Flowing in Slightly Tapered Pipes

This paper presents a quasi-analytical framework for modeling the flow of Carreau-Yasuda-like shear-thinning fluids in slightly tapered pipes, deriving expressions via order-of-magnitude analysis and validating them through numerical simulations applied to extrusion bioprinting.

The Gist

The Big Picture: The "Traffic Jam" in a Funnel

Imagine you are squeezing a tube of toothpaste or honey through a funnel.

• The Fluid: Most people think of fluids like water, which flows the same way no matter how fast you push them. But many real-world fluids (like ketchup, paint, blood, or the "bio-inks" used to 3D print human tissue) are shear-thinning. This means they are thick and gooey when sitting still, but if you push them hard or fast, they suddenly become runny and thin, like water.
• The Pipe: The paper looks at fluids flowing through a tapered pipe—basically a cone that gets narrower at the end, like a nozzle on a spray bottle or a syringe.
• The Problem: Scientists have known how to calculate the flow for simple fluids (like water) or simple "thickening" fluids for a long time. But for these complex "thinning" fluids (specifically those following the Carreau-Yasuda model), there was no easy math formula. You usually had to use super-computers to guess the answer, which is slow and expensive.

The Solution: The "Three-Zone" Map

The authors of this paper created a quasi-analytical solution. In plain English, they built a "smart map" that predicts exactly how the fluid moves without needing a super-computer.

They did this by simplifying the fluid's behavior into three distinct zones, like a traffic system:

1. The "Calm" Zone (Low Speed): Near the center of the pipe, where the fluid moves slowly, it acts like thick honey. It's stubborn and doesn't want to move.
2. The "Slippery" Zone (Medium Speed): As you get closer to the walls, the fluid gets squeezed and speeds up. Here, it starts to thin out and slide easily, acting like a power-law fluid (a specific mathematical rule for thinning).
3. The "Runaway" Zone (High Speed): Right at the very edge of the pipe, where the fluid is moving fastest, it becomes as thin as water.

The Analogy: Imagine a crowded hallway (the pipe).

• In the middle, people are walking slowly and bumping into each other (thick viscosity).
• As they get closer to the exit, they start running and weaving through gaps (shear-thinning).
• Right at the door, they are sprinting in a single file line (thin viscosity).

The paper's math figures out exactly where these "zones" start and stop at every point along the pipe.

How They Did It: The "Lego" Approach

The real math for these fluids is incredibly messy and hard to solve. The authors' trick was to approximate the fluid's behavior using a Piecewise Approximation (PWA).

Think of the fluid's viscosity curve as a smooth, curvy hill.

• Old way: Trying to solve the math for the whole curvy hill at once is like trying to climb a jagged mountain with your eyes closed.
• New way: The authors broke that curvy hill into three flat, straight steps (like a Lego staircase).
  • Step 1: Flat plateau (Thick).
  • Step 2: A ramp (Thinning).
  • Step 3: Flat plateau (Thin).

By treating the fluid as these three simple steps, they could use standard, easy math to solve the problem. They then stitched these three solutions together to get a complete picture of the flow.

Why Does This Matter? (The "Bio-Printer" Test)

The authors tested their new math by simulating Extrusion Bioprinting. This is a high-tech process where robots print living tissue (like skin or cartilage) by squeezing bio-ink through a tiny needle.

• The Risk: If you squeeze the ink too hard, the shear stress (the "rubbing" force) at the wall of the needle gets too high. This can kill the cells inside the ink. If you squeeze too gently, the print doesn't hold its shape.
• The Result: Their new "smart map" allows engineers to instantly calculate: "If I push the plunger at this speed, will the cells die?"

They compared their math to a full computer simulation (the "super-computer" method) and found their solution was 97-99% accurate but calculated instantly.

The Takeaway

This paper gives engineers a fast, easy-to-use calculator for designing nozzles for complex fluids.

• Before: You had to run a slow, complex computer simulation to see if your design would work.
• Now: You can use this formula to instantly check if your pipe is too narrow, if your fluid is too thick, or if you are pushing too hard and killing the cells.

It's like going from having to build a full-scale wind tunnel to test a car, to just having a simple app on your phone that tells you if the car will fly. This is a huge step forward for making better 3D printed organs, better food processing, and more efficient plastic manufacturing.

Technical Summary

1. Problem Statement

The paper addresses the mathematical modeling of non-Newtonian inelastic fluids (specifically "Carreau-Yasuda-like" fluids) flowing through slightly tapered pipes (conical ducts).

• Context: This flow configuration is ubiquitous in industrial and research applications, including polymer manufacturing, food processing, and biomedical extrusion bioprinting.
• The Gap: While analytical solutions exist for Newtonian fluids and simpler non-Newtonian models (like Power-Law or Bingham) in tapered pipes, a closed-form or quasi-analytical solution for fluids exhibiting the complex rheology of the Carreau-Yasuda model has been lacking. The Carreau-Yasuda model is characterized by two viscosity plateaus (at low and high shear rates) connected by a shear-thinning branch, making exact analytical integration difficult due to the nonlinearity of the viscosity function.

2. Methodology

The authors developed a quasi-analytical framework based on the following steps:

A. Rheological Approximation (PWA)

To overcome the mathematical intractability of the continuous Carreau-Yasuda equation, the authors introduced a Power-Law-based Piecewise Approximation (PWA). This model approximates the viscosity curve $\mu(\dot{\gamma})$ as three distinct linear regions in a log-log space:

1. Low-Shear Rate (LSR) Region: Constant Newtonian viscosity ($\mu_0$) for $\dot{\gamma} \leq 1/\lambda_0$.
2. Medium-Shear Rate (MSR) Region: Power-law behavior ($K\dot{\gamma}^{n-1}$) for $1/\lambda_0 < \dot{\gamma} \leq 1/\lambda_\infty$.
3. High-Shear Rate (HSR) Region: Constant Newtonian viscosity ($\mu_\infty$) for $\dot{\gamma} > 1/\lambda_\infty$.

B. Mathematical Formulation

• Governing Equations: The study assumes axisymmetric, steady, creeping flow (neglecting inertial terms in the Navier-Stokes equations), which is valid for low Reynolds numbers common in bioprinting and polymer extrusion.
• Flow Regimes: Depending on the local wall shear stress, the flow in any cross-section can be classified into three regimes:
  • LSR Flow: Entirely Newtonian ($\mu_0$).
  • MSR Flow: A core of Newtonian fluid ($\mu_0$) surrounded by a shear-thinning annulus ($K, n$).
  • HSR Flow: A core of Newtonian fluid ($\mu_0$), a middle shear-thinning annulus, and an outer Newtonian annulus ($\mu_\infty$).
• Derivation: By applying order-of-magnitude analysis and integrating the momentum balance equation, the authors derived expressions for axial velocity ($v_z$), radial velocity ($v_r$), shear rate, and flow rate ($Q$) for each regime.
• Iterative Algorithm: Since the pressure gradient $p'(z)$ varies along the tapered pipe and determines the boundaries between viscosity regions ($R_0(z)$ and $R_\infty(z)$), a closed-form relationship between $Q$ and total pressure drop is impossible. The authors developed Algorithm 1, an iterative semi-analytical procedure to compute the axial pressure gradient distribution $p'(z)$ by satisfying mass conservation and matching the flow rate $Q$.

3. Key Contributions

1. First Quasi-Analytical Solution for Carreau-Yasuda in Tapered Pipes: The paper provides the first analytical framework capable of handling the specific rheological complexity of Carreau-Yasuda fluids in converging/diverging geometries.
2. Piecewise Approximation Strategy: The introduction of the PWA allows for the derivation of analytical expressions for velocity and pressure fields that would otherwise require full numerical simulation.
3. Algorithm for Pressure Gradient: The development of a robust iterative algorithm (Algorithm 1) that determines the pressure gradient profile along the pipe by dynamically switching between Newtonian and Power-Law sub-solutions based on local shear stress.
4. Application to Bioprinting: The framework is specifically tailored for extrusion bioprinting, offering a tool to predict pressure drops and shear stresses to prevent cell damage.

4. Results and Verification

The proposed solution was validated through two main approaches:

A. Numerical Verification (CFD)

• Setup: The quasi-analytical results were compared against Computational Fluid Dynamics (CFD) simulations (Ansys Fluent) using the full, non-approximated Carreau-Yasuda (SRB) model.
• Fluids: Two bio-ink-like fluids (Fluid A and Fluid B) with different shear-thinning indices were tested.
• Accuracy:
  • Pressure Drop: The quasi-analytical solution matched CFD results with high accuracy. Deviations were minimal (approx. 2.5% at the inlet for Fluid A, increasing to ~4.1% for larger taper angles).
  • Velocity Profiles: Excellent agreement was observed for both axial and radial velocity components.
  • Parametric Sensitivity: Varying rheological parameters ($\mu_0, \lambda_0, n, \lambda_\infty$) by $\pm 20\%$ resulted in mean extrusion velocity errors generally below 2.6%, confirming the robustness of the PWA.

B. Application Analysis

• Flow Regime Transitions: The model successfully visualized how increasing the pressure drop transitions the flow from LSR to MSR and finally to HSR regimes along the pipe axis.
• Cell Damage Prediction: The framework was used to calculate the maximum allowable extrusion velocity ($\bar{V}_{out, dam}$) before the wall shear stress exceeds a critical threshold for cell damage.
  • Results showed that reducing zero-shear viscosity ($\mu_0$) or the power-law index ($n$) allows for higher extrusion velocities without damaging cells.
  • The model demonstrated that for a higher damage threshold (500 Pa vs. 100 Pa), the allowable extrusion velocity could increase by a factor of ten.

5. Significance and Conclusion

• Engineering Utility: The solution provides a fast, ready-to-use screening tool for engineers and researchers. Unlike full CFD, which is computationally expensive, this quasi-analytical approach allows for rapid parametric studies to optimize nozzle geometry and operating conditions.
• Biomedical Impact: In extrusion bioprinting, where cell viability is critical, this method offers a direct way to calculate safe operating windows (pressure vs. flow rate) to avoid shear-induced cell death.
• Limitations: The method relies on the "slightly tapered" assumption (lubrication theory), meaning it is valid for small taper angles where radial velocity components are negligible compared to axial ones. For sharp contractions or large taper angles, full numerical methods remain necessary.

In summary, this paper bridges the gap between complex rheological modeling and practical engineering design for tapered pipe flows, offering a validated, efficient analytical alternative to numerical simulations for Carreau-Yasuda fluids.