Capillary Rheometry Shear Rate Correction Protocols for Bimodal High Density Polyethylene
Correcting capillary rheometry data via Bagley and Rabinowitsch protocols is mandatory to prevent up to 50 percent viscosity errors in bimodal HDPE die design.
Barrel
Capillary testing of dual-peak polyethylene reveals severe melt distortions when apparent values pass unadjusted through process design calculations. High-density polyethylene synthesized with a bimodal molecular weight distribution incorporates two distinct polymer populations within a single resin pellet. The low molecular weight component provides melt fluidity during processing, while the ultra-high molecular weight fraction supplies environmental stress crack resistance and long-term hydrostatic strength in blown film, blow molding, and pressure pipe applications.
Standard laboratory capillary rheometers measure force and volumetric flow through a constrained orifice to calculate wall shear stress and apparent wall shear rate. Raw values calculated directly from piston speed and capillary die geometry systematically overestimate true viscosity across processing shear rates ranging from 100 to 10,000 inverse seconds.
High Pressure Extrusion Dynamics of Bimodal Resins
Combining low and high molecular weight fractions creates a distinct dual-regime flow in narrow channels. At low deformation rates, long chains in the high molecular weight tail remain entangled, dominating elastic response and generating significant normal stress differences. As shear rates elevate inside the rheometer channel, these elongated molecules align parallel to flow vectors, triggering a sudden shear-thinning transition.
The low molecular weight fraction acts as an internal lubricant, shifting non-Newtonian behavior to lower shear rates than in unimodal polyethylenes with an equivalent average melt flow rate. Pressure transducers mounted along the rheometer tube record non-linear pressure gradients when high molecular weight entanglements resist entrance deformation at the die convergence zone.
Discrepancies between Melt Flow Rate and Capillary Rheology
Standard single-point index values measured under static mass loads mask the structural shear-thinning response required for high-throughput die design. ASTM D1238 and ISO 1133 standard measurements at 190 degrees Celsius using a 2.16 kilogram load yield identical nominal melt flow rates for resins with radically different molecular architectures. A bimodal pipe grade PE100 resin displaying a melt flow rate of 0.25 grams per ten minutes can exhibit an apparent shear viscosity forty percent lower than a unimodal grade of the same melt index when evaluated at a shear rate of 1,000 inverse seconds.
Bimodal high density polyethylene exhibits dual-stage shear thinning that causes uncorrected melt viscosity data to deviate by up to fifty percent from true in-cavity rheology.
Process simulation software fed with uncorrected capillary rheology underpredicts extrusion head pressures, causing motor overload failures and inaccurate die sizing during tooling fabrication. Tooling adjustments based on raw test figures produce severe wall thickness variations and unexpected melt fracture in high-speed blow molding lines.
Orifice
Pressure losses occurring at the inlet and outlet of a capillary die corrupt apparent shear stress values calculated from total extrusion pressure. As the bimodal melt transitions from the large reservoir into the narrow capillary bore, streamlines converge rapidly, generating severe extensional deformation. This entrance flow stores elastic energy in ultra-high molecular weight polymer chains, creating an extra pressure drop that does not contribute to overcoming fully developed viscous resistance along the capillary length.
Exit pressure drops caused by elastic recovery upon discharge further inflate total measured differential pressure.
The Bagley End Correction Mathematical Formulation
Quantifying entrance and exit pressure drops demands testing the resin through multiple capillaries of identical diameter but varying capillary length-to-diameter ratios. ISO 11443 and ASTM D3835 mandate using at least three dies, typically with length-to-diameter ratios of 10:1, 20:1, and 30:1, alongside a zero-length orifice die possessing a length-to-diameter ratio below 0.2. Total extrusion pressure is plotted against the capillary length-to-diameter ratio at constant apparent shear rates.
Linear extrapolation of these Bagley lines to a hypothetical length-to-diameter ratio of zero yields the end-correction pressure loss, or alternatively, the dimensionless entrance correction factor e. True shear stress at the capillary wall is derived by subtracting the end pressure drop from total measured pressure, or by dividing total pressure by four times the sum of the actual length-to-diameter ratio and the correction factor e.
- Mount a tungsten carbide capillary die with a length-to-diameter ratio of 10:1 and a diameter of 1.0 millimeter into the heated rheometer barrel maintained at 190 degrees Celsius.
- Pack bimodal high density polyethylene pellets into the barrel in four increments, tamping between charges to purge air pockets and preheating the charge for six minutes.
- Execute a programmed speed profile generating apparent shear rates at 50, 100, 250, 500, 1000, 2500, and 5000 inverse seconds while recording steady-state barrel pressures.
- Clean the barrel thoroughly and repeat the preheat and piston speed schedule using dies with length-to-diameter ratios of 20:1 and 30:1 at identical test temperatures.
- Plot measured total pressure against length-to-diameter ratio for each shear rate, performing linear regression to establish the zero-length intercept pressure.
- Calculate true wall shear stress by applying the derived zero-length pressure drop to the raw capillary pressure data for every target shear rate setting.
Non-Linear Bagley Response in Bimodal Systems
Extrapolating pressure data across broad shear rate bands frequently reveals non-linear Bagley plots in resins containing ultra-high molecular weight fractions. At high piston speeds, hydrostatic pressure within the rheometer barrel reaches levels exceeding 50 megapascals, increasing local melt viscosity and causing the pressure gradient along the capillary to curve upward. Entrance pressure losses elevate non-linearly due to flow-induced crystallization or severe extensional strain hardening of long-chain entanglements.
Applying a linear fit to upward-curving Bagley data introduces systematic errors, underestimating entrance losses at elevated shear rates and overestimating true shear stress.
An ISO 11443 Bagley correction performed at 190 degrees Celsius on a bimodal pipe resin yields end-correction factors between 4.2 and 8.7 die diameters across commercial extrusion shear rates.
| Polymer Grade Architecture | Apparent Shear Rate (1/s) | Total Pressure L/D 30 (MPa) | Zero-Length Pressure Drop (MPa) | Bagley Correction Factor e | True Wall Shear Stress (kPa) |
|---|---|---|---|---|---|
| Unimodal Film Grade (MFR 0.05) | 100 | 18.4 | 1.8 | 1.47 | 138.3 |
| Unimodal Film Grade (MFR 0.05) | 1,000 | 38.2 | 4.2 | 1.65 | 283.3 |
| Bimodal Pipe Grade PE100 | 100 | 22.1 | 3.6 | 2.44 | 154.2 |
| Bimodal Pipe Grade PE100 | 1,000 | 46.8 | 10.8 | 3.46 | 300.0 |
| Bimodal Pipe Grade PE100 | 5,000 | 74.5 | 22.4 | 4.51 | 434.2 |
| Bimodal Blow Molding Grade | 1,000 | 42.1 | 8.9 | 3.17 | 276.7 |
| Data measured per ISO 11443 using 1.0 mm diameter capillary dies with entrance angle of 90 degrees under nitrogen purge. | |||||
Extrusion flow instabilities are routinely attributed to raw material contamination when off-spec shear viscosity charts are derived without multi-die entrance corrections.
Gradient
Calculating shear rates inside a capillary die requires correcting for the non-Newtonian velocity profile of the polymer melt. The apparent wall shear rate formula assumes a parabolic velocity profile characteristic of Newtonian fluids, where fluid velocity scales quadratically from zero at the wall to a maximum at the channel centerline. High density polyethylene exhibits pronounced non-Newtonian shear-thinning, causing the velocity profile to flatten near the center and develop steep velocity gradients adjacent to the die wall.
The actual shear rate experienced by the melt at the wall exceeds the apparent shear rate calculated from volumetric throughput.
Weissenberg-Rabinowitsch Correction Mechanics
True wall shear rate determination relies on the Weissenberg-Rabinowitsch protocol, which applies a mathematical correction factor based on the local derivative of logarithmic shear stress with respect to logarithmic apparent shear rate. This derivative represents the power-law index slope n prime at the specific stress level. True wall shear rate equals the apparent wall shear rate multiplied by a factor of three-fourths plus one-fourth of the inverse power-law slope.
For a Newtonian fluid, the slope n prime equals one, reducing the factor to unity. For strongly shear-thinning bimodal resins, n prime drops to values between 0.25 and 0.40, increasing the true wall shear rate by thirty to fifty percent above apparent values.
Are Inflection Points Diagnostic for Bimodal Polymer Tails?
Plotting the local power-law slope n prime against apparent shear rate for bimodal high density polyethylene yields a distinct non-monotonic curve containing two inflection points. The first minimum in n prime occurs at moderate shear rates between 100 and 500 inverse seconds, corresponding to the disentanglement transition of the ultra-high molecular weight tail. The second regime develops at shear rates above 2,000 inverse seconds as the lower molecular weight matrix reaches maximum shear thinning.
Unimodal resins display a smooth, continuous decrease in n prime without intermediate inflection points. Tracking these slope transitions confirms the physical presence of low-frequency high-molecular-weight fractions without performing gel permeation chromatography.
True wall shear rates for bimodal high density polyethylene surpass apparent shear rates by up to forty-five percent due to extreme velocity profile flattening near the die channel boundary.
Accurate curve fitting of the double logarithmic shear stress versus apparent shear rate relationship is vital for precise derivative calculation. Standard polynomial fits often introduce mathematical oscillations that distort the local slope n prime, generating artificial viscosity spikes. Fitting data with an adjusted low-order Cross model or local sliding-window polynomial smoothing prevents mathematical artifacts while capturing genuine structural inflections.
| Apparent Shear Rate (1/s) | True Shear Stress (kPa) | Local Power-Law Slope (n’) | Rabinowitsch Correction Factor | True Shear Rate (1/s) | True Dynamic Viscosity (Pa·s) |
|---|---|---|---|---|---|
| 50 | 112.0 | 0.52 | 1.231 | 61.5 | 1,821.1 |
| 100 | 154.2 | 0.41 | 1.360 | 136.0 | 1,133.8 |
| 250 | 208.5 | 0.33 | 1.508 | 376.9 | 553.2 |
| 500 | 252.1 | 0.29 | 1.612 | 806.0 | 312.8 |
| 1,000 | 300.0 | 0.31 | 1.556 | 1,556.0 | 192.8 |
| 2,500 | 368.4 | 0.26 | 1.712 | 4,278.8 | 86.1 |
| 5,000 | 434.2 | 0.24 | 1.792 | 8,958.3 | 48.5 |
- Polynomial Over-Fitting creates false oscillations in the power-law derivative, generating artificial steps in corrected viscosity curves.
- Truncated Data Ranges obscure the low-shear relaxation spectrum, leading to inaccurate zero-shear viscosity extrapolations.
- Uncorrected Temperature Drifts reduce local melt viscosity during high-shear runs, artificially steepening the power-law slope.
- Capillary Aspect Ratios below 20:1 amplify entrance turbulence, distorting the local shear stress measurements required for derivative estimation.
The structural relationship between the double-peak relaxation spectrum and local power-law slope inflections remains an active field of investigation among polymer physicists attempting to predict melt fracture thresholds directly from molecular weight distributions.
Slippage
Standard continuum fluid mechanics assumes zero fluid velocity at the solid boundary wall. High molecular weight polyethylenes violate this boundary condition when local wall shear stresses exceed critical thresholds, typically between 0.10 and 0.18 megapascals. At these stress levels, ultra-high molecular weight chains adsorbed onto the metallic die surface disentangle abruptly from the bulk melt, causing the fluid to slide along the channel boundary.
Mooney Slip Correction Protocol
Quantifying genuine melt viscosity in the presence of boundary slip requires applying the Mooney correction method. Capillary rheometry testing is performed across multiple die diameters while maintaining constant capillary length-to-diameter ratios, typically using dies with internal diameters of 0.5, 1.0, and 2.0 millimeters. Apparent wall shear rate is plotted against inverse die radius at constant true wall shear stress levels.
The slope of the resulting linear plots yields slip velocity, while the vertical intercept defines the true bulk wall shear rate corrected for slip effects.
Distinguishing wall slip from volume shear-thinning requires inspecting the extrusate surface appearance and monitoring pressure transducer stability. A visual identification protocol guides compounding lines through instability classification:
- Surface Sharkskin manifests as high-frequency, low-amplitude periodic roughness caused by local extensional tensile failure of the melt surface skin as it exits the die die land.
- Stick-Slip Instability displays periodic pressure oscillations and alternating bands of smooth and rough extrusate resulting from repetitive cycle transitions between slip and non-slip boundary conditions.
- Gross Melt Fracture exhibits chaotic, non-periodic extrusate distortion occurring when bulk elastic strain energy exceeds the cohesive strength of the polymer melt melt matrix.
- Adhesive Wall Slip produces a sudden decrease in measured apparent viscosity without visual extrusate distortion, verified by die diameter dependent flow curves.
Eliminating boundary slip artifacts requires applying fluoroelastomer processing aids during rheological characterization or restricting test stress levels below the slip activation threshold. Dynamic mechanical analysis in the linear viscoelastic regime serves as a non-destructive alternative for characterizing base resin viscosity without slip interference.
Capillary dies with ground internal surface finishes below 0.1 micrometers Ra prevent mechanical anchor points that artificially delay slip onset.
Invoice
Procuring commercial bimodal high density polyethylene requires translating rheological corrections into enforceable quality specifications and landed cost calculations. Polymer producers supply resin data sheets carrying single-point melt flow index values that fail to capture batch-to-batch variations in molecular weight distribution broadness. A shift in the ratio between low and high molecular weight reactor steps alters high-shear processability while maintaining static melt index numbers.
Converting incoming raw material specifications to multi-point capillary viscosity protocols protects converter margins against unexpected scrap rate spikes during high-speed extrusion processing.
Commercial Specification Framework for Bimodal Pipe Grades
Purchase contracts for PE100 pipe compounds must incorporate corrected capillary viscosity limits alongside standard density and melt flow rate windows. Raw material acceptance standards mandate providing Bagley-corrected and Rabinowitsch-corrected viscosity curves at 190 degrees Celsius across shear rates from 100 to 5,000 inverse seconds. A commercial batch tolerance window limits true shear viscosity variations to within plus or minus five percent of the target reference curve across the entire shear rate spectrum.
| Rheological Parameter | Target Specification | Test Standard & Conditions | Deviation Limit | Processing Impact | Landed Cost Penalty |
|---|---|---|---|---|---|
| True Viscosity (100 s⁻¹) | 1,130 Pa·s | ISO 11443, 190 °C, L/D 30:1 | ±5.0% | Pipe wall sag and parison drop | $45/tonne sorting fee |
| True Viscosity (1,000 s⁻¹) | 192 Pa·s | ISO 11443, 190 °C, L/D 30:1 | ±5.0% | Extruder head pressure spikes | $120/tonne regrind discount |
| Bagley Correction e (1,000 s⁻¹) | 3.45 | ISO 11443 multi-die matrix | ±10.0% | Die entrance turbulence and swell | Rejection of railcar lot |
| Stick-Slip Stress Threshold | > 0.14 MPa | ASTM D3835 pressure scan | Minimum value | Premature surface sharkskin onset | $80/tonne processing aid cost |
Extrusion processors purchasing resin on spot markets face substantial commercial risks when raw material suppliers substitute unimodal or broad-distribution grades under generic bimodal trade names.
A two percent shift in true dynamic viscosity at 1,000 inverse seconds alters extrusate swell by four percent, driving pipe wall thickness outside dimensional tolerance limits.
A worked financial scenario demonstrates the financial consequence of processing unverified bimodal resin lots. Assume a pipe extrusion facility processes 40,000 metric tonnes of PE100 resin annually at a delivered pellet cost of $1,450 per tonne. Processing a resin batch with an uncorrected viscosity discrepancy that reduces high-shear flowability causes a three percent drop in line speed and increases wall thickness variation scrap from 1.5% to 4.2%.
Over a 1,000-tonne production run, the excess scrap generates 27 tonnes of off-spec pipe requiring reprocessing or re-sale at a $500 per tonne downgrade penalty, incurring direct losses of $13,500 per lot alongside $18,000 in lost machine throughput value.
Standard raw material purchase contracts incorporate quality assurance clause 14.3, specifying that resin deliveries failing to meet true Rabinowitsch-corrected viscosity profiles within specified tolerance bands trigger immediate seller indemnification for processing scrap and line downtime penalties.