Laser welding of aluminum battery tab to variable Al/Cu busbars in Li-ion battery joint

Ehsan Harati; Paul Kah; Ehsan Harati; Paul Kah

doi:10.3934/matersci.2022053

AIMS Materials Science

2022, Volume 9, Issue 6: 884-918. doi: 10.3934/matersci.2022053

Previous Article Next Article

Review Topical Sections

Laser welding of aluminum battery tab to variable Al/Cu busbars in Li-ion battery joint

Ehsan Harati ^,,
Paul Kah

Department of Engineering Science, University West, SE-46186 Trollhättan, Sweden

Received: 09 June 2022 Revised: 29 August 2022 Accepted: 10 September 2022 Published: 07 November 2022

The trend is shifting from internal combustion engines (ICEs) to battery electric vehicles (BEVs). One of the important battery joints is battery tabs to the busbar connection. Aluminum (Al) and copper (Cu) are among the common materials for busbar and battery tab manufacturing. A wide range of research shows that the laser welding of busbar to battery tabs is a very promising technique. It can enhance the battery module's safety and reliability owing to its unique properties. The desired strength, ductility, fatigue life as well as electrical resistivity are crucial to attain in laser welding of dissimilar materials aluminum and copper in busbar to battery tab in BEVs. Therefore, an adequate understating of the principal factors influencing the Al–Cu busbar to battery tabs joint properties are of prime importance. The current review paper provides information on laser welding and laser brazing of dissimilar Al–Cu with thin thicknesses. Also, the common defects, the effect of materials properties on laser joining, and laser-materials interaction during the laser welding process are discussed. Laser process parameters adjustment (e.g., laser power or speed), laser operational mode, and proper choice of materials (e.g., base metals, alloying elements, filler metals, etc.) may enhance the joint properties in terms of mechanical and electrical properties.

Keywords:

Citation: Ehsan Harati, Paul Kah. Laser welding of aluminum battery tab to variable Al/Cu busbars in Li-ion battery joint[J]. AIMS Materials Science, 2022, 9(6): 884-918. doi: 10.3934/matersci.2022053

Related Papers:

[1]	Xiao Xu, Hong Liao, Xu Yang . An automatic density peaks clustering based on a density-distance clustering index. AIMS Mathematics, 2023, 8(12): 28926-28950. doi: 10.3934/math.20231482
[2]	Fatimah Alshahrani, Wahiba Bouabsa, Ibrahim M. Almanjahie, Mohammed Kadi Attouch . $k$ NN local linear estimation of the conditional density and mode for functional spatial high dimensional data. AIMS Mathematics, 2023, 8(7): 15844-15875. doi: 10.3934/math.2023809
[3]	Shabbar Naqvi, Muhammad Salman, Muhammad Ehtisham, Muhammad Fazil, Masood Ur Rehman . On the neighbor-distinguishing in generalized Petersen graphs. AIMS Mathematics, 2021, 6(12): 13734-13745. doi: 10.3934/math.2021797
[4]	Yuming Chen, Jie Gao, Luan Teng . Distributed adaptive event-triggered control for general linear singular multi-agent systems. AIMS Mathematics, 2023, 8(7): 15536-15552. doi: 10.3934/math.2023792
[5]	Luca Bisconti, Paolo Maria Mariano . Global existence and regularity for the dynamics of viscous oriented fluids. AIMS Mathematics, 2020, 5(1): 79-95. doi: 10.3934/math.2020006
[6]	Hsin-Lun Li . Leader–follower dynamics: stability and consensus in a socially structured population. AIMS Mathematics, 2025, 10(2): 3652-3671. doi: 10.3934/math.2025169
[7]	Hongjie Li . Event-triggered bipartite consensus of multi-agent systems in signed networks. AIMS Mathematics, 2022, 7(4): 5499-5526. doi: 10.3934/math.2022305
[8]	Neelakandan Subramani, Abbas Mardani, Prakash Mohan, Arunodaya Raj Mishra, Ezhumalai P . A fuzzy logic and DEEC protocol-based clustering routing method for wireless sensor networks. AIMS Mathematics, 2023, 8(4): 8310-8331. doi: 10.3934/math.2023419
[9]	Yinwan Cheng, Chao Yang, Bing Yao, Yaqin Luo . Neighbor full sum distinguishing total coloring of Halin graphs. AIMS Mathematics, 2022, 7(4): 6959-6970. doi: 10.3934/math.2022386
[10]	Jin Cai, Shuangliang Tian, Lizhen Peng . On star and acyclic coloring of generalized lexicographic product of graphs. AIMS Mathematics, 2022, 7(8): 14270-14281. doi: 10.3934/math.2022786

Abstract

1. Introduction

Mathematics was initially utilized in biology in the $13^{th}$ century by Fibonacci, who developed the famed Fibonacci series to explain an increasing population. Daniel Bernoulli utilized mathematics to describe the impact on tiny shapes. Johannes Reinke coined the phrase "bio maths" in 1901. Biomathematics involves the theoretical examination of mathematical models to understand the principles governing the formation and functioning of biological systems.

Mathematical models are utilized to investigate specific questions related to the studied disease. For example, epidemiological models are important for predicting how infectious diseases spread, and help to control them by identifying important factors in the community. In this case, we want to analyze a particular model to understand how the COVID-19 virus behaves. This virus appeared in late 2019, and continues to be a global challenge. To gain a deeper insight into the underlying physical processes, we delve into the realm of fractional calculus. Previous literature has introduced various operators through the framework of fractional calculus ^[1,2]. In the field of C $^\eta$ -Calculus, Golmankhaneh et al. ^[3] provided an explanation of the Sumudu transform and Laplace transform. Additionally, in 2019, Goyal ^[4] proposed a fractional model that demonstrated the potential to manage the Lassa hemorrhagic fever disease.

Advancements in technology have led to significant progress in the field of epidemiology, enabling the examination of various infectious diseases for treatment, control, and cure ^[5]. It is essential to underscore that mathematical biology plays a pivotal role in investigating numerous diseases. Significant strides have been taken in the mathematical modeling of infectious diseases in recent decades, as evidenced by various studies ^[6,7]. Over the last thirty years, mathematical modeling has gained prominence in research, making substantial contributions to the development of effective public health strategies for disease control ^[8,9]. Mathematical models serve as invaluable tools for analyzing spatiotemporal patterns and the dynamic behavior of infections. Acknowledging their significance, researchers have approached the study of COVID-19 from diverse perspectives in the last three years ^[10,11].

Various methodologies have been employed by researchers in this field to devise successful techniques for managing this condition, with recent studies offering additional insights ^[12,13]. For example, a recent investigation utilized a mathematical model to evaluate the impacts of immunization in nursing homes ^[14]. Furthermore, researchers have explored mathematical modeling and effective intervention strategies for controlling the COVID-19 outbreak ^[15]. Additionally, some studies have delved into COVID-19 mathematical models using stochastic differential equations and environmental white noise ^[16].

In 2019, China experienced a notable outbreak of the coronavirus disease 2019 (COVID-19), prompting concerns about its potential to escalate into a worldwide pandemic ^[17]. Researchers from China, particularly Zhao et al., made significant contributions in addressing the challenges posed by COVID-19. This disease is attributed to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a viral infection. The initial verified case was documented in Wuhan, China, in December 2019 ^[18]. The infection rapidly spread worldwide, leading to the declaration of a COVID-19 pandemic.

Transmission occurs through various means, including respiratory droplets from coughing, sneezing, close contact, and touching contaminated surfaces. Key preventive measures include consistent mask usage, frequent hand washing, and maintaining safe interpersonal distances ^[19]. Effective interventions and real-time data play a crucial role in managing the coronavirus outbreak ^[20]. Prior studies have employed real-time analysis to comprehend the transmission of the virus among individuals, the severity of the disease, and the early stages of the pathogen, particularly in the initial week of the outbreak ^[21].

In December 2019, an outbreak of pneumonia cases was reported in Wuhan, initially with unidentified origins. Some cases were associated with exposure to wet markets and seafood. Chinese health authorities, in collaboration with the Chinese Center for Disease Control and Prevention (China CDC), initiated an investigation into the cause and spread of the disease on December 31, 2019 ^[22]. We conducted an analysis of temporal changes in the outbreak by examining the time interval between hospital admission dates and fatalities. Clinical studies on COVID-19 have indicated that symptoms typically manifest around 7 days after the onset of illness ^[23]. It is important to consider the duration between hospitalization and death in order to accurately assess the risk of mortality ^[24]. The information regarding the incubation period of COVID-19 and epidemiological data was sourced from publicly available records of confirmed cases ^[25].

An established method for the fractional-order model is elaborated upon in ^[26]. Recent contributions include various fractional models related to COVID-19, such as the analysis by Atangana and Khan focusing on the pandemic's impact on China ^[27]. Additionally, the COVID-19 model's dynamical aspects were explored using fuzzy Caputo and ABC derivatives, as demonstrated in ^[28]. A similar type of approach using fractional operator techniques are given in ^[29,30,31]. Different author's have investigated the transmission of different infectious disease like COVID-19 with symptomatic and asymptomatic effects in the community by using the fractal-fractional definition ^[32,33].

In light of the aforementioned significance, we aim is to address fundamental issues by concentrating on the distinctive challenges posed by the dynamics of COVID-19. To achieve this, we employ a model tailored to accurately capture the characteristics of COVID-19 dynamics and account for the limitations in our response to the pandemic. Initially, we examine the epidemic dynamics within a specific community characterized by a unique social pattern. For this analysis, we adopt a conventional SEIR design that accommodates prolonged incubation periods.

Here, the previous model is given in ^[34] as follows:

$\begin{eqnarray} \frac{dS}{dt}& = &\textbf{a}-\frac{\rho_{1} S I}{1+\gamma I}-(\delta+\rho)S,\\ \frac{dE}{dt}& = &\rho S-\delta E-\rho_{2}\alpha E I,\\ \frac{dI}{dt}& = &\frac{\rho_{1} S I}{1+\gamma I}\,\,+\rho_{2}\,\alpha\, E I-(\delta+\mu_{0}+\omega-\textbf{b})\,I,\\ \frac{dR}{dt}& = &\omega I-\delta R. \end{eqnarray}$

(1.1)

Initial conditions corresponds to the aforementioned system:

$\begin{equation*} S^{0}(t) \, = \, S_{0},\, \, E^{0}(t) \, = \, E_{0}, \, \,I^{0}(t) \, = \, I_{0},\,\, R^{0}(t)\, = \, R_{0}. \end{equation*}$

The primary goal of this research is to employ novel fractional derivatives within mathematical analysis and simulation to enhance the COVID-19 model. COVID-19 is a highly dangerous disease that presents a significant risk to human life. Verification of the existence and distinct characteristics of the solution system is undertaken, coupled with a qualitative assessment of the system. So, we introduce vaccination measures for low immune individuals. We developed new mathematical model by taking vaccination measures which helps to control COVID-19 early which we shall observe on simulation easily. The research involves confirming the presence of a solution system with unique characteristics and conducting a qualitative evaluation of this system. Furthermore, the fractal-fractional derivative is utilized to investigate the real-world behavior of the newly developed mathematical model. Finally, numerical simulations are used to reinforce and authenticate the biological findings.

Definition 1.1. If 0 $< \xi \leq$ 1 and 0 $< \lambda \leq$ 1, then the Riemann-Liouville operator for the fractal-fractional Operator (FFO) with a Mittag-Leffler (ML) kernel is defined as $U(t)$ ^[35].

$\begin{eqnarray} ^{FFM}_{0} D_{t}^{\xi, \lambda} U(t) = \frac{A B(\xi)}{1- \xi} \int^{t}_{0} E_{\xi} \frac{d \,U(\Omega)}{dt^{\lambda}}\left[- \frac{\xi}{1- \xi} (t-\Omega)^{\xi}\right] {d \Omega}\,\,, \end{eqnarray}$

involving $0 \; < \; \xi, \; \lambda \; \leq \; 1$ , and $AB(\xi) = 1 - \xi + \frac{\xi}{\Gamma(\xi)}.$

Therefore, the function $U(t)$ , which has an order of $(\xi, \lambda)$ and a Mittag-Leffler (ML) kernel, is given as follows:

$\begin{eqnarray} ^{FFM}_{0} D_{t}^{\xi, \lambda} U(t) = \frac{\lambda(1-\xi) t^{\lambda -1} U(t)}{A B(\xi)} + \frac{\xi \lambda}{A B(\xi)} \int^{t}_{0} {\Omega}^{\xi -1} (t-\Omega) U(\Omega) {d \Omega} \,\,. \end{eqnarray}$

2. Formulation of the SEVIR model

A newly developed model for SARS-COVID-19 includes the vaccinated effect, whereas the previous model used the SEIR framework. In this new model, we introduce a new variable called "Vaccinated." With the addition of Vaccinated, the new model is referred to as SEVIR, where "S" represents the Susceptible class, "E" represents the Exposed class, "V" represents the Vaccinated class, "I" represents the Infected class, and "R" represents the Recovered class.

We define several parameters in this model: "a" represents the recruitment rate, " $\delta$ " represents the death rate due to natural causes, " $\rho+\rho_{1}$ " represents the contact rate from the Susceptible class to the Exposed class, " $\rho_{2}$ " represents the vaccination rate, " $\alpha$ " represents the rate at which the infection is reducing due to vaccination effects, " $\mu_{0}$ " represents the infection death rate, " $\omega$ " represents the rate at which an individual recovers from vaccination and becomes recovered, " $\textbf{b}$ " represents the recruitment rate to the Vaccinated class, " $\phi$ " represents the contact rate from the Vaccinated class to the Infected class, and " $\psi$ " represents the recovery rate.

We want to investigate spread of the SEIR model for SARS-COVID-19 with the vaccinated effect.

So, the flow chart for newly developed model SEVIR is given as Figure 1.

Figure 1. The flow chart illustrates the newly developed model.

DownLoad: Full-Size Img PowerPoint

The model that was developed based on the generalized hypothesis with the vaccinated effect is presented as follows:

$\begin{eqnarray} \frac{d S}{dt}& = &\textbf{a}-(\delta+\rho+\rho_{1})S,\\ \frac{d E}{dt}& = &(\rho+\rho_{1})S-\delta E-\rho_{2}\alpha E V,\\ \frac{d V}{dt}& = &\rho_{2}\alpha E V-(\delta+\mu_{0}+\omega-\textbf{b}+\phi)V,\\ \frac{d I}{dt}& = &\phi V-(\mu_{0}+\delta+\psi)I,\\ \frac{d R}{dt}& = &\omega V+\psi I -\delta R. \end{eqnarray}$

(2.1)

The following are initial conditions linked with the described system:

$\begin{eqnarray} S^{0}(t) \, = \, S_{0},\, \, E^{0}(t) \, = \, E_{0}, \, \,V^{0}(t) \, = \, V_{0}, \, \,I^{0}(t) \, = \, I_{0},\,\, R^{0}(t)\, = \, R_{0}. \end{eqnarray}$

Using the fractal-fractional order (FFO) with a Mittag-Leffler (ML) definition, the above model becomes

$\begin{eqnarray} ^{FFM}_{0}D^{\xi,\lambda}_{t}S(t)& = &\textbf{a}-(\delta+\rho+\rho_{1})S,\\ ^{FFM}_{0}D^{\xi,\lambda}_{t}E(t)& = &(\rho+\rho_{1})S-\delta E-\rho_{2}\alpha E V,\\ ^{FFM}_{0}D^{\xi,\lambda}_{t}V(t)& = &\rho_{2}\alpha E V-(\delta+\mu_{0}+\omega-\textbf{b}+\phi)V,\\ ^{FFM}_{0}D^{\xi,\lambda}_{t}I(t)& = &\phi V-(\mu_{0}+\delta+\psi)I,\\ ^{FFM}_{0}D^{\xi,\lambda}_{t}R(t)& = &\omega V+\psi I-\delta R. \end{eqnarray}$

(2.2)

Here, $^{FFM}_{0}D^{\xi, \lambda}_{t}$ is the fractal-fractional operator with Mittag-Leffler (ML), where $0 < \xi \leq 1$ and $0 < \lambda \leq\; 1$ .

The following are initial conditions linked with the described system:

$\begin{eqnarray} S^{0}(t) \, = \, S_{0},\, \, E^{0}(t) \, = \, E_{0}, \, \,V^{0}(t) \, = \, V_{0}, \, \,I^{0}(t) \, = \, I_{0},\,\, R^{0}(t)\, = \, R_{0}. \end{eqnarray}$

Parameter descriptions are given in the following table.

Parameters	Representation	Reference
$\text{a}$	$\text{contact rate from the susceptible to the exposed class}$	^[34]
$\delta$	$\text{death rate due to natural causes}$	^[34]
$\rho+\rho_{1}$	$\text{contact rate from the susceptible to the exposed class}$	^[34]
$\alpha$	$\text{rate at which the infection is reducing due to vaccination}$	^[34]
$\rho_{2}$	$\text{vaccination rate}$	^[34]
$\mu_{0}$	death rate due to infection	^[34]
$\omega$	rate at which an individuals becomes recovered	^[34]
$\text{b}$	recruitment rate to the vaccinated class	^[34]
$\phi$	contact rate from the vaccinated to the infected class	Assumed
$\psi$	recovery rate	Assumed

2.1. Equilibrium point and reproductive number

For this model, the point of equilibrium without disease (disease free) is

$\begin{eqnarray} D_{1}(S,\,E,\,V,\,I,\,R) = \left(\frac{\textbf{a}}{\delta+\rho+\rho_{1}},\frac{\textbf{a}(\rho+\rho_{1})}{\delta(\delta+\rho+\rho_{1})},0,0,0\right), \end{eqnarray}$

as well as the endemic points of equilibrium $D_{2}(S^{*}, E^{*}, V^{*}, I^{*}, R^{*})$ , where

$\begin{eqnarray} S^{*}& = &\frac{\textbf{a}}{\delta+\rho+\rho_{1}},\\ E^{*}& = &\frac{-\textbf{b}+\delta+\phi+\omega+\mu_{0}}{\alpha\rho_{2}},\\ V^{*}& = &\frac{(\delta+\psi+\mu_{0})\textbf{A}-\textbf{b}\delta\rho_{1}+\delta^{2}\rho_{1}+\delta\phi\rho_{1}+\delta\omega\rho_{1}+\delta\mu_{0}\rho_{1}-\textbf{a}\alpha\rho\rho_{2}-\textbf{a}\alpha\rho_{1}\rho_{2}}{\alpha \phi(\textbf{b}-\delta-\phi-\omega-\mu_{0})(\delta+\rho+\rho_{1})\rho_{2}},\\ I^{*}& = &\phi\Bigg(\frac{\textbf{A}-\textbf{b}\delta\rho_{1}+\delta^{2}\rho_{1}+\delta\phi\rho_{1}+\delta\omega\rho_{1}+\delta\mu_{0}\rho_{1}-\textbf{a}\alpha\rho\rho_{2}-\textbf{a}\alpha\rho_{1}\rho_{2}}{\alpha(\textbf{b}-\delta-\phi-\omega-\mu_{0})(\delta+\rho+\rho_{1})\rho_{2}}\Bigg),\\ R^{*}& = &\big(\phi\psi+\delta\omega+\psi\omega+\omega\mu_{0}\big)\times\Bigg(\frac{\textbf{A}-\textbf{b}\delta\rho_{1}+\delta^{2}\rho_{1}+\delta\phi\rho_{1}+\delta\omega\rho_{1}+\delta\mu_{0}\rho_{1}-\textbf{a}\alpha\rho\rho_{2}-\textbf{a}\alpha\rho_{1}\rho_{2}}{\alpha \delta \phi(\delta+\psi+\mu_{0})(\textbf{b}-\delta-\phi-\omega-\mu_{0})(\delta+\rho+\rho_{1})\rho_{2}}\Bigg), \end{eqnarray}$

where

$\begin{eqnarray} \textbf{A} = -\textbf{b}\delta^{2}+\delta^{3}-\textbf{b}\delta\rho+\delta^{2}\rho+\delta^{2}\phi +\delta\rho\phi+\delta^{2}\omega+\delta\rho\omega+\delta^{2}\mu_{0}+\delta\rho\mu_{0}. \end{eqnarray}$

The reproductive number for the newly developed system by using the next generation method is

$\begin{eqnarray} R_{0} = \frac{\textbf{b}\delta(\delta+\rho+\rho_{1})+\textbf{a}\alpha\rho_{2}(\rho+\rho_{1})}{\delta(\delta+\rho+\rho_{1})(\delta+\phi+\omega+\mu_{0})}. \end{eqnarray}$

3. Bounded and positive solutions

In this section, we demonstrate the boundedness and positivity of the developed model.

Theorem 3.1. The considered initial condition is

$\begin{eqnarray} \{S_{0},E_{0}, V_{0},I_{0},R_{0}\}\subset\Upsilon, \end{eqnarray}$

and therefore the solutions $\{ S, E, V, I, R \}$ will be positive $\forall$ t $\geq$ 0.

Proof. We will begin the primary analysis to show the improved quality of the solutions. These solutions effectively address real-world issues and have positive outcomes. We will follow the methodology provided in references ^[36,37,38]. In this segment, we will examine the conditions required to ensure positive outcomes from the proposed model. To accomplish this, we will establish a standard.

$\begin{eqnarray} \parallel \beta \parallel_{\infty} = \sup\limits_{t\in D_{\beta}}\mid\beta(t)\mid, \end{eqnarray}$

where " $D_{\beta}$ " represents the $\beta$ domain. Now, we continue with $S(t)$ .

$\begin{eqnarray} ^{FFM}_{0}D^{\xi,\lambda}_{t}S(t)& = &\textbf{a}-(\delta+\rho+\rho_{1})S, \quad \forall t\geq0,\\ &\geq&-\left(\delta+\rho+\rho_{1}\right)S, \quad \forall t\geq0. \end{eqnarray}$

This yield

$\begin{eqnarray} S(t)\geq S(0)E_{\xi}\bigg[-\frac{c^{1-\lambda}{\xi}{(\delta+\rho+\rho_{1})}{t^{\xi}}}{{A B}(\xi)-(1-\xi)(\delta+\rho+\rho_{1})}\bigg], \quad \forall t\geq0, \end{eqnarray}$

where "c" represents the time element. This demonstrates that the $S(t)$ individuals must be positive $\forall \, \, t \geq$ 0. Now, we have the $E(t)$ individuals as follows:

$\begin{eqnarray} ^{FFM}_{0}D^{\xi,\lambda}_{t}E(t)& = &(\rho+\rho_{1})S-\delta E-\rho_{2}\alpha E V, \quad \forall \,\, t\geq0,\\ &\geq&-(\delta+\rho_{2}\alpha\mid V \mid)E, \quad \forall \,\,t\geq0,\\ &\geq&-(\delta+\rho_{2}\alpha\sup\limits_{t\in D_{V}}\mid V \mid )E, \quad \forall \,\,t\geq0,\\ &\geq&-(\delta+\rho_{2}\alpha\parallel V \parallel_{\infty})E, \quad \forall \,\,t\geq 0. \end{eqnarray}$

This yield

$\begin{eqnarray} E(t)\geq E(0)E_{\xi}\bigg[-\frac{c^{1-\lambda}{\xi}{(\delta+\rho_{2}\alpha\parallel V\parallel_{\infty})}{t^{\xi}}}{AB(\xi)-(1-\xi)(\delta+\rho_{2}\alpha\parallel V\parallel_{\infty})}\bigg], \quad \forall \,\,t\geq0, \end{eqnarray}$

where "c" represents the time element. This demonstrates that the $E(t)$ individuals must be positive $\forall \, \, t \geq$ 0. Now, we have the $V(t)$ individuals as follows:

$\begin{eqnarray} ^{FFM}_{0}D^{\xi,\lambda}_{t}V(t)& = &\rho_{2}\alpha E I-(\delta+\mu_{0}+\omega-\textbf{b}+\phi)V, \quad \forall t\geq0,\\ &\geq&-(-\rho_{2}\alpha \mid E \mid+\mu_{0}+\delta+\omega-\textbf{b}+\phi)V, \quad \forall t\geq0,\\ &\geq&-(-\rho_{2}\alpha \sup\limits_{t\in D_{E}}\mid E \mid+\mu_{0}+\delta+\omega-\textbf{b}+\phi)V, \quad \forall t\geq0,\\ &\geq&-(-\rho_{2}\alpha \parallel E \parallel_{\infty}+\mu_{0}+\delta+\omega-\textbf{b}+\phi)V, \quad \forall t\geq0. \end{eqnarray}$

This yield

$\begin{eqnarray} V(t)\geq V(0)E_{\xi}\bigg[-\frac{c^{1-\lambda{\xi}}{(-\rho_{2}\alpha \parallel E \parallel_{\infty}+\mu_{0}+\delta+\omega-\textbf{b}+\phi)}{t^{\xi}}}{A B(\xi)-(1-\xi)\left(-\rho_{2}\alpha\parallel E \parallel_{\infty}+\mu_{0}+\delta+\omega-\textbf{b}+\phi\right)}\bigg], \quad \forall t\geq0, \end{eqnarray}$

where "c" represents the time element. This demonstrates that the $V(t)$ individuals must be positive $\forall \, \, t \geq$ 0. Now, we have the $I(t)$ individuals as follows:

$\begin{eqnarray} ^{FFM}_{0}D^{\xi,\lambda}_{t}I(t)& = &\phi V-(\mu_{0}+\delta+\psi)I, \quad \forall t\geq0,\\ &\geq&-(\mu_{0}+\delta+\psi)I(t), \quad \forall t\geq0. \end{eqnarray}$

This yield

$\begin{eqnarray} I(t)\geq I(0)E_{\xi}\bigg[-\frac{c^{1-\lambda}{\xi}{(\mu_{0}+\delta+\psi)}{t^{\xi}}}{A B(\xi)-(1-\xi)(\mu_{0}+\delta+\psi)}\bigg], \quad \forall t\geq0, \end{eqnarray}$

where "c" represents the time element. This demonstrates that the $I(t)$ individuals must be positive $\forall \, \, t \geq$ 0. Now, we have the $R(t)$ individuals as follows:

$\begin{eqnarray} ^{FFM}_{0}D^{\xi,\lambda}_{t}R(t)& = &\omega V+\psi I-\delta R , \quad \forall \,t\,\geq\,0,\\ &\geq&-(\delta)R(t), \quad \forall \,t\,\geq\,0. \end{eqnarray}$

This yield

$\begin{eqnarray} R(t)\geq R(0){E_{\xi}}\bigg[-\frac{c^{1-\lambda}\xi{(\delta)}{t^{\xi}}}{{A B}(\xi)-(1-\xi)(\delta)}\bigg], \quad \forall t\geq0, \end{eqnarray}$

where "c" represents the time element. This demonstrates that the $R(t)$ individuals must be positive $\forall \, \, t \geq$ 0.

Theorem 3.2. Solutions of our developed model given in Eq (2.2) with positive initial values are all bounded.

Proof. The above theorem demonstrates that the solutions of our developed model must be positive $\forall \, \, t\, \geq\, 0$ , and the strategies are described in ^[39]. Because $X = S\, +\, E\, +\, V$ , then

$\begin{eqnarray} ^{FFM}_{0}D^{\xi,\lambda}_{t}X(t) = \textbf{a}-\delta X-(\mu_{0}+\omega-\textbf{b}+\phi)V. \end{eqnarray}$

We achieved as follows:

$\begin{eqnarray} \Psi_{p} = \{S,\quad E,\quad V \in R_{+}^{3}\mid S + V\leq X\} \quad \forall \,\,t\geq0. \end{eqnarray}$

Further we have $X_{\upsilon}$ = I+R. So, we have

$\begin{eqnarray} ^{FFM}_{0}D^{\xi,\lambda}_{t}X_{\upsilon}(t) = (\phi+\omega)V-\mu_{0}I-X_{\upsilon}\delta. \end{eqnarray}$

Upon solving the above equation and taking $t\rightarrow \infty$ , we get

$\begin{eqnarray} X_{\upsilon}\leq\frac{(\phi+\omega)V-\mu_{0}I}{\delta}. \end{eqnarray}$

Thus,

$\begin{eqnarray} \Psi_{\upsilon} = \bigg\{I,\,\,R\in R_{+}^{2}\mid X_{\upsilon}\leq\frac{(\phi+\omega)I-\mu_{0}Q}{\delta}\bigg\} \quad \forall \,\,t \geq 0. \end{eqnarray}$

The model's mathematical solutions (2.2) are confined to the region $\Psi$ .

$\begin{eqnarray} \Psi = \bigg\{S,\,\,E,\,\,V,\,\,I,\,\,R \in R_{+}^{5}\mid S+V\leq X,X_{\upsilon}\leq\frac{(\phi+\omega)V-\mu_{0}I}{\delta}\bigg\} \quad \forall \,\,t\geq0. \end{eqnarray}$

This demonstrates that for every $t \geq$ 0, all solutions remain positive, consistent with the provided initial conditions in the domain $\Psi$ .

Theorem 3.3. The proposed coronavirus model (2.2) in $R_{+}^{5}$ has positive invariant solutions, in addition to the initial conditions.

Proof. In this particular scenario, we applied the procedure described in ^[40]. We have

$\begin{eqnarray} ^{FFM}_{0}D^{\xi,\lambda}_{t}(S(t))_{S = 0}& = &\textbf{a} \quad \geq0,\\ ^{FFM}_{0}D^{\xi,\lambda}_{t}(E(t))_{E = 0}& = &(\rho+\rho_{1})S \quad \geq0,\\ ^{FFM}_{0}D^{\xi,\lambda}_{t}(V(t))_{V = 0}& = &\rho_{2}\alpha E V+\textbf{b}V \quad \geq0,\\ ^{FFM}_{0}D^{\xi,\lambda}_{t}(I(t))_{I = 0}& = &\phi V \quad \geq0,\\ ^{FFM}_{0}D^{\xi,\lambda}_{t}(R(t))_{R = 0}& = &\omega V+\psi I \quad \geq0. \end{eqnarray}$

(3.1)

If $(S_{0}, E_{0}, V_{0}, I_{0}, R_{0})$ $\in$ $R_{+}^{5}$ , then our obtained solution is unable to escape from the hyperplane, as stated in Eq (3.1). This proves that the $R_{+}^{5}$ domain is positive invariant.

4. Impact of global derivatives for uniqueness and existence of solutions

The Riemann-Stieltjes integral has been widely recognized in the literature as the most commonly used integral. If

$\begin{eqnarray} Y(x) = \int y(x){d}x, \end{eqnarray}$

then the Riemann-Stieltjes integral is given as follows:

$\begin{eqnarray} Y_{w}(x) = \int y(x){d}{w}(x), \end{eqnarray}$

where the $y(x)$ global derivative with respect to $w(x)$ is

$\begin{eqnarray} D_{w}y(x) = \lim\limits_{h\rightarrow0}\frac{y(x+h)-y(x)}{w(x+h)-w(x)}. \end{eqnarray}$

If the above function's numerator and denominator are differentiated, we get

$\begin{eqnarray} D_{w}y(x) = \frac{y^{'}(x)}{w^{'}(x)}, \end{eqnarray}$

assuming that $w^{'}(x)\, \, \neq$ 0, $\forall x\in D_{w^{'}}$ . Now, we will test the impact on the coronavirus by using the global derivative instead of the classical derivative:

$\begin{eqnarray} D_{w}S& = &\textbf{a}-(\delta+\rho+\rho_{1})S,\\ D_{w}E& = &(\rho+\rho_{1})S-\delta E-\rho_{2}\alpha E V,\\ D_{w}V& = &\rho_{2}\alpha E V-(\delta+\mu_{0}+\omega-\textbf{b}+\phi)V,\\ D_{w}I& = &\phi V-(\mu_{0}+\delta+\psi)I,\\ D_{w}R& = &\omega V+\psi I-\delta R. \end{eqnarray}$

For the sake of clean notation, we shall suppose that w is differentiable.

$\begin{eqnarray} S^{'}& = &w^{'}[\textbf{a}-(\delta+\rho+\rho_{1})S],\\ E^{'}& = &w^{'}[(\rho+\rho_{1})S-\delta E-\rho_{2}\alpha E V],\\ V^{'}& = &w^{'}[\rho_{2}\alpha E V-(\delta+\mu_{0}+\omega-\textbf{b}+\phi)V],\\ I^{'}& = &w^{'}[\phi V-(\mu_{0}+\delta+\psi)I],\\ R^{'}& = &w^{'}[\omega V+\psi I-\delta R]. \end{eqnarray}$

An appropriate choice of the function $w(t)$ will lead to a specific outcome. For instance, if $w(t) = t^{\alpha}$ , where $\alpha$ is a real number, we will observe fractal movement. We had to take action due to the circumstances that

$\begin{eqnarray} \parallel w^{'}\parallel_{\infty} = \sup\limits_{t\in D_{w^{'}}}\mid w^{'}(t)\mid < N. \end{eqnarray}$

The below example demonstrates the unique solution for the developed system:

$\begin{eqnarray} S^{'}& = &w^{'}[\textbf{a}-(\delta+\rho+\rho_{1})S] = Z_{1}(t, S, G),\\ E^{'}& = &w^{'}[(\rho+\rho_{1})S-\delta E-\rho_{2}\alpha E V] = Z_{2}(t, S, G),\\ V^{'}& = &w^{'}[\rho_{2}\alpha E V-(\delta+\mu_{0}+\omega- \textbf{b}+\phi)V] = Z_{3}(t, S, G), \\ I^{'}& = &w^{'}[\phi V-(\mu_{0}+\delta+\psi)I] = Z_{4}(t, S, G),\\ R^{'}& = &w^{'}[\omega V+\psi I-\delta R] = Z_{5}(t, S, G), \end{eqnarray}$

where $G = E, V, I, R$ .

We need to confirm the first two requirements as follows:

(1) $\mid Z(t, S, G)\mid^{2} < K(1+\mid S\mid^{2}$ ,

(2) $\forall \quad S_{1}, S_{2}$ , we have, $\parallel Z(t, S_{1}, G)-Z(t, S_{2}, G)\parallel^{2} < \bar{K}\parallel S_{1}-S_{2}\parallel_{\infty}^{2}$ .

Initially,

$\begin{eqnarray} \mid Z_{1}(t, S, G)\mid^{2}& = &\mid w^{'}\big[\textbf{a}-(\delta+\rho+\rho_{1})S\big]\mid^{2},\\ & = &\mid w^{'}\big[\textbf{a}+\left(-\delta-\rho-\rho_{1}\right)S\big]\mid^{2},\\ &\leq&2\mid w^{'}\mid^{2}\left(\textbf{a}^{2}+\mid(-\delta-\rho-\rho_{1})S\mid^{2}\right),\\ &\leq&2\sup\limits_{t\in D_{w^{'}}}\mid w^{'}\mid^{2}\textbf{a}^{2}+6\sup\limits_{t\in D_{w^{'}}}\mid w^{'}\mid^{2}\left(\rho^{2}+\delta^{2}+\rho_{1}\right)\mid S\mid^{2},\\ &\leq&2\parallel w^{'}\parallel^{2}_{\infty}\textbf{a}^{2}+6\parallel w^{'}\parallel^{2}_{\infty}\left(\rho^{2}+\delta^{2}+\rho_{1}^{2}\right)\mid S\mid^{2},\\ &\leq&2\parallel w^{'}\parallel^{2}_{\infty}\textbf{a}^{2}\left(1+\frac{3}{\textbf{a}^{2}}(\rho^{2}+\delta^{2}+\rho_{1}^{2})\mid S\mid^{2}\right),\\ &\leq&K_{1}(1+\mid S\mid^{2}), \end{eqnarray}$

under the condition

$\begin{eqnarray} \frac{3}{\textbf{a}^{2}}(\rho^{2}+\delta^{2}+\rho_{1}^{2}) < 1, \end{eqnarray}$

involving

$\begin{eqnarray} K_{1} = 2\parallel w^{'}\parallel^{2}_{\infty}\textbf{a}^{2}. \end{eqnarray}$

$\begin{eqnarray} \mid Z_{2}(t, S, G)\mid^{2}& = &\mid w^{'}\big[(\rho+\rho_{1})S-\delta E-\rho_{2}\alpha E V\big]\mid^{2},\\ & = &\mid w^{'}\big[(\rho+\rho_{1})S+(-\delta-\rho_{2}\alpha V)E\big]\mid^{2},\\ &\leq&2\mid w^{'}\mid^{2}\left(\mid(\rho+\rho_{1})S\mid^{2}+(-\delta-\rho_{2}\alpha V)E\mid^{2}\right),\\ &\leq&4\sup\limits_{t\in D_{w^{'}}}\mid w^{'}\mid^{2}(\rho^{2}+\rho_{1}^{2})\sup\limits_{t\in D_{S}}\mid S\mid^{2}+4\sup\limits_{t\in D_{w^{'}}}\mid w^{'}\mid^{2}(\delta^{2}+\rho^{2}_{2}\alpha^{2}\sup\limits_{t\in D_{V}}\mid V \mid^{2})\mid E\mid^{2},\\ &\leq&4\parallel w^{'}\parallel^{2}_{\infty}(\rho^{2}+\rho_{1}^{2})\parallel S\parallel^{2}_{\infty}+4\parallel w^{'}\parallel^{2}_{\infty}(\delta^{2}+\rho^{2}_{2}\alpha^{2}\parallel V \parallel^{2}_{\infty})\mid E\mid^{2},\\ &\leq&4\parallel w^{'}\parallel^{2}_{\infty}(\rho^{2}+\rho_{1}^{2})\parallel S\parallel^{2}_{\infty}\left(1+\frac{(\delta^{2}+\rho^{2}_{2}\alpha^{2}\parallel V \parallel^{2}_{\infty})\mid E\mid^{2}}{(\rho^{2}+\rho_{1}^{2})\parallel S\parallel^{2}_{\infty}}\right),\\ &\leq&K_{2}(1+\mid E\mid^{2}), \end{eqnarray}$

under the condition

$\begin{eqnarray} \frac{(\delta^{2}+\rho^{2}_{2}\alpha^{2}\parallel V \parallel^{2}_{\infty})}{(\rho^{2}+\rho_{1}^{2})\parallel S\parallel^{2}_{\infty}} < 1, \end{eqnarray}$

where

$\begin{eqnarray} K_{2} = 2\parallel w^{'}\parallel^{2}_{\infty}(\rho^{2}+\rho_{1}^{2})\parallel S\parallel^{2}_{\infty}. \end{eqnarray}$

$\begin{eqnarray} \mid Z_{3}(t, S, G)\mid^{2}& = &\mid w^{'}\big[\rho_{2}\alpha E V-(\delta+\mu_{0}+\omega-\textbf{b}+\phi)V\big]\mid^{2},\\ & = &\mid w^{'}\big[\rho_{2}\alpha E V+\textbf{b}V+(-\delta-\mu_{0}-\omega-\phi)V\big]\mid^{2},\\ &\leq&2\mid w^{'}\mid^{2}\left(\mid\rho_{2}\alpha E+\textbf{b}\mid^{2}+\mid(-\delta-\mu_{0}-\omega-\phi)\mid^{2}\right)\mid V \mid^{2},\\ &\leq&4\sup\limits_{t\in D_{w^{'}}}\mid w^{'}\mid^{2}\big[(\rho_{2}^{2}\alpha^{2}\sup\limits_{t\in D_{E}}\mid E\mid^{2}+\textbf{b}^{2})+2(\delta^{2}+\mu_{0}^{2}+\omega^{2}+\phi^{2})\big]\mid V \mid^{2},\\ &\leq&4\parallel\mid w^{'}\parallel^{2}_{\infty}\bigg[(\rho_{2}^{2}\alpha^{2}\parallel E\parallel^{2}_{\infty}+\textbf{b}^{2})+2(\delta^{2}+\mu_{0}^{2}+\omega^{2}+\phi^{2})\bigg]\mid V \mid^{2},\\ &\leq&4\parallel\mid w^{'}\parallel^{2}_{\infty}(\rho_{2}^{2}\alpha^{2}\parallel E\parallel^{2}_{\infty}+\textbf{b}^{2})\bigg[1+\frac{2(\delta^{2}+\mu_{0}^{2}+\omega^{2}+\phi^{2})}{\left(\rho_{2}^{2}\alpha^{2}\parallel E\parallel^{2}_{\infty}+\textbf{b}^{2}\right)}\bigg]\mid V \mid^{2},\\ &\leq&K_{3}(2\mid V \mid^{2}), \end{eqnarray}$

under the condition

$\begin{eqnarray} \frac{2(\delta^{2}+\mu_{0}^{2}+\omega^{2}+\phi^{2})}{\left(\rho_{2}^{2}\alpha^{2}\parallel E\parallel^{2}_{\infty}+\textbf{b}^{2}\right)}\leq1, \end{eqnarray}$

where

$\begin{eqnarray} K_{3} = 4\parallel\mid w^{'}\parallel^{2}_{\infty}(\rho_{2}^{2}\alpha^{2}\parallel E\parallel^{2}_{\infty}+\textbf{b}^{2}). \end{eqnarray}$

$\begin{eqnarray} \mid Z_{4}(t, S, G)\mid^{2}& = &\mid w^{'}\big[\phi V-(\mu_{0}+\delta+\psi)I\big]\mid^{2},\\ & = &\mid w^{'}\big[\phi V+(-\mu_{0}-\delta-\psi)I\big]\mid^{2},\\ &\leq&2\mid w^{'}\mid^{2}\left(\mid\phi V \mid^{2}+\mid(-\mu_{0}-\delta-\psi)I\mid^{2}\right),\\ &\leq&2\sup\limits_{t\in D_{w^{'}}}\mid w^{'}\mid^{2}\phi^{2}\sup\limits_{t\in D_{V}}\mid V \mid^{2}+6\sup\limits_{t\in D_{w^{'}}}\mid w^{'}\mid^{2}(\mu_{0}^{2}+\delta^{2}+\psi^{2})\mid I \mid^{2},\\ &\leq&2\parallel w^{'}\parallel^{2}_{\infty}\phi^{2}\parallel V \parallel^{2}_{\infty}+6\parallel w^{'}\parallel^{2}_{\infty}(\mu_{0}^{2}+\delta^{2}+\psi^{2})\mid I\mid^{2},\\ &\leq&2\parallel w^{'}\parallel^{2}_{\infty}\phi^{2}\parallel I\parallel^{2}_{\infty}\left(1+\frac{3(\mu_{0}^{2}+\delta^{2}+\psi^{2})\mid I \mid^{2}}{\phi^{2}\parallel I\parallel^{2}_{\infty}}\right),\\ &\leq&K_{4}(1+\mid I \mid^{2}), \end{eqnarray}$

under the condition

$\begin{eqnarray} \frac{3(\mu_{0}^{2}+\delta^{2}+\psi^{2})}{\phi^{2}\parallel V \parallel^{2}_{\infty}} < 1, \end{eqnarray}$

where

$\begin{eqnarray} K_{4} = 2\parallel w^{'}\parallel^{2}_{\infty}\phi^{2}\parallel V \parallel^{2}_{\infty}. \end{eqnarray}$

$\begin{eqnarray} \mid Z_{5}(t, S, G)\mid^{2}& = &\mid w^{'}\big[\omega V+\psi I-\delta R\big]\mid^{2},\\ & = &\mid w^{'}\big[(\omega V+\psi I)+(-\delta R)\big]\mid^{2},\\ &\leq&2\mid w^{'}\mid^{2}\left(\mid(\omega V+\psi I)\mid^{2}+\mid-\delta R\mid^{2}\right),\\ &\leq&4\sup\limits_{t\in D_{w^{'}}}\mid w^{'}\mid^{2}\left(\omega^{2}\sup\limits_{t\in D_{V}}\mid V \mid^{2}+\psi\sup\limits_{t\in D_{I}}\mid I \mid^{2}\right)+2\sup\limits_{t\in D_{w^{'}}}\mid w^{'}\mid^{2}\delta^{2}\mid R\mid^{2},\\ &\leq&4\parallel w^{'}\parallel^{2}_{\infty}\left(\omega^{2}\parallel V \parallel^{2}_{\infty}+\psi\parallel I \parallel^{2}_{\infty}\right)+2\parallel w^{'}\parallel^{2}_{\infty}\delta^{2}\mid R\mid^{2}),\\ &\leq&4\parallel w^{'}\parallel^{2}_{\infty}(\omega^{2}\parallel V \parallel^{2}_{\infty}+\psi^{2}\parallel I\parallel^{2}_{\infty})\left(1+\frac{\delta^{2}\mid R\mid^{2}}{2(\omega^{2}\parallel V \parallel^{2}_{\infty}+\psi^{2}\parallel I \parallel^{2}_{\infty})}\right),\\ &\leq&K_{5}(1+\mid R\mid^{2}), \end{eqnarray}$

under the condition

$\begin{eqnarray} \frac{\delta^{2}}{2(\omega^{2}\parallel V \parallel^{2}_{\infty}+\psi^{2}\parallel I \parallel^{2}_{\infty})} < 1, \end{eqnarray}$

where

$\begin{eqnarray} K_{5} = 4\parallel w^{'}\parallel^{2}_{\infty}(\omega^{2}\parallel V \parallel^{2}_{\infty}+\psi^{2}\parallel I \parallel^{2}_{\infty}). \end{eqnarray}$

Hence the linear growth condition is satisfied.

Further, we validate the Lipschitz condition.

$\begin{eqnarray} \mid Z_{1}(t, S_{1}, G)-Z_{1}(t, S_{2}, G)\mid^{2}& = &\mid w^{'}(-\delta-\rho-\rho_{1})(S_{1}-S_{2})\mid^{2},\\ \mid Z_{1}(t, S_{1}, G)-Z_{1}(t, S_{2}, G)\mid^{2}&\leq&\mid w^{'}\mid^{2}(3\delta^{2}+3\rho^{2}+3\rho_{1}^{2})\mid S_{1}-S_{2}\mid^{2},\\ \sup\limits_{t\in D_{S}}\mid Z_{1}(t, S_{1}, G)-Z_{1}(t, S_{2}, G)\mid^{2}&\leq&\sup\limits_{t\in D_{w^{'}}}\mid w^{'}\mid^{2}\left(3\delta^{2}+3\rho^{2}+3\rho_{1}^{2}\right)\sup\limits_{t\in D_{S}}\mid S_{1}-S_{2}\mid^{2},\\ \parallel Z_{1}(t, S_{1}, G)-Z_{1}(t, S_{2}, G)\parallel_{\infty}^{2}&\leq&\parallel w^{'}\parallel^{2}_{\infty}(3\delta^{2}+3\rho^{2}+3\rho_{1}^{2})\parallel S_{1}-S_{2}\parallel^{2}_{\infty},\\ \parallel Z_{1}(t, S_{1}, G)-Z_{1}(t, S_{2}, G)\parallel_{\infty}^{2}&\leq&\bar{K_{1}}\parallel S_{1}-S_{2}\parallel^{2}_{\infty}, \end{eqnarray}$

where

$\begin{eqnarray} \bar{K_{1}} = \parallel w^{'}\parallel^{2}_{\infty}\left(3\delta^{2}+3\rho^{2}+3\rho_{1}^{2}\right). \end{eqnarray}$

$\begin{eqnarray} \mid Z_{2}(t,\,\,S,\,\,E_{1},\,\,V,\,\,I,\,\,R)-Z_{2}(t,\,\,S,\,\,E_{2},\,\,V,\,\,I,\,\,R)\mid^{2}& = &\mid w^{'}(-\delta-\rho_{2}\alpha V)(E_{1}-E_{2})\mid^{2},\\ \mid Z_{2}(t,\,\,S,\,\,E_{1},\,\,V,\,\,I,\,\,R)-Z_{2}(t,\,\,S,\,\,E_{2},\,\,V,\,\,I,\,\,R)\mid^{2}&\leq&\mid w^{'}\mid^{2}(2\delta^{2}+2\rho_{2}^{2}\alpha^{2}\mid V \mid^{2})\mid E_{1}-E_{2}\mid^{2},\\ \sup\limits_{t\in D_{E}}\mid Z_{2}(t,\,\,S,\,\,E_{1},\,\,V,\,\,I,\,\,R)-Z_{2}(t,\,\,S,\,\,E_{2},\,\,V,\,\,I,\,\,R)\mid^{2}&\leq&\sup\limits_{t\in D_{w^{'}}}\mid w^{'}\mid^{2}(2\delta^{2}+2\rho_{2}^{2}\alpha^{2}\sup\limits_{t\in D_{V}}\mid V \mid^{2})\sup\limits_{t\in D_{E}}\\&&\times\mid E_{1}-E_{2}\mid^{2},\\ \parallel Z_{2}(t,\,\,S,\,\,E_{1},\,\,V,\,\,I,\,\,R)-Z_{2}(t,\,\,S,\,\,E_{2},\,\,V,\,\,I,\,\,R)\parallel_{\infty}^{2}&\leq&\parallel w^{'}\parallel^{2}_{\infty}(2\delta^{2}+2\rho_{2}^{2}\alpha^{2}\parallel V \parallel^{2}_{\infty}),\parallel E_{1}-E_{2}\parallel^{2}_{\infty},\\ \parallel Z_{2}(t,\,\,S,\,\,E_{1},\,\,V,\,\,I,\,\,R)-Z_{2}(t,\,\,S,\,\,E_{2},\,\,V,\,\,I,\,\,R)\parallel_{\infty}^{2}&\leq&\bar{K_{2}}\parallel E_{1}-E_{2}\parallel^{2}_{\infty}, \end{eqnarray}$

where

$\begin{eqnarray} \bar{K_{2}} = \parallel w^{'}\parallel^{2}_{\infty}(2\delta^{2}+2\rho_{2}^{2}\alpha^{2}\parallel V \parallel^{2}_{\infty}). \end{eqnarray}$

$\begin{eqnarray} &&\mid Z_{3}(t,\,\,S,\,\,E,\,\,V_{1},\,\,I,\,\,R)-Z_{3}(t,\,\,S,\,\,E,\,\,V_{2},\,\,I,\,\,R)\mid^{2}\\& = &\mid w^{'}\left(\rho_{2}\alpha E+\textbf{b}+(-\delta-\mu_{0}-\omega-\phi)\right)(V_{1}-V_{2})\mid^{2}, \end{eqnarray}$

$\begin{eqnarray} &&\mid Z_{3}(t,\,\,S,\,\,E,\,\,V_{1},\,\,I,\,\,R)-Z_{3}(t,\,\,S,\,\,E,\,\,V_{2},\,\,I,\,\,R)\mid^{2}\\&\leq&2\mid w^{'}\mid^{2}\left(\mid\rho_{2}\alpha E+\textbf{b}\mid^{2}+\mid(-\delta-\mu_{0}-\omega-\phi)\mid^{2}\right)\\&&\times\mid V_{1}-V_{2}\mid^{2}, \end{eqnarray}$

$\begin{eqnarray} \sup\limits_{t\in D_{V}}\mid Z_{3}(t,\,\,S,\,\,E,\,\,V_{1},\,\,I,\,\,R)-Z_{3}(t,\,\,S,\,\,E,\,\,V_{2},\,\,I,\,\,R)\mid^{2}&\leq&4\sup\limits_{t\in D_{w^{'}}}\mid w^{'}\mid^{2}\bigg(\rho_{2}^{2}\alpha^{2}\sup\limits_{t\in D_{E}}\mid E\mid^{2}+\textbf{b}^{2}\\&&+2(\delta^{2}+\mu_{0}^{2}+\omega^{2}+\phi^{2})\bigg)\sup\limits_{t\in D_{V}}\mid V_{1}-V_{2}\mid^{2},\\ \parallel Z_{3}(t,\,\,S,\,\,E,\,\,V_{1},\,\,I,\,\,R)-Z_{3}(t,\,\,S,\,\,E,\,\,V_{2},\,\,I,\,\,R)\parallel_{\infty}^{2}&\leq&4\parallel w^{'}\parallel^{2}_{\infty}\bigg(\rho_{2}^{2}\alpha^{2}\parallel E\parallel^{2}_{\infty}+\textbf{b}^{2}\\&&+2(\delta^{2}+\mu_{0}^{2}+\omega^{2}+\phi^{2})\bigg)\parallel V_{1}-V_{2}\parallel^{2}_{\infty},\\ \parallel Z_{3}(t,\,\,S,\,\,E,\,\,V_{1},\,\,I,\,\,R)-Z_{3}(t,\,\,S,\,\,E,\,\,V_{2},\,\,I,\,\,R)\parallel_{\infty}^{2}&\leq&\bar{K_{3}}\parallel V_{1}-V_{2}\parallel^{2}_{\infty}, \end{eqnarray}$

where

$\begin{eqnarray} \bar{K_{3}} = 4\parallel w^{'}\parallel^{2}_{\infty}\bigg(\rho_{2}^{2}\alpha^{2}\parallel E\parallel^{2}_{\infty}+\textbf{b}^{2}+2(\delta^{2}+\mu_{0}^{2}+\omega^{2}+\phi^{2})\bigg). \end{eqnarray}$

$\begin{eqnarray} \mid Z_{4}(t,\,\,S,\,\,E,\,\,V,\,\,I_{1},\,\,R)-Z_{4}(t,\,\,S,\,\,E,\,\,V,\,\,I_{2},\,\,R)\mid^{2}& = &\mid w^{'}(-\mu_{0}-\delta-\psi)\big](Q_{1}-Q_{2})\mid^{2},\\ \mid Z_{4}(t,\,\,S,\,\,E,\,\,V,\,\,I_{1},\,\,R)-Z_{4}(t,\,\,S,\,\,E,\,\,V,\,\,I_{2},\,\,R)\mid^{2}& = &\mid w^{'}\mid^{2}(3\mu_{0}^{2}+3\delta^{2}+3\psi^{2})\mid I_{1}-I_{2}\mid^{2},\\ \sup\limits_{t\in D_{I}}Z_{4}(t,\,\,S,\,\,E,\,\,V,\,\,I_{1},\,\,R)-Z_{4}(t,\,\,S,\,\,E,\,\,V,\,\,I_{2},\,\,R)\mid^{2}& = &\sup\limits_{t\in D_{w^{'}}}\mid w^{'}\mid^{2}(3\mu_{0}^{2}+3\delta^{2}+3\psi^{2})\sup\limits_{t\in D_{I}}\mid I_{1}-I_{2}\mid^{2},\\ \parallel Z_{4}(t,\,\,S,\,\,E,\,\,V,\,\,I_{1},\,\,R)-Z_{4}(t,\,\,S,\,\,E,\,\,V,\,\,I_{2},\,\,R)\parallel_{\infty}^{2}&\leq&\parallel w^{'}\parallel^{2}_{\infty}(3\mu_{0}^{2}+3\delta^{2}+3\psi^{2})\parallel I_{1}-I_{2}\parallel^{2}_{\infty},\\ \parallel Z_{4}(t,\,\,S,\,\,E,\,\,V,\,\,I_{1},\,\,R)-Z_{4}(t,\,\,S,\,\,E,\,\,V,\,\,I_{2},\,\,R)\parallel_{\infty}^{2}&\leq&\bar{K_{4}}\parallel I_{1}-I_{2}\parallel^{2}_{\infty}, \end{eqnarray}$

where

$\begin{eqnarray} \bar{K_{4}} = \parallel w^{'}\parallel^{2}_{\infty}(3\mu_{0}^{2}+3\delta^{2}+3\psi^{2}). \end{eqnarray}$

$\begin{eqnarray} \mid Z_{5}(t,\,\,S,\,\,E,\,\,V,\,\,I,\,\,R_{1})-Z_{5}(t,\,\,S,\,\,E,\,\,V,\,\,I,\,\,R_{2})\mid^{2}& = &\mid w^{'}(-\delta)(R_{1}-R_{2})\mid^{2},\\ \mid Z_{5}(t,\,\,S,\,\,E,\,\,V,\,\,I,\,\,R_{1})-Z_{5}(t,\,\,S,\,\,E,\,\,V,\,\,I,\,\,R_{2})\mid^{2}&\leq&\mid w^{'}\mid^{2}\delta^{2}\mid(R_{1}-R_{2})\mid^{2},\\ \sup\limits_{t\in D_{R}}\mid Z_{5}(t,\,\,S,\,\,E,\,\,V,\,\,I,\,\,R_{1})-Z_{5}(t,\,\,S,\,\,E,\,\,V,\,\,I,\,\,R_{2})\mid^{2}&\leq&\sup\limits_{t\in D_{w^{'}}}\mid w^{'}\mid^{2}\delta^{2}\sup\limits_{t\in D_{R}}\mid R_{1}-R_{2}\mid^{2},\\ \parallel Z_{5}(t,\,\,S,\,\,E,\,\,V,\,\,I,\,\,R_{1})-Z_{5}(t,\,\,S,\,\,E,\,\,V,\,\,I,\,\,R_{2})\parallel_{\infty}^{2}&\leq&\parallel w^{'}\parallel^{2}_{\infty}\delta^{2}\parallel R_{1}-R_{2}\parallel_{\infty}^{2},\\ \parallel Z_{5}(t,\,\,S,\,\,E,\,\,V,\,\,I,\,\,R_{1})-Z_{5}(t,\,\,S,\,\,E,\,\,V,\,\,I,\,\,R_{2})\parallel_{\infty}^{2}&\leq&\bar{K_{5}}\parallel R_{1}-R_{2}\parallel^{2}_{\infty}, \end{eqnarray}$

involving

$\begin{eqnarray} \bar{K_{5}} = \parallel w^{'}\parallel^{2}_{\infty}\delta^{2}. \end{eqnarray}$

Then, given the condition, system (2.2) has a particular solution.

$\begin{eqnarray} \max\Big[\frac{3}{\textbf{a}^{2}}(\rho^{2}+\delta^{2}+\rho_{1}^{2}),\frac{(\delta^{2}+\rho^{2}_{2}\alpha^{2}\parallel V \parallel^{2}_{\infty})}{(\rho^{2}+\rho_{1}^{2})\parallel S\parallel^{2}_{\infty}},\frac{2(\delta^{2}+\mu_{0}^{2}+\omega^{2}+\phi^{2})}{\left(\rho_{2}^{2}\alpha^{2}\parallel E\parallel^{2}_{\infty}+\textbf{b}^{2}\right)},\frac{3(\mu_{0}^{2}+\delta^{2}+\psi^{2})}{\phi^{2}\parallel V \parallel^{2}_{\infty}},\\ \frac{\delta^{2}}{2(\omega^{2}\parallel V \parallel^{2}_{\infty}+\psi^{2}\parallel I \parallel^{2}_{\infty})}\Big] < 1. \end{eqnarray}$

5. Global stability for developed system

We use Lyapunov's approach and LaSalle's concept of invariance to analyze global stability and determine the conditions for eliminating diseases.

5.1. Lyapunov's first derivative

Theorem 5.1. ^[41] When the reproductive number $R_{0} >$ 1, the endemic equilibrium points of the SEVIR model are globally asymptotically stable.

Proof. The Lyapunov function can be expressed in the following manner:

$\begin{eqnarray} L(S^{*},E^{*},V^{*},I^{*},R^{*})& = &\left(S-S^{*}-S^{*}\log\frac{S}{S^{*}}\right)+\left(E-E^{*}-E^{*}\log\frac{E}{E^{*}}\right) +\left(V-V^{*}-V^{*}\log\frac{V}{V^{*}}\right)\\ &&+\left(I-I^{*}-I^{*}\log\frac{I}{I^{*}}\right)+\left(R-R^{*}-R^{*}\log\frac{R}{R^{*}}\right). \end{eqnarray}$

By applying a derivative on both sides,

$\begin{eqnarray} D^{\xi,\lambda}_{t} L = \dot{L} = \left(\frac{S-S^{*}}{S}\right)D^{\xi,\lambda}_{t}S+\left(\frac{E-E^{*}}{E}\right)D^{\xi,\lambda}_{t}E+\left(\frac{V-V^{*}}{V}\right)D^{\xi,\lambda}_{t}V+\left(\frac{I-I^{*}}{I}\right)D^{\xi,\lambda}_{t}I+\left(\frac{R-R^{*}}{R}\right)D^{\xi,\lambda}_{t}R, \end{eqnarray}$

we get

$\begin{eqnarray} D^{\xi,\lambda}_{t} L & = &\left(\frac{S-S^{*}}{S}\right)\left(\textbf{a}-(\delta+\rho+\rho_{1})S\right) +\left(\frac{E-E^{*}}{E}\right)\left((\rho+\rho_{1})S-\delta E-\rho_{2}\alpha E V\right)+\left(\frac{V-V^{*}}{V}\right)\\&&\times\left(\rho_{2}\alpha E V -(\delta+\mu_{0}+\omega-\textbf{b}+\phi)V\right)+\left(\frac{I-I^{*}}{I}\right)\left(\phi V-(\mu_{0}+\delta+\psi)I\right)+\left(\frac{R-R^{*}}{R}\right)\\&&\times\left(\omega V+\psi I-\delta R\right), \end{eqnarray}$

and setting $S = S-S^{*}, E = E-E^{*}, V = V-V^{*}, I = I-I^{*}$ and $R = R-R^{*}$ results in

$\begin{eqnarray} D^{\xi,\lambda}_{t} L& = &\textbf{a}-\textbf{a}\frac{S^{*}}{S}-(\delta+\rho+\rho_{1})\frac{(S-S^{*})^{2}}{S} +(\rho+\rho_{1})S-(\rho+\rho_{1})S^{*}-(\rho+\rho_{1})S\frac{E^{*}}{E}+(\rho+\rho_{1})S^{*}\frac{E^{*}}{E} \\&&-\delta\frac{(E-E^{*})^{2}}{E}-\rho_{2}\alpha V\frac{(E-E^{*})^{2}}{E}+\rho_{2}\alpha V^{*}\frac{(E-E^{*})^{2}}{E}+\rho_{2}\alpha E\frac{(V-V^{*})^{2}}{V}-\rho_{2}\alpha E^{*}\frac{(V-V^{*})^{2}}{V} \\&&+\textbf{b}\frac{(V-V^{*})^{2}}{V}-(\delta+\mu_{0}+\omega+\phi)\frac{(V-V^{*})^{2}}{V}+\phi V-\phi V^{*}-\phi V\frac{I^{*}}{I}+\phi V^{*}\frac{I^{*}}{I}\\&&-(\mu_{0}+\delta+\psi)\frac{(I-I^{*})^{2}}{I}+\omega V -\omega V^{*}-\omega V\frac{R^{*}}{R}+\omega V^{*}\frac{R^{*}}{R}+\psi I-\psi I^{*}\\&&-\psi I\frac{R^{*}}{R}+\psi I^{*}\frac{R^{*}}{R}-\delta\frac{(R-R^{*})^{2}}{R}. \end{eqnarray}$

We can write $D^{\xi, \lambda}_{t} L = \Sigma-\Omega$ , where

$\begin{eqnarray} \Sigma& = &\textbf{a}+(\rho+\rho_{1})S+(\rho+\rho_{1})S^{*}\frac{E^{*}}{E}+\rho_{2}\alpha V^{*}\frac{(E-E^{*})^{2}}{E}+\rho_{2}\alpha E\frac{(V-V^{*})^{2}}{V}+\textbf{b}\frac{(V-V^{*})^{2}}{V} +\phi V\\&&+\phi V^{*}\frac{I^{*}}{I}+\omega V+\omega V^{*}\frac{R^{*}}{R}+\psi I+\psi I^{*}\frac{R^{*}}{R}, \end{eqnarray}$

and

$\begin{eqnarray} \Omega& = &\textbf{a}\frac{S^{*}}{S}+(\delta+\rho+\rho_{1})\frac{(S-S^{*})^{2}}{S}+(\rho+\rho_{1})S^{*} +(\rho+\rho_{1})S\frac{E^{*}}{E}+\delta\frac{(E-E^{*})^{2}}{E}+\rho_{2}\alpha V\frac{(E-E^{*})^{2}}{E}\\&&+\rho_{2}\alpha E^{*}\frac{(V-V^{*})^{2}}{V}+(\delta+\mu_{0}+\omega+\phi)\frac{(V-V^{*})^{2}}{V}+\phi V^{*}+\phi V\frac{I^{*}}{I}+(\mu_{0}+\delta+\psi)\frac{(I-I^{*})^{2}}{I} +\omega V^{*}\\&&+\omega V\frac{R^{*}}{R}+\psi V^{*}+\psi I\frac{R^{*}}{R}+\delta\frac{(R-R^{*})^{2}}{R}. \end{eqnarray}$

We conclude that if $\Sigma < \Omega$ , this yields $D^{\xi, \lambda}_{t} L < 0$ , however when $S = S^{*}, E = E^{*}, V = V^{*}, I = I^{*}$ and $R = R^{*}$ , $\Sigma-\Omega = 0 \quad \Rightarrow D^{\xi, \lambda}_{t} L = 0$ .

We can observe that $\left\{(S^{*}, E^{*}, V^{*}, I^{*}, R^{*})\in\Gamma : D^{\xi, \lambda}_{t} L = 0\right\}$ represents the point $D_2$ for the developed model.

According to Lasalles' concept of invariance, $D_2$ is globally uniformly stable in $\Gamma$ if $\Sigma-\Omega = 0$ .

6. Solutions by fractal-fractional operator

Now, we will develop a solution using a numerical approach for our newly developed model given in Eq (2.2). We use the ML kernel in the current scenario instead of the classical derivative operator.

Furthermore, we will use the variable order version.

$\begin{eqnarray} ^{FFM}_{0}D^{\xi,\lambda}_{t} S(t)& = &\textbf{a}-(\delta+\rho+\rho_{1})S,\\ ^{FFM}_{0}D^{\xi,\lambda}_{t} E(t)& = &(\rho+\rho_{1})S-\delta E-\rho_{2}\alpha E V,\\ ^{FFM}_{0}D^{\xi,\lambda}_{t} V(t)& = &\rho_{2}\alpha E V-(\delta+\mu_{0}+\omega-\textbf{b}+\phi)V,\\ ^{FFM}_{0}D^{\xi,\lambda}_{t} I(t)& = &\phi V-(\mu_{0}+\delta+\psi)I,\\ ^{FFM}_{0}D^{\xi,\lambda}_{t} R(t)& = &\omega V+\psi I-\delta R. \end{eqnarray}$

For clarity, we express the above equation as follows:

$\begin{eqnarray} ^{FFM}_{0}D^{\xi,\lambda}_{t} S(t)& = &S_{1}(t, S, G),\\ ^{FFM}_{0}D^{\xi,\lambda}_{t} E(t)& = &E_{1}(t, S, G),\\ ^{FFM}_{0}D^{\xi,\lambda}_{t} V(t)& = &V_{1}(t, S, G),\\ ^{FFM}_{0}D^{\xi,\lambda}_{t} I(t)& = &I_{1}(t, S, G),\\ ^{FFM}_{0}D^{\xi,\lambda}_{t} R(t)& = &R_{1}(t, S, G). \end{eqnarray}$

Where

$\begin{eqnarray} S_{1}(t, S, G)& = &\textbf{a}-(\delta+\rho+\rho_{1})S,\\ E_{1}(t, S, G)& = &(\rho+\rho_{1})S-\delta E-\rho_{2}\alpha E V,\\ V_{1}(t, S, G)& = &\rho_{2}\alpha E V-(\delta+\mu_{0}+\omega-\textbf{b}+\phi)V,\\ I_{1}(t, S, G)& = &\phi V-(\mu_{0}+\delta+\psi)I,\\ R_{1}(t, S, G)& = &\omega V+\psi I-\delta R. \end{eqnarray}$

After using the fractal-fractional integral with the ML kernel, we obtain the following results:

$\begin{eqnarray} S(t_{\eta}+1)& = &\frac{\lambda(1-\xi)}{A B(\xi)}t_{\eta}^{\lambda-1}S_{1}\left(t_{\eta},S(t_{\eta}),G(t_{\eta})\right)\\&&+ \frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}\int_{t_{\nu}}^{t_{\nu+1}}S_{1}(t,S,G)\tau^{\xi-1}(t_{\eta+1} -\tau)^{\xi-1}{d}\tau,\\ E(t_{\eta}+1)& = &\frac{\lambda(1-\xi)}{A B(\xi)}t_{\eta}^{\lambda-1}E_{1}\left(t_{\eta},S(t_{\eta}),G(t_{\eta})\right) \\&&+\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}\int_{t_{\nu}}^{t_{\nu+1}}E_{1}(t,S,G)\tau^{\xi-1}(t_{\eta+1} -\tau)^{\xi-1}{d}\tau,\\ V(t_{\eta}+1)& = &\frac{\lambda(1-\xi)}{A B(\xi)}t_{\eta}^{\lambda-1}V_{1}\left(t_{\eta},S(t_{\eta}),G(t_{\eta})\right) \\&&+\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}\int_{t_{\nu}}^{t_{\nu+1}}V_{1}(t,S,G)\tau^{\xi-1}(t_{\eta+1} -\tau)^{\xi-1}{d}\tau,\\ I(t_{\eta}+1)& = &\frac{\lambda(1-\xi)}{A B(\xi)}t_{\eta}^{\lambda-1}I_{1}\left(t_{\eta},S(t_{\eta}),G(t_{\eta})\right) \\&&+\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}\int_{t_{\nu}}^{t_{\nu+1}}I_{1}(t,S,G)\tau^{\xi-1}(t_{\eta+1} -\tau)^{\xi-1}{d}\tau,\\ R(t_{\eta}+1)& = &\frac{\lambda(1-\xi)}{A B(\xi)}t_{\eta}^{\lambda-1}R_{1}\left(t_{\eta},S(t_{\eta}),G(t_{\eta})\right) \\&&+\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}\int_{t_{\nu}}^{t_{\nu+1}}R_{1}(t,S,G)\tau^{\xi-1}(t_{\eta+1} -\tau)^{\xi-1}{d}\tau, \end{eqnarray}$

(6.1)

where $G(t_{\eta}) = E(t_{\eta}), \, \, V(t_{\eta}), \, \, I(t_{\eta}), \, \, R(t_{\eta})$ .

Remember that the Newton polynomial can be obtained by using the Newton interpolation formula.

$\begin{eqnarray} N(t, S, G)&\simeq&N(t_{\eta-2},S_{\eta-2},G_{\eta-2})+\frac{1}{\Delta t}\Big[N(t_{\eta-1},S_{\eta-1},G_{\eta-1})\\ &-&N(t_{\eta-2},S_{\eta-2},G_{\eta-2})\Big]\left(\tau-t_{\eta-2}\right)+\frac{1}{2\Delta t^{2}}\Big[N(t_{\eta},S_{\eta},E_{\eta},I_{\eta},Q_{\eta},R_{\eta})\\ &-&2N(t_{\eta-1},S_{\eta-1},G_{\eta-1})-N(t_{\eta-2},S_{\eta-2},G_{\eta-2})\Big] \left(\tau-t_{\eta-2}\right)\left(\tau-t_{\eta-1}\right), \end{eqnarray}$

where $G_{\eta-2} = E_{\eta-2}, V_{\eta-2}, I_{\eta-2}, R_{\eta-2}$ , $G_{\eta-1} = E_{\eta-1}, V_{\eta-1}, I_{\eta-1}, R_{\eta-1}$ .

When we substitute the Newton polynomial into the system of Eqs (6.1), we obtain the following:

$\begin{eqnarray} S^{\eta+1}& = &\frac{\lambda(1-\xi)}{A B(\xi)}t_{\eta}^{\lambda-1}S_{1}(t_{\eta},S(t_{\eta}),G(t_{\eta}))+\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}S_{1} (t_{\nu-2},S^{\nu-2},G^{\nu-2})\\&&\times t^{\lambda-1}_{\nu-2} \int_{t_{\nu}}^{t_{\nu+1}}(t_{\eta+1}-\tau)^{\xi-1}{d}\tau +\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}\frac{1}{\Delta t}\Big[t^{\lambda-1}_{\nu-1}S_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1}) \\&&-t^{\lambda-1}_{\nu-2} S_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \int_{t_{\nu}}^{t_{\nu+1}}(\tau-t_{\nu-2})(t_{\eta+1}-\tau)^{\xi-1}{d}\tau\\&&+\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}\frac{1}{2\Delta t^{2}}\Big[t^{\lambda-1}_{\nu}S_{1}(t_{\nu},S^{\nu},G^{\nu}) -2t^{\lambda-1}_{\nu-1}S_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\&&+t^{\lambda-1}_{\nu-2} S_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \int_{t_{\nu}}^{t_{\nu+1}}(\tau-t_{\nu-2})(\tau-t_{\nu-1})(t_{\eta+1}-\tau)^{\xi-1}{d}\tau.\\ E^{\eta+1}& = &\frac{\lambda(1-\xi)}{A B(\xi)}t_{\eta}^{\lambda-1}E_{1}(t_{\eta},S(t_{\eta}),G(t_{\eta}))+\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}E_{1} (t_{\nu-2},S^{\nu-2},G^{\nu-2})\\&&\times t^{\lambda-1}_{\nu-2}\int_{t_{\nu}}^{t_{\nu+1}}(t_{\eta+1}-\tau)^{\xi-1}{d}\tau\ +\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}\frac{1}{\Delta t}\Big[t^{\lambda-1}_{\nu-1}E_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\ &&-t^{\lambda-1}_{\nu-2} E_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \int_{t_{\nu}}^{t_{\nu+1}}(\tau-t_{\nu-2})(t_{\eta+1}-\tau)^{\xi-1}{d}\tau\\ &&+\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}\frac{1}{2\Delta t^{2}}\Big[t^{\lambda-1}_{\nu} E_{1}(t_{\nu},S^{\nu},G^{\nu}) -2t^{\lambda-1}_{\nu-1}E_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\&&+t^{\lambda-1}_{\nu-2} E_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \int_{t_{\nu}}^{t_{\nu+1}}(\tau-t_{\nu-2})(\tau-t_{\nu-1})(t_{\eta+1}-\tau)^{\xi-1}{d}\tau.\\ V^{\eta+1}& = &\frac{\lambda(1-\xi)}{A B(\xi)}t_{\eta}^{\lambda-1}V_{1}(t_{\eta},S(t_{\eta}),G(t_{\eta}))+\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}V_{1} (t_{\nu-2},S^{\nu-2},G^{\nu-2})\\&&\times t^{\lambda-1}_{\nu-2}\int_{t_{\nu}}^{t_{\nu+1}}(t_{\eta+1}-\tau)^{\xi-1}{d}\tau +\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}\frac{1}{\Delta t} \Big[t^{\lambda-1}_{\nu-1} V_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\ &&-t^{\lambda-1}_{\nu-2} V_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \int_{t_{\nu}}^{t_{\nu+1}}(\tau-t_{\nu-2})(t_{\eta+1}-\tau)^{\xi-1}{d}\tau\\&&+\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}\frac{1}{2\Delta t^{2}}\Big[t^{\lambda-1}_{\nu}V_{1}(t_{\nu},S^{\nu},G^{\nu}) -2t^{\lambda-1}_{\nu-1}V_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\&&+t^{\lambda-1}_{\nu-2} V_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \int_{t_{\nu}}^{t_{\nu+1}}(\tau-t_{\nu-2})(\tau-t_{\nu-1})(t_{\eta+1}-\tau)^{\xi-1}{d}\tau.\\ I^{\eta+1}& = &\frac{\lambda(1-\xi)}{A B(\xi)}t_{\eta}^{\lambda-1}I_{1}(t_{\eta},S(t_{\eta}),G(t_{\eta}))+\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}I_{1} (t_{\nu-2},S^{\nu-2},G^{\nu-2})\\ &&\times t^{\lambda-1}_{\nu-2}\int_{t_{\nu}}^{t_{\nu+1}}(t_{\eta+1}-\tau)^{\xi-1}{d}\tau +\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}\frac{1}{\Delta t}\Big[t^{\lambda-1}_{\nu-1} I_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\&&-t^{\lambda-1}_{\nu-2} I_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \int_{t_{\nu}}^{t_{\nu+1}}(\tau-t_{\nu-2})(t_{\eta+1}-\tau)^{\xi-1}{d}\tau\\ &&+\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}\frac{1}{2\Delta t^{2}}\Big[t^{\lambda-1}_{\nu}I_{1}(t_{\nu},S^{\nu},G^{\nu}) -2t^{\lambda-1}_{\nu-1}Q_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\&&+t^{\lambda-1}_{\nu-2} I_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \int_{t_{\nu}}^{t_{\nu+1}}(\tau-t_{\nu-2})(\tau-t_{\nu-1})(t_{\eta+1}-\tau)^{\xi-1}{d}\tau.\\ R^{\eta+1}& = &\frac{\lambda(1-\xi)}{A B(\xi)}t_{\eta}^{\lambda-1}R_{1}(t_{\eta},S(t_{\eta}),G(t_{\eta}))+\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}R_{1} (t_{\nu-2},S^{\nu-2},G^{\nu-2})\\&&\times t^{\lambda-1}_{\nu-2}\int_{t_{\nu}}^{t_{\nu+1}}(t_{\eta+1} -\tau)^{\xi-1}{d}\tau\ +\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}\frac{1}{\Delta t}\Big[t^{\lambda-1}_{\nu-1}R_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\ &&-t^{\lambda-1}_{\nu-2}R_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \int_{t_{\nu}}^{t_{\nu+1}}(\tau-t_{\nu-2})(t_{\eta+1}-\tau)^{\xi-1}{d}\tau\\ &&+\frac{\xi \lambda}{A B(\xi)\Gamma(\xi)}\sum^{\eta}_{\nu = 2}\frac{1}{2\Delta t^{2}}\Big[t^{\lambda-1}_{\nu}R_{1}(t_{\nu},S^{\nu},G^{\nu}) -2t^{\lambda-1}_{\nu-1}R_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\&&+t^{\lambda-1}_{\nu-2} R_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \int_{t_{\nu}}^{t_{\nu+1}}(\tau-t_{\nu-2})(\tau-t_{\nu-1})(t_{\eta+1}-\tau)^{\xi-1}{d}\tau, \end{eqnarray}$

(6.2)

where $G^{\nu-2} = E^{\nu-2}, V^{\nu-2}, I^{\nu-2}, R^{\nu-2}$ , $G^{\nu-1} = E^{\nu-1}, V^{\nu-1}, I^{\nu-1}, R^{\nu-1}$ , $G^{\nu} = E^{\nu}, V^{\nu}, I^{\nu}, R^{\nu}$ , and $G(t_{\eta}) = E(t_{\eta}), V(t_{\eta}), I(t_{\eta}), R(t_{\eta})$ .

We can perform the following calculations for the integral in Eq (6.2):

$\begin{eqnarray} \int_{t_{\nu}}^{t_{\nu+1}}(t_{\eta+1}-\tau)^{\xi-1}{d}\tau& = &\frac{(\Delta t)^{\xi}}{\xi}\Big[(\eta-\nu+1)^{\xi}-(\eta-\nu)^{\xi}\Big],\\ \int_{t_{\nu}}^{t_{\nu+1}}(\tau-t_{\nu-2})(t_{\eta+1}-\tau)^{\xi-1}{d}\tau& = &\frac{(\Delta t)^{\xi+1}}{\xi(\xi+1)}\Big[(\eta-\nu+1)^{\xi}\\&&\times (\eta-\nu+3+2\xi)-(\eta-\nu)^{\xi}(\eta-\nu+3+3\xi)\Big].\\ \int_{t_{\nu}}^{t_{\nu+1}}(\tau-t_{\nu-2})(\tau-t_{\nu-1})(t_{\eta+1}-\tau)^{\xi-1}{d}\tau& = &\frac{(\Delta t)^{\xi+2}}{\xi(\xi+1)(\xi+2)}\Big[(\eta-\nu+1)^{\xi}\Big\{2(\eta-\nu)^{2}+(3\xi+10)\\&&\times(\eta-\nu)+2\xi^{2}+9\xi+12\Big\}-(\eta-\nu)^{\xi}\Big\{2(\eta-\nu)^{2}\\ &&+(5\xi+10)(\eta-\nu)+6\xi^{2}+18\xi+12\Big\}\Big], \end{eqnarray}$

(6.3)

substituting integral calculation values into Eq (6.2).

We acquire the numerical solutions $S(t), E(t), V(t), I(t)\, \, and \, \, R(t)$ :

$\begin{eqnarray} S^{\eta+1}& = &\frac{\lambda(1-\xi)}{A B(\xi)}t_{\eta}^{\lambda-1}S_{1}(t_{\eta},S(t_{\eta}),G(t_{\eta})) +\frac{\xi \lambda(\Delta t)^{\xi}}{A B(\xi)\Gamma(\xi+1)}\sum^{\eta}_{\nu = 2} S_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\\&&\times t^{\lambda-1}_{\nu-2} \Big[(\eta-\nu+1)^{\xi}-(\eta-\nu)^{\xi}\Big]+\frac{\xi \lambda(\Delta t)^{\xi}}{A B(\xi)\Gamma(\xi+2)}\sum^{\eta}_{\nu = 2} \Big[t^{\lambda-1}_{\nu-1}S_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\&&- t^{\lambda-1}_{\nu-2}S_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \Big[(\eta-\nu+1)^{\xi}(\eta-\nu+3+2\xi)-(\eta-\nu)^{\xi}(\eta-\nu+3+3\xi)\Big]\\&&+\frac{\xi \lambda(\Delta t)^{\xi}}{2A B(\xi)\Gamma(\xi+3)}\sum^{\eta}_{\nu = 2}\Big[t^{\lambda-1}_{\nu}S_{1} (t_{\nu},S^{\nu},G^{\nu}) -2t^{\lambda-1}_{\nu-1}S_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\ &&+t^{\lambda-1}_{\nu-2}S_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \Big[(\eta-\nu+1)^{\xi}\Big\{2(\eta-\nu)^{2}+(3\xi+10)(\eta-\nu)+2\xi^{2}+9\xi+12\Big\}\\ &&-(\eta-\nu)^{\xi}\times\Big\{2(\eta-\nu)^{2}+(5\xi+10)(\eta-\nu)+6\xi^{2}+18\xi+12\Big\}\Big],\\ E^{\eta+1}& = &\frac{\lambda(1-\xi)}{A B(\xi)}t_{\eta}^{\lambda-1}E_{1}(t_{\eta},S(t_{\eta}),G(t_{\eta})) +\frac{\xi \lambda(\Delta t)^{\xi}}{A B(\xi)\Gamma(\xi+1)}\sum^{\eta}_{\nu = 2} E_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\\&&\times t^{\lambda-1}_{\nu-2} \Big[(\eta-\nu+1)^{\xi}-(\eta-\nu)^{\xi}\Big]+\frac{\xi \lambda(\Delta t)^{\xi}}{A B(\xi)\Gamma(\xi+2)}\sum^{\eta}_{\nu = 2} \Big[t^{\lambda-1}_{\nu-1}E_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\ &&- t^{\lambda-1}_{\nu-2}E_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \Big[(\eta-\nu+1)^{\xi}(\eta-\nu+3+2\xi)-(\eta-\nu)^{\xi}(\eta-\nu+3+3\xi)\Big]\\ &&+\frac{\xi \lambda(\Delta t)^{\xi}}{2A B(\xi)\Gamma(\xi+3)}\sum^{\eta}_{\nu = 2}\Big[t^{\lambda-1}_{\nu}E_{1} (t_{\nu},S^{\nu},G^{\nu}) -2t^{\lambda-1}_{\nu-1}E_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\ &&+t^{\lambda-1}_{\nu-2}E_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \Big[(\eta-\nu+1)^{\xi}\Big\{2(\eta-\nu)^{2}+(3\xi+10)(\eta-\nu)+2\xi^{2}+9\xi+12\Big\}\\&&- (\eta-\nu)^{\xi}\times\Big\{2(\eta-\nu)^{2}+(5\xi+10)(\eta-\nu)+6\xi^{2}+18\xi+12\Big\}\Big],\\ V^{\eta+1}& = &\frac{\lambda(1-\xi)}{A B(\xi)}t_{\eta}^{\lambda-1}V_{1}(t_{\eta},S(t_{\eta}),G(t_{\eta}))+\frac{\xi \lambda(\Delta t)^{\xi}}{A B(\xi)\Gamma(\xi+1)}\sum^{\eta}_{\nu = 2} V_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\\&&\times t^{\lambda-1}_{\nu-2} \Big[(\eta-\nu+1)^{\xi}-(\eta-\nu)^{\xi}\Big]+\frac{\xi \lambda(\Delta t)^{\xi}}{A B(\xi)\Gamma(\xi+2)}\sum^{\eta}_{\nu = 2} \Big[t^{\lambda-1}_{\nu-1}V_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\&&- t^{\lambda-1}_{\nu-2}V_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \Big[(\eta-\nu+1)^{\xi}(\eta-\nu+3+2\xi)-(\eta-\nu)^{\xi}(\eta-\nu+3+3\xi)\Big]\\ &&+\frac{\xi \lambda(\Delta t)^{\xi}}{2A B(\xi)\Gamma(\xi+3)}\sum^{\eta}_{\nu = 2}\Big[t^{\lambda-1}_{\nu}V_{1} (t_{\nu},S^{\nu},G^{\nu}) -2t^{\lambda-1}_{\nu-1}V_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\ &&+t^{\lambda-1}_{\nu-2}V_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \Big[(\eta-\nu+1)^{\xi}\Big\{2(\eta-\nu)^{2}+(3\xi+10)(\eta-\nu)+2\xi^{2}+9\xi+12\Big\}\\&&- (\eta-\nu)^{\xi}\times\Big\{2(\eta-\nu)^{2}+(5\xi+10)(\eta-\nu)+6\xi^{2}+18\xi+12\Big\}\Big],\\ I^{\eta+1}& = &\frac{\lambda(1-\xi)}{A B(\xi)}t_{\eta}^{\lambda-1}I_{1}(t_{\eta},S(t_{\eta}),G(t_{\eta})) +\frac{\xi \lambda(\Delta t)^{\xi}}{A B(\xi)\Gamma(\xi+1)}\sum^{\eta}_{\nu = 2} I_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\\&&\times t^{\lambda-1}_{\nu-2} \Big[(\eta-\nu+1)^{\xi}-(\eta-\nu)^{\xi}\Big]+\frac{\xi \lambda(\Delta t)^{\xi}}{A B(\xi)\Gamma(\xi+2)}\sum^{\eta}_{\nu = 2} \Big[t^{\lambda-1}_{\nu-1}I_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\ &&-t^{\lambda-1}_{\nu-2}I_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \Big[(\eta-\nu+1)^{\xi}(\eta-\nu+3+2\xi)-(\eta-\nu)^{\xi}(\eta-\nu+3+3\xi)\Big]\\&&+\frac{\xi \lambda(\Delta t)^{\xi}}{2A B(\xi)\Gamma(\xi+3)}\sum^{\eta}_{\nu = 2}\Big[t^{\lambda-1}_{\nu}I_{1} (t_{\nu},S^{\nu},G^{\nu}) -2t^{\lambda-1}_{\nu-1}I_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\ &&+t^{\lambda-1}_{\nu-2}I_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \Big[(\eta-\nu+1)^{\xi}\Big\{2(\eta-\nu)^{2}+(3\xi+10)(\eta-\nu)+2\xi^{2}+9\xi+12\Big\}\\ &&-(\eta-\nu)^{\xi}\times\Big\{2(\eta-\nu)^{2}+(5\xi+10)(\eta-\nu)+6\xi^{2}+18\xi+12\Big\}\Big],\\ R^{\eta+1}& = &\frac{\lambda(1-\xi)}{A B(\xi)}t_{\eta}^{\lambda-1}R_{1}(t_{\eta},S(t_{\eta}),G(t_{\eta}))+\frac{\xi \lambda(\Delta t)^{\xi}}{A B(\xi)\Gamma(\xi+1)}\sum^{\eta}_{\nu = 2} R_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\\&&\times t^{\lambda-1}_{\nu-2} \Big[(\eta-\nu+1)^{\xi}-(\eta-\nu)^{\xi}\Big]+\frac{\xi \lambda(\Delta t)^{\xi}}{A B(\xi)\Gamma(\xi+2)}\sum^{\eta}_{\nu = 2} \Big[t^{\lambda-1}_{\nu-1}R_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\&&- t^{\lambda-1}_{\nu-2}R_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \Big[(\eta-\nu+1)^{\xi}(\eta-\nu+3+2\xi)-(\eta-\nu)^{\xi}(\eta-\nu+3+3\xi)\Big]\\ &&+\frac{\xi \lambda(\Delta t)^{\xi}}{2 A B(\xi)\Gamma(\xi+3)}\sum^{\eta}_{\nu = 2}\Big[t^{\lambda-1}_{\nu}R_{1} (t_{\nu},S^{\nu},G^{\nu}) -2t^{\lambda-1}_{\nu-1}R_{1}(t_{\nu-1},S^{\nu-1},G^{\nu-1})\\ &&+t^{\lambda-1}_{\nu-2}R_{1}(t_{\nu-2},S^{\nu-2},G^{\nu-2})\Big] \Big[(\eta-\nu+1)^{\xi}\Big\{2(\eta-\nu)^{2}+(3\xi+10)(\eta-\nu)+2\xi^{2}+9\xi+12\Big\}\\&&-(\eta-\nu)^{\xi}\times\Big\{2(\eta-\nu)^{2}+(5\xi+10)(\eta-\nu)+6\xi^{2}+18\xi+12\Big\}\Big]. \end{eqnarray}$

7. Simulation explanation

In this section, we utilized an advanced technique to obtain theoretical outcomes and assess their effectiveness. The newly developed SEVIR system was analyzed through simulation. By applying non-integer parametric values in the SARS-COVID-19 model, we obtained interesting findings. – display the solutions for $S(t)$ , $E(t)$ , $V(t)$ , $I(t)$ , and $R(t)$ by reducing the fractional values to the desired level. To validate the efficiency of the theoretical outcomes, we provide the following examples. Numerical simulations for the SARS-COVID-19 model were performed using MATLAB. The initial conditions used in the newly developed model are $S(0) = 217.342565, E(0) = 100, V(0) = 1.386348, I(0) = 1.1$ and $R(0) = 1.271087$ . The parameter values used in the developed system are as follows: $\textbf{a} = 1.43$ , $\delta = 0.000065$ , $\rho = 0.45$ , $\rho_{1} = 0.10$ , $\rho_{2} = 0.020$ , $\alpha = 0.0008601$ , $\mu_{0} = 0.19$ , $\omega = 0.98$ , $\textbf{b} = 0.135$ , $\phi = 0.0001$ , and $\psi = 0.0001$ . Figures 2, 4, and 5 illustrate the changes in susceptible, vaccinated, and infected individuals respectively, showing a sharp decrease before reaching a stable position. Meanwhile, Figures 3 and 6 demonstrate the dynamics of exposed and recovered individuals respectively at different fractional orders in which both individuals increases and after certain time the number of individuals approaches stable state using different dimensions. The research predicts future infection rates and suggests ways to decrease the spread of infection units more effectively. By utilizing a fractal-fractional approach, the study yields reliable and accurate results for all compartments at non-integer order derivatives, which are more trustworthy when fractional values are reduced as well a by reducing its dimensions. The findings suggest that the number of infected individuals decreases significantly due to vaccination measures, while the number of recovered individuals increases due to a decline in infected individuals and the effect of vaccination.

Figure 2. The value of

$S(t)$ using the fractal-fractional operator with various fractional values at different dimensions.

DownLoad: Full-Size Img PowerPoint

Figure 3. The value of

$E(t)$ using the fractal-fractional operator with various fractional values at different dimensions.

DownLoad: Full-Size Img PowerPoint

Figure 4. The value of

$V(t)$ using the fractal-fractional operator with various fractional values at different dimensions.

DownLoad: Full-Size Img PowerPoint

Figure 5. The value of

$I(t)$ using the fractal-fractional operator with various fractional values at different dimensions.

DownLoad: Full-Size Img PowerPoint

Figure 6. The value of

$R(t)$ using the fractal-fractional operator with various fractional values at different dimensions.

DownLoad: Full-Size Img PowerPoint

8. Conclusions

This article employs a fractional order SEVIR model for SARS-COVID-19 with vaccinated effects using an FFO to find reliable solutions. We provide advice on controlling this virus to help our community overcome the pandemic by implementing vaccinated measures for low immune individuals. We analyze the dangerous coronavirus disease with the effect of vaccination to understand its real impact on the community. Qualitative and quantitative analyses are conducted to verify its stable position in a continuous dynamical system. We also verify that the fractional order coronavirus disease model has bounded and unique solutions. We examine the impact of global measures to control the spread of the coronavirus disease. Also, analyses are performed to see how the rate of infection changes after the implementation of vaccination measures. We ensure that our findings are reliable and realistic. FFO is used for continuously monitoring the spread as well as control of the disease in society after vaccination measures. It was observed that infected individuals recover quickly due to the vaccinated strategy. The fractal-fractional operator (FFO) is used for continuously monitoring for the spread of the diseases using different fractional values as well as reliable solutions. In fractal-fractional operators, fractal represents the dimensions of the spread of the disease, and fractional represents the fractional ordered derivative operator which provides the real behavior of the spread as well as control of COVID-19 with different dimensions and continuous monitoring, which can be observed in simulation. We conduct numerical simulations to observe how the disease controlled in the community after the implementation of vaccination measures using different fractional values with different dimensions. Additionally, future estimates are provided based on our findings, which can help in mitigating the risk of the disease spreading in the environment.

Use of AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

The authors would like to extend their sincere appreciation to Researchers Supporting Project number (RSP2024R472), King Saud University, Riyadh, Saudi Arabia.

Conflict of interest

All authors declare no conflicts of interest in this paper.

References

[1]	Curry C, Bloomberg New Energy Finance, 2017. Lithium-ion battery costs and market.
[2]	Bloomberg New Energy Finance, 2020. Electric Vehicle Outlook. Available from: https://about.bnef.com/electric-vehicle-outlook/.
[3]	HS Markit, 2020. The battery electric vehicle: Why mass adoption is inevitable, yet elusive.
[4]	Samsung SDI. The Composition of EV Batteries: Cells? Modules? Packs? Let's Understand Properly! Available from: https://www.samsungsdi.com/column/all/detail/54344.html?listType=&pageIndex=1&searchCondition=&searchKeyword=&idx=54344.
[5]	Babbitt CW (2020) Sustainability perspectives on lithium-ion batteries. Clean Technol Envir 22: 1213–1214. https://doi.org/10.1007/s10098-020-01890-3 doi: 10.1007/s10098-020-01890-3
[6]	Richa K, Babbitt CW, Gaustad G (2017) Eco-efficiency analysis of a lithium-ion battery waste hierarchy inspired by circular economy. J Ind Ecol 21: 715–730. https://doi.org/10.1111/jiec.12607 doi: 10.1111/jiec.12607
[7]	Liu WW, Zhang H, Liu LH, et al. (2016) Remanufacturing cathode from end-of-life of lithium-ion secondary batteries by Nd:YAG laser radiation. Clean Technol Envir 18: 231–243. https://doi.org/10.1007/s10098-015-1010-1 doi: 10.1007/s10098-015-1010-1
[8]	Zwicker MFR, Moghadam M, Zhang W, et al. (2020) Automotive battery pack manufacturing-a review of battery to tab joining. J Adv Joi Process 1: 100017. https://doi.org/10.1016/j.jajp.2020.100017 doi: 10.1016/j.jajp.2020.100017
[9]	Das A, Li D, Williams D, et al. (2018) Joining technologies for automotive battery systems manufacturing. World Electr Veh J 9: 22. https://doi.org/10.3390/wevj9020022 doi: 10.3390/wevj9020022
[10]	Shui L, Chen F, Garg A, et al. (2018) Design optimization of battery pack enclosure for electric vehicle. Struct Multidisc Optim 58: 331–347. https://doi.org/10.1007/s00158-018-1901-y doi: 10.1007/s00158-018-1901-y
[11]	Lee JH, Lee HM, Ahn S (2003) Battery dimensional changes occurring during charge/discharge cycles-thin rectangular lithium ion and polymer cells. J Power Sources 119–121: 833–837. https://doi.org/10.1016/S0378-7753(03)00281-7 doi: 10.1016/S0378-7753(03)00281-7
[12]	Zhao N, Zhao D, Xu L, et al. (2014) A multimode responsive aptasensor for adenosine detection. J Nanomater 2014: 1–7. https://doi.org/10.1155/2014/360347 doi: 10.1155/2014/360347
[13]	Fleckenstein M, Bohlen O, Roscher MA, et al. (2011) Current density and state of charge inhomogeneities in Li-ion battery cells with LiFePO₄ as cathode material due to temperature gradients. J Power Sources 196: 4769–4778. https://doi.org/10.1016/j.jpowsour.2011.01.043 doi: 10.1016/j.jpowsour.2011.01.043
[14]	Yang N, Zhang X, Shang B, et al. (2016) Unbalanced discharging and aging due to temperature differences among the cells in a lithium-ion battery pack with parallel combination. J Power Sources 306: 733–741. https://doi.org/10.1016/j.jpowsour.2015.12.079 doi: 10.1016/j.jpowsour.2015.12.079
[15]	Wang Q, Ping P, Zhao X, et al. (2012) Thermal runaway caused fire and explosion of lithium ion battery. J Power Sources 208: 210–224. https://doi.org/10.1016/j.jpowsour.2012.02.038 doi: 10.1016/j.jpowsour.2012.02.038
[16]	Baumann M, Wildfeuer L, Rohr S, et al. (2018) Parameter variations within Li-ion battery packs-Theoretical investigations and experimental quantification. J Energy Storage 18: 295–307. https://doi.org/10.1016/j.est.2018.04.031 doi: 10.1016/j.est.2018.04.031
[17]	Baumhöfer T, Brühl M, Rothgang S, et al. (2014) Production caused variation in capacity aging trend and correlation to initial cell performance. J Power Sources 247: 332–338. https://doi.org/10.1016/j.jpowsour.2013.08.108 doi: 10.1016/j.jpowsour.2013.08.108
[18]	Brand MJ, Schmidt PA, Zaeh MF, et al. (2015) Welding techniques for battery cells and resulting electrical contact resistances. J Energy Storage 1: 7–14. https://doi.org/10.1016/j.est.2015.04.001 doi: 10.1016/j.est.2015.04.001
[19]	Hesse T, Zaske S (2018) A "World Champion" in Energy Efficiency: Environmental audit sees laser welding ahead of resistance spot welding. Laser Tech J 15: 56–58. https://doi.org/10.1002/latj.201800014 doi: 10.1002/latj.201800014
[20]	Kohn S, Berdichevsky G, Hewett BC (2011) Tunable frangible battery pack system. US Patent 7923144B2.
[21]	Shin HS, de Leon M (2017) Mechanical performance and electrical resistance of ultrasonic welded multiple Cu–Al layers. J Mater Process Tech 241: 141–153. https://doi.org/10.1016/j.jmatprotec.2016.11.004 doi: 10.1016/j.jmatprotec.2016.11.004
[22]	Lee SS, Kim TH, Hu SJ, et al. (2010) Joining technologies for automotive lithium-ion battery manufacturing: A review. ASME 2010 International Manufacturing Science and Engineering Conference (MSEC 2010), Erie, Pennsylvania, USA, 1: 541–549. https://doi.org/10.1115/MSEC2010-34168
[23]	Wagner G, Balle F, Eifler D (2012) Ultrasonic welding of hybrid joints. JOM 64: 401–406. https://doi.org/10.1007/s11837-012-0269-5 doi: 10.1007/s11837-012-0269-5
[24]	Chamberlain A, Larsson H, Nilsson L, et al. (2019) Welding methods for electrical connections in battery systems. Uppsala University.
[25]	Das A, Ashwin TR, Barai A (2019) Modelling and characterisation of ultrasonic joints for Li-ion batteries to evaluate the impact on electrical resistance and temperature raise. J Energy Storage 22: 239–248. https://doi.org/10.1016/j.est.2019.02.017 doi: 10.1016/j.est.2019.02.017
[26]	McGovern ME, Rinker TJ, Sekol RC (2019) Assessment of ultrasonic welds using pulsed infrared thermography. JNDE 2: 011009. https://doi.org/10.1115/1.4042260 doi: 10.1115/1.4042260
[27]	Mohan Raj N, Kumaraswamidhas LA, Nalajam PK, et al. (2018) Studies on electro mechanical aspects in ultrasonically welded Al/Cu joints. Trans Indian Inst Met 71: 107–116. https://doi.org/10.1007/s12666-017-1140-8 doi: 10.1007/s12666-017-1140-8
[28]	Choi S, Fuhlbrigge T, Nidamarthi S (2012) Vibration analysis in robotic ultrasonic welding for battery assembly. 2012 IEEE International Conference on Automation Science and Engineering (CASE), Seoul, Korea, 550–554. https://doi.org/10.1109/CoASE.2012.6386349
[29]	Kang B, Cai W, Tan CA (2014) Dynamic stress analysis of battery tabs under ultrasonic welding. J Manuf Sci Eng 136: 041011. https://doi.org/10.1115/1.4026990 doi: 10.1115/1.4026990
[30]	Zhao J, Li H, Choi H, et al. (2013) Insertable thin film thermocouples for in situ transient temperature monitoring in ultrasonic metal welding of battery tabs. J Manuf Process 15: 136–140. https://doi.org/10.1016/j.jmapro.2012.10.002 doi: 10.1016/j.jmapro.2012.10.002
[31]	Li H, Choi H, Ma C, et al. (2013) Transient temperature and heat flux measurement in ultrasonic joining of battery tabs using thin-film microsensors. J Manuf Sci Eng 135: 051015. https://doi.org/10.1115/1.4024816 doi: 10.1115/1.4024816
[32]	Zhang H, Senkara J (2011) Resistance Welding: Fundamentals and Applications, 2 Eds., CRC Press. https://doi.org/10.1201/b11752
[33]	Saariluoma H, Piiroinen A, Unt A, et al. (2020) Overview of optical digital measuring challenges and technologies in laser welded components in EV battery module design and manufacturing. Batteries 6: 47. https://doi.org/10.3390/batteries6030047 doi: 10.3390/batteries6030047
[34]	Mehta KP, Badheka VJ (2016) A review on dissimilar friction stir welding of copper to aluminum: Process, properties, and variants. Mater Manuf Process 31: 233–254. https://doi.org/10.1080/10426914.2015.1025971 doi: 10.1080/10426914.2015.1025971
[35]	Mypati O, Mishra D, Sahu S, et al. (2020) A study on electrical and electrochemical characteristics of friction stir welded lithium-ion battery tabs for electric vehicles. J Electron Mater 49: 72–87. https://doi.org/10.1007/s11664-019-07711-8 doi: 10.1007/s11664-019-07711-8
[36]	Li W, Wen Q, Yang X, et al. (2017) Interface microstructure evolution and mechanical properties of Al/Cu bimetallic tubes fabricated by a novel friction-based welding technology. Mater Design 134: 383-393. https://doi.org/10.1016/j.matdes.2017.08.065 doi: 10.1016/j.matdes.2017.08.065
[37]	Brand MJ, Berg P, Kolp EI, et al. (2016) Detachable electrical connection of battery cells by press contacts. J Energy Storage 8: 69–77. https://doi.org/10.1016/j.est.2016.09.011 doi: 10.1016/j.est.2016.09.011
[38]	Bolsinger C, Zorn M, Birke KP (2017) Electrical contact resistance measurements of clamped battery cell connectors for cylindrical 18650 battery cells. J Energy Storage 12: 29–36. https://doi.org/10.1016/j.est.2017.04.001 doi: 10.1016/j.est.2017.04.001
[39]	Taheri P, Hsieh S, Bahrami M (2011) Investigating electrical contact resistance losses in lithium-ion battery assemblies for hybrid and electric vehicles. J Power Sources 196: 6525–6533. https://doi.org/10.1016/j.jpowsour.2011.03.056 doi: 10.1016/j.jpowsour.2011.03.056
[40]	Fu R, Choe SW, Jackson RL, et al. (2012) Modeling and analysis of vibration-induced changes in connector resistance of high power electrical connectors for hybrid vehicles. Mech Based Des Struc 40: 349–365. https://doi.org/10.1080/15397734.2012.670098 doi: 10.1080/15397734.2012.670098
[41]	American Welding Society (1982) Welding Handbook: Metals and Their Weldability, 7 Eds. https://www.amazon.com/Welding-Handbook-Metals-Their-Weldability/dp/0871712180.
[42]	Davis JR, Ferutz K, Wheaton N (1980) Asm Handbook: Welding, Brazing, and Soldering, ASM International, 6.
[43]	American Welding Society (1991) Brazing Handbook, 4 Eds..
[44]	Hayashi J, Miyazawa Y (2014) Spot brazing of aluminum to copper with a cover plate. IOP Conf Ser-Mater Sci Eng 61: 012015. https://doi.org/10.1088/1757-899X/61/1/012015 doi: 10.1088/1757-899X/61/1/012015
[45]	Solchenbach T, Plapper P, Cai W (2014) Electrical performance of laser braze-welded aluminum-copper interconnects. J Manuf Process 16: 183–189. https://doi.org/10.1016/j.jmapro.2013.12.002 doi: 10.1016/j.jmapro.2013.12.002
[46]	Brand MJ, Kolp EI, Berg P, et al. (2017) Electrical resistances of soldered battery cell connections. J Energy Storage 12: 45–54. https://doi.org/10.1016/j.est.2017.03.019 doi: 10.1016/j.est.2017.03.019
[47]	Sommer AK, Türpe M, Füssel U, et al. (2020) Considerations on a new brazing concept for vacuum brazing of aluminium and copper. Weld World 64: 1589–1596. https://doi.org/10.1007/s40194-020-00935-x doi: 10.1007/s40194-020-00935-x
[48]	Krishnaja D, Cheepu M, Venkateswarlu D (2018) A review of research progress on dissimilar laser weld-brazing of automotive applications. IOP Conf Ser-Mater Sci Eng 330: 012073. https://doi.org/10.1088/1757-899X/330/1/012073 doi: 10.1088/1757-899X/330/1/012073
[49]	Markovits T, Jászberényi A (2018) Investigation of laser-material interaction in case of aluminium brazing process. Period Polytech Transp Eng 46: 59–62. https://doi.org/10.3311/PPtr.11528 doi: 10.3311/PPtr.11528
[50]	Li Y, Chen C, Yi R, et al. (2020) Review: Special brazing and soldering. J Manuf Process 60: 608–635. https://doi.org/10.1016/j.jmapro.2020.10.049 doi: 10.1016/j.jmapro.2020.10.049
[51]	Xiao R, Dong P (2011) Laser penetration brazing of dissimilar alloys. International Congress on Applications of Lasers & Electro-Optics, Orlando, Florida, USA, 2011: 308–314. https://doi.org/10.2351/1.5062251
[52]	Deyev G, Deyev D (2005) Surface Phenomena in Fusion Welding Processes, 1 Ed., CRC Press. https://doi.org/10.1201/9781420036299-1
[53]	Zhou L, Luo LY, Tan CW, et al. (2018) Effect of welding speed on microstructural evolution and mechanical properties of laser welded-brazed Al/brass dissimilar joints. Opt Laser Technol 98: 234–246. https://doi.org/10.1016/j.optlastec.2017.08.004 doi: 10.1016/j.optlastec.2017.08.004
[54]	Xia C, Li Y, Puchkov UA, et al. (2008) Microstructure and phase constitution near the interface of Cu/Al vacuum brazing using Al–Si filler metal. Vacuum 82: 799–804. https://doi.org/10.1016/j.vacuum.2007.11.007 doi: 10.1016/j.vacuum.2007.11.007
[55]	Furuya HS, Sato YS, Kokawa H, et al. (2018) Improvement of interfacial strength with the addition of Ni in Al/Cu dissimilar joints produced via laser brazing. Metall Mater Trans A 49: 6215–6223. https://doi.org/10.1007/s11661-018-4938-z doi: 10.1007/s11661-018-4938-z
[56]	Sun Z, Ion JC (1995) Laser welding of dissimilar metal combinations. J Mater Sci 30: 4205–4214. https://doi.org/10.1007/BF00361499 doi: 10.1007/BF00361499
[57]	Steen WM (1987) Lasers in Manufacturing: Proceedings of the 4th International Conference, Birmingham, UK, Kempston, Bedford: IFS Publications; Berlin, New York: Springer-Verlag, 12–14.
[58]	Weigl M, Albert F, Schmidt M (2011) Enhancing the ductility of laser-welded copper-aluminum connections by using adapted filler materials. Physics Procedia 12: 332–338. https://doi.org/10.1016/j.phpro.2011.03.141 doi: 10.1016/j.phpro.2011.03.141
[59]	Sun Z, Moisio T (1993) Laser beam welding of austenitic/ferritic dissimilar steel joints using nickel based filler wire. Mater Sci Tech 9: 603–608. https://doi.org/10.1179/mst.1993.9.7.603 doi: 10.1179/mst.1993.9.7.603
[60]	60. Kancharla V, Mendes M, Grupp M, et al. (2018) Recent advances in fiber laser welding. Biul Inst Spawal Gliwicach 62: 175–181. https://doi.org/10.17729/ebis.2018.5/20
[61]	61. Blecher JJ, Galbraith CM, Van Vlack C, et al. (2014) Real time monitoring of laser beam welding keyhole depth by laser interferometry. Sci Technol Weld Joi 19: 560–564. https://doi.org/10.1179/1362171814Y.0000000225 doi: 10.1179/1362171814Y.0000000225
[62]	Vasudevan AK, Doherty RD (1989) Aluminum Alloys-Contemporary Research and Applications, Boston: Academic Press.
[63]	Katgerman L, Soetens F (2010) New frontiers in light metals. Proceedings of the 11th International Aluminium Conference INALCO 2010, the Netherlands, 23–25.
[64]	Kah P, Rajan R, Martikainen J, et al. (2015) Investigation of weld defects in friction-stir welding and fusion welding of aluminium alloys. Int J Mech Mater Eng 10: 26. https://doi.org/10.1186/s40712-015-0053-8 doi: 10.1186/s40712-015-0053-8
[65]	Sun T, Franciosa P, Ceglarek D (2021) Effect of focal position offset on joint integrity of AA1050 battery busbar assembly during remote laser welding. J Mater Res Technol 14: 2715–2726. https://doi.org/10.1016/j.jmrt.2021.08.002 doi: 10.1016/j.jmrt.2021.08.002
[66]	Kirchhoff M (2021) Welding: E-mobility, Trumpf Laser- und Systentechnik GmbH.
[67]	Laserline GmbH, 2021. Laserline LDblue welding of electrical copper connectors.
[68]	Kraetzsch M, Standfuss J, Klotzbach A, et al. (2011) Laser beam welding with high-frequency beam oscillation: Welding of dissimilar materials with brilliant fiber lasers. Physics Procedia 12: 142–149. https://doi.org/10.1016/j.phpro.2011.03.018 doi: 10.1016/j.phpro.2011.03.018
[69]	Dimatteo V, Ascari A, Fortunato A (2019) Continuous laser welding with spatial beam oscillation of dissimilar thin sheet materials (Al–Cu and Cu–Al): Process optimization and characterization. J Manuf Process 44: 158–165. https://doi.org/10.1016/j.jmapro.2019.06.002 doi: 10.1016/j.jmapro.2019.06.002
[70]	Lerra F, Ascari A, Fortunato A (2019) The influence of laser pulse shape and separation distance on dissimilar welding of Al and Cu films. J Manuf Process 45: 331–339. https://doi.org/10.1016/j.jmapro.2019.07.015 doi: 10.1016/j.jmapro.2019.07.015
[71]	Olowinsky A (2021) How Lasers improve battery production. Available from: https://www.ilt.fraunhofer.de/de/mediathek/vortraege.html.
[72]	Huang W, Wang H, Rinker T, et al. (2020) Investigation of metal mixing in laser keyhole welding of dissimilar metals. Mater Design 195: 109056. https://doi.org/10.1016/j.matdes.2020.109056 doi: 10.1016/j.matdes.2020.109056
[73]	Schmalen P, Plapper P, Cai W (2016) Process robustness of laser braze-welded Al/Cu connectors. SAE Int J Alt Power 5: 195–204. https://doi.org/10.4271/2016-01-1198 doi: 10.4271/2016-01-1198
[74]	Schmalen P, Plapper P (2016) Evaluation of laser braze-welded dissimilar Al–Cu joints. Physics Procedia 83: 506–514. https://doi.org/10.1016/j.phpro.2016.08.052 doi: 10.1016/j.phpro.2016.08.052
[75]	Reisgen U, Olschok S, Jakobs S, et al. (2018) Influence of the degree of dilution with laser beam vacuum-welded Cu–Al mixed joints on the electrical properties. Procedia CIRP 74: 23–26. https://doi.org/10.1016/j.procir.2018.08.022 doi: 10.1016/j.procir.2018.08.022
[76]	Solchenbach T, Plapper P (2013) Mechanical characteristics of laser braze-welded aluminium-copper connections. Optics Laser Technol 54: 249–256. https://doi.org/10.1016/j.optlastec.2013.06.003 doi: 10.1016/j.optlastec.2013.06.003
[77]	Britten SW, Schmid L, Molitor T, et al. (2020) Blue high-power laser sources for processing solutions in e-mobility and beyond. Procedia CIRP 94: 592–595. https://doi.org/10.1016/j.procir.2020.09.082 doi: 10.1016/j.procir.2020.09.082
[78]	Zediker MS, Fritz RD, Finuf MJ, et al. (2020) Laser welding components for electric vehicles with a high-power blue laser system. J Laser Appl 32: 022038. https://doi.org/10.2351/7.0000054 doi: 10.2351/7.0000054
[79]	Mathivanan K, Plapper P (2019) Laser welding of dissimilar copper and aluminum sheets by shaping the laser pulses. Procedia Manuf 36: 154–162. https://doi.org/10.1016/j.promfg.2019.08.021 doi: 10.1016/j.promfg.2019.08.021
[80]	Schmalen P, Plapper P (2017) Resistance measurement of laser welded dissimilar Al/Cu joints. J Laser Micro Nanoen 12: 189–194. https://doi.org/10.2961/jlmn.2017.03.0003 doi: 10.2961/jlmn.2017.03.0003
[81]	Schmalen P, Plapper P, Peral I, et al. (2018) Composition and phases in laser welded Al–Cu joints by synchrotron X-ray microdiffraction. Procedia CIRP 74: 27–32. https://doi.org/10.1016/j.procir.2018.08.006 doi: 10.1016/j.procir.2018.08.006
[82]	Solchenbach T, Plapper P, Greger M, et al. (2014) Thermal and electrical aging of laser braze-welded aluminum-copper interconnects. Transl Mater Res 1: 015001. https://doi.org/10.1088/2053-1613/1/1/015001 doi: 10.1088/2053-1613/1/1/015001
[83]	Mathivanan K, Plapper PW (2019) Laser overlap joining from copper to aluminum and analysis of failure zone. Lasers in Manufacturing Conference 2019, Munich, Germany.
[84]	Hailat MM, Mian A, Chaudhury ZA, et al. (2012) Laser micro-welding of aluminum and copper with and without tin foil alloy. Microsyst Technol 18: 103–112. https://doi.org/10.1007/s00542-011-1378-8 doi: 10.1007/s00542-011-1378-8
[85]	Helm J, Dietz von Bayer I, Olowinsky A, et al. (2019) Influence of the surface properties of the connector material on the reliable and reproducible contacting of battery cells with a laser beam welding process. Weld World 63: 1221–1228. https://doi.org/10.1007/s40194-019-00727-y doi: 10.1007/s40194-019-00727-y
[86]	Lee SJ, Nakamura H, Kawahito Y, et al. (2014) Effect of welding speed on microstructural and mechanical properties of laser lap weld joints in dissimilar Al and Cu sheets. Sci Technol Weld Joi 19: 111–118. https://doi.org/10.1179/1362171813Y.0000000168 doi: 10.1179/1362171813Y.0000000168
[87]	Risdiyanto A, Arifin M, Khayam U (2013) Study on temperature distribution at busbar connection based on contact resistance of different plating contact surface. 2013 Joint International Conference on Rural Information & Communication Technology and Electric-Vehicle Technology (rICT & ICeV-T), Bandung, Indonesia. https://doi.org/10.1109/rICT-ICeVT.2013.6741537
[88]	Farahat MA, Gockenbach E, El-Alaily AA, et al. (1996) Effect of coating materials on the electrical performance of copper joints. Proceedings of the Forty-Second IEEE Holm Conference on Electrical Contacts, 472–478. https://doi.org/10.1109/HOLM.1996.557229
[89]	Wong TT, Liang GY, An G, et al. (2005) The electrical conductivity of laser-remelted and plasma-sprayed Ni and Cr coatings. J Mater Process Tech 159: 265–271. https://doi.org/10.1016/j.jmatprotec.2004.05.012 doi: 10.1016/j.jmatprotec.2004.05.012
[90]	Swedish Institute of Standards (2016) Aluminum and aluminum alloy-Sheet metal and strip-Part 2: Mechanical properties. SS-EN 485-2.
[91]	ASM Handbook Committee (1990) Properties and selection: Nonferrous alloys and special-purpose materials, ASM Handbook, ASM International, 2.
[92]	Swedish Institute of Standards (2019) Aluminium and aluminium alloys—Chemical composition and form of wrought products-Part 3: Chemical composition and form of products. SS-EN 573-3.
[93]	Solchenbach T, Plapper P (2013) Combined laser beam braze-welding process for fluxless Al–Cu connections. International Conference on competitive Manufacturing (COMA), 13: 131–136.
[94]	Kah P, Vimalraj C, Martikainen J, et al. (2015) Factors influencing Al–Cu weld properties by intermetallic compound formation. Int J Mech Mater Eng 10: 10. https://doi.org/10.1186/s40712-015-0037-8 doi: 10.1186/s40712-015-0037-8
[95]	Lee YG, Duh JG (1998) Characterizing the formation and growth of intermetallic compound in the solder joint. J Mater Sci 33: 5569–5572. https://doi.org/10.1023/A:1004499728840 doi: 10.1023/A:1004499728840
[96]	Mys I, Schmidt M (2006) Laser micro welding of copper and aluminum. Lasers and Applications in Science and Engineering, San Jose, California, United States, 610703. https://doi.org/10.1117/12.648376
[97]	Watanabe W, Onda S, Tamaki T, et al. (2006) Space-selective laser joining of dissimilar transparent materials using femtosecond laser pulses. Appl Phys Lett 89: 021106. https://doi.org/10.1063/1.2221393 doi: 10.1063/1.2221393
[98]	Braunovic M (2007) Reliability of power connections. J Zhejiang Univ-Sc A 8: 343–356. https://doi.org/10.1631/jzus.2007.A0343 doi: 10.1631/jzus.2007.A0343
[99]	Kaplan AFH (2005) Model of the absorption variation during pulsed laser heating applied to welding of electronic Au/Ni-coated Cu-leadframes. Appl Surf Sci 241: 362–370. https://doi.org/10.1016/j.apsusc.2004.07.035 doi: 10.1016/j.apsusc.2004.07.035
[100]	Huntington CA, Eagar TW (1983) Laser welding of aluminum and aluminum alloys. Weld Res Suppl 62: 105–107.
[101]	Kannatey-Asibu Jr E (2009) Principles of Laser Materials Processing, John Wiley & Sons. https://doi.org/10.1002/9780470459300
[102]	Delgado-Pinar M, Zalvidea D, Diez A, et al. (2006) Q-switching of an all-fiber laser by acousto-optic modulation of a fiber Bragg grating. Opt Express 14: 1106. https://doi.org/10.1364/OE.14.001106 doi: 10.1364/OE.14.001106
[103]	Heider A, Weber R, Herrmann D, et al. (2015) Power modulation to stabilize laser welding of copper. J Laser Appl 27: 022003. https://doi.org/10.2351/1.4906127 doi: 10.2351/1.4906127
[104]	Minasian RA (2005) MODULATORS\|modulation and demodulation of optical signals, In: Guenther CD, Encyclopedia of Modern Optics, Elsevier, 129–138. https://doi.org/10.1016/B0-12-369395-0/00667-9
[105]	Purdue University, Modulation of light. Available from: https://www.physics.purdue.edu/webapps/index.php/course_document/index/phys570P/1684/25/14317
[106]	Markettech, 2008. Overview of modulated and pulsed diode laser systems. Available from: https://www.markettechinc.net/wp-content/uploads/2018/05/Modulated-and-Pulsed-Diode-Lasers.pdf.
[107]	Schiry M, Plapper P (2019) Novel opportunities by laser welding of dissimilar materials. Proceedings of International Conference on Competitive Manufacturing COMA'19, South Africa, 293–300.
[108]	Müller A, Goecke SF, Rethmeier M (2018) Laser beam oscillation welding for automotive applications. Weld World 62: 1039–1047. https://doi.org/10.1007/s40194-018-0625-3 doi: 10.1007/s40194-018-0625-3
[109]	Xue Z, Hu S, Zuo D, et al. (2013) Molten pool characterization of laser lap welded copper and aluminum. J Phys D Appl Phys 46: 495501. https://doi.org/10.1088/0022-3727/46/49/495501 doi: 10.1088/0022-3727/46/49/495501
[110]	Deyneka-Dupriez N, Denkel A (2020) For a better world: OCT unlocks real-time quality control during oscillation welding. Lessmüller Lasertechnik. Available from: https://www.lasersystemseurope.com/analysis-opinion/better-world-oct-unlocks-real-time-quality-control-during-oscillation-welding.
[111]	Müller A, Goecke SF, Sievi P, et al. (2014) Laser beam oscillation strategies for fillet welds in lap joints. Physics Procedia 56: 458–466. https://doi.org/10.1016/j.phpro.2014.08.149 doi: 10.1016/j.phpro.2014.08.149
[112]	Dimatteo V, Ascari A, Fortunato A (2021) Dissimilar laser welding of copper and aluminum alloys in multilayer configuration for battery applications. J Laser Appl 33: 042028. https://doi.org/10.2351/7.0000476 doi: 10.2351/7.0000476
[113]	Wahba M, Kawahito Y, Katayama S (2011) Laser direct joining of AZ91D thixomolded Mg alloy and amorphous polyethylene terephthalate. J Mater Process Tech 211: 1166–1174. https://doi.org/10.1016/j.jmatprotec.2011.01.021 doi: 10.1016/j.jmatprotec.2011.01.021
[114]	Schmidt PA, Pauleser T, Zaeh MF (2014) Optimisation of weld seam configurations using a genetic algorithm. Procedia CIRP 25: 393–399. https://doi.org/10.1016/j.procir.2014.10.054 doi: 10.1016/j.procir.2014.10.054
[115]	Mai TA, Spowage AC (2004) Characterisation of dissimilar joints in laser welding of steel-kovar, copper-steel and copper-aluminium. Mater Sci Eng-A Struct 374: 224–233. https://doi.org/10.1016/j.msea.2004.02.025 doi: 10.1016/j.msea.2004.02.025
[116]	Cai ZP, Ai BQ, Cao R, et al. (2016) Microstructure and properties of aluminum AA6061-T6 to copper (Cu)-T2 joints by cold metal transfer joining technology. J Mater Res 31: 2876–2887. https://doi.org/10.1557/jmr.2016.295 doi: 10.1557/jmr.2016.295
[117]	Mehta KP, Badheka VJ (2017) Hybrid approaches of assisted heating and cooling for friction stir welding of copper to aluminum joints. J Mater Process Tech 239: 336–345. https://doi.org/10.1016/j.jmatprotec.2016.08.037 doi: 10.1016/j.jmatprotec.2016.08.037
[118]	Esser G, Mys I, Schmidt MH (2004) Laser micro welding of copper and aluminium using filler materials. Fifth International Symposium on Laser Precision Microfabrication 5662: 337–342. https://doi.org/10.1117/12.596330 doi: 10.1117/12.596330
[119]	Tavassoli S, Abbasi M, Tahavvori R (2016) Controlling of IMCs layers formation sequence, bond strength and electrical resistance in Al–Cu bimetal compound casting process. Mater Design 108: 343–353. https://doi.org/10.1016/j.matdes.2016.06.076 doi: 10.1016/j.matdes.2016.06.076
[120]	Shi Y, Zhou XL, Zhu M, et al. (2017) Effect of filler wires on brazing interface microstructure and mechanical properties of Al/Cu dissimilar metals welding-brazing joint. Mater Rep 31: 61–64 (in Chinese).

Reader Comments

Your name:*

Email:*
© 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)