Next Article in Journal
Algorithm for the Accelerated Calculation of Conceptual Distances in Large Knowledge Graphs
Previous Article in Journal
Modeling of Nonlinear Sea Wave Modulation in the Presence of Ice Coverage
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In-Plane Wave Propagation Analysis of Human Breast Lesions Using a Higher-Order Nonlocal Model and Deep Learning

1
Adelaide Medical School, University of Adelaide, The Queen Elizabeth Hospital, Woodville South, SA 5011, Australia
2
Robinson Research Institute, University of Adelaide, Adelaide, SA 5006, Australia
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(23), 4804; https://doi.org/10.3390/math11234804
Submission received: 27 October 2023 / Revised: 23 November 2023 / Accepted: 24 November 2023 / Published: 28 November 2023
(This article belongs to the Section Engineering Mathematics)

Abstract

:
The wave propagation characteristics of biological tissues are of high importance in improving healthcare technologies and can be used as an early clinical indicator of many diseases. However, the current mathematical models that describe the mechanical properties of biological tissues do not account for the difference in softening and hardening observed at different scales and this limits their utility in biomedical imaging. In this paper, a higher-order nonlocal model is developed to study in-plane wave propagation in healthy, benign, and cancerous breast tissues. To verify the mathematical approach, finite element simulations are conducted. Furthermore, a sequential deep neural network model of feedforward type with multiple hidden layers is developed to understand the intrinsic in-plane wave characteristics of breast tissues. The deep learning algorithm shows potential in accurately extracting the frequencies and phase velocities of breast lesions under in-plane waves even when there is a limited number of clinical samples. Using the higher-order nonlocal model, significant differences between healthy fibroglandular tissue and early breast cancer in the form of ductal carcinoma in situ have been found. The combination of nonlocal and strain gradient parameters allows for the concurrent incorporation of stiffness hardening and softening, solving the rigid-tumour–soft-cell paradox of cancer biomechanics.

1. Introduction

Understanding the biomechanical features and electromechanical response of biological tissues and biomedical components is of significance in many medical science and technology fields including implantable devices, tissue-engineered constructs, drug delivery systems and disease detection [1,2]. To increase the efficiency and durability of biomedical implantable devices such as knee joint replacement tools, cardioverter defibrillators, and dental implants, mechanical features including elasticity properties and vibrational characteristics have vital roles since they are highly associated with the overall performance of the device [3,4]. Biomechanical properties are also important in the design and manufacture of biomaterials used as the building blocks of artificial organs and tissue-engineering scaffolds [5,6]. Optimised mechanical properties lead to the smooth and robust integration of biomaterials into target tissues, reducing inflammation and enhancing the efficacy of the artificial organ or scaffold [5,7]. In drug delivery systems, acoustic analysis and fluid flow studies help researchers and clinicians gain a better understanding of therapeutic particle motion and improve drug delivery rates [8,9]. Moreover, biomechanical features are used as indicators of disease for detection purposes and in treatment monitoring to assess therapy effectiveness [10,11,12].
In clinical applications, ultrasound transducers are commonly used to create waves within tissue and the patient’s body [13]. Depending on the direction in which the ultrasound transducer is applied by the professional nurse or general practitioner, it can create in-plane, shear, and flexural waves [14,15]. If the transducer as the initial source of the waves is applied perpendicular to the human tissue, it leads to flexural waves; meanwhile, when the ultrasound transducer is utilised along the side of the tissue or human body, the in-plane wave propagation assumption is valid [14,15]. Each direction of wave propagation provides useful information about potential diseases or abnormalities. A combination of them is used in some applications for better spatial visualisation [13,15].
Elastography is a promising non-invasive biomedical imaging technique in which the biomechanical features of the target biological tissue are visualised on an image called an elastogram for disease detection or treatment efficacy assessments [15,16]. Depending on the source of the mechanical stimulation, there are several types of elastography imaging including ultrasound [17], magnetic resonance [18], and optical [12]. In all of these applications, mechanical deformation is created using a safe and non-invasive approach and then the corresponding displacements per unit length (strain) are measured. The obtained stress and strain are utilised in a mathematical model to estimate the biomechanical features of the target tissue [19]. Elastography imaging has been used in the rapid and quantitative detection of many diseases and abnormalities, especially those in which fibrosis, stiffness hardening, or/and softening are involved. These biomedical imaging applications include—but are not limited to—breast cancer [20], liver fibrosis [21], polycystic ovary syndrome [22], and prostate cancer [23].
Accurate and computationally fast mathematical models are essential in elastography imaging in order to calculate and plot intrinsic biomechanical features such as stiffness [24], Poisson’s ratio [25], and poroelastic [26] and viscoelastic properties [27]. To provide elastography mathematical models of deformation for clinical applications, different continuum mechanics-based models have been developed in recent years. Samani et al. [28] introduced a finite element approach for studying the elasticity behaviour of human breast samples for cancer detection. They conducted finite element simulations and experimentally obtained the elasticity modulus of a wide list of different pathological breast tissues including both malignant and benign conditions. It was found that the presence of cancer in human breast tissue leads to a significant increase in the elasticity modulus, which could be clinically beneficial in the detection of breast cancer. Kheirkhah et al. [29] used mathematical constraints obtained based on the principles of tissue mechanics for the regularization of displacements and strains in ultrasound elastography, incorporating factors such as tissue incompressibility and deformation compatibility. Poul et al. [30] studied different rheological models of viscoelasticity for elastography imaging applications including the Kelvin–Voigt, standard linear solid, and Kelvin–Voigt fractional derivative models. They concluded that the Kelvin–Voigt fractional derivative model performs consistently in terms of stress relaxation. Moreover, hyperelastic constitutive models [31] and porohyperviscoelastic approaches [32] have been developed for shear-wave elastography and have been investigated for tissue-mimicking materials and human liver, respectively. A doublet mechanics model [33], a convolutional neural network technique [34], and an analytical poroelastic approach [35] have also been proposed to improve tissue mechanical property visualisation. For more details about the state-of-the-art methods and techniques used to obtain the biomechanical features of biological tissues, an interested reader is referred to a recent review paper by Alekya et al. [36].
In recent years, several higher-order biomechanical models have been developed to study the deformation behaviour of biological components and tissues under mechanical forces. Xiang and Liew [37] introduced a highly accurate computational platform to analyse the transverse compression of protein microtubules using a higher-order Cauchy–Born law. Moreover, Xiang et al. [38] developed a Cauchy–Born rule of third order in conjunction with an element-free computational framework to simulate the biomechanical behaviour of microtubules. More recently, Dwairy et al. [39] applied the biphasic theory to predict interstitial fluid pressure and stress within solid tumours. They found that mechanical stress in solid tumours highly depends on the nonlinearity and constitutive model of deformation while the interstitial fluid pressure is relatively independent of the fundamental constitutive equation.
However, the current mathematical models that describe the mechanical properties of biological tissues do not account for the fundamental differences in tissue softening and hardening observed at different scales. Recent experimental research studies have revealed that the biomechanical features of biological tissues depend on the scale (i.e., size) at which the investigation is conducted [40]. Fuhs et al. studied the stiffness of breast cancer cells as well as the elasticity modulus of human breast tissues and found that breast cancer cells are softer than healthy cells while cancer tissue is more rigid compared to the healthy tissue [40]. This paradox cannot be explained in the context of classical continuum mechanics models as they are limited to large-scale dimensions and can only predict stiffness hardening at the tissue level [41].
To the best of our knowledge, no scale-dependent mathematical models have been proposed to solve this paradox (i.e., the rigid-tumours–soft-cells paradox). To address this, a scale-dependent modified continuum model incorporating stress nonlocality and strain gradient effects for the in-plane wave propagation analysis of biological tissues is developed. Particular attention is paid to human breast lesions to visualise the results of the scale-dependent model. The higher-order nonlocal strain gradient continuum mechanics theory, in conjunction with Hamilton’s law, is employed to derive the differential equations of the breast tissue. An analytical solution is obtained for the in-plane wave propagation in biological tissues. A deep learning model made of one flatten layer, several hidden layers, and a few dense layers is also developed and trained to improve the accuracy and allow for extra flexibility to better fit the clinical data. In addition, other machine learning models such as linear, Ridge, Lasso, Random Forest, and ElasticNet regression models are developed for comparison purposes. To verify the accuracy of the analytical approach, a finite element method is also implemented to extract the in-plane wave characteristics of healthy and diseased breast tissue including adipose and fibroglandular tissue, fibroadenoma, invasive ductal carcinoma, infiltrating lobular carcinoma, ductal carcinoma in situ, fibrocystic disease, fat necrosis, and invasive mucinous carcinoma. Figure 1 shows the flowchart of this research study with the four main steps, namely, (1) problem specification, (2) development of the higher-order nonlocal model, (3) machine learning modelling, and (4) finite element verification study. The new higher-order mathematical formulation integrated with machine learning allows for the introduction of three new parameters in biological tissue modelling: (1) zero-order nonlocal parameter, (2) first-order nonlocal parameter, and (3) strain gradient parameter. The combination of these parameters enables the description of both stiffness softening and hardening, which have been observed at cell and tissue levels, respectively, and solves the rigid-tumours–soft-cells paradox. The proposed mathematical model and simulation data have applications in biomedical elastography imaging for the early and accurate detection of breast cancer.

2. Theoretical Modelling

2.1. A Higher-Order Nonlocal Model of Biological Tissue

The main novel aspect of this research work is the development of the first scale-dependent, higher-order nonlocal model for analysing in-plane elastic waves propagated within human breast tissue. In addition, the higher-order model of nonlocality is integrated with a deep learning neural network model, which is later described in Section 2.3.
To the best of our knowledge, this is the first refined combination of a higher-order nonlocal strain gradient model and a deep learning model for the extraction of the wave propagation characteristics of biological samples at different scales.
In Eringen’s nonlocal elasticity theory [42], stress at a given point of the tissue is a function of the strain at all points of the domain (stress nonlocality). This is the basic assumption of the nonlocal continuum mechanics [42,43] and is initially stated in the form of an integro-partial differential equation between stress and strain. However, as this integro-partial differential equation is challenging and difficult to use for modelling purposes, Eringen introduced a simpler differential form and proved that it could incorporate stress nonlocality in the mechanical analysis [42,43]. Moreover, he found that experimental data and lattice dynamics simulations supported this differential form very nicely [42,43]. In this section, a more advanced version of the nonlocal continuum mechanics known as the higher-order nonlocal strain gradient theory [44] is used, which can incorporate not only the zeroth-order stress nonlocality but also the first-order one.
The classical elasticity theory assumes that the stress at a given point is only a function of the strain at that point. This does not allow for the incorporation of stress nonlocality, which has been widely reported to be associated with scale effects and is responsible for stiffness softening. Moreover, in the classical local models, the influence of the strain gradient is neglected, which is a complementary scale parameter associated with stiffness hardening at small scales. The combination of these effects within the framework of a newly introduced mathematical model solves the rigid-tumour–soft-cell paradox as it allows for both a reduction and increase in the stiffness due to the cancer’s formation and growth. Based on the higher-order version of the nonlocal strain gradient model (HNSM), the stress is related to strain as [44,45,46]
Γ 10 σ i j = Γ 1 C i j k l ε k l Γ 0 C i j k l Θ s t g 2 2 ε k l ,
where
Γ 10 f = f λ 0 n l s + λ 1 n l s 2 f + λ 0 n l s λ 1 n l s 4 f , Γ 1 f = f λ 1 n l s 2 f , Γ 0 f = f λ 0 n l s 2 f ,
and
2 f = 2 f x 2 + 2 f y 2 , 4 f = 2 2 f = 4 f x 4 + 2 4 f y 2 x 2 + 4 f y 4 ,
where σ i j , C i j k l , ε k l , and Θ s t g are the stress tensor, elastic property matrix, strain, and strain gradient parameter, respectively. Furthermore, the coefficients in Equation (2) are given by λ 1 n l s = κ 1 l c e l l 2 and λ 0 n l s = κ 0 l c e l l 2 , which are, respectively, the first and second nonlocal stress parameters [46,47,48,49]. Γ 10 f , Γ 1 f , and Γ 0 f are the combined, first-order, and zeroth-order nonlocal operators, respectively. When the zeroth-order nonlocal operator is equal to the first-order one (i.e., Γ 0 = Γ 1 ), the HNSM is reduced to the lower-order nonlocal strain gradient model (LNSM) as in [50,51].
Γ 0 σ i j = C i j k l ε k l C i j k l Θ s t g 2 2 ε k l .
On the other hand, the classical nonlocal model (CLNM) is recovered when the nonlocal operators are identical and the strain gradient operator is zero, namely, Γ 0 = Γ 1   and   Θ s t g = 0 :
Γ 0 σ i j = C i j k l ε k l .
The CLNM has been successfully utilised to analyse the biomechanics of protein microtubules [52,53], providing evidence on the scale-dependency behaviour of biological components and the importance of nonlocal effects [54,55]. If only the strain gradient effects are included, both nonlocal operators are neglected, and the strain gradient operator is maintained ( Γ 0 = 0   and   Θ s t g 0 ), leading to the following constitutive equation of biological tissue.
σ i j = C i j k l ε k l C i j k l Θ s t g 2 2 ε k l .
Equation (6) describes the fundamental stress–strain relation of the classical strain gradient model (CLSM) [56,57,58]. Moreover, the constitutive equation of the HNSM can be reduced to the classical elasticity model (CLEM) when Γ 0 = Γ 1 = 0   and   Θ s t g = 0 , which is mathematically described as:
σ i j = C i j k l ε k l .
In this section, a higher-order model of elasticity with strain gradient and stress nonlocality effects is developed for the in-plane wave propagations in rectangular-shaped, thin biological tissue sections (Figure 2). Using Equations (1) and (2) in the two-dimensional Cartesian coordinate system, the stress components are linked to the strain components as
σ x x λ 0 n l s + λ 1 n l s 2 σ x x + λ 0 n l s λ 1 n l s 4 σ x x = E 11 1 ν 12 ν 21 Γ 1 ε x x Θ s t g 2 Γ 0 2 ε x x + ν 12 E 22 1 ν 12 ν 21 Γ 1 ε y y Θ s t g 2 Γ 0 2 ε y y , σ y y λ 0 n l s + λ 1 n l s 2 σ y y + λ 0 n l s λ 1 n l s 4 σ y y = ν 12 E 22 1 ν 12 ν 21 Γ 1 ε x x Θ s t g 2 Γ 0 2 ε x x + E 22 1 ν 12 ν 21 Γ 1 ε y y Θ s t g 2 Γ 0 2 ε y y , σ x y λ 0 n l s + λ 1 n l s 2 σ x y + λ 0 n l s λ 1 n l s 4 σ x y = 2 G 12 Γ 1 ε x y 2 G 12 Θ s t g 2 Γ 0 2 ε x y ,
ν i j and E i j are the Poisson’s ratio and elasticity modulus of the biological tissue.
For small deformation induced by the in-plane wave propagation, the strain components are given by
ε x x ( x , y , z ) = x u ( x , y ) , ε y y ( x , y , z ) = y v ( x , y ) , ε x y ( x , y , z ) = 1 2 y u ( x , y ) + x v ( x , y ) ,
where u and v are the in-plane displacements along the x and y axes, respectively. The force stress resultants can be obtained from Equations (8) and (9) as
T x x λ 0 n l s + λ 1 n l s 2 T x x + λ 0 n l s λ 1 n l s 4 T x x = A 11 Γ 1 u x Θ s t g 2 Γ 0 2 u x + A 12 Γ 1 v y Θ s t g 2 Γ 0 2 v y , T y y λ 0 n l s + λ 1 n l s 2 T y y + λ 0 n l s λ 1 n l s 4 T y y = A 21 Γ 1 u x Θ s t g 2 Γ 0 2 u x + A 22 Γ 1 v y Θ s t g 2 Γ 0 2 v y , T x y λ 0 n l s + λ 1 n l s 2 T x y + λ 0 n l s λ 1 n l s 4 T x y = A 33 Γ 1 u y + v x Θ s t g 2 Γ 0 2 u y + v x ,
where the force stress resultants (Tij) are defined by
T x x = h / 2 h / 2 σ x x d z ,   T x y = h / 2 h / 2 σ x y d z ,   T y y = h / 2 h / 2 σ y y d z .
Here, h denotes the thickness of the tissue section. Aij are the in-plane stiffness components of the tissue, which are calculated using
A 11 = E 11 h 1 ν 21 ν 12 ,   A 22 = E 22 h 1 ν 21 ν 12 ,   A 12 = A 21 = ν 12 E 22 h 1 ν 21 ν 12 ,   A 33 = G 12 h ,
where G12 is the shear elasticity modulus of the biological tissue. Using Hamilton’s law for motions in the x-y plane, one can obtain
T x x x + T x y y = m 0 2 u t 2 , T x y x + T y y y = m 0 2 v t 2 ,
where
m 0 = h / 2 h / 2 ρ b i o d z .  
Here, ρ b i o is the biological tissue mass density. It is worth mentioning that mass density is entirely different from mammographic breast density. In the context of breast tissue, mammographic density refers to the relative amount of fibroglandular tissue to fat tissue, which is an independent risk factor for breast cancer and is also a significant problem in mammography imaging as it masks cancer [59,60,61]. By contrast, mass density in Equation (14) refers to the mass per unit volume of the tissue, which is measured in terms of kg/m3. Using Equations (10) and (13), the governing differential equations for predicting the deformation behaviour of biological tissues under in-plane wave propagation are derived, which are given in Appendix A as Equations (A1) and (A2). It should be noted that, in the present analysis, the effect of energy loss is not considered. In terms of the mechanical response of breast tissues, it has been shown that, after an initial transient period, the system reaches a steady state, and the assumption of negligible energy loss becomes valid [35].
To determine the in-plane wave propagation characteristics of the breast tissue, the solution of Equations (A1) and (A2) can be assumed as
u = U exp k y y + k x x ω t i , v = V exp k y y + k x x ω t i ,
where V and U represent the in-plane vibration amplitudes along the y and x axes, respectively. The in-plane wave numbers along the y and x axes are ky and kx, respectively. Moreover, ω is the frequency of vibrations induced by in-plane waves. Substituting Equation (15) into the in-plane motion equations, namely, Equations (A1) and (A2), leads to
Θ i , j U V ω i p 2 Π i , j U V = 0 0 ,
where Θ i , j , Π i , j , and ω i p represent the total stiffness matrix, total mass matrix, and in-plane vibration frequency, respectively. The square root of the eigenvalues of Equation (16) are the frequencies of the in-plane vibrations within the tissue induced by in-plane waves. The elements of matrix Θ i , j and Π i , j are listed in Appendix A.

2.2. Finite Element Simulation

To verify the accuracy of the higher-order nonlocal modelling, finite element simulations have been conducted using COMSOL Multiphysics Simulation software version 5.5. The Solid Mechanics (solid) module in conjunction with a linear elastic material model was used for the analysis. The general structural analysis of three-dimensional bodies was adopted for a more accurate analysis. The eigenfrequency solver technique was the multifrontal massively parallel sparse direct (MUMPS), which is a powerful rapid numerical method for solving large linear systems using a parallel computational approach. Three different breast tissue samples were considered for the finite element study. First, fibroglandular tissue was tested as the main source of mammographic breast density, which causes a masking effect problem in mammography imaging [62,63]. The second case study has been performed on ductal carcinoma in situ (DCIS) samples, which are widely known as the early form of human breast cancer [64]. DCIS cancers are noninvasive, meaning that they have not invaded surrounding breast tissue or spread into other parts of the body. Detecting DCIS enables more effective and faster treatment, improving the survival rates of breast cancer patients. The third samples, which have been considered for the finite element study, were fat necrosis as a benign condition that might be misdiagnosed as cancerous tumours. They are commonly formed after a trauma or injury [65]. The finite element simulations were conducted for the first ten mode shapes of tissue samples. A rectangular-shaped geometry with a thickness of 1 cm, length of 10 cm, and width of 10 cm was created in a three-dimensional coordinate system as the tissue section. Scale effects including the stress nonlocality and strain gradients were neglected as the available commercial finite element software packages lack the ability to simulate these effects. After the geometry creation, the elastic properties such as Poisson’s ratio and Young’s modulus were assigned to the material (tissue sections) based on the available clinical data in the literature [28]. The displacement and velocity components were initialised to zero at the start of the analysis, indicating a zero or at-rest initial condition. The boundary conditions of the samples were assumed to be the following: the displacement components along the x and z axes were zero at x = 0 and x = a, the displacement components in the y and z directions were set to zero at y = 0 and y = b, and a free type of boundary condition was adopted for the top and bottom surfaces of the rectangular sample. A combination of tetrahedra and triangles (element type) was used to create a smooth and well-distributed mesh pattern with an average extra-fine mesh size. An eigenfrequency technique was utilised as the solver, which is suitable for eigenvalue problems with eigenmodes.

2.3. Artificial Neural Network Simulation

To capture the different mechanical behaviours of breast lesions at different scales and include both stiffness hardening and softening, a machine learning model that can account for nonlocality, complexity, and nonlinear patterns is required. In addition, the model needs to have the capability of dealing with a limited number of samples and small datasets while preventing overfitting. Deep neural network models, especially with rectified linear unit (ReLU) activation functions, provide the ability to incorporate the complexity of the in-plane wave characteristics of breast lesions and their nonlocal deformation behaviour at small scales. Furthermore, neural network models allow for flexibility in the architecture by offering different hyperparameter adjustment options including the change in the number of hidden layers and neurons as well as the choice of activation function, optimiser, and learning rate. This flexibility in the fine-tuning process is essential for optimising the performance of the machine learning model and for preventing overfitting when the training dataset is relatively small, as it is in many clinical diagnosis applications. Overall, the benefits of neural network models include—but are not limited to—a fast training process, flexibility in fine-tuning, the ability to predict complexity, and the capability of mitigating overfitting.
TensorFlow, NumPy, and pandas were used to develop the deep neural network model, to work with multi-dimensional matrixes, and to create the database, respectively [66,67]. Keras, as a high-level application programming interface for conveniently working with TensorFlow 2.14.0 machine learning software, was adopted [68]. Moreover, to divide the dataset into 70% training and 30% testing subsets, the train_test_split method of the model_selection module of the sklearn python library was used. The dataset used for training and testing contained 100 human samples diagnosed with DCIS. To build the deep learning model, several architectures were considered using the sequential method of Keras. Overall, the neural network was composed of one flatten layer, a couple of hidden dense layers, and an output layer as shown in Figure 3a. In clinical applications, it is often challenging to find patients of all different pathological conditions, leading to relatively small datasets. Thus, a simple neural network can be adopted to prevent overfitting and ensure the generalisation capability of the model despite the restricted human breast samples. In addition, simpler neural network algorithms are computationally inexpensive and faster to implement as they require less training time as well as a faster hyperparameter adjustment process. A sequential deep neural network model of feedforward type was utilised in this study for prediction, which was made of several hidden layers. A rectified linear unit (ReLU) was employed as the activation function of dense layers, which allows for the incorporation of nonlinearity and complexity in real-world problems. An Adam optimiser with a learning rate of 0.001 and mean-square-error loss function was used to compile the machine learning model. Data scaling was performed prior to the training process to ensure consistency with the scale of all features. The training process was conducted for 500 epochs. During each epoch, the entire training dataset was passed forward and backwards within the artificial neural network. The rationale for choosing a relatively high epoch number is the limited sample size, which is usually the case in clinical studies where providing human tissue samples is not straightforward. This ensures that the neural network model has been sufficiently exposed to the available data for training and generalisation purposes. Figure 3b shows a flowchart to visualise the different steps and processes involved in the machine learning modelling of the in-plane wave propagation analysis of human breast lesions. The model evaluation, which is essential for hyperparameter adjustment, was conducted using the mean square error of both training and test datasets. The learning rate, number of neurons within a dense layer, number of dense layers, activation function, number of epochs, and the optimiser technique have been adjusted during the hyperparameter adjustment process to obtain a highly accurate neural network model.
In this research work, three different architectures of neural networks were considered: (1) FL-3HDL(512,128,10): one flatten layer (FL), three hidden layers (HDL) with 512, 128, and 10 neurons as well as one dense layer with one neuron as output layer (5 layers and 651 neurons in total); (2) FL-2HDL(128,10): one FL, two HDLs with 512 and 10 neurons, as well as one dense layer with one neuron; (3) FL-1HDL(10): one FL, one HDL with 10 neurons, and one output dense layer. In all three different architectures, the activation function of the hidden and dense layers was set to the rectified linear unit (ReLU) to incorporate nonlocality and nonlinearity and to capture complex patterns within the data. The Adam optimiser with a learning rate of 0.001 was used to compile the machine learning model. The results of these three neural network models will be compared and discussed in the next section.

3. Results

To verify the analytical higher-order nonlocal strain gradient modelling, a finite element study was conducted on three different types of breast tissue, including healthy fibroglandular tissue (FGT), ductal carcinoma in situ (DCIS), and fat necrosis, which is a benign condition. The results (dimensionless phase velocities) are listed in Table 1 for the first three mode numbers (n = 1, 2, 3). The elasticity moduli of different human breast tissue were adopted from those reported by Samani et al. [28]. They conducted experimental mechanical testing on 169 human samples to extract the elasticity properties of pathological breast tissues. The scale parameters were set to zero as these parameters are not incorporated in commercial finite element software packages like Comsol Multiphysics version 5.5. The thickness of the tissue samples was set to 1 cm while the tissue side length was 10 cm in this study. The Poisson’s ratio of all breast lesions was taken as 0.495 [28] while the physical mass density of breast tissue was set to 945 kg/m3 [69,70]. The patterns of in-plane wave propagation within the breast tissue, obtained using the finite element simulations, are also depicted in Figure 4 for the first four in-plane modes. The visualisation was performed on the DCIS. Similar patterns were found for other breast lesion types, which are not included here for the sake of brevity.
Table 2 lists the mean square errors of training and testing datasets for the neural network of deep learning models studied in this paper. The outcome is indicated for three different deep learning architectures including (1) FL-3HDL(512,128,10), (2) FL-2HDL(128,10), and (3) FL-1HDL(10). Here, FL and HDL stand for the flatten layer and hidden dense layer, respectively. The integer preceding the ‘HDL’ indicates the quantity of hidden dense layers in the architecture of the neural network, and the corresponding numbers of neurons in each layer are denoted in parentheses in an ordered manner. For example, FL-2HDL(128,10) refers to a neural network with one flatten layer and two hidden dense layers, including 10 and 128 units within the second and first hidden layers, respectively. The mean square errors of these models are listed for various numbers of epochs from 1 to 500. The first model with one FL and three HDLs demonstrated the minimum mean square error, reaching around 4.71 × 10−7 and 1.06 × 10−6 after 500 epochs for the train and test datasets, respectively. In all cases, the mean square error of the deep learning model significantly decreased as the epoch number increased, indicating the importance of an appropriate epoch number in the accuracy of the machine learning modelling. Table 3 compares the predicted dimensionless wave velocity within the early breast cancer of DCIS obtained using the deep learning model with those of the analytical approach (actual data) for different dimensionless (normalised) wave numbers. Scaling (normalisation) was conducted by dividing each feature by its maximum value to ensure that all features were in the same range.
The results of a number of different machine learning models are listed and compared in Table 4. The machine learning models were used to predict the dimensionless wave velocity of in-plane elastic waves propagated within the breast lesions of women with high mammographic density. All models have been trained on 100 data points obtained using the higher-order nonlocal model. In addition to the neural network algorithm, Random Forest, Lasso, linear, Ridge, and ElasticNet regression models are presented in this table for comparison purposes. In this study, the neural network model consisted of one flatten layer and three hidden layers with 512, 128, and 10 neurons as well as one dense layer at the end with one neuron. The rectified linear unit (ReLU) function was used as the activation function to incorporate complexity in the form of nonlinearity and nonlocality into the model. The RandomForestRegressor from the open-source scikit-learn library was utilised with 100 estimators (decision trees). Each estimator made a prediction, and the final output was the average of these predictions for the regression problems. To adjust the machine learning models of regression and obtain the best hyperparameters, the GridSearchCV method of the sklearn.model_selection was used. The optimised values of the regularization parameter alpha for ElasticNet, Lasso, and Ridge were obtained as 0.005, 0.001, and 0.2, respectively. From this table, it is found that the neural network and Random Forest led to more accurate predictions compared to the other machine learning models. The results calculated using the neural network were slightly better than those of the Random Forest model (closer to the actual data) for various amounts of in-plane dimensionless wave numbers.
Figure 5 and Figure 6 show the variation of phase velocity and frequency with respect to the in-plane wave number for the first and second modes, respectively. Healthy fibroglandular tissue, which is the main component of breast tissue, with high mammographic density and DCIS are compared in these figures. The first and second nonlocal parameters were set to 2 and 1 cm, respectively, and the strain gradient parameter was 0.5 cm. It was assumed that the thickness of the tissue was 1 cm. The phase velocity of the dense breast tissue was significantly lower than that of DCIS in both the first and second in-plane wave modes. In addition, the frequency of in-plane wave propagation within the DCIS was higher than that of the dense breast tissue due to the fact that the DCIS exhibited more rigidity under wave propagation and possessed a higher stiffness. In Figure 5c,d, the influence of the length of the breast sample on the in-plane phase velocity curves is shown for human dense breast and DCIS, respectively. Various non-classical continuum models, together with the traditional elasticity model, were developed and used to plot these figures. In the higher-order nonlocal strain gradient model, the first and second nonlocal parameters were, respectively, set to 2 and 1 microns, and the strain gradient parameter was 0.5 microns. In the nonlocal model of the breast samples, the strain gradient effect was ignored while the first and second nonlocal parameters were 2 and 1 micron, respectively. By contrast, both the first and second nonlocal parameters were set to zero while the strain gradient parameter was taken as 0.5 micron in the strain gradient model. In the traditional classical model, the effects of all scale parameters were neglected, leading to a scale-free continuum model. The phase velocity estimated using the nonlocal model was lower than that of the classical model of breast tissues due to the stiffness-softening effect of stress nonlocality. By contrast, the strain gradient effect was associated with stiffness hardening and, hence, higher phase velocities. The nonlocal model incorporated stiffness softening at small scales while the strain gradient model could estimate structural stiffness hardening of the breast tissue. A refined combination of the models developed in this study simultaneously simulated both softening and hardening responses. At larger scales, as anticipated, all models exhibited convergence toward the classical model.
To study the influence of various scale parameters including first nonlocal, second nonlocal, and strain gradient on the in-plane wave propagation characteristics of breast lesions, phase velocity versus the in-plane wave number was plotted in Figure 7a–c, respectively. The breast lesion was assumed to be of intermediate invasive ductal carcinoma (IDC) type, which accounts for around 75–80% of all invasive breast cancers [71]. The elasticity and Poisson’s ratio of the intermediate IDC were taken as 19.99 kPa and 0.495, respectively [28]. From Figure 7a, it can be concluded that, when the first nonlocal parameter was increased from 0 to 2 cm, the phase velocity experienced a significant reduction. This change in phase velocity became higher by increasing the in-plane wave number. A slight reduction in the phase velocity of in-plane waves within the intermediate IDC by increasing the second nonlocal parameter was found (Figure 7b). At small in-plane wave numbers (k < 30 1/m), the influence of the first-order stress nonlocality was negligible and the zeroth-order nonlocality was dominant. Figure 7c depicts the variation of phase velocity with the wave number for different strain gradient parameters and for in-plane wave propagation. The phase velocity was very sensitive to the amount of strain gradient, especially at higher in-plane wave numbers. A small increase in the strain gradient parameter led to a sharp increase in phase velocity. Strain gradient effects were associated with stiffness hardening, leading to an increase in the frequency and phase velocity of in-plane waves within the human breast lesion.
To investigate and compare the wave propagation characteristics of various breast lesions, the frequency and phase velocity were plotted in Figure 8 and Figure 9 for nonlocal and strain gradient models, respectively. Poisson’s ratio of human breast tissue was considered to be 0.495. Elasticity moduli of adipose, fibroglandular, fibroadenoma, DCIS, low-grade, intermediate-grade, and high-grade IDC were 3.25 ± 0.91, 3.24 ± 0.61, 6.41 ± 2.86, 16.38 ± 1.55, 10.40 ± 2.60, 19.99 ± 4.2, and 42.52 ± 12.47, respectively [28]. Furthermore, the elasticity moduli of the infiltrating lobular carcinoma (ILC), fibrocystic disease, invasive mucinous carcinoma (IMC), and fat necrosis were, respectively, 15.62 ± 2.64, 17.11 ± 7.35, 20.21 ± 4.2, and 4.45 ± 0.61. The in-plane wave number was set to 20 for this investigation. Based on the frequency and phase velocity, the differentiation between the common benign conditions such as the fibroadenoma and fat necrosis and cancerous breast lesions including IMC and IDC would be feasible in elastography biomedical imaging using in-plane elastic wave propagation. A remarkable difference was also observed between the in-plane wave response of healthy breast tissue and that of early breast cancer in the form of DCIS.

4. Discussion

From Table 1, it can be concluded that there is an excellent match between the results of finite element simulation and those of the analytical model in all cases. The analytical model was developed based on the higher-order nonlocal theory, in which the gradients of strain were also incorporated. However, as the scale of the biological component increases, strain gradient and nonlocal effects become less important and can be ignored at a certain scale. In this case, the higher-order model reduces to the classical model, which can be validated using finite element simulations. The validation study was performed on three important breast tissues, including DCIS, healthy fibroglandular tissue, and fat necrosis. The DCIS and fibroglandular tissue are, respectively, taken into consideration as early breast cancers and a type of healthy tissue. In the mammography imaging of mammographically dense breasts, in which there is a significant amount of fibroglandular tissue compared to fat tissue, the DCIS is likely to be hidden by the fibroglandular tissue, making the diagnosis process challenging. In addition, fat necrosis was considered in this finite element study as a common type of benign condition. Figure 4 shows that, at higher in-plane modes, the change in the tissue displacement occurs in smaller regions, making the deformation field more complex but more section-sensitive. Exciting higher vibration modes could be an effective way to improve the resolution and sensitivity of wave propagation-based detection tools like elastography imaging, as higher modes are associated with more ups and downs in the tissue, and, consequently, finer spatial details about the tissue’s biomechanical features.
Table 2 indicates the importance of an appropriate architecture design for the deep learning model. When the number of hidden layers is increased from one to three, the MSE of both test and train data significantly decreases. This significant reduction in the MSE of the model is not only related to the number of hidden layers but also is associated with the number of neurons in the architecture. The third model of deep learning, which is called FL-1HDL(10) in Table 2, contains 10 neurons, while the first model with the name FL-3HDL(512,128,10) incorporates 650 neurons. Table 3 further confirms this finding by comparing the results of the deep learning model with those of actual data obtained by using the higher-order nonlocal model of elasticity. It is observed that a deep learning algorithm with one flatten layer and two dense layers containing around several hundred neurons can highly accurately predict the results of the nonlocal elasticity model without having any elasticity theory knowledge.
The rigid-tumour–soft-cell paradox is associated with stiffness hardening at the tissue scale level while the stiffness of a clump of cells experiences a softening behaviour. The classical model of elasticity can only predict stiffness hardening and lacks the simulation of stiffness softening at cell scales. The nonlocal elasticity model with higher-order scale and strain parameters can predict both stiffness hardening and softening depending on the scale of the biological component or tissue. At large scales, the strain gradient effect is dominant, which is highly linked with stiffness hardening, as can be concluded from Figure 5, Figure 6 and Figure 7. However, at certain cell levels, the stress nonlocalities of zeroth and first orders account for the stiffness softening, as observed in Figure 7a,b. Specifically, in terms of wave propagation characteristics, stress nonlocalities lead to a reduction in stiffness, and this results in lower frequencies and phase velocities, which is supported by the data in Figure 8 and Figure 9. By contrast, strain gradients are linked to a significant enhancement in the total stiffness, leading to higher in-plane wave frequencies and phase velocities. Another important finding is that, when the wave mode number increases, the difference between the malignant and healthy breast tissue significantly increases at both tissue and cellular levels, making the frequency-based diagnosis process easier and more accurate. Moreover, stress nonlocalities become more prominent at higher in-plane modes since the interaction between different sections of the breast tissue is greater when the in-plane wave number increases. Based on the comparison study made in this research, a significant difference between the in-plane wave propagation characteristics of healthy fibroglandular/adipose tissues and that of early breast cancer (DCIS) was observed, which makes wave-propagation-based imaging tools like elastography promising and feasible for the early detection of breast cancer. Unlike X-ray mammography, there is no masking effect when the density (i.e., the relative amount of fibroglandular tissue) of the human breast increases, making in-plane wave propagation-based elastography imaging a promising complementary imaging technology for women with dense breasts.

5. Conclusions

In this study, a higher-order model of nonlocal elasticity in conjunction with a deep learning algorithm was developed and used to accurately obtain the in-plane wave propagation characteristics of human breast tissues. Eleven different healthy, benign, and pathological conditions including adipose, fibroglandular, DCIS, ILC, IDC (intermediate, low- and high-grade), IMC, fat necrosis, and fibrocystic and fibroadenoma were taken into account. Using Hamilton’s law, the differential equations of the in-plane motion of breast tissue were derived. Analytical solutions were presented for the frequency and phase velocity of the elastic waves propagated within the thin breast tissue sections. To show the validity of the model, three finite element studies were conducted on DCIS as the very early form of breast cancer, using fibroglandular tissue as the healthy mammographically dense breast and fat necrosis as a common benign condition. The results of the finite element studies were in excellent agreement with those of higher-order nonlocal elasticity modelling. Furthermore, an artificial deep neural network was designed, created, and adjusted for the in-plane wave propagation analysis of breast lesions. The deep learning model can discover the mathematical formulas that govern the in-plane wave propagation response of breast lesions and precisely predict the phase velocity with training on only a few hundred samples. For the deep learning model with FL-3HDL(512,128,10) architecture, the mean square error of training and testing datasets were, respectively, reduced to 4.7 × 10−7 and 1.06 × 10−6 after 500 epochs. The results of the deep learning model were also compared to those of other machine learning models, including linear, Ridge, Lasso, Random Forest, and ElasticNet regression models. An excellent match was found between the predictions made by the neural network algorithm and that of the Random Forest model. The phase velocity of in-plane waves within the DCIS was, on average, more than 2.5 times higher than that of healthy fibroglandular tissue. Moreover, a significant difference (p-value < 0.01) was found between breast cancers such as high-grade IDC and IMC and benign conditions like fat necrosis and fibroadenoma. Both zeroth- and first-order stress nonlocalities were associated with a reduction in the phase velocity of in-plane waves since they were linked to tissue stiffness softening. By contrast, the strain gradient parameter enhanced the in-plane phase velocity of breast tissue as it was correlated with stiffness hardening. The combination of nonlocal and strain gradient parameters in the higher-order elasticity model allowed for the concurrent incorporation of stiffness hardening and softening, solving the rigid-tumour–soft-cell paradox of breast biomechanics.

6. Patents

Ali Farajpour, Wendy V. Ingman, “Scale-dependent elastography method for detection of small inclusions in biological tissue”. Applicant: The University of Adelaide, Application number: PCT/AU2023/050855. Patent Cooperation Treaty (PCT) Submission Date: 1 September 2023.

Author Contributions

Conceptualization, A.F. and W.V.I.; methodology, A.F. and W.V.I.; software, A.F.; validation, A.F. and W.V.I.; formal analysis, A.F.; investigation, A.F.; resources, A.F. and W.V.I.; data curation, A.F.; writing—original draft preparation, A.F.; writing—review and editing, W.V.I.; visualization, A.F.; supervision, W.V.I.; project administration, W.V.I.; funding acquisition, A.F. and W.V.I. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Robinson Research Institute’s Innovation Seed Funding and the Faculty of Health and Medical Sciences (Adelaide Medical School) Building Research Leaders Award given to A.F.

Institutional Review Board Statement

Ethical review and approval were waived by the University of Adelaide Human Research Ethics Committee as the research involved the use of existing collections of data that contain only non-identifiable data about human beings.

Informed Consent Statement

Patient consent was waived because the research involved the use of existing collections of data that contain only non-identifiable data about human beings. Permission to reuse data (elastic moduli of human breast tissues) from Table 1 of Ref. [28] has been formally granted through the Copyright Clearance Center (CCC) from IOP Publishing, Ltd.

Data Availability Statement

Data are available upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The governing equations of the breast tissue under the assumption of in-plane wave propagation using the higher-order nonlocal strain gradient model are derived as
A 11 2 u x 2 A 11 λ 1 n l s + Θ s t g 2 4 u x 4 + 4 u y 2 x 2 + A 11 Θ s t g 2 λ 0 n l s 6 u x 6 + 2 6 u y 2 x 4 + 6 u y 4 x 2 + A 12 2 v x y A 12 λ 1 n l s + Θ s t g 2 4 v x 3 y + 4 v x y 3 + A 12 Θ s t g 2 λ 0 n l s 6 v x 5 y + 2 6 v y 3 x 3 + 6 v x y 5 + A 33 2 u y 2 + 2 v x y A 33 λ 1 n l s + Θ s t g 2 4 u x 2 y 2 + 4 u y 4 + 4 v y x 3 + 4 v y 3 x + A 33 Θ s t g 2 λ 0 n l s 6 u x 4 y 2 + 2 6 u x 2 y 4 + 6 u y 6 + 6 v y x 5 + 2 6 v y 3 x 3 + 6 v y 5 x = m 0 2 u t 2 m 0 λ 0 n l s + λ 1 n l s 4 u x 2 t 2 + 4 u y 2 t 2 + m 0 λ 0 n l s λ 1 n l s 6 u x 4 t 2 + 2 6 u y 2 x 2 t 2 + 6 u y 4 t 2 ,
and
A 33 2 u x y + 2 v x 2 A 33 λ 1 n l s + Θ s t g 2 4 u x 3 y + 4 u x y 3 + 4 v x 4 + 4 v y 2 x 2 + A 33 Θ s t g 2 λ 0 n l s 6 u x 5 y + 2 6 u y 3 x 3 + 6 u x y 5 + 6 v x 6 + 2 6 v y 2 x 4 + 6 v y 4 x 2 + A 21 2 u y x A 21 λ 1 n l s + Θ s t g 2 4 u y x 3 + 4 u y 3 x + A 21 Θ s t g 2 λ 0 n l s 6 u y x 5 + 2 6 u y 3 x 3 + 6 u y 5 x + A 22 2 v y 2 A 22 λ 1 n l s + Θ s t g 2 4 v x 2 y 2 + 4 v y 4 + A 22 Θ s t g 2 λ 0 n l s 6 v x 4 y 2 + 2 6 v y 4 x 2 + 6 v y 6 = m 0 2 v t 2 m 0 λ 0 n l s + λ 1 n l s 4 v x 2 t 2 + 4 v y 2 t 2 + m 0 λ 0 n l s λ 1 n l s 6 v x 4 t 2 + 2 6 v y 2 x 2 t 2 + 6 v y 4 t 2 .
The elements of the matrix components of Equation (16), namely, Θ i , j and Π i , j , are calculated using
Θ 1 , 1 = A 11 i k x 2 A 11 λ 1 n l s + Θ s t g 2 i k y 2 i k x 2 + i k x 4 + A 11 Θ s t g 2 λ 0 n l s i k x 6 + 2 i k x 4 i k y 2 + i k y 4 i k x 2 + A 33 i k y i k y A 33 λ 1 n l s + Θ s t g 2 i k y 2 i k x 2 + i k y 2 + A 33 Θ s t g 2 λ 0 n l s i k y 2 i k x 4 + 2 i k y 4 i k x 2 + i k y 6 , Θ 1 , 2 = A 12 i k y i k x 1 λ 1 n l s + Θ s t g 2 i k x 2 + i k y 2 + A 12 Θ s t g 2 λ 0 n l s i k y i k x i k x 4 + 2 i k y 2 i k x 2 + i k y 4 + A 33 i k y i k x A 33 λ 1 n l s + Θ s t g 2 i k y i k x i k x 2 + i k y 2 + A 33 Θ s t g 2 λ 0 n l s i k y i k x 5 + 2 i k y 3 i k x 3 + i k y 5 i k x , Θ 2 , 1 = A 33 i k y i k x A 33 λ 1 n l s + Θ s t g 2 i k y i k x 3 + i k y 3 i k x + A 33 Θ s t g 2 λ 0 n l s i k y i k x 5 + 2 i k y 3 i k x 3 + i k y 5 i k x + A 21 i k y i k x A 21 λ 1 n l s + Θ s t g 2 i k y i k x 3 + i k y 3 i k x + A 21 Θ s t g 2 λ 0 n l s i k y i k x 5 + 2 i k y 3 i k x 3 + i k y 5 i k x , Θ 2 , 2 = A 33 i k x 2 A 33 λ 1 n l s + Θ s t g 2 i k x 4 + i k y 2 i k x 2 + A 33 Θ s t g 2 λ 0 n l s i k x 6 + 2 i k y 2 i k x 4 + i k y 4 i k x 2 + A 22 i k y 2 A 22 λ 1 n l s + Θ s t g 2 i k y 2 i k x 2 + i k y 4 + A 22 Θ s t g 2 λ 0 n l s i k y 2 i k x 4 + 2 i k y 4 i k x 2 + i k y 6 .
and
Π 1 , 1 = m 0 m 0 λ 0 n l s + λ 1 n l s i k x 2 + i k y 2 + m 0 λ 0 n l s λ 1 n l s i k x 4 + 2 i k y 2 i k x 2 + i k y 4 , Π 1 , 2 = 0 ,   Π 2 , 1 = 0 , Π 2 , 2 = m 0 m 0 λ 0 n l s + λ 1 n l s i k x 2 + i k y 2 + m 0 λ 0 n l s λ 1 n l s i k x 4 + 2 i k y 2 i k x 2 + i k y 4 .

References

  1. Joldes, G.R.; Doyle, B.; Wittek, A.; Nielsen, P.M.; Miller, K. Computational Biomechanics for Medicine: Imaging, Modeling and Computing; Springer: Berlin/Heidelberg, Germany, 2016; pp. 1–203. [Google Scholar]
  2. Stashko, C.; Hayward, M.-K.; Northey, J.J.; Pearson, N.; Ironside, A.J.; Lakins, J.N.; Oria, R.; Goyette, M.-A.; Mayo, L.; Russnes, H.G. A convolutional neural network STIFMap reveals associations between stromal stiffness and EMT in breast cancer. Nat. Commun. 2023, 14, 3561. [Google Scholar] [CrossRef]
  3. Knudson, D.V.; Knudson, D. Fundamentals of Biomechanics; Springer: Berlin/Heidelberg, Germany, 2007; Volume 183. [Google Scholar]
  4. Gao, X.; Fraulob, M.; Haïat, G. Biomechanical behaviours of the bone–implant interface: A review. J. R. Soc. Interface 2019, 16, 20190259. [Google Scholar] [CrossRef] [PubMed]
  5. Ajalloueian, F.; Lemon, G.; Hilborn, J.; Chronakis, I.S.; Fossum, M. Bladder biomechanics and the use of scaffolds for regenerative medicine in the urinary bladder. Nat. Rev. Urol. 2018, 15, 155–174. [Google Scholar] [CrossRef] [PubMed]
  6. Noroozi, R.; Arif, Z.U.; Taghvaei, H.; Khalid, M.Y.; Sahbafar, H.; Hadi, A.; Sadeghianmaryan, A.; Chen, X. 3D and 4D Bioprinting Technologies: A Game Changer for the Biomedical Sector? Ann. Biomed. Eng. 2023, 51, 1683–1712. [Google Scholar] [CrossRef] [PubMed]
  7. Yang, L.; Li, X.; Wang, D.; Mu, S.; Lv, W.; Hao, Y.; Lu, X.; Zhang, G.; Nan, W.; Chen, H. Improved mechanical properties by modifying fibrin scaffold with PCL and its biocompatibility evaluation. J. Biomater. Sci. Polym. Ed. 2020, 31, 658–678. [Google Scholar] [CrossRef] [PubMed]
  8. Pourmehran, O.; Cazzolato, B.; Tian, Z.; Arjomandi, M. The effect of inlet flow profile and nozzle diameter on drug delivery to the maxillary sinus. Biomech. Model. Mechanobiol. 2022, 21, 849–870. [Google Scholar] [CrossRef]
  9. Ferreira, J.A.; Jordão, D.; Pinto, L. Drug delivery enhanced by ultrasound: Mathematical modeling and simulation. Comput. Math. Appl. 2022, 107, 57–69. [Google Scholar] [CrossRef]
  10. Lan, G.; Larin, K.V.; Aglyamov, S.; Twa, M.D. Characterization of natural frequencies from nanoscale tissue oscillations using dynamic optical coherence elastography. Biomed. Opt. Express 2020, 11, 3301–3318. [Google Scholar] [CrossRef]
  11. Lan, G.; Aglyamov, S.; Larin, K.V.; Twa, M.D. In vivo human corneal natural frequency quantification using dynamic optical coherence elastography: Repeatability and reproducibility. J. Biomech. 2021, 121, 110427. [Google Scholar] [CrossRef]
  12. Kennedy, B.F.; Wijesinghe, P.; Sampson, D.D. The emergence of optical elastography in biomedicine. Nat. Photonics 2017, 11, 215–221. [Google Scholar] [CrossRef]
  13. Cui, X.-W.; Li, K.-N.; Yi, A.-J.; Wang, B.; Wei, Q.; Wu, G.-G.; Dietrich, C.F. Ultrasound elastography. Endosc. Ultrasound 2022, 11, 252. [Google Scholar]
  14. Korukonda, S.; Nayak, R.; Carson, N.; Schifitto, G.; Dogra, V.; Doyley, M.M. Noninvasive vascular elastography using plane-wave and sparse-array imaging. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2013, 60, 332–342. [Google Scholar] [CrossRef]
  15. Gennisson, J.-L.; Deffieux, T.; Fink, M.; Tanter, M. Ultrasound elastography: Principles and techniques. Diagn. Interv. Imaging 2013, 94, 487–495. [Google Scholar] [CrossRef]
  16. Wong, V.W.S.; Chan, H.L.Y. Transient elastography. J. Gastroenterol. Hepatol. 2010, 25, 1726–1731. [Google Scholar] [CrossRef]
  17. Sigrist, R.M.; Liau, J.; El Kaffas, A.; Chammas, M.C.; Willmann, J.K. Ultrasound elastography: Review of techniques and clinical applications. Theranostics 2017, 7, 1303. [Google Scholar] [CrossRef]
  18. Sack, I. Magnetic resonance elastography from fundamental soft-tissue mechanics to diagnostic imaging. Nat. Rev. Phys. 2023, 5, 25–42. [Google Scholar] [CrossRef]
  19. Goodbrake, C.; Li, D.S.; Aghakhani, H.; Contreras, A.; Reece, G.P.; Markey, M.K.; Sacks, M.S. On the three-dimensional mechanical behavior of human breast tissue. Ann. Biomed. Eng. 2022, 50, 601–613. [Google Scholar] [CrossRef] [PubMed]
  20. Yao, Z.; Luo, T.; Dong, Y.; Jia, X.; Deng, Y.; Wu, G.; Zhu, Y.; Zhang, J.; Liu, J.; Yang, L. Virtual elastography ultrasound via generative adversarial network for breast cancer diagnosis. Nat. Commun. 2023, 14, 788. [Google Scholar] [CrossRef]
  21. Idilman, I.S.; Li, J.; Yin, M.; Venkatesh, S.K. MR elastography of liver: Current status and future perspectives. Abdom. Radiol. 2020, 45, 3444–3462. [Google Scholar] [CrossRef]
  22. Sumbul, H.E.; Avci, B.S.; Bankir, M.; Pekoz, B.C.; Gulumsek, E.; Koc, A.S. Ovarian stiffness is significantly increased in polycystic ovary syndrome and related with anti-mullerian hormone: A point shear wave elastography study. Ultrasound Q. 2022, 38, 83–88. [Google Scholar] [CrossRef]
  23. Yoo, J.W.; Koo, K.C.; Chung, B.H.; Lee, K.S. Role of the elastography strain ratio using transrectal ultrasonography in the diagnosis of prostate cancer and clinically significant prostate cancer. Sci. Rep. 2022, 12, 21171. [Google Scholar] [CrossRef]
  24. Caenen, A.; Pernot, M.; Nightingale, K.R.; Voigt, J.-U.; Vos, H.J.; Segers, P.; D’hooge, J. Assessing cardiac stiffness using ultrasound shear wave elastography. Phys. Med. Biol. 2022, 67, 02TR01. [Google Scholar] [CrossRef] [PubMed]
  25. Islam, M.T.; Tang, S.; Liverani, C.; Saha, S.; Tasciotti, E.; Righetti, R. Non-invasive imaging of Young’s modulus and Poisson’s ratio in cancers in vivo. Sci. Rep. 2020, 10, 7266. [Google Scholar] [CrossRef]
  26. Sowinski, D.R.; McGarry, M.D.; Van Houten, E.E.; Gordon-Wylie, S.; Weaver, J.B.; Paulsen, K.D. Poroelasticity as a model of soft tissue structure: Hydraulic permeability reconstruction for magnetic resonance elastography in silico. Front. Phys. 2021, 8, 617582. [Google Scholar] [CrossRef] [PubMed]
  27. Jin, Z.; Zhou, Y.; Shen, M.; Wang, Y.; Lu, F.; Zhu, D. Assessment of corneal viscoelasticity using elastic wave optical coherence elastography. J. Biophotonics 2020, 13, e201960074. [Google Scholar] [CrossRef]
  28. Samani, A.; Zubovits, J.; Plewes, D. Elastic moduli of normal and pathological human breast tissues: An inversion-technique-based investigation of 169 samples. Phys. Med. Biol. 2007, 52, 1565. [Google Scholar] [CrossRef] [PubMed]
  29. Kheirkhah, N.; Dempsey, S.; Sadeghi-Naini, A.; Samani, A. A novel tissue mechanics-based method for improved motion tracking in quasi-static ultrasound elastography. Med. Phys. 2023, 50, 2176–2194. [Google Scholar] [CrossRef]
  30. Poul, S.S.; Ormachea, J.; Gary, R.G.; Parker, K.J. Comprehensive experimental assessments of rheological models’ performance in elastography of soft tissues. Acta Biomater. 2022, 146, 259–273. [Google Scholar] [CrossRef]
  31. Rosen, D.P.; Jiang, J. A comparison of hyperelastic constitutive models applicable to shear wave elastography (SWE) data in tissue-mimicking materials. Phys. Med. Biol. 2019, 64, 055014. [Google Scholar] [CrossRef]
  32. Zheng, Y.; Jiang, Y.; Cao, Y. A porohyperviscoelastic model for the shear wave elastography of the liver. J. Mech. Phys. Solids 2021, 150, 104339. [Google Scholar] [CrossRef]
  33. Gentile, F.; Sakamoto, J.; Righetti, R.; Decuzzi, P.; Ferrari, M. A doublet mechanics model for the ultrasound characterization of malignant tissues. J. Biomed. Sci. Eng. 2011, 4, 362. [Google Scholar] [CrossRef]
  34. Tang, S.; Yang, X.; Shajudeen, P.; Sears, C.; Taraballi, F.; Weiner, B.; Tasciotti, E.; Dollahon, D.; Park, H.; Righetti, R. A CNN-based method to reconstruct 3-D spine surfaces from US images in vivo. Med. Image Anal. 2021, 74, 102221. [Google Scholar] [CrossRef] [PubMed]
  35. Islam, M.T.; Righetti, R. An analytical poroelastic model of a spherical tumor embedded in normal tissue under creep compression. J. Biomech. 2019, 89, 48–56. [Google Scholar] [CrossRef]
  36. Alekya, B.; Rao, S.; Pandya, H.J. Engineering approaches for characterizing soft tissue mechanical properties: A review. Clin. Biomech. 2019, 69, 127–140. [Google Scholar]
  37. Xiang, P.; Liew, K.M. A computational framework for transverse compression of microtubules based on a higher-order Cauchy–Born rule. Comput. Methods Appl. Mech. Eng. 2013, 254, 14–30. [Google Scholar] [CrossRef]
  38. Xiang, P.; Zhang, L.; Liew, K. A third-order Cauchy-Born rule for modeling of microtubules based on the element-free framework. Compos. Struct. 2017, 161, 215–226. [Google Scholar] [CrossRef]
  39. Dwairy, M.; Reddy, J.; Righetti, R. Predicting stress and interstitial fluid pressure in tumors based on biphasic theory. Comput. Biol. Med. 2023, 167, 107651. [Google Scholar] [CrossRef]
  40. Fuhs, T.; Wetzel, F.; Fritsch, A.W.; Li, X.; Stange, R.; Pawlizak, S.; Kießling, T.R.; Morawetz, E.; Grosser, S.; Sauer, F. Rigid tumours contain soft cancer cells. Nat. Phys. 2022, 18, 1510–1519. [Google Scholar] [CrossRef]
  41. Ghayesh, M.H.; Farajpour, A. A review on the mechanics of functionally graded nanoscale and microscale structures. Int. J. Eng. Sci. 2019, 137, 8–36. [Google Scholar] [CrossRef]
  42. Eringen, A.C.; Wegner, J. Nonlocal continuum field theories. Appl. Mech. Rev. 2003, 56, B20–B22. [Google Scholar] [CrossRef]
  43. Eringen, A.C. On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves. J. Appl. Phys. 1983, 54, 4703–4710. [Google Scholar] [CrossRef]
  44. Lim, C.; Zhang, G.; Reddy, J. A higher-order nonlocal elasticity and strain gradient theory and its applications in wave propagation. J. Mech. Phys. Solids 2015, 78, 298–313. [Google Scholar] [CrossRef]
  45. Farajpour, A.; Yazdi, M.H.; Rastgoo, A.; Mohammadi, M. A higher-order nonlocal strain gradient plate model for buckling of orthotropic nanoplates in thermal environment. Acta Mech. 2016, 227, 1849–1867. [Google Scholar] [CrossRef]
  46. Ghandourah, E.E.; Daikh, A.A.; Khatir, S.; Alhawsawi, A.M.; Banoqitah, E.M.; Eltaher, M.A. A Dynamic Analysis of Porous Coated Functionally Graded Nanoshells Rested on Viscoelastic Medium. Mathematics 2023, 11, 2407. [Google Scholar] [CrossRef]
  47. Barati, A.; Hadi, A.; Nejad, M.Z.; Noroozi, R. On vibration of bi-directional functionally graded nanobeams under magnetic field. Mech. Based Des. Struct. Mach. 2022, 50, 468–485. [Google Scholar] [CrossRef]
  48. Zarezadeh, E.; Hosseini, V.; Hadi, A. Torsional vibration of functionally graded nano-rod under magnetic field supported by a generalized torsional foundation based on nonlocal elasticity theory. Mech. Based Des. Struct. Mach. 2020, 48, 480–495. [Google Scholar] [CrossRef]
  49. Emadi, M.; Nejad, M.Z.; Ziaee, S.; Hadi, A. Buckling analysis of arbitrary two-directional functionally graded nano-plate based on nonlocal elasticity theory using generalized differential quadrature method. Steel Compos. Struct. Int. J. 2021, 39, 565–581. [Google Scholar]
  50. Li, X.; Li, L.; Hu, Y.; Ding, Z.; Deng, W. Bending, buckling and vibration of axially functionally graded beams based on nonlocal strain gradient theory. Compos. Struct. 2017, 165, 250–265. [Google Scholar] [CrossRef]
  51. Melaibari, A.; Abdelrahman, A.A.; Hamed, M.A.; Abdalla, A.W.; Eltaher, M.A. Dynamic analysis of a piezoelectrically layered perforated nonlocal strain gradient nanobeam with flexoelectricity. Mathematics 2022, 10, 2614. [Google Scholar] [CrossRef]
  52. Civalek, Ö.; Demir, C. A simple mathematical model of microtubules surrounded by an elastic matrix by nonlocal finite element method. Appl. Math. Comput. 2016, 289, 335–352. [Google Scholar] [CrossRef]
  53. Demir, C.; Civalek, Ö. Torsional and longitudinal frequency and wave response of microtubules based on the nonlocal continuum and nonlocal discrete models. Appl. Math. Model. 2013, 37, 9355–9367. [Google Scholar] [CrossRef]
  54. DEMİR, Ç.; Civalek, Ö. Nonlocal deflection of microtubules under point load. Int. J. Eng. Appl. Sci. 2015, 7, 33–39. [Google Scholar] [CrossRef]
  55. Farajpour, A.; Rastgoo, A. Influence of carbon nanotubes on the buckling of microtubule bundles in viscoelastic cytoplasm using nonlocal strain gradient theory. Results Phys. 2017, 7, 1367–1375. [Google Scholar] [CrossRef]
  56. Karamanli, A.; Vo, T.P.; Civalek, O. Higher order finite element models for transient analysis of strain gradient functionally graded microplates. Eur. J. Mech.-A/Solids 2023, 99, 104933. [Google Scholar] [CrossRef]
  57. Uzun, B.; Civalek, Ö.; Yaylı, M.Ö. Torsional and axial vibration of restrained saturated nanorods via strain gradient elasticity. Arch. Appl. Mech. 2023, 93, 1605–1630. [Google Scholar] [CrossRef]
  58. Tocci Monaco, G.; Fantuzzi, N.; Fabbrocino, F.; Luciano, R. Trigonometric solution for the bending analysis of magneto-electro-elastic strain gradient nonlocal nanoplates in hygro-thermal environment. Mathematics 2021, 9, 567. [Google Scholar] [CrossRef]
  59. Archer, M.; Dasari, P.; Evdokiou, A.; Ingman, W.V. Biological mechanisms and therapeutic opportunities in mammographic density and breast cancer risk. Cancers 2021, 13, 5391. [Google Scholar] [CrossRef] [PubMed]
  60. Ingman, W.V.; Richards, B.; Street, J.M.; Carter, D.; Rickard, M.; Stone, J.; Dasari, P. Breast density notification: An Australian perspective. J. Clin. Med. 2020, 9, 681. [Google Scholar] [CrossRef]
  61. McCormack, V.A.; dos Santos Silva, I. Breast density and parenchymal patterns as markers of breast cancer risk: A meta-analysis. Cancer Epidemiol. Biomark. Prev. 2006, 15, 1159–1169. [Google Scholar] [CrossRef]
  62. Huo, C.; Chew, G.; Britt, K.; Ingman, W.; Henderson, M.; Hopper, J.; Thompson, E. Mammographic density—A review on the current understanding of its association with breast cancer. Breast Cancer Res. Treat. 2014, 144, 479–502. [Google Scholar] [CrossRef]
  63. Archer, M.; Dasari, P.; Walsh, D.; Britt, K.L.; Evdokiou, A.; Ingman, W.V. Immune regulation of mammary fibroblasts and the impact of mammographic density. J. Clin. Med. 2022, 11, 799. [Google Scholar] [CrossRef]
  64. Litin, S.C. Mayo Clinic Family Health Book; HarperCollins: New York, NY, USA, 2005. [Google Scholar]
  65. Vasei, N.; Shishegar, A.; Ghalkhani, F.; Darvishi, M. Fat necrosis in the Breast: A systematic review of clinical. Lipids Health Dis. 2019, 18, 139. [Google Scholar] [CrossRef] [PubMed]
  66. Singh, P.; Manure, A.; Singh, P.; Manure, A. Introduction to tensorflow 2.0. In Learn TensorFlow 2.0: Implement Machine Learning and Deep Learning Models with Python; Springer: Berlin/Heidelberg, Germany, 2020; pp. 1–24. [Google Scholar]
  67. Zouch, W.; Sagga, D.; Echtioui, A.; Khemakhem, R.; Ghorbel, M.; Mhiri, C.; Hamida, A.B. Detection of COVID-19 from CT and chest X-ray images using deep learning models. Ann. Biomed. Eng. 2022, 50, 825–835. [Google Scholar] [CrossRef] [PubMed]
  68. Joseph, F.J.J.; Nonsiri, S.; Monsakul, A. Keras and TensorFlow: A hands-on experience. Advanced Deep Learning for Engineers and Scientists: A Practical Approach; Springer: Berlin/Heidelberg, Germany, 2021; pp. 85–111. [Google Scholar]
  69. Sanchez, A.Z.; Mills, C.; Scurr, J.C. Estimating breast mass-density: A retrospective analysis of radiological data. Breast J. 2017, 23, 237–239. [Google Scholar] [CrossRef] [PubMed]
  70. Mills, C.; Sanchez, A.; Scurr, J. Estimating the gravity induced three dimensional deformation of the breast. J. Biomech. 2016, 49, 4134–4137. [Google Scholar] [CrossRef] [PubMed]
  71. Bharat, A.; Gao, F.; Margenthaler, J.A. Tumor characteristics and patient outcomes are similar between invasive lobular and mixed invasive ductal/lobular breast cancers but differ from pure invasive ductal breast cancers. Am. J. Surg. 2009, 198, 516–519. [Google Scholar] [CrossRef]
Figure 1. A flowchart to depict the present research workflow: (1) in-plane waves are propagated within various types of healthy, benign, and malignant breast conditions; (2) a higher-order nonlocal model is developed to extract the in-plane wave characteristics; (3) the mathematical model is integrated with a neural network algorithm to improve the accuracy and flexibility of the model for clinical applications; (4) a finite element comparison study is conducted to verify the validity of the model.
Figure 1. A flowchart to depict the present research workflow: (1) in-plane waves are propagated within various types of healthy, benign, and malignant breast conditions; (2) a higher-order nonlocal model is developed to extract the in-plane wave characteristics; (3) the mathematical model is integrated with a neural network algorithm to improve the accuracy and flexibility of the model for clinical applications; (4) a finite element comparison study is conducted to verify the validity of the model.
Mathematics 11 04804 g001
Figure 2. A schematic representation of human breast sections under in-plane elastic wave propagation. An ultrasound transducer can be used to create in-plane waves. The direction of vibrations caused by the transducer is placed within the x-y plane. There is no significant lateral displacement (i.e., perpendicular to the plane of the breast sections).
Figure 2. A schematic representation of human breast sections under in-plane elastic wave propagation. An ultrasound transducer can be used to create in-plane waves. The direction of vibrations caused by the transducer is placed within the x-y plane. There is no significant lateral displacement (i.e., perpendicular to the plane of the breast sections).
Mathematics 11 04804 g002
Figure 3. (a) Machine learning overview. There are five main processes involved in the deep learning algorithm: (1) scaling process: the dataset is scaled to ensure that all data features are in the same range; (2) train/test split: the scaled dataset is divided into two subsets as training (70%) and test (30%); (3) model fitting: training process; (4) model evaluation: the mean square errors of both test and training data are monitored; (5) hyperparameter adjustment: the model parameters are adjusted to enhance accuracy; (b) architecture of the deep neural network with one flatten layer (FL) at the beginning, a couple of hidden dense layers (HDLs) in the middle, and an output layer at the end. nk indicates the number of neurons in the kth HDL. In each neuron, the input quantities are first aggregated through summation and then an activation function ( f A ) is applied.
Figure 3. (a) Machine learning overview. There are five main processes involved in the deep learning algorithm: (1) scaling process: the dataset is scaled to ensure that all data features are in the same range; (2) train/test split: the scaled dataset is divided into two subsets as training (70%) and test (30%); (3) model fitting: training process; (4) model evaluation: the mean square errors of both test and training data are monitored; (5) hyperparameter adjustment: the model parameters are adjusted to enhance accuracy; (b) architecture of the deep neural network with one flatten layer (FL) at the beginning, a couple of hidden dense layers (HDLs) in the middle, and an output layer at the end. nk indicates the number of neurons in the kth HDL. In each neuron, the input quantities are first aggregated through summation and then an activation function ( f A ) is applied.
Mathematics 11 04804 g003
Figure 4. The first four mode shapes of vibrations induced by in-plane wave propagation within the human breast tissue. All patterns of in-plane elastic waves were plotted for the DCIS using the finite element simulations. Similar patterns were found for other breast lesions. However, the magnitudes of vibration amplitudes were different. The total in-plane displacement fields in terms of m within the breast tissue sections were also plotted in this figure using colour bars. As the mode number of the wave propagation increased, the pattern became more sophisticated with more local sharp changes in the in-plane displacement.
Figure 4. The first four mode shapes of vibrations induced by in-plane wave propagation within the human breast tissue. All patterns of in-plane elastic waves were plotted for the DCIS using the finite element simulations. Similar patterns were found for other breast lesions. However, the magnitudes of vibration amplitudes were different. The total in-plane displacement fields in terms of m within the breast tissue sections were also plotted in this figure using colour bars. As the mode number of the wave propagation increased, the pattern became more sophisticated with more local sharp changes in the in-plane displacement.
Mathematics 11 04804 g004
Figure 5. (a) Phase velocity of in-plane waves propagated within the human dense breast and DCIS versus the wave number; (b) frequency of in-plane vibrations induced by the wave propagation within the human dense breast and DCIS. The first mode of in-plane vibrations was excited to extract the wave propagation characteristics of the breast lesion. A scale-dependent model of higher-order nonlocal elasticity was used for calculations. Phase velocity of in-plane elastic waves within (c) a human dense breast and (d) DCIS versus sample size for different models including nonlocal, strain gradient, classical and a combination of them.
Figure 5. (a) Phase velocity of in-plane waves propagated within the human dense breast and DCIS versus the wave number; (b) frequency of in-plane vibrations induced by the wave propagation within the human dense breast and DCIS. The first mode of in-plane vibrations was excited to extract the wave propagation characteristics of the breast lesion. A scale-dependent model of higher-order nonlocal elasticity was used for calculations. Phase velocity of in-plane elastic waves within (c) a human dense breast and (d) DCIS versus sample size for different models including nonlocal, strain gradient, classical and a combination of them.
Mathematics 11 04804 g005
Figure 6. (a) Phase velocity of in-plane waves propagated within the human dense breast and DCIS versus the wave number; (b) frequency of in-plane vibrations induced by the wave propagation within the human dense breast and DCIS. The fibroglandular tissue was assumed as the main building block of the dense breast while early breast cancer was mainly found in the form of the DCIS. To estimate the wave propagation characteristics of breast lesions, the second mode of in-plane vibrations was excited.
Figure 6. (a) Phase velocity of in-plane waves propagated within the human dense breast and DCIS versus the wave number; (b) frequency of in-plane vibrations induced by the wave propagation within the human dense breast and DCIS. The fibroglandular tissue was assumed as the main building block of the dense breast while early breast cancer was mainly found in the form of the DCIS. To estimate the wave propagation characteristics of breast lesions, the second mode of in-plane vibrations was excited.
Mathematics 11 04804 g006
Figure 7. (a) Role of zeroth-order nonlocal parameter (NLP0) in the phase velocity of in-plane waves propagated within the intermediate-grade IDC as the most common type of breast cancer, which starts from the breast milk ducts and spreads to other tissues. The elastic properties of the intermediate IDC were adopted from the available clinical data [28]. (b) Role of first-order nonlocal parameter (NLP1) in the phase velocity of in-plane waves propagated within the intermediate-grade IDC. (c) Role of the strain gradient parameter (SGP) in the phase velocity of in-plane waves propagated within the intermediate-grade IDC tissue. The IDC is known as the most frequently encountered breast cancer, which accounts for about 75–80% of all breast cancer types that exhibit invasiveness [71]. The thickness of the tissue sections was 1 cm while the tissue length and width were 10 cm.
Figure 7. (a) Role of zeroth-order nonlocal parameter (NLP0) in the phase velocity of in-plane waves propagated within the intermediate-grade IDC as the most common type of breast cancer, which starts from the breast milk ducts and spreads to other tissues. The elastic properties of the intermediate IDC were adopted from the available clinical data [28]. (b) Role of first-order nonlocal parameter (NLP1) in the phase velocity of in-plane waves propagated within the intermediate-grade IDC. (c) Role of the strain gradient parameter (SGP) in the phase velocity of in-plane waves propagated within the intermediate-grade IDC tissue. The IDC is known as the most frequently encountered breast cancer, which accounts for about 75–80% of all breast cancer types that exhibit invasiveness [71]. The thickness of the tissue sections was 1 cm while the tissue length and width were 10 cm.
Mathematics 11 04804 g007
Figure 8. (a) Frequency and (b) phase velocity of in-plane waves propagated within various breast lesions estimated based on the nonlocal elasticity model, in which the zeroth-order nonlocal parameter was 2 cm while the remaining scale parameters were ignored. Here ‘**’ means p ≤ 0.01, ‘***’ represents p ≤ 0.001, and ‘****’ indicates p ≤ 0.0001.
Figure 8. (a) Frequency and (b) phase velocity of in-plane waves propagated within various breast lesions estimated based on the nonlocal elasticity model, in which the zeroth-order nonlocal parameter was 2 cm while the remaining scale parameters were ignored. Here ‘**’ means p ≤ 0.01, ‘***’ represents p ≤ 0.001, and ‘****’ indicates p ≤ 0.0001.
Mathematics 11 04804 g008
Figure 9. (a) Frequency and (b) phase velocity of in-plane waves propagated within various breast lesions estimated based on the strain gradient elasticity model, in which both zeroth- and first-order nonlocal parameters were zero while the remaining scale parameter was 2 cm. Here ‘**’ means p ≤ 0.01, ‘***’ represents p ≤ 0.001, and ‘****’ indicates p ≤ 0.0001.
Figure 9. (a) Frequency and (b) phase velocity of in-plane waves propagated within various breast lesions estimated based on the strain gradient elasticity model, in which both zeroth- and first-order nonlocal parameters were zero while the remaining scale parameter was 2 cm. Here ‘**’ means p ≤ 0.01, ‘***’ represents p ≤ 0.001, and ‘****’ indicates p ≤ 0.0001.
Mathematics 11 04804 g009
Table 1. Finite element results of dimensionless phase velocities of breast tissues. Three different human breast samples from benign, healthy, and cancerous conditions are considered. The dimensionless phase velocity is defined as (phase velocity) × k/(6.7454 ω 1 ), where ω 1 is the first vibration frequency of FGT and 6.7454 is the maximum observed dimensionless phase velocity.
Table 1. Finite element results of dimensionless phase velocities of breast tissues. Three different human breast samples from benign, healthy, and cancerous conditions are considered. The dimensionless phase velocity is defined as (phase velocity) × k/(6.7454 ω 1 ), where ω 1 is the first vibration frequency of FGT and 6.7454 is the maximum observed dimensionless phase velocity.
Human
Breast Lesion
Elasticity Modulus (KPa) [28]Mode
Number
In-Plane
Wave Number
Finite Element SimulationPresent
Analytical Model
Fibroglandular tissue3.24n = 10.33330.14800.1482
n = 20.66670.29670.2965
n = 31.000.44500.4447
DCIS16.38n = 10.33330.33310.3333
n = 20.66670.66700.6667
n = 31.001.001.00
Fat necrosis4.45n = 10.33330.17410.1737
n = 20.66670.34770.3475
n = 31.000.52100.5212
Table 2. The mean square errors (MSEs) of both training and testing datasets for three different architectures of deep learning algorithms: (1) FL-3HDL(512,128,10): one flatten layer (FL) and three hidden dense layers (HDLs) with 512, 128, and 10 neurons; (2) FL-2HDL(128,10): one FL and two HDLs with 128 and 10 neurons; (3) FL-1HDL(10): one FL and one HDL with 10 neurons.
Table 2. The mean square errors (MSEs) of both training and testing datasets for three different architectures of deep learning algorithms: (1) FL-3HDL(512,128,10): one flatten layer (FL) and three hidden dense layers (HDLs) with 512, 128, and 10 neurons; (2) FL-2HDL(128,10): one FL and two HDLs with 128 and 10 neurons; (3) FL-1HDL(10): one FL and one HDL with 10 neurons.
Deep Learning ModelArchitecture Design
(FL-nHDL) *
Number of EpochsMSE
(Train Data)
MSE (Test Data)
1FL-3HDL(512,128,10)10.43980.2520
100.10700.0621
503.3644 × 10−62.2851 × 10−6
1002.0467 × 10−63.1722 × 10−6
5004.7098 × 10−71.0561 × 10−6
2FL-2HDL(128,10)10.46800.3417
100.15830.0991
500.01950.0106
1000.00150.0015
5006.4016 × 10−58.3080 × 10−5
3FL-1HDL(10)10.47320.3897
100.36270.2816
500.14480.0907
1000.09130.0596
5000.00190.0019
* FL-nHDL (l1, l2, l3, …, ln): a deep learning model with one flatten layer (FL) and n hidden dense layers (nHDLs) in which li indicates the number of units (neurons) in the ith hidden dense layer.
Table 3. The predicted dimensionless wave velocities obtained using the deep learning algorithm compared to the actual data calculated using the analytical model of the higher-order nonlocal elasticity. Three different architectures including (1) FL-3HDL(512,128,10), (2) FL-2HDL(128,10), and (3) FL-1HDL(10) were built. To standardise tissue features, a scaling pre-processing method, which included dividing the quantities by their maximum values, was used.
Table 3. The predicted dimensionless wave velocities obtained using the deep learning algorithm compared to the actual data calculated using the analytical model of the higher-order nonlocal elasticity. Three different architectures including (1) FL-3HDL(512,128,10), (2) FL-2HDL(128,10), and (3) FL-1HDL(10) were built. To standardise tissue features, a scaling pre-processing method, which included dividing the quantities by their maximum values, was used.
Deep Learning ModelArchitecture
(FL-nHDL) *
Dimensionless Wave NumberPredicted
Dimensionless Wave Velocity
Actual
Dimensionless Wave Velocity
1FL-3HDL(512,128,10)0.810.55060.5497
0.390.71340.7140
0.180.89990.90
0.470.66310.6621
0.280.80700.8056
2FL-2HDL(128,10)0.810.55640.5497
0.390.72570.7140
0.180.89520.90
0.470.66110.6621
0.280.81450.8056
3FL-1HDL(10)0.810.53710.5497
0.390.76090.7140
0.180.87270.90
0.470.71830.6621
0.280.81950.8056
* FL-nHDL(l1, l2, l3, …, ln): a deep learning model with one flatten layer (FL) and n hidden dense layers (nHDLs) in which li indicates the number of units (neurons) in the ith hidden dense layer.
Table 4. A comparison study conducted on different machine learning models including the deep neural network algorithm, Random Forest model, and Lasso, linear, Ridge, and ElasticNet regression models. The neural network model was made of one flatten layer, three hidden layers with 512, 128, and 10 neurons, as well as one dense layer at the end with one neuron.
Table 4. A comparison study conducted on different machine learning models including the deep neural network algorithm, Random Forest model, and Lasso, linear, Ridge, and ElasticNet regression models. The neural network model was made of one flatten layer, three hidden layers with 512, 128, and 10 neurons, as well as one dense layer at the end with one neuron.
Dimensionless Wave NumberActual
Data
Machine Learning Model
Neural
Network
Random ForestLinearLassoRidgeElasticNet
0.030.99700.99710.99770.95600.95070.94820.9427
0.110.95830.95840.96170.91300.90870.90660.9021
0.20.88130.88010.87600.86480.86150.85990.8564
0.330.76130.76110.75780.79500.79320.79240.7905
0.370.7290.72820.72400.77360.77220.77160.7702
0.450.67390.67350.67980.73060.73020.73000.7296
0.520.63580.63550.63990.69310.69350.69370.6941
0.570.61370.61320.61490.66630.66720.66770.6687
0.660.58280.58220.58170.61800.62000.62090.6230
0.740.56280.56190.56360.57510.57800.57940.5825
0.830.54650.54780.54630.52680.53070.53260.5368
0.850.54360.54480.54240.51600.52020.52220.5267
0.930.53380.53290.53480.47310.47820.48070.4861
0.950.53170.52990.53110.46240.46770.47030.4759
0.980.52890.52540.52930.44630.45200.45470.4607
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.

Share and Cite

MDPI and ACS Style

Farajpour, A.; Ingman, W.V. In-Plane Wave Propagation Analysis of Human Breast Lesions Using a Higher-Order Nonlocal Model and Deep Learning. Mathematics 2023, 11, 4804. https://doi.org/10.3390/math11234804

AMA Style

Farajpour A, Ingman WV. In-Plane Wave Propagation Analysis of Human Breast Lesions Using a Higher-Order Nonlocal Model and Deep Learning. Mathematics. 2023; 11(23):4804. https://doi.org/10.3390/math11234804

Chicago/Turabian Style

Farajpour, Ali, and Wendy V. Ingman. 2023. "In-Plane Wave Propagation Analysis of Human Breast Lesions Using a Higher-Order Nonlocal Model and Deep Learning" Mathematics 11, no. 23: 4804. https://doi.org/10.3390/math11234804

APA Style

Farajpour, A., & Ingman, W. V. (2023). In-Plane Wave Propagation Analysis of Human Breast Lesions Using a Higher-Order Nonlocal Model and Deep Learning. Mathematics, 11(23), 4804. https://doi.org/10.3390/math11234804

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop