
This article presents the results of research on the effects of electrolyte plasma hardening on the structure, phase composition, tribological, and mechanical properties of medium-carbon structural steel 45, which is widely used in the manufacturing of tools and machine parts. Hardening experiments were conducted using an electrolyte plasma hardening setup with electrolytes varying in sodium carbonate (Na2CO3) concentration in distilled water (15%, 20%, and 25%). With a consistent heating duration of 4 s during quenching, significant phase changes in the steel's microstructure were observed, enhancing hardness and wear resistance. The transformation of the initial structure of steel 45, which consists of ferrite and pearlite into martensite on the surface of the samples, led to an increase in microhardness up to 506–690 HV01. This value is 2.5–3.5 times higher compared to the untreated sample, and the thickness of the hardened layer reached up to 3.2 mm. Additionally, wear volume measurements showed that after electrolyte plasma hardening, the wear resistance of the samples increased by 1.3–1.5 times (2.01 × 10−4, 2.26 × 10−4 m3). The obtained results on the changes in microstructure and mechano-tribological properties of steel 45 confirm the potential of electrolyte plasma hardening technology for improving operational characteristics and extending the service life of heavily loaded and critical machine parts.
Citation: Bauyrzhan Rakhadilov, Rinat Kussainov, Aisulu Kalitova, Zarina Satbayeva, Aibek Shynarbek. The impact of technological parameters of electrolytic-plasma treatment on the changes in the mechano-tribological properties of steel 45[J]. AIMS Materials Science, 2024, 11(4): 666-683. doi: 10.3934/matersci.2024034
[1] | Md. Mehedi Hassan, Rezuana Haque, Sheikh Mohammed Shariful Islam, Hossam Meshref, Roobaea Alroobaea, Mehedi Masud, Anupam Kumar Bairagi . NeuroWave-Net: Enhancing epileptic seizure detection from EEG brain signals via advanced convolutional and long short-term memory networks. AIMS Bioengineering, 2024, 11(1): 85-109. doi: 10.3934/bioeng.2024006 |
[2] | Yi Deng, Zhiguo Wang, Xiaohui Li, Yu Lei, Owen Omalley . Advancing biomedical engineering through a multi-modal sensor fusion system for enhanced physical training. AIMS Bioengineering, 2023, 10(4): 364-383. doi: 10.3934/bioeng.2023022 |
[3] | Kayode Oshinubi, Augustina Amakor, Olumuyiwa James Peter, Mustapha Rachdi, Jacques Demongeot . Approach to COVID-19 time series data using deep learning and spectral analysis methods. AIMS Bioengineering, 2022, 9(1): 1-21. doi: 10.3934/bioeng.2022001 |
[4] | Kuna Dhananjay Rao, Mudunuru Satya Dev Kumar, Paidi Pavani, Darapureddy Akshitha, Kagitha Nagamaleswara Rao, Hafiz Tayyab Rauf, Mohamed Sharaf . Cardiovascular disease prediction using hyperparameters-tuned LSTM considering COVID-19 with experimental validation. AIMS Bioengineering, 2023, 10(3): 265-282. doi: 10.3934/bioeng.2023017 |
[5] | Abdulmajeed Alsufyani . Performance comparison of deep learning models for MRI-based brain tumor detection. AIMS Bioengineering, 2025, 12(1): 1-21. doi: 10.3934/bioeng.2025001 |
[6] | Mohammad Z. Ansari, Anil K. Nirala . Assessment of Fevicol (adhesive) Drying Process through Dynamic Speckle Techniques. AIMS Bioengineering, 2015, 2(2): 49-59. doi: 10.3934/bioeng.2015.2.49 |
[7] | Torsten Schüling, Alina Eilers, Thomas Scheper, Johanna Walter . Aptamer-based lateral flow assays. AIMS Bioengineering, 2018, 5(2): 78-102. doi: 10.3934/bioeng.2018.2.78 |
[8] | Abhimanu Singh, Smita Jain . Enhanced brain tumor detection from brain MRI images using convolutional neural networks. AIMS Bioengineering, 2025, 12(2): 215-224. doi: 10.3934/bioeng.2025010 |
[9] | Poovayar Priya Mohan, Ezhilarasan Murugesan . A novel approach for detecting the presence of a lymphatic rosary in the iris using a 3-stage algorithm. AIMS Bioengineering, 2024, 11(4): 506-526. doi: 10.3934/bioeng.2024023 |
[10] | Lamia Fatiha KAZI TANI, Mohammed Yassine KAZI TANI, Benamar KADRI . Gas-Net: A deep neural network for gastric tumor semantic segmentation. AIMS Bioengineering, 2022, 9(3): 266-282. doi: 10.3934/bioeng.2022018 |
This article presents the results of research on the effects of electrolyte plasma hardening on the structure, phase composition, tribological, and mechanical properties of medium-carbon structural steel 45, which is widely used in the manufacturing of tools and machine parts. Hardening experiments were conducted using an electrolyte plasma hardening setup with electrolytes varying in sodium carbonate (Na2CO3) concentration in distilled water (15%, 20%, and 25%). With a consistent heating duration of 4 s during quenching, significant phase changes in the steel's microstructure were observed, enhancing hardness and wear resistance. The transformation of the initial structure of steel 45, which consists of ferrite and pearlite into martensite on the surface of the samples, led to an increase in microhardness up to 506–690 HV01. This value is 2.5–3.5 times higher compared to the untreated sample, and the thickness of the hardened layer reached up to 3.2 mm. Additionally, wear volume measurements showed that after electrolyte plasma hardening, the wear resistance of the samples increased by 1.3–1.5 times (2.01 × 10−4, 2.26 × 10−4 m3). The obtained results on the changes in microstructure and mechano-tribological properties of steel 45 confirm the potential of electrolyte plasma hardening technology for improving operational characteristics and extending the service life of heavily loaded and critical machine parts.
A large variety of models of interactions between agents have been developed within the framework of thermostatted kinetic theory [1,2] and notably applied to traffic [3], swarming [4] and cell competition [5,6,7,8]. The interactions between immune system and cancer are known to be complex. The competition between immune system cells and cancer cells can be compared to the interactions between immune system cells and viruses. In both cases, pathogenic cells induce immune system responses [9,10,11,12]. Different cell types are involved in the defense mechanism. Three types of cells are commonly recognized as being involved in the immune response, antigen-presenting cells, T cells and B cells. Dentritic cells ingest cancer cells, isolate antigens and present them to T cells, which leads to their activation and proliferation [13]. Activated T cells mutate into killer T cells, which destroy cancer cells. Like innate dendritic cells, B cells are antigen-presenting cells developed by the adaptive immune system. B cells produce antibodies, specific to each pathogen, which bind to the surface of an invading cell and mark it for destruction by T cells. When a cancer cell is killed, additional antigens are released, which stimulates further activation, i.e. learning of new T and B cells. In addition, regulatory T cells are produced to control the response of immune system. They adjust the level of antibody secretion, i.e. the regulation of B and T cell production and activation. While T cells are activated, cancer cells mutate and acquire the ability to avoid detection and destruction by the immune system. Specifically, cancer cells express proteins on their surface that induce the inactivation of immune cells, thus promoting the proliferation of cancer cells [14,15]. Through cellular interactions, cancer cells learn not to be detected, misleading the regulatory T cells, which in turn limits the production of immune system cells and let the cancer cells proliferate [16,17]. In the model, the learning processes and activation of both immune system cells and cancer cells is reproduced by the increase of a real quantity called activity and carried by each cell [2].
Stimulating the immune system of a patient can lead to disappointing results and conversely, chemotherapy may boost immune system response [18]. For some cancers, immunotherapy leads to a long remission followed by a fatal outcome, reflecting the complex relationship between the immune system and the cancer. This phenomenon, known as the three Es characterizes the following three phases commonly observed in a patient, elimination, equilibrium and escape of cancer, introduced in the theory of immunosurveillance as early as the mid-20th century [9,19]. The initial elimination of cancer cells is only apparent and corresponds to the attack by killer T cells, before the cancer cells have learned to mislead the immune system. During the pseudo-equilibrium phase, a small number of cancer cells remain and learn from the immune system cells. Eventually, the cancer cells are able to proliferate and the cancer escapes the control of the immune system. C-Reactive Protein (CRP) levels have been used as markers of inflammation in response to antigen exposure in relation to the progression of many cancers, such as melanoma and lung cancer [20]. Oscillations of CRP concentration in advanced cancer patients have been observed and interpreted in terms of immune regulatory cycles [21,22]. The oscillations of CRP concentration reveal the antagonistic effects of immune system stimulation on a tumor. Many macroscopic models of interactions between cancer and immune system leading to oscillating solutions for concentrations after a Hopf bifurcation can be found in the recent literature [23,24,25,26], in particular after introduction of one or more delays in the equations [27,28], but without direct connection with a description at cell scale.
We recently designed a minimal model of interactions between cancer and immune system cells in a homogeneous system, ignoring the spatial description of a tumor [7,8]. The spatio-temporal description of the immune system [29] or tumors [30] often involves partial differential equations for macroscopic quantities, with nevertheless an emphasis on the importance of multiscale modeling [31]. In this work we introduce a space-velocity description of tumor growth at cell scale in the framework of thermostatted kinetic theory. The homogeneous model includes autocatalytic cell production and increased cell activity during favorable interactions. In chemical kinetics, autocatalysis refers to the formation of a product which is also a reactant. Hence, the production rate of the autocatalytically formed species depends on its own concentration, leading to long induction times, when the species is present in low concentration at the beginning of the reaction, and then to explosive behaviors, when the concentration exceeds a threshold. Autocatalysis is a crucial ingredient in the destabilization of a steady state and the emergence of more complex behaviors such as oscillations [32,33]. In the model, the fluctuations of cell activity are regulated by a so-called thermostat [2]. The introduced interaction processes only consider the death of a cell when its activity is not favorable. However, all cells [34], including memory T cells involved in the long-lasting protection acquired after vaccination [35,36,37,38], have a finite life span, whether or not they are very active. Thus, the spontaneous death of educated cells leads to a loss of information. In order to take this phenomenon into account, the model introduces the regulation of activity fluctuations, by imposing that its second moment remains constant. By analogy with the dissipation of energy in a system maintained at constant temperature [1], the activity regulator is also called a thermostat.
We design a kinetic Monte Carlo algorithm inspired by the direct simulation Monte Carlo (DSMC) method initially devoted to dilute gases [39,40] and simulate the evolution of cell position, velocity, nature, and activity. Contrary to the standard DSMC method in which the velocities of the dilute gas particles are updated due to binary collisions, jump processes are introduced to randomize the direction of cell velocities in the present inhomogeneous model involving dense soft matter. Specifically, the algorithm that we develop harnesses the discretization of space into boxes and the implementation of interactions only between cells randomly chosen inside a same spatial box, as suggested by standard DSMC. In spite of this powerful procedure, which earns the name of Monte Carlo to the algorithm, the simulations remain expensive in computation time, because of the large number of cancer cells often created in the case of cancer escape. Our aim is to check if the inhomogeneous model contains enough biological ingredients to reproduce the complex antagonistic effects observed in cancer evolution, including the "three Es" (elimination, equilibrium, and escape of cancer from immunosurveillance) [19,26] and spontaneous oscillations in the number of immune system cells observed during tumor growth [21,22] but not included in the homogeneous model [7,8]. In addition, the simulations of an inhomogeneous system reported in the present paper give access to the evolution of the geometry of the tumor, its possible deformation, shifting, and splitting that were clearly not affordable in the framework of a homogeneous description.
The paper is organized as follows. In section 2, we present the inhomogeneous model and the kinetic equations governing the evolution of the probability distributions of cancer, immune system and normal cells. The kinetic Monte Carlo algorithm is made precise. Simulation results for an initially localized tumor are presented in section 3. Specifically, we examine if, for a poorly effective thermalization, the number of cancer cells rapidly drops, then slightly increases during a long quasi-steady regime, and eventually explodes, although the autocatalytic production of immune system cells is much more favored than cancer cell formation. The existence of oscillations for the total number of immune system cells is also discussed. A sensitivity analysis is performed and the different behaviors observed as the field that regulates activity and the rate constant related to the production of immune system cells vary are shown. The variations of the pseudo-period of oscillations with the field and the rate constant are determined. Section 4 contains conclusions.
The competition between immune system cells and cancer cells can be compared to the interactions between immune system cells and viruses. In both cases, the foreign cells induce responses from the immune system [10,11,12]. Different types of cells are involved in the defense mechanism [41,42]. Dentritic cells ingest cancer cells, isolate antigens and present them to T cells, leading to their activation and proliferation [13]. Killer T cells are able to destroy cancer cells and regulatory T cells control the number of lymphocytes and prevent their overproduction [18,43]. In parallel with the learning of the immune system cells, the learning of the cancer cells takes place, which may lead to cancer escape. Through cell interactions, the cancer cells learn not to be detected, misleading the regulatory T cells which in turn limit the production of immune system cells and let cancer cells proliferate [16,17].
Recently, we have introduced a homogeneous model of competition between cancer cells and immune system cells which reproduces cell interactions, activation, and learning processes [7,8] and accounts for the three Es of immunotherapy. We propose to develop a space-velocity description of tumor growth, which requires to model the updating of cell position x and velocity v in addition to their nature j and activity u. As in the homogeneous model, we introduce normal cells n, cancer cells c, and immune system cells i as a whole with the aim of introducing as few variables and parameters as possible and providing qualitative interpretations of clinical observations.
Inhomogeneous initial conditions with tumor cells located in the middle of the system are preferably considered. A Gaussian distribution Pu(u) of mean μ and variance σ is used to generate the initial activities of the cells, regardless of their nature. Possible heterogeneity in initial activity distribution may be considered. All the cells are supposed to have the same speed which remains constant. The direction of cell velocity v is defined by an angle θv∈]0,2π] randomly chosen using a uniform distribution Pv(θv).
Cell activity u and nature j=c,i,n, where c stands for cancer, i for immune, and n for normal, may change during local interactions that are assumed to leave position x and velocity v of cancer cells and immune system cells unchanged. We consider only three processes which include interaction, activation, proliferation or death, in a similar manner as autocatalytic chemical processes
{c(v,u)+n(v′,u′)κcn(v,u,v′,u′)→c(v,u+ε)+c(v′,u′)R⟶n(v′′,u′′) | (2.1) |
i(v,u)+c(v′,u′)κic(v,u,v′,u′)→i(v,u+ε)+i(v′,u′) | (2.2) |
c(v,u)+i(v′,u′)κci(v,u,v′,u′)→c(v,u+ε)+c(v′,u′) | (2.3) |
where the rate constants κcn(v,u,v′,u′), κic(v,u,v′,u′), and κci(v,u,v′,u′) have the following nontrivial dependence on the activities and the velocities of the interacting couple:
κcj(v,u,v′,u′)=kcj(u−u′)H(u−u′)H(θinter−θ(v,v′))forj=n,i | (2.4) |
κic(v,u,v′,u′)=kic(u−u′)H(u−u′)H(θinter−θ(v,v′)) | (2.5) |
where kcn, kic and kci are constant, θinter∈]0,π] is the maximum angle between the velocity directions permitting interaction and θ(v,v′) is defined as follows:
θ(v,v′)={∣θv−θv′∣if∣θv−θv′∣≤π∣∣θv−θv′∣−2π∣if∣θv−θv′∣>π | (2.6) |
where θv∈]0,2π] is the angle defining the direction of v. The condition on the directions of the velocities of the interacting pair imposed by the Heaviside function H(θinter−θ(v,v′)) ensures that the two cells remain close during a sufficient time for activation and learning to take place.
The learning process is reproduced by the small increase ε of the activity of the cell with the already larger activity before the interaction. The process described in Eq. (2.1) has the form of an autocatalytic production of cancer cells, the presence of a cancer cell being a pre-requisite for the formation of another cancer cell, as in a chemical process catalysed by a product of the reaction.
Equation (2.1) accounts for both the mutation of a normal cell into a cancer cell and the division of a cancer cell. The specific rate constant expression ensures that proliferation of cancer cells is more favorable when they have already learned to blend into their environment, i.e. have increased their activity u. The reservoir R of normal cells maintains their number constant and plays also the role of reservoir of activity. The notion of reservoir is useful in far-from-equilibrium systems able to self-organize. The formation of complex structures requires that the system is maintained far from equilibrium, which occurs through the consumption of energy or matter. A reservoir is an ideal device which is supposed to be able to maintain the concentration of a given species constant by instantaneously absorbing incoming fluxes or injecting the species to compensate for losses [32]. In the present model, the condition on the activities of the interacting cells imposed by the Heaviside function H(u−u′) introduces a selection. The normal cell which disappears has a smaller activity than the interacting cancer cell, which could induce a non physical increase in the mean activity of the remaining normal cells. The model provides for reinjecting a normal cell with a randomly chosen activity u″ using the initial Gaussian distribution Pu(u) of fixed mean value μ, so that the mean activity of normal cells remains constant, consistent with biology. The velocity v″ of the reinjected normal cell is randomly chosen using the initial uniform distribution Pv(θv).
The process given in Eq. (2.2) corresponds to the autocatalytic production of immune system cells and accounts for both the death of a cancer cell attacked by a killer T cell and the proliferation of immune system cells with increased activity, i.e. the knowledge of the presence of cancer cells in their environment. As in Eq. (2.1), the rate constant depends on a Heaviside step function H(u−u′) implying that the increase of the number of immune cells is locally more probable when they have a larger activity than the neighboring cancer cells.
The regulation of the number of immune cells by T-regulatory cells and the learning process of cancer cells are taken into account in the third process given in Eq. (2.3). This process accounts for both the death of an immune system cell and the division of a cancer cell. The rate constant being proportional to H(u−u′), the proliferation of cancer cells locally occurs if their activity is larger than the one of the surrounding immune system cells. Hence, the model accounts for the ability of trained cancer cells to mislead the immune system and trigger the detrimental response of T-regulatory cells.
In order to introduce a regulation of the explosive cell production due to autocatalytic processes, we impose that the rate of each process is proportional to the difference of activities of the interacting pair. Specifically, according to Eq. (2.3), further growth of a cluster of cancer cells with a smaller activity than the surrounding immune system cells is not likely to be observed. Following Eq. (2.2), an analogous phenomenon is expected for immune system cells: A local increase in the number of immune system cells will not be sustained if their activity is smaller that the one of the neighboring cancer cells. According to Eqs. (2.4) and (2.5), the expression of the rate constants, proportional to the relative activity of the interacting cells introduces a remote analogy to the collision term of the Boltzmann equation, which depends on the relative velocity of the encounters [44].
In order to mimic the preferential increase in the number of cancer cells due to the deleterious response of T-regulatory cells rather than by mutation of normal cells, we assign values of kci in Eq. (2.3) larger than the values of kcn in Eq. (2.1).
The processes given in Eqs. (2.1-2.3) associate learning process with systematic increased cell activity. However, cell interactions that are not explicitly taken into account in the scheme make learning less effective. In the model the regulation of the activity fluctuations is ensured by a thermostat mimicking information loss due to dissipation. Thermalization, characterized by a field E and a friction coefficient α, acts on the activity uk of each cell k according to [1,5]
dukdt=E−αuk | (2.7) |
Choosing
α=⟨u⟩E⟨u2⟩ | (2.8) |
where ⟨u⟩ is the mean activity of the entire system, ensures that the second moment ⟨u2⟩ remains constant.
All cells are supposed to move at constant speed in all directions. We introduce a stochastic turning operator that modifies the direction of cell velocities inside a cone of apex angle θjump with 0<θjump≤π/2. Hence, the rate constant associated with the velocity jumps obeys
κv(v,v′)=kvH(θjump−∣θv−θv′∣) | (2.9) |
where kv is constant.
The distribution function fj for cells of type j=c,i,n obeys the following kinetic equation
(∂t+v⋅∇x)fj(t,x,v,u)+∂u((E−αu)fj(t,x,v,u))=Ij+Vj | (2.10) |
where v⋅∇x is the standard advection operator, ∂u((E−αu)fj) is the thermalization operator, Ij is the interaction operator and Vj is the velocity randomization operator. Following Eqs. (2.1-2.3), the interaction operator Ic associated with cancer cells reads
Ic=∫R3∫R+κcn(v,u−ε,v′,u′)fc(t,x,v,u−ε)fn(t,x,v′,u′)dv′du′+∫R3∫R+κcn(v′,u′,v,u)fc(t,x,v′,u′)fn(t,x,v,u)dv′du′−∫R3∫R+κcn(v,u,v′,u′)fc(t,x,v,u)fn(t,x,v′,u′)dv′du′−∫R3∫R+κic(v′,u′,v,u)fi(t,x,v′,u′)fc(t,x,v,u)dv′du′+∫R3∫R+κci(v,u−ε,v′,u′)fc(t,x,v,u−ε)fi(t,x,v′,u′)dv′du′+∫R3∫R+κci(v′,u′,v,u)fc(t,x,v′,u′)fi(t,x,v,u)dv′du′−∫R3∫R+κci(v,u,v′,u′)fc(t,x,v,u)fi(t,x,v′,u′)dv′du′ | (2.11) |
with κcj(v,u,v′,u′) and κic(v,u,v′,u′) given in Eqs. (2.4) and (2.5). The two first terms of the right-hand side account for the positive contributions to the distribution function of cancer cells of velocity v and activity u after the interaction of a cancer cell and a normal cell through the process given in Eq. (2.1). Specifically, the first term corresponds to the formation of a cancer cell with the desired properties from a cancer cell of initial activity u−ε after update by the increment ε. All contributions are obtained by integrating over all acceptable velocities v′ and activities u′ of the interacting normal cell. The second term is obtained when assigning the velocity v and activity u of the disappearing normal cell to a newly formed cancer cell and integrating over all acceptable velocities v′ and activities u′ of the interacting cancer cell. The fourth term is related to the process given in Eq. (2.2) and the fifth to seventh terms are associated with Eq. (2.3). The interaction operator Ii for immune system cells is formed following analogous rules as for cancer cells:
Ii=∫R3∫R+κic(v,u−ε,v′,u′)fi(t,x,v,u−ε)fc(t,x,v′,u′)dv′du′+∫R3∫R+κic(v′,u′,v,u)fi(t,x,v′,u′)fc(t,x,v,u)dv′du′−∫R3∫R+κic(v,u,v′,u′)fi(t,x,v,u)fc(t,x,v′,u′)dv′du′−∫R3∫R+κci(v′,u′,v,u)fc(t,x,v′,u′)fi(t,x,v,u)dv′du′ | (2.12) |
For normal cells, the interaction operator In reads
In=−∫R3∫R+κcn(v′,u′,v,u)fc(t,x,v′,u′)fn(t,x,v,u)dv′du′+Pu(u)Pv(θv)∬ | (2.13) |
The second term of the right-hand side is related to the action of the reservoir, which injects a normal cell of velocity {\bf v} and activity u exactly at the same rate as a normal cell of velocity {\bf v''} and activity u'' just disappeared through Eq. (2.1). The contribution to the evolution of f_n(t, {\bf x}, {\bf v}, u) includes integration over all velocities {\bf v''} and weighting by P_{\bf v}({\theta_v}) , integration over all activities u'' and weighting by P_u(u) , as well as integration over all acceptable velocities {\bf v'} and activities u' of the interacting cancer cell.
Finally, the velocity randomization operator V_j which models the velocity-jump process inside a cone reads
\begin{eqnarray} V_j & = & \int_{{\bf{R}}^3} \kappa_v({\bf v}, {\bf v'}) \left( f_j(t, {\bf x}, {\bf v'}, u)-f_j(t, {\bf x}, {\bf v}, u) \right){\rm d}{\bf v'} \end{eqnarray} | (2.14) |
where the rate constant \kappa_v({\bf v}, {\bf v'}) of clockwise or anticlockwise jump of velocity direction is given in Eq. (2.9).
Graeme Bird [39,40] introduced an efficient algorithm, the direct simulation Monte Carlo (DSMC) method, to simulate the Boltzmann equations associated with a dilute gas [45] and successfully adapted to concentrated solutions [46]. The performance of DSMC with respect to molecular dynamics is related to considering collisions only for particles belonging to a same spatial box and following an acceptance-rejection technique to treat the collisions. Our aim is to show that these hypotheses can be satisfactorily extended to cellular dynamics in order to directly solve the three kinetic equations given in Eq. (2.10) for the different cell types j = c, i, n . We propose the following algorithm. During a time step \Delta t , cell natures, activities, positions and velocities are updated according to the model of interactions, thermalization and velocity jumps described in the previous section. Specifically, for given values of the constants k_{cn} , k_{ic} , and k_{ci} , the maximum number of interactions between pairs is determined in each spatial box and the corresponding interactions between randomly chosen pairs are actually accepted provided that the conditions on activities and velocities expressed in Eqs. (2.4) and (2.5) are met. The nature and the activity of the interacting cells are then updated. After the interactions in all spatial boxes have been considered, the activity of each cell is thermalized according to Eqs. (2.7). Cell positions are modified according to their velocities, which may result in changing of spatial box. Finally, the directions of cell velocities are updated following the stochastic deflection process of rate constant given in Eq. (2.9).
In the next section we present the simulation results obtained for a two-dimensional inhomogeneous system.
Two-dimensional simulations are performed in a square of side l_x = l_y = 30 divided in boxes of side \Delta x = \Delta y = 1 . Periodic boundary conditions are applied. Initially, 20 normal cells and 20 immune system cells per box are homogeneously spread in the entire system. In order to study the growth of an initially localized tumor, we start from 20 cancer cells per box in a central discretized disc of diameter D = 5 and no cancer cells elsewhere. The effects of the rate constant values on the system behavior have been studied in detail in a homogeneous system [7,8] and we have checked that they are not sensitively different in the inhomogeneous case. The absolute values of the rate constants are arbitrary. We choose to present the results obtained in a boosted immune system such that k_{ic} > k_{ci} . This choice results in the more favorable autocatalytic formation of immune system cells than cancer cells for populations with similar activities. We impose k_{ci} > k_{cn} in order to comply with the much more probable division of cancer cells under the pressure exerted by regulatory T cells through Eq. (2.3) than the mutation of normal cells into cancer cells reproduced by Eq. (2.1). We impose k_{ci}/k_{cn} = 10^3 . We refer to the sensitivity analysis performed in a homogeneous system for a detailed discussion of the effect of the parameters introduced in the model [8]. In order to speed up the simulation, we accept the interaction between two cells without condition on the direction of their velocities, which amounts to choosing \theta_{\rm inter} = \pi . The rate constant for velocity deflection k_v and the maximum angle of velocity jump \theta_{\rm jump} are chosen sufficiently large for the randomization of velocity direction to be efficient. We have k_v = 0.1 and \theta_{\rm jump} = \pi/6 .
Before performing a sensitivity analysis and showing the different behaviors encountered as essential parameters of the model vary, we focus on two nontrivial behaviors reproducing clinical observations, the three Es of immunotherapy [19] and oscillations of the number of immune system cells [21,22]. The same parameter values are chosen to illustrate the two typical behaviors, except the ratio k_{ic}/k_{ci} , equal to 20 for the 3Es and 10 for the oscillations.
Figure 1 shows the evolution of the total number N_c of cancer cells and the total number N_i of immune system cells in the entire system. Due to the action of the reservoir, the number of normal cells is constant by construction in agreement with biological requirements. In the simulations, the total number of normal cells is set at 18000 in a system of 30\times30 spatial boxes. A phenomenon analogous to the three Es (elimination, equilibrium, and escape of cancer) observed in immunotherapy [19] is reproduced by the inhomogeneous model for a critical value of the field E depending on the values of the other parameters. The initial elimination phase observed in a homogeneous system [7,8] is not observed and the number N_c of cancer cells first slightly increases. Then, N_c remains nearly constant during a long pseudo-equilibrium phase. The elimination phase would be recovered for an initial number of cancer cells larger than the value of the plateau reached during the equilibrium phase. Finally, N_c explodes, becomes greater than the number of cells of the immune system, N_i , which in turn fluctuates with a large amplitude. Shortly after a very large fluctuation, N_i cancels out and the cancer escapes from immunosurveillance. In the inhomogeneous system, we rather observe the 2Es, for equilibrium and escape. If one excepts the very beginning of the evolution, the behaviors of the homogeneous and inhomogeneous systems are similar and reproduce a very long induction period before explosion, obtained for sufficiently inefficient thermalization in a small interval of field values.
Cancer cells invade the entire system despite k_{ci} < k_{ic} : Inefficient control of activity fluctuations severely limits the chances of success of immunotherapy. Following the expression of the rate constant \kappa_{ci} given in Eq. (2.4), the local important production of cancer cells is related to the learning of these cells, i.e. to the increase of their activity with respect to the neighboring immune system cells. The learning process implemented in the model consists of a small activity increase through cell interactions and a rate dependence on the relative activity of the interacting pairs. This simple procedure is refined enough to reproduce the complex learning process of cancer cells, which can mislead the immune system and induce a harmful regulatory T cell response that wrongly orders a decrease in the production of immune cells.
The spatial description gives access to the evolution of the clusters of cancer cells and immune system cells. A cluster of cancer cell is defined as a continuous set of spatial boxes with neighboring boxes containing at least 20 cancer cells, having a common side and not only a common corner. The value 20 corresponds to the initial number of cancer cells per box in the small, localized tumor. An analogous definition is used for a cluster of immune system cells. The size of a cluster is given by the number of boxes it contains. We use the Hoshen-Kopelman algorithm to build the clusters [47,48]. Figure 2 gives the evolutions of the averaged number of clusters and the size of the biggest cluster for each cell type for the same conditions as in Figure 1.
During the short transient initial period preceding the pseudo-equilibrium phase, the number of clusters of immune system cells rapidly increases while the size of the biggest cluster decreases from the total system size 30 \times 30 to a small value around 10 , revealing the rapid fragmentation of the initial big cluster. During the same period, the number of clusters of cancer cells slightly fluctuates between 1 and 0 and the size of the biggest cluster instantaneously drops from 21 to 0 , meaning that the number of cancer cells per box rapidly becomes smaller than 20 . During the long pseudo-equilibrium phase, the number of clusters of immune system cells linearly decreases due to two phenomena, the coalescence of growing clusters and the decrease of the number of immune system cells, in agreement with the slight increase of the size of the biggest cluster. During the same phase, the number of clusters and the size of the biggest cluster of cancer cells vanish, as if cancer has been eliminated. Then, the number of clusters of cancer cells exponentially grows implying a global, nonlocal escape of cancer while the number of clusters of immune system cells drops. Eventually, the clusters of cancer cells coalesce and a single cluster invades the entire system while the immune system disappears: Cancer becomes out of control. System size has an impact on the time at which an initially localized tumor invades the entire system but does not affect the qualitative behavior of the system.
The main contribution of the inhomogeneous description is shown in Figure 3. The total numbers N_i of immune system cells and the mean activities, u_c , and u_i , of both types of cells display pseudo-oscillations of quite well-defined period and increasing amplitude. The total number N_c of cancer cells rapidly becomes much larger than the total number N_i of immune system cells and the chosen scale prevents from showing the end of the evolution of N_c . The extrema of the mean activities u_c and u_i are associated with typical spatial configurations as shown in the different snapshots of Figures 3-5. A minimum of u_c is rapidly followed by a maximum of u_i and then a maximum of u_c .
As shown in Figure 3, a minimum of u_c corresponds to a nearly symmetric disc of cancer cells attacked by a small cluster of immune system cells due to a local favorable fluctuation of u_i . Rapidly, the cluster of immune system cells grows and, after a small delay, the mean activity u_i reaches a maximum, as displayed in Figure 6. The cluster of immune system cells has an exotic shape which corresponds to a domain of the central tumor in which the number of cancer cells is depleted. According to Figure 5, the symmetry of the entire system is rapidly restored: due to their larger activities in average, the cancer cells deplete the immune system cells in the central part of the tumor, where they are more abundant. As a result, the immune system cells are pushed outwards the disc, forming an external ring. Over a longer period of time, the ring of immune system cells will dissolve, leading to a minimum of N_i , a nearly homogeneous distribution of immune system cells, and a nearly perfect disc of cancer cells. A new local fluctuation of u_i will appear and trigger the cycle again.
As exemplified by the Brusselator model [32], a chemical model presenting time oscillations contains the autocatalytic production of one species X to the detriment of another species Y as well as a feedback mechanism which regenerates Y and consumes X . In the model of cancer/immune system interaction, the conditions on the relative activity of an interacting pair of cells given in Eqs. (2.2) and (2.3) introduce such a feedback mechanism. In a small, localized cluster containing a high number of immune system cells with high activity, the autocatalytic production of immune cells by the process given in Eq. (2.2) is first favored. Then, the total number of cancer cells N_c decreases. The effect on the average activity u_i is non-trivial, the activity of the already present immune system cell being increased by \varepsilon and the newly formed immune system cell inheriting the lower activity of the interacting cancer cell. Simultaneously, the interaction tends to deplete the cancer cell population of its lowest activity elements. As a result, the number of cancer cells decreases but their mean activity u_i increases. The autocatalytic production of cancer cells is then promoted by the process given in Eq. (2.3), which in turn tends to decrease the total number of immune cells N_i but increase their average activity u_i . The system is again in the situation where the process given in Eq. (2.2) favors the autocatalytic production of immune cells and a new cycle begins. Hence, activity-weighted rate constants for the autocatalytic productions of immune system cells and cancer cells is able to induce an oscillating behavior and reproduce immune regulatory cycles observed in advanced cancer patients and revealed by oscillations of C-Reactive Protein concentration [20,21].
The interaction processes do not only select cells depending on their activity, they also introduce criteria on their velocity direction. Small values of the angle of interaction \theta_{\rm inter} coupled to inefficient velocity randomization, i.e. small values of both the rate constant k_v and the maximum angle of velocity jump \theta_{\rm jump} , could induce a bias in the results for large values of v/\Delta t with respect to box size \Delta x . Specifically, a cluster of autocatalytically formed immune system cells with increased activity would move in a similar direction during a certain time, inducing a persistent anisotropy in cluster shape, which would induce a persistent deformation of tumor shape. Adapted parameter values could be chosen to simulate specific tumor shapes. Conversely, efficient velocity randomization is essential to restore the initial tumor symmetry in the presence of velocity selection during cell interactions. The simulation results shown in Figures 3-5 are given for \theta_{\rm inter} = \pi , i.e. in the absence of selection of the angle of the interacting pair, which prevents the formation of clusters with a selected velocity. Cell speed and time step obey v = 10^{-3} and \Delta t = 1 , respectively. The choice of k_v = 0.1 for the rate constant of velocity deflection leads to about 100 velocity jumps before a cell leaves a box, which consists of an efficient velocity randomization. Hence, a randomly occurring cluster of immune system cells surrounded by cancer cells gives rise to the isotropic propagation of a radial front of immune system cells. Then, the cluster depletes from the center. Indeed, the cancer cells in the boxes occupied by the initial cluster of immune system cells which have not been consumed have an increased activity and are thus able to inhibit the immune system cells, leading to the formation of an external ring of immune system cells with radial velocity, similar to the one observed in Figure 4.
The evolution of the averaged number of clusters and size of the biggest cluster are given in Figure 6 for each cell type and for the parameter values of Figures 3-5 leading to oscillations. The number of clusters of cancer cells slightly increases during the entire duration of the simulation, the biggest cluster presenting a first fluctuating growing period followed by a faster increase accompanied by more regular oscillations. The behavior of the size of the biggest cluster is close to the one of the mean number N_c of cancer cells, which shows that the tumor swarms little, remains localized and grows around the initial cluster of cancer cells. The size of the biggest cluster of cancer cells does not exceed 600 boxes and remains smaller than system size equal to 30 \times 30 cells, making explicit that the oscillating behavior shown in Figures 3-5 is not affected by small size effects and boundary conditions. The dependence of the pseudo-period of oscillations on system size is checked in the next subsection.
The first increase of the number of clusters of immune system cells reveals the fragmentation of the big cluster initially occupying the entire system. Then the number of clusters of immune system cells decreases while the size of the biggest cluster remains constant in average, proving that clusters disappear due to immune system cells destruction and not to coalescence. As shown in the snapshot of Figure 4, the maxima of N_i correspond to the periodic appearance of a small cluster containing a large number of immune system cells in the central disc occupied by the large killing population of cancer cells, in which the immune system cells will be soon consumed before a new cycle.
Qualitatively different behaviors are observed as some relevant parameters of the model vary. In a homogeneous system, the field E controlling activity fluctuations and one of the rate constants, k_{ic} or k_{ci} , related to immune system cell or cancer cell production proved to have a significant impact on the phenomena [8]. In the case of an inhomogeneous system, we perform systematic analyses for increasing values of either E or k_{ic} with the other parameters set to the values shown in the caption of Figure 1.
Figures 7 and 8 illustrate typical evolutions of the numbers of cells and mean activities, respectively, as the field E varies for k_{ic} = 0.05 , chosen smaller than k_{ci} . The same vertical scale is chosen for all subfigures of a given figure. For a sufficiently small value of the field, typically E \lesssim 10^{-5} , corresponding to an inefficient control of activity fluctuations, a fast monotonous convergence toward elimination of the immune system and cancer growth is obtained. As E increases, the evolutions of both cell numbers and activities become oscillating, in a chaotic manner for E \simeq 5 \times 10^{-5} and a more regular way illustrated in Figures 3-5 for E \simeq 10^{-4} . For a critical value of the field, E \simeq 1.4 \times 10^{-4} , associated with a bifurcation, irregular fluctuations appear after a long induction period and the behavior illustrates the 3Es, with the final escape of cancer from immunosurveillance. For a slightly larger value of the field, E \simeq 1.5 \times 10^{-4} , the behavior is completely different and the system monotonically converges toward a steady state in which the activities of the two cell types are equal: Cancer is controlled, exactly as in the case of a sufficiently well thermalized homogeneous system [8]. The choice k_{ic} > k_{ci} implies that the final steady number of immune system cells is larger than the final steady number of cancer cells.
The results shown in Figures 7 and 8 correspond to simulations stopped at different times such that the immune system cells have disappeared, the number of cells in a spatial box exceeds a threshold set to 3000 as obtained for the critical value E \simeq 1.4 \times 10^{-4} , or when the mean activities of the cancer cells or the immune system cells are equal. The two first cases obtained for E \leq 1.4 \times 10^{-4} correspond to cancer escape. The last case observed for E \ge 1.5 \times 10^{-4} is associated with cancer control and the simulation time is then supposed to mimic the life expectancy of a patient dying of a disease other than cancer. The characteristic times of the evolution are completely different depending on the field values, in particular because the 3Es are associated with a bifurcation, i.e. a divergence of the induction time before explosion [32]. Sufficiently far from the critical field value, the relaxation toward the stationary state is much faster.
Figure 9 represents the variation of the mean pseudo-period of oscillations of the total number of immune system cells N_i versus E in the field range where they are observed. The uncertainty about their determination is minimum for E \simeq 10^{-4} associated with the most regular oscillations of both N_i , u_i , and u_c and increases as the boundaries of the domain of existence get closer. Even for E \simeq 10^{-4} the time between two maxima of N_i decreases over the evolution, leading to rather large error bars for the pseudo-period. The averaged pseudo-period \langle T \rangle slightly decreases as E increases. The value of \langle T \rangle in a system of 30 \times 30 boxes is satisfactorily compared to the result obtained in a system of 20 \times 20 boxes. The difference between the two results is much smaller than the variability of the pseudo-period during time, proving that the simulation results deduced from a 30 \times 30 system or even a 20 \times 20 system are not sensitive to boundary effects and reliable.
The different typical behaviors encountered as the rate constant k_{ic} controlling the autocatalytic formation of immune system cells increases are represented in Figures 10 and 11, all other parameters being fixed.
As shown in Figure 11 for sufficiently small values of k_{ic} obeying k_{ic} \lesssim 0.007 , the system rapidly and monotonously converges to a steady solution with the same final value u_c^{\rm end} = u_i^{\rm end} of the mean activities of the two cell types c and i . The final cell numbers obey N_c^{\rm end} > N_i^{\rm end} for k_{ic} < k_{ci} , they coincide for k_{ic} = k_{ci} , and obey N_c^{\rm end} < N_i^{\rm end} for k_{ci} < k_{ic} \lesssim 0.007 , as shown in Figure 10. The symmetric roles played by the constants k_{ic} and k_{ci} clearly appear. In all these cases, cancer is controlled. On the contrary, for sufficiently large k_{ic} , such that k_{ic} \gtrsim 0.01 , cancer escapes from the surveillance of the immune system. The case obtained for k_{ic} = 0.01 , already shown in Figure 1, illustrates the 3Es, with a long induction period and the eventual explosion of cancer cell number and both activities. For k_{ic} \simeq 0.05 , rather regular oscillations of cell numbers and activities are obtained. Before disappearing for larger values of k_{ic} , the oscillations become erratic for k_{ic} \simeq 0.25 . The value k_{ic} = 0.01 can be considered as critical, smaller values leading to monotonous cancer control and similar mean activities of both cell types, larger values to oscillating cancer escape and larger mean activity for cancer cells. The results presented in Figures 10 and 11 are nonintuitive since increasing k_{ic} , which should reinforce the production of immune system cells, induces the exit from the domain where cancer is controlled. This result is explained by the feedback effect associated with the activity-dependence of the rate constant \kappa_{ic} given in Eq. (2.5). At constant field E , increasing k_{ic} locally depletes the population of cancer cells from its most active representatives, which induces a local increase of cancer cell activity. This fluctuation of cancer activity is then responsible for cancer cell production through the autocatalytic process given in Eq. (2.3), which explains the eventual explosion of the number of cancer cells and cancer escape from immunosurveillance. In other words, a larger value of k_{ic} leads to a larger critical value E_c of the field associated with the 3Es and delimiting the domains between cancer escape and cancer control: As shown in Figures 7 and 8, for k_{ic} = 0.05 we find E_c\simeq 1.4 \times 10^{-4} and in Figures 10 and 11, for k_{ic} = 0.01 the bifurcation arises for a smaller critical field value E_c\simeq 10^{-4} . A stronger stimulation of the immune system requires a more efficient regulation of activity fluctuations for cancer to be controlled.
Figure 12 shows the variation of the mean pseudo-period of oscillations \langle T \rangle versus the rate constant k_{ic} in the range in which they exist. The uncertainty is much smaller than the variation of the averaged pseudo-period and we conclude without ambiguity that \langle T \rangle increases as k_{ic} increases. It was not easy to anticipate such a result, due, once again, to the antagonistic effects included in the model. At first glance, increasing k_{ic} should lead to the faster production of immune system cells according to the autocatalytic process given in Eq. (2.2), which should decrease the time between two maxima of N_i or u_i . However, as shown in Figures 4 and 5, the formation of a cluster of immune system cells associated with a maximum of mean activity u_i is rapidly followed by a maximum for the mean activity u_c of the cancer cells, due to the depletion of small activity cancer cells induced by the process given in Eq. (2.2). Then, the very active cancer cells consume the immune system cells according to the autocatalytic process given in Eq. (2.3). As a result, the feedback effect on the activities is stronger than the increase of k_{ic} and the time between two maxima of N_i or u_i is longer for larger k_{ic} .
We end the sensitivity analysis of the system to the field E and the rate constant k_{ic} with a few remarks about the effect of the other parameters. As explicitly shown in the case of a homogeneous system [8], the role of k_{ci} is symmetrical to the one of k_{ic} . The results displayed in Figures 10 and 11 for k_{ic} < k_{ci} , k_{ic} = k_{ci} , and k_{ic} > k_{ci} illustrate this property. A small change of cell speed v modifies the critical values of the field and the constant k_{ic} for which the 3Es are observed but does not qualitatively alter the behavior of the system. However, another issue, far beyond the scope of this article, could be addressed in a different domain of parameter space. The model could be extended to the study of the anisotropic growth of a tumor by choosing smaller interaction angles \theta_{\rm inter} , larger cell speeds v , smaller rate constants k_v of velocity randomisation and smaller deflection angles \theta_{\rm jump} in order to favor the formation of clusters of cancer cells with a high activity moving in the same direction during a sufficiently long time to distort the tumor.
The interactions between cancer and the immune system are complex and nontrivial antagonistic effects have been observed in patients with long-standing tumors and revealed by oscillations of some markers of inflammation [21]. We extend a model of cancer and immune system competition [7,8] to an inhomogeneous system in order to follow the dynamics of tumor growth at cell scale within the framework of thermostatted kinetic theory. The model includes cell-cell interactions mimicking mutation, death, division, regulation, and learning associated with increased activity. Activity fluctuations are thermostatted to account for dissipation of information through nonspecific cell interactions. We implement a kinetic Monte Carlo algorithm to directly simulate the kinetic equations in the spirit of DSMC method [39]. The impact of thermalization on the behavior of the system is important and a bifurcation occurs as the field controlling activity fluctuations increases. For small field values, cancer invades the entire system whereas, for large field values, the system reaches a steady state with non vanishing numbers of cancer cells and immune system cells. Cancer is then controlled and the tumor reaches a stationary size. For a critical value of the field for which large activity fluctuations are eventually generated, the model reproduces the three Es of immunoediting [19], leading to the unexpected escape of cancer from immune system control after the elimination and equilibrium phases. For a field value slightly smaller than the critical value, the model is sufficiently complex to reproduce an observed feature of tumor growth associated with rather regular oscillations of the immune response.
We did not observe time oscillations in simulations of a homogeneous system [7,8]. Their stimulation requires non trivial initial conditions that are spontaneously generated in the simulations of the inhomogeneous system. In particular oscillations may begin when a cluster of immune system cells have a sufficiently higher activity than the neighboring cancer cells. The simulations have been performed for an efficient cell velocity randomization which ensures the rapid recovery of the initial symmetry of the tumor. Following the evolution of the mean number of clusters as well as the size of the biggest cluster reveals that the tumor remains a connected space apart from the loss of very small clusters on its periphery. The model involves two autocatalytic reactions, the one producing cancer cells, the other, immune system cells, with a regulation through activity-dependent rate constants. In the present case, large-activity immune system cells kill low-activity cancer cells, leading to increased activity for the remaining cancer cells. These educated cancer cells are then able to kill low-activity immune system cells. Consequently, the mean activity of the remaining immune system cells increases and a new cycle may begin.
It is worth noting that decreasing the rate constant k_{ic} controlling the autocatalytic production of immune system cells has a similar effect on the behavior of the system as increasing the field E . The same succession of bifurcations from cancer proliferation associated with chaotic oscillations, periodic oscillations, long induction related to the 3Es, and cancer control is observed as shown when comparing Figures 7 and 10. Counterintuitively, larger rates of autocatalytic production of immune system cells require more efficient thermalization for cancer to be controlled, due to the feedback effect induced by activity-dependent rate constants.
To our knowledge, the results of the model we developed are the first to account for the 3Es of immunotherapy and oscillations of immune system response within the framework of thermostatted kinetic theory and, in particular, at cell scale. As a perspective, macroscopic equations for concentrations and activities could be derived from the kinetic equations in the limit of appropriate space and time scaling [5,6]. This derivation could give a microscopic interpretation at the scale of cell interactions to mean-field models [23,24,25,26] and to the delays [27,28] introduced in the equations for macroscopic quantities such as the concentrations of different cell types.
The results suggest improved treatment protocols, combining vaccination and chemotherapy at proper times. Reducing the dominance of regulatory T cells over killer T cells may be crucial in immunotherapy: Whereas killer T cells destroy cancer cells, regulatory T cells are susceptible to improperly reduce the number of immune system cells, due to misleading information conveyed by cancer cells. Although the model introduces a single type of immune system cells, it takes into account the role of killer T cells through the process (2.2) and the role of regulatory T cells through the process (2.3). Vaccination can be used to maximize killer T cell response and chemotherapy or radiotherapy, to deplete, inter alia, regulatory T cell population. Hence, the choice of the timing of vaccination, on the one hand, and chemotherapy, on the other hand, with respect to the immune regulatory cycle revealed for example by C-Reactive Protein oscillations, could be crucial in the treatment of some tumors. The model, which reproduces intrinsic oscillations could be used to minimize the averaged number of cancer cells after periodic series of immune system boosting followed by immune system cell depletion at different times of the cycle. The results could give some hints on how a tumor reacts to periodic double disturbances mimicking vaccination and chemotherapy.
All authors declare no conflicts of interest in this paper.
[1] | Sharaya O, Vodolazskaya N (2019) Technological aspects of modification of surface layer of agricultural machines parts. Inn Agric Compl Probl Perspect 3: 82–92 (in Russian). Available from: https://www.bsaa.edu.ru/InfResource/library/Journal3(23)2019.pdf. |
[2] | Aulov V, Rozhkov Y (2016) On the issue of combining electrospark and thermodiffusion methods for hardening the machine parts surfaces in the agribusiness. Mach Equip Rural Area 2: 28–32 (in Russian). Available from: https://doi.org/10.33267/2072-9642-2022-2-28-32. |
[3] | Zobnev V, Markov A, Ivanov S, et al. (2014) Wear resistance of multicomponent diffusion boride coatings on working organs of agricultural machines. Mater Sci Mach Build 1: 435–439 (in Russian). Available from: https://elibrary.ru/download/elibrary_22610860_80996072.pdf. |
[4] |
Martínez-Vázquez J, Rodríguez-Ortiz G, Hortelano-Capetillo J, et al. (2021) Effect of induction heating on Vickers and Knoop hardness of 1045 steel heat treated. J Mech Eng 5: 8–15. https://doi.org/10.35429/JME.2021.15.5.8.15 doi: 10.35429/JME.2021.15.5.8.15
![]() |
[5] |
Dudnikov I (2011) Ensure the safety properties of parts that define security of agricultural machinery. Tech Audit Prod Res 1: 33–36 (in Russian). https://doi.org/10.15587/2312-8372.2011.4853 doi: 10.15587/2312-8372.2011.4853
![]() |
[6] |
Morshed-Behbahani K, Farhat Z, Nasiri A (2024) Effect of surface nanocrystallization on wear behavior of steels: A review. Materials 17: 1618. https://doi.org/10.3390/ma17071618 doi: 10.3390/ma17071618
![]() |
[7] |
Pour-Ali S, Kiani-Rashid A, Babakhani A, et al. (2018) Severe shot peening of AISI 321 with 1000% and 1300% coverages: A comparative study on the surface nanocrystallization, phase transformation, sub-surface microcracks, and microhardness. Int J Mater Res 109: 451–459. https://doi.org/10.3139/146.111622 doi: 10.3139/146.111622
![]() |
[8] |
Roy R, Ghosh S, Kaisar T, et al. (2022) Multi-response optimization of surface grinding process parameters of AISI 4140 alloy steel using response surface methodology and desirability function under dry and wet conditions. Coatings 12: 104. https://doi.org/10.3390/coatings12010104 doi: 10.3390/coatings12010104
![]() |
[9] | Stepanova T (2009) Technologies of Surface Hardening of Machine Parts, Ivanovo: Ivanovo State Chemical-Technological University, 64p (in Russian). Available from: https://www.isuct.ru/sites/default/files/department/ightu/ktmio/08.pdf. |
[10] | Ostromenskiy P, Aksenov V, Korotaev B, et al. (2001) Prospects for the use of high-energy technologies to increase the lateral wear resistance of rails. Curr Probl Transp Asian Part Russ 2001: 92–98 (in Russian). Available from: https://is.gd/dUCNDS. |
[11] | Maisuradze M (2022) Induction and Laser Thermal Treatment of Steel Products: A Textbook, Yekaterinburg: Ural University Publishing, 96p (in Russian). Available from: https://elar.urfu.ru/bitstream/10995/117127/1/978-5-7996-3544-2_2022.pdf. |
[12] |
Khisamutdinov R, Zvezdin V, Israfilov I, et al. (2016) Study of processes of steels surfaces modification with highly concentrated energy flows. J Phys Conf Ser 669: 012024. https://doi.org/10.1088/1742-6596/669/1/012024 doi: 10.1088/1742-6596/669/1/012024
![]() |
[13] | Barmin A, Rizakhanov R, Rudstein R (2012) Optimization of quenching regimes for carbon steels by electron beam treatment. Proc Interuniv Sci Sch Young Spec 1: 62–67 (in Russian). Available from: http://nuclphys.sinp.msu.ru/school/s12/12_15.pdf. |
[14] | Belinin D, Shchitsyn Y (2012) Features of structurization at plasma surface hardening on big depth of products from 40Cr13. Proc Samara Sci Cent Russ Acad Sci 4: 1200–1205 (in Russian). Available from: https://elar.urfu.ru/bitstream/10995/30926/1/sid_2014_04.pdf. |
[15] | Belinin D, Verkhorubov V, Kuchev P, et al. (2011) Plasma surface hardening of hard loading constructions made of steel 40Kh13. Bulletin Pnrpu 2: 12–18 (in Russian). Available from: https://cyberleninka.ru/article/n/plazmennaya-zakalka-tyazhelonagruzhennyh-detaley-iz-stali-40h13. |
[16] | Sidorov S (1998) Technical level and resource of working organs of agricultural machinery. Tract Agric Mach 3: 29 (in Russian). Available from: https://rusneb.ru/catalog/000200_000018_RU_NLR_PER_B_2381825_1998_3/. |
[17] | Kravchenko N (2013) Plasma Methods of Hardening and Restoration of Working Bodies Road-Building and Soil-Cultivating Machines, Moscow: Eco-Press, 328 p (in Russian). Available from: https://expo-books.ru/category/book?id = 12005. |
[18] | Kanayev A (2015) Modernization of the surface layer structure of structural steel. Herald Scin Seifullin Kaz Agro Tehn Univ 3: 78–86 (in Russian). Available from: https://bulletinofscience.kazatu.edu.kz/index.php/bulletinofscience/article/view/701. |
[19] |
Korotkov V (2011) Wear resistance of plasma-hardened materials. J Frict Wear 32: 17–22. https://doi.org/10.3103/S1068366611010077 doi: 10.3103/S1068366611010077
![]() |
[20] |
Jumbad V, Chel A, Verma U (2020) Application of electrolytic plasma process in surface improvement of metals: A review. Letters Appl NanoBioSci 9: 1249–1262. https://doi.org/10.33263/lianbs93.12491262 doi: 10.33263/lianbs93.12491262
![]() |
[21] |
Gupta P, Tenhundfeld G, Daigle E (2007) Electrolytic plasma technology: Science and engineering—An overview. Surf Coat Tech 201: 8746–8760. https://doi.org/10.1016/j.surfcoat.2006.11.023 doi: 10.1016/j.surfcoat.2006.11.023
![]() |
[22] |
Rakhadilov B, Bayatanova L, Kurbanbekov S, et al. (2023) Investigation on the effect of technological parameters of electrolyte-plasma cementation method on phase structure and mechanical properties of structural steel 20X. AIMS Mater Sci 10: 934–947. https://doi.org/10.3934/matersci.2023050 doi: 10.3934/matersci.2023050
![]() |
[23] |
Meletis E, Nie X, Wang F (2002) Electrolytic plasma processing for cleaning and metal-coating of steel surfaces. Surf Coat Tech 150: 246–256. https://doi.org/10.1016/S0257-8972(01)01521-3 doi: 10.1016/S0257-8972(01)01521-3
![]() |
[24] | Belkin P (2013) Electrolytic-plasma modification of metals and alloys. Bull Kostroma State Univ 5: 5–11 (in Russian). Available from: https://cyberleninka.ru/article/n/elektrolitno-plazmennaya-modifikatsiya-metallov-i-splavov-1. |
[25] |
Jiang Y, Bao Y, Wang M (2017) Kinetic analysis of additive on plasma electrolytic boriding. Coatings 7: 61. https://doi.org/10.3390/coatings7050061 doi: 10.3390/coatings7050061
![]() |
[26] |
Belkin V, Belkin P, Krit B, et al. (2019) Increasing wear resistance of low-carbon steel by anodic plasma-electrolytic nitroboriding. J Mater Eng Perform 29: 564–572. https://doi.org/10.1007/s11665-019-04521-1 doi: 10.1007/s11665-019-04521-1
![]() |
[27] |
Taheri P, Dehghanian C, Aliofkhazraei M, et al. (2007) Nanocrystalline structure produced by complex surface treatments: plasma electrolytic nitrocarburizing, boronitriding, borocarburizing, and borocarbonitriding. Plasma Process Polym 4: 721–727. https://doi.org/10.1002/ppap.200731805 doi: 10.1002/ppap.200731805
![]() |
[28] |
Skakov M, Rakhadilov B, Sheffler M (2013) Influence of electrolyte plasma treatment on structure, phase composition and microhardness of steel Р6М5. Key Eng Mater 531–532: 627–631. https://doi.org/10.4028/www.scientific.net/KEM.531-532.627 doi: 10.4028/www.scientific.net/KEM.531-532.627
![]() |
[29] |
Luk S, Leung T, Miu W (1999) A study of the effect of average preset voltage on effective case depth during electrolytic surface-hardening. Mater Charact 42: 65–71. https://doi.org/10.1016/S1044-5803(98)00044-8 doi: 10.1016/S1044-5803(98)00044-8
![]() |
[30] |
Tarakci M, Korkmaz K, Gencer Y (2005) Plasma electrolytic surface carburizing and hardening of pure iron. Surf Coat Tech 199: 205–212. https://doi.org/10.1016/j.surfcoat.2005.02.117 doi: 10.1016/j.surfcoat.2005.02.117
![]() |
[31] | Satbayeva Z (2022) Structure Formation in Alloyed Steels During Electrolytic-Plasma Surface Hardening, Oskemen: Sarsen Amanzholov East Kazakhstan University, 160p (in Russian). Available from: https://nabrk.kz/ru/e-catalog?catalog=4&language=rus&page=8&sphere=3&topic=0&publication_type=13. |
[32] | Suminov I (2011) Plasma Electrolytic Modification of the Surface of Metals and Alloys, Moscow: Tekhnosfera, 464 p (in Russian). Available from: https://f.eruditor.link/file/2686301/. |
[33] |
Cenk Kumruoğlu L, Özel A (2010) Surface modification of AISI 4140 steel using electrolytic plasma thermocyclic treatment. Mater Manuf Process 25: 923–931. https://doi.org/10.1080/10426911003720839 doi: 10.1080/10426911003720839
![]() |
[34] |
Sagdoldina Z, Zhurerova L, Tyurin Y, et al. (2022) Modification of the surface of 40Kh steel by electrolytic plasma hardening. Metals 12: 2071. https://doi.org/10.3390/met12122071 doi: 10.3390/met12122071
![]() |
[35] | Smirnov M (1999) Fundamentals of Heat Treatment of Steel, Yekaterinburg: Ural Branch of the Russian Academy of Sciences, 494p (in Russian). Available from: http://www.materialscience.ru/shared_folder/matved/books/Smirnov_Osnovy_TO_stali.djvu. |
[36] | Kartonova L (2020) Theory and Technology of Heat Treatment, Vladimir: Publishing house of Vladimir State University, 128p (in Russian). Available from: https://dspace.www1.vlsu.ru/bitstream/123456789/8725/1/02082.pdf. |
[37] |
Rasouli D, Khameneh Asl S, Akbarzadeh A, et al. (2008) Effect of cooling rate on the microstructure and mechanical properties of microalloyed forging steel. J Mater Process Technol 206: 92–98. https://doi.org/10.1016/j.jmatprotec.2007.12.006 doi: 10.1016/j.jmatprotec.2007.12.006
![]() |
[38] |
Sunardi S, Lusiani R, Listijorini E, et al. (2021) The effect of airflow speed as cooling media in the hardening process to the hardness, corrosion rate and fatigue life of medium carbon steel. Mater Sci Forum 1045: 40–49. https://doi.org/10.4028/www.scientific.net/msf.1045.40 doi: 10.4028/www.scientific.net/msf.1045.40
![]() |
[39] |
Jo H, Kang M, Park G, et al. (2020) Effects of cooling rate during quenching and tempering conditions on microstructures and mechanical properties of carbon steel flange. Materials 13: 4186. https://doi.org/10.3390/ma13184186 doi: 10.3390/ma13184186
![]() |
[40] |
Basori I, Pratiwi W, Dwiyati S (2019) Effect of salt quenching on the microstructures and mechanical properties of AISI 1045 steel. J Phys Conf Ser 5: 055102. https://doi.org/10.1088/1742-6596/1402/5/055102 doi: 10.1088/1742-6596/1402/5/055102
![]() |
[41] | Pérez R, Llano M, Ravagli R, et al. (2018) Effect of machining fluid like quenching media on the friction and wear behavior of AISI 1045 steel. Int J Mech Eng Technol 9: 146–154. Available from: https://iaeme.com/MasterAdmin/Journal_uploads/IJMET/VOLUME_9_ISSUE_7/IJMET_09_07_017.pdf. |
[42] |
Vieira E, Biehl L, Medeiros J, et al. (2021) Evaluation of the characteristics of an AISI 1045 steel quenched in different concentration of polymer solutions of polyvinylpyrrolidone. Sci Rep 11: 1313. https://doi.org/10.1038/s41598-020-79060-0 doi: 10.1038/s41598-020-79060-0
![]() |
[43] |
Dayanç A, Karaca B, Kumruoğlu L (2017) The cathodic electrolytic plasma hardening of steel and cast iron based automotive camshafts. Acta Phys Pol A 131: 374–378. https://doi.org/10.12693/APhysPolA.131.374 doi: 10.12693/APhysPolA.131.374
![]() |
[44] |
Belkin P, Kusmanov S (2016) Plasma electrolytic hardening of steels: review. Surf Engin Appl Electrochem 52: 531–546. https://doi.org/10.3103/S106837551606003X doi: 10.3103/S106837551606003X
![]() |
[45] |
Bayati M, Molaei R, Janghorban K (2011) Surface alloying of carbon steels from electrolytic plasma. Met Sci Heat Treat 53: 91–94. https://doi.org/10.1007/s11041-011-9347-5 doi: 10.1007/s11041-011-9347-5
![]() |
[46] | Rakhadilov B, Satbayeva Z, Bayatanova L, et al. (2019) Influence of electrolyte-plasma surface hardening on the structure and properties of steel 40KhN. J Phys Conf Ser 1393: 012119. https://doi.org/101088/1742-6596/1393/1/012119 |
[47] |
Ayday A, Durman M (2013) Surface hardening of ductile cast iron by electrolytic plasma technology. Acta Phys Pol A 123: 291–293. https://doi.org/10.12693/APhysPolA.123.291 doi: 10.12693/APhysPolA.123.291
![]() |
[48] | Gorelik S (2002) X-ray and Electron-optical Analysis, 4 Eds., Moscow: Publishing house MISIS, 360 p (in Russian). Available from: https://www.geokniga.org/bookfiles/geokniga-rentgenograficheskiy-analiz.pdf. |
[49] |
Rastkar A, Shokri B (2012) Surface modification and wear test of carbon steel by plasma electrolytic nitrocarburizing. Surf Interface Anal 44: 342–351. https://doi.org/10.1002/sia.3808 doi: 10.1002/sia.3808
![]() |
[50] |
Zhang L, Jin Y, Wang X, et al. (2019) Surface alloys of 0.45 C carbon steel produced by high current pulsed electron beam. High Temp Mater Process 38: 444–451. https://doi.org/10.1515/htmp-2018-0065 doi: 10.1515/htmp-2018-0065
![]() |
[51] |
Zheng B, Huang Z, Xing J, et al. (2016) Three-body abrasive behavior of cementite–iron composite with different cementite volume fractions. Tribol Lett 62: 1–11. https://doi.org/10.1007/s11249-016-0683-x doi: 10.1007/s11249-016-0683-x
![]() |
[52] |
Yaghmazadeh M, Dehghanian C (2009) Surface hardening of AISI H13 steel using pulsed plasma electrolytic carburizing (PPEC). Plasma Processes Polym 6: S168–S172. https://doi.org/10.1002/ppap.200930410 doi: 10.1002/ppap.200930410
![]() |
[53] |
Kurbanbekov S, Skakov M, Baklanov V, et al. (2017) Changes in mechanical properties and structure of electrolytic plasma treated×12 CrNi18 10 Ti stainless steel. Mater Test 59: 361–365. https://doi.org/10.3139/120.111014 doi: 10.3139/120.111014
![]() |
[54] |
Skakov M, Zhurerova L, Scheffler M (2013) Influence of regimes electrolytic-plasma processing on phase structure and hardening of steel 30CrMnSi. Adv Mat Res 601: 79–83. https://doi.org/10.4028/www.scientific.net/AMR.601.79 doi: 10.4028/www.scientific.net/AMR.601.79
![]() |
[55] |
Tabieva E, Zhurerova L, Baizhan D (2020) Influence of electrolyte-plasma hardening technological parameters on the structure and properties of banding steel 2. Key Eng Mater 839: 57–62. https://doi.org/10.4028/www.scientific.net/KEM.839.57 doi: 10.4028/www.scientific.net/KEM.839.57
![]() |
[56] |
Bayati M, Molaei R, Janghorban K (2010) Surface modification of AISI 1045 carbon steel by the electrolytic plasma process. Metall Mater Trans A 41: 906–911. https://doi.org/10.1007/s11661-009-0165-y doi: 10.1007/s11661-009-0165-y
![]() |
[57] |
Nunura C, Santos C, Spim J (2015) Numerical–experimental correlation of microstructures, cooling rates and mechanical properties of AISI 1045 steel during the Jominy end-quench test. Mater Des 76: 230–243. https://doi.org/10.1016/j.matdes.2015.03.031 doi: 10.1016/j.matdes.2015.03.031
![]() |
[58] |
Li Q, Xiao K, Lu Z, et al. (2024) Mechanism of laser carbonitriding enhancing the wear resistance of 45# steel. Mater Today Commun 38: 108327. https://doi.org/10.1016/j.mtcomm.2024.108327 doi: 10.1016/j.mtcomm.2024.108327
![]() |
[59] |
Korotkov V (2015) Wear resistance of carbon steel with different types of hardening. J Frict Wear 36: 149–152. https://doi.org/10.3103/S1068366615020105 doi: 10.3103/S1068366615020105
![]() |
[60] |
Algodi S, Salman A, Al-Helli A (2023) Microstructure and mechanical properties of AISI 1106/AISI 1045 steels drawn arc stud welded joints. Adv Sci Technol Res J 17: 48–55. https://doi.org/10.12913/22998624/171020 doi: 10.12913/22998624/171020
![]() |