Squeeze flow

Squeeze flow (also called squeezing flow, squeezing film flow, or squeeze flow theory) is a type of flow in which a material is pressed out or deformed between two parallel plates or objects. First explored in 1874 by Josef Stefan, squeeze flow describes the outward movement of a droplet of material, its area of contact with the plate surfaces, and the effects of internal and external factors such as temperature, viscoelasticity, and heterogeneity of the material. Several squeeze flow models exist to describe Newtonian and non-Newtonian fluids undergoing squeeze flow under various geometries and conditions. Numerous applications across scientific and engineering disciplines including rheometry, welding engineering, and materials science provide examples of squeeze flow in practical use.

Basic Assumptions

Conservation of mass (expressed as a continuity equation), the Navier-Stokes equations for conservation of momentum, and the Reynolds number provide the foundations for calculating and modeling squeeze flow. Boundary conditions for such calculations include assumptions of an incompressible fluid, a two-dimensional system, neglecting of body forces, and neglecting of inertial forces.

Relating applied force to material thickness:

$F=-{\frac {4*L^{3}*\eta *W}{h^{3}}}{dh \over dt}$

Where $F$ is the applied squeezing force, $2L$ is the initial length of the droplet, $\eta$ is the fluid viscosity, $W$ is the width of the assumed rectangular plate, $2h$ is the final height of the droplet, and ${dh \over dt}$ is the change in droplet height over time. To simplify most calculations, the applied force is assumed to be constant.

Visualization of surface-to-surface contact; asperity highlighted in red.

Newtonian fluids

Several equations accurately model Newtonian droplet sizes under different initial conditions.

Consideration of a single asperity, or surface protrusion, allows for measurement of a very specific cross-section of a droplet. To measure macroscopic squeeze flow effects, models exist for two the most common surfaces: circular and rectangular plate squeeze flows.

Single asperity squeeze flow diagram at initial and follow-on conditions; plates (assumed to be semi-infinite, in gray), droplet (green).

Single asperity

For single asperity squeeze flow:

${\frac {h_{0}}{h}}=\left(1+{\frac {5*F*t*h_{0}^{2}}{4*\eta *W*L_{0}^{3}}}\right)^{1/5}$

Where $2h_{0}$ is the initial height of the droplet, $2h$ is the final height of the droplet, $F$ is the applied squeezing force, $t$ is the squeezing time, $\eta$ is the fluid viscosity, $W$ is the width of the assumed rectangular plate, and $2L_{0}$ is the initial length of the droplet.

Based on conservation of mass calculations, the droplet width is inversely proportional to droplet height; as the width increases, the height decreases in response to squeezing forces.

Circular plate squeeze flow diagram at initial and follow-on conditions; circular plates (gray), droplet (hash marked ellipse between both plates).

Circular plate

For circular plate squeeze flow:

${\frac {h_{0}}{h}}=\left(1+{\frac {16*F*t*h_{0}^{2}}{3*\pi *\eta *R^{4}}}\right)^{1/2}$

$R$ is the radius of the circular plate.

Rectangular plate

Rectangular plate squeeze flow diagram at initial and follow-on conditions; rectangular plates (gray), droplet (green).

For rectangular plate squeeze flow:

${\frac {h_{0}}{h}}=\left(1+{\frac {F*t*h_{0}^{2}}{2*\mu *W*L^{3}}}\right)^{1/2}$

These calculations assume a melt layer that has a length much larger than the sample width and thickness.

Non-Newtonian fluids

Simplifying calculations for Newtonian fluids allows for basic analysis of squeeze flow, but many polymers can exhibit properties of non-Newtonian fluids, such as viscoelastic characteristics, under deformation. The power law fluid model is sufficient to describe behaviors above the melting temperature for semicrystalline thermoplastics or the glass transition temperature for amorphous thermoplastics, and the Bingham fluid model provides calculations based on variations in yield stress calculations.

Power law fluid

For squeeze flow in a power law fluid:

${\frac {h_{0}}{h}}=\left(1+t*({\frac {2n+3}{4n+2}})({\frac {(4*h_{0}*L_{0})^{n+1}*F*(n+2)}{(2*L_{0})^{2n+3}*W*m}})^{1/n}\right)^{n/2n+3}$

Where $m$ (or $K$) is the flow consistency index and $n$ is the dimensionless flow behavior index.

$m=m_{0}*exp\left({\frac {-E_{a}}{R*T}}\right)$

Where $m$ is the flow consistency index, $m_{0}$is the initial flow consistency index, $E_{a}$is the activation energy, $R$is the universal gas constant, and $T$is the absolute temperature.

During experimentation to determine the accuracy of the power law fluid model, observations showed that modeling slow squeeze flow generated inaccurate power law constants ($m$ and $n$) using a standard viscometer, and fast squeeze flow demonstrated that polymers may exhibit better lubrication than current constitutive models will predict. The current empirical model for power law fluids is relatively accurate for modeling inelastic flows, but certain kinematic flow assumptions and incomplete understanding of polymeric lubrication properties tend to provide inaccurate modeling of power law fluids.

Bingham fluid

Bingham fluids exhibit uncommon characteristics during squeeze flow. While undergoing compression, Bingham fluids should fail to move and act as a solid until achieving a yield stress; however, as the parallel plates move closer together, the fluid shows some radial movement. One study proposes a “biviscosity” model where the Bingham fluid retains some unyielded regions that maintain solid-like properties, while other regions yield and allow for some compression and outward movement.

$\tau ={\begin{cases}\eta _{2}*{du \over dy}+\tau _{1},&{\text{if }}\tau \geq \tau _{1}\\\eta _{1}*{du \over dy},&{\text{if }}\tau <\tau _{1}\end{cases}}$

Where $\eta _{2}$ is the known viscosity of the Bingham fluid, $\eta _{1}$ is the "paradoxical" viscosity of the solid-like state, and $\tau _{1}$ is the biviscosity region stress. To determine this new stress:

$\tau _{0}=\tau _{1}(1-\epsilon )$

Where $\tau _{0}$ is the yield stress and $\epsilon ={\frac {\eta _{2}}{\eta _{1}}}$ is the dimensionless viscosity ratio. If $\epsilon =1$, the fluid exhibits Newtonian behavior; as $\epsilon \rightarrow 0$, the Bingham model applies.

Applications

Squeeze flow application is prevalent in several science and engineering fields. Modeling and experimentation assist with understanding the complexities of squeeze flow during processes such as rheological testing, hot plate welding, and composite material joining.

Rheological testing

Squeeze flow rheometry allows for evaluation of polymers under wide ranges of temperatures, shear rates, and flow indexes. Parallel plate plastometers provide analysis for high viscosity materials such as rubber and glass, cure times for epoxy resins, and fiber-filled suspension flows. While viscometers provide useful results for squeeze flow measurements, testing conditions such as applied rotation rates, material composition, and fluid flow behaviors under shear may require the use of rheometers or other novel setups to obtain accurate data.

Hot plate welding; (a) Heating Phase, (b) Change-over Phase, (c) Joining Phase. Polymer materials (gray), heating element (red), melt layers (orange), squeeze out regions (orange, hash mark borders).

Hot plate welding

During conventional hot plate welding, a successful joining phase depends on proper maintenance of squeeze flow to ensure that pressure and temperature create an ideal weld. Excessive pressure causes squeeze out of valuable material and weakens the bond due to fiber realignment in the melt layer, while failure to allow cooling to room temperature creates weak, brittle welds that crack or break completely during use.

Composite material joining

Prevalent in the aerospace and automotive industries, composites serve as expensive, yet mechanically strong, materials in the construction of several types of aircraft and vehicles. While aircraft parts are typically composed of thermosetting polymers, thermoplastics may become an analog to permit increased manufacturing of these stronger materials through their melting abilities and relatively inexpensive raw materials. Characterization and testing of thermoplastic composites experiencing squeeze flow allow for study of fiber orientations within the melt and final products to determine weld strength. Fiber strand length and size show significant effects on material strength, and squeeze flow causes fibers to orient along the load direction while being perpendicular to the joining direction to achieve the same final properties as thermosetting composites.