1. Introduction
Research studies show that the morphology of nonwoven filter material has an important influence on its filtration performance [
1,
2]. Therefore, it is of great significance to measure the morphological parameters of fiber in nonwoven filter material to clarify the relationship between the structure and filtration performance and realize structural optimization and performance improvements.
Currently, the finite element method (FEM) has been widely used to predict the mechanical and filtration properties of fiber filtration materials. Elnaz Saberi et al. [
3] simulated the puncture behavior of nonwoven meso-structures with finite element methods based on the hyperelastic model, and they observed that the nonwoven fabric with higher weights showed higher puncture stress. Artem Kulachenko et al. [
4] established a three-dimensional random fiber mesh finite element model using deposition technologies and simulated the deformation and failure behavior of the fiber mesh. Alvaro Ridruejo et al. [
5] used the finite element method to simulate the mechanical behavior of glass fiber felt before fracture. N Chen et al. [
6] modeled the three-dimensional geometric model of polypropylene nonwoven meso-structures and predicted the bond performance of the model using finite element simulation. Hou X.N et al. [
7] established a discontinuous finite element model, studied the deformation mechanism of nonwoven samples in two main directions, and numerically analyzed the influence of discontinuous and nonuniform fiber mesh and different joint point arrangements. In finite element analysis, the effectiveness of the simulation results heavily depends on the established meso-structure model. The model that is more similar to the real internal structure of the nonwoven filter material is higher in reference value. Therefore, a real 3D model of nonwoven filter materials is rather crucial.
However, due to the complex internal structure and disordered spatial arrangement of fibers within nonwoven fabrics, accurately extracting each fiber’s morphology for modeling is difficult. Optical microscopy is one of the most important methods in fiber detection because it can measure fiber morphological characteristics quickly and accurately without introducing damage at a certain magnification [
8,
9]. Image processing technology based on an optical microscope is very important for establishing a 3D model of nonwoven fabric. Jing et al. [
10] calculated fiber orientations and length distributions by using the 3D model of a fiber net of wool nonwoven material established by micro-CT. Huang et al. [
11] explored the 3D skeleton image technology for establishing metal fiber mesh and calculated the morphological parameters such as the curvature and orientation of fiber fragments using skeleton tracking. Yu et al. [
12] obtained a microscopic image of a continuous section of the prefabricated fiber net of the needled carbon-carbon composite by using a layer removal method, and they established a 3D microscopic finite element model after image segmentation and noise reduction. Although the above method retains the real shape data of fiber in the nonwoven material fiber net, it can obtain an accurate model and has great development potential [
12]. However, due to the large fiber slenderness ratio and fiber bending, fibers within nonwoven fabrics inevitably involve lamination, interlacing, and entanglement. Therefore, the following two problems will be caused in the image acquisition process. Firstly, the three-dimensional stacking of fibers will cause the depth of the nonwoven fabric to exceed the microscopic system’s depth of field, which is the so-called “multi-focal plane” phenomenon [
13]. Secondly, due to interlacing and entanglement among fibers, some fibers parts are blocked, which causes difficulty in accurately extracting their complete morphology, thus reducing the accuracy of subsequent image detection results [
11].
The reported studies show that image fusion technology can effectively expand the depth of field of the microscopic system [
14,
15]. Image fusion utilizes the focusing layers of different depths, which can solve the “multi-focal plane” problem. However, the obtained 3D image still exhibits the occlusion of targets since depth information is still lost [
16]. Therefore, to establish the real 3D meso-structure of nonwoven fabrics, the problems overlapping the occlusion of the fiber system need to be completely solved.
The optical microscope cannot directly obtain the depth information of the target. Some scholars collect multiple sequence images by controlling the continuous movement of the object platform in the vertical direction and reconstructing the 3D surface of the target by focusing on the depth information [
17,
18,
19]. Based on this idea, this paper explores the method of reconstructing the surface depth image of the nonwoven fabric, thus extracting 3D fiber information. Due to the mutual occlusion of fibers, an optical microscope is limited to acquiring surface depth. It cannot directly obtain the spatial morphology and three-dimensional coordinates of the shielded part of fibers. Therefore, a “point-line-volume” modeling method based on 2D sequence images is proposed in this paper. A 3D model similar to the real meso-structure of nonwoven filter materials is established based on a series of 2D images (as shown in
Figure 1c).
This work aims to develop an effective method to establish the real meso-structure model of nonwoven fabric, thereby predicting the filtration performance with the help of FEM. The processing steps of sequential images are discussed in
Section 2. Then, the meso-structure of the nonwoven fabric is reconstructed in
Section 3.
Section 4 introduces research on finite element simulations. In
Section 5, the filtration experiments were carried out. Finally, the conclusions were obtained in
Section 6.
2. Image Processing Method
The sample used in this paper comprises a piece of polypropylene spun-bond nonwoven fabric with a weight of 18 g/m
3 and a thickness of 0.08 mm. Seventy-one images were collected by an optical microscope under different focal positions. Parts of the captured images are shown in
Figure 1a.
2.1. Image Processing
Before the image fusion process, a bilateral filtering model was adopted to denoise the captured images.
The bilateral filter is a nonlinear filter with edge-preserving denoising properties. The bilateral filter optimizes the weighted coefficient of the Gaussian filter and convolves with the gray information of the image’s pixels. In the bilateral filter, the pixel values in the region near the pixel point are weighted by a certain combination. The calculation formula is shown in Equation (1):
Here, (k, l) represents the central coordinates of the filter window, and (i, j) indicates the coordinates of other coefficients of the filter window. f (i, j) denotes the pixel value of the image at the point (i, j), and f (k, l) is the pixel value of the point (k, l); σ is the standard deviation of the Gaussian function. The weighted coefficient u (i, j, k, l) is the product of proximity factor d and brightness factor similarity factor r in pixel space.
2.2. Image Sharpness Evaluation
Since the spun-bonded nonwoven images selected in this paper have rich details and edge information, to obtain the 3D coordinate information of fibers in the nonwoven fiber, a sharpness evaluation algorithm a regional gradient variance algorithm is adopted [
20].
The algorithm uses a 7 × 7 sharpness window to evaluate the sharpness of the central pixel centering upon the pixel for which its sharpness needs to be measured, which can be expressed as Equation (2):
where
Tn (
a,
b) represents the gradient of pixels (
a,
b) in the 7 × 7 sharpness window,
Un (
i,
j) represents the average gradient of all pixels in the 7 × 7 sharpness window, and
Vn (
i,
j) represents the sharpness of the center pixel of the 7 × 7 sharpness window in the
ith focused image layer.
Gradient
Tn (
a,
b) is calculated by the Prewitt operator, and the vertical and horizontal directions of the Prewitt operator are first-order differential operators, as shown in Equation (3):
where |
a′
b′
c′| and |
d′
d′
d′| are Prewitt operators, and
represents convolution calculations.
2.3. Image Fusion
In the image sequence, the focus degree of the single fiber in different spatial positions is different, so extracting the single fiber in one captured image layer is difficult. Therefore, it is necessary to fuse the image sequence to obtain a full clear image of the single fiber.
For sample images, the sharpness of each pixel point in each focused image is calculated using Equation (4) and recorded in the 3D coordinate matrix
M(
i,
j,
n), where
i is the image’s pixel width,
j is the image’s height, and
n is the number of the image. This method compares all image frames (
n = 1, 2, …, 71) in the same position (
i,
j), calculates the focus image frame
k with the highest resolution (
i,
j), which is denoted as
H(
i,
j,
k), and stores (
i,
j,
k) coordinate information (
x1,
y1,
z1), as shown in Equation (4):
As shown in
Figure 1a, the fiber focus degree in different image sequences varies, and the fibers interlock with one another, making it difficult for a single fiber to be completely segmented. The multi-focus image fusion algorithm adopted in this paper is the multi-wavelet algorithm. This algorithm carries out a fusion method for high-frequency components based on regional direction and energy weighting, and the fusion rules are described in Equation (5):
In the formula, and correspond to the maximum and minimum values of the energy matching degree, respectively. corresponds to the three directions of 0, 45, and 90. and correspond to the j-grade decomposition coefficients of the original image A and B, respectively. is the final fuse image.
2.4. Fiber Fineness Measurement
In nonwovens, fibers are overlapped and distributed randomly in three dimensions, which makes it difficult to measure the fiber’s diameter in nonwovens materials. This paper uses digital image processing technology to measure fiber diameter in fiber mesh. The fused image in
Figure 1a shows all fibers in the net mesh. The pixel scale in the figure is 0.01 mm for every 7 pixels. Ten fiber segments were randomly selected from the fusion image to calculate the diameter of a single fiber, and diameter lines were drawn across the parallel edges of each fiber segment. The fiber diameter was obtained by measuring the diameter line of each fiber segment and taking the average value. The calculation is shown in Equation (6):
where
N represents the number of pixels and
ε is the pixel scale.
The diameter of the fiber in the spun-bonded nonwoven fabric is randomly distributed. To reduce the calculation’s cost, the diameters of all fibers within one nonwoven fabric are assumed to be constant and thus set as
d, as obtained from Equation (6). Finally, we calculated that the diameter of the fibers in the spun-bonded nonwovens is 0.025 mm.
Figure 1 shows the image fusion process and 3D coordinate extraction.
4. CFD Simulation Research
4.1. Model Meshing and Conditional Setting
This paper used computational fluid dynamics to analyze the 3D model and establish the numerical model. Before the built 3D model is imported into Ansys Fluent for numerical simulation, the 3D model needs to be meshed. The ICEMCFD components were chosen to be contained in the Ansys workbench for 3D model mesh processing. Grid processing was selected to discretize the area of the simulation calculation, which is also the basic link of numerical simulation calculations. The local mesh is shown in
Figure 3a.
When the meshed files were imported into Fluent for numerical analysis, the boundary conditions of the calculation region needed to be set, as shown in
Figure 3b. In the CFD-coupled calculation, the airflow inlet adopted the velocity-inlet condition, and the outlet adopted the pressure-outlet condition. In the text, the particle size was set between 0.5 μm and 0.7 μm. To ensure the stability of airflow entering the computational domain, the distance between the model’s particle injection surface and the upper surface was set to 2, which is the same as the distance between the particle’s escape surface and the lower surface. When the gas phase calculation reached a stable state, the particles are injected into the calculation domain along with the airflow’s direction. If the particles moved beyond the calculation area or encountered the boundary of the calculation domain. It would be recorded as an escape.
Figure 3c shows the particle trajectory in the model.
4.2. Governing Equation
The fluid control equation follows the law of mass conservation and energy conservation of the local mean variable. Gidaspow et al. [
22]. proposed two models, A and B, which are suitable for the CFD-DEM simulation of gas-solid two-phase flows. Model A assumes that the pressure drop of the flow is borne by both the solid and gas phases, while model B assumes that the pressure drop is borne only by the gas phase. Since the gas passing through the fibrous medium exhibits laminar flows, model A with laminar flow is selected in this paper to solve the equation, and the continuity equation and momentum equation are expressed as Equations (8) and (9):
where
ε,
u,
ρf,
P, and
τ represent the tensors of porosity, average fluid velocity, three-dimensional density, pressure, and viscous fluid stress, respectively.
Ffp is the total interaction forces between particles in grid cells, including the object resistance, pressure gradient force, etc.
where
fd,i is the force acting on particles, and ∆
V and
n denote the volume of grid computing cells and the number of particles within the volume of each grid cell.
In the process of particle filtering, particles and fibers produce an interaction between particles and gases so that momentum and energy exchange occur between the particles. This paper studies the contact model used in the particle: The solid nonlinear softball model. In the solid two-phase flow, according to Newton’s second law, the translation and rotation motion equation of the particle can be obtained, and the expression is as follows:
where
mp,
Up,i,
fn,i, and
Ti represent the mass, velocity and angular velocity, tangential impact force, and the impact moment of the particulate matter, respectively.
Ip denotes the inertial moment of the particle, and
F′
fp,i and
Fnc,i denote fluid and non-contact forces.
To improve the accuracy of the numerical simulation results, the SIMPLE algorithm was selected in Ansys Fluent to solve the coupling process between the velocity and pressure fields.
4.3. Numerical Simulation
The simulated particles move through the 3D model under the action flow field and pressure in the numerical simulation process. The simulation calculation will judge whether the particles interact with the fiber or are captured by the fiber. In this paper, a Lagrange multiphase discrete phase model (DPM) in Fluent was selected to track and solve the motion trajectory of the particles.
First, the material property of particles is set as anthracite particles, the velocity inlet of the fluid is set as the injection surface of the particles, and the diameter of the particles is set as five gradients: 0.5 μm, 1 μm, 2.5 μm, 5 μm, and 7 μm. The velocity of the inlet is set according to the flow of fluid, and it is divided into five gradients: 0.2063 m/s, 0.2652 m/s, 0.3242 m/s, 0.3831 m/s, and 0.4421 m/s. The collision between the emitted particles and the fibers in the 3D model is considered when the quantitative particles are launched at the same inlet velocity. The difference between particles at the launch inlet and particles escaping from the pressure outlet is the particle numbers captured by the fibers; the particles touching the boundary of the computational domain are also regarded as an escape. During the simulation calculation, if the distance between the particle and the fiber surface in the 3D model is less than or equal to the particle’s size, the particle is considered to be captured by the fiber.
Figure 3d shows the deposition process of the particles on the fiber. The movement track of the particles can be observed from the side view. It can be observed from the detailed figure on the upper right that the particle’s sizes vary. Particles with larger particle sizes are easily captured by fibers, while particles with smaller particle sizes easily escape. In addition, the motion of individual particles is shown on the lower right side of the graph.
Figure 3e shows the movement track of the particle in the model. It can be observed that the direction changes during the moving process in the flow field, which is more consistent with the movement law of the particle in the actual filtration experiment.
Figure 3b shows the filtration efficiency of the model for particles with 0.5 μm, 1 μm, 2.5 μm, 5 μm, and 7 μm particle sizes at the inlet velocity of 0.2063 m/s, 0.2652 m/s, 0.3242 m/s, 0.3831 m/s, and 0.4421 m/s. The filtering efficiency,
η, is defined as follows:
where
Cout represents the number of particles before filtering, and
Cin represents the number of particles after filtering.