Mathematical Diffusion Modeling of Saturated Hydrocarbon Oligomers in Multi-Layer Packaging

Flux

Polyolefin packaging releases saturated hydrocarbon oligomers into food contact matrices through concentration-driven diffusion. These migrants ~ linear and cyclic alkanes ranging from ten to thirty-five carbon units ~ stem from incomplete polymerization, residual wax processing aids, and Ziegler-Natta catalyst carrier diluents. Under Regulation (EU) 10/2011, validated numerical modeling serves as an accepted alternative to direct migration testing for compliance work.

Calculating mass transfer requires solving unsteady-state partial differential equations under defined contact temperatures and layer geometries, using Fickian second-law models across mono- and co-extruded films.

Planar geometry dictates the governing rate of one-dimensional mass transport within an isotropic polymer phase:

∂C/∂t = D × (∂²C/∂x²)

Variable C denotes oligomer concentration at position x across film thickness L at contact elapsed duration t. Parameter D represents the diffusion coefficient of the saturated hydrocarbon species within the host polymer. Molecular weight dictates the relative velocity of individual hydrocarbon chains through the tangled polymer network.

Low molecular weight oligomers beneath 300 Daltons migrate rapidly, whereas species exceeding 600 Daltons encounter severe steric hindrance from polymer crystal lamellae.

Saturated hydrocarbon oligomers beneath three hundred Daltons penetrate polyolefin networks during early ambient storage.

Polymer morphology sets migration rates through tortuosity and free volume distribution. The branched chains in low-density polyethylene yield higher diffusion coefficients than the tighter structures of high-density polyethylene. For regulatory submissions, the Piringer model offers a standardized, conservative estimate of diffusion coefficients:

D = D₀ × exp(Ap – 0.55 × M_r^(2/3) – 0.059 × M_r^(1/3) + 10454 / T_ref – 10454 / T)

Parameter Ap describes the polymer-specific diffusion parameter, while tau represents the matrix activation parameter. Parameter M_r defines migrant relative molecular mass, and T defines absolute thermodynamic contact temperature in Kelvin. Regulatory authorities accept Ap values that overestimate measured physical transfer to preserve safety margins.

Polymer chain mobility and backbone architecture directly determine migrant diffusion rates through the film.

Layer

Flexible packaging commonly uses co-extruded or laminated multi-layer films to combine sealability, mechanical stiffness, and barrier performance. A standard structure might pair a low-density polyethylene sealant, a maleic anhydride-grafted tie layer, an ethylene vinyl alcohol barrier core, and an oriented polyamide outer ply. Hydrocarbon oligomers in outer plies migrate inward unless intercepted by dense intermediate layers.

Modeling migration through multi-layer laminates requires enforcing continuous mass flux and thermodynamic partitioning at every internal boundary.

May Multi-Layer Functional Barriers Halt Oligomer Transfer?

Ethylene vinyl alcohol and metallized polyolefin sheets act as functional barrier plies against hydrocarbon transmission. The boundary condition linking ply i and ply i+1 at interface coordinate x_int satisfies continuous mass conservation:

D_i × (∂C_i/∂x) = D_(i+1) × (∂C_(i+1)/∂x)

Interfacial concentrations demonstrate a thermodynamic step discontinuity governed by the polymer-to-polymer partition coefficient K_(i,i+1):

C_i(x_int, t) = K_(i,i+1) × C_(i+1)(x_int, t)

A partition coefficient deviating from unity shifts the equilibrium concentration gradient across the laminate joint. Saturated hydrocarbon oligomers partition preferentially into polyolefin phases over polar polyamide or ethylene vinyl alcohol cores. Polar polymers exhibit low equilibrium solubility for non-polar aliphatic hydrocarbon chains, suppressing transient migration into the core.

Gaps or pinholes in the core layer allow oligomers to bypass the barrier entirely. Adhesive lamination lines applied below nominal thickness targets leave micro-channels across tie boundaries, while conversion temperatures alter regional crystallinity enough to increase diffusion coefficients by an order of magnitude.

Thin barrier skins do not automatically remove the requirement for batch-specific extraction testing under strict audit conditions.

Grid

Analytical solutions to transient multi-layer diffusion equations are limited to infinite boundary assumptions and single-layer structures. Practical compliance assessments use numerical discretization across discrete spatial and temporal points, dividing the cross-section into non-uniform meshes with refined node spacing at polymer interfaces.

Will Finite Difference Meshes Overestimate Interfacial Flux?

Explicit Euler methods risk numeric instability whenever time increments exceed thresholds set by the spatial grid spacing. Crank-Nicolson implicit discretization avoids this stability limit by averaging spatial derivatives between the current and subsequent time steps across each internal node j in ply m.

Discretization equates the temporal concentration change to the averaged second-order central spatial differences across consecutive time intervals.

Constructing tridiagonal matrices allows rapid inversion along every time step n. Interfacial boundary nodes incorporate flux conservation equations through three-point asymmetric central differences. When spatial increments Δx drop below 0.1 microns, numerical truncation errors vanish entirely.

A ten-day exposure at forty degrees Celsius in vegetable oil simulant anchors standard European compliance dossiers for long-term shelf storage.

Initial migrant profiles reflect conversion and storage history. Polyolefin rolls stored under tension undergo contact transfer before pouch forming or filling, as printed outer plies transfer oligomer fractions to the inner food-contact surface through roll set-off. Numerical models account for this by setting non-zero initial boundary concentrations C(x, 0) across the sealant layer.

Finite difference modeling accounts for specific laminate structural configurations:

• Total migrant mass balance tracks initial compound load across combined polymer layers against cumulative release into the contacting medium to eliminate non-physical mass generation.
• Partition step corrections maintain equilibrium concentration ratios across adjacent polymer nodes at each advancing time increment.
• Dynamic food boundary constraints update migrant solubility limits within finite food volumes as migrant saturation conditions develop.
• Solvent swelling factors adjust polymer diffusion parameters when aggressive fatty food simulants plasticize polyolefin sealant layers during contact.

Laminators meeting verified numerical modeling protocols avoid repeating costly solvent extraction runs across every film conversion batch under Regulation (EU) 10/2011 Annex XVI provisions.

Simulant

Validating numerical migration models requires benchmarking predicted transfer curves against standardized laboratory extractions. Regulation (EU) 10/2011 establishes specific food simulants to replicate real food matrices: Simulant D1 (50 percent aqueous ethanol) represents dairy products and oil-in-water emulsions, Simulant D2 (refined vegetable oil) covers fatty food contact, and Simulant E (poly(2,6-diphenyl-p-phenylene oxide), or Tenax) simulates dry food contact.

Fatty simulants swell polyolefins by penetrating the amorphous regions. Vegetable oil sorption into low-density polyethylene lowers the matrix glass transition temperature and accelerates oligomer diffusion. The equilibrium food-to-polymer partition coefficient K_fp governs this distribution:

K_fp = C_polymer,eq / C_food,eq

For lipophilic saturated hydrocarbons in contact with vegetable oil, K_fp values range between 0.1 and 1.0, promoting near-complete extraction of accessible migrants. In contrast, aqueous media such as Simulant A yield K_fp values above 1000, holding the hydrocarbons inside the polymer network.

Lipophilic saturated hydrocarbon oligomers transfer quantitatively into vegetable oil simulants while remaining largely immobilized against aqueous test media.

Consider a worked compliance verification calculation for a three-layer co-extruded film: 50 microns of low-density polyethylene sealant containing 1200 mg/kg initial C18-C24 saturated hydrocarbon oligomers, 15 microns of ethylene vinyl alcohol barrier, and 20 microns of oriented polypropylene backing. The film contacts Simulant D2 with a surface-to-volume ratio of 6 dm² per 1 kg of food under OM2 conditions (10 days at 40 °C). Numerical integration calculates an initial rapid release from the inner 12 microns of the sealant layer, delivering 1.8 mg/kg of oligomers into the simulant within 48 hours.

The middle ethylene vinyl alcohol barrier halts outward migration from the outer polypropylene layer, restricting total ten-day migration to 2.4 mg/kg against a calculated infinite single-ply migration scenario of 8.9 mg/kg.

Submitting an unverified model that fails laboratory audit forces immediate customs impoundment and product withdrawal at the port of entry.

Dispute

Enforcement authorities and packaging converters continue to dispute toxicological thresholds and chromatographic classification for saturated hydrocarbons. Testing protocols group saturated fractions into Mineral Oil Saturated Hydrocarbons (MOSH) and Polyolefin Oligomeric Saturated Hydrocarbons (POSH). On gas chromatography with flame ionization detection (GC-FID), both fractions form overlapping unresolved complex humps, making precise speciation difficult.

Regulatory authorities scrutinize migration dossiers based on whether modeling parameters reflect realistic or worst-case physics:

1. Piringer parameter selection generates debate because generic Ap coefficients inflate estimated migration rates by a factor of two to five above experimental values.
2. Barrier activation energies vary between fresh extruded films and converted pouches subjected to flex cracking during transport.
3. Solvent absorption corrections remain unstandardized across commercial modeling platforms when calculating fatty food exposure.
4. Toxicological threshold allocation divides regulators who enforce a 0.05 mg/kg limit for unlisted oligomers against converters applying 10 mg/kg general thresholds.

Converters favor using experimentally measured diffusion coefficients tailored to specific resin grades and run conditions. Regulators, however, routinely reject proprietary constants unless verified under standardized EN 13130 protocols. High-density polyethylene masterbatches formulated with synthetic waxes draw particular scrutiny in fatty food applications.

Separating synthetic polyolefin oligomers from recycled paperboard mineral oils poses severe analytical challenges. While offline liquid chromatography separates aliphatics from aromatics, branched polyolefin oligomers elute in the same retention windows as natural isoprenoid hydrocarbons, leaving residual uncertainty during cross-laboratory audits.

Scientific consensus remains absent regarding whether mathematical diffusion equations can fully account for polymer aging and structural relaxation during multi-year storage cycles.