Published: 25 October 2022

The dual-phase-lag bioheat transfer of a skin tissue subjected to thermo-electrical shock

Hamdy M. Youssef1
Raafat A. A. Salem2
1Mechanical Engineering Department-College of Engineering and Islamic Architecture, Umm Al, Qura University, Makkah, Saudi Arabia
2Mathematical Sciences Department, College of Applied Sciences, Umm Al, Qura University, Makkah, Saudi Arabia
Corresponding Author:
Hamdy M. Youssef
Views 193
Reads 102
Downloads 443

Abstract

The current paper is dealing with the thermal reaction and response of skin tissue subjected to a constant heat flux due to thermo-electrical shock on the bounding plane. The dual-phase-lag bio heat conduction equation based on the Tzou model has been applied and solved in the Laplace transform domain by using a direct method. The numerical inversions of the Laplace transform, and numerical solution has been obtained. The thermos-electrical shock for a small value of time has been considered with constant voltage and resistance. The temperature increment responses have been obtained and discussed with various values of voltage, resistance, and dual-lag value times. Moreover, a comparison between the well-known three models of bioheat conduction of Pennes, Ventott-Cattaneo, and Tzou has been shown in the figure. The results show that the voltage, resistance, and dual-lag value of times have significant effects on the temperature increment distribution and the thermal wave propagation through the skin tissue.

The dual-phase-lag bioheat transfer of a skin tissue subjected to thermo-electrical shock

Highlights

  • This research examined how skin tissue responds to a steady heat flux caused by a thermo-electrical shock on the bounding plane in terms of temperature reaction and response.
  • A direct technique has been used to apply and solve the dual-phase-lag bio heat conduction equation based on the Tzou model in the domain of the Laplace transform.
  • The voltage, resistance, and dual-lag time values significantly affect the distribution of temperature increments and the propagation of thermal waves through skin tissue.
  • The Tzou model is the most effective model that accurately predicts the physical behaviour of skin tissue.
  • The Tzou and Ventott-Cattaneo models also provide a thermal wave that propagates with a limited speed and is extremely similar to one another.

1. Introduction

Many treatments and their applications depend heavily on temperature transmission via live skin tissue. Microwaves, lasers, and other technological advancements have also expanded the field of bioheat transfer work [1]. Researcher Pennes, who improved the parabolic differential equation of the biological tissues for the first bioheat transfer, provides commentary on biothermal transmission analyses in real biological tissues [2]. Additionally, the consistency and nature of thermal activity in biological tissues have been modelled using the Pennes biological transmission approach. Similarities to such erratic outcomes that exhibit non-Fourier or hyperbolic conduction behaviour on the outside. To characterize the kind of hyperbolic equation, Cattaneo and Vernotte gave several modifications of heat conduction of Fourier's law in its linearity form [3, 4]. They proposed a thermal wave model to evaluate the impact of microwave and flux thermal activity. To treat diverse skin tissue disorders without harming the healthy tissue surrounding them, many sorts of treatments have been attempted. The Pennes bioheat transfer equation (PBTE) was analytically solved by Xu et al. They also studied the skin tissue bioheat transmission, skin thermal damage, tissue structure, skin biomechanics, and biothermomechanics [5] The dual-phase-lag (DPL), parabolic models of biomass transport, and hyperbolic were used to examine the non-Fourier thermomechanical activity of skin tissues in the context of various thermal loading limitations on the surface. Xu and colleagues also observed significant differences between Pennes, thermal waves, and dual-phase-lag (DPL) [6]. Additionally, Rossmann et al. analyzed how biological tissues’ perfusion, thermal characteristics, and dielectric properties change with the temperature at hyperthermic and ablation temperatures [7]. The dual-phase-lag (DPL) approach has been used to expand Tzou’s model while considering the relaxation activity at a high rate of reaction. Despite the slow procedure, the small-scale reaction is detected in time [8, 9]. Tzou introduced a new parameter time of phase lag for a temperature gradient [10]. Dual-phase-lag (DPL) heat conduction model has been used by Askarizadeh et al. to address transient heat transmission issues in skin tissue [11].

Before analyzing the general bioheat transfer paradigm in live tissue, Liu et al. [12]. The two-stage model for heat transmission issues in biological tissue uses the hyperthermia process to assess Liu et al. [13, 14]. The Liu and Chen thermal conductor model (DPL) employed non-Fourier thermal activity to characterize the thermal transit taking place in biological tissue and to diagnose hyperthermia [15].

Liu and Xu developed a closed-form analytical method to solve the Pennes bioheat equation for fluctuations in skin tissue temperature caused by sinusoidal heat flow [16]. Ahmadika et al. will have the analytical solution for the Pennes equations’ bioheat (parabolize) and hyperbolized bioheat transfer models for periodic constant and pulsed train heat flow boundary conditions [17]. Youssef and Alghamdi investigated a one-dimensional application of thermoelastic dual-phase-lag (DPL) in skin tissue under specific heat stress [18].

The body resists the flow of current. The skin is where the body experiences more than 99 % of its total electrical resistance. Ohms are used to measuring resistance. A calloused, dry hand may contain more than 100,000 due to the stratum corneum’s thick outer layer of dead cells. Due to the damp, moderately salty tissues beneath the skin, the internal body resistance is around 300. If there is skin breakdown due to high voltage, a cut, a severe abrasion, or submersion in water, the skin resistance can be efficiently overcome [19-21]. The skin behaves similarly to a capacitor or other electrical device in that it permits greater current to pass when voltage changes quickly. If one's hand is holding a metal instrument that abruptly comes into contact with a source of voltage, a quickly varying voltage will be applied to the palm and fingers. This kind of interaction will result in a substantially higher body current amplitude than would otherwise happen [19-21].

A quick, intense reaction to an electric current flowing through any portion of a person’s body is known as an electric shock. Electric shock results in electrocution, which is death. Tissue damage caused by electrical current, or voltage is referred to as primary electrical injury. Falls are a frequent source of secondary injuries. High resistance in the skin’s outer layer disintegrates above 500 V or greater. This significantly reduces the body’s resistance to current flow [19-21]. As a result, the quantity of current that flows at any given voltage increases. Skin breakdown can occasionally appear as tiny wounds the size of a pinhead that is simple to ignore. They frequently indicate that a significant quantity of current could enter the body. Muscles, nerves, and other tissues may sustain profound tissue damage because of this current. This is one of the reasons why high-voltage accidents frequently result in considerable deep tissue harm but little surface burns [19-21]. Youssef and El-Bary studied the thermal damage of the tumour tissue due to the thermoelectrical effect based on the Thomson effect. Different voltage and resistance effects have been studied and discussed in this work for three different types of tumour tissues [22].

The present paper’s goal is to build a mathematical model of bioheat skin tissue within the framework of a dual-phase-lag heat conduction model (Tzou model). The skin tissue’s surface is treated with an electrical shock with a constant voltage for a short period. The work is unique in that it investigates the impacts of dual-phase-lag time parameters, voltage value, and skin resistance on the temperature increment distribution using three distinct bioheat models: Pennes, Ventott-Cattaneo, and Tzou.

1.1. Pennes model of bioheat transfer for biological tissues

The bioheat transmission model was developed to examine the time-dependent temperature change as a heat response generated by either a heat source or thermal heating. Pennes developed the first model of biological tissues based on Fourier's law of heat conduction, which is as follows [18, 23, 24]:

1
K 2T=ρCTt-ρbwbCbTb-T-Qmet+Qext,

where ρ, C and K are the density, specific heat, and thermal conductivity of the tissue, respectively. Cb, wb, ρb and Tb denote the specific heat of the blood, blood perfusion rate, blood density, and blood temperature, respectively. T is the absolute temperature function. Qmet denotes the metabolic heat generated by the chemical reaction inside the tissue, and it is constant, and the external heat source is given by Qext. Moreover, 2 is the Laplace operator.

1.2. The Ventott-Cattaneo model of bioheat conduction (Modified Fourier law)

Ventott-Cattaneo (V-C) amended the conventional Fourier thermal drive law by postulating the assumption of limited thermal wave propagation speed, resulting in the following hyperbolic form of the thermal wave [18, 23, 24]:

2
K 2T=ρC1+τqtTt+wbCbρb1+τqtT-Tb
-1+τqtQmet+Qext,

where τq=αc02>0 is a material property and is called the relaxation time parameter, and α is the thermal diffusivity while c0 is the speed of the thermal wave inside the medium.

1.3. Dual-phase-lag model (DPL) of bioheat transfer (Tzou model)

The dual-phase-lag (DPL) model based on the dual reaction between the gradient of the temperature T and the heat flux q has been modified by Tzou to be in the following heat conduction equation [18, 23, 24]:

3
K1+τTt2T=ρC1+τqtTt+wbCbρb1+τqtT-Tb
-1+τqtQmet+Qext,

where τT0 is the second parameter of the relaxation time, which is the phase lag of the temperature gradient passing through the medium.

We assume the temperature increment as follows:

4
θ=T-Tb.

Hence, we have:

5
K1+τTt2θ=ρC1+τqtθt+wbCbρb1+τqtθ
-1+τqtQmet+Qext.

2. Basic formulation

We proposed the medium 0xL has been filled with skin tissue as in figure 1 and obeys the dual-phase-lag (DPL) as in Eq. (5) [18, 24].

Initially, the medium is dormant and has no external heat source Qext=0while the metabolic heat source Qmet is constant. The bounding plane of the outer surface of the skin tissue x=0 is subjected to constant heat flux q0, while the bounding plane of the inner surface x=L has no heat flux (see Fig. 1).

Hence, the heat conduction equations take the forms:

6
τTt+12θx2=τqρCK2t2+ρC+τqwbCbρbKt+wbCbρbKθ-QmetK.

Fig. 1The skin tissue subjected to electrical shock with constant voltage

The skin tissue subjected to electrical shock with constant voltage

The initial conditions are:

7
θx,tt=0=θx,t tt=0=0.

The boundary conditions are:

8
θx,t xx=0=-q0K, θx,t xx=L=0.

The solution in the Laplace transform domain:

9
F-x,s=Lfx,t=0fx,te-stdt , t0,

where the inverse of the Laplace transform takes the form:

10
L-1f-s=fteκtt12f-κ+Ren=1N-1nf-κ+i nπt,

where Re means the real part, while i is the imaginary number unit, for faster convergence, many experiments confirmed that the value of the parameter κ could satisfy the relation κt 4.7 [8, 25].

Then, we have:

11
d2θ-x,sdx2-λ2sθ-x,s=-γs,

where:

12
λ2s=τqρCs2+ρC+τqwbCbρbs+wbCbρbτTs+1K, γs=Qmets K,
dθ-x,sd xx=0=-q0sK, dθ-x,sd xx=L=0.

The skin tissue is subjected to thermal shock that comes from the thermal effect due to electrical shock with constant voltage (Joule’s equation of electrical heating) VV and constant resistance ReΩ when the current has a vibration of 60 Hz, the mean value of the skin resistance will be considered in the range 300-400 Ω [7, 19-22].

Then, Joule’s equation of electrical heating is given by [7, 19, 21, 22]:

13
q0=V2t0Re,

where t0t0>0 is a parameter which gives the time interval of the electrical shock.

The general solution of the Eq. (11) is given by:

14
θ-x,s=θ1se-λx+θ2seλx+γsλ2s,

where θ1s and θ2s are some parameters which must be determined by using the boundary condition Eq. (12).

Hence, we have:

15
-θ1+θ2=-α, -θ1e-λL+θ2eλL=0,

where α=q0/sKλ.

By solving the above system, we obtain:

16
θ1=αeλLeλL-e-λL, θ2=αe-λLeλL-e-λL.

Then, we have the solution in the form:

17
θ-x,s=αcoshλL-xsinhLλ+γλ2,

which completes the solution in the Laplace transform domain where the solution in the physical time domain can be obtained by applying iteration Eq. (10).

3. Numerical results

The temperature distribution through skin tissue is explored in this work using three bioheat transfer models (Pennes, Ventott-Cattaneo, and Tzou) under constant heat flow circumstances on the skin surface. Table 1 shows the values of the essential thermal parameters employed in the current computations [5, 18, 23, 24, 26].

Table 1The material properties of the skin tissue

Parameter
Value
Unit
Parameter
Value
Unit
K
0.215
W/m°C
Wb
0.00052
ml/C m
ρ
1000
kg/m3
Tb
37
°C
ρb
1060
kg/m3
τT,τq
20.0
s
C
4187
J/kg°C
t0
1.0
s
Cb
3800
J/kg°C
L
0.01
m

4. Discussions

Fig. 2 shows the temperature increment distribution for variance values of voltage V= (100, 110, 120) V when t= 120 s, Re= 300 Ω, τT=τq= 20 s, and t0= 1.0 s. It is noted that the value of the voltage of the thermo-electrical shock has a significant effect on the thermal wave. The value of the temperature increment increases when the value of the voltage increases. The values of the temperature increment on the boundary of the skin tissue x=0 take the following order:

18
θV=100 Vx=0<θV=110 Vx=0<θV=120 Vx=0.

Fig. 3 shows the temperature increment distribution for variance values of skin tissue resistance Re= (300, 350, 400) Ω when t= 120 s, V= 100 V, τT=τq= 20 s, and t0= 1.0 s. It is noted that the value of the voltage of the skin tissue resistance has a significant effect on the thermal wave. The value of the temperature increment decreases when the value of the skin tissue resistance increases. The values of the temperature increment on the boundary of the skin tissue x=0 take the following order:

19
θRe=400 Ωx=0<θRe=350 Ωx=0<θRe=300 Ωx=0.

Fig. 2The temperature increment for variance values of voltage when t= 120 s, Re= 300 Ω, τT=τq= 20 s, and t0= 1.0 s.

The temperature increment for variance values of voltage  when t= 120 s, Re= 300 Ω, τT=τq= 20 s, and t0= 1.0 s.

Fig. 3The temperature increment for variance values of resistance when t= 120 s, V= 100 V, τT=τq= 20 s, and t0= 1.0 s

The temperature increment for variance values of resistance  when t= 120 s, V= 100 V, τT=τq= 20 s, and t0= 1.0 s

Fig. 4 shows the temperature increment distribution for variance values of electrical shock time t0= (1.0, 1.5, 2.0) s when t= 120 s, V= 100 V, τT=τq= 20 s, and Re= 300 Ω. It is noted that the value of the electrical shock time has a significant effect on the thermal wave. The value of the temperature increment increases when the value of the electrical shock time increases. The values of the temperature increment on the boundary of the skin tissue x=0 take the following order:

20
θt0=1.0 sx=0<θt0=1.5 sx=0<θt0=2.0 sx=0.

Fig. 5 shows the temperature increment distribution for variance situations of the dual-lag times parameters τq<τT,τq=τT,τq>τT when t= 120 s, V= 100 V, t0= 1.0 s and Re= 300 Ω. It is noted that the value of the dual-lag times parameters has a significant effect on the thermal wave. The values of the temperature increment on the boundary of the skin tissue x=0 take the following order:

21
θτT<τqx=0<θτT=τqx=0<θτT>τqx=0.

We can see that the situation when the value of the lag-time of the temperature is smaller than the value of the lag-time parameter of the temperature gradient, generates a thermal wave which has a finite speed of propagation where in this case, the temperature increment goes to zero before the other two cases.

Fig. 4The temperature increment for variance values of the electrical shock time when t= 120 s, V= 100 V, τT=τq= 20 s, and Re= 300 Ω

The temperature increment for variance values of the electrical shock time  when t= 120 s, V= 100 V, τT=τq= 20 s, and Re= 300 Ω

Fig. 5The temperature increment for variance situations of the dual-lag times parameters when t= 120 s, V= 100 V, t0= 1.0 s and Re= 300 Ω

The temperature increment for variance situations of the dual-lag times parameters  when t= 120 s, V= 100 V, t0= 1.0 s and Re= 300 Ω

Fig. 6 shows the temperature increment distribution for Pennes, Ventott-Cattaneo, and Tzou models when t= 120 s, V= 100 V, t0= 1.0 s, and Re= 300 Ω. It is noted that the temperature increment distribution has different values based on the three different models. The values of the temperature increment take the following order:

22
θTzou<θVentott-Cattaneo<θPennes.

The results of this figure proved that the Tzou model is the most successful model which agrees with the physical behaviour of the skin tissue. Moreover, Tzou and Ventott-Cattaneo models are very close and offer a thermal wave propagating with a finite speed.

Fig. 6The temperature increment for Pennes, Ventott-Cattaneo, and Tzou models when t= 120 s, V= 100 V, t0= 1.0 s, and Re= 300 Ω

The temperature increment for Pennes, Ventott-Cattaneo, and Tzou models  when t= 120 s, V= 100 V, t0= 1.0 s, and Re= 300 Ω

5. Conclusions

This research examined how skin tissue responds to a steady heat flux caused by a thermo-electrical shock on the bounding plane in terms of temperature reaction and response. A direct technique has been used to apply and solve the dual-phase-lag bio heat conduction equation based on the Tzou model in the domain of the Laplace transform. The Laplace transforms’ numerical inversions and solutions have been discovered. With constant voltage and resistance, the thermos-electrical shock for a short period has been explored.

With different voltage, resistance, and dual-lag time values, the temperature increment responses have been measured and discussed. Additionally, a comparison of the well-known Pennes, Ventott-Cattaneo, and Tzou bioheat conduction models has been illustrated.

The findings demonstrate that the voltage, resistance, and dual-lag time values significantly affect the distribution of temperature increments and the propagation of thermal waves through skin tissue.

The outcomes of this research demonstrated that the Tzou model is the most effective model that accurately predicts the physical behaviour of skin tissue. The Tzou and Ventott-Cattaneo models also provide a thermal wave that propagates with a limited speed and is extremely similar to one another.

References

  • F. Xu, T. J. Lu, K. A. Seffen, and E. Y. K. Ng, “Mathematical modeling of skin bioheat transfer,” Applied Mechanics Reviews, Vol. 62, No. 5, p. 05080, Sep. 2009, https://doi.org/10.1115/1.3124646
  • H. H. Pennes, “Analysis of Tissue and Arterial Blood Temperatures in the Resting Human Forearm,” Journal of Applied Physiology, Vol. 1, No. 2, pp. 93–122, Aug. 1948, https://doi.org/10.1152/jappl.1948.1.2.93
  • C. Cattaneo, “A form of heat-conduction equations which eliminates the paradox of instantaneous propagation,” Comptes Rendus, Vol. 247, p. 431, 1958.
  • P. Vernotte, “Les paradoxes de la theorie continue de l’equation de la chaleur,” in Comptes Rendus, Vol. 246, p. 3154, 1958.
  • F. Xu, T. J. Lu, and K. A. Seffen, “Biothermomechanical behavior of skin tissue,” Acta Mechanica Sinica, Vol. 24, No. 1, pp. 1–23, Feb. 2008, https://doi.org/10.1007/s10409-007-0128-8
  • F. Xu, K. A. Seffen, and T. J. Lu, “Non-Fourier analysis of skin biothermomechanics,” International Journal of Heat and Mass Transfer, Vol. 51, No. 9-10, pp. 2237–2259, May 2008, https://doi.org/10.1016/j.ijheatmasstransfer.2007.10.024
  • C. Rossmanna and D. Haemmerich, “Review of temperature dependence of thermal properties, dielectric properties, and perfusion of biological tissues at hyperthermic and ablation temperatures,” Critical Reviews in Biomedical Engineering, Vol. 42, No. 6, pp. 467–492, 2014, https://doi.org/10.1615/critrevbiomedeng.2015012486
  • D. Y. Tzou, “A unified field approach for heat conduction from macro – to micro-scales,” Journal of Heat Transfer, Vol. 117, No. 1, pp. 8–16, Feb. 1995, https://doi.org/10.1115/1.2822329
  • D. Jou, “DY Tzou: macro-to microscale heat transfer. the lagging behaviour. series in chemical and mechanical engineering,” Journal of Non Equilibrium Thermodynamics, Vol. 23, No. 2, pp. 192–193, 1998.
  • M. N. O. ¨Zis¸Ik and D. Y. Tzou, “On the wave theory in heat conduction,” Journal of Heat Transfer, Vol. 116, No. 3, pp. 526–535, Aug. 1994, https://doi.org/10.1115/1.2910903
  • H. Askarizadeh and H. Ahmadikia, “Analytical analysis of the dual-phase-lag model of bioheat transfer equation during transient heating of skin tissue,” Heat and Mass Transfer, Vol. 50, No. 12, pp. 1673–1684, Dec. 2014, https://doi.org/10.1007/s00231-014-1373-6
  • K.-C. Liu, Y.-N. Wang, and Y.-S. Chen, “Investigation on the bio-heat transfer with the dual-phase-lag effect,” International Journal of Thermal Sciences, Vol. 58, pp. 29–35, Aug. 2012, https://doi.org/10.1016/j.ijthermalsci.2012.02.026
  • K.-C. Liu and H.-T. Chen, “Analysis for the dual-phase-lag bio-heat transfer during magnetic hyperthermia treatment,” International Journal of Heat and Mass Transfer, Vol. 52, No. 5-6, pp. 1185–1192, Feb. 2009, https://doi.org/10.1016/j.ijheatmasstransfer.2008.08.025
  • K.-C. Liu, H.-T. Chen, and P.-J. Cheng, “Inverse investigation of non-Fourier heat conduction in tissue,” Journal of Thermal Biology, Vol. 62, pp. 123–128, Dec. 2016, https://doi.org/10.1016/j.jtherbio.2016.07.002
  • K.-C. Liu and H.-T. Chen, “Investigation for the dual phase lag behavior of bio-heat transfer,” International Journal of Thermal Sciences, Vol. 49, No. 7, pp. 1138–1146, Jul. 2010, https://doi.org/10.1016/j.ijthermalsci.2010.02.007
  • Jing Liu and L. X. Xu, “Estimation of blood perfusion using phase shift in temperature response to sinusoidal heating at the skin surface,” IEEE Transactions on Biomedical Engineering, Vol. 46, No. 9, pp. 1037–1043, 1999, https://doi.org/10.1109/10.784134
  • H. Ahmadikia, R. Fazlali, and A. Moradi, “Analytical solution of the parabolic and hyperbolic heat transfer equations with constant and transient heat flux conditions on skin tissue,” International Communications in Heat and Mass Transfer, Vol. 39, No. 1, pp. 121–130, Jan. 2012, https://doi.org/10.1016/j.icheatmasstransfer.2011.09.016
  • H. M. Youssef and N. A. Alghamdi, “Modeling of one-dimensional thermoelastic dual-phase-lag skin tissue subjected to different types of thermal loading,” Scientific Reports, Vol. 10, No. 1, pp. 1–12, Dec. 2020, https://doi.org/10.1038/s41598-020-60342-6
  • R. M. Fish and L. A. Geddes, “Conduction of electrical current to and through the human body: a review,” Eplasty, Vol. 9, 2009.
  • R. M. Fish and L. A. Geddes, “Electrophysiology of connection current spikes,” Cardiovascular Engineering, Vol. 8, No. 4, pp. 219–224, Dec. 2008, https://doi.org/10.1007/s10558-008-9059-7
  • S. Grimnes, “Dielectric breakdown of human skinin vivo,” Medical and Biological Engineering and Computing, Vol. 21, No. 3, pp. 379–381, May 1983, https://doi.org/10.1007/bf02478510
  • H. M. Youssef and A. A. El-Bary, “Voltage and time required for irreversible thermal damage of tumor tissues during electrochemotherapy under thomson effect,” Mathematics, Vol. 8, No. 9, p. 1488, Sep. 2020, https://doi.org/10.3390/math8091488
  • H. M. Youssef and N. A. Alghamdi, “Characterization of thermal damage due to two-temperature high-order thermal lagging in a three-dimensional biological tissue subjected to a rectangular laser pulse,” Polymers, Vol. 12, No. 4, p. 922, Apr. 2020, https://doi.org/10.3390/polym12040922
  • H. M. Youssef and N. A. Alghamdi, “Three-dimensional biological tissue under high-order effect of two-temperature thermal lagging to thermal responses due to a laser irradiation,” 36th International Conference on Vibroengineering in Dubai, United Arab Emirates, March 15-17th, 2019, Vol. 22, pp. 112–117, Mar. 2019, https://doi.org/10.21595/vp.2019.20521
  • D. Tzou and M.-T. M. H. Transfer, The Lagging Behavior. Washington: Taylor Francis, 1997.
  • Y. Zhang, “Generalized dual-phase lag bioheat equations based on nonequilibrium heat transfer in living biological tissues,” International Journal of Heat and Mass Transfer, Vol. 52, No. 21-22, pp. 4829–4834, Oct. 2009, https://doi.org/10.1016/j.ijheatmasstransfer.2009.06.007

About this article

Received
21 September 2022
Accepted
20 October 2022
Published
25 October 2022
Keywords
phase-lag
skin tissue
heat flux
Fourier series
thermo-electrical effect
Acknowledgements

The authors have not disclosed any funding.

Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflict of interest

The authors declare that they have no conflict of interest.