Model | Surface | Input | Output | ||||
NN-ANN Model-I | (γp=0.0) | (εT) | (Hp) | (Pr) | (βp) | (Rp) | NN |
NN-ANN Model-II | (γp=0.5) | (εT) | (Hp) | (Pr) | (βp) | (Rp) | NN |
Cancer is a complex group of diseases characterized by uncontrolled cell growth that can spread throughout the body, leading to serious health issues. Traditional treatments mainly include chemotherapy, surgery, and radiotherapy. Although combining different therapies is becoming more common, predicting how these treatments will interact and what side effects they may cause, such as gastrointestinal or neurological problems, can be challenging. This research applies optimal control theory (OCT) to create precise and personalized treatment plans for cancer patients. OCT helps identify the most effective doses of chemotherapy and immunotherapy by forecasting how various treatment combinations will impact tumor growth and the immune response over time. It optimizes the integration of chemotherapy with immunotherapy to minimize side effects while maximizing therapeutic benefits. The study proposes a model for managing malignant tumors using a mix of immunotherapy, vaccines, and chemotherapy. The aim is to develop the best treatment plan that reduces new tumor growth while keeping healthy cells stable. It also takes into account individual differences among patients, including variations in tumor biology and immune responses in both younger and older individuals. To do this, we compared different optimal control strategies: interior point optimization (IPOPT), an open-source tool for nonlinear optimization; state-dependent Riccati equation (SDRE), which adapts linear control methods for nonlinear situations; and approximate sequence Riccati equation (ASRE), a globally optimal feedback control approach for nonlinear systems. The optimization criterion showed that the proposed work achieved a cost value of 52.3573 for IPOPT, compared with 52.424 for both SDRE and ASRE. For CD8+ T cells, the proposed method maintained a consistent value of 1.6499 for continuous (C) and dosed (D) across all techniques. Tumor cell counts had a C value of 0.0007 for IPOPT, compared with 0.0006 for ISDRE and ASRE, with D values remaining at 0 across all methods. This comparison demonstrates the successful use of control theory techniques and highlights their potential for developing personalized and effective treatment strategies for complex cancer cases. By optimizing treatment schedules and dosages, OCT can help minimize the side effects of cancer therapies, thereby enhancing patients' overall quality of life.
Citation: Ahmed J. Abougarair, Mohsen Bakouri, Abdulrahman Alduraywish, Omar G. Mrehel, Abdulrahman Alqahtani, Tariq Alqahtani, Yousef Alharbi, Md Samsuzzaman. Optimizing cancer treatment using optimal control theory[J]. AIMS Mathematics, 2024, 9(11): 31740-31769. doi: 10.3934/math.20241526
[1] | Xiaolei Dong, Yuming Qin . Finite fractal dimension of pullback attractors for a nonclassical diffusion equation. AIMS Mathematics, 2022, 7(5): 8064-8079. doi: 10.3934/math.2022449 |
[2] | Li Yang . Pullback random attractors of stochastic strongly damped wave equations with variable delays on unbounded domains. AIMS Mathematics, 2021, 6(12): 13634-13664. doi: 10.3934/math.2021793 |
[3] | Xin Liu . Stability of random attractors for non-autonomous stochastic p-Laplacian lattice equations with random viscosity. AIMS Mathematics, 2025, 10(3): 7396-7413. doi: 10.3934/math.2025339 |
[4] | Xiaobin Yao . Random attractors for non-autonomous stochastic plate equations with multiplicative noise and nonlinear damping. AIMS Mathematics, 2020, 5(3): 2577-2607. doi: 10.3934/math.2020169 |
[5] | Jean De Dieu Mangoubi, Mayeul Evrard Isseret Goyaud, Daniel Moukoko . Pullback attractor for a nonautonomous parabolic Cahn-Hilliard phase-field system. AIMS Mathematics, 2023, 8(9): 22037-22066. doi: 10.3934/math.20231123 |
[6] | Qinghua Zhou, Li Wan, Hongbo Fu, Qunjiao Zhang . Pullback attractor of Hopfield neural networks with multiple time-varying delays. AIMS Mathematics, 2021, 6(7): 7441-7455. doi: 10.3934/math.2021435 |
[7] | Binbin Miao, Chongbin Xu, Caidi Zhao . Statistical solution and piecewise Liouville theorem for the impulsive discrete Zakharov equations. AIMS Mathematics, 2022, 7(5): 9089-9116. doi: 10.3934/math.2022505 |
[8] | Jianlong Wu . Regularity criteria for the 3D generalized Navier-Stokes equations with nonlinear damping term. AIMS Mathematics, 2024, 9(6): 16250-16259. doi: 10.3934/math.2024786 |
[9] | Chunting Ji, Hui Liu, Jie Xin . Random attractors of the stochastic extended Brusselator system with a multiplicative noise. AIMS Mathematics, 2020, 5(4): 3584-3611. doi: 10.3934/math.2020233 |
[10] | Shaoliang Yuan, Lin Cheng, Liangyong Lin . Existence and uniqueness of solutions for the two-dimensional Euler and Navier-Stokes equations with initial data in H1. AIMS Mathematics, 2025, 10(4): 9310-9321. doi: 10.3934/math.2025428 |
Cancer is a complex group of diseases characterized by uncontrolled cell growth that can spread throughout the body, leading to serious health issues. Traditional treatments mainly include chemotherapy, surgery, and radiotherapy. Although combining different therapies is becoming more common, predicting how these treatments will interact and what side effects they may cause, such as gastrointestinal or neurological problems, can be challenging. This research applies optimal control theory (OCT) to create precise and personalized treatment plans for cancer patients. OCT helps identify the most effective doses of chemotherapy and immunotherapy by forecasting how various treatment combinations will impact tumor growth and the immune response over time. It optimizes the integration of chemotherapy with immunotherapy to minimize side effects while maximizing therapeutic benefits. The study proposes a model for managing malignant tumors using a mix of immunotherapy, vaccines, and chemotherapy. The aim is to develop the best treatment plan that reduces new tumor growth while keeping healthy cells stable. It also takes into account individual differences among patients, including variations in tumor biology and immune responses in both younger and older individuals. To do this, we compared different optimal control strategies: interior point optimization (IPOPT), an open-source tool for nonlinear optimization; state-dependent Riccati equation (SDRE), which adapts linear control methods for nonlinear situations; and approximate sequence Riccati equation (ASRE), a globally optimal feedback control approach for nonlinear systems. The optimization criterion showed that the proposed work achieved a cost value of 52.3573 for IPOPT, compared with 52.424 for both SDRE and ASRE. For CD8+ T cells, the proposed method maintained a consistent value of 1.6499 for continuous (C) and dosed (D) across all techniques. Tumor cell counts had a C value of 0.0007 for IPOPT, compared with 0.0006 for ISDRE and ASRE, with D values remaining at 0 across all methods. This comparison demonstrates the successful use of control theory techniques and highlights their potential for developing personalized and effective treatment strategies for complex cancer cases. By optimizing treatment schedules and dosages, OCT can help minimize the side effects of cancer therapies, thereby enhancing patients' overall quality of life.
Non-Newtonian fluids are used in a variety of products, including culinary items, personal protective equipment, printing technologies, damping and braking systems, and drag-reducing agents [1]. Fluids that are not Newtonian exhibit complex rheological behavior [2]. Unlike Newtonian fluids, these materials do not behave linearly between the shear rate or shear and viscosity. Non-Newtonian fluids [3,4] are categorized into many groups based on how they react to changes in shear stress or rate. Based on their behavior, non-Newtonian fluids can be classified as shear thinning [5,6], shear thickening [7,8], thixotropic, visco-plastic, and rheopectic fluids. Following such characteristics, various studies were given by researchers to explore the type of non-Newtonian fluid models and the corresponding daily life applications. In this direction, the Casson fluid model is one of the non-Newtonian fluid models that best capture yield stress and has noteworthy applications in the biomechanics and polymer sectors. Several different goods, including synthetic lubricants, medicinal compounds, paints, coal, china clay, and many others, can be prepared using Casson fluid. Human blood can also be referred to as Casson fluid since it contains a variety of components including human red blood cells, fibrinogen, protein, and globulin in aqueous base plasma. The operative use of Casson fluid in biological treatments, drilling processes, bio-engineering and food processing has captured the attention of researchers like the analysis of diffusion phenomena in stenosed capillary-tissue exchange systems provided by Siddiqui and Mishra [9]. In order to depict the blood flow they utilize a modified Casson's fluid description. Through the alteration of a measure known as the retention parameter, the severity of the condition could be evaluated. Another significant aspect of this analysis was the increase of these variables together with the progression of the stenosis. For the investigation of the normal and sick states, the concentration profiles were identified. In the recent past, according to Abolbashari et al. [10], engineering and industrial applications have found usage for nano-fluids because of their fantastic thermal conductivity increase. Due to their significance, analyze the fluid non-Newtonian flow with nanoparticles towards the sheet. In special circumstances, the current analytical solution shows a very strong correlation with those of the earlier research that was published. The impact of various flow physical factors were thoroughly explored. Bhattacharyya et al. [11] conducted an examination of Casson fluid flow past a shrinking surface with different effects. Self-similar nonlinear ordinary differential equations (ODEs) were obtained. The translated energy and velocity equations both had two exact answers. According to the depicted data, the temperature drops for greater values of radiation parameter, Prandtl number, and power-law exponent. Additionally, the thickness of the thermal boundary layer increases with direct variation in wall temperature. In some instances, temperature overshoot was seen in the temperature field's graphical depiction. Therefore, under some circumstances, heat transfer from the surface happens instead of heat absorption at the surface. In the presence of fluctuating suction, Abd El-Aziz and Yahya [12] investigated the impact of Hall current on the Casson fluid flow over a sheet. To obtain non-dimensional flow fields, the resulting equations were analytically solved. The temperature and velocity profiles are studied for various parameter values that entered the problem and were displayed graphically. Physically relevant Nusselt number was sorted out and displayed as graphics. For an axisymmetric Casson stagnation point, Nawaz et al. [13] looked at the effects of Joule heating and viscous dissipation. The ensuing boundary value issue was resolved using the homotopy analysis method. The homotopy solutions and the numerical outcomes were compared. The behavior of various parameters was examined. The Casson fluid parameter's resultant effect seems to accelerate the fluid's velocity. The actions of dissipation and Joule heating cause the system to warm up. The thermal and velocity slip effects on a Magnetohydrodynamics (MHD) Casson fluid were studied by Usman et al. [14] over a cylinder with Buongiorno's model. The Casson nanofluid model's equations were constructed and reduced by using similarity transformation. The collocation approach was used to arrive at the numerical solution. The graphs that compare the nanofluid fluid flow to the emerging physical parameters provide an analysis of the physical quantities of interest. The effects of Newtonian heating on the Casson fluid flow were examined by Tassaddiq et al. [15]. In this study, porous effects for such fluids and MHD effects were also taken into account. Equations involving partial differentials were used to model the fundamental issue. The analytical tool is known as the Laplace transform was used to achieve the "Velocity" and "Temperature" functions. In order to analyze the modeling parameters that were used, graphical representations were used. To verify the facts, numerical calculations were made. The graphical results show that velocity clearly declines for magnetic parameter is intensified and increases as the porosity parameter is increased (conjugate parameter). For all conceivable Casson parameter values, the fluid flow might be controlled. A comparative examination of multiphase flow over a steep channel was carried out by Hussain et al. [16]. The combination of microscopic gold particles with Casson fluid creates an intuitive non-Newtonian particulate suspension. Consideration of a uniformly slanted conduit allowed for the consideration of gravitational and magnetic influences. The shear-thinning phenomena is further influenced by heating effects at the channel's edge. For two-phase fluid flow with heat transfer, an exaction was found. In addition, Newtonian-gold particulate flow and Casson-gold particulate flow, which was a non-Newtonian suspension of particle flow, are contrasted. It was discovered that the shear-thinning effects of magnetized Newtonian particle suspension across the inclined channel are more pronounced and experience reduced skin friction. However, it was discovered that Casson fluid was a very helpful suspension for the fabrication and coatings processes, among other things. Additionally, the magnetic field affects the mobility of both phases as a resistive force. Due to the system's significant viscous dissipation, more energy was added. Akolade and Tijani [17] discussed the effect of various transport parameters on Casson- nanofluids through a Riga surface. After imposing sufficient similarity factors, the flow equations were derived. To account for all relevant flow characteristics, the spectral quasi-linearization technique (SQLM) was used. Aside from great agreement with previously published work, the approach performance (numerical estimation tools) was examined in terms of CPU time and mistakes. For the current challenge, engineering dimensionless parameters were also supplied. The outcomes demonstrated that the Casson fluid had a lower flow resistance than the Williamson fluid. Furthermore, Casson fluid diffuses and conducts less than Williamson fluid. Siddiqui et al. [18] investigated the Casson nanofluid between two spinning and stretchy disks in this study. The effects of thermal radiation, magnetic field, Darcy Forchheimer flow, and heat source were investigated. The shooting strategy was owned to solve the existing equations with the help of MATLAB packages. The temperature profile when the Prandtl number was varied was also discussed. Furthermore, the thermal profile towards the thermophoresis, radiation, and Brownian parameter was investigated. Graphs were used to evaluate the impact of nanoparticle concentration on microorganisms with rising Peclet and Lewis numbers. References [19,20,21,22] provide an evaluation of the most recent attempts at Casson fluid flow in a variety of configurations.
In this paper, we used the prediction application of artificial intelligence to examine the thermal stagnation point magnetized flow field of non-Newtonian fluid flow due to inclined heated stretched surfaces with thermos-physical effects namely viscous dissipation, mixed convection, heat generation, and variable thermal conductivity. The flow is mathematically formulated for a cylindrical surface and later mathematical reduction is done for thermal flow over a flat surface. The ultimate results are obtained by using the shooting method and are offered as line graphs. The quantity of interest includes the temperature and NN at both flat and cylindrical surfaces. Two different ANN models are constructed to forecast the heat transfer coefficient at surfaces. We are confident that the thermal findings by the use of ANN models will helpful for investigators affiliated with the area of thermal engineering.
Stretching inclined surfaces, such as the cylinder and plate, are used to introduce non-Newtonian fluid flow. By taking into account the magnetic field [23,24,25], mixed convection, stagnation point flow [26], temperature-dependent conductivity [27], viscous dissipation [28], and heat generation [29] the novelty of the flow field is strengthened. By utilizing the concentration equation, the mass transfer is taken into account. It is assumed that the strength of temperature and concentration is considered higher at surfaces in comparison with strength far away from the surface. The concluding Casson fluid flow equations [26] for the given problem are expressed as
∂(ˆrˆu)∂ˆx+∂(ˆrˆv)∂ˆr=0, | (1) |
ˆu∂ˆu∂ˆx+ˆv∂ˆu∂ˆr=ν(1+1β)(∂2ˆu∂ˆr2+1ˆr∂ˆu∂ˆr)+ˆueˆue∂ˆx−σB02ρ(ˆu−ˆue)+g0(ˆT−ˆT∞)BTcos(α)+g0(ˆC−ˆC∞)BCcos(α), | (2) |
ρcρ(ˆu∂ˆT∂ˆx+ˆv∂ˆT∂ˆr)=1ˆr∂∂ˆr(κ∂ˆT∂ˆr)+ˉμ(1+1β)(∂ˆu∂ˆr)2−1ˆr∂∂ˆr(ˆrˆq)+Q0(ˆT−ˆT∞). | (3) |
ˆu∂ˆC∂ˆx+ˆv∂ˆC∂ˆr=Dm∂2ˆC∂ˆr2. | (4) |
The radioactive flux and thermal conductivity relations are given as:
ˆq=−16σ∗T3∞3k∗∂ˆT∂ˆr, |
κ(ˆT)=κ∞(1+εˆT−ˆT∞ˆTw−ˆT∞), | (5) |
The flow conditions are:
ˆu=ˆUw=aˆx,,ˆC=ˆCw,ˆv=0ˆT=ˆTw,at ˆr=Rc,ˆu=ˆue=dˆx,ˆC→ˆC∞,ˆT→ˆT∞,as r→∞. | (6) |
Here, ν is kinematic viscosity, Rc stands for radius of cylinder, BT is the thermal expansion coefficient, BC denote the coefficient of solutal expansion, ρ,ˆT,ˆC,Q0,ε,κ∞ and Dm are fluid density, temperature, concentration, heat generation coefficient, small parameter, heat conductivity far from the surface and mass diffusivity respectively. By using NN-ANN model, we are interested in comparing how flow parameters affect the Nusselt number and temperature of Casson fluid flow toward both surfaces. Consequently, we must convert these partial differential equations (PDEs) into corresponding ODEs in order to solve Eq (1) through (4). For this reason, we have own the following set of variables [26,27]:
η=ˆr2−Rc22Rc√U0νL,ˆu=ˆxU0LF′(η),ˆv=−Rcˆr√νU0LF(η), |
ϕn(η)=ˆC−ˆC∞ˆCw−ˆC∞,θn(η)=ˆT−ˆT∞ˆTw−ˆT∞. | (7) |
Owning Eq (7), Eqs (1)-(3) turns into
(F‴(1+2γpη)+2γpF″)(1+1/βp)+FF″−F′2+GTθncos(α)+GCϕncos(α)−Mp2(F′−Ap)+Ap2=0, | (8) |
(θ″n(1+2γpη)+2γpθ′n)(1+43Rp)+εT((θn′2+θnθ″n)(1+2ηγp)+2θnγpθ′n)+PrEn(1+2ηγp)F′′2(1+1βp)+PrHpθ+Prθ′nF=0, | (9) |
ϕ″n(1+2ηγp)+2γpϕ′n+ScFϕ′n=0. | (10) |
and flow conditions reduced as:
F′=1,F=0,θn=1,ϕn=1,at η=0, |
F′=Ap,θn=0,ϕn=0,at η→∞. | (11) |
The parameters are given as:
βp=ˉμ√2πcτr,Rp=4σ∗ T3∞κk∗,γp=√υLc2U0, |
Ap=da,Mp=√σB20LρU0,Pr=ˉμcpκ,En=U20(ˆx/L)2cp(ˆTw−ˆT∞),GT=g0βT(ˆTw−ˆT∞)L2U0ˆx,GC= |
g0βC(ˆCw−ˆC∞)L2U0ˆx,Sc=νDm,Hp=LQ0U0ρcp. | (12) |
The Casson fluid parameter, Prandtl number, Eckert number, thermal radiation parameter, Solutal Grashof number, thermal Grashof number, velocities ratio parameter, and Schmidt number are symbolized by βp,Pr,En,Rp,GC,GT,ApandSc. The heat transfer at inclined surfaces is considered a key physical effect and to evaluate heat transfer we considered the Nusselt number. The mathematical relationship [26] of NN is:
Nuˆx=ˆxqwκ(ˆTw−ˆT∞),qw=−κ(1+16σ∗ˆT3∞3k∗κ)(∂ˆT∂ˆr)ˆr=Rc,Nuˆx√Reˆx=−(1+43Rp)θ′n(0)}. | (13) |
The lower order equations, Eqs (8)-(10), can be used to examine the non-Newtonian fluid flow regime. Due to their coupling and nonlinearity, they are challenging to solve analytically. The flow differential equations are reported using a number of various solution approaches [30,31,32,33,34,35], however, the shooting method will be employed to get numerical solution for the present research problem. The name of the shooting technique is taken from an analogy with target shooting: we shoot at the target and watch where it lands; based on the misses, we can modify our aim and shoot again in the hopes that it would hit close to the target. To convert boundary value problems into similar initial value problems, the shooting method was developed. The supportive material in this regard can be accessed in Refs. [25,26,27]. Having an initial value system is a vital first step in this regard. We can achieve it by presuming:
Y1=F(η),Y2=F′(η),Y3=F″(η), |
Y4=θn(η),Y5=θ′n(η), |
Y6=ϕn(η),Y7=ϕ′n(η), | (14) |
under Eq (14), Eqs (8)-(10) reduces to
Y′1=Y2, |
Y′2=Y3, |
Y′3=1(1+1βp)(1+2ηγp)[−2γpY3(1+1βp)+Y22−Y1Y3−GTY4cosα−GCY6cosα−Mp2(Y2−Ap)−Ap2], |
Y′4=Y5, |
Y′5=−1(1+43Rp)(1+2ηγp)+εT(1+2ηγpY4[(1+43Rp)(2γpY5)+εT((1+2ηγp)Y25+2γpY4Y5)+PrY1Y5+PrHpY4+PrEn(1+2ηγp)(1+1βp)Y23] |
Y′6=Y7, |
Y′7=−ScY1Y6−2γpY7(1+2ηγp). | (15) |
The flow equations turns down to:
Y1=0,Y2=1,Y6=1,Y4=1,at η=0, |
Y2→Ap,Y6→0,Y4→0,as η=∞. | (16) |
The self-coding scheme of shooting is carried out in Matlab for the present problem and outcomes are shared in terms of graphs and numerical data.
Mathematical modeling is done for the Casson flow over heated stretched surfaces. The equations for non-Newtonian flow are numerically solved. Both a flat plate and a cylindrical surface are used to evaluate the Nusselt number. The Nusselt number at surfaces is predicted using ANN models in both situations. Researchers firmly think that the multilayer perceptron (MLP) and ANN model [36,37,38,39,40] have excellent learning capabilities and can be utilized to anticipate a variety of physical occurrences. In MLP networks, three separate layers are utilized. The first layer uses the inputs, while the centrally important layer is referred to as the hidden layer. The output layer, which is the final layer, is where the prediction data is stored. The NN-ANN-I and II are two ANN models that we have created. Whereas NN-ANN-II is designed to provide predictions of Nusselt number values at cylindrical surfaces, NN-ANN-I is built to forecast the Nusselt number at flat surfaces. The five parameters are employed as inputs and the NN is taken into account as an output for both NN-ANN models. Variable thermal conductivity, heat generation, Prandtl number, the Casson fluid, and the thermal radiation parameters are the flow parameters.
Table 1 provides the symbolic information for models. Network architecture is given in Figure 1. We have gathered 80 samples for the five inputs, leading to 80 sample values for the Nusselt number. 15% of the data is used for validation while the remaining 70% is used to train the network. Testing also makes advantage of the 15%. Table 2 provides all available details in this regard. Ten (10) neurons are taken into account in the hidden layer, and the network is trained using the Levenberg-Marquardt method. The transfer functions used in hidden and output have the following expressions:
Model | Surface | Input | Output | ||||
NN-ANN Model-I | (γp=0.0) | (εT) | (Hp) | (Pr) | (βp) | (Rp) | NN |
NN-ANN Model-II | (γp=0.5) | (εT) | (Hp) | (Pr) | (βp) | (Rp) | NN |
Samples | NN-ANN Model-I | NN-ANN Model-II |
Total samples | 80 | 80 |
Training | 56 | 56 |
Validation | 12 | 12 |
Testing | 12 | 12 |
FT(x)=11+e−x, | (17) |
Pureline(x)=x. | (18) |
We have taken into consideration the mean squared error (MSE) and coefficient of determination (R) for performance analysis of ANN models. The mathematical formulas are expressed as:
MSE=1N∑Ni=1(Xnum(i)−XANN(i))2. | (19) |
R=√1−∑Ni=1(Xnum(i)−XANN(i))2∑Ni=1(Xnum(i))2. | (20) |
Using ANN models, the Casson flow toward two distinct heated surfaces is studied. The heat transfer coefficient and Casson fluid temperature are included in the quantity of interest. A mathematical formulation and shooting method is used to solve flow equations.
By taking the curvature parameter's zero value into consideration, the Casson fluid flow over a flat surface is assessed. The most important flow parameters for studying the flow field are the thermal radiation parameter, Prandt number, heat generation parameter, and variable thermal conductivity. We have assessed these parameters' effects on the NN at flat and cylindrical surfaces see Tables 3-10. Table 3 shows the variation in the NN in the presence of positive values for the variable thermal conductivity parameter. Both plates and cylinders are subject to this kind of inspection. The heat transfer normal to the surface is measured as a Nusselt number. The NN is seen to decrease for variations in the variable thermal conductivity parameter. For cylindrical surfaces, the Nusselt number is larger in magnitude. Table 4 shows the influence of a varied thermal conductivity parameter on the NN in the absence of heat generation.
εT |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −0.3868 | −0.5200 | 0.5415 | 0.7280 |
0.2 | −0.3782 | −0.5152 | 0.5295 | 0.7213 |
0.3 | −0.3702 | −0.5107 | 0.5183 | 0.7150 |
0.4 | −0.3625 | −0.5064 | 0.5075 | 0.7090 |
0.5 | −0.3553 | −0.5023 | 0.4974 | 0.7032 |
0.6 | −0.3484 | −0.4985 | 0.4878 | 0.6979 |
0.7 | −0.3419 | −0.4947 | 0.4787 | 0.6926 |
0.8 | −0.3357 | −0.4912 | 0.4700 | 0.6877 |
0.9 | −0.3298 | −0.4877 | 0.4617 | 0.6828 |
1.0 | −0.3187 | −0.4813 | 0.4462 | 0.6738 |
εT |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −1.3699 | −0.9213 | 1.9179 | 1.2898 |
0.2 | −1.3414 | −0.9107 | 1.8780 | 1.2750 |
0.3 | −1.3145 | −0.9005 | 1.8403 | 1.2607 |
0.4 | −1.2890 | −0.8908 | 1.8046 | 1.2471 |
0.5 | −1.2649 | −0.8815 | 1.7709 | 1.2341 |
0.6 | −1.2420 | −0.8726 | 1.7388 | 1.2216 |
0.7 | −1.2203 | −0.8641 | 1.7084 | 1.2097 |
0.8 | −1.1995 | −0.8559 | 1.6793 | 1.1983 |
0.9 | −1.1798 | −0.8480 | 1.6517 | 1.1872 |
1.0 | −1.1609 | −0.8404 | 1.6253 | 1.1766 |
The NN permits the inverse relation with varied heat conductivity for both surfaces. It is significant to notice that the Nusselt number's strength is higher for flat surfaces than for cylindrical surfaces. At the plate and cylinder, the effect of the Prandtl number on the NN is explored (see Table 5). It is obvious that the NN displays levels of stimulation for higher Prandtl numbers. When the influence of heat generation is taken into account, the effect of the Prandtl number on NN is explored (see Table 6). We have seen, the NN rises as the Prandtl number does for both plates and cylinders. Both plate and cylindrical surfaces are considered in the evaluation of the Casson fluid parameter's impact on the Nusselt number. About the scenario of heat generation, Table 7 is provided. We found that the Nusselt number declines normally when Casson fluid parameter values increase. Moreover, cylindrical surfaces have a somewhat greater Nusselt number magnitude.
In Table 8, the NN at both surfaces are analyzed in relation to the Casson fluid parameter for the case when heat generation has no impact. We found that heat transfer normal to the surface tends to incline for larger Casson fluid parameter values.
Pr |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.7 | −0.4278 | −0.5235 | 0.5989 | 0.7329 |
0.8 | −0.4588 | −0.5283 | 0.6423 | 0.7396 |
0.9 | −0.4879 | −0.5333 | 0.6831 | 0.7466 |
1.0 | −0.5152 | −0.5385 | 0.7213 | 0.7539 |
2.0 | −0.5652 | −0.5492 | 0.7913 | 0.7689 |
3.0 | −0.6114 | −0.5601 | 0.8560 | 0.7841 |
4.0 | −0.6537 | −0.5708 | 0.9152 | 0.7991 |
5.0 | −0.6931 | −0.5815 | 0.9703 | 0.8141 |
6.0 | −0.7301 | −0.5914 | 1.0221 | 0.8280 |
7.0 | −0.8893 | −0.6333 | 1.245 | 0.8866 |
Pr |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.7 | −0.4356 | −0.5686 | 0.6098 | 0.7960 |
0.8 | −0.4669 | −0.5781 | 0.6537 | 0.8093 |
0.9 | −0.4964 | −0.5875 | 0.6950 | 0.8225 |
1.0 | −0.5242 | −0.5969 | 0.7339 | 0.8357 |
2.0 | −0.7450 | −0.6839 | 1.0430 | 0.9575 |
3.0 | −0.9102 | −0.7546 | 1.2743 | 1.0564 |
4.0 | −1.0471 | −0.8102 | 1.4659 | 1.1343 |
5.0 | −1.1662 | −0.8547 | 1.6327 | 1.1966 |
6.0 | −1.2728 | −0.8911 | 1.7819 | 1.2475 |
7.0 | −1.3699 | −0.9213 | 1.9179 | 1.2898 |
βp |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −0.4052 | −0.6093 | 0.5673 | 0.8530 |
0.2 | −0.3953 | −0.6095 | 0.5534 | 0.8533 |
0.3 | −0.3911 | −0.6089 | 0.5475 | 0.8525 |
0.4 | −0.3888 | −0.6085 | 0.5443 | 0.8519 |
0.5 | −0.3863 | −0.6082 | 0.5408 | 0.8515 |
0.6 | −0.3855 | −0.6080 | 0.5397 | 0.8512 |
0.7 | −0.3850 | −0.6078 | 0.5390 | 0.8509 |
0.8 | −0.3845 | −0.6077 | 0.5383 | 0.8508 |
0.9 | −0.3840 | −0.6076 | 0.5376 | 0.8506 |
1.0 | −0.3838 | −0.6075 | 0.5373 | 0.8505 |
βp |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −0.2982 | −0.2125 | 0.4175 | 0.2975 |
0.2 | −0.3281 | −0.4017 | 0.4593 | 0.5624 |
0.3 | −0.3412 | −0.4635 | 0.4777 | 0.6489 |
0.4 | −0.3489 | −0.4940 | 0.4885 | 0.6916 |
0.5 | −0.3534 | −0.5120 | 0.4948 | 0.7168 |
0.6 | −0.3562 | −0.5238 | 0.4987 | 0.7333 |
0.7 | −0.3582 | −0.5283 | 0.5015 | 0.7396 |
0.8 | −0.3592 | −0.5384 | 0.5029 | 0.7538 |
0.9 | −0.3605 | −0.5432 | 0.5047 | 0.7605 |
1.0 | −0.4250 | −0.5470 | 0.5950 | 0.7658 |
Table 9 shows how the thermal radiation parameter affects the NN. NN increases as thermal radiation strength does. Also, in the case of cylindrical surfaces, the NN is larger. When the heat generation effect is absent, Table 10 shows how the thermal radiation parameter affects the NN at both surfaces. For both frames, we saw that greater radiation parameter values lead to larger NN values. It is significant to notice that in the case of flat plates, the Nusselt number strength is higher. Our interest is to predict the values of the NN by using ANN models. For the prediction of the Nusselt number at both surfaces, we have built two separate ANN models (NN-ANN-I and NN-ANN-II). At a flat surface, NN-ANN-I owns the Nusselt number values, whereas NN-ANN-II owns the Nusselt number values at a cylindrical surface. One crucial step in the development of an ANN model is network training. The network was trained using the Levenberg-Marquardt method. Figures 2 and 3 show the training results for NN-ANN models. Specifically, Figure 2 presents the training of the ANN model for the forecast of NN at a flat surface, whereas Figure 3 for the prediction of NN at a cylindrical surface. For NN-ANN-I, the finest validation (FV) is 0.029329 at epoch 5 and for NN-ANN-II, the FV is 0.0.0099961 at epoch 29. We noticed that the MSE values for both graphs start out high and then start to decline as the stages go-on. The error histogram (EH) of the NN-ANN models being built to forecast Nusselt number values is shown in Figures 4 and 5. The EH of model-I is shown in Figure 4, whereas the EH of the model-II for a cylindrical surface is shown in Figure 5. The EHs in both cases have 20 Bins.
Rp |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −0.5235 | −0.4278 | 0.5933 | 0.4848 |
0.2 | −0.5215 | −0.4058 | 0.6605 | 0.5140 |
0.3 | −0.5200 | −0.3868 | 0.7280 | 0.5415 |
0.4 | −0.5188 | −0.3703 | 0.7954 | 0.5678 |
0.5 | −0.5179 | −0.3557 | 0.8631 | 0.5929 |
0.6 | −0.5171 | −0.3428 | 0.9308 | 0.6170 |
0.7 | −0.5166 | −0.3314 | 0.9988 | 0.6407 |
0.8 | −0.5161 | −0.3212 | 1.0667 | 0.6638 |
0.9 | −0.5157 | −0.3121 | 1.1345 | 0.6867 |
1.0 | −0.5154 | −0.3038 | 1.2026 | 0.7089 |
Rp |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −1.4994 | −0.9564 | 1.6993 | 1.0839 |
0.2 | −1.4302 | −0.9383 | 1.8116 | 1.1885 |
0.3 | −1.3699 | −0.9213 | 1.9179 | 1.2898 |
0.4 | −1.3167 | −0.9054 | 2.0189 | 1.3883 |
0.5 | −1.2693 | −0.8905 | 2.1155 | 1.4842 |
0.6 | −1.2267 | −0.8765 | 2.2081 | 1.5777 |
0.7 | −1.1882 | −0.8633 | 2.2972 | 1.669 |
0.8 | −1.1530 | −0.8509 | 2.3829 | 1.7585 |
0.9 | −1.1208 | −0.8392 | 2.4658 | 1.8462 |
1.0 | −1.0912 | −0.8281 | 2.5461 | 1.9322 |
Since both NN-ANN models' error values are quite low, the training stages of NN-ANN-I and NN-ANN-II to estimate Nusselt numbers at both surfaces are successfully finished. The mean square error graphs for the Nusselt number are shown in Figures 6 and 7. Here, Figure 6 shows the MSE for the NN at the plat surface, and Figure 7 shows the MSE for the NN at cylinder. Both figures show that NN-ANN-I and NN-ANN-II have successfully completed their learning stages. For use of 80 data samples, the MSE values are displayed. The MSE values' proximity to zero indicates that the NN-ANN models were trained with fewer mistakes to predict the Nusselt number. While MSE is recorded as 0.03261 for flat surface, and the average MSE for cylindrical surfaces is 0.008562777. The regression outcomes of Nusselt number data for both surfaces are given in Figures 8 and 9. It is key to remember that the regression value provides a correlation between the objectives and forecasts values; if this value is near to one (1), we can state there is a close link.
The developed NN-ANN-I and NN-ANN-II are the best models to forecast the NN at both surfaces, according to Figures 8 and 9.
Figures 10 and 11 compare the Nusselt number's actual values with those predicted by an artificial neural networking model. Figure 10 compares NN-ANN and the target data set for plate whereas Figure 11 compares NN-ANN and the target data set for cylinder. We can note from figures that the bulk of NN-ANN model outputs is in good agreement with the Nusselt number target values. Due to this overlap, we may conclude that the generated models, NN-ANN-I and NN-ANN-II, are highly accurate in predicting Nusselt numbers.
The effect of a positive modification in the Hp on temperature is shown in Figures 12 and 13. Figure 12 shows the impact of heat generation on temperature for a plate, while Figure 13 shows the impact of heat generation on temperature for a cylindrical surface. It can be demonstrated that the temperature of the Casson fluid and the heat generation parameter are directly related. Positive values of the Hp cause energy generation, which raises the temperature. The effect of the Pr on the temperature is seen in Figures 14 and 15. In more detail, Figures 14 and 15 show how the Prandt number affects the temperature of a fluid over a flat surface and a cylindrical surface, respectively. We have seen that the Prandt number enables an indirect relation to Casson fluid temperature.
As the Prandtl number admits an inverse relationship to thermal conductivity, larger Prandtl number values result in a decrease in thermal conductivity. The fluid temperature drops as a result. For both flat and cylindrical surfaces, the effect of varying thermal conductivity on Casson fluid temperature is examined. Figures 16 and 17 provide evidence in this regard. Figure 16 describes the effect of variable thermal conductivity on the temperature in more detail, and we found that temperature rises as the εT rises. Similar to Figure 16, Figure 17 demonstrates how εT affects temperature over a cylinder. We found that fluid temperature exhibits tendencies that are indicative of variations in εT. Table 11 provides performance details for both ANN models, and it is clear that the constructed models are capable of accurately forecasting the values of the heat transfer coefficient at both flat and cylindrical surfaces.
MSE | R | ||||||
Model | Surface | Training | Validation | Test | Training | Validation | Test |
NN-Model-I | (γp=0.0) | 0.0350 | 0.0293 | 0.0246 | 0.9547 | 0.9727 | 0.9619 |
NN-Model-II | (γp=0.5) | 0.0094 | 0.0100 | 0.0032 | 0.9549 | 0.9442 | 0.9605 |
Thermally magnetized Casson stagnation flow over an inclined surfaces is studied numerically. Mixed convection, thermal radiations, heat generation, viscous dissipation, and temperature-dependent thermal conductivity all contribute to the flow. Using artificial neural networking models, the NN is evaluated at both surfaces toward five various flow factors. The following are the key outcomes:
• For an inclined cylindrical surface, an average Mean Square Error (MSE) is recorded as 0.008562777 while for a flat surface it is recorded as 0.03261.
• The values of the coefficient of determination suggest artificial neural networks as the best method for predicting the NN at both surfaces.
• For both surfaces, the temperature admits direct relation for heat generation and thermal conductivity parameters.
• For higher Prandtl number, the Casson fluid temperature shows decline values. The such trend holds for both surfaces.
• According to the predictions made by NN-ANN models, we found that for cylindrical surface, the strength of the NN is greater than the flat surface.
• The Nusselt number admits indirect relation for variable thermal conductivity, and Casson parameters while the opposite is the case for thermal radiation parameter and Prandtl number.
• A neural networking approach can be used to predict the surface quantities being involved for the fluid flow over non-linear stretched surfaces having engineering standpoints.
The authors would like to thank Prince Sultan University, Saudi Arabia, for the technical support through the TAS research lab. Author (Haitham M. S. Bahaidara) is thankful for the financial support provided by the Interdisciplinary Research Center for Renewable Energy and Power Systems, King Fahd University of Petroleum and Minerals, for funding this research through Project no: INRE2212.
The authors declare no conflict of interest.
[1] |
J. Z. Shing, J. Corbin, A. R. Kreimer, L. J. Carvajal, K. Taparra, M. S. Shiels, et al., Human papillomavirus–associated cancer incidence by disaggregated Asian American, Native Hawaiian, and other Pacific Islander ethnicity, JNCI Cancer Spectrum, 7 (2023), 1–9. https://doi.org/10.1093/jncics/pkad012 doi: 10.1093/jncics/pkad012
![]() |
[2] |
J. Galon, D. Daniela, Tumor immunology and tumor evolution: intertwined histories, Immunity, 52 (2020), 55–81. https://doi.org/10.1016/j.immuni.2019.12.018 doi: 10.1016/j.immuni.2019.12.018
![]() |
[3] |
W. H. Fridman, F. Pagès, C. Sautès-Fridman, J. Galon, The immune contexture in human tumors: impact on clinical outcome, Nat. Rev. Cancer, 12 (2012), 298–306. https://doi.org/10.1038/nrc3245 doi: 10.1038/nrc3245
![]() |
[4] |
H. Song, C. Ruan, Y. Xu, T. Xu, R. Fan, T. Jiang, et al., Survival stratification for colorectal cancer via multi-omics integration using an autoencoder-based model, Exp. Biol. Med., 247 (2022), 898–909. https://doi.org/10.1177/15353702211065010 doi: 10.1177/15353702211065010
![]() |
[5] |
G. Rajput, S. Agrawal, K. Biyani, S. Vishvakarma, Early breast cancer diagnosis using cogent activation function-based deep learning implementation on screened mammograms, Int. J. Imaging Syst. Technol., 32 (2022), 1101–1118. https://doi.org/10.1002/ima.22701 doi: 10.1002/ima.22701
![]() |
[6] |
A. Abougarair, A. Oun, S. Sawan, A. Ma'arif, Deep learning-based automated approach for classifying bacterial images, Int. J. Robotics Control Syst., 4 (2024), 849–876. https://doi.org/10.31763/ijrcs.v4i2.1423 doi: 10.31763/ijrcs.v4i2.1423
![]() |
[7] |
M. Itik, M. U. Salamci, S. P. Banks, Optimal control of drug therapy in cancer treatment, Nonlinear Anal., 71 (2009), e1473–e1486. https://doi.org/10.1016/j.na.2009.01.214 doi: 10.1016/j.na.2009.01.214
![]() |
[8] |
T. Çimen, Systematic and effective design of nonlinear feedback controllers via the state-dependent Riccati equation (SDRE) method, Annu. Rev. Control, 34 (2010), 32–51. https://doi.org/10.1016/j.arcontrol.2010.03.001 doi: 10.1016/j.arcontrol.2010.03.001
![]() |
[9] |
Y. Batmani, H. Khaloozadeh, Optimal chemotherapy in cancer treatment: state dependent Riccati equation control and extended Kalman filter, Optim. Contr. Appl. Met., 34 (2013), 562–577. https://doi.org/10.1002/oca.2039 doi: 10.1002/oca.2039
![]() |
[10] |
L. De Pillis, A. Radunskaya, The dynamics of an optimally controlled tumor model: a case study. Math. Comput. Model., 37 (2003), 1221–1244. https://doi.org/10.1016/S0895-7177(03)00133-X doi: 10.1016/S0895-7177(03)00133-X
![]() |
[11] |
F. A. Rihan, N. F. Rihan, Dynamics of cancer-immune system with external treatment and optimal control, J. Cancer Sci. Ther., 8 (2016), 257–261. https://doi.org/10.4172/1948-5956.1000423 doi: 10.4172/1948-5956.1000423
![]() |
[12] |
M. Sarhaddi, M. Yaghoobi, A new approach in cancer treatment regimen using adaptive fuzzy back-stepping sliding mode control and tumor-immunity fractional order model, Biocybern. Biomed. Eng., 40 (2020), 1654–1665. https://doi.org/10.1016/j.bbe.2020.09.003 doi: 10.1016/j.bbe.2020.09.003
![]() |
[13] |
V. A. Kuznetsov, I. A. Makalkin, M. A. Taylor, A. S. Perelson, Nonlinear dynamics of immunogenic tumors: parameter estimation and global bifurcation analysis, Bull. Math. Biol., 56 (1994), 295–321. https://doi.org/10.1016/S0092-8240(05)80260-5 doi: 10.1016/S0092-8240(05)80260-5
![]() |
[14] |
F. Angaroni, A. Graudenzi, M. Rossignolo, D. Maspero, T. Calarco, R. Piazza, et al., An optimal control framework for the automated design of personalized cancer treatments, Front. Bioeng. Biotechnol, 8 (2020), 523. https://doi.org/10.3389/fbioe.2020.00523 doi: 10.3389/fbioe.2020.00523
![]() |
[15] |
J. J. Cunningham, J. S. Brown, R. A. Gatenby, K. Staňková, Optimal control to develop therapeutic strategies for metastatic castrate resistant prostate cancer, J. Theor. Biol., 459 (2018), 67–78. https://doi.org/10.1016/j.jtbi.2018.09.022 doi: 10.1016/j.jtbi.2018.09.022
![]() |
[16] |
Y. Yang, C. Y. Shiranthika, C. Y. Wang, K. W. Chen, S. Sumathipala, Reinforcement learning strategies in cancer chemotherapy treatments: a review, Comput. Meth. Prog. Bio., 229 (2023), 107280. https://doi.org/10.1016/j.cmpb.2022.107280 doi: 10.1016/j.cmpb.2022.107280
![]() |
[17] |
P. Samadi, P. Sarvarian, E. Gholipour, K. S. Asenjan, L. Aghebati-Maleki, R. Motavalli, et al., Berberine: a novel therapeutic strategy for cancer, IUBMB Life, 72 (2020), 2065–2079. https://doi.org/10.1002/iub.2350 doi: 10.1002/iub.2350
![]() |
[18] |
A. M. Jarrett, D. Faghihi, D. A. Hormuth, E. A. B. F. Lima, J. Virostko, G. Biros, et al., Optimal control theory for personalized therapeutic regimens in oncology: Background, history, challenges, and opportunities, J. Clin. Med., 9 (2020), 1314. https://doi.org/10.3390/jcm9051314 doi: 10.3390/jcm9051314
![]() |
[19] | H. Schättler, U. Ledzewicz, Optimal control for mathematical models of cancer therapies: an application of geometric methods, Vol. 42, Springer, 2015. https://doi.org/10.1007/978-1-4939-2972-6 |
[20] |
S. Oke, M. B. Matadi, S. Xulu, Optimal control analysis of a mathematical model for breast cancer, Math. Comput. Appl., 23 (2018), 21. https://doi.org/10.3390/mca23020021 doi: 10.3390/mca23020021
![]() |
[21] |
M. Engelhart, D. Lebiedz, S. Sager, Optimal control for selected cancer chemotherapy ODE models: a view on the potential of optimal schedules and choice of objective function, Math. Biosci., 229 (2011), 123–134. https://doi.org/10.1016/j.mbs.2010.11.007 doi: 10.1016/j.mbs.2010.11.007
![]() |
[22] |
A. Cappuccio, F. Castiglione, B. Piccoli, Determination of the optimal therapeutic protocols in cancer immunotherapy, Math. Biosci., 209 (2007), 1–13. https://doi.org/10.1016/j.mbs.2007.02.009 doi: 10.1016/j.mbs.2007.02.009
![]() |
[23] |
M. Gluzman, J. G. Scott, A. Vladimirsky, Optimizing adaptive cancer therapy: dynamic programming and evolutionary game theory, Proc. Biol. Sci., 287 (2020), 20192454. https://doi.org/10.1098/rspb.2019.2454 doi: 10.1098/rspb.2019.2454
![]() |
[24] | J. A. Adam, N. Bellomo, A survey of models for tumor-immune system dynamics, Springer Science & Business Media, 2012. |
[25] | J. C. Arciero, T. L. Jackson, D. E. Kirschner, A mathematical model of tumor-immune evasion and siRNA treatment, Discrete Cont. Dyn. Syst.-B, 4 (2004), 39–58. |
[26] | M. Chaplain, A. Matzavinos, Mathematical modelling of spatio-temporal phenomena in tumour immunology, In: A. Friedman, Tutorials in mathematical biosciences III, Lecture Notes in Mathematics, Springer, 1872 (2006), 131–183. https://doi.org/10.1007/11561606_4 |
[27] |
R. J. De Boer, P. Hogeweg, H. F. Dullens, R. A. De Weger, W. Den Otter, Macrophage T lymphocyte interactions in the anti-tumor immune response: a mathematical model, J. Immunol., 134 (1985), 2748–2758. https://doi.org/10.4049/jimmunol.134.4.2748 doi: 10.4049/jimmunol.134.4.2748
![]() |
[28] |
L. G. de Pillis, W. Gu, A. E. Radunskaya, Mixed immunotherapy and chemotherapy of tumors: modeling, applications and biological interpretations, J. Theor. Biol., 238 (2006), 841–862. https://doi.org/10.1016/j.jtbi.2005.06.037 doi: 10.1016/j.jtbi.2005.06.037
![]() |
[29] |
L. de Pillis, K. R. Fister, W. Gu, C. Collins, M. Daub, D. Gross, et al., Mathematical model creation for cancer chemo-immunotherapy, Comput. Math. Methods Med., 10 (2009), 165–184, 2009. https://doi.org/10.1080/17486700802216301 doi: 10.1080/17486700802216301
![]() |
[30] | D. Kirk, Optimal control theory: an introduction, Dover Publications, 2004. |
[31] |
D. Jacobson, M. Lele, A transformation technique for optimal control problems with a state variable inequality constraint, IEEE Trans. Automat. Control, 14 (1969), 457–464. https://doi.org/10.1109/TAC.1969.1099283 doi: 10.1109/TAC.1969.1099283
![]() |
[32] | D. A. Redfern, C. J. Goh, Feedback control of state constrained optimal control problems, In: J. Doležal, J. Fidler, System modelling and optimization, IFIP-The International Federation for Information Processing, Springer, 1996,442–449. https://doi.org/10.1007/978-0-387-34897-1_53 |
[33] | L. S. Pontryagin, Mathematical theory of optimal processes, 1 Ed., CRC press, 1987. https://doi.org/10.1201/9780203749319 |
[34] | H. Khalil, Nonlinear systems, 3 Eds., Prentice Hall, 2002. |
[35] |
J. A. E. Andersson, J. Gillis, G. Horn, J. B. Rawlings, M. Diehl, CasADi: a software framework for nonlinear optimization and optimal control, Math. Prog. Comput., 11 (2019), 1–36. https://doi.org/10.1007/s12532-018-0139-4 doi: 10.1007/s12532-018-0139-4
![]() |
[36] |
A. Wächter, L. Biegler, On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming, Math. Program., 106 (2006), 25–57. https://doi.org/10.1007/s10107-004-0559-y doi: 10.1007/s10107-004-0559-y
![]() |
[37] |
E. V. Kumar, J. Jerome, G. Raaja, State dependent Riccati equation based nonlinear controller design for ball and beam system, Procedia Eng., 97 (2014), 1896–1905. https://doi.org/10.1016/j.proeng.2014.12.343 doi: 10.1016/j.proeng.2014.12.343
![]() |
[38] |
K. Hałas, E. Krysiak, T. Hałas, S. Stępień, Numerical solution of SDRE control problem–comparison of the selected methods, Found. Comput. Decis. Sci., 45 (2020), 79–95. https://doi.org/10.2478/fcds-2020-0006 doi: 10.2478/fcds-2020-0006
![]() |
[39] |
A. J. Abougarair, S. E. Elwefati, Identification and control of epidemic disease based neural networks and optimization technique, Int. J. Rob. Control Syst., 3 (2023), 780–803. https://doi.org/10.31763/ijrcs.v3i4.1151 doi: 10.31763/ijrcs.v3i4.1151
![]() |
[40] | S. E. Elwefati, A. J. Abougarair, M. M. Bakush, Control of epidemic disease based optimization technique, 2021 IEEE 1st International Maghreb Meeting of the Conference on Sciences and Techniques of Automatic Control and Computer Engineering MI-STA, 2021, 52–57. https://doi.org/10.1109/MI-STA52233.2021.9464453 |
[41] | A. J. Abougarair, A. A. Oun, S. I. Sawan, T. Abougard, H. Maghfiroh, Cancer treatment precision strategies through optimal control theory, J. Rob. Control, 5 (2024), 1261–1290. |
[42] |
T. Yuan, G. Guan, S. Shen, L. Zhu, Stability analysis and optimal control of epidemic-like transmission model with nonlinear inhibition mechanism and time delay in both homogeneous and heterogeneous networks, J. Math. Anal. Appl., 526 (2023), 127273. https://doi.org/10.1016/j.jmaa.2023.127273 doi: 10.1016/j.jmaa.2023.127273
![]() |
[43] |
L. Zhu, T. Yuan, Optimal control and parameter identification of a reaction–diffusion network propagation model, Nonlinear Dyn., 111 (2023), 21707–21733. https://doi.org/10.1007/s11071-023-08949-y doi: 10.1007/s11071-023-08949-y
![]() |
[44] |
Y. Ke, L. Zhu, P. Wu, L. Shi, Dynamics of a reaction-diffusion rumor propagation model with non-smooth control, Appl. Math. Comput., 435 (2022), 127478. https://doi.org/10.1016/j.amc.2022.127478 doi: 10.1016/j.amc.2022.127478
![]() |
[45] |
L. Zhu, X. Tao, S. Shen, Pattern dynamics in a reaction-diffusion predator–prey model with Allee effect based on network and non-network environments, Eng. Appl. Artif. Intel., 128 (2024), 107491. https://doi.org/10.1016/j.engappai.2023.107491 doi: 10.1016/j.engappai.2023.107491
![]() |
[46] |
B. Li, L. Zhu, Turing instability analysis of a reaction–diffusion system for rumor propagation in continuous space and complex networks, Inform. Process. Manag., 61 (2024), 103621. https://doi.org/10.1016/j.ipm.2023.103621 doi: 10.1016/j.ipm.2023.103621
![]() |
Model | Surface | Input | Output | ||||
NN-ANN Model-I | (γp=0.0) | (εT) | (Hp) | (Pr) | (βp) | (Rp) | NN |
NN-ANN Model-II | (γp=0.5) | (εT) | (Hp) | (Pr) | (βp) | (Rp) | NN |
Samples | NN-ANN Model-I | NN-ANN Model-II |
Total samples | 80 | 80 |
Training | 56 | 56 |
Validation | 12 | 12 |
Testing | 12 | 12 |
εT |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −0.3868 | −0.5200 | 0.5415 | 0.7280 |
0.2 | −0.3782 | −0.5152 | 0.5295 | 0.7213 |
0.3 | −0.3702 | −0.5107 | 0.5183 | 0.7150 |
0.4 | −0.3625 | −0.5064 | 0.5075 | 0.7090 |
0.5 | −0.3553 | −0.5023 | 0.4974 | 0.7032 |
0.6 | −0.3484 | −0.4985 | 0.4878 | 0.6979 |
0.7 | −0.3419 | −0.4947 | 0.4787 | 0.6926 |
0.8 | −0.3357 | −0.4912 | 0.4700 | 0.6877 |
0.9 | −0.3298 | −0.4877 | 0.4617 | 0.6828 |
1.0 | −0.3187 | −0.4813 | 0.4462 | 0.6738 |
εT |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −1.3699 | −0.9213 | 1.9179 | 1.2898 |
0.2 | −1.3414 | −0.9107 | 1.8780 | 1.2750 |
0.3 | −1.3145 | −0.9005 | 1.8403 | 1.2607 |
0.4 | −1.2890 | −0.8908 | 1.8046 | 1.2471 |
0.5 | −1.2649 | −0.8815 | 1.7709 | 1.2341 |
0.6 | −1.2420 | −0.8726 | 1.7388 | 1.2216 |
0.7 | −1.2203 | −0.8641 | 1.7084 | 1.2097 |
0.8 | −1.1995 | −0.8559 | 1.6793 | 1.1983 |
0.9 | −1.1798 | −0.8480 | 1.6517 | 1.1872 |
1.0 | −1.1609 | −0.8404 | 1.6253 | 1.1766 |
Pr |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.7 | −0.4278 | −0.5235 | 0.5989 | 0.7329 |
0.8 | −0.4588 | −0.5283 | 0.6423 | 0.7396 |
0.9 | −0.4879 | −0.5333 | 0.6831 | 0.7466 |
1.0 | −0.5152 | −0.5385 | 0.7213 | 0.7539 |
2.0 | −0.5652 | −0.5492 | 0.7913 | 0.7689 |
3.0 | −0.6114 | −0.5601 | 0.8560 | 0.7841 |
4.0 | −0.6537 | −0.5708 | 0.9152 | 0.7991 |
5.0 | −0.6931 | −0.5815 | 0.9703 | 0.8141 |
6.0 | −0.7301 | −0.5914 | 1.0221 | 0.8280 |
7.0 | −0.8893 | −0.6333 | 1.245 | 0.8866 |
Pr |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.7 | −0.4356 | −0.5686 | 0.6098 | 0.7960 |
0.8 | −0.4669 | −0.5781 | 0.6537 | 0.8093 |
0.9 | −0.4964 | −0.5875 | 0.6950 | 0.8225 |
1.0 | −0.5242 | −0.5969 | 0.7339 | 0.8357 |
2.0 | −0.7450 | −0.6839 | 1.0430 | 0.9575 |
3.0 | −0.9102 | −0.7546 | 1.2743 | 1.0564 |
4.0 | −1.0471 | −0.8102 | 1.4659 | 1.1343 |
5.0 | −1.1662 | −0.8547 | 1.6327 | 1.1966 |
6.0 | −1.2728 | −0.8911 | 1.7819 | 1.2475 |
7.0 | −1.3699 | −0.9213 | 1.9179 | 1.2898 |
βp |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −0.4052 | −0.6093 | 0.5673 | 0.8530 |
0.2 | −0.3953 | −0.6095 | 0.5534 | 0.8533 |
0.3 | −0.3911 | −0.6089 | 0.5475 | 0.8525 |
0.4 | −0.3888 | −0.6085 | 0.5443 | 0.8519 |
0.5 | −0.3863 | −0.6082 | 0.5408 | 0.8515 |
0.6 | −0.3855 | −0.6080 | 0.5397 | 0.8512 |
0.7 | −0.3850 | −0.6078 | 0.5390 | 0.8509 |
0.8 | −0.3845 | −0.6077 | 0.5383 | 0.8508 |
0.9 | −0.3840 | −0.6076 | 0.5376 | 0.8506 |
1.0 | −0.3838 | −0.6075 | 0.5373 | 0.8505 |
βp |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −0.2982 | −0.2125 | 0.4175 | 0.2975 |
0.2 | −0.3281 | −0.4017 | 0.4593 | 0.5624 |
0.3 | −0.3412 | −0.4635 | 0.4777 | 0.6489 |
0.4 | −0.3489 | −0.4940 | 0.4885 | 0.6916 |
0.5 | −0.3534 | −0.5120 | 0.4948 | 0.7168 |
0.6 | −0.3562 | −0.5238 | 0.4987 | 0.7333 |
0.7 | −0.3582 | −0.5283 | 0.5015 | 0.7396 |
0.8 | −0.3592 | −0.5384 | 0.5029 | 0.7538 |
0.9 | −0.3605 | −0.5432 | 0.5047 | 0.7605 |
1.0 | −0.4250 | −0.5470 | 0.5950 | 0.7658 |
Rp |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −0.5235 | −0.4278 | 0.5933 | 0.4848 |
0.2 | −0.5215 | −0.4058 | 0.6605 | 0.5140 |
0.3 | −0.5200 | −0.3868 | 0.7280 | 0.5415 |
0.4 | −0.5188 | −0.3703 | 0.7954 | 0.5678 |
0.5 | −0.5179 | −0.3557 | 0.8631 | 0.5929 |
0.6 | −0.5171 | −0.3428 | 0.9308 | 0.6170 |
0.7 | −0.5166 | −0.3314 | 0.9988 | 0.6407 |
0.8 | −0.5161 | −0.3212 | 1.0667 | 0.6638 |
0.9 | −0.5157 | −0.3121 | 1.1345 | 0.6867 |
1.0 | −0.5154 | −0.3038 | 1.2026 | 0.7089 |
Rp |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −1.4994 | −0.9564 | 1.6993 | 1.0839 |
0.2 | −1.4302 | −0.9383 | 1.8116 | 1.1885 |
0.3 | −1.3699 | −0.9213 | 1.9179 | 1.2898 |
0.4 | −1.3167 | −0.9054 | 2.0189 | 1.3883 |
0.5 | −1.2693 | −0.8905 | 2.1155 | 1.4842 |
0.6 | −1.2267 | −0.8765 | 2.2081 | 1.5777 |
0.7 | −1.1882 | −0.8633 | 2.2972 | 1.669 |
0.8 | −1.1530 | −0.8509 | 2.3829 | 1.7585 |
0.9 | −1.1208 | −0.8392 | 2.4658 | 1.8462 |
1.0 | −1.0912 | −0.8281 | 2.5461 | 1.9322 |
MSE | R | ||||||
Model | Surface | Training | Validation | Test | Training | Validation | Test |
NN-Model-I | (γp=0.0) | 0.0350 | 0.0293 | 0.0246 | 0.9547 | 0.9727 | 0.9619 |
NN-Model-II | (γp=0.5) | 0.0094 | 0.0100 | 0.0032 | 0.9549 | 0.9442 | 0.9605 |
Model | Surface | Input | Output | ||||
NN-ANN Model-I | (γp=0.0) | (εT) | (Hp) | (Pr) | (βp) | (Rp) | NN |
NN-ANN Model-II | (γp=0.5) | (εT) | (Hp) | (Pr) | (βp) | (Rp) | NN |
Samples | NN-ANN Model-I | NN-ANN Model-II |
Total samples | 80 | 80 |
Training | 56 | 56 |
Validation | 12 | 12 |
Testing | 12 | 12 |
εT |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −0.3868 | −0.5200 | 0.5415 | 0.7280 |
0.2 | −0.3782 | −0.5152 | 0.5295 | 0.7213 |
0.3 | −0.3702 | −0.5107 | 0.5183 | 0.7150 |
0.4 | −0.3625 | −0.5064 | 0.5075 | 0.7090 |
0.5 | −0.3553 | −0.5023 | 0.4974 | 0.7032 |
0.6 | −0.3484 | −0.4985 | 0.4878 | 0.6979 |
0.7 | −0.3419 | −0.4947 | 0.4787 | 0.6926 |
0.8 | −0.3357 | −0.4912 | 0.4700 | 0.6877 |
0.9 | −0.3298 | −0.4877 | 0.4617 | 0.6828 |
1.0 | −0.3187 | −0.4813 | 0.4462 | 0.6738 |
εT |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −1.3699 | −0.9213 | 1.9179 | 1.2898 |
0.2 | −1.3414 | −0.9107 | 1.8780 | 1.2750 |
0.3 | −1.3145 | −0.9005 | 1.8403 | 1.2607 |
0.4 | −1.2890 | −0.8908 | 1.8046 | 1.2471 |
0.5 | −1.2649 | −0.8815 | 1.7709 | 1.2341 |
0.6 | −1.2420 | −0.8726 | 1.7388 | 1.2216 |
0.7 | −1.2203 | −0.8641 | 1.7084 | 1.2097 |
0.8 | −1.1995 | −0.8559 | 1.6793 | 1.1983 |
0.9 | −1.1798 | −0.8480 | 1.6517 | 1.1872 |
1.0 | −1.1609 | −0.8404 | 1.6253 | 1.1766 |
Pr |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.7 | −0.4278 | −0.5235 | 0.5989 | 0.7329 |
0.8 | −0.4588 | −0.5283 | 0.6423 | 0.7396 |
0.9 | −0.4879 | −0.5333 | 0.6831 | 0.7466 |
1.0 | −0.5152 | −0.5385 | 0.7213 | 0.7539 |
2.0 | −0.5652 | −0.5492 | 0.7913 | 0.7689 |
3.0 | −0.6114 | −0.5601 | 0.8560 | 0.7841 |
4.0 | −0.6537 | −0.5708 | 0.9152 | 0.7991 |
5.0 | −0.6931 | −0.5815 | 0.9703 | 0.8141 |
6.0 | −0.7301 | −0.5914 | 1.0221 | 0.8280 |
7.0 | −0.8893 | −0.6333 | 1.245 | 0.8866 |
Pr |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.7 | −0.4356 | −0.5686 | 0.6098 | 0.7960 |
0.8 | −0.4669 | −0.5781 | 0.6537 | 0.8093 |
0.9 | −0.4964 | −0.5875 | 0.6950 | 0.8225 |
1.0 | −0.5242 | −0.5969 | 0.7339 | 0.8357 |
2.0 | −0.7450 | −0.6839 | 1.0430 | 0.9575 |
3.0 | −0.9102 | −0.7546 | 1.2743 | 1.0564 |
4.0 | −1.0471 | −0.8102 | 1.4659 | 1.1343 |
5.0 | −1.1662 | −0.8547 | 1.6327 | 1.1966 |
6.0 | −1.2728 | −0.8911 | 1.7819 | 1.2475 |
7.0 | −1.3699 | −0.9213 | 1.9179 | 1.2898 |
βp |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −0.4052 | −0.6093 | 0.5673 | 0.8530 |
0.2 | −0.3953 | −0.6095 | 0.5534 | 0.8533 |
0.3 | −0.3911 | −0.6089 | 0.5475 | 0.8525 |
0.4 | −0.3888 | −0.6085 | 0.5443 | 0.8519 |
0.5 | −0.3863 | −0.6082 | 0.5408 | 0.8515 |
0.6 | −0.3855 | −0.6080 | 0.5397 | 0.8512 |
0.7 | −0.3850 | −0.6078 | 0.5390 | 0.8509 |
0.8 | −0.3845 | −0.6077 | 0.5383 | 0.8508 |
0.9 | −0.3840 | −0.6076 | 0.5376 | 0.8506 |
1.0 | −0.3838 | −0.6075 | 0.5373 | 0.8505 |
βp |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −0.2982 | −0.2125 | 0.4175 | 0.2975 |
0.2 | −0.3281 | −0.4017 | 0.4593 | 0.5624 |
0.3 | −0.3412 | −0.4635 | 0.4777 | 0.6489 |
0.4 | −0.3489 | −0.4940 | 0.4885 | 0.6916 |
0.5 | −0.3534 | −0.5120 | 0.4948 | 0.7168 |
0.6 | −0.3562 | −0.5238 | 0.4987 | 0.7333 |
0.7 | −0.3582 | −0.5283 | 0.5015 | 0.7396 |
0.8 | −0.3592 | −0.5384 | 0.5029 | 0.7538 |
0.9 | −0.3605 | −0.5432 | 0.5047 | 0.7605 |
1.0 | −0.4250 | −0.5470 | 0.5950 | 0.7658 |
Rp |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −0.5235 | −0.4278 | 0.5933 | 0.4848 |
0.2 | −0.5215 | −0.4058 | 0.6605 | 0.5140 |
0.3 | −0.5200 | −0.3868 | 0.7280 | 0.5415 |
0.4 | −0.5188 | −0.3703 | 0.7954 | 0.5678 |
0.5 | −0.5179 | −0.3557 | 0.8631 | 0.5929 |
0.6 | −0.5171 | −0.3428 | 0.9308 | 0.6170 |
0.7 | −0.5166 | −0.3314 | 0.9988 | 0.6407 |
0.8 | −0.5161 | −0.3212 | 1.0667 | 0.6638 |
0.9 | −0.5157 | −0.3121 | 1.1345 | 0.6867 |
1.0 | −0.5154 | −0.3038 | 1.2026 | 0.7089 |
Rp |
θn′(0) | −(1+4/3Rp) | ||
γp=0 (Plate) | γp=0.5(Cylinder) | γp=0 (Plate) | γp=0.5 (Cylinder) | |
0.1 | −1.4994 | −0.9564 | 1.6993 | 1.0839 |
0.2 | −1.4302 | −0.9383 | 1.8116 | 1.1885 |
0.3 | −1.3699 | −0.9213 | 1.9179 | 1.2898 |
0.4 | −1.3167 | −0.9054 | 2.0189 | 1.3883 |
0.5 | −1.2693 | −0.8905 | 2.1155 | 1.4842 |
0.6 | −1.2267 | −0.8765 | 2.2081 | 1.5777 |
0.7 | −1.1882 | −0.8633 | 2.2972 | 1.669 |
0.8 | −1.1530 | −0.8509 | 2.3829 | 1.7585 |
0.9 | −1.1208 | −0.8392 | 2.4658 | 1.8462 |
1.0 | −1.0912 | −0.8281 | 2.5461 | 1.9322 |
MSE | R | ||||||
Model | Surface | Training | Validation | Test | Training | Validation | Test |
NN-Model-I | (γp=0.0) | 0.0350 | 0.0293 | 0.0246 | 0.9547 | 0.9727 | 0.9619 |
NN-Model-II | (γp=0.5) | 0.0094 | 0.0100 | 0.0032 | 0.9549 | 0.9442 | 0.9605 |