# Investigation on the Coalescence of Filaments in Fused Filament Fabrication Process for Amorphous Polyether Ketone Ketone for Varying Nozzle Diameters and Chamber Temperatures

Antoine Runacher, Thomas Joffre, Germain Fauny, Claudia Salvan, Nils Marchal, Francisco Chinesta

<figure><img src="image_1.png" /></figure>

## To cite this version:

Antoine Runacher, Thomas Joffre, Germain Fauny, Claudia Salvan, Nils Marchal, et al. Investigation on the Coalescence of Filaments in Fused Filament Fabrication Process for Amorphous Polyether Ketone Ketone for Varying Nozzle Diameters and Chamber Temperatures. Polymers, 2024, 16 (19), pp.2732. 10.3390/polym16192732. hal-05357092

HAL Id: hal-05357092

https://hal.science/hal-05357092v1

Submitted on 10 Nov 2025

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
Article

# Investigation on the Coalescence of Filaments in Fused Filament Fabrication Process for Amorphous Polyether Ketone Ketone for Varying Nozzle Diameters and Chamber Temperatures

Antoine Runacher <fnref n="1" />,<fnref n="2" />,* , Thomas Joffre <fnref n="1" />, , Germain Fauny <fnref n="1" />, Claudia Salvan <fnref n="1" />, Nils Marchal <fnref n="1" /> and Francisco Chinesta <fnref n="2" />

**Abstract:** The use of Fused Filament Fabrication (FFF) of high-performance polymers is becoming increasingly prevalent, leading to the exploration of new applications. The use of such materials in critical cases for aerospace applications necessitates the verification of industry standards, particularly with regard to the requirements for part porosity. The authors investigate the effect of nozzle diameter and cooling temperature printing parameters on the porosity of the part by using existing modelling methods based on the sintering of cylinders and spheres and comparing the results to microscope snapshots of sections of parts. The models are able to be used as limits for predicting the longitudinal neck growth of the part. The authors show through experiments that the value of the cooling temperature of the deposited filament has a minimal effect on the outcome, while nozzle diameter has a strong impact on the resulting porosity. The modelling results show that there is a significant impact of both the nozzle diameter and cooling temperature on the porosity of the part. This implies that further refinement of the models is needed for the resulting parts to be applied in critical structures.

**Keywords:** fused filament fabrication; high-performance polymer; filament modelling; sintering; coalescence; microscope imaging; nozzle diameter; cooling temperature; mesostructure; porosity

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Runacher, A.; Joffre, T.; Fauny, G.; Salvan, C.; Marchal, N.; Chinesta, F. Investigation on the Coalescence of Filaments in Fused Filament Fabrication Process for Amorphous Polyether Ketone Ketone for Varying Nozzle Diameters and Chamber Temperatures. *Polymers* **2024**, *16*, 2732. https://doi.org/10.3390/polym16192732

Academic Editor: Yang Li

Received: 29 August 2024

Revised: 16 September 2024

Accepted: 19 September 2024

Published: 27 September 2024

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2024 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

## 1. Introduction

Since its birth in the 1960s, Additive Manufacturing (AM) methods have been rapidly developing. While not as efficient as traditional subtractive and formative methods, the quick design-to-part time and the customisation possibilities of AM have contributed to its development, first for prototyping endeavours and now even for structural uses of parts [1].

A large number of such methods were developed over the years for metals and polymer materials, each having its own drawbacks and perks. This multitude of methods has resulted in a rapid increase in scientific publications regarding this set of methods and also contributed to the expansion of the field [2].

One of these fabrication processes is Fused Filament Fabrication (FFF) [3]. By feeding a polymer through a heating nozzle, a melted filament of material is extruded following a specific path [4]. Adjacent filaments will fuse to form a solid ensemble after the cooling of the material [5]. Successive layers of material deposited on top of each other will then form a finished part. This process can be applied to a wide array of materials [6–8] to produce parts of complex geometry, without moulds or welds. It is, however, slower than traditional manufacturing, the manufactured parts are limited in size, and their mechanical properties are often worse than parts fabricated with the injection process [9].

<footnote n="1">¹ Centre Technique Innovation Plasturgie et Composites, 2 Rue Pierre et Marie Curie, 01100 Bellignat, France; nils.marchal@ct-ipc.com (N.M.); germain.fauny@ct-ipc.com (G.F.); claudia.salvan@ct-ipc.com (C.S.); thomas.joffre@ct-ipc.com (T.J.); francisco.chinesta@ensam.eu (F.C.)</footnote>

<footnote n="2">² Laboratoire PIMM, Ecole Nationale des Arts et Métiers, 151 Boulevard de l'Hôpital, 75013 Paris, France</footnote>

<footnote>* Correspondence: antoine.runacher@ensam.eu</footnote>
There exists a market for FFF-manufactured parts: while the manufacturing times are high, there are no added costs of moulds. For a low volume of parts, these characteristics can be more advantageous than injection or welding. This is particularly the case in the aviation industry, as the parts made using FFF are lighter than their aluminium counterparts used in planes today and boast very low excess material consumption. Still, the high safety requirements of critical-use parts need to be addressed [10].

To satisfy this drive for FFF-processed parts in primary aerospace structures, stringent safety requirements must be implemented. The first is the need for fire-resistant materials that do not produce smoke when burned, produce fumes that are safe for users, and have a high melting temperature. The second requirement is to keep part porosity below 2% [11].

The PAEK family of thermoplastic materials (PAEK, PEKK, PEEK, ...) complies with the specified safety regulations but is notoriously hard to use in Fused Filament Fabrication [12], requiring a heated chamber for the coalescence between filaments to occur. By increasing this temperature, the capillary time during which the coalescence takes place is longer, resulting in higher part strength [13].

During the deposition process, two adjacent beads of material are placed in intimate contact with each other. The high temperature of printing starts the bond formation process between filaments, which is the fusing of polymer chains at a molecular scale. When the cooling step of the process is finished, the healing of the part is characterized by the interdiffusion of polymer chains, from one bead to the next [14].

This adhesion mechanism is shown in Figure 1.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Bond formation process between two filaments: (1) surface contact; (2) neck growth; (3) molecular diffusion at interface and randomisation—inspired from [15].</figcaption></figure>

The finished parts have varying levels of coalescence, causing interrogations regarding the structural soundness of the part and hence safety requirements regarding part porosity.

This phenomenon was first modelled by Frenkel and Eshelby in 1945 and 1949 [16]. They described the coalescence of two spherical beads of material modelled as an incompressible and viscous Newtonian fluid in isothermal conditions with a uniform stress tensor. This problem was solved with a small angle approximation.

Pokluda et al., in 1997 [17], expended the Frenkel–Eshelby model for all angles with a Runge–Kuta numerical method to solve the differential heat equation. This model provides a dimensionless solution to the problem of neck growth of PLA beads. From this, the expected porosity should not vary with a change in the radius of the considered beads.

As shown by Bakrari Balani [18], the surface tension of the material changes during the cooling process and so does the viscosity. An accurate characterisation of both is needed to model the behaviour of the porosity of a printed part.

Lepoivre et al., in 2021 [19], characterized the evolution of the surface tension and viscosity of PEKK with regard to the temperature and then used a numerical model to describe the evolution of the degree of coalescence with varying surface tension and viscosity. Their findings are used to accurately calculate the surface tension and viscosity of the PEKK for the presented model.
In 2023, Jiang et al. [20] investigated the coalescence of PLA in FFF by modelling the behaviour of the coalescence and testing the results on printed filaments of varying sizes, but the long diameter of the filaments studied was limited to a nozzle diameter of 0.4 mm.

There is a need to determine the behaviour of larger filaments, as the increase in cooling time of the filaments increases with each increment in diameter. The effects of the chamber temperature, a parameter with a sure impact on a part's Young modulus, are not often seen in the literature [21].

In this publication, a framework based on Jiang et al.'s investigations, with Lepoivre et al.'s data on the viscosity and surface tension of PEKK, explores the neck growth of PEKK filaments. This study will show the impact of nozzle diameter and cooling temperature on the longitudinal neck growth between filaments. The authors have not found other examples in the literature of similar studies for the considered nozzle diameter.

These investigations on the neck growth will also be conducted in different chamber temperatures to quantify the effect of the parameter on the porosity of the part.

## 2. Materials and Methods

This section pertains to the methodology used to analyse section cuts of printed PEKK test pieces.

### 2.1. Modelling of the Temperature of the Bead

The cooling of a single filament is simplified to a one-dimensional transient heat transfer model. This method works under the assumption of a uniform temperature distribution throughout the cross section of the deposited filament [16]. The nozzle moves at a constant speed $v$. The temperature of the substrate is considered constant throughout the deposition and cooling process, from which the conduction of heat between layers is considered convection. Figure 2 shows the energy interaction on a finite element $dx$ of a heated filament.

<figure><img src="image_6.png" /><figcaption>Figure 2. Energy interaction on finite element $dx$, inspired by [22].</figcaption></figure>

From this, the heat differential equation of the system can be written.

$$ \rho C A \frac{\partial T}{\partial t} = A \frac{\partial(k \frac{\partial T}{\partial x})}{\partial x} - h(T - T_{\infty}) \quad (1) $$

Where $T = T_0$ at $x = 0$, $t \ge 0$ and $T = T_\infty$ at $x = \infty$ and $t \ge 0$, and $C$ is the specific heat coefficient, $\rho$ is the density, $A$ is the cross section area of the filament, $k$ is the heat conductivity coefficient and $h$ is the convection coefficient of the system.

The coefficient $h$ governs the effects of both heat convection with air and conduction with the substrate.
The nozzle is moving at a velocity $v$ and $x = vt$. The time-dependent term $\frac{\partial T}{\partial t}$ can be rewritten as

$$ \frac{\partial T}{\partial t} = \frac{\partial T}{\partial x} \frac{\partial x}{\partial t} = \frac{\partial T}{\partial x} v \qquad (2) $$

The previous equation then becomes:

$$ \rho C A v \frac{\partial T}{\partial x} = A \frac{\partial(k \frac{\partial T}{\partial x})}{\partial x} - h(T - T_{\infty}) \qquad (3) $$

And has the solution:

$$ T(x) = T_{\infty} + (T_0 - T_{\infty})e^{-mx} \qquad (4) $$

Where $m = \frac{\sqrt{1+4\alpha\beta-1}}{2\alpha}$, $x = vt$, $\alpha = \frac{k}{\rho C v}$, and $\beta = \frac{hP}{\rho C v A}$. $P = \pi(a+b)(\frac{64-3\lambda^4}{64-16\lambda^2})$ and $A = \pi ab$, with a as the major axis and b the minor axis of the cross section of the deposited filament, considered elliptic, and where $\lambda = \frac{a-b}{a+b}$.

## 2.2. Modelling of Coalescence between Two Beads of Material

The porosity of the part is related to the sintering of two adjacent deposited filaments. The polymer filaments are printed at a temperature higher than the glass transition temperature. At those temperatures, the polymer is considered a Newtonian fluid.

Polychronopoulos and Vlachopoulos [23], in 2020, studied the role of heating and cooling in the sintering of both spheres and cylinders for additive manufacturing, detailing the steps to solve the coalescence problem of Fused Filament Fabrication.

The sintering process or coalescence is illustrated in Figure 3.

The sintering is the behaviour that two adjacent liquid beads follow to form one unique bead of a longer radius. It is a phenomenon that is dependent on the surface tension of the fluid, $\Gamma$, and the viscosity of the fluid $\eta$. The sintering process occurs until the temperature of the material crosses below the glass transition temperature as highlighted in the previous section. The material is then considered solid and will not be subject to further deformations related to the sintering of the beads. The sintering process is modelled from the dynamic equilibrium between the work of the surface tension forces and the work of the viscosity forces in the given fluid.

In the following sections, the equations leading to the expression of the evolution of the half angle of coalescence $\theta$ are reported. It is obtained in the form of a differential equation.

From $\theta$, geometry relations between quantities give the theoretical longitudinal neck growth—the quantity of interest.

<figure><img src="image_7.png" /><figcaption>Figure 3. Sintering process of two cylinders or two spheres, inspired by [23].</figcaption></figure>

### 2.2.1. Sintering of Spheres

If two spheres of radius $a_0$ intersect at one point of contact, they will gradually coalesce into one sphere of radius $a_f$ with $a_f = a_0\sqrt[3]{2}$, with that point of contact as the centre of the newly formed sphere.
As shown in Figure 3, at a given time t, the values of the radius of the bead a(t), the half angle of coalescence θ(t), and the longitudinal neck growth x(t) will be referred to as a, θ and x, respectively, in the following.

Pokluda et al. [17] established a framework to model the behaviour of the sintering of beads. The radius a of the beads is given by the conservation of the mass of the polymer, with the density of a bead constant with the temperature.

$$a = a_0 \left( \frac{4}{[1 + \cos \theta]^2 [2 - \cos \theta]} \right)^{1/3} \qquad (5)$$

The area of a section S of two beads at a given time is given by Equation (6).

$$S = 4\pi a^2 [1 + \cos \theta] \qquad (6)$$

The flow of material for spherical bead sintering is given by a biaxial extensional flow field. This leads to solving the dynamic equilibrium between the work of viscous forces and surface tension, the expression of a differential equation θ, and the half angle of the coalescence of the beads.

Polychronopoulos et al. [23] show the steps to the solution, and the expressions of the quantities of interest θ, the half angle of coalescence, and y, the longitudinal neck growth between filaments, as follows. The differential equation in θ is solved using a Runge-Kuta algorithm with an initial value of θ, θ₀ = 0.01 to avoid numerical instabilities.

$$\dot{\theta} = \frac{\Gamma 2^{-5/3} \cos \theta \sin \theta (2 - \cos \theta)^{1/3}}{a_0 \eta (1 - \cos \theta) (1 + \cos \theta)^{1/3}} \qquad (7)$$

Finally, for spherical sintering, the longitudinal neck growth y is given by:

$$y = a_0 \sin \theta \left( \frac{4}{(1 + \cos \theta)^2 (2 - \cos \theta)} \right)^{1/3} \qquad (8)$$

In the following sections, the half angle of coalescence and the longitudinal neck growth for spherical beads of material will, respectively, be noted θsph and ysph

### 2.2.2. Sintering of Cylinders

In the case of fused filament fabrication, the sintering process may also be modelled as two elliptic cylinders fusing, which can be characterized by two circular cylinders of a unique radius $r_0 = b^2/a$ [20], with a the semi-major axis and b the semi-minor axis of the ellipse, as described in Figure 4. The ultimate result of such coalescence is a unique cylinder in which the temperature of the polymer is maintained high enough for the sintering to continue. The variation in the length of the cylinder will be ignored as it is negligible before the variation in the radius of the cylinder. For ease of comprehension and for the rest of the section, the functions r(t) (the radius of the two cylinders), θ(t) (the half angle of coalescence, as shown in Figure 4), and y(t) (the longitudinal neck growth of the process) will be referred to as r, θ, and y, respectively.

<figure><img src="image_8.png" /><figcaption>Figure 4. Diagram of bonding, inspired by [20].</figcaption></figure>
The area of a circular segment of a disk formed by the section of two cylinders is as follows [24] and is shown in Figure 5, where R is the radius, Θ is the central angle in radians equal to 2θ, c is the cord length, s is the arc length, h is the sagitta, d is the apothem, and Aseg is the area of the segment (orange area in Figure 5):

<figure><img src="image_9.png" /><figcaption>Figure 5. Area of a segment of a circle.</figcaption></figure>

$$A_{seg} = \frac{R^2}{2}(\Theta - \sin\Theta) = R^2(\theta - \cos\theta \sin\theta) \quad (9)$$

From which the area of the section of two cylinders during the sintering process A_c shown in orange in Figure 4 is given by:

$$A_c = 2(A - A_{seg}) = 2\pi r^2 - 2r^2(\theta - \cos\theta \sin\theta) = 2r^2(\pi - \theta + \cos\theta \sin\theta) \quad (10)$$

From the conservation of mass between the steps of the sintering, the radius r of the circular cylinder is given [23]:

$$r = \frac{\sqrt{\pi r_0}}{\sqrt{\pi - \theta + \cos \theta \sin \theta}} \quad (11)$$

Where $r_0 = \frac{b^2}{a}$ and θ is the half angle of the coalescence, as shown in Figure 4.

From the given expressions of nominal radius r and area A_c and the dynamic equilibrium between the work of viscous forces and the surface tension, a differential equation with the half angle of coalescence θ is given.

The flow of material induced by the viscous forces and the surface tension of the fluid is a planar extensional flow when the sintering process is approaching the sintering of cylinders, meaning that the work of viscous forces is expressed as a function of the strain rate of the system. By definition, the strain rate and the velocity of the deformation are linked, and so the dynamic equilibrium between viscous forces and surface tension that exists in the system when the gravity is neglected results in a differential expression of θ, i.e., the half angle of coalescence.

Polycronopoulos et al. [23] describe the steps to the solution of this problem, and the quantities of interest θ and y are given by Equations (12) and (13). The differential equation in θ is solved using a Runge-Kuta algorithm with an initial value of θ, where θ₀ = 0.01 to avoid numerical instabilities.

$$\dot{\theta} = \frac{\Gamma}{2\sqrt{\pi\eta r_0}} \frac{[(\pi - \theta)(\cos\theta) + \sin\theta][\pi - \theta + \cos\theta \sin\theta]^{\frac{1}{2}}}{(\pi - \theta)^2 \sin\theta \tan\theta} \quad (12)$$

$$y = r_0 \frac{\sin\theta}{\sqrt{\pi - \theta + \cos\theta \sin\theta}} \quad (13)$$

In the following sections, θ and y for cylinder sintering will be referred to as θcyl and ycyl.
## 2.3. Material

The material is the 60–40 amorphous ThermaX Polyether Ketone Ketone polymer (PEKK-A) formulated by 3DXTech in Grand Rapids, MI, USA and made using Arkema's Kepstan [60/40 copolymer] fabricated in Mobile, AL, USA. It is amorphous and contains only PEKK. Following the manufacturer processing parameters, the printed temperature is 325–360 °C, and the chamber temperature is 70 to 150 °C. It has a glass transition temperature of 162 °C and a melting temperature of 335 °C.

## 2.4. Specimen Printing

The printed test specimens are created according to the 1BA Standard EN ISO 527-2 [25] dimensions. The .STL file of the corresponding geometry was sliced into a readable G-code for a funmat 610 HT printer fabricated by INTAMSYS in Shanghai, China.

The G-code path consists of superposed deposited filaments, as shown in Figure 6.

<figure><img src="image_10.png" /><figcaption>**Figure 6.** G-code path.</figcaption></figure>

The values of the varying parameters of the factorial Design Of Experiment (DOE) used are given in Table 1. The DOE has three levels and two factors, where cooling temperature $T_{inf}$ has the values 110, 130 and 150 °C, and the nozzle diameter $D$ has the values 0.4, 0.6 and 0.8 mm. Printing failed for cooling temperatures below 110 °C. For each specified nozzle diameter $D$ (mm) and chamber temperature $T_{inf}$, a set of test pieces will be produced, resulting in nine parameter couples and samples. The layer thickness $H$ of the part is given by $H = D/2$, as the optimum for print quality in the field is given for this value of $H$ [26].

<table><caption>**Table 1.** Test piece DOE.</caption><thead><tr><th>Test Piece</th><th>X<sup>0.4</sup><sub>110</sub></th><th>X<sup>0.4</sup><sub>130</sub></th><th>X<sup>0.4</sup><sub>150</sub></th><th>X<sup>0.6</sup><sub>110</sub></th><th>X<sup>0.6</sup><sub>130</sub></th><th>X<sup>0.6</sup><sub>150</sub></th><th>X<sup>0.8</sup><sub>110</sub></th><th>X<sup>0.8</sup><sub>130</sub></th><th>X<sup>0.8</sup><sub>150</sub></th></tr></thead><tbody><tr><td>T<sub>inf</sub> (°C)</td><td>110</td><td>130</td><td>150</td><td>110</td><td>130</td><td>150</td><td>110</td><td>130</td><td>150</td></tr><tr><td>D (mm)</td><td>0.4</td><td>0.4</td><td>0.4</td><td>0.6</td><td>0.6</td><td>0.6</td><td>0.8</td><td>0.8</td><td>0.8</td></tr><tr><td>H (mm)</td><td>0.2</td><td>0.2</td><td>0.2</td><td>0.3</td><td>0.3</td><td>0.3</td><td>0.4</td><td>0.4</td><td>0.4</td></tr></tbody></table>

Other process parameters, such as printing temperature, also have an impact on the porosity of the part and are fixed to remain constant in the DOE. Test pieces were printed for each varying parameter couple. For $D$, the given nozzle diameter, and $T_{inf}$, the given chamber temperature of a test piece, the denomination $X_{T_{inf}}^D$ will be given to the test piece.

## 2.5. Specimen Preparation

The test specimens were then cut in the YZ axis using an Isomet 1000 precision cutter into 6 observable surfaces.

The cut pieces of the testing specimen were cured inside a transparent KM-U resin from Presi, made in Eybens, France.

The resin cylinders were then polished using a MECAPOL P200 polishing machine made by Presi in Eybens. The abrasive disc P180 was used to denude the surfaces of interest; then, the abrasive discs P360, P840, P1000, P2400 and P4000 were used successively for 5 min each, with water as a suspension. Subsequently, a TOP polishing disk was applied for 5 min with a 6 µm Reflex LDM suspension before a RAM polishing disk with a 3 µm Reflex LDM suspension.
The polishing machine head has a speed of 600 laps per minute with all disks and applies a pressure of 1.1 daN onto the polishing head. The resin cylinders were maintained as immobile during the process.

## 2.6. Mesostructure Observation

The samples were observed using a VHX-7100 microscope, made by Keyence in Osaka, Japan, where two captures of each surface of interest were taken. The lighting of the scene was made by coaxial lights and then optimized thanks to the microscope software. The images were then binarised using image processing tools.

There are 6 such mesotructures available for each pair of chamber temperature and nozzle diameter.

The mesostructures were then cut by layers, where the longitudinal length of the voids was measured and subtracted to the height of the layer, as shown in Figure 7, to obtain the longitudinal neck growth of this particular element in the layer. This operation is repeated for all voids in the layer.

In Figure 7, the longitudinal neck growth is obtained by subtracting the longitudinal void length (represented in green) from the layer height (represented in orange).

<figure><img src="image_11.png" /><figcaption>**Figure 7.** Measurement of longitudinal neck growth of one coalescence occurrence in a designated layer. The colour lengths are the layer height (orange) and the longitudinal void length (green).</figcaption></figure>

The results of all the measurements are then plotted in the next section.

# 3. Results

The results of the various investigations are shown and are divided into three parts. The first section presents the results of the theoretical model used to predict the length of the neck between filaments, the second displays the experimental data measured using the microscope, and the third compares the results of the two previous sections. All the results are then discussed in the following section.

## 3.1. *Theoretical Longitudinal Neck Length*

Using the introduced parameters, the theoretical thermal temperature of the filament is given in Figure 8:
<figure><img src="image_12.png" /><figcaption>**Figure 8.** Theoretical temperature of the deposited filament.</figcaption></figure>

From this, the coalescence time $t_c$ given to the polymer to form bonds in each case is obtained. It is shown as the time at which the temperature of the filament crosses the glass transition temperature in Figure 8. This coalescence time is then used as a limit in the time integration to obtain the theoretical half angle $\theta$ by solving the differential Equation (7) for spherical sintering and (12) cylinder sintering. By definition, an amorphous polymer under its glass transition temperature is considered solid, but the conformal changes in segments are considered infinitely slow at the glass transition temperature. By offering the possibility to deform until the glass transition temperature, the assumption is that the deformation modelled may be superior to the experimental deformation. In this study, the viscosity of the polymer changes with the temperature of the polymer. Below 200 °C, the model shows a negligible increase in longitudinal neck growth, which is in agreement with the expected behaviour of the polymer, with the majority of the deformation taking place in the first instants after the deposition.

From Equations (8) and (13), the theoretical neck growth is then obtained. All the theoretical results have been compiled in Table 2.

<table><caption>**Table 2.** Theoretical results of capillary time $t_c$, theoretical half angle of coalescence $\theta_{cyl}$ and $\theta_{sph}$, theoretical longitudinal neck growth $y_{cyl}$ and $y_{sph}$ for both spherical and cylinder sintering.</caption><thead><tr><th>Test Piece</th><th>X<sup>0.4</sup><sub>110</sub></th><th>X<sup>0.4</sup><sub>130</sub></th><th>X<sup>0.4</sup><sub>150</sub></th><th>X<sup>0.6</sup><sub>110</sub></th><th>X<sup>0.6</sup><sub>130</sub></th><th>X<sup>0.6</sup><sub>150</sub></th><th>X<sup>0.8</sup><sub>110</sub></th><th>X<sup>0.8</sup><sub>130</sub></th><th>X<sup>0.8</sup><sub>150</sub></th></tr></thead><tbody><tr><td>T<sub>inf</sub> (°C)</td><td>110</td><td>130</td><td>150</td><td>110</td><td>130</td><td>150</td><td>110</td><td>130</td><td>150</td></tr><tr><td>D (mm)</td><td>0.4</td><td>0.4</td><td>0.4</td><td>0.6</td><td>0.6</td><td>0.6</td><td>0.8</td><td>0.8</td><td>0.8</td></tr><tr><td>H (mm)</td><td>0.2</td><td>0.2</td><td>0.2</td><td>0.3</td><td>0.3</td><td>0.3</td><td>0.4</td><td>0.4</td><td>0.4</td></tr><tr><td>t<sub>c</sub> (s)</td><td>5.82</td><td>7.22</td><td>10.43</td><td>8.73</td><td>10.83</td><td>15.64</td><td>11.53</td><td>14.44</td><td>20.85</td></tr><tr><td>θ<sub>cyl</sub> (rad)</td><td>0.26</td><td>0.3</td><td>0.35</td><td>0.35</td><td>0.39</td><td>0.43</td><td>0.41</td><td>0.45</td><td>0.49</td></tr><tr><td>θ<sub>sph</sub> (rad)</td><td>0.18</td><td>0.22</td><td>0.27</td><td>0.28</td><td>0.33</td><td>0.38</td><td>0.36</td><td>0.4</td><td>0.45</td></tr><tr><td>y<sub>cyl</sub> (mm)</td><td>0.13</td><td>0.15</td><td>0.17</td><td>0.26</td><td>0.29</td><td>0.32</td><td>0.41</td><td>0.44</td><td>0.48</td></tr><tr><td>y<sub>sph</sub> (mm)</td><td>0.09</td><td>0.11</td><td>0.14</td><td>0.21</td><td>0.25</td><td>0.29</td><td>0.36</td><td>0.41</td><td>0.47</td></tr></tbody></table>

## 3.2. *Experimental Results*

The pictures of the cross section of test pieces for different nozzle diameters were taken via microscope. An example of such mesostructure is shown in Figure 9.
<figure><img src="image_13.png" /><figcaption>**Figure 9.** Microscope image for 0.4 mm nozzle, 110 °C chamber temperature: $X_{110}^{0.4}$.</figcaption></figure>

After image binarisation and cropping, Figure 10 shows the mesostructure of parts made with nozzle diameters of 0.4 mm, 0.6 mm and 0.8 mm at different chamber temperatures.

<figure><img src="image_14.png" /><figcaption>**Figure 10.** Examples of mesostructures of printed test pieces at different chamber temperatures $T_{\text{inf}}$ and nozzle diameters $D$ noted $X_{T_{\text{inf}}}^{D}$.</figcaption></figure>

This process resulted in more than 60 measures of longitudinal neck growth per DOE set.
### 3.3. *Experimental and Theoretical Comparison*

From the mesostructure, the longitudinal neck growth values are extracted, as described previously, and compiled by layer and chamber temperature in Figure 11 for nozzle diameter *D* = 0.4 mm and Figure 12 for nozzle diameter *D* = 0.6 mm.

<figure><img src="image_15.png" /><figcaption>**Figure 11.** Neck growth for nozzle diameter *D* = 0.4 mm, for all layers and cooling temperatures (*T*<sub>inf</sub> = 110 °C in red, *T*<sub>inf</sub> = 130 °C in blue, *T*<sub>inf</sub> = 150 °C in green), for experimental data (boxplot), and models (spherical model as dotted line and cylinder model as line).</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>**Figure 12.** Neck growth for nozzle diameter *D* = 0.6 mm, for all layers and cooling temperatures (*T*<sub>inf</sub> = 110 °C in red, *T*<sub>inf</sub> = 130 °C in blue, *T*<sub>inf</sub> = 150 °C in green), for experimental data (boxplot), and models (spherical model as dotted line and cylinder model as line).</figcaption></figure>

To help visualize the data, it was plotted in the form of a box plot, where the limits of the box represent the first and third quartiles of the data, and the whiskers represent the limits of the data. Outliers were plotted separately beyond the whiskers of the plot. The mean value of the data for each layer and DOE set was plotted as the black line in the box. The standard deviation of the data *s* can be approximated through the following equation:

$$s = \frac{3}{4}(Q_3 - Q_1), \text{ with } Q_1 \text{ the first quartile and } Q_3 \text{ the third quartile.}$$

From the shape of the porosity shown in the experimental results displayed for nozzle diameter *D* = 0.8 mm, the sintering of the part is not the principal phenomenon behind the porosity of the part, so the comparison with the theoretical model is not conclusive.

In the abscissa of the graphs, *L* is the abbreviation of layer. It is followed by the number of the considered layer, starting with 1 for the bottom layer in contact with the raster. *T* is the abbreviation of temperature and is followed by the value of the chamber temperature corresponding to that particular print.
In the following section, the term “theoretical longitudinal neck growth” will refer
to both the cylinder- and spherical-modelled neck growths, representing an interval of
possible neck growth measurements. The spherical calculation serves as the lower limit of
the interval, and the cylinder calculation serves as the upper limit.

As shown in the previous section, the longitudinal neck growth is not the driving
force behind the porosity observed for *D* = 0.8 mm, and the results of the experimental
measurements are not the longitudinal neck growth of the filaments. The image analysis of
the printed parts with a nozzle diameter of *D* = 0.8 mm is not relevant to the comparison
of the model with experimental results.

For D = 0.4 mm, for the first few layers, the coalescence measured is less than the theoretical value: this is due to two phenomena. Firstly, the absence of contact between filaments is observed, as shown in Figure 10c, leading to measured longitudinal neck growth equal to zero, unlike the predictions of the model. This behaviour is observed for elements printed with all chamber temperatures, but mostly for T<sub>inf</sub> = 110 °C and T<sub>inf</sub> = 150 °C. Secondly, the values of the non-zero elements is also lower than expected, which is due to the proximity to the heating bed. The assumption that the conduction of heat through the material does not change its temperature is false: the substrate has a similar mass to the deposited filament, and hence, changes in the temperature in the substrate are induced by the presence of the filament. Therefore, the cooling of the material is faster than expected, leading to a lower experimental longitudinal neck growth measured than the theoretical value.

For the middle to top layers, the theoretical and experimental estimations of the
neck growth results show an overlap as the majority of experimental measurements are
contained between the two boundaries defined by the models, but a large number of
experimental measurements show a quasi-total neck growth. This is due to the change
in geometry observed in the mesostructure: the void shape going from diamonds to
triangles. In the literature, occurrences of this are linked with over-extrusion, as shown in
Ghorbani et al. [27].

For *D* = 0.6 mm, there is contact between filaments in the first layers for all temper-
atures. As observed for *D* = 0.4 mm, the neck growth for the first deposited layers is
reportedly smaller than that of the middle or top of the printed test piece. The observed
voids are almost all triangle-shaped. A complete absence of voids is also observed in
several areas of the section. Large defects are also starting to appear in the top layer, leading
to poor surface quality of the part for chamber temperature *T*<sub>inf</sub> = 150 °C. The filaments
deposited in the last layers have a rounder appearance than those deposited in the first
layer, while the extrusion speed of the material has not changed.

4. Discussion

4.1. *Physics Analysis*

This section pertains to the analysis of the physics governing the results previously
displayed. Two phenomena are discussed: the first is the filament shape, and the second is
the coalescence between filaments.

4.1.1. Filament Shape

The model considers a perfect ellipse or sphere on which the equations of the system are solved. The experiments show that the shape of the filament is not so simple to describe. First, the nozzle is a main contributor to the shape of the deposited filament. Second, the asymmetry of the deposited filament is due to small variations in the angle of the nozzle during the printing process.

* For a nozzle diameter *D* = 0.4 mm, and while the material is printed, the nozzle constrains the expansion of the material on the Z-axis by its presence. The excess material is pushed to the sides of the nozzle. As observed in Figure 10, the centre of mass of the filament is higher than the middle of the filament on the Z-axis, and the
skinny part of the filament is at the bottom. This particular deformation is a result of
the path of the nozzle, restraining the degrees of freedom of the deposited material.

* For nozzle diameter *D* = 0.6 mm, the shape of the filament shows a centre of mass lower than the middle of the filament on the Z-axis, and the skinny part is at the top. As stated in the previous section, the aspect of the filament is rounder for filaments at the top of the print for cooling temperature *T*<sub>inf</sub> = 150 °C.

* For nozzle diameter *D* = 0.8 mm, the filaments in the first layer show little to no voids, and the flow of material induced by the high printing temperature was enough to change the behaviour of the deposited filament. Higher up in the test piece, the aspect ratio of the observable filaments is near 1, meaning that the nozzle does not constrain the expansion of material on the Z-axis: the substrate of the considered layer is lower than expected. A direct consequence is that the geometry of the part shows defects, and the mesostructure is no longer easily observable.

* Small variations in the angle of the nozzle during the deposition process result in a favoured direction of flow of material. The resulting filaments show asymmetry, with the centre of mass of the filament skewing to either right or left. This small change in shape may result in changes in the mesostructure of the complete test piece. This phenomenon is observable for all nozzle diameters and cooling temperatures.

4.1.2. Coalescence between Filaments

The coalescence between filaments varies a lot with the process parameters and, during
the process, with the number of layers deposited. Here, the different types of coalescence
observed are shown: absence of coalescence, coalescence resulting in diamond-shaped voids,
coalescence resulting in triangular-shaped voids, and coalescence resulting in no voids.

* For nozzle diameter *D* = 0.4 mm and cooling temperature *T* = 110 °C or *T* = 150 °C, the first layers of deposited filaments display instances of the absence of coalescence. The asymmetry observed in the previous section, resulting in an absence of contact between filaments, stalls the coalescence between neighbouring filaments.

* For nozzle diameter *D* = 0.4 mm, the voids can be in the shape of diamonds, formed by the coalescence of two elliptic filaments from the previous layer, and two from the current layer, meaning that the neck growth is as modelled in the theoretical part. It is, nevertheless, not the only type of void shape observed.

* For nozzle diameters *D* = 0.4 mm and *D* = 0.6 mm, the voids are often in the shape of triangles, meaning that the bottom left with the top right filament or top left with bottom right filament have started the sintering process, leading to two voids of smaller size and of triangular shape, or just one triangular-shaped void, as one of the two disappears during the sintering process. This is due, in part, to the asymmetry of the deposited filaments but also to the flow of material observed for *D* = 0.6 mm filaments. This is seen in Figure 10d–f as the filament bottom part is sintering with its neighbours, but the top part is not, creating the aforementioned voids, when voids are created.

* For nozzle diameter *D* = 0.8 mm, the analysis of the coalescence is not relevant, as the main contributor of voids is not the coalescence but the printing defects observed that resulted from the flow of material.

4.2. *Theoretical and Experimental Results Discussion*

The theoretical models describe two flow states between spherical, corresponding
to a local sintering of the filament with its neighbour, and cylinder, corresponding to the
global sintering of the whole filament length with its neighbour. The size of the zone heated
by the deposition of a filament in its neighbour is the parameter that governs the flow
of material from one filament to the next. If the heated zone is small before the length
of the filament, then the flow state will approach a spherical sintering (the contact with
the neighbouring filament changes the kinematics of the problem, and the flow state is
not strictly spherical). If the heated zone is comparable to the length of the filament, then
the flow state will correspond to a cylinder-shaped sintering. The reality is supposedly a combination of the two.

The given models may produce limits for the porosity inside a printed part, as the results for *D* = 0.4 mm and *D* = 0.6 mm show a significant overlap between the theoretical and experimental results, with both the spherical and cylinder sintering predictions acting as limits to the experimental results. There are cases where the prediction is not as precise:

* In the first layers of the print, this behaviour is changed as the conduction of the heat between the filament and the substrate is faster than expected, leading to smaller longitudinal neck growths measured.

* For *D* = 0.6 mm, the Biot number of the filament *Bi* = *h* *V* / *k* *A*, is very similar to 0.1. The approximation giving the uniformity of the temperature in the section of the considered filament may not be verified and the actual temperature of the point of contact between filaments may be lower than modelled, leading to the overestimation of the longitudinal neck growth by the considered models.

The models provided are also constant throughout the printing process and do not take into account the layer number of the print, while the experimental results show a clear influence of the layer number and especially the first two deposited layers. There are also deposited filaments that have a very different aspect ratio than other filaments in the part. They are found in the higher layers of the print.

The theoretical longitudinal neck growth may be equal or greater than the layer height *H* = 0.3 mm. While the experimental neck growth is bounded by the process parameters, as the method used to measure the experimental neck growth requires a void length to be subtracted from the layer height, it is not the case for the theoretical longitudinal neck growth. The cases in which the coalescence is such that there is no void between filaments were simply not measured. The rest are represented in Figure 12.

### Impact of the Cooling Temperature on Porosity

While the models show a clear dependence on the cooling temperature of the process, the experimental results show similar levels of mean longitudinal neck growth for all cooling temperatures. The impact of the cooling temperature on the part is mainly seen in the variations in neck growth measured between layers.

The porosity of the printed parts shows little to no variation with regards to the cooling temperature, but the models do. While it is often true that the lower the cooling temperature, the smaller the neck growth, the variation in neck growth is small. Still, the spherical model is better at approximating the experimental data than the cylinder model, especially for *T*<sub>inf</sub> = 150 °C, where both models predict higher values than the experimental data.

### 4.3. *Hypothesis for the Observed Differences*

In this section, the differences in results are explored. First, the differences between experimental results for different parameters are hypothesised. Secondly, the differences between the theoretical and experimental results are explained.

#### 4.3.1. Differences between Experimental Results for the Given Parameters

The differences described in the previous sections can be explained thanks to a few hypotheses. The differences mentioned can be grouped by the phenomena responsible for their existence. The first group of differences is as follows:

* The changes in aspect ratio observed for *D* = 0.6 mm and *D* = 0.8 mm but not for *D* = 0.4 mm. The aspect ratio for the last layer for *D* = 0.6 mm and the last few layers for *D* = 0.8 mm is much closer to 1 than expected.

* When the observed aspect ratio of the filament, given by the ratio of *H*/*D*, is equal to 0.5, the position of the centre of mass of the filament changes with the nozzle diameter: the higher the nozzle diameter, the lower the centre of mass.
An explanation for these few differences might be the effect of gravity on the filament and the flow of material, causing the filament to “sag”. The changes in the position of the centre of mass of the filament show that a force is applied to the filament and acts on the direction of the flow of material. This force is negligible for *D* = 0.4 mm but starts to affect the filament shape for *D* = 0.6 mm.

A slight variation in layer height induced by gravity that is undetected by the nozzle and compounds at each layer, may also explain the variation in aspect ratio for such filaments. With the substrate at a lower position than expected by the nozzle, the plastic deformation induced by the nozzle will be lower than for the previous layers, leading to higher chances of defects and changes in the aspect ratio of the filament.

Another set of differences can be explained thanks to the preferred direction of the flow of material from the nozzle, induced by a slight angling defect of the nozzle:

* The coalescence of filaments from different layers that are not directly above one another that occurs systematically for $D \ge 0.6$ mm and often for $D = 0.4$ mm.

* The lack of coalescence between filaments observed in the first few layers for $D = 0.4$ mm.

As the nozzle is not completely parallel to the build plate, a small angle defect will impact the direction of the flow of filaments, skewing the position of filament centre toward the sides and creating an absence of coalescence in-between neighbouring filaments. In the case for coalescence between filaments that should not be in contact, the changes in geometry induced by this slight defect create now-neighbouring filaments to coalesce with, impacting the voids' shape and length.

The last set of differences seen pertains to:

* The value of the longitudinal neck growth for non-zero elements of the first few layers being lower than expected.

* The changes observed for the longitudinal neck growth for different layers of the printed test piece.

This last set of differences is explained thanks to the changes in conduction that occur during the build of the part. In the first layers, the heated zone of the substrate by the filament may be longer than the substrate itself. However, the printing bed cools the material to a given value, making the filament cool faster than expected, as the heat given to the substrate is absorbed by the printing bed. After the first few layers, a plateau is reached and the longitudinal neck growth of the following layers shows less variation once the zone heated by the filament is farther away from the printing bed.

### 4.3.2. Limits of the Models

The models, either describing the sintering of spheres or cylinders, are able to give an approximation of the values of the measured longitudinal neck growth of a print when compared with compatible data, but they are not able to account for all experimental variations.

This approximation is sometimes too high. This can be explained by the combination of two factors:

* For $D = 0.6$ mm, the variation in temperature resulting in the variation in the model prediction to values around and above the layer height of the print $H = 0.3$ mm. It is not represented in the experimental data because of the absence of voids for a large part of the observed section, resulting in no longitudinal neck growth to compare the model prediction to.

* The temperature inside a filament is not entirely constant, so the temperature at the points for contact may be lower than modelled, leading to smaller longitudinal neck growth measured.

Still, the majority of measurements of the neck growth are found between the upper limit given by the modelling of the sintering of cylinders and the lower limit given by the modelling of the sintering of spheres.

The variations in the longitudinal neck growth between layers are also not shown by the models, as the variations induced by the proximity of the printing bed and the detected
aspect ratio variations are not studied. In general, unexpected experimental variations are
not predicted by the model.

The variations in theoretical longitudinal neck growth recorded when varying the cooling
temperature of the print were not measured experimentally. The experimental longitudinal
neck growth shows little variation in absolute value with the cooling temperature, but there is
variation in the layer number where the maximum longitudinal neck growth is recorded.

During the deposition, two phenomena occur at the same time:

* The heating of the neighbouring filament by the deposited filament.

* The cooling of both filaments by conduction and convection.

These two phenomena delimit a heated zone in which the sintering between filaments
can occur. The shape of the heated zone as a result of the deposition of a filament is a
strong indicator of which spherical or cylindrical model to use. If the cooling of the part is
the dominating phenomenon, the shape of the heated zone will tend towards the sphere.
Otherwise, the shape of the heated zone will tend towards the cylinder. Here, the spherical
model shows a better approximation of the longitudinal neck growth than the cylindrical
model, and the shape of the heated zone is closer to the sphere than to the cylinder.

While the used models provide an interesting approximation towards the solution, here
are a few phenomena not considered that may have an impact on the longitudinal neck growth:

* The gradient of the temperature inside the filament and the dynamic changes in temperature resulting from such considerations.

* The effect of the previous layers on the actual layer.

* The consideration of voids of various shapes inside the mesostructure.

**5. Conclusions**

A numerical model of the coalescence between filaments was applied to PEKK for larger
filament diameters than available in the literature. The impact of the nozzle diameter and
chamber temperature on the longitudinal neck growth was measured using image analysis.
The mesostructure variations observed were discussed and the differences in the numerical
models explained, furthering the knowledge of the behaviour of high-performance polymer
parts made by Fused Filament Fabrication. The following conclusions are given:

* While the order of magnitude of the model is correct, the experimental nature of fused filament fabrication leads to a high disparity in neck growths inside a single layer, leading also to a disparity in void size between layers.

* When the shape of the void is in accordance with the model, the model can be used to obtain a quantitative approximation of the longitudinal neck growth when knowing the printing parameters. The instances where the shape of the void is not in accordance with the model are when the absence of contact between filaments is observed or when the voids are not in the shape of diamonds, as shown in Figure 10.

* In a perfect environment, the impact of the nozzle diameter will easily be explained thanks to Equations (7) and (12). The plastic deformation of the filament by the nozzle changes the geometry of the filament, leading to disparities between the theoretical model and the experimental measures.

* The spherical model acts as the upper boundary of the porosity inside a part. As the porosity of a part is limited, a conservative approach to the calculation of the neck growth results in the higher porosity of a part, leading to stricter security margins, as discussed previously in the comparison between experimental and theoretical neck growth and shown in Figures 11 and 12.

To further the prediction of the effect of chamber temperature and nozzle diameter on
the longitudinal neck growth, a complete study on computational fluid dynamics should be
carried out with the goal of solving the Navier–Stokes equation of the coalescence of two
adjacent filaments to compare to the experimental data and validating the approximations used
in the model.
A model of the neck growth of triangle-shaped voids should be added to the model to predict the neck growth in the test piece when relevant. A diagnostic of when to use diamond- or triangle-shaped voids should be added to the model.

To link the porosity of the part to its mechanical behaviour, mechanical testing was conducted on the printed parts, and further work to link printing parameters and mechanical behaviour is in progress. Then, a framework for printing parameters, porosity, and mechanical behaviour will be available for use.

The road to the zero-waste fabrication objective of additive manufacturing for high-end polymers starts with our understanding of the porosity of the part. Computational fluid dynamics models, coupled with high-resolution cameras, automatic defect detection algorithms and machine learning may result in instant parameter optimisation during the fabrication process, leading to fewer defects and less waste.

**Author Contributions:** Conceptualisation, A.R., N.M. and C.S.; methodology, A.R. and N.M.; software, A.R.; validation, A.R. and N.M.; formal analysis, A.R. and N.M.; investigation, A.R. and G.F.; resources, A.R. and G.F.; data curation, A.R. and G.F.; writing—original draft preparation, A.R.; writing—review and editing, N.M. and F.C.; visualisation, A.R.; supervision, F.C. and N.M.; project administration, T.J.; funding acquisition, T.J. All authors have read and agreed to the published version of the manuscript.

**Funding:** This project has received funding from the European Union Horizon 2020 research and innovation program under grant agreement No 101007022.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data available upon request.

**Acknowledgments:** The authors acknowledge the contribution and support of the Horizon 2020 European Project DOMMINIO.

**Conflicts of Interest:** The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

**Abbreviations**

The following abbreviations are used in this manuscript:

<table><tr><td>FFF</td><td>Fused Filament Fabrication</td></tr><tr><td>PAEK</td><td>polyaryletherketone</td></tr><tr><td>PEKK</td><td>polyether ketone ketone</td></tr><tr><td>PEEK</td><td>polyether ether ketone</td></tr><tr><td>Bi</td><td>biot number</td></tr><tr><td>H</td><td>layer height</td></tr><tr><td>D</td><td>nozzle diameter</td></tr><tr><td>T<sub>inf</sub></td><td>cooling temperature (bed and ambient)</td></tr><tr><td>t<sub>c</sub></td><td>coalescence time</td></tr><tr><td>θ<sub>cyl</sub></td><td>theoretical half angle of coalescence for cylinder-based modelling</td></tr><tr><td>θ<sub>sph</sub></td><td>theoretical half angle of coalescence for spherical modelling</td></tr><tr><td>y<sub>cyl</sub></td><td>theoretical longitudinal neck growth for cylinder-based modelling</td></tr><tr><td>y<sub>sph</sub></td><td>theoretical longitudinal neck growth for spherical modelling</td></tr><tr><td>T<sub>g</sub></td><td>glass transition temperature of PEKK-A</td></tr><tr><td>PEKK-A</td><td>amorphous PolyEther Ketone Ketone polymer</td></tr><tr><td>°C</td><td>degree Celsius</td></tr><tr><td>a<sub>0</sub></td><td>initial radius of sphere</td></tr><tr><td>a</td><td>in Section 2.2.1: radius of sphere during coalescence, in Section 2.2.2: semi-major axis of elliptic cylinder</td></tr></table>
<table><tr><td>a<sub>f</sub></td><td>radius of sphere after coalescence</td></tr><tr><td>S</td><td>area of section of two spheres during coalescence</td></tr><tr><td>W<sub>V</sub></td><td>work of viscosity forces</td></tr><tr><td>W<sub>S</sub></td><td>work of surface tension</td></tr><tr><td>r<sub>0</sub></td><td>initial radius of circular cylinder built from elliptic cylinder</td></tr><tr><td>b</td><td>semi-minor axis of elliptic cylinder</td></tr><tr><td>s</td><td>standard deviation</td></tr><tr><td>Q<sub>1</sub></td><td>first quartile of longitudinal neck growth</td></tr><tr><td>Q<sub>3</sub></td><td>third quartile of longitudinal neck growth</td></tr></table>

References

1. Pereira, T.; Kennedy, J.V.; Potgieter, J. A comparison of traditional manufacturing vs additive manufacturing, the best method for the job. *Procedia Manuf.* **2019**, *30*, 11–18. [CrossRef]

2. Rajaguru, K.; Karthikeyan, T.; Vijayan, V. Additive manufacturing—State of art. *Mater. Today Proc.* **2020**, *21*, 628–633. [CrossRef]

3. Ahmadifar, M.; Benfriha, K.; Shirinbayan, M.; Tcharkhtchi, A. Additive Manufacturing of Polymer-Based Composites Using Fused Filament Fabrication (FFF): A Review. *Appl. Compos. Mater.* **2021**, *28*, 1335–1380. [CrossRef]

4. Bikas, H.; Stavropoulos, P.; Chryssolouris, G. Additive manufacturing methods and modelling approaches: A critical review. *Int. J. Adv. Manuf. Technol.* **2016**, *83*, 389–405. [CrossRef]

5. Popescu, D.; Zapciu, A.; Amza, C.; Baciu, F.; Marinescu, R. FDM process parameters influence over the mechanical properties of polymer specimens: A review. *Polym. Test.* **2018**, *69*, 157–166. [CrossRef]

6. Nachal, N.; Moses, J.A.; Karthik, P.; Anandharamakrishnan, C. Applications of 3D Printing in Food Processing. *Food Eng. Rev.* **2019**, *11*, 123–141. [CrossRef]

7. Ranasinghe, K.; Sabatini, R.; Gardi, A.; Bikkahalli, S.; Kapoor, R.; Fahey, T.; Thangavel, K. Advances in Integrated System Health Management for mission-essential and safety-critical aerospace applications. *Prog. Aerosp. Sci.* **2022**, *128*, 100758. [CrossRef]

8. Huang, S.H.; Liu, P.; Mokasdar, A.; Hou, L. Additive manufacturing and its societal impact: A literature review. *Int. J. Adv. Manuf. Technol.* **2013**, *67*, 1191–1203. [CrossRef]

9. Wong, K.V.; Hernandez, A. A Review of Additive Manufacturing. *Int. Sch. Res. Not.* **2012**, 2012, e208760. [CrossRef]

10. Fu, X.; Lin, Y.; Yue, X.J.; Xun, M.; Hur, B.; Yue, X.Z. A Review of Additive Manufacturing (3D Printing) in Aerospace: Technology, Materials, Applications, and Challenges. In *Proceedings of the Mobile Wireless Middleware, Operating Systems and Applications*; Tang, D., Zhong, J., Zhou, D., Eds.; Innovations in Communication and Computing; Springer: Cham, Switzerland, 2022; pp. 73–98. [CrossRef]

11. Tao, Y.; Kong, F.; Li, Z.; Zhang, J.; Zhao, X.; Yin, Q.; Xing, D.; Li, P. A review on voids of 3D printed parts by fused filament fabrication. *J. Mater. Res. Technol.* **2021**, *15*, 4860–4879. [CrossRef]

12. Das, A.; Chatham, C.A.; Fallon, J.J.; Zawaski, C.E.; Gilmer, E.L.; Williams, C.B.; Bortner, M.J. Current understanding and challenges in high temperature additive manufacturing of engineering thermoplastic polymers. *Addit. Manuf.* **2020**, *34*, 101218. [CrossRef]

13. Zawaski, C.; Williams, C. Design of a low-cost, high-temperature inverted build environment to enable desktop-scale additive manufacturing of performance polymers. *Addit. Manuf.* **2020**, *33*, 101111. [CrossRef]

14. Rashid, A.A.; Koç, M. Fused Filament Fabrication Process: A Review of Numerical Simulation Techniques. *Polymers* **2021**, *13*, 3534. [CrossRef]

15. Sun, Q.; Rizvi, G.; Bellehumeur, C.; Gu, P. Effect of processing conditions on the bonding quality of FDM polymer filaments. *Rapid Prototyp. J.* **2008**, *14*, 72–80. [CrossRef]

16. Bellehumeur, C.; Li, L.; Sun, Q.; Gu, P. Modeling of Bond Formation Between Polymer Filaments in the Fused Deposition Modeling Process. *J. Manuf. Process.* **2004**, *6*, 170–178. [CrossRef]

17. Pokluda, O.; Bellehumeur, C.T.; Vlachopoulos, J. Modification of Frenkel’s model for sintering. *AIChE J.* **1997**, *43*, 3253–3256. [CrossRef]

18. Shahriar, B.B.; France, C.; Valerie, N.; Arthur, C.; Christian, G. Toward improvement of the properties of parts manufactured by FFF (fused filament fabrication) through understanding the influence of temperature and rheological behaviour on the coalescence phenomenon. *AIP Conf. Proc.* **2017**, *1896*, 040008. [CrossRef]

19. Lepoivre, A.; Levy, A.; Boyard, N.; Gaudefroy, V.; Sobotka, V. Coalescence in fused filament fabrication process: Thermodependent characterization of high-performance polymer properties. *Polym. Test.* **2021**, *98*, 107096. [CrossRef]

20. Jiang, S.; Yun, C.; Huang, X.; Zhao, C. Investigation on the Forming Process of Polylactic Acid in Material Extrusion Additive Manufacturing Technique. *SSRN* **2023**. [CrossRef]

21. Fang, L.; Yan, Y.; Agarwal, O.; Yao, S.; Seppala, J.E.; Kang, S.H. Effects of Environmental Temperature and Humidity on the Geometry and Strength of Polycarbonate Specimens Prepared by Fused Filament Fabrication. *Materials* **2020**, *13*, 4414. [CrossRef]

22. Bhalodi, D.; Zalavadiya, K.; Gurrala, P.K. Influence of temperature on polymer parts manufactured by fused deposition modeling process. *J. Braz. Soc. Mech. Sci. Eng.* **2019**, *41*, 113. [CrossRef]

23. Polychronopoulos, N.D.; Vlachopoulos, J. The role of heating and cooling in viscous sintering of pairs of spheres and pairs of cylinders. *Rapid Prototyp. J.* **2020**, *26*, 719–726. [CrossRef]
24. Weisstein, E.W. Circular Segment. From MathWorld—A Wolfram Web Resource. Available online: https://mathworld.wolfram.com/CircularSegment.html (accessed on 24 July 2024).

25. ISO 527-2:2012; Plastics—Determination of Tensile Properties. ISO: Geneva, Switzerland, 2012.

26. Zemcik, O.; Sedlak, J. Application of Linear Optimization on Parameters of 3D FDM Print. *Teh. Vjesn.* **2019**, *26*, 1164–1170. [CrossRef]

27. Ghorbani, J.; Koirala, P.; Shen, Y.L.; Tehrani, M. Eliminating voids and reducing mechanical anisotropy in fused filament fabrication parts by adjusting the filament extrusion rate. *J. Manuf. Process.* **2022**, *80*, 651–658. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.