Anisotropic Extrudate Swell from a Slit Die: A Velocity-Centre Hypothesis and Numerical Verification

TL;DR

This paper introduces a new theoretical model to understand and predict anisotropic extrudate swell in polymer processing using a velocity-centre hypothesis and numerical simulations.

Contribution

The paper proposes a novel semi-empirical framework to isolate geometric and kinematic effects in extrudate swell using a multi-velocity-centre hypothesis.

Findings

Analytical predictions align well with CFD data (error < 5%) in the core zone of high-aspect-ratio dies.

Velocity decay near die edges causes local deviations due to 3D flow kinematics.

The model provides a mathematical basis for inverse die-profile design and shape compensation.

Abstract

While anisotropic extrudate swell in polymer processing is fundamentally driven by physical viscoelastic recovery, this paper proposes a theoretical framework to explicitly isolate and map the purely geometric and kinematic components of this phenomenon. Serving as a mathematical proof-of-concept, a multi-velocity-centre hypothesis is proposed. By introducing a semi-empirical, lumped material-flow calibration parameter, the macroscopic diameter swell ratio is mathematically extended to the discrete local flow field of a rectangular slit die. To evaluate its validity, the analytical framework is subjected to a numerical test for kinematic consistency utilizing isothermal, inelastic power-law fluid CFD simulations, thereby separating geometric mapping from complex viscoelastic stress relaxation. Results indicate that analytical predictions show good agreement with CFD data (error < 5%) strictly within the core zone of high-aspect-ratio dies. However, due to the infinite-slit assumption, 3D flow kinematics near die edges induce velocity decay, leading to local deviations that require future empirical corrections. Although comprehensive physical extrusion experiments and non-isothermal viscoelastic coupling are required for industrial deployment, this semi-empirical kinematic mapping provides a foundational mathematical basis that could potentially inform future inverse die-profile design and shape distortion compensation.

Figures

Full text

1. Introduction

Polymer extrusion is a fundamental manufacturing process for polymeric products. Dimensional accuracy and geometric stability of the extrudate are critically influenced by extrudate swell (the Barus effect) upon exiting the die [1,2]. This phenomenon results from the release of viscoelastic stresses accumulated in the melt during flow through the confined die channel. For axisymmetric flows (e.g., through circular dies), extrudate swell is relatively uniform, and well-established predictive models exist [2,3]. In contrast, for profile extrusion through slit dies, the swell exhibits pronounced anisotropy, with systematic directional variations in the swell ratio [4,5]. However, in profile extrusion—particularly through slit dies—the swell manifests pronounced anisotropy, characterized by systematic variations in the swell ratio along different directions [6,7,8]. This non-uniform swelling critically distorts the cross-sectional profile of extrudates, posing a fundamental challenge for precision die design and product quality control [9,10,11,12]. Consequently, elucidating physical origins of anisotropic swell and establishing accurate predictive frameworks remain persistent and industrially significant scientific pursuits in polymer processing.

The phenomenon of polymer extrudate swell significantly influences extrusion die design [13,14,15,16]. Significant scholarly attention has focused on three core aspects: (1) theoretical prediction of extrudate swell ratios, (2) techniques for tracking and capturing free surface evolution during swelling [17,18], and (3) characterization of non-uniform swell behaviour. Tanner [1] established foundational formulas quantifying swell ratios for both circular and slit die geometries. Pearson et al. [5] provided comparative validation of Tanner’s formulations, extending their application to long- and short-slit dies to derive generalized computational frameworks. Tang et al. [6,7] established a tripartite methodology—encompassing rheological theory, three-dimensional viscoelastic flow simulations, and digital holography experiments—to comprehensively characterize swelling mechanisms in slit die extrudates. Researchers have conducted complementary theoretical and experimental investigations. These collective advances establish an effective methodology for swell prediction in axisymmetric configurations. Nevertheless, theoretical determination of extrudate swelling remains challenging for geometrically complex dies, necessitating the adoption of numerical simulations as the prevailing predictive approach.

The underlying mechanisms governing anisotropic extrudate swell in viscoelastic polymer flows remain incompletely characterized. While extensive studies employing theoretical frameworks, experimental techniques, and computational simulations have addressed the general phenomenon of extrudate swell, the specific origins of its non-uniform manifestations remain elusive to rigorous theoretical elucidation. This knowledge gap persists despite the established significance of anisotropic swelling in the fidelity of extrusion die design and dimensional accuracy of the product. Consequently, a systematic investigation is warranted to: (1) establish quantitative descriptors for anisotropic swell phenomena across slit dies; (2) develop predictive models for slit-flow specific anisotropy; and (3) explore the geometric and kinematic mapping rules that link local flow fields to directional swell heterogeneity.

While classical theories, notably Tanner’s elastic recovery framework [1,5], have provided foundational predictions for axisymmetric extrudate swell, they fundamentally attribute swelling to the release of stored elastic energy upon stress cessation—a mechanism implicitly assuming isotropic material response and stress-equilibrated flow at the die exit. This perspective, though successful for capillary geometries, encounters significant limitations in explaining anisotropic swelling in slit dies, where the swelling ratio varies by over 30% in different directions [6]. The long-standing dichotomy between elastic recovery and viscous redistribution mechanisms remains unresolved in such confined, non-axisymmetric flows. Proponents of the elastic recovery paradigm argue that normal-stress-driven recoil dominates swelling [1,19], whereas advocates of viscous redistribution emphasize the role of streamline reorientation and shear-induced velocity rearrangement at the free surface [7,20]. While Tanner’s stress-release approach provides a macro-scale thermodynamic perspective, the present velocity-centre hypothesis complements this by focusing on localized flow kinematics. It posits that the macroscopic anisotropic swelling profile can be mathematically correlated with multifocal reorganization of the velocity gradient field immediately downstream of the die lip. This hypothesis aims to link the spatial distribution of velocity cores within the die land to directional swell heterogeneity, thereby isolating the kinematic mapping component from complex coupled viscoelastic recovery.

To address this, the present investigation develops a theoretical framework predicated upon a multi-velocity-centre hypothesis to mathematically describe flow rearrangement processes responsible for non-uniform swelling phenomena. While the anisotropic extrudate swell of real polymers is fundamentally driven by physical viscoelastic recovery, the primary objective of this theoretical study is to explicitly isolate and map purely geometric and kinematic components of this complex phenomenon. Therefore, this paper is intended as a semi-empirical mathematical proof-of-concept for the multi-velocity-centre hypothesis.

2. Extrudate Swell Behaviour

Upon exiting the extrusion die, polymer melt exhibits extrudate swell (the Barus effect). For products with axisymmetric cross-sections (e.g., from circular dies), the flow is essentially one-dimensional, as the velocity profile depends only on radial position. Consequently, the resulting extrudate swell is usually uniform and isotropic. However, for profile extrusion products with non-circular cross-sections (e.g., from slit or complex dies), the swell becomes significantly anisotropic, meaning the degree of swelling varies with direction [1,7]. This anisotropy stems from several factors: (1) geometric irregularities in the flow channel, (2) distortions in the flow field caused by complex cross-sectional geometries, and (3) material nonlinearities such as viscoelastic memory effects. Furthermore, processing parameters like flow rate and temperature gradients can exacerbate anisotropic swelling.

Extrudate swell is quantitatively characterized by the extrudate swell ratio (B), defined as:

$$[eqn]$$

where Se is the stabilized cross-sectional area of extrudate in the free-swelling zone, and S represents the geometric cross-sectional area at the die exit.

For circular die extrusion systems, extrudate swell is alternatively quantified using the diameter swell ratio ( $[eqn]$ ), defined as:

$$[eqn]$$

where De is the stabilized diameter measured at the equilibrium swell state, and D denotes the nominal geometric diameter at the die exit plane.

For axisymmetric circular dies, the extrudate swell ratio (B) reduces to a function of diameter swell ratio ( $[eqn]$ ):

$$[eqn]$$

Building upon Tanner’s elastic recovery theory for circular extrudate swell [1], Kauzlarich [21] derived a low-shear-stress approximation for the relative extrudate swell ratio $[eqn]$ :

$$[eqn]$$

where $[eqn]$ is the Elastic parameter ( $[eqn]$ [s^2^/Pa], material constant), and $[eqn]$ is wall shear stress (Pa).

For fully developed laminar flow of a power-law fluid in a circular tube, the constitutive relationship between wall shear stress ( $[eqn]$ ) and wall shear rate ( $[eqn]$ ) is given by:

$$[eqn]$$

where K is the consistency index, zero-shear viscosity parameter, and n is the flow behaviour index (dimensionless).

Relative extrudate swell ratio $[eqn]$ defined in Equation (4) is derived from Equation (5) through viscoelastic energy recovery theory:

$$[eqn]$$

For steady laminar flow of an incompressible power-law fluid in a circular channel, the dimensionless velocity distribution is given by [9]:

$$[eqn]$$

where $[eqn]$ is local velocity at radial position r, $[eqn]$ is maximum velocity at the centreline, and R is the tube radius.

Wall shear rate ( $[eqn]$ ) for steady laminar flow of an incompressible power-law fluid is derived from the velocity gradient at the wall:

$$[eqn]$$

The relative extrudate swell ratio ( $[eqn]$ ) quantifying elastic recovery energy is derived from fundamental rheological parameters:

$$[eqn]$$

The relative extrudate swell ratio ( $[eqn]$ ) can be expressed through the material velocity constant $[eqn]$ and geometric scaling:

$$[eqn]$$

where $[eqn]$ is the material velocity constant (Pa^2^·s^2n^) which combines elastic recovery ( $[eqn]$ ), viscous dissipation ( $[eqn]$ ), kinetic energy, and power-law index (n).

To bridge fundamental rheological parameters with the geometric scaling of extrudate swell, a lumped material flow constant $[eqn]$ is introduced. This constant links the viscoelastic memory of the polymer melt—quantified by the elastic parameter $[eqn]$ —with its steady-state viscous dissipation, characterized by the power-law consistency index $[eqn]$ and the shear-thinning exponent $[eqn]$ . Furthermore, $[eqn]$ incorporates the kinetic energy of the flow via maximum velocity, thereby coupling the intrinsic material response with specific processing conditions. Physically, $[eqn]$ represents the recoverable viscoelastic energy density per unit flow history that drives the swelling upon stress cessation. Its derivation is based on low-shear-stress approximation of elastic recovery theory [1] and the fully developed velocity profile of a power-law fluid. Consequently, $[eqn]$ serves as a scaling parameter that combines the combined effects of fluid elasticity, shear-dependent viscosity, and flow kinematics into a single constant, enabling the subsequent geometric scaling of the relative swell ratio $[eqn]$ with the characteristic flow dimension $[eqn]$ .

Note that the parameter $[eqn]$ is not a fundamental material property that can be derived solely from standard rheometry. Instead, it serves as a semi-empirical, lumped calibration parameter that implicitly captures coupling between material flow behaviour and local geometry. In practice, $[eqn]$ must be implicitly fitted using a baseline numerical simulation or a single-point experimental measurement.

Equation (10) provides a key scaling relationship for the relative extrudate swell ratio ( $[eqn]$ ) under controlled flow conditions. For axisymmetric flow with constant maximum velocity ( $[eqn]$ ) and low wall shear stress ( $[eqn]$ ), relative extrudate swell ratio $[eqn]$ is inversely proportional to $[eqn]$ [22].

3. Non-Uniform Extrudate Swell Behaviour

The derivation of the multi-velocity-centre hypothesis is based on the following governing assumptions to ensure scientific transparency and mathematical reproducibility:

• (1)The fluid is considered incompressible, ensuring volumetric flow conservation.
• (2)The flow is fully developed within the slit die and isothermal conditions are maintained.
• (3)The local swelling contribution from any discrete velocity element decays radially according to an inverse square distance relationship.

For slit-extruded profiles, both experiments and simulations reveal pronounced anisotropy [7,8], where the swell ratio differs significantly depending on the measurement direction (Figure 1). Unlike circular dies which exhibit point-symmetric flow with a singular velocity maximum at the centreline, slit dies (aspect ratio > 10) develop a planar-symmetric flow. This symmetry, defined by the plane ( $[eqn]$ ) where shear stress $[eqn]$ , leads to anisotropic swelling through differential stress recovery, expressed as:

$$[eqn]$$

where $[eqn]$ and $[eqn]$ are the width and thickness swell ratio, respectively.

To investigate this correlation, a theoretical model for non-uniform behaviour of extrudate swell in slit extrusion is developed, as shown in Figure 1. The X coordinates of points A and B are 0 mm and 20 mm, respectively. The extrudate swell of slit extrusion products is symmetrical about the axis, so half of it can be studied. Leveraging the symmetry condition ( $[eqn]$ ), a computationally efficient half-domain modelling is employed.

Figure 1 illustrates the theoretical model. Line segment AB represents the symmetry plane (y = 0) of the fluid flow within the straight slit die channel. Point X denotes an arbitrary position along this symmetry plane, while point C is located on the die wall surface. The channel height is defined as $[eqn]$ , exhibiting symmetric geometry about the symmetry plane. As illustrated, the radius R here is the distance between points X and C, which follows the geometric relationship:

$$[eqn]$$

where H is the half-channel height, and $[eqn]$ is the angle between line XC and the slit channel’s vertical axis.

At any arbitrary point X along the symmetry plane (y = 0), the diameter swell ratio $[eqn]$ at point C on the extrudate surface is defined by the geometric deformation following extrusion:

$$[eqn]$$

where Re is the distance R after elastic recovery in the extrudate.

Derived from Equation (13), the radial displacement vector magnitude $[eqn]$ resulting from extrusion swell at material point C is expressed as:

$$[eqn]$$

The swell-induced displacement of material point C is resolved into orthogonal components relative to the die coordinate system:

$$[eqn]$$ $$[eqn]$$

where $[eqn]$ is lateral displacement along the die width dimension (X-axis), and $[eqn]$ is the thickness-direction displacement (Y-axis).

Equation (17) resolves lateral displacement through geometric parametrization of Figure 1:

$$[eqn]$$

where $[eqn]$ is the normal projection coordinate of point C on the die centreline (X-axis), and $[eqn]$ identifies the velocity centreline position of material point X.

The derivation follows the kinematic relationship $[eqn]$ per Figure 1, confirming a linear dependence on coordinate separation.

Displacement components $[eqn]$ and $[eqn]$ arise from cumulative contributions of all velocity elements along the centreline. For any material point X on this centreline (spanning $[eqn]$ ), total displacements at observation point C are obtained through integration over the entire domain:

$$[eqn]$$ $$[eqn]$$

where $[eqn]$ , $[eqn]$ : total displacement vectors in the width (X) and thickness (Y) directions at C; $[eqn]$ , $[eqn]$ : boundary coordinates defining the velocity centreline segment (Figure 1); integral operators: integrate the localized kinematic swelling contributions along the symmetry plane.

From Equation (10), defining the radial swell ratio, the displacement integrals in the x- and y-directions are derived as follows:

$$[eqn]$$ $$[eqn]$$

Total displacement in the x- and y-directions due to non-uniform extrudate swell is quantified by the integrals:

$$[eqn]$$ $$[eqn]$$

where $[eqn]$ is the semi-empirical calibration parameter, $[eqn]$ and $[eqn]$ represent the geometric integral functions in the x- and y-directions, respectively.

$$[eqn]$$ $$[eqn]$$

To prevent integration ambiguity, let $[eqn]$ denote the fixed geometric evaluation point on the cross-section, and let $[eqn]$ act strictly as the dummy integration variable along the die width. Distance $[eqn]$ (as depicted in Figure 1) between the flow element and the evaluation point is consistently defined as

$$[eqn]$$

For non-Newtonian fluids, the power-law index (n) ranges from 0 to 1. The total swell ratio in the x-direction at point C is derived from the integral expression governing swell behaviour:

$$[eqn]$$

For power-law fluids with an index of n = 0.5, the total swell ratio in the x-direction at any coordinate $[eqn]$ on the extrudate profile due to anisotropic extrudate swell is defined by integral expression:

$$[eqn]$$

To characterize the geometric dependence of anisotropic extrudate swell, slit channel half-heights (H) of 0.5, 1.0, and 2.0 mm were simulated with a fixed width (W) of 20 mm and a power-law index n = 0.5. Equation (28) governs the predicted x-direction swell ratios. As shown in Figure 2, swell ratios scale with H due to its direct proportionality to $[eqn]$ , while the distribution pattern remains height-invariant. Notably, the displacement profile exhibits antisymmetric behaviour about the slit centreline ( $[eqn]$ ).

The total swell ratio in the y-direction at point C is derived from the integral expression governing swell behaviour:

$$[eqn]$$

For power-law fluids with an index of n = 0.5, the total swell ratio in the y-direction at any coordinate $[eqn]$ on the extrudate profile due to anisotropic extrudate swell is defined by integral expression:

$$[eqn]$$

Function $[eqn]$ attains its maximum at the slit centreline, $[eqn]$ , where its first derivative is zero. At this critical coordinate, $[eqn]$ attains its global maximum value, corresponding to the peak swelling displacement at the slit die centreline. Functional behaviour exhibits strict monotonicity.

The transverse strain evolution along the streamwise coordinate is illustrated in Figure 3. As derived from Equation (30), the y-direction swell profile exhibits characteristic parabolic symmetry with a maximum at x0 = 10 mm, indicating strain localization at the flow centreline. The reduction in expansion toward the edges (x0 = 0 and x0 = 20) is a direct consequence of 3D edge effects and the presence of side-wall constraints. Unlike the central region which approximates 2D planar Poiseuille flow, the corners are subjected to multi-axial shear that modifies local velocity gradients and restricts purely kinematic rearrangement.

The resultant extrudate swelling emerges from coupled swell effects along orthogonal coordinates. Continuum mechanics analysis of slot extrusion phenomena reveals a fundamental mapping relationship between die geometry coordinates $[eqn]$ and extrudate profiles $[eqn]$ :

$$[eqn]$$

where $[eqn]$ and $[eqn]$ denote total displacement tensors and $[eqn]$ represents the strain transformation operator. The invertibility of $[eqn]$ establishes:

• (i)A predictive framework for anisotropic swelling via forward strain-field computation.
• (ii)An inverse rheological design methodology for die lip contour synthesis.

4. Numerical Verification

In this study, the simulated slit die features a specific aspect ratio of 10. Three-dimensional isothermal flow simulations of power-law fluids (K = 1000 Pa·s, n = 0.5) through confined slit dies were conducted using the ANSYS Polyflow 2022 R1 finite-element framework. These isothermal simulations employing an inelastic power-law model were designed to isolate geometric and kinematic aspects of the multi-velocity-centre hypothesis, thereby circumventing the complexities of viscoelastic stress relaxation and thermal effects inherent to real polymer melts. The computational domain (Figure 4) comprised a 20 mm die land ( $[eqn]$ ) and 40 mm free-jet section ( $[eqn]$ ). Slit dimensions were a gap height (2H) of 2 mm and a width (W) of 20 mm, giving an aspect ratio (W/2H) of 10. Leveraging geometric symmetry, the model employed a quarter-domain formulation to optimize computational efficiency. The spatial discretization employed graded elements across key dimensions: the height direction used 10 elements (bias factor 4), the width direction used 40 elements with end refinement (bias 4), the die land was discretized with 30 axial elements refined at the inlet and outlet (bias 6), and the free-jet section used 35 elements with clustering near the exit zone (bias 13).

A mesh sensitivity analysis was performed to ensure grid-independent results. The baseline mesh was iteratively refined and coarsened, varying the global node count by approximately ±5%. Mesh convergence was deemed achieved when variation in the primary extrudate swell ratios between successive mesh generations was less than 1%. The mesh satisfying this criterion was adopted for all subsequent simulations. The final mesh selected for all production runs satisfied this criterion. For the steady-state solver, a stringent convergence criterion required that the normalized residuals for mass, momentum, and constitutive equations all fall below a threshold of 10^−6^, ensuring sufficient solution accuracy for the nonlinear flow fields. These measures collectively confirm that the presented simulation results are quantitatively robust and independent of numerical discretization parameters.

Boundary conditions were applied as follows:

(1)Inlet ( $[eqn]$ ): A fully developed flow profile with a volumetric flow rate of Q = 200 mm^3^/s.(2)Walls ( $[eqn]$ ): No-slip condition ( $[eqn]$ ).(3)Outlet ( $[eqn]$ ): Stress-free boundary ( $[eqn]$ ).(4)Free surfaces ( $[eqn]$ ): Kinematic condition ( $[eqn]$ ) with zero traction force.(5)Symmetry planes ( $[eqn]$ ): Mirror symmetry ( $[eqn]$ , $[eqn]$ ).

Analysis of the numerical results (Figure 5) reveals that geometric symmetry about the x-axis permits evaluation of half-domain extrusion behaviour; along this coordinate, the x-direction swelling displacement—defined as the spatial offset between die exit and extrudate profiles—displays a monotonic increase aligned with the predictive curve, while y-direction displacement fields exhibit parabolic spatial progression consistent with theoretical distribution, confirming non-uniform extrudate swell kinematics.

Velocity distribution at the slit die exit is illustrated in Figure 6, revealing that while a velocity core along the x-axis exhibits uniform magnitude (with isovels parallel to this axis) at central regions—consistent with infinitely broad slit assumptions for 1D flow—lateral decay initiates proximal to |x| = 6 mm, manifesting as curved isovels indicative of transitional 2D flow behaviour. Consequently, quantitative comparison between theoretical frameworks (relying on multi-centre velocity homogeneity) and numerical simulations is confined to the symmetry-preserving domain x ∈ [−6, 6] mm, where centreline velocity invariance remains experimentally validated.

A comparison of extrudate swell between numerical simulations and theoretical predictions is shown in Table 1. Quantitative agreement is achieved between numerical simulations and theoretical predictions at $[eqn]$ .

As shown in Figure 7, the x-direction swelling displacement profiles from both the theoretical model and numerical simulations exhibit a consistent progression. However, simulated values (see Table 1 for details) show a slightly steeper gradient than the theoretical curve. This discrepancy indicates amplified edge effects in the near-wall regions, likely due to localized viscous dissipation and velocity decay near the boundaries, phenomena not fully captured by the theoretical velocity core assumption.

Figure 8 reveals systemic alignment between theoretical and simulated y-direction swelling displacement profiles along the x-axis (Table 1), yet divergent behaviour emerges in boundary regions (|x|→6 mm). This divergence primarily originates from the infinite-slit assumption in theoretical derivation, wherein uniform velocity along the centreline ( $[eqn]$ ) contrasts with experimental velocity decay. Consequently, while x-directional displacement—dominated by far-field inertial effects—exhibits <5% deviation, y-directional displacement shows greater variance due to localized viscous dissipation in low-velocity boundary layers, consistent with non-Newtonian stress relaxation characteristics at n = 0.5.

The present theoretical model is grounded in the “infinitely wide slit” assumption, which postulates a uniform velocity profile in the width direction (i.e., ∂v/∂x ≡ 0). This simplification reduces the flow to a one-dimensional problem, enabling an analytical formulation consistent with classical die-swell theories that attribute swelling primarily to elastic recovery upon stress release [1,5]. However, as evidenced by numerical simulations (Figure 7), a discernible velocity decay emerges in the boundary regions (|x|→6 mm), signifying a transition from idealized one-dimensional flow to a more complex two-dimensional behaviour near the edges. This deviation directly underpins the observed discrepancies in the y-direction swelling displacement (Figure 8), where theoretical predictions diverge from simulated values at larger |x|. The velocity decay in these regions intensifies localized viscous dissipation and alters the stress relaxation dynamics—effects that are inherently omitted in the infinite-slit framework, particularly for shear-thinning power-law fluids (e.g., n = 0.5) [6]. As demonstrated, the high-accuracy agreement (errors < 5%) is firmly restricted to the core zone ( $[eqn]$ mm). Beyond this region, the infinite-slit assumption breaks down due to significant velocity decay near the side walls, causing the analytical predictions to deviate from CFD results with local edge errors exceeding 30%.

Despite this limitation, the model remains notably applicable and accurate for dies with high aspect ratios (W/2H > 10). In such configurations, the core flow domain (approximately |x| < 3 mm) retains a nearly uniform velocity, thereby upholding the key model assumption. Consequently, theoretical predictions align well with simulations in the central zone, with minimal deviation (e.g., <5% in the x-direction displacement). For comprehensive design applications, especially where edge precision is critical, future model enhancements could integrate edge-effect corrections derived from 3D flow analyses [7].

5. Scope and Limitations

The multi-velocity-centre hypothesis provides a mathematical mapping for anisotropic swell. However, it is crucial to acknowledge the fundamental limitations of the current theoretical framework and its numerical verification, which point to necessary directions for future research.

First, numerical verification is confined to an inelastic, isothermal power-law fluid (n = 0.5). Given that real polymer melts exhibit strong viscoelastic memory, this CFD validation primarily tests geometric and kinematic consistency of the proposed mapping, not the full physics of viscoelastic stress relaxation.

Second, the isothermal assumption neglects practical thermal effects. Temperature and viscosity gradients near the die walls can shift the velocity core and alter the anisotropic swell profile.

Third, the model’s infinite-slit assumption ( $[eqn]$ ) fails near the side walls, as evidenced by our results. Consequently, accurate predictions are confined to the core zone of high-aspect-ratio dies. For industrial dies with finite width, future work must incorporate edge-correction functions or hybrid CFD strategies.

Therefore, to transform this mathematical proof-of-concept into a predictive engineering tool, comprehensive physical extrusion experiments and non-isothermal viscoelastic simulations are essential next steps.

6. Conclusions

In this study, a theoretical mathematical framework is developed to geometrically describe the anisotropic extrudate swell profile of polymer melts in rectangular slit dies. Based on a detailed geometric and kinematic analysis, the following main conclusions are drawn:

• (1)Proof-of-Concept for Kinematic Mapping: Recognizing that actual polymer extrudate swell is fundamentally driven by viscoelastic recovery, this study explored a multi-velocity-centre hypothesis to approximate geometric and kinematic components of this phenomenon. This approach serves as a mathematical proof-of-concept, adapting the macroscopic diameter swell ratio to estimate the localized swelling profile.
• (2)Semi-Empirical Mathematical Formulation: By eliminating redundant geometric decompositions, a simplified analytical integration scheme was formulated. The model relies on the introduction of χ, which is strictly defined as a semi-empirical, lumped calibration parameter rather than a fundamental material property. This parameter serves as a proxy to capture the apparent local coupling between flow kinematics and macroscopic swelling.
• (3)Core-Zone Accuracy and Edge Limitations: The framework was subjected to a numerical test for kinematic consistency using isothermal, inelastic power-law fluid (n = 0.5) CFD simulations. The analytical predictions show good agreement (errors < 5%) strictly within the core zone of high-aspect-ratio dies. However, the infinite-slit assumption causes the model to break down in the edge regions, where 3D velocity decay leads to localized deviations exceeding 30%.
• (4)Future Prospects and Engineering Application: To transition this purely kinematic mapping into a generalized predictive model, future research must incorporate comprehensive physical extrusion experiments, non-isothermal viscoelastic constitutive modelling, and empirical edge-correction functions. Despite current limitations, this semi-empirical framework offers an initial geometric approach that may assist in the future development of inverse die-profile optimization and shape distortion compensation algorithms.