Wall slip effects on the fiber orientation of a short-fiber suspension in hyperbolic channel flow

This study numerically investigates how increasing wall slip in a hyperbolic channel flow enhances the alignment of short rigid fibers by extending the region of high orientation from the midplane toward the walls, utilizing an analytical velocity field and a second-order tensor formulation with hybrid closure.

The Gist

Imagine you are pouring a thick milkshake filled with tiny, rigid sticks (like toothpicks) through a funnel. As the milkshake flows, those sticks don't just float randomly; they start to line up, pointing in the direction the liquid is moving. This alignment is crucial because if you were to freeze that milkshake into a solid block (like a plastic part in a car), the strength of that block would depend entirely on how those sticks are arranged.

This paper is a scientific study about how to control that alignment, specifically by looking at what happens when the walls of the funnel are "slippery."

Here is the breakdown of the research using simple analogies:

1. The Setup: The "Hyperbolic Funnel"

The researchers are studying a specific shape of channel called a hyperbolic channel.

• The Analogy: Imagine a hallway that starts wide and gets narrower and narrower as you walk down it, but the walls curve inward smoothly (like a funnel or an hourglass).
• The Flow: As the fluid (the milkshake) moves from the wide part to the narrow part, it has to speed up.
  • Near the walls: The fluid rubs against the sides, creating a "shear" effect (like sliding a deck of cards). The sticks here tend to spin and align at an angle.
  • In the middle: The fluid is being stretched out like taffy. This is called "extensional flow." Here, the sticks get pulled straight and line up perfectly with the flow.

2. The Problem: The "Sticky" vs. "Slippery" Walls

In most standard physics models, we assume the fluid sticks perfectly to the wall (like tape on a surface). This is called the no-slip condition.

• The Reality: However, in many modern manufacturing processes (like 3D printing or injection molding), the fluid might be slippery, or the walls might be treated to be smooth. This creates wall slip.
• The Question: The researchers wanted to know: If the walls are slippery, does the fluid slide past them faster? And if so, how does that change the way the tiny sticks inside line up?

3. The Discovery: The "Slippery Slide" Effect

The team used complex math (which they call "tensor equations") to simulate this flow. Here is what they found, translated into everyday terms:

• The "Slippery" Effect on Speed: When the walls are slippery, the fluid near the edges doesn't get stuck. It slides along. This makes the flow more uniform. Instead of the fluid moving super fast in the middle and barely moving at the edges, the whole stream moves more evenly.
• The "Slippery" Effect on Alignment:
  • In the Middle: Surprisingly, even though the flow changes, the sticks in the very center still line up perfectly straight, just like before. The "stretching" force is so strong there that it overpowers the slip.
  • Near the Walls (The Big Change): This is where the magic happens. When the walls are slippery, the "zone of perfect alignment" (where the sticks are straight) expands.
  • The Metaphor: Imagine a crowd of people walking through a narrowing hallway.
    • No Slip (Sticky Walls): People near the walls get stuck and have to shuffle sideways. Only the people in the very center can walk straight.
    • Slip (Slippery Walls): The people near the walls can glide past the walls easily. Now, not just the center people, but also the people standing closer to the walls can walk straight and align with the crowd. The "straight-line zone" gets wider.

4. Why Does This Matter?

The researchers found that by making the walls slippery (or by designing the channel to encourage slip), you can force more of the fibers to align in the direction of the flow.

• The Result: You end up with a material that is stronger and more uniform.
• The Application: This is huge for industries like 3D printing or making carbon-fiber car parts. If you know how to tweak the "slipperiness" of the mold, you can engineer the final product to be stronger in the specific direction you need it to be.

Summary

Think of this paper as a guide on how to use slippery walls to turn a chaotic crowd of tiny sticks into a well-organized army. By letting the fluid slide along the edges, the "army" (the fibers) aligns more perfectly throughout the entire channel, not just in the center. This allows engineers to build stronger, better materials by simply changing how the fluid interacts with the walls.

Technical Summary

1. Problem Statement

The study investigates how wall slip influences the orientation of non-Brownian, short, rigid, high-aspect-ratio cylindrical fibers suspended in a Newtonian fluid flowing through a symmetric hyperbolic planar channel.

• Context: Fiber orientation is critical for determining the mechanical properties of final composite materials (e.g., in injection molding and 3D printing).
• Flow Geometry: The channel features an entrance region with parallel walls that transitions smoothly into a hyperbolic contraction. This geometry creates a flow field that transitions from pure shear near the walls to pure extension along the midplane.
• Gap in Literature: While previous work (Housiadas et al., 2025) analyzed this flow under standard no-slip boundary conditions, the effect of wall slip (common in polymeric fluids and treated surfaces) on fiber orientation in such mixed kinematic flows had not been systematically quantified.

2. Methodology

The authors employed a semi-analytical and numerical approach based on the following framework:

• Governing Equations:
  • Flow Field: The fluid is treated as a Newtonian matrix. The fiber extra-stress contribution to the total stress is neglected (decoupled approach), valid for dilute-to-semi-dilute suspensions with high aspect ratios ($r > 10$). The velocity field is derived analytically using extended lubrication theory (Sialmas & Housiadas, 2024) with Navier's slip law applied at the walls.
  • Fiber Orientation: The evolution of the second-order orientation tensor ($\mathbf{a}$) is modeled using the Advani & Tucker (1987) formulation.
    • Closure: A hybrid closure (convex combination of linear and quadratic closures) is used to approximate the fourth-order tensor, balancing accuracy and computational efficiency.
    • Interactions: Fiber-fiber interactions are modeled via a phenomenological rotary diffusion term proportional to the local rate-of-deformation magnitude.
• Coordinate Transformation: To handle the varying channel height, the authors mapped the physical domain to a fixed rectangular domain using transformed coordinates $(Y, Z)$. They also introduced a rotated coordinate system $(s, n, x)$ aligned with the local flow direction to simplify the boundary conditions and tensor evolution equations.
• Numerical Solution:
  • The orientation tensor equations were solved using a fully implicit finite difference method.
  • Integration along the flow direction ($Z$) was treated as an initial value problem, starting from the analytical solution of fully developed Poiseuille flow in the entrance region.
  • A Newton iterative scheme was used to solve the resulting non-linear algebraic equations at each cross-section.

3. Key Contributions

1. Extension of No-Slip Models: The work extends previous no-slip analyses to include wall slip, quantifying how slip modifies the kinematic environment and subsequent fiber alignment.
2. Decoupled Kinematic Analysis: It isolates the effect of slip on orientation by utilizing a pre-calculated Newtonian velocity field, demonstrating that slip significantly alters the shear/extension balance without requiring complex coupled flow-stress simulations for this regime.
3. Validation of Infinite Aspect Ratio: The study justifies the use of the infinite-aspect-ratio limit for fibers with aspect ratios $r \ge 20$, showing that finite aspect ratio effects are negligible for practical industrial fibers in this flow regime (Appendix B).
4. Analytical-Numerical Hybrid: The paper successfully combines high-order analytical solutions for the velocity field with robust numerical integration for the orientation tensor, ensuring high accuracy in capturing boundary layer effects.

4. Key Results

• Effect of Slip on Flow Kinematics:
  • Increasing the slip coefficient ($K$) reduces the magnitude of the rate-of-deformation tensor ($\dot{\gamma}$) throughout the domain.
  • Unlike the no-slip case where the midplane strain rate is constant, slip introduces a smooth spatial variation in the strain rate, though the flow remains purely extensional at the midplane.
• Fiber Orientation Evolution:
  • Walls: Fibers near the walls exhibit a "boundary layer" effect where they align tangentially to the wall immediately upon entering the hyperbolic section. Higher slip coefficients lead to stronger initial alignment due to increased slip velocity.
  • Midplane: The flow remains purely extensional. Surprisingly, despite a significant reduction in the strain rate magnitude due to slip, the degree of fiber alignment at the midplane remains practically unchanged. This is because the governing ODE depends on the ratio of velocity to strain rate, which varies linearly with $Z$ and is insensitive to $K$.
• Overall Alignment Trends:
  • As the slip coefficient increases, the region of high fiber alignment (dominated by extensional flow) expands from the midplane toward the walls.
  • The transition from wall-aligned (shear-dominated) to midplane-aligned (extension-dominated) profiles becomes more gradual and monotonic with higher slip.
  • Eigenvector Analysis: Increasing slip reduces the angle between the major fiber orientation eigenvector and the flow direction, indicating faster and stronger alignment with the main flow across the entire channel cross-section.

5. Significance and Implications

• Process Control: The study demonstrates that wall slip is not merely a boundary condition artifact but a controllable parameter for tailoring fiber orientation in composite manufacturing. By adjusting surface treatments or fluid properties to induce slip, manufacturers can enhance fiber alignment in the final product.
• Design Optimization: The findings suggest that hyperbolic channel designs can be optimized in conjunction with slip conditions to achieve desired anisotropic mechanical properties in extruded or molded parts.
• Theoretical Advancement: The work provides a rigorous benchmark for future studies involving viscoelastic matrices or fully coupled fiber-stress interactions, establishing a baseline for slip effects in mixed shear-extension flows.

In conclusion, the paper establishes that while wall slip reduces the local deformation rates, it paradoxically enhances the overall efficiency of fiber alignment by expanding the extensional flow regime and reducing the shear-dominated boundary layer thickness.