Next Article in Journal
Comparative Analysis of the Chaotic Behavior of a Five-Dimensional Fractional Hyperchaotic System with Constant and Variable Order
Next Article in Special Issue
Artificial Intelligence in Chromatin Analysis: A Random Forest Model Enhanced by Fractal and Wavelet Features
Previous Article in Journal
Enhancing Transient Stability in Multi-Machine Power Systems through a Model-Free Fractional-Order Excitation Stabilizer
Previous Article in Special Issue
On the Application of Mann-Iterative Scheme with h-Convexity in the Generation of Fractals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractal Operators Abstracted from Arterial Blood Flow

1
Department of Engineering Mechanics, Tsinghua University, Beijing 100084, China
2
Department of Civil Engineering, Chengdu Technological University, Chengdu 611730, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(7), 420; https://doi.org/10.3390/fractalfract8070420
Submission received: 11 May 2024 / Revised: 16 July 2024 / Accepted: 16 July 2024 / Published: 18 July 2024
(This article belongs to the Special Issue Fractals in Biophysics and Their Applications)

Abstract

:
In this paper, we present a re-established functional fractal circuit model of arterial blood flow that incorporates the shunt effect of the branch vessels. Under the background of hemodynamics, we abstracted a family of fractal operators and investigate the kernel function and properties thereof. Based on fractal operators, the intrinsic relation between Bessel function and Struve function was revealed, and some new special functions were found. The results provide mathematical tools for biomechanics and automatic control.

1. Introduction

The functional circuit model (FCM) [1,2] is extensively employed in the field of hemodynamics. The arterial lumen can be considered as an elastic cavity with periodic contraction and dilation, known as the Windkessel model, introduced by Frank [2]. Based on the electrical-mechanical analogy, the blood flow movement within the elastic cavity can be simulated by a functional circuit model composed of basic electrical components. Frank’s functional circuit model only involves a single elastic cavity and two components (a capacitor and a resistor), with limited simulation capabilities. Goldwyn et al. [3] proposed a multi-level elastic chamber model that incorporates the blood acceleration in space. Burattini et al. [4] and Stergiopulos et al. [5] introduced the inductive component to simulate the inertia of blood. With the increasement of the variety of electrical components and the number of cavities, more and more complex circuits were used to simulate blood flow movements more accurately [6,7,8,9,10,11,12,13,14].
Generally, the utilization of fluid dynamics methods based on the Navier–Stokes equations for simulating blood flow demands a substantial amount of computational resources and time due to high computation costs. In contrast, Windkessel models and FCM, which only require solving circuit equations without dealing with complex models and equations, can effectively simulate the pressure-flow relationship of hemodynamics at a relatively low computational cost, so they have been widely utilized and applied in medical diagnostics [15,16].
The concepts of fractal cells and fractal operators have been developed in recent years and successfully applied to various disciplines such as tendon and ligament, nerve, bone and hemodynamics [17,18,19,20]. The results indicate a close correlation between functional fractals and fractional dynamics. Based on the constitutive established by functional fractals, both fractional-order operators and fractional-order responses can be derived [21,22], which significantly differs from merely introducing fractional-order elements phenomenally [23,24,25].
In 2020, through the infinite-level division of the elastic cavities, Peng et al. [20] achieved the fractal of blood flow image and designed the functional fractal circuit of aorta and arterioles, developing the hemodynamics in the functional fractal space. Recently, we realized that Peng’s model ignored the shunt effect of branch vessels. Therefore, this paper enhances the model, establishing a new functional fractal circuit, and a family of fractal operators are abstracted from the new model. The properties of the operator family are studied. The study of operator properties is not limited to hemodynamics. In fact, since the circuits we have given are constructible with the components and parameters designable, the properties of the operators we have studied can serve as a foundation for automatic control and a better description of biological phenomena.
The present paper primarily encompasses the subsequent sections. In Section 2, the functional fractal circuit is established. In Section 3, the operator algebraic equations are established and solved. In Section 4 and Section 5, the kernel functions of the fractal operator family are investigated.

2. Fractal Circuit Model

In the classical elastic cavity model (Windkessel model), the blood pressure P is compared to the voltage U, and the flow rate Q is compared to the current I. For the aorta, the longitudinal inertia of the blood flow is compared to the inductance L, causing the longitudinal flow; the wall compliance is compared to the capacitance C, causing the transverse flow. For the arteriole, transverse flow is also caused by the wall compliance, while the longitudinal flow is dominated by the blood viscosity, compared to the resistance R. the relationship between longitudinal and transverse flow can be likened to that of two branches in parallel of a circuit.
The comparative relationship between hemodynamics and the electric parameters is shown in Table 1. c L and c R are constant coefficients. r is the vessel radius, and h is the thickness of the vessel.
Peng et al. [20] differentiated the classical single segment elastic cavity into a combination of infinite micro-elastic ones, which is inherited by this article as a basic thought. However, we notice that both the single elastic cavity model and Peng’s infinite micro-elastic cavity model took no account of the branch vessels and therefore ignored the shunt effect of the branch vessels.
Nevertheless, the real artery has a multilevel branching topology, called “arterial vascular tree” in textbooks. Taking the aorta as an example, due to its numerous branch vessels, blood flowing through the aortic lumen continuously exits along these branches, forming tributaries. In this case, if we establish the elastic cavity model of the aorta, the shunt effect of tributaries should be taken into account.
When the shunt vessels are taken into account, the basic flow pattern of arterial blood flow can be described as follows, as depicted in Figure 1.
When the aortic wall undergoes contraction and dilation (indicated by black arrows in Figure 1), the blood flow within the vascular lumen exhibits three distinct patterns of flow (indicated by red arrows in Figure 1): axial (longitudinal) flow, radial (transverse) flow, and lateral branch flow. The axial flow is normal, that is, irrespective of whether the vessel wall dilates or contracts, axial flow always exists and maintains a consistent direction towards the distal end. There is a cooperative relationship between radial flow and branch flow, both driven by the motion of the vessel wall.
When the vessel wall contracts (indicated by the black arrows in Figure 1a), radial blood flow is directed from the wall towards the center (indicated by the red arrows in Figure 1a). Simultaneously, lateral branch flow is initiated, causing intracavity blood to exit the main cavity through branch blood vessels (indicated by the red arrows in Figure 1a). When the wall undergoes dilation (indicated by the black arrows in Figure 1b), radial blood flow is directed from the central region towards the wall (indicated by the red arrows in Figure 1b). Simultaneously, closure of the branch vessel channel occurs, effectively arresting branch flow and preventing any backflow of blood from the branch vessel into the main cavity.
According to the above flow patterns, the principal topological characteristics of arterial flow images can be summarized: for each micro-elastic cavity unit, the longitudinal (axial) flow in the main cavity is in parallel with the transverse (radial) flow; the blood flow in branch blood vessels runs parallel to that in the main vessel, and the branch blood flow of each micro-elastic cavity unit leads to infinity and no longer flows into the main blood vessel. By combining an infinite number of micro-elastic cavities, we can construct the functional fractal circuit shown in Figure 2. In Figure 2, M 1 , M 2 , and M 3 are basic electrical components, which could be inductors, capacitors, or resistors.
It could be seen that the functional fractal circuit in Figure 2 fully conforms to the mentioned topological characteristics of the blood flow image. Specifically, the parallel operation of M 1 and M 2 mirrors the parallel operation of the longitudinal and transverse blood flow; the branches of M 3 extend infinitely, corresponding to the branching outflow towards infinity, thereby becoming disconnected with the main vessel’s blood flow.
Moreover, the functional fractal circuit in Figure 2 has the property of infinite self-similarity. The right side of Figure 2 gives the smallest repeating unit in the circuit, called the fractal cell. While classical fractals exhibit “self-similarity” [27,28], the functional fractal circuit demonstrates “self-congruence”.
It is noted that the flow patterns depicted in Figure 1 are universally observed in living organisms. Blood flow, lymph flow, and fluid flow in tissues can all be reduced to the movement patterns shown in Figure 1. In other words, the functional fractal shown in Figure 2 represents a spatial pattern with universal significance in biological systems. Therefore, the primary objective of this paper is not to investigate hemodynamics, but rather to examine the fundamental characteristics of functional fractal spaces utilizing blood flow images as a contextual backdrop. This approach offers a fresh perspective and establishes a novel mathematical foundation for studying fluid dynamics within living organisms.

3. Operator Algebraic Equations of the Fractal Circuit Model and Their Radical Solutions

According to Mikusinski’s work [29], the basic differential and integral operators are defined as follows. If the function f ( t ) has continuous derivatives at t 0 , then the definition of the differential operator acting on the function is
p f ( t ) = d f ( t ) d t + f ( 0 ) δ ( t ) ,
where δ ( t ) is the Dirac pulse function. The definition of integral operator l is that
l f ( t ) = 0 t f ( τ ) d τ .
The definitions of the basic differential and integral operators p and l in this paper strictly adhere to the explicit forms presented in Mikusinski’s work [29]. Based on basic differential operators and integral operators, Mikusinski constructed the operator domain, which is referred to as the Mikusinski operator field in this paper, and its internal structure has been described in detail in [29,30]. Mikusinski employed rigorous mathematical logic within his book to prove that the basic differential operator p and integral operator l defined by (1) and (2) are mutually inverse, also demonstrating that their combined operation yields results independent of the order of action. In the Mikusinski operator field, Equation (3) is valid, which serves as a crucial logical foundation for operator algebra theory.
l = 1 p .
Based on Equation (3), both differential operators and integral operators can be regarded as functions of the basic differential operator p defined in Equation (1).
Similar to Peng et al. [20], the admittance operator T is used to denote the map from voltage u ( t ) to current i ( t ) , that is:
i ( t ) = T u ( t ) ,
The concept of admittance is widely utilized in the field of electricity, and it serves as the reciprocal of impedance. Each electrical component has its corresponding admittance operator. Admittance operators for the three types of basic electrical components are presented in Table 2.
The three admittance operators listed in Table 2 are all basic operators in the Mikusinski operator field.
In the functional fractal circuit shown in Figure 2, the admittance operators corresponding to the electrical components M 1 , M 2 , and M 3 are denoted as T 1 , T 2 , and T 3 , respectively. The admittance operator of the fractal element M is denoted as T, referred to as the fractal operator. T 1 , T 2 , and T 3 are determined based on the selection of M 1 , M 2 , and M 3 , respectively, according to Table 2. Therefore, T 1 , T 2 , and T 3 are operators in the Mikusinski operator field, which is a prerequisite for performing following algebraic operations on them.
The fractal cell in Figure 2 is equivalent to the entire fractal circuit, and the fractal circuit entirety is equivalent to a fractal element. Thus, it is immediately deduced that the fractal cell is equivalent to a fractal element. For the convenience of discussion, the equivalence is depicted in Figure 3, the admittance of each component or element denoted in parentheses.
Based on such equivalence and combined with operator algebra theory [29,30,31,32], the operator algebraic equation can be established:
( T 1 + T 2 ) T T 1 + T 2 + T + T 3 = T .
The process of deriving Equation (5) is as follows: The parallel combination of M 1 and M 2 can be treated as a composite component with an admittance of ( T 1 + T 2 ) . This composite component is then in series with the fractal element M, whose admittance is T, resulting in a total admittance of ( T 1 + T 2 ) T ( T 1 + T 2 ) + T . Subsequently, the entire branch is in parallel with M 3 , and their admittances are summed to obtain the overall admittance of the fractal cell, that is the left side of Equation (5). As the fractal cell is equivalent to the fractal element, their admittances are equal, thus ensuring that Equation (5) holds.
By rearranging Equation (5) we can obtain
T 2 T 3 T T 3 ( T 1 + T 2 ) = 0 .
Equation (6) is the algebraic equation of the fractal operator T. It is a quadratic equation with one unknown, which has the radical solution
T = T 3 ± T 3 2 + 4 T 3 ( T 1 + T 2 ) 2 .
Since T 1 , T 2 , and T 3 are all functions of the basic differential operator p defined by Equation (1), the fractal operator T is also a function of p, i.e.,
T = T ( p ) .
As introduced at the beginning of Section 2, the arteries are divided into the aorta and the arteriole, whose longitudinal flows are, respectively, dominated by inertia and viscosity. According to Table 1, the longitudinal flow of the aorta and the arteries are respectively denoted by the inductive component L and the resistive component R. The elasticity of the artery wall is represented by the capacitive component C. Therefore, as a direct application, M 1 and M 3 should be selected from the inductive L and the resistive R according to the type of main and branch vessels described, i.e., T 1 and T 3 should be selected from T L and T R ; M 2 should be taken as the capacitive C, i.e., T 2 should be taken as T C .
However, the discussion in this paper is not limited to this. The functional fractal circuit shown in Figure 2 provides a completely different topological structure than that of previous models, and the topological structure is universal in living bodies. Therefore, we hope to leverage the topological structure and optimize the control characteristics through different element matching methods, thereby obtaining valuable control results for the biomechanics and automatic control, not limited to hemodynamics. In the scope of this paper, the types of basic components are finite, so we can take an ergodic approach and discuss all possible cases where T 1 , T 2 and T 3 take on admittance operators of various basic electronic components. These cases are constructible, and the components and parameter range are designable. See the analysis below for details.
Notice that there is a radical ( ) in the fractal operator. In the previous articles [21,22], we referred to the fractal operators with radical ( ) forms as fractional operators of apparent 1/2 order. Apparently, they are irrational fractional operators.
The fractal operator mentioned in Equation (7) is rooted in hemodynamics. It should be emphasized that the subjects and physical mechanisms studied in hemodynamics, ligament fiber mechanics, and nerve signal release dynamics are fundamentally distinct [33,34,35,36,37]. However, disciplines with such great differences have similar functional fractals, similar fractal operator algebraic equations, and similar radical operator solutions for algebraic equations. This means that it is imperative to study the radical operators uniformly.
The fractal operators in Equation (7) constitute a “family of fractal operators” of radical type. This raises a question: If we abandon the hemodynamic background and focus solely on investigating the fractal operator family, can we obtain results that hold universal significance? Specifically, we can ask the following question: How many basic or independent operators are there in the family of fractal operators? What are their properties? What are the similarities and differences between these properties? Clearly, by clarifying these questions, the toolbox of biomechanics and control theory can be greatly enriched.
There are three types of basic electrical components, resistance, inductance, and capacitance, as shown in Table 2. In Equation (7), there are exactly three operators, T 1 , T 2 , and T 3 . Therefore, T 1 , T 2 , and T 3 can be selected as admittance operators of the three basic components respectively. By combining the three basic admittance operators, nine fractal operators can be obtained, as shown in Table 3.
In Table 3, the powers of the basic differential operator p in the radical ( ) , from −2 to 2, all appear. It can be said that Table 3 gives the basic graph of the 1/2 order fractal operator.
All constants in Table 3 have specific physical significance: RC, LC, L R , determining the characteristic frequency or characteristic time of a functional fractal system. In order to ensure the universality of the research results, we discard the specific meanings of the physical parameters, and then abstract the operators with radical ( ) into the following more general form: p 2 + α 2 , 1 p p 2 + α 2 ; p + α , 1 p p + α , 1 p p + α ; ( p + β ) 2 α 2 , 1 p ( p + β ) 2 α 2 , 1 p ( p + β ) 2 α 2 . Without special instructions, we set α , β > 0 . α and β are physical constants that generally depends on the characteristic frequency or characteristic time.

4. Definition of Kernel Functions of Fractal Operators

The convolution of two functions f ( t ) and g ( t ) are defined as follows.
f ( t ) g ( t ) 0 t f ( t τ ) g ( τ ) d τ .
According to Mikusinski’s operator theory [29], as functions of the basic differential operator p, operators T ( p ) can be defined by their kernel functions K ( t ) . Specifically, the action of operator T ( p ) on a function f ( t ) is equal to the convolution between the kernel function K ( t ) of the operator T ( p ) and the function f ( t ) , that is:
T ( p ) f ( t ) = K ( t ) f ( t ) = 0 t K ( t τ ) f ( τ ) d τ .
If we take f ( t ) as the Dirac pulse function δ ( t ) from Equation (10), we can derive Equation (11).
T ( p ) δ ( t ) = K ( t ) .
Equation (11) means that, the kernel function (also called the integral kernel) of operator T ( p ) is equal to the result of the action of T ( p ) on the Dirac pulse function δ ( t ) . As a direct definition of the kernel function, Equation (10) repeatedly is used in the remain part of the article.
Mikusinski did not provide a universal method for explicitly obtaining the kernel function of a general operator. Yu et al. [30], based on an analysis of the properties of the Minkusinsiki operator field and kernel function space, rigorously derived the following conclusion: when the definition of the basic differential operator p follows Equation (1), the kernel function K ( t ) of the operator T ( p ) can be obtained by taking the Laplace inverse transform of the operator function T ( p ) , that is:
1 [ T ( p ) ] = K ( t ) ,
i.e.,
T ( p ) δ ( t ) = K ( t ) = 1 [ T ( p ) ] ,
Based on the above theory, the kernel function of the fractal operator family can be clarified.
To solve differential equations, the Laplace transform on the differential equation is a well-known operation. Yu’s work [30] rigorously demonstrated that within Mikusinski’s operator algebra system, the integral kernel of an operator can be obtained by directly applying the Laplace inverse transform to the operator itself, eliminating the need for performing both Laplace transform and inverse transform on the entire differential equation. In fact, the operator solely encompasses the intrinsic properties of the system, whereas the differential equation incorporates both intrinsic information and external excitation (as depicted in (4)). Transforming the operator directly to obtain the integral kernel is more conducive to comprehending the inherent structure and properties of the system, which is exactly what we want in this paper.

5. Kernel Functions of the Fractal Operator Family

The present section focuses on the kernel functions of the exemplary fractal operators listed in Table 3.

5.1. Kernel Functions of the Fractal Operators p 2 + α 2 and 1 p p 2 + α 2

Peng et al. [20] used to derive the operator p 2 + α 2 p in hemodynamics. Mikusinski also discussed the operator in his work [29]. The inverse Laplace transformation of the fractal operator p 2 + α 2 p is:
1 ( p 2 + α 2 p ) = α t J 1 ( α t ) ,
where J 1 is the Bessel function of the first kind. As a direct application of Equation (12), the right side of Equation (14) gives the kernel function of the fractal operator p 2 + α 2 p . Thus, according to Equation (11) we have
[ p 2 + α 2 p ] δ ( t ) = α t J 1 ( α t ) .
The kernel function results given by Equations (14) and (15) are consistent with those obtained by Mikusinski [20] using a series expansion. Note that p 2 + α 2 p is an irrational operator; therefore, its kernel is a non-elementary function, specifically a weighted Bessel function α t J 1 ( α t ) .
Based on Equation (15), we can analyze the operator 1 p p 2 + α 2 . Equation (13), combined with the inverse Laplace transformation, derives that
1 ( p 2 + α 2 p p ) = 0 t α τ J 1 ( α τ ) d τ ,
So, the kernel function is
[ p 2 + α 2 p p ] δ ( t ) = 0 t α τ J 1 ( α τ ) d τ = α 2 { 2 [ α t J 0 ( α t ) J 1 ( α t ) ] + π α t [ J 1 ( α t ) H 0 ( α t ) J 0 ( α t ) H 1 ( α t ) ] } .
The proof of the second equation of Equation (17) can be seen in Appendix A. H 0 and H 1 are zero-order and first-order Struve functions, respectively. p 2 + α 2 p p is an irrational operator, and its kernel function is a non-elementary function 0 t α τ J 1 ( α τ ) d τ . Curiously, the non-elementary function can be characterized by the Bessel function and the Struve function.
Figure 4 involves Bessel function J 1 ( t ) , weighted Bessel function J 1 ( t ) t , and its integral 0 t 1 τ J 1 ( τ ) d τ . The graphs in Figure 4 demonstrate that all the three types of functions exhibit volatility, which evidently comes from the inherent properties of the Bessel function.
As time tends to infinity, all three types of functions tend to converge towards a constant value, that is:
lim t J 1 ( t ) = lim t J 1 ( t ) t = 0 .
lim t 0 t 1 τ J 1 ( τ ) d τ = 1 .
Compared with the Bessel function J 1 ( t ) , the weighted Bessel function J 1 ( t ) t and its integral 0 t 1 τ J 1 ( τ ) d τ exhibit a shorter convergence time to a constant value. Such superior convergence properties hold potential applications in the fields of biomechanics, biomimetic mechanics and automatic control.
Equation (17) is instructive. If the operator term on the left side of Equation (16) is discarded and only the second equation remains, it can be obtained that (see Appendix A for detailed proof):
0 t 1 τ J 1 ( α τ ) d τ = 1 2 { 2 [ α t J 0 ( α t ) J 1 ( α t ) ] + π α t [ J 1 ( α t ) H 0 ( α t ) J 0 ( α t ) H 1 ( α t ) ] } .
The Bessel function and the Struve function are associated by Equation (19). The Struve function arises from the nonhomogeneous Bessel equation given below [38]:
x 2 d 2 y d x 2 + x d y d x + ( x 2 m 2 ) y = 4 ( x 2 ) m + 1 π Γ ( m + 1 2 ) .
A set of solutions to the Equation (21) is
H m ( z ) = ( z 2 ) m + 1 n = 0 ( 1 ) n ( z 2 ) 2 n Γ ( n + 3 2 ) Γ ( n + m + 3 2 ) ,
where H m ( z ) is the Struve function. Thus, Equation (19) is not particularly unexpected. Although there is a correlation between the Bessel function and the Struve function, what is really surprising is that they are correlated in such a complex way. It should be emphasized that they are very different special functions. Obviously, the irrational fractional operator p 2 + α 2 p p serves as an exceptional bridge, connecting two distinct special functions in a remarkably ingenious manner.
The fractal operators have a clear biomechanical background, indicating that it is insufficient to describe the mechanical behavior of biological materials or structures solely through elementary functions. In practice, it is difficult to characterize the long-term behavior of organisms using elementary functions alone; non-elementary functions are necessary in such cases [39,40]. To a certain extent, it could be said that biomechanics is mechanics characterized by fractal operators and non-elementary functions.
As depicted in Figure 4, the kernel function of the fractal operator p 2 + α 2 p p converges very fast and almost converges to a stable value after two oscillations. These desirable characteristics align perfectly with the requirements of automatic control theory. In particular, the parameters in the operator p 2 + α 2 p p are designable. In other words, by designing the component parameters in the fractal circuit, we can obtain useful control characteristics.

5.2. Kernel Functions of the Fractal Operators p 2 α 2 and 1 p p 2 α 2

Although operators p 2 α 2 and 1 p p 2 α 2 do not appear in Table 3, they are deeply intrinsically related to the fractal operators p 2 + α 2 and 1 p p 2 + α 2 in the above section, so they are also included in the family of radical type operators.
Introducing the fractal operator p 2 α 2 p , there are two ways to obtain its kernel function. The first path is the inverse Laplace transformation. The inverse Laplace transformation of operator p 2 α 2 p gives that
1 ( p 2 α 2 p ) = α t J ¯ 1 ( α t ) .
The function α t J ¯ 1 ( α t ) at the right side of Equation (23) is the kernel function of the operator p 2 α 2 p , where J ¯ 1 ( α t ) is the modified Bessel function of the first kind.
The second path is to consider a variable substitution. Extend the parameter α to a complex and introduce the complex transformation
α i α .
Then, the fractal operator p 2 α 2 p has the following transformation:
p 2 + α 2 p p 2 α 2 p .
Thus, Equation (15) is transformed into
[ p 2 α 2 p ] δ ( t ) = i α t J 1 ( i α t ) .
The transformation between the Bessel function J m ( α t ) and the modified Bessel function J ¯ m ( α t ) is that
J m ( i α t ) = i m J ¯ m ( α t ) .
The combination of Equations (22) and (23) immediately obtains Equation (26).
[ p 2 α 2 p ] δ ( t ) = α t J ¯ 1 ( α t ) .
Similarly, Equation (16) is transformed into Equation (27).
[ p 2 α 2 p p ] δ ( t ) = 0 t i α τ J 1 ( i α τ ) d τ = i α 2 { 2 [ i α t J 0 ( i α t ) J 1 ( i α t ) ] + π i α t [ J 1 ( i α t ) H 0 ( i α t ) J 0 ( i α t ) H 1 ( i α t ) ] } .
With the help of Equation (25), Equation (27) can be rewritten as Equation (28).
[ p 2 α 2 p p ] δ ( t ) = 0 t α τ J ¯ 1 ( α τ ) d τ = α 2 { 2 [ α t J ¯ 0 ( α t ) J ¯ 1 ( α t ) ] + π α t [ J ¯ 1 ( α t ) i H 0 ( i α t ) J ¯ 0 ( α t ) H 1 ( i α t ) ] } .
On the surface, Equation (27) is strange: the first equation exists in the field of real numbers, but the imaginary number appears in the last term of the second equation. However, further analysis reveals that this imaginary number can be eliminated. Based on the transformation between modified Struve function H ¯ m ( t ) and Struve function H m ( t ) :
H m ( i t ) = i m + 1 H ¯ m ( t ) ,
H ¯ m ( t ) ( t 2 ) m + 1 n = 0 ( t 2 ) 2 n Γ ( n + 3 2 ) Γ ( n + m + 3 2 ) ,
we can obtain
H 0 ( i α t ) = i H ¯ 0 ( α t ) ,
H 1 ( i α t ) = H ¯ 1 ( α t ) .
Then Equation (28) can be changed into
[ p 2 α 2 p p ] δ ( t ) = 0 t α τ J ¯ 1 ( α τ ) d τ = α 2 { 2 [ α t J ¯ 0 ( α t ) J ¯ 1 ( α t ) ] + π α t [ J ¯ 0 ( α t ) H ¯ 1 ( α t ) J ¯ 1 ( α t ) H ¯ 0 ( α t ) ] } .
That is, the kernel function of the fractal operator p 2 α 2 p p can be characterized by modified Bessel function and modified Struve function. Comparing Equations (33) and (17), fantastic symmetry can be seen. Discard the operator term and only retain the second equation of Equation (33), we obtain
0 t 1 τ J ¯ 1 ( α τ ) d τ = 1 2 { 2 [ α t J ¯ 0 ( α t ) J ¯ 1 ( α t ) ] + π α t [ J ¯ 0 ( α t ) H ¯ 1 ( α t ) J ¯ 1 ( α t ) H ¯ 0 ( α t ) ] }
The modified Bessel function J ¯ m ( α t ) exhibits a profound intrinsic correlation with the modified Struve function H ¯ m ( α t ) . Equations (34) and (20) still exhibit an elegant symmetry, which can be further enhanced as follows. The physical parameter α is present in both Equations (20) and (34). As mentioned earlier, α is the characteristic frequency of a functional fractal. However, mathematical theorems are typically independent of physics. To reveal the independence, the following variable substitution is introduced:
α t t , α τ τ .
Then, Equations (20) and (34) are transformed into the following forms, respectively:
0 t 1 τ J 1 ( τ ) d τ = 1 2 { 2 [ t J 0 ( t ) J 1 ( t ) ] + π t [ J 1 ( t ) H 0 ( t ) J 0 ( t ) H 1 ( t ) ] } .
0 t 1 τ J ¯ 1 ( τ ) d τ = 1 2 { 2 [ t J ¯ 0 ( t ) J ¯ 1 ( t ) ] + π t [ J ¯ 0 ( t ) H ¯ 1 ( t ) J ¯ 1 ( t ) H ¯ 0 ( t ) ] } .
The physical parameters disappear, and the symmetry is enhanced.
Comparing the content of this section with that of the previous section, it becomes evident that the entire theoretical system exhibits remarkable symmetry. The Bessel function J m ( α t ) and the modified Bessel function J ¯ m ( α t ) demonstrate symmetrical characteristics. Similarly, the Struve function H m ( α t ) and the modified Struve function H ¯ m ( α t ) also possess inherent symmetry. The symmetry arises from the symmetry between operators, that is, the structure of operator p 2 + α 2 and operator p 2 α 2 are completely symmetric.
The graphs of modified Bessel function J ¯ 1 ( t ) , weighted modified Bessel function J ¯ 1 ( t ) t , and its integral 0 t 1 τ J ¯ 1 ( τ ) d τ are shown in Figure 5. When time tends to zero, we obtain Equation (37).
lim t 0 J ¯ 1 ( t ) = 0 , lim t 0 J ¯ 1 ( t ) t = 1 2 , lim t 0 0 t 1 τ J ¯ 1 ( τ ) d τ = 0
That is, the values of the three types of special functions J ¯ 1 ( t ) , J ¯ 1 ( t ) t , and 0 t 1 τ J ¯ 1 ( τ ) d τ remain finite at t = 0 . When time tends to infinity, we obtain Equation (38).
lim t J ¯ 1 ( t ) , lim t J ¯ 1 ( t ) t , lim t 0 t 1 τ J ¯ 1 ( τ ) d τ .
That is, the three types of special functions all diverge, in which 0 t 1 τ J ¯ 1 ( τ ) d τ demonstrates the fastest divergence rate, followed by J ¯ 1 ( t ) and then lim t J ¯ 1 ( t ) t . It is evident that the divergence observed in the three types of functions stems from the inherent properties of the modified Bessel function.
Since the three kinds of special functions all diverge, the fractal operators [ p 2 α 2 p ] and [ p 2 α 2 p p ] and their kernel functions do not converge.
The theoretical and practical value of operators and functions that fail to converge is significantly constrained. How can we overcome these limitations? An effective method is as follows: mathematically, decay terms can be introduced into the kernel function; physically, the parameters of the physical components can be changed, introducing a translation into the operator. See the Section 5.3 for details.

5.3. Kernel Functions of the Fractal Operators ( p + β ) 2 α 2 , 1 p ( p + β ) 2 α 2 and 1 p ( p + β ) 2 α 2

In this section, the operators ( p + β ) 2 α 2 and 1 p ( p + β ) 2 α 2 are considered. If we set a = 0 , the two operators are reduced to the two operators p 2 α 2 and 1 p p 2 α 2 in the previous section.
The direct acquisition of the kernel function of operators ( p + β ) 2 α 2 and 1 p ( p + β ) 2 α 2 poses a challenge. However, we can take a detour. By combining operator p 2 α 2 from the previous section with the translation transformation, operator ( p + β ) 2 α 2 and its kernel function can be obtained.
By translation transformation, the operator p in Equation (23) is transformed into p + a:
p p + β ,
Then the operator p 2 α 2 p is transformed into:
p 2 α 2 p ( p + β ) 2 α 2 ( p + β ) .
Combined with the translation property of the Laplace transformation, we obtain
1 [ ( p + β ) 2 α 2 ( p + β ) ] = α t J ¯ 1 ( α t ) e β t .
Then the kernel function of the fractal operator ( p + β ) 2 α 2 ( p + β ) is
[ ( p + β ) 2 α 2 ( p + β ) ] δ ( t ) = α t J ¯ 1 ( α t ) e β t .
From Equation (40), the kernel function of operator ( p + β ) 2 α 2 could be derived. Since ( p + β ) 2 α 2 is an irrational radical operator, its kernel is a non-elementary function. Notice that the decay term appears in the kernel function, specifically as the negative power exponential function e β t . If we let β 0 by translation transform, Equation (40) could be reduced to Equation (26).
When t 0 , the limit of the kernel function 1 t J ¯ 1 ( α t ) e β t is
lim t 0 [ 1 t J ¯ 1 ( α t ) e β t ] = α 2 .
When t , the limit of the kernel function 1 t J ¯ 1 ( α t ) e β t should be discussed in two cases. In the case of 0 < β α < 1 , there is
lim t [ 1 t J ¯ 1 ( α t ) e β t ] = .
While in the case of β α 1 , there is
lim t [ 1 t J ¯ 1 ( α t ) e β t ] = 0 .
Obviously, β α = 1 is exactly the turning point of the monotonicity, the detailed analysis as follows. Suppose α is fixed as α = 1 , and the graphs of the kernel function 1 t J ¯ 1 ( t ) e β t with different values of a are shown in Figure 6. As shown in Figure 6, the convergence of the operator ( p + β ) 2 α 2 ( p + β ) when t is determined by the value of β α . In the case of 0 < β α < 1 , the kernel function exhibits a decrement followed by an increment, ultimately tending towards infinity. In the case of β α = 1 , the function exactly loses its capability to increase. Evidently, in the case of β α 1 , the convergence of the kernel function 1 t J ¯ 1 ( α t ) e β t in this section is much better than that of the non-convergent kernel function 1 t J ¯ 1 ( α t ) in the previous section. It can be seen that in this case, the modulation effect of the negative power exponential function e β t is highly significant, facilitating the transformation of a non-convergent kernel function 1 t J ¯ 1 ( α t ) into a rapidly convergent one 1 t J ¯ 1 ( α t ) e β t .
For the operator ( p + β ) 2 α 2 ( p + β ) p , Equation (38) combined with the inverse Laplace transform obtains
1 [ ( p + β ) 2 α 2 ( p + β ) p ] = 0 t α τ J ¯ 1 ( α τ ) e β τ d τ .
Therefore, the kernel function of operator ( p + β ) 2 α 2 ( p + β ) p could be given:
[ ( p + β ) 2 α 2 ( p + β ) p ] δ ( t ) = α 0 t 1 τ J ¯ 1 ( α τ ) e β τ d τ .
From Equation (45), the kernel function of the operator could be 1 p ( p + β ) 2 α 2 derived. Since 1 p ( p + β ) 2 α 2 is an irrational operator, its kernel function is a non-elementary function, given by the integral of the weighted modified Bessel function. If we set β 0 by translation transform, Equation (45) could be reduced to Equation (33).
The convergence of the integral 0 t 1 τ J ¯ 1 ( α τ ) e β τ d τ is ensured due to the boundedness of the function 1 τ J ¯ 1 ( α τ ) e β τ within the interval [ 0 , t ] for any t > 0 . For the convenience of discussion, we define
0 t 1 τ J ¯ 1 ( α τ ) e β τ d τ J ¯ 1 e ( β , α t ) .
The modified Bessel function J ¯ 1 ( α τ ) is a special function. After applying weighting, 1 τ J ¯ 1 ( α τ ) e β τ is still a special function; therefore, it is reasonable to anticipate that the integral J ¯ 1 e ( β , α t ) should be a new special function.
At last, consider the operator 1 p ( p + β ) 2 α 2 . The inverse Laplace transformation gives that
1 ( 1 p ) = 1 π t .
Equations (39) and (47), combined with the convolution theorem of the Laplace transformation, obtain that
1 [ ( p + β ) 2 α 2 ( p + β ) p ] = α π 0 t 1 t τ 1 τ J ¯ 1 ( α τ ) e β τ d τ .
The kernel function of the operator in Equation (48) is given as
[ ( p + β ) 2 α 2 ( p + β ) p ] δ ( t ) = α π 0 t 1 t τ 1 τ J ¯ 1 ( α τ ) e β τ d τ .
The proof of the convergence of the integral 0 t 1 t τ 1 τ J ¯ 1 ( α τ ) e β τ d τ could be seen in Appendix B. In Equation (49), the kernel function of operator 1 p is given by Equation (47), and the property of operator p could be referred to in Courant’s work [31]. Therefore, the calculability of the operator 1 p ( p + β ) 2 α 2 is ensured.
If we let β 0 by translation transform, Equation (49) could be reduced to
[ p 2 α 2 p p ] δ ( t ) = α π 0 t 1 t τ 1 τ J ¯ 1 ( α τ ) d τ .
For the convenience of discussion, we define
0 t 1 t τ 1 τ J ¯ 1 ( α τ ) e β τ d τ J ¯ ˜ 1 e ( β , α t ) .
After applying weighting on the Bessel function J ¯ 1 ( α τ ) , 1 t τ 1 τ J ¯ 1 ( α τ ) e β τ is still a special function. It is reasonable to anticipate that the integral J ¯ ˜ 1 e ( β , α t ) should be a new special function.
In the case of 0 < β α < 1 , due to the divergence of the function 1 t J ¯ 1 ( t ) e β t , the values of the integrals J ¯ 1 e ( β , α t ) and J ¯ ˜ 1 e ( β , α t ) tends to infinity when t , indicating poor convergence of the operator and limited practical applicability. Therefore, the following content focuses on the case of β α 1 .
Suppose that α is fixed as α = 1 , β is set as β = 1 and β = 2 , respectively, and the graphs of functions J ¯ 1 e ( β , t ) and J ¯ ˜ 1 e ( β , t ) are shown in Figure 7. The two functions exhibit different monotonicity. The function J ¯ 1 e ( β , t ) tends to a constant value and saturated evolution after monotonically increasing. When t , there is
lim t J ¯ 1 e ( 1 , t ) = 1.000 ,   lim t J ¯ 1 e ( 2 , t ) = 0.268 .
The function J ¯ ˜ 1 e ( β , t ) first increases for a short time and then slowly decreases, tending to a constant value. The function J ¯ ˜ 1 e ( 1 , t ) attains its maximum value 0.568 at amount t = 1.1 , and the function J ¯ ˜ 1 e ( 2 , t ) attains its maximum value 0.387 at amount t = 0.45 .
Following the aforementioned analysis, we can obtain that all the kernel functions of the whole operator family in Table 3 exist. It can also be observed that certain operators’ kernel functions cannot be directly expressed through established special functions, necessitating the introduction of novel special functions in such cases.
In disciplines such as rock mechanics, the pressure difference often exhibits a gradual decay over time. While exponential and other elementary functions can illustrate the convergence of physical quantities, the convergence is frequently too rapid to accurately represent this slow decay. The gradual decline depicted in Figure 7 suggests that the operator 1 p ( p + β ) 2 α 2 can effectively characterize this phenomenon.

5.4. Kernel Functions of the Fractal Operators ( p + β ) 2 + α 2 , 1 p ( p + β ) 2 + α 2 and 1 p ( p + β ) 2 + α 2

Section 5.1 is symmetric with Section 5.2, and the symmetry can be used for reference. Now that we have Section 5.3, can we construct content that is symmetric with it? The answer is yes. We can construct Section 5.4, the contents of this section.
Extend the parameter α to complex field and introduce the complex transformation
α i α .
Then, the fractal operator ( p + β ) 2 α 2 ( p + β ) has the following transformation
( p + β ) 2 α 2 ( p + β ) ( p + β ) 2 + α 2 ( p + β ) .
Based on the above parameter transformation and operator transformation, the content of Section 5.3 can be transformed into the content of this section.
Equation (40) is transformed into
[ ( p + β ) 2 + α 2 ( p + β ) ] δ ( t ) = i α t J ¯ 1 ( i α t ) e β t .
From the transformation between the Bessel function and the modified Bessel function (Equation (25)), we can immediately obtain that
[ ( p + β ) 2 + α 2 ( p + β ) ] δ ( t ) = α t J 1 ( α t ) e β t .
From Equation (54), the kernel function of the operator ( p + β ) 2 + α 2 could be derived. The comparison between Equation (54) and Equation (40) reveals their complete symmetry, as both the algebraic structure of the operators and the kernel functions exhibit symmetrical characteristics.
When t 0 , there is
lim t 0 [ 1 t J 1 ( α t ) e β t ] = α 2 .
When t , the limit of the kernel function is
lim t [ 1 t J 1 ( α t ) e β t ] = 0 .
Unlike in Section 5.3, the limit of the kernel function 1 t J 1 ( α t ) e β t is not affected by the value of β α , as shown in Figure 8. This advantage stems from the inherent oscillation and convergence properties of the Bessel function.
In the case of β α 1 , both functions 1 t J 1 ( α t ) e β t and 1 t J ¯ 1 ( α t ) e β t converge rapidly to 0, indicating a high degree of similarity in their graphs. This similarity highlights the inherent symmetry in the fractal operators ( p + β ) 2 + α 2 ( p + β ) and ( p + β ) 2 + α 2 ( p + β ) .
Similarly, Equations (45) and (49) are transformed into
[ ( p + β ) 2 + α 2 ( p + β ) p ] δ ( t ) = α 0 t 1 τ J 1 ( α τ ) e β τ d τ ,
[ ( p + β ) 2 + α 2 ( p + β ) p ] δ ( t ) = α π 0 t 1 t τ 1 τ J 1 ( α τ ) e β τ d τ .
The proof of the convergence of the integrals at the right side of Equations (57) and (58) is identical to that in Section 5.3. Thus, the calculability of the operators ( p + β ) 2 + α 2 p and ( p + β ) 2 + α 2 p is ensured.
If we set β 0 by translation transform, Equation (57) could be reduced to Equation (17), and Equation (58) could be reduced to
[ p 2 + α 2 p p ] δ ( t ) = α π 0 t 1 t τ 1 τ J 1 ( α τ ) d τ .
Corresponding to Section 5.3, we have the definitions
0 t 1 τ J 1 ( α τ ) e β τ d τ J 1 e ( β , α t ) ,
0 t 1 t τ 1 τ J 1 ( α τ ) e β τ d τ J ˜ 1 e ( β , α t ) .
Suppose α is fixed as α = 1 and β is set to β = 1 and β = 2 , respectively, and the graphs of functions J 1 e ( β , t ) and J ˜ 1 e ( β , t ) are shown in Figure 9.
By comparing Figure 9 and Figure 7, it can be observed that the graph trends of J 1 e ( β , t ) and J ¯ 1 e ( β , t ) are similar, the similarity increasing with the increase of the value of a. The similarity can be attributed to the regulation imposed by the exponential function. This shows that the fractal operators ( p + β ) 2 + α 2 ( p + β ) p and ( p + β ) 2 α 2 ( p + β ) p are symmetric, the symmetry increasing with the increase of the value of β . When t , there is
lim t J 1 e ( 1 , t ) = 414 ,   lim t J 1 e ( 2 , t ) = 0.236 .
It can be seen that in the case of β = 2 , the difference between functions J 1 e ( 2 , t ) and J ¯ 1 e ( 2 , t ) is already minimal, and their graph contours are exactly the same.
Similarly, the graph trends of J ˜ 1 e ( β , t ) and J ¯ ˜ 1 e ( β , t ) in Figure 9 and Figure 7, respectively, are also similar, which means that the fractal operators ( p + β ) 2 + α 2 ( p + β ) p and ( p + β ) 2 α 2 ( p + β ) p also exhibit a pleasing symmetry. The function J ˜ 1 e ( 1 , t ) attains its maximum value 0.522 at amount t = 0.75 , and the function J ˜ 1 e ( 2 , t ) attains its maximum value 0.379 at amount t = 0.4 . Also, we can see that in the case of β = 2 , the difference between functions J ˜ 1 e ( 2 , t ) and J ¯ ˜ 1 e ( 2 , t ) is already minimal, and their graph contours are exactly the same.
The logical basis of the above symmetries is the symmetry between the Bessel function and the modified Bessel function, as well as the symmetry between the Struve function and the modified Struve function. It is quite an unexpected result.
After comparing Figure 8, Figure 9, and Figure 4, it becomes evident that the oscillation characteristics induced by the Bessel function J 1 ( t ) are completely eradicated following the modulation of 1 t J 1 ( t ) e 2 t , leaving no discernible trace of oscillation on the curve, which highlights the efficiency of eliminating oscillations by designing parameters with the help of fractal operators.

5.5. Kernel Functions of the Fractal Operators p + β , 1 p p + β and 1 p p + β

p + β is a classical operator. Jian et al. [19] conducted a detailed analysis and derived the following expression:
p + α 2 η ( t ) = α erf ( α 2 t ) + e α 2 t π t ,
where η ( t ) is the Heaviside step function. Setting β = α 2 and noting that
p η ( t ) = δ ( t ) ,
we can give that
1 p p + β δ ( t ) = β erf ( β t ) + e β t π t ,
where e r f ( ) is the error function.
The operator 1 p p + β will now be considered from a different perspective. When β > 0 , the inverse Laplace transformation gives that
1 ( 1 p p + β ) = e β t π t + β erf ( β t ) .
Therefore, the kernel function of the fractal operator is
1 p p + β δ ( t ) = e β t π t + β erf ( β t ) .
Equation (67) is consistent with the result derived by Jian et al. [19]. Since 1 p p + β is an irrational fractional-type operator, the kernel function contains the non-elementary function erf ( ) .
Then, consider the operator 1 p p + β . For the Laplace transformation of the modified Bessel function of the first kind J ¯ 0 and J ¯ 1 , the following equation is valid:
{ β 2 [ J ¯ 0 ( β t 2 ) + J ¯ 1 ( β t 2 ) ] } = p + β 2 p β 2 1 .
Observing the right side of Equation (68), if we introduce the translation transformation, we can obtain that:
p p + β 2 ,
Then, we have
p + β 2 p β 2 p + β p .
That is, when the operator p is translated to p + β 2 , the operator is exactly translated to p + β p . Therefore, Equation (68) combined with the translation property of the Laplace transformation obtains that
{ β 2 e β t 2 [ J ¯ 0 ( β t 2 ) + J ¯ 1 ( β t 2 ) ] } = p + β p 1 .
It can be derived further that
[ 1 p p + β 1 ] δ ( t ) = β 2 e β t 2 [ J ¯ 0 ( β t 2 ) + J ¯ 1 ( β t 2 ) ] .
1 p p + β 1 is an irrational fractional-type operator, so its kernel function contains the non-elementary function, specifically the modified Bessel function.

6. Conclusions

In this paper, the functional fractal circuit model of arterial blood flow was reconstructed, the fractal operators were abstracted, and their kernel functions were studied. The main conclusions drawn from this are as follows:
  • The structure and motion of living organisms offer abundant resources for exploring new spatial forms and fractal operators. After taking the shunt effect of the branch vessels into account, the reconstructed functional circuit model exhibited a distinct topological structure compared to that of previous models. Such topological structure can indeed yield insightful results.
  • By selecting basic admittance operators, the fractal operator family abstracted from the fractal circuit model constituted the fundamental graph of the apparent 1/2-order operators. The kernel functions and properties of these operators offer novel perspectives on the fields of automatic control and biomechanics.
  • The Bessel function and Struve function were related through explicit integration; the modified Bessel function and modified Struve function also exhibited analogous connection. The connections could be employed for constructing the kernel functions of operators 1 p p 2 + α 2 and 1 p p 2 α 2 .
  • By translation transformation, the divergent operators p 2 α 2 and 1 p p 2 α 2 were transformed into convergent operators ( p + β ) 2 α 2 and 1 p ( p + β ) 2 α 2 , greatly broadening the scope of the operators.
  • Under certain conditions, operators ( p + β ) 2 α 2 , 1 p ( p + β ) 2 α 2 , and 1 p ( p + β ) 2 α 2 demonstrate notable symmetry with operators ( p + β ) 2 + α 2 , 1 p ( p + β ) 2 + α 2 , and 1 p ( p + β ) 2 + α 2 respectively. The constructions of the kernel functions for those operators necessitates the introduction of novel special functions, which possess properties challenging to articulate by elementary functions.
  • These results provide a powerful mathematical lexicon and toolkit for the biomechanics.
The results presented in this paper possess broad applicability and can be effectively employed in the fields of biomechanics and automatic control. The scope of this paper is limited to the combination of the fundamental electrical admittance operators: inductance, capacitance, and resistance. To enhance our understanding of biological phenomena, future research will incorporate more intricate electrical components (including fractional-order components) for a more comprehensive depiction.

Author Contributions

Conceptualization, T.Z., Y.Y., G.P., C.L. and Z.J.; methodology, T.Z., Y.Y., G.P., C.L. and Z.J.; software, T.Z. and Y.Y.; validation, T.Z., Y.Y. and C.L.; formal analysis, T.Z., Y.Y. and G.P.; investigation, T.Z., Y.Y. and G.P.; writing—original draft preparation, T.Z. and Y.Y.; writing—review and editing, T.Z. and Y.Y.; funding acquisition, Y.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the National Natural Science Foundation of China, grant number 12050001.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. Proof of Equation (19)

The following equation without parameter α is considered first:
{ 1 2 { 2 [ t J 0 ( t ) J 1 ( t ) ] + π t [ J 1 ( t ) H 0 ( t ) J 0 ( t ) H 1 ( t ) ] } } = 1 t J 1 ( t ) .
The Bessel function exhibits recursion properties:
J n ( t ) = 1 2 J n 1 ( t ) 1 2 J n + 1 ( t ) ,
n J n ( t ) = t 2 J n 1 ( t ) + t 2 J n + 1 ( t ) .
From Equations (A2) and (A3) we can derive
J n ( t ) = n t J n ( t ) J n + 1 ( t ) .
The Struve function exhibits recursion properties:
H n 1 ( t ) H n + 1 ( t ) = 2 H n ( t ) 1 π Γ ( n + 3 2 ) ( t 2 ) n ,
H n 1 ( t ) + H n + 1 ( t ) = 2 n t H n ( t ) + 1 π Γ ( n + 3 2 ) ( t 2 ) n .
From Equations (A2), (A4) and (A5), we can derive that
J 0 ( t ) = J 1 ( t ) ,
J 1 ( t ) = 1 2 J 0 ( t ) 1 2 J 2 ( t ) ,
H 0 ( t ) = 1 2 [ H 1 ( t ) + H 1 ( t ) ] + 4 π ,
H 1 ( t ) = 1 2 [ H 0 ( t ) + H 2 ( t ) ] + t 3 π .
Substituting Equations (A7)–(A10) into the left side of Equation (A1) obtains
{ 1 2 { J 1 ( t ) [ 2 + π t H 0 ( t ) ] + t J 0 ( t ) [ 2 π H 1 ( t ) ] } } = 1 2 J 0 ( t ) 1 6 t 2 J 0 ( t ) 1 2 t J 1 ( t ) + 1 2 J 2 ( t ) + 1 4 π t J 1 ( t ) H 1 ( t ) + 1 2 π J 1 ( t ) H 0 ( t ) 1 4 π t J 2 ( t ) H 0 ( t ) 1 2 π J 0 ( t ) H 1 ( t ) + 1 4 π t J 1 ( t ) H 1 ( t ) + 1 4 π t J 0 ( t ) H 2 ( t ) .
From Equation (A3), we can obtain that the right side of Equation (A1) satisfies
1 t J 1 ( t ) = 1 2 [ J 0 ( t ) + J 2 ( t ) ] .
From Equations (A11) and (A12), we can obtain that Equation (A1) is equivalent to the following equation:
1 4 π J 1 ( t ) H 1 ( t ) + 1 2 t π J 1 ( t ) H 0 ( t ) 1 4 π J 2 ( t ) H 0 ( t ) 1 2 t π J 0 ( t ) H 1 ( t ) + 1 4 π J 1 ( t ) H 1 ( t ) + 1 4 π J 0 ( t ) H 2 ( t ) = t 6 J 0 ( t ) + 1 2 J 1 ( t ) .
Then, we only need to prove Equation (A13). Substituting Equation (A12) into the left side of the Equation (A13); eliminating 1 t J 1 ( t ) , (A13) changes into
1 4 π J 1 ( t ) H 1 ( t ) + 1 4 π J 0 ( t ) H 0 ( t ) 1 2 t π J 0 ( t ) H 1 ( t ) + 1 4 π J 1 ( t ) H 1 ( t ) + 1 4 π J 0 ( t ) H 2 ( t ) = t 6 J 0 ( t ) + 1 2 J 1 ( t ) .
From Equation (A6) we can obtain that
H 0 ( t ) + H 2 ( t ) = 2 t H 1 ( t ) + 2 t 3 π ,
H 1 ( t ) + H 1 ( t ) = 2 π .
Substituting Equation (A15) into Equation (A14), eliminating H 0 ( t ) + H 2 ( t ) , Equation (A14) changes into
1 4 π J 1 ( t ) H 1 ( t ) + 1 4 π J 1 ( t ) H 1 ( t ) + 1 6 t J 0 ( t ) = t 6 J 0 ( t ) + 1 2 J 1 ( t ) .
Substituting Equation (A16) into the left side of Equation (A17) directly obtains the right side of Equation (A17); therefore, Equation (A17) is valid, i.e., the Equation (A1) is valid.
In Equation (A1), substituting the variable t into α t , i.e.,
t α t ,
we can obtain that
d d ( α t ) { 1 2 { 2 [ α t J 0 ( α t ) J 1 ( α t ) ] + π α t [ J 1 ( α t ) H 0 ( α t ) J 0 ( α t ) H 1 ( α t ) ] } } = 1 α t J 1 ( α t ) ,
i.e.,
d d t { 1 2 { 2 [ α t J 0 ( α t ) J 1 ( α t ) ] + π α t [ J 1 ( α t ) H 0 ( α t ) J 0 ( α t ) H 1 ( α t ) ] } } = 1 t J 1 ( α t ) .
The Equation (A19) is obviously equivalent to Equation (19); therefore, Equation (19) is proved.

Appendix B. Proof of the Convergence of the Integral 0 t 1 τ t τ J ¯ 1 ( α τ ) e β τ d τ

The theorem of convergence for improper integrals is introduced:
Theorem A1. 
Suppose that the function f ( x ) is continuous over the interval ( a , b ] , and there is f ( x ) 0 , and x = a is a flaw of f ( x ) . If there exists 0 < q < 1 , which makes the limit lim x a + ( x a ) q f ( x ) exist; then, the improper integral a b f ( x ) d x is convergent.
In the integral 0 t 1 τ t τ J ¯ 1 ( α τ ) e β τ d τ , introduce the variable substitution ζ = τ ; then, we have
0 t 1 t τ 1 τ J ¯ 1 ( α τ ) e β τ d τ = 0 t 1 t + ( τ ) 1 τ J ¯ 1 ( α τ ) e β ( τ ) d τ = 0 t 1 t + ζ 1 ( ζ ) J ¯ 1 ( α ζ ) e β ζ d ( ζ ) = t 0 1 t + ζ 1 ζ J ¯ 1 ( α ζ ) e β ζ d ζ .
We have α , β > 0 , and t > 0 , so there is 1 ζ t + ζ J ¯ 1 ( α ζ ) e β ζ 0 over the interval t < ζ < 0 . Obviously, ζ = t is the flaw of the integrand. At ζ = 0 , the limit of the integrand is lim ζ 0 1 t + ζ 1 ζ J ¯ 1 ( α ζ ) e β ζ = α 2 t , i.e., the limit exists. Therefore, the case at ζ = 0 does not affect the convergence of the entire integral. Based on the fact that
lim ζ ( t ) + ( ζ + t ) 1 2 1 t + ζ 1 ζ J ¯ 1 ( α ζ ) e β ζ = lim ζ ( t ) + 1 ζ J ¯ 1 ( α ζ ) e β ζ = 1 t J ¯ 1 ( α t ) e β t ,
the limit lim ζ ( t ) + ( ζ + t ) 1 2 1 t + ζ 1 ζ J ¯ 1 ( α ζ ) e β ζ exists. Then according to the Theorem, the integral t 0 1 t + ζ 1 ζ J ¯ 1 ( α ζ ) e β ζ d ζ is convergent, i.e., the integral 0 t 1 τ t τ J ¯ 1 ( α τ ) e β τ d τ is convergent.

References

  1. Hales, S. Statical Essays: Containing Haemostaticks; Innys and Manby: London, UK, 1733. [Google Scholar]
  2. Frank, O. Die grundform des arteriellen pulses. Ztg. Biol. 1899, 37, 483–586. [Google Scholar]
  3. Goldwyn, R.M.; Watt, T.B. Arterial pressure pulse contour analysis via a mathematical model for the clinical quantification of human vascular properties. IEEE Trans. Biomed. Eng. 1967, BME-14, 11–17. [Google Scholar] [CrossRef]
  4. Burattini, R.; Gnudi, G. Computer identification of models for the arterial tree input impedance: Comparison between two new simple models and first experimental results. Med. Biol. Eng. Comput. 1982, 20, 134–144. [Google Scholar] [CrossRef] [PubMed]
  5. Stergiopulos, N.; Westerhof, B.E.; Westerhof, N. Total arterial inertance as the fourth element of the windkessel model. Am. J. Physiol.-Heart C 1999, 276, H81–H88. [Google Scholar] [CrossRef] [PubMed]
  6. Jager, G.N.; Westerhof, N.; Noordergraaf, A. Oscillatory flow impedance in electrical analog of arterial system: Representation of sleeve effect and non-Newtonian properties of blood. Circ. Res. 1965, 16, 121–133. [Google Scholar] [CrossRef] [PubMed]
  7. Westerhof, N.; Elzinga, G.; Sipkema, P. An artificial arterial system for pumping hearts. J. Appl. Physiol. 1971, 31, 776–781. [Google Scholar] [CrossRef] [PubMed]
  8. Shim, Y.T.; Pasipoularides, A.; Straley, C.A.; Hampton, T.G.; Soto, P.F.; Owen, C.H.; Davis, J.W.; Glower, D.D. Arterial windkessel parameter estimation: A new time-domain method. Ann. Biomed. Eng. 1994, 22, 66–77. [Google Scholar] [CrossRef] [PubMed]
  9. Mandeville, J.B.; Marota, J.J.A.; Ayata, C.; Zaharchuk, G.; Moskowitz, M.A.; Rosen, B.R.; Weisskoff, R.M. Evidence of a cerebrovascular postarteriole Windkessel with delayed compliance. J. Cereb. Blood Flow Metab. 1999, 19, 679–689. [Google Scholar] [CrossRef]
  10. Rupnic, M.; Runvovc, F. Simulation of steady state and transient phenomena by using the equivalent electronic circuit. Comput. Meth. Programs Biomed. 2002, 67, 1–12. [Google Scholar] [CrossRef]
  11. Abdolrazaghi, M.; Navidbakhsh, M.; Hassani, K. Mathematical modelling of intra-aortic balloon pump. Comput. Methods Biomech. Biomed. Eng. 2010, 13, 567–576. [Google Scholar] [CrossRef]
  12. Gul, R.; Schütte, C.; Bernhard, S. Mathematical modeling and sensitivity analysis of arterial anastomosis in the arm. Appl. Math. Model. 2016, 40, 7724–7738. [Google Scholar] [CrossRef]
  13. Baker, W.B.; Parthasarathy, A.B.; Gannon, K.P.; Kavuri, V.C.; Busch, D.R.; Abramson, K.; He, L.; Mesquita, R.C.; Mullen, M.T.; Detre, J.A.; et al. Noninvasive optical monitoring of critical closing pressure and arteriole compliance in human subjects. J. Cereb. Blood Flow Metab. 2017, 37, 2691–2705. [Google Scholar] [CrossRef]
  14. Li, B.; Mao, B.Y.; Feng, Y.; Liu, J.C.; Zhao, Z.; Duan, M.Y.; Liu, Y.J. The hemodynamic mechanism of FFR-guided coronary arterybypass grafting. Front. Physiol. 2021, 12, 8. [Google Scholar]
  15. Nichols, W.W.; Conti, C.R.; Walker, W.E.; Milnor, W.R. Input impedance of the systemic circulation in man. Circ. Res. 1971, 40, 451–458. [Google Scholar] [CrossRef]
  16. Packard, R.R.S.; Luo, Y.; Abiri, P.; Jen, N.; Aksoy, O.; Suh, W.M.; Tai, Y.C.; Hsiai, T.K. 3-D Electrochemical Impedance Spectroscopy Mapping of Arteries to Detect Metabolically Active but Angiographically Invisible Atherosclerotic Lesions. Theranostics 2017, 7, 2431–2442. [Google Scholar] [CrossRef] [PubMed]
  17. Guo, J.Q.; Yin, Y.J.; Ren, G.X. Abstraction and operator characterization of fractal ladder viscoelastic hyper-cell for ligaments and tendons. Appl. Math. Mech. 2019, 40, 1429–1448. [Google Scholar] [CrossRef]
  18. Guo, J.Q.; Yin, Y.J.; Hu, X.L.; Ren, G.X. Self-similar network model for fractional-order neuronal spiking: Implications of dendritic spine functions. Nonlinear Dynam 2020, 100, 921–935. [Google Scholar] [CrossRef]
  19. Jian, Z.M.; Peng, G.; Li, D.A.; Yu, X.B.; Yin, Y.J. Correlation between Convolution Kernel Function and Error Function of Bone Fractal Operators. Fractal Fract. 2023, 7, 707. [Google Scholar] [CrossRef]
  20. Peng, G.; Guo, J.Q.; Yin, Y.J. Self-Similar Functional Circuit Models of Arteries and Deterministic Fractal Operators: Theoretical Revelation for Biomimetic Materials. Int. J. Mol. Sci. 2021, 22, 12897. [Google Scholar] [CrossRef]
  21. Yin, Y.J.; Guo, J.Q.; Peng, G.; Yu, X.B.; Kong, Y.Y. Fractal Operators and Fractional Dynamics with 1/2 Order in Biological Systems. Fractal Fract. 2022, 6, 378. [Google Scholar] [CrossRef]
  22. Yin, Y.J.; Peng, G.; Yu, X.B. Algebraic equations and non-integer orders of fractal operators abstracted from biomechanics. Acta Mech. Sin. 2022, 38, 521488. [Google Scholar]
  23. Colombaro, I.; Garra, R.; Giusti, A.; Mainardi, F. Scott-Blair models with time-varying viscosity. Appl. Math. Lett. 2018, 86, 57–63. [Google Scholar] [CrossRef]
  24. Zhou, H.W.; Wang, C.P.; Han, B.B.; Duan, Z.Q. A creep constitutive model for salt rock based on fractional derivatives. Int. J. Rock Mech. Min. Sci. 2011, 48, 116–121. [Google Scholar] [CrossRef]
  25. Yin, D.S.; Zhang, W.; Cheng, C.; Li, Y.Q. Fractional time-dependent Bingham model for muddy clay. J. Non-Newton. Fluid Mech. 2012, 187, 32–35. [Google Scholar] [CrossRef]
  26. Westerhof, N. Snapshots of Hemodynamics an Aid for Clinical Research and Graduate Education, 3rd ed.; Springer International Publishing: Cham, Switzerland, 2019. [Google Scholar]
  27. Mandelbrot, B.B.; Aizenman, M. Fractals: Form, Chance, and Dimension. Phys. Today 1979, 32, 65–66. [Google Scholar] [CrossRef]
  28. Feder, J. Fractals; Plenum Press: New York, NY, USA, 1988. [Google Scholar]
  29. Mikusinski, J. Operational Calculus, 2nd ed.; Pergamon Press: Oxford, UK, 1983. [Google Scholar]
  30. Yu, X.B.; Yin, Y.J. Operator Kernel Functions in Operational Calculus and Applications in Fractals with Fractional Operators. Fractal Fract. 2023, 7, 755. [Google Scholar] [CrossRef]
  31. Heaviside, O. On Operators in Physical Mathematics, Part I. Proc. R. Soc. Lond. Ser. A-Math. Phys. Eng. Sci. 1893, 52, 315–320. [Google Scholar]
  32. Courant, R.; Hilbert, D. Methods of Mathematical Physics; Interscience Publishers: New York, NY, USA, 1962. [Google Scholar]
  33. Milnor, W.R. Hemodynamics, 2nd ed.; William & Wilkins: Baltimore, MD, USA, 1989. [Google Scholar]
  34. Screen, H.R.C. Investigating load relaxation mechanics in tendon. J. Mech. Behav. Biomed. Mater. 2008, 1, 51–58. [Google Scholar] [CrossRef]
  35. Shen, Z.L.; Kahn, H.; Ballarini, R.; Eppell, S.J. Viscoelastic properties of isolated collagen fibrils. Biophys. J. 2011, 100, 3008–3015. [Google Scholar] [CrossRef]
  36. Joglekar, M.R.; Mejias, J.F.; Yang, G.R.; Wang, X.J. Inter-areal balanced amplification enhances signal propagation in a large-scale circuit model of the primate cortex. Neuron 2018, 98, 222–234.e8. [Google Scholar] [CrossRef]
  37. Chaudhuri, R.; Knoblauch, K.; Gariel, M.A.; Kennedy, H.; Wang, X.J. A large-scale circuit mechanism for hierarchical dynamical processing in the primate cortex. Neuron 2015, 88, 419–431. [Google Scholar] [CrossRef] [PubMed]
  38. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables; Dover Publications: New York, NY, USA, 1972. [Google Scholar]
  39. Dai, Z.; Peng, Y.; Mansy, H.A.; Sandler, R.H.; Royston, T.J. A model of lung parenchyma stress relaxation using fractional viscoelasticity. Med. Eng. Phys. 2015, 37, 752–758. [Google Scholar] [CrossRef] [PubMed]
  40. Craiem, D.; Rojo, F.J.; Atienza, J.M.; Armentano, R.L.; Guinea, G.V. Fractional-order viscoelasticity applied to describe uniaxial stress relaxation of human arteries. Phys. Med. Biol. 2008, 53, 4543. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The basic pattern of arterial flow: (a) The arterial wall contracts. (b) The arterial wall dilates. There are three distinct patterns of flow: axial (longitudinal) flow, radial (transverse) flow, and the lateral branch flow, indicated by the red arrows. When the arterial wall contracts and dilates, the directions of the transverse flow are opposite. The direction of the longitudinal flow is constant. The branch flow exists when the arterial wall contracts. More details are explained in the content.
Figure 1. The basic pattern of arterial flow: (a) The arterial wall contracts. (b) The arterial wall dilates. There are three distinct patterns of flow: axial (longitudinal) flow, radial (transverse) flow, and the lateral branch flow, indicated by the red arrows. When the arterial wall contracts and dilates, the directions of the transverse flow are opposite. The direction of the longitudinal flow is constant. The branch flow exists when the arterial wall contracts. More details are explained in the content.
Fractalfract 08 00420 g001
Figure 2. Functional fractal circuit and fractal cell. In the functional fractal circuit, electronic components M 1 and M 2 denote longitudinal and transverse flow in the main arterial cavity respectively, while M 3 denotes the flow in branch vessels. The right side gives the smallest repeating unit of the circuit, called the fractal cell.
Figure 2. Functional fractal circuit and fractal cell. In the functional fractal circuit, electronic components M 1 and M 2 denote longitudinal and transverse flow in the main arterial cavity respectively, while M 3 denotes the flow in branch vessels. The right side gives the smallest repeating unit of the circuit, called the fractal cell.
Fractalfract 08 00420 g002
Figure 3. The equivalence of the fractal cell and the fractal element.
Figure 3. The equivalence of the fractal cell and the fractal element.
Fractalfract 08 00420 g003
Figure 4. The graphs of three special classes of functions determined by the Bessel function.
Figure 4. The graphs of three special classes of functions determined by the Bessel function.
Fractalfract 08 00420 g004
Figure 5. The graphs of three special classes of functions determined by the modified Bessel function.
Figure 5. The graphs of three special classes of functions determined by the modified Bessel function.
Fractalfract 08 00420 g005
Figure 6. The graphs of the function 1 t J ¯ 1 ( t ) e β t with different values of β .
Figure 6. The graphs of the function 1 t J ¯ 1 ( t ) e β t with different values of β .
Fractalfract 08 00420 g006
Figure 7. The graphs of functions J ¯ 1 e ( β , t ) and J ¯ ˜ 1 e ( β , t ) with different values of β .
Figure 7. The graphs of functions J ¯ 1 e ( β , t ) and J ¯ ˜ 1 e ( β , t ) with different values of β .
Fractalfract 08 00420 g007
Figure 8. The graphs of the function 1 t J 1 ( t ) e β t with different values of β .
Figure 8. The graphs of the function 1 t J 1 ( t ) e β t with different values of β .
Fractalfract 08 00420 g008
Figure 9. The graphs of functions J 1 e ( β , t ) and J ˜ 1 e ( β , t ) with different values of β .
Figure 9. The graphs of functions J 1 e ( β , t ) and J ˜ 1 e ( β , t ) with different values of β .
Fractalfract 08 00420 g009
Table 1. The comparative correlation between parameters of hemodynamics and functional circuits [20,26].
Table 1. The comparative correlation between parameters of hemodynamics and functional circuits [20,26].
HemodynamicElectricQuantitative Correlation
Pressure PVoltage U P = U
Flow rate QCurrent I Q = I
Blood inertia ρ Inductance L L = c L ρ / ( π r 2 )
Wall Compliance ECapacitance C C = 3 π r 2 / ( 2 E h )
Blood Viscosity μ Resistance R R = 8 c R μ / ( π r 4 )
Table 2. Admittance operators for the three types of basic components [20].
Table 2. Admittance operators for the three types of basic components [20].
ComponentAdmittance OperatorNote
Inductance L T L = 1 L p First-order integral operator
Capacitance C T C = C p First-order derivative operator
Resistance R T R = 1 R Zero-order operator
Table 3. The family of fractal operators derived by the functional fractal circuit.
Table 3. The family of fractal operators derived by the functional fractal circuit.
T 1 , T 2 = T C , T L T 1 , T 2 = T C , T R T 1 , T 2 = T L , T R
T 3
T R 1 2 R [ 1 ± 1 + 4 R C ( p + 1 L C 1 p ) ] 1 2 R [ 1 ± 5 + 4 R C p ] 1 2 R [ 1 ± 5 + 4 R L 1 p ]
T L 1 2 L [ 1 p ± 5 p 2 + 4 L C ] 1 2 L [ 1 p ± 1 p 2 + 4 L R 1 p + 4 L C ] 1 2 L [ 1 p ± 5 p 2 + 4 L R 1 p ]
T C C 2 [ p ± 5 p 2 + 4 L C ] C 2 [ p ± 5 p 2 + 4 R C p ] C 2 [ p ± p 2 + 4 R C p + 4 L C ]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhou, T.; Yin, Y.; Peng, G.; Luo, C.; Jian, Z. Fractal Operators Abstracted from Arterial Blood Flow. Fractal Fract. 2024, 8, 420. https://doi.org/10.3390/fractalfract8070420

AMA Style

Zhou T, Yin Y, Peng G, Luo C, Jian Z. Fractal Operators Abstracted from Arterial Blood Flow. Fractal and Fractional. 2024; 8(7):420. https://doi.org/10.3390/fractalfract8070420

Chicago/Turabian Style

Zhou, Tianyi, Yajun Yin, Gang Peng, Chaoqian Luo, and Zhimo Jian. 2024. "Fractal Operators Abstracted from Arterial Blood Flow" Fractal and Fractional 8, no. 7: 420. https://doi.org/10.3390/fractalfract8070420

APA Style

Zhou, T., Yin, Y., Peng, G., Luo, C., & Jian, Z. (2024). Fractal Operators Abstracted from Arterial Blood Flow. Fractal and Fractional, 8(7), 420. https://doi.org/10.3390/fractalfract8070420

Article Metrics

Back to TopTop