# Describing and Modeling Rough Composites Surfaces by Using Topological Data Analysis and Fractional Brownian Motion

Antoine Runacher, Mohammad-Javad Kazemzadeh-Parsi, Daniele Di Lorenzo, Victor Champaney, Nicolas Hascoet, Amine Ammar, Francisco Chinesta

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

## To cite this version:

Antoine Runacher, Mohammad-Javad Kazemzadeh-Parsi, Daniele Di Lorenzo, Victor Champaney, Nicolas Hascoet, et al. Describing and Modeling Rough Composites Surfaces by Using Topological Data Analysis and Fractional Brownian Motion. Polymers, 2023, 15 (6), pp.1449. 10.3390/polym15061449. hal-04093566

HAL Id: hal-04093566

https://cnam.hal.science/hal-04093566

Submitted on 10 May 2023

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>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Describing and Modeling Rough Composites Surfaces by Using Topological Data Analysis and Fractional Brownian Motion

Antoine Runacher<fnref n="1" />,<fnref n="2" /> ID, Mohammad-Javad Kazemzadeh-Parsi<fnref n="3" /> ID, Daniele Di Lorenzo<fnref n="1" />,<fnref n="4" />, Victor Champaney<fnref n="1" />, Nicolas Hascoet<fnref n="1" />, Amine Ammar<fnref n="3" /> ID and Francisco Chinesta<fnref n="1" />,<fnref n="4" />,*

**Abstract:** Many composite manufacturing processes employ the consolidation of pre-impregnated preforms. However, in order to obtain adequate performance of the formed part, intimate contact and molecular diffusion across the different composites' preform layers must be ensured. The latter takes place as soon as the intimate contact occurs and the temperature remains high enough during the molecular reptation characteristic time. The former, in turn, depends on the applied compression force, the temperature and the composite rheology, which, during the processing, induce the flow of asperities, promoting the intimate contact. Thus, the initial roughness and its evolution during the process, become critical factors in the composite consolidation. Processing optimization and control are needed for an adequate model, enabling it to infer the consolidation degree from the material and process features. The parameters associated with the process are easily identifiable and measurable (e.g., temperature, compression force, process time, ...). The ones concerning the materials are also accessible; however, describing the surface roughness remains an issue. Usual statistical descriptors are too poor and, moreover, they are too far from the involved physics. The present paper focuses on the use of advanced descriptors out-performing usual statistical descriptors, in particular those based on the use of homology persistence (at the heart of the so-called topological data analysis—TDA), and their connection with fractional Brownian surfaces. The latter constitutes a performance surface generator able to represent the surface evolution all along the consolidation process, as the present paper emphasizes.

**Keywords:** topological data analysis—TDA; composite consolidation; rough surfaces; fractional Brownian surfaces

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

check for
updates

**Citation:** Runacher, A.;
Kazemzadeh-Parsi, M.-J.; Di Lorenzo,
D.; Champaney, V.; Hascoet, N.;
Ammar, A.; Chinesta, F. Describing
and Modeling Rough Composites
Surfaces by Using Topological Data
Analysis and Fractional Brownian
Motion. *Polymers* **2023**, *15*, 1449.
https://doi.org/10.3390/
polym15061449

Academic Editors: Abderrahmane
Ayadi, Patricia Krawczak and Chung
Hae Park

Received: 28 January 2023

Revised: 9 March 2023

Accepted: 13 March 2023

Published: 14 March 2023

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

**Copyright:** © 2023 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

Many manufacturing processes are based on the consolidation of composites' pre-
impregnated preforms that are available in the form of sheets or tapes. In both cases, the
consolidation requires putting a new sheet (or tape) in contact with the one already laid,
which now constitutes the so-called substrate [1]. If the temperature at the sheet-substrate
contact level is high enough, the polymer viscosity becomes low enough to ensure the
asperities flow under the externally applied pressure (from a press or from a roller in the
case of automated tape laying as sketched in Figure 1). Thus, the initial asperities flow and
spread, inducing the roughness smoothing and the increase in the intimate contact. As soon
as the intimate contact occurs, at a temperature ensuring molecular mobility and during a
time period long enough to ensure molecular reptation across the interface, consolidation
is attained, and ideally bulk properties are recovered at the sheet interfaces.

<footnote n="1">¹ PIMM Laboratory and ESI Group Chair, Arts et Metiers Institute of Technology, 151 Boulevard de Hopital, 75013 Paris, France</footnote>

<footnote n="2">² IPC, 2 Rue Pierre et Marie Curie, 01100 Bellignat, France</footnote>

<footnote n="3">³ LAMPA Laboratory and ESI Group Chair, Arts et Metiers Institute of Technology, 2 bd du Ronceray, 49035 Angers, France</footnote>

<footnote n="4">⁴ ESI Group, 3 Rue Saarinen, 94150 Rungis, France</footnote>

<footnote>* Correspondence: francisco.chinesta@ensam.eu</footnote>
Analysis and modeling of the roughness evolution are critical points for evaluating the quality of parts [2], to predict the remaining roughness that will result in a residual porosity in the composited laminate [3], and to evaluate tribological [4] or mechanical properties [5]. Other processes, such as machining, can induce noticeable roughness in the parts' surfaces [6].

<figure><img src="image_5.png" /><figcaption>**Figure 1.** The automated tape placement (ATP) composite manufacturing: (a) physical phenomena; (b) processing.</figcaption></figure>

The design, optimization and control of the process involves description, modeling and simulation issues [7], revisited below.

## 1.1. Rough Surface Representation

Different approaches have been proposed to represent a rough surface in view of consolidation analyses. The simplest consists of representing the asperities as rectangular elements of a given height and width [8,9]. To consider more realistic representations, fractal surfaces have also been employed, where several characteristic dimensions are hierarchically present [10–15]. Finally, the most realistic option consists of acquiring the real surface itself, on which physics-based models can be applied. The use of wavelet representations is a valuable route, because of the fact that they favor a multi-resolution description, and that a simple description consisting of rectangles is easily manipulable when, for example, the Haar wavelet is considered [7].

All these descriptions enable surface representation for solving realistic physical models on them; however, as discussed later, they do not allow to easily perform surface comparisons or classifications. The latter, in general, require the use of adequate metrics, different to the usually considered Euclidean metrics. The use of dynamic time wrapping [16] could help in some cases, but it remains too macroscopic when different description scales co-exist (fine and coarse surface features), as is the case when considering micro and macro roughness.

When looking for more intrinsic descriptors, extracting roughness features that could serve to compare and classify surfaces, statistical descriptors [17–21] represent a first natural choice. However, in many cases such descriptors are quite far from the physics involved in the consolidation process, and, even if they are adequate for describing roughness from a geometrical point of view, they seem limited in representing the physics-based surface properties as well as the time-evolution during processing.

Indeed, there exist many advanced and richer descriptors of microstructures, such as the ones reported in a study by Torquato [22]. Richer descriptions should retain much more information but at the same time be more intrinsic than when one proceeds on the real surfaces, as previously mentioned. Thus, real surfaces could be expressed in alternative spaces by using well-experienced transformations, such as Fourier, discrete cosine transform, wavelets, ..., or by learning the appropriate transformation such as *Code2Vector* [23] performs, for instance.

Recently, alternative transformations emphasizing topological persistence, which are at the heart of the so-called TDA—topological data analysis (revisited later), are emerging and proving their value [24–29]. The so-called persistence diagram *PD* and its associated
persistence image *PI* describe the topology of data in a very convenient way, such that surfaces exhibiting the same topology result in a similar *PI*. Persistence images, being defined in a vector space, can be easily compared by using adequate metrics (e.g., Euclidean, Wasserstein, ...) or manipulated by using standard machine learning technologies such as convolutional neural networks—CNN [30].

In [31] we proved that TDA can be used for classifying surfaces according to their topological content, out-performing previous studies making use of statistical descriptors [17], and proved it is possible to use that information to infer process-induced properties. By analyzing usual statistical roughness descriptors and TDA-based descriptors, respectively, in [17,31], the superiority of the latter was observed.

Another possible route consists of extracting the main modes representing these images by performing linear (principal component analysis—PCA) or nonlinear dimensionality reduction based on manifold learning [32], as successfully accomplished in [33], or by employing auto-encoders [34].

The main issue related to the use of statistical descriptors and transformation procedures is the possibility of generating synthetic surfaces according to these descriptors.

A valuable surface generator consists of exploiting the link between the random nature of roughness and Brownian motion. Generating time series, curves or surfaces from Brownian motion has been widely considered. However, Brownian motion seems, in some cases, too limited, in particular when processing occurs and the surface is strongly modified. In those cases, fractional Brownian motion seems a valuable route. Fractional motion is closely related to fractional diffusion, which, in turn, makes use of fractional calculus. Anomalous diffusion, also called non-Brownian diffusion or fractional diffusion, occurs when the mean square displacement scale with a non-integer power of time, this depending on the exponent value, defines surdiffusion (long-jumps) or subdiffusion (long-rest) [35].

Generating random surfaces by using fractional Brownian motion becomes a valuable route, and, even more interesting, this motion (as well as the resulting surfaces) is characterized from a single scalar, the so-called Hurst index. Fractional Brownian surfaces represent a timely research topic widely considered [36–40].

Thus, one could expect that modeling the Hurst index from some features, enables constructing models and makes it possible to reconstruct surfaces compatible with the features, a key route for material and process optimization and for processing control. This point represents the main goal of the present paper.

## 1.2. Original Contribution and Paper Outline

The present paper represents a step forward in the description of machine learning-based process modeling and process-induced properties.

Its main original contributions concern: (i) proving that TDA is able to discriminate fractional Brownian surfaces, the latter characterized by the Hurst index, in a very efficient manner, as soon as a regression model is developed for predicting the Hurst index from TDA-based persistence images; (ii) evaluating the effect of the composite compression on the persistence image, proving the ability of TDA to discriminate original and compressed surfaces; (iii) evidencing that the surface compression increases the Hurst index; and (iv) emphasizing the ability of fractional Brownian motion to generate synthetic surfaces with controlled topology, enabling further statistical analyses.

The paper is organized as follows. After the present introduction, the next section revisits the main methodologies employed in the present research work, in particular topological data analysis and fractional Brownian surface generation. Then, Section 3 addresses numerical analyses proving the ability of TDA to infer the Hurst index of fractional Brownian surfaces, as well as the evolution of topology and Hurst index during surface compression.
## 2. Methods

### 2.1. *Topological Data Analysis*

When TDA is applied on a time series, it transforms the time series into an image defined in a vector space, the so-called persistence image, that could be considered as a topological descriptor, enabling topological metrics to serve in comparing curves, or even to be employed as inputs in learned models.

For that purpose, and as explained in [31,41], when using a level-set-based filtration, neighboring local minima/local maxima are paired, leading to the so-called min-max pairs. Each min-max constitutes a point of the persistence diagram. The persistence diagram consists of a two-dimensional representation, where, on its horizontal axis, the minimum value of each data-pair is reported, while the associated maxima are displayed on the vertical axis. Because the maximum is always higher than the associated minimum, the axes of the persistence diagram are usually labeled as birth and death axes, respectively.

For the random profile depicted in Figure 2, (11, 14), (7, 9) and (9, 10) represent the three neighboring local minima/local maxima.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Random profile.</figcaption></figure>

Because of the nature of maximum and minimum (or birth and death) the points in the persistence diagram are located in the upper bisector defining the space $y \ge x$, with $x$ and $y$ referring to birth and death (horizontal and vertical), respectively. The persistence diagram related to the profile depicted in Figure 2 is shown in Figure 3.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Typical persistence diagram *PD*.</figcaption></figure>

The higher the distance of a point to the diagonal $x = y$, the more persistent the topology of the considered point (maximum to minimum difference). Points located close to the diagonal contain a small persistent topology, which in many cases can be associated to noise (or finer physical scales).
Another important point concerns the invariance properties that the topology provides. If we consider a point on the persistence diagram depicted in Figure 3, for example the (11, 14)-pair, this point is associated with local-minimum/local-maximum neighboring data points. All the pairs with a minimum value equal to 11, and the neighboring local maximum equal to 14, are associated with the referred (11, 14)-pair, all of them represented by the same data-point in the persistence diagram, independently of their location on the time axis, or the time delay between the minimum and maximum constituting each pair.

Thus, two identical time series, but one of them translated in the time axis with respect to the second one, are represented by the same persistence diagram. Two time series, one of them consisting of a dilation of the other, have the same persistence diagram. These invariance properties are very important when the topology is more important than the location at which these topological events take place.

When addressing heat and momentum transfer on rough interfaces, TDA seems a very promising descriptor that extracts geometrical features closely related to the physics operating on it.

When topologies are close but not identical, their associated persistence diagrams consist of two set of points. To evaluate the distance between them, appropriate metrics must be employed for comparing the distance between both point clouds. Thus, the Wasserstein distance consists of a simple and efficient choice for measuring the distance between two persistence diagrams.

Another derived representation consists of the so-called persistence image. Prior to its introduction, an intermediate transformation is needed. Instead of representing birth–death, an equivalent representation consists of representing birth–lifetime, the latter defined as the difference between the death and the birth. The main advantage of using the lifetime diagram instead of the usual persistence diagram, is that in the former points fill the whole 2D space instead of only the upper bisector ($y \ge x$). The lifetime diagram $LT$ associated with the persistence diagram $PD$ shown in Figure 3 is given in Figure 4.

<figure><img src="image_8.png" /><figcaption>Figure 4. Lifetime diagram $LT$.</figcaption></figure>

Now, the last step consists of smoothing the point representation involved in the lifetime diagram, to transform a point into a 2D function. For that purpose, and as described in [31], points in the lifetime diagram are replaced by a Gaussian function centered at each point of the lifetime diagram $LT$.

Next, the so-called persistent image $PI$, defined in a vector space is constructed. For that purpose we consider a continuous piece-wise differentiable non-negative weight function $w(x, y)$ and a bivariate normal distribution $g_{x,y}(u, v)$ centered at each point $(x, y)$, from which a function $\rho(u, v)$ can be defined at each point $(u, v)$ in the lifetime domain [31]

$$ \rho(u, v) = \sum_{(x,y) \in LT} w(x, y) g_{(x,y)}(u, v). \qquad (1) $$
Now, the domain is partitioned in a series of non-overlapping subdomains covering it, the so-called pixels $P_i$, $i = 1, 2, \dots$, where function $\rho(u, v)$ is averaged to define the persistence image $\mathcal{PI}$, according to

$$ \mathcal{PI}(P_i) = \iint_{P_i} \rho(u, v) du dv. \qquad (2) $$

A typical persistence image related to a profile larger than the one shown in Figure 2 is illustrated in Figure 5. Now, this being a quite standard image, one is tempted to use it as input in learning procedures, in particular those making use of convolutional neural networks—CNN [33], even if many other options exist. Other possibilities consist of applying dimensionality reduction, either linear (e.g., PCA) or nonlinear (e.g., auto-encoders), to facilitate the neural network training, that now makes use of smaller input data, or of applying the weights of PCA modes when using the PCA or the data mapped into the latent reduced space when considering auto-encoders.

<figure><img src="image_9.png" /><figcaption>Figure 5. Typical persistence image.</figcaption></figure>

## 2.2. From Brownian Diffusion to Anomalous Diffusion

This section revisits the derivation of Brownian diffusion and its connection with Brownian motion, viewed as a random walk. Then, in the case of anomalous diffusion, fractional Brownian motion will lead to the so-called fractional Brownian trajectories and surfaces.

### 2.2.1. Brownian Diffusion

Following Einstein's works, the particle displacement $\Delta$ in the unbounded one-dimensional axis $x$ is represented by a random variable whose probability density reads $\phi(\Delta)$.

The particle's balance considers both, the local time evolution

$$ \rho(x, t + \tau) = \rho(x, t) + \frac{\partial \rho(x, t)}{\partial t} \tau + \Theta(\tau^2), \qquad (3) $$

and the transport in space

$$ \rho(x, t + \tau) = \int_{\mathbb{R}} \rho(x + \Delta, t) \phi(\Delta) d\Delta. \qquad (4) $$

Developing $\rho(x + \Delta, t)$

$$ \rho(x + \Delta, t) = \rho(x, t) + \frac{\partial \rho(x, t)}{\partial x} \Delta + \frac{1}{2} \frac{\partial^2 \rho(x, t)}{\partial x^2} \Delta^2 + \Theta(\Delta^3) \qquad (5) $$
and taking into account the normality $\int_{\mathbb{R}} \phi(\Delta)d\Delta = 1$ and symmetry $\int_{\mathbb{R}} \Delta\phi(\Delta)d\Delta = 0$ conditions, the final result is

$$ \rho(x,t) + \frac{\partial\rho(x,t)}{\partial t} \tau = \rho(x,t) + \frac{1}{2} \frac{\partial^2\rho(x,t)}{\partial x^2} \int_{\mathbb{R}} \Delta^2\phi(\Delta)d\Delta, \quad (6) $$

or

$$ \frac{\partial\rho(x,t)}{\partial t} \tau = D \frac{\partial^2\rho(x,t)}{\partial x^2}, \quad (7) $$

with

$$ D = \frac{1}{2\tau} \int_{\mathbb{R}} \Delta^2 \phi(\Delta) d\Delta. \quad (8) $$

The integration of the diffusion equation from the localized (at $x = 0$) initial condition $\rho(x, t = 0) = \delta_0(x)$, leads to the Gaussian distribution

$$ \rho(x,t) = \frac{1}{\sqrt{4\pi Dt}} e^{-\frac{x^2}{4Dt}}, \quad (9) $$

whose main square displacement scales linearly with the time

$$ \langle x^2 \rangle = 2Dt. \quad (10) $$

### 2.2.2. Brownian Motion

The diffusion equation can be also derived from a random walk. For illustrating the derivation, we consider a grid on the x-axis, with step $\Delta x$ being the time step $\Delta t$. At each time step, each particle is assumed to jump to one of its two neighboring sites with the same probability.

Thus, the balance at the j-site (grid node) reads

$$ W_j(t + \Delta t) = \frac{1}{2} W_{j+1}(t) + \frac{1}{2} W_{j-1}(t), \quad (11) $$

where $W_j(t)$ is the probability of having the particle at site j at time t.

By developing $W_i(t + \Delta t)$, $W_{j+1}(t)$ and $W_{j-1}(t)$ according to

$$ \left\{ \begin{array}{l} W_i(t + \Delta t) = W_j(t) + \left. \frac{\partial W_j(t)}{\partial t} \right|_t \Delta t + \Theta(\Delta t^2) \\ W_{j+1}(t) = W_j(t) + \left. \frac{\partial W_j(t)}{\partial x} \right|_j \Delta x + \frac{1}{2} \left. \frac{\partial^2 W_j(t)}{\partial x^2} \right|_j \Delta x^2 + \Theta(\Delta x^3), \\ W_{j-1}(t) = W_j(t) - \left. \frac{\partial W_j(t)}{\partial x} \right|_j \Delta x + \frac{1}{2} \left. \frac{\partial^2 W_j(t)}{\partial x^2} \right|_j \Delta x^2 - \Theta(\Delta x^3) \end{array} \right. \quad (12) $$

the balance results in the usual diffusion equation

$$ \frac{\partial W}{\partial t} = D \frac{\partial^2 W}{\partial x^2}, \quad (13) $$

with D defined in the limit of infinitesimal $\Delta x$ and $\Delta t$ by

$$ D = \frac{\Delta x^2}{2\Delta t}. \quad (14) $$

### 2.3. Fractional Brownian Motion

The equivalence between Brownian diffusion and Brownian motion motivates the modeling of anomalous diffusion (which results in the fractional diffusion equation [35]) by using fractional Brownian motion.

In fractional Brownian motion, also known as fractal Brownian motion, the increments are not independent, in contrast to Brownian motion.
For that purpose, we consider the continuous Gaussian process $B_H(t)$ (the index ●$_H$ refers to the Hurst index), with $B_H(t = 0) = 0$, having a null expectation at each time, and defined by the covariance

$$ \mathbb{E}(B_H(r), B_H(s)) = \frac{1}{2} (|r|^{2H} + |s|^{2H} - |r-s|^{2H}), \qquad (15) $$

with $H = 0.5$ in the case of Brownian motion.

For generating a random trajectory we consider a set of time instants $t_1, t_2, \dots, t_n$ and construct the matrix $\mathbf{R}$ (symmetric positive-definite), whose component $R_{ij}$ given by

$$ R_{ij} = \frac{1}{2} (|t_i|^{2H} + |t_j|^{2H} - |t_i - t_j|^{2H}). \qquad (16) $$

The singular value decomposition—SVD of $\mathbf{R}$ writes $\mathbf{R} = \mathbf{V}\mathbf{D}\mathbf{V}^T$. Thus, we can define the matrix $\mathbf{S} = \mathbf{V}\mathbf{D}^{1/2}\mathbf{V}^T$.

Now, we generate $n$ values following a standard Gaussian distribution (zero mean and unit variance), grouped in vector $\mathbf{G}$. Finally, this results in the associated fractional values in vector $\mathbf{F}_H$:

$$ \mathbf{F}_H = \mathbf{S}\mathbf{G}. \qquad (17) $$

### 3. Numerical Tests

#### 3.1. Topological Data Analysis

Different pre-impregnated composite sheets, delivered by the same provider and prepared with the same materials and under the same process conditions, were analyzed, to check the existence of topological variability on their respective surfaces. A Bruker optical profilometer was used to extract the surface roughness profiles along different directions on the sheet surfaces, in particular along the direction of fibers, along the perpendicular direction to those fibers, and along long fibers that constitute the pre-impregnated composite sheets reinforcement.

Roughness profiles were measured at 17 different locations on both surfaces of each composite sheet. The locations of those measured zones and their dimensions are illustrated in Figure 6, where the fibers' orientation is also indicated.

<figure><img src="image_10.png" /><figcaption>Figure 6. Measurement locations on the composite sheets' surfaces.</figcaption></figure>

The persistence images associated with the different profiles measured along a fiber's direction and along its perpendicular direction, were obtained. A supervised random forest classifier involving 400 trees was trained to discriminate the profiles related to each of the three available composite sheets. A total of 80% of the 1020 available profiles was considered in the training, whereas the remaining 20% served to test the trained classifier performances. Figure 7 shows the associated confusion matrix related to the test set. This figure proves that the three sheet surfaces exhibited similar roughness from a topological viewpoint, and consequently the trained classifier failed to discriminate to which sheet each profile in the test set belongs.
Different classifiers were tested, all them exhibiting similar performances. To conclude on the reasons of the poor classifier performance, different persistence images extracted from the three sheets were analyzed and it was noticed that the differences between the images coming from the same sheets were similar to the differences between the images of different sheets. Thus, it can be concluded that the poor classification capabilities are due to the nature of the surfaces, instead of being a consequence of the classifiers' performance. In these analyses, the Wasserstein distance was considered for quantifying the differences between the topological description of the surfaces.

<figure><img src="image_11.png" /><figcaption>**Figure 7.** Confusion matrix related to the profiles of the three composite sheets.</figcaption></figure>

However, during composite processing, the sheet surfaces will experience a compres-
sion affecting the roughness. Even if, as commented, no difference is expected between
the different sheets, a difference is expected concerning the roughness topology before
and after compression. To prove that discrimination ability, a sample of each sheet was
heated to 200 degrees Celsius and compressed by using a pneumatic press, during 5 min,
for ensuring the flow of surface asperities.

Figure 8 compares two characteristic surface topographies before and after compres-
sion, from which a noticeably higher smoothness can be seen on the compressed surface.
Then, different rough profiles were measured by again using the optical profilometer in
both directions (with respect to the fiber orientation).

As the roughness is significantly higher in the perpendicular direction to the fiber orientation, our analyses mainly focus on profiles along that perpendicular direction.

<figure><img src="image_12.png" /><figcaption>**Figure 8.** Surface topography before (left) and after (right) compression.</figcaption></figure>

Figure 9 shows the confusion matrix associated with a random forest classifier (in-
volving 400 trees), trained from 1016 profiles and tested on 254 profiles that were excluded
from the training set. It can be concluded that TDA perfectly discriminates original and
compressed surfaces.
<figure><img src="image_13.png" /><figcaption>**Figure 9.** Confusion matrix related to the profiles of the three composite sheets.</figcaption></figure>

To better visualize the roughness-topology change induced by the compression, Figure 10 shows three persistence images, one extracted from each of the three composite sheets along the perpendicular direction to the long fibers, where no major differences are noticed between them.

These figures also reveal an almost mono-disperse topological content, which explains the localized pattern present in the persistence images.

On the other side, Figure 11 compares the roughness topologies before and after compression, where a significant difference can clearly be noticed; the topology is more poly-disperse after compression. On the original surface the asperities' height seems quite mono-disperse, with the associated localization in the persistence image. When compressing the surface, the highest asperities start to be compressed and the asperities' height starts exhibiting a larger spectrum, with the associated effect on the persistence image.

<figure><img src="image_14.png" /><figcaption>**Figure 10.** Persistence images related to three roughness profiles, each extracted from one of the three composite sheets.</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Figure 11.** Characteristic persistence images before (left) and after (right) compression.</figcaption></figure>

### 3.2. Fractional Brownian Surfaces

This section aims at generating different roughness profiles by using fractional Brownian motion, characterized by different Hurst indexes. Then, a classifier is trained with
the associated TDA surfaces, and the classifier's performance is tested on the surfaces composing the test set.

Figure 12 shows three surfaces generated from a fractional Brownian motion with $H = 0.25$, $H = 0.5$ and $H = 0.75$. No major differences are noticeable at first glance; however, a topological footprint is expected to exist, enabling an efficient classification.

<figure><img src="image_16.png" /><figcaption>**Figure 12.** Profiles generated by a fractional Brownian motion characterized by $H = 0.25$ (**left**), $H = 0.5$ (**center**) and $H = 0.75$ (**right**).</figcaption></figure>

To prove the Brownian behavior of profiles generated with $H = 0.5$, and the anomalous diffusion for $H < 0.5$ and $H > 0.5$, we generated 1000 trajectories and evaluated their mean square displacement evolution. Figure 13 superposes 30 different profiles generated with $H = 0.25$, $H = 0.5$ and $H = 0.75$. Figure 14 depicts the evolution of the mean square displacement. As expected, it can be noticed that for $H = 0.5$ (Brownian motion) the mean square displacement evolves linearly.

<figure><img src="image_17.png" /><figcaption>**Figure 13.** Population of 30 profiles generated by a fractional Brownian motion characterized by $H = 0.25$ (**left**), $H = 0.5$ (**center**) and $H = 0.75$ (**right**).</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>**Figure 14.** Mean square displacement evolution for the population of profiles generated by a fractional Brownian motion characterized by $H = 0.25$ (**left**), $H = 0.5$ (**center**) and $H = 0.75$ (**right**).</figcaption></figure>

Now, the persistence images associated with each profile are obtained, and a part of them (2400 images) used for training a random forest classifier (involving 400 trees), which is then applied for classifying 600 persistence images composing the test set.

The confusion matrix associated with the classifier, again a random forest involving 400 trees, trained from 8000 images, applied to the test set (composed of 2000 images) is shown in Figure 15 and proves the ability of the TDA to capture the different data topologies associated with the different Hurst indexes.
<figure><img src="image_19.png" /><figcaption>**Figure 15.** Confusion matrix related to the trained classifier applied to the three fractional Brownian surfaces characterized by *H* = 0.25, *H* = 0.5 and *H* = 0.75.</figcaption></figure>

In order to prove the proposed model's robustness, the procedure was repeated, but now considered 10 different Hurst indexes. Again, the confusion matrix associated with the trained classifier, shown in Figure 16, proves the TDA's robustness for representing fractional Brownian surfaces.

<figure><img src="image_20.png" /><figcaption>**Figure 16.** Confusion matrix related to the trained classifier applied to ten fractional Brownian surfaces.</figcaption></figure>

### 3.3. Hurst Index Evolution during the Surface Compression

To better understand the effect of surface compression on the Hurst index, we considered many Brownian surfaces generated by considering *H* = 0.5. Then, we removed all the positive peaks, that is, the positive part of all the profiles was removed and replaced by a null value, for emulating the flatness induced by the compression. The asperities' flattening is compulsory for consolidation to occur, because of the fact that molecular reparation needs the intimate contact that flatness facilitates. Then, a regressor (random forest involving 400 trees) operating on the persistence images of 4800 profiles, with Hurst indexes ranging in the interval (0.05, 0.75), and tested on the 1200 images composing the test set, was applied to the persistence images associated with flattened profiles. It was observed that the inferred Hurst index ranged in the interval (0.55, 0.6), proving that, as expected, higher correlations are induced in the data representing the flattened profiles.

Now, the same rationale was applied to the real composites' profiles, for comparing the Hurst index distribution before and after compression. For that purpose 250 non-compressed and other 250 compressed profiles, randomly chosen from the 1270 available, were analyzed, and their Hurst indexes inferred by using the trained regressor just introduced. The histograms related to the inferred Hurst indexes are compared in Figure 17, evidencing that the compression induces a noticeable increase in the Hurst index distribution.
<figure><img src="image_21.png" /><figcaption>**Figure 17.** Inferred Hurst index histograms for 250 non-compressed (red) and 250 compressed (blue) profiles.</figcaption></figure>

**4. Conclusions**

The present paper first proved the ability of topological data analysis to quantify the roughness evolution during composite processing involving compression, and at the same time characterize the roughness of pre-impregnated composite sheets.

Then, the link between TDA and fractional Brownian surfaces was evidenced. TDA enables classifying those surfaces according to their associated Hurst indexes.

Finally, compressed and non-compressed surfaces were characterized using a Hurst index evolution, which could represent a valuable descriptor in processing optimization and control.

**Author Contributions:** Conceptualization, A.R., F.C., A.A.; methodology, A.R., D.D.L., V.C.; experiments and data analysis, M.-J.K.-P.; software, N.H., A.R. 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 programme under grant agreement No 101007022.

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

**Informed Consent Statement:** Not applicable.

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

**Acknowledgments:** The authors acknowledge the contribution and support of: (i) the CREATE-ID ESI-ENSAM Chair; (ii) the European Union Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 956401 (XS-Meta); and (iii) the H2020 DOMMINIO project.

**Conflicts of Interest:** The authors declare no conflict of interest.
## References

1. Chinesta, F.; Leygue, A.; Bognet, B.; Ghnatios, C.; Poulhaon, F.; Bordeu, F.; Barasinski, A.; Poitou, A.; Chatel, S.; Maison-Le-Poec, S. First steps towards an advanced simulation of composites manufacturing by automated tape placement. *Int. J. Mater. Form.* **2014**, *7*, 81–92. [CrossRef]

2. Boon, Y.D.; Joshi, S.C.; Bhudolia, S.K. Filament Winding and Automated Fiber Placement with In Situ Consolidation for Fiber Reinforced Thermoplastic Polymer Composites. *Polymers* **2021**, *13*, 1951. [CrossRef] [PubMed]

3. Song, Q.; Liu, W.; Chen, J.; Zhao, D.; Yi, C.; Liu, R.; Geng, Y.; Yang, Y.; Zheng, Y.; Yuan, Y. Research on Void Dynamics during In Situ Consolidation of CF/High-Performance Thermoplastic Composite. *Polymers* **2022**, *14*, 1401. [CrossRef]

4. Pierik, E.R.; Grouve, W.J.B.; Wijskamp, S.; Akkerman, R. Prediction of the peak and steady-state ply-ply friction response for UDC/PAEK tapes. *Compos. Part A* **2022**, *163*, 107185. [CrossRef]

5. Li, C.; Fei, J.; Zhang, T.; Zhao, S.; Qi, L. Relationship between surface characteristics and properties of fiber-reinforced resin-based composites. *Compos. Part B Eng.* **2023**, *249*, 110422. [CrossRef]

6. Rajasekaran, T.; Palanikumar, K.; Latha, B. Investigation and analysis of surface roughness in machining carbon fiber reinforced polymer composites using artificial intelligence techniques. *Carbon Lett.* **2022**, *32*, 615–627. [CrossRef]

7. Leon, A.; Argerich, C.; Barasinski, A.; Soccard, E.; Chinesta, F. Effects of material and process parameters on in-situ consolidation. *Int. J. Mater. Form.* **2019**, *12*, 491–503. [CrossRef]

8. Lee, W.I.; Springer, G.S. A model of the manufacturing process of thermoplastic matrix composites. *J. Compos. Mater.* **1987**, *21*, 1057–1082.

9. Levy, A.; Heider, D.; Tierney, J.; Gillespie, J. Inter-layer thermal contact resistance evolution with the degree of intimate contact in the processing of thermoplastic composite laminates. *J. Compos. Mater.* **2014**, *48*, 491–503. [CrossRef]

10. Borodich, F.; Mosolov, A. Fractal roughness in contact problems. *J. Appl. Math. Mech.* **1992**, *56*, 786–795. [CrossRef]

11. Ganti, S.; Bhushan, B. Generalized fractal analysis and its applications to engineering surfaces. *Wear* **1995**, *180*, 17–34. [CrossRef]

12. Leon, A.; Barasinski, A.; Chinesta, F. Microstructural analysis of pre-impregnated tapes consolidation. *Int. J. Mater. Form.* **2017**, *10*, 369–378. [CrossRef]

13. Majumdar, A.; Tien, C. Fractal Characterization and simulation of rough surfaces. *Wear* **1990**, *136*, 313–327. [CrossRef]

14. Mandelbrot, B.; Passoja, D.; Paullay, A. Fractal character of fracture surfaces of metals. *Nature* **1984**, *308*, 721–722. [CrossRef]

15. Yang, F.; Pitchumani, R. A fractal Cantor set based description of interlaminar contact evolution during thermoplastic composites processing. *J. Mater. Sci.* **2001**, *36*, 4661–4671. [CrossRef]

16. Senin, P. *Dynamic Time Warping Algorithm Review*; Technical Report; University of Hawaii at Manoa: Honolulu, HI, USA, 2008.

17. Argerich, C.; Ibáñez, R.; León, A.; Abisset-Chavanne, E.; Chinesta, F. Tape surface characterization and classification in automated tape placement processability: Modeling and numerical analysis. *AIMS Mater. Sci.* **2018**, *5*, 870–888. [CrossRef]

18. Longuet-Higgins, M. Statistical properties of an isotropic random surface. *Ser. A-Math. Phys. Sci.* **1957**, *250*, 157–174.

19. Longuet-Higgins, M. The Statistical Analysis of a Random, moving surface. *Ser. A-Math. Phys. Sci.* **1957**, *249*, 321–387.

20. Sayles, R.; Thomas, T. The spatial representation of surface roughness by means of the structure function: A practical alternative to correlation. *Wear* **1977**, *42*, 263–276. [CrossRef]

21. Yaglom, A. Volume I—Basic Results. In *Correlation Theory of Stationary and Related Random Function*; Springer: New York, NY, USA, 1987.

22. Torquato, S. Statistical Description of Microstructures. *Annu. Rev. Mater. Res.* **2002**, *32*, 77–111. [CrossRef]

23. Argerich, C.; Ibanez, R.; Leon, A.; Barasinski, A.; Chinesta, F. Code2vect: An efficient heterogenous data classifier and nonlinear regression technique. *C. R. Mécanique* **2019**, *347*, 754–761. [CrossRef]

24. Carlsson, G.; Zomorodian, A.; Colling, A.; Guibas, L. Persistence Barcodes for Shapes. In Proceedings of the 2004 Eurographics/ACM SIGGRAPH Symposium on Geometry Processing, Nice, France, 8–10 July 2004.

25. Carlsson, G. Topology and Data. *Bull. Am. Math. Soc.* **2009**, *46*, 255–308. [CrossRef]

26. Chazal, F.; Michel, B. An introduction to Topological Data Analysis: Fundamental and practical aspects for data scientists. *arXiv* **2017**, arXiv:1710.04019.

27. Oudot, S.Y. *Persistence Theory: From Quiver Representation to Data Analysis*; Mathematical Surveys and Monographs; American Mathematical Society: Providence, RI, USA, 2010; Volume 209.

28. Rabadan, R.; Blumberg, A. *Topological Data Analysis For Genomics And Evolution*; Cambridge University Press: Cambridge, UK, 2020.

29. Saul, N.; Tralie, C. *Scikit-TDA: Topological Data Analysis for Python*. 2019. Available online: https://github.com/scikit-tda/scikit-tda (accessed on 12 March 2023).

30. Venkatesan, R.; Li, B. *Convolutional Neural Networks in Visual Computing: A Concise Guide*; CRC Press: Boca Raton, FL, USA, 2017.

31. Frahi, T.; Yun, M.; Argerich, C.; Falco, A.; Chinesta, F. Tape Surfaces Characterization with Persistence Images. *AIMS Mater. Sci.* **2020**, *7*, 364–380. [CrossRef]

32. Lee, J.A.; Verleysen, M. *Nonlinear Dimensionality Reduction*; Springer: New York, NY, USA, 2007.

33. Yun, M.; Argerich, C.; Cueto, E.; Duval, J.L.; Chinesta, F. Nonlinear regression operating on microstructures described from Topological Data Analysis for the real-time prediction of effective properties. *Materials* **2020**, *13*, 2335. [CrossRef] [PubMed]

34. Hinton, G.E.; Zemel, R.S. Autoencoders, minimum description length and Helmholtz free energy. In *Advances in Neural Information Processing Systems 6 (NISP 1993)*; Morgan-Kaufmann: Burlington, MA, USA, 1993; pp. 3–10.
35. Chinesta, F.; Abisset, E. *A Journey Around the Different Scales Involved in the Description of Matter and Complex Systems*; SpringerBrief: Cham, Germany, 2017.

36. Bardet, J.M.; Surgailis, D. Measuring the roughness of random paths by increment ratios. *Bernoulli* **2011**, *17*, 749–780. [CrossRef]

37. Gelbaum, Z.; Titus, M. Simulation of Fractional Brownian Surfaces via Spectral Synthesis on Manifolds. *arXiv* **2013**, arXiv:1303.6377v1.

38. Kroese, D.P.; Botev, Z.I. Spatial Process Generation. *arXiv* **2013**, arXiv:1308.0399v1.

39. Rabiei, H.; Coulon, O.; Lefevre, J.; Richard, F. Surface regularity via the estimation of fractional Brownian motion index. *IEEE Trans. Image Process.* **2020**, *30*, 1453–1460. [CrossRef]

40. Stein, M.L. Fast and Exact Simulation of Fractional Brownian Surfaces. *J. Comput. Graph. Stat.* **2002**, *11*, 587–599. [CrossRef]

41. Frahi, T.; Chinesta, F.; Falco, A.; Badias, A.; Cueto, E.; Choi, H.Y.; Han, M.; Duval, J.L. Empowering Advanced Driver-Assistance Systems from Topological Data Analysis. *Mathematics* **2021**, *9*, 634. [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.