Open access peer-reviewed chapter

Non-Linear Behaviours in the Dynamics of Some Biostructures

By Emil Anton, Anna Gavrilut, Maricel Agop and Daniel Timofte

Submitted: November 15th 2017Reviewed: January 25th 2018Published: July 18th 2018

DOI: 10.5772/intechopen.74435

Downloaded: 273

Abstract

Various differentiable models are frequently used to describe the dynamics of complex systems (see the kinetic models, fluid models, etc.). Given the complexity of all the physical phenomena involved in the dynamics of such systems, it is required to introduce the dynamic variable dependencies both on the space-time coordinates and on the scale resolutions. Therefore, in this case an adequate theoretical approach may be the use of non-linear physical models either in the form of the Scale Relativity Theory or of the Extended Scale Relativity Theory, i.e., the Scale Relativity Theory with an arbitrary constant fractal dimension. In the framework of the Extended Scale Relativity Theory, fractal velocity field is described both by topological solitons of kink type and by non-topological soliton varieties of breather type. Applications for the blood flow are proposed. The results revealed the directional flow toward the walls, which can explain the thickening effect which is one of the source of arteriosclerosis.

Keywords

  • complex system
  • fractal
  • non-linear
  • bio-structure
  • non-differentiable

1. Introduction

Bio-structure is a complex non-Newtonian fluid made of: plasma and formed cells, cholesterol vesicles and other suspended elements [1]; thus, the laws of fractal physics are completely applicable to sanguine circulation. For conformity, the perfect Newtonian fluid is a fluid in which viscosity is independent of the shear stress, thus having no relation to the sanguine fluid.

However, not only the complex structure of bio-structure justifies the using of fractality, but also the complex structure of the arterial system, with its multiple ramifications, which generate turbulence areas and interruptions of the linear flowing that make classical physics not applicable in this context. We can thus discuss about multi-fractality: a morphological one due to complex structure of the arterial tree as well as a functional one due to bio-structure flow “regimes” [2]. Also, the stress of a visco-elastic fluid, unlike the Newtonian fluid, depends not only on the actually stress applied but on the one applied during previous deformation of the fluid [3, 4].

Standard theoretical models usually used in complex fluid dynamics and particularly that of flow through bio-structure vessels are ambiguous [5]. The assessment that the entities contained by bio-structure move along continuous and differentiable courses proves to be false, as it cannot comprise the entire variety of dynamics that are induced by the flowing of bio-structure through the vessel system (from the separation of bio-structure components through turbulence regimes to bio-structure-bio-structure vessel interactions).

In this context, we propose a new hypothesis according to which bio-structure structural unit move on continuous but non-differentiable curves (and particularly on fractal curves). We cannot predict the entirety of the bio-structure-vessel system, bio-structure-organic tissues, etc., or at elementary level, bio-structure entity-bio-structure structural unit (i.e., lymphocyte – granulocyte, lymphocyte – platelet and others) interactions. This is why accepting the above stated hypothesis comes up as a simple, elegant and efficient solution, the impossibility of predicting all these interactions that take place being substituted by the use of fractality [6].

We are led to the dynamics of a special type of fluid, free of interactions, in which the stream lines are continuous and non-differentiable curves.

Multiple physical models have been developed in the attempt to explain the dynamics of bio-structure flow and its physiological and pathological changes on the course of the entire arterial trunk, starting from the big elastic arteries and continuing with the small arteries and arterioles – resistance vessels -, following with the bio-structure capillaries (with arterial and venous components) and, backwards, with the post capillary veins, then with the middle and large veins – capacitance vessels – ensuring the anti-gravitational mobilization of bio-structure.

Thus, the hypothesis of the geometric risk factor in the development of the circulatory system’s suffering has been proposed. This initially promising theory [7] proved its drawbacks that derive from a non-differentiable, Euclidean approach to the dynamics of bio-structure flow and its effects on the vessels’ wall. The proposed counterbalance is represented by fractal physics laws [4, 8] that offer a great amount of freedom due to accepting the relativity of the complex fluids’ behavior (bio-structure belongs to this category).

Correspondingly, the theoretical models that describe the complex fluids’ dynamics are sophisticated [4, 8]. However, these models can be simplified since the complexity of the interaction process imposes various temporal resolution scales. Also, one should take into account the fact that the pattern evolution imposes different freedom degrees [3, 28].

One could develop new theoretical models assuming the fact that the complex fluids displaying chaotic behavior are recognized to acquire self-similarity in association with strong fluctuations at all possible space-time scales [4, 9]. One can replace the deterministic trajectories by collections of potential trajectories for temporal scales which are large with respect to the inverse of the highest Lyapunov exponent (see, e.g., [10, 11]). Also, the concept of definite positions can be replaced by that of probability density. An interesting example is represented by that of collisions processes in complex system, where the dynamics of the particles can be described by fractal (non-differentiable) curves.

Since fractality (non-differentiability) is a universal property of complex fluids, one should build a fractal (non-differentiable) physics. In this context, replacing the complexity of the interactions processes by non-differentiability, it is not necessary anymore to use the classical “arsenal” of quantities from the differentiable physics. This was developed in the Scale Relativity Theory (SRT) [12, 13] and in the non-standard Scale Relativity Theory (NSRT), i.e., Scale Relativity Theory with arbitrary constant fractal dimension [14]. In the framework of SRT or NSRT let us suppose that the complex system structural unit motions take place on continuous but non-differentiable curves (called fractal curves). In this way, all physical phenomena that are involved in the dynamics depend on the space-time coordinates and also on the space-time scales resolution. In this context, the physical quantities describing the complex systems dynamics could be considered to be fractal functions [13, 14]. Additionally, the complex system entities could be reduced to and identified with their own trajectories. In this way, the complex system behaves as a special interaction-less “fluid” by means of its geodesics in a non-differentiable (fractal) space.

In the present chapter, various bio-structure flow dynamics are analyzed aiming to propose new mechanisms for the genesis and evolution of different bio-structure-related pathologies (arterial occlusion, cholesterol deposition, etc.).

2. Material and method

Assuming that the motions of bio-structure’s structural units take place on continuous but non-differentiable curves (fractal curves), the following consequences emerge [13, 14]:

  1. (i) Any continuous but non-differentiable curve of bio-structure’s structural units (bio-structure fractal curve) is scale resolution δtdependent. This means that when δttends to zero, its length tends to infinity.

Let us recall that a curve is non-differentiable if it satisfies the Lebesgue theorem [9]. This means that when the scale resolution goes to zero, its length becomes infinite. In consequence, in this limit, a curve has a zigzag form and consequently it has the property of self-similarity in every one of its points. Since every part reflects the whole, this can be translated into a holography property [9];

  1. (ii) The physics of bio-structure phenomena are related to the behavior of a set of functions during the zoom operation of the resolution scale δt. Through the substitution principle, δtcan be identified with dt, that is, δtdt. Consequently, it will play the role of an independent variable. We shall use the notation dtfor the usual time as in the Hamiltonian bio-structure dynamics;

  2. (iii) The dynamics of bio-structure’s structural units are described by means of fractal variables. Since the differential time reflection invariance of any dynamical variable is broken, these functions depend on the space-time coordinates and also on the resolution scale. In consequence, one can define two derivatives of the variable field Qtdtat any point of the bio-structure fractal curve:

d+Qtdtdt=limΔt0+Qt+ΔtΔtQtΔtΔtdQtdtdt=limΔt0QtΔtQtΔtΔtΔtE1

The “+” sign corresponds to forward processes of bio-structure’s structural units, while the “−” sign correspond to the backwards ones;

  1. (iv) The differential of the spatial coordinate field dXitdtby means of which we can describe the bio-structure dynamics, is expressed as the sum of the two differentials, one of them being scale resolution independent (differential part d±xitand the other one being scale resolution dependent (fractal part d+ζit)i.e.,

d±Xit=d±xit+d±ζitE2
  1. (v) The non-differentiable part of the spatial coordinate field, by means of which we can describe the bio-structure dynamics, satisfies the fractal equation [14]:

d±ζitdt=λ±idt1/DFE3

where λ±iare constant coefficients that help to specify the fractalization type which describes the bio-structure dynamics. Also, DFdefines the fractal dimension of the bio-structure non-differentiable curve.

In this way, the bio-structure processes imply dynamics on geodesics having different fractal dimensions. This variety of fractal dimensions of the bio-structure geodesics is a result of the bio-structure’s structure. For DF=2, quantum type processes are generated in bio-structure dynamics [15]. For DF<2, correlative type processes are induced and for DF>2non-correlative type processes are generated [6, 12, 13].

  1. (vi) The differential time reflection invariance of any dynamical variable of the bio-structure is recovered by combining the derivatives d+/dtand d/dtin the non-differentiable operator

d̂dt=12d++ddti2d+ddtE4

This is a natural result of the complex prolongation procedure applied to bio-structure dynamics [14, 16]. Applying now the non-differentiable operator to the spatial coordinate field, by means of which we can describe the bio-structure dynamics, yields the complex velocity field of the bio-structure.

V̂i=d̂Xidt=ViiUi,i=1E5

with

Vi=12d+Xi+dXidt,Ui=12d+XidXidtE6

The real part Viof the bio-structure complex velocity field is differentiable and scale resolution independent (differentiable velocity field). The imaginary part Uiis non-differentiable and scale resolution dependent (fractal velocity field).

  1. (vii) If we have no external constraint, one can find an infinite number of fractal curves (geodesics) relating any pair of points. This happens on all scales of bio-structure dynamics. Then, in the fractal space of the bio-structure, all the structural units are substituted with the geodesics themselves so that any external constraint can be interpreted as a selection of bio-structure geodesics. The infinity of geodesics in the bundle, their non-differentiability and the two values of the derivative imply a generalized statistical fluid-like description of the bio-structure dynamics. Then, the average values of the bio-structure variables must be considered in the previously mentioned sense, so the average of d±Xiis

d±Xid±xiE7

with

d±ζi=0E8

The previous relation (8) implies that the average of the fractal fluctuations is null.

  1. (viii) One can describe the bio-structure dynamics by means of a scale covariant derivative. Its explicit form can be obtained as follows. We assume that the bio-structure fractal curves are immersed in a 3-dimensional space. We also suppose that Xiis the spatial coordinate field of a point on this fractal curve. Let us also consider a variable field QXitand the following Taylor expansion up to the second order

d±QXit=tQdt+iQd±Xi+12lkQd±Xld±XkE9

These relations are valid at any point and more for the points Xi on the bio-structure fractal curve which we have selected in Eq. (9). From here, forward and backward average values for bio-structure variables from Eq. (9) become

d±Q=tQdt+iQd±Xi+12<lkQd±Xld±Xk>E10

We suppose that the average values of the all variable field Q and its derivatives coincide with themselves and the differentials d±Xiand dtare independent. Therefore, the average of their products coincides with the product of averages. Consequently, Eq. (10) becomes

d±Q=tQdt+iQ<d±Xi>+12lkQd±Xld±XkE11

Even the average value of d±ζiis null, for the higher order of d±ζithe situation can still be different. Let us focus on the averages d±ζld±ζk. Using Eq. (3) we can write

d±ζld±ζk=±λ±lλ±kdt2DF1dtE12

where the sign + corresponds to dt>0and the sign – corresponds to dt<0.

Then, Eq. (11) takes the form:

d±Q=tQdt+iQ<d±Xi>+12lkQd±xld±xk±12lkQλ±lλ±kdt2DF1dtE13

If we divide by dtand neglect the terms that contain differential factors (for details, see the method from [13, 14]) we obtain:

d±Qdt=tQ+ν±iiQ±12λ±lλ±kdt2DF1lkQE14

where ν+i=d+xidt,νi=dxidt.

These relations also allow us to define the operators

d±dt=t+ν±ii±12λ±lλ±kdt2DF1lkE15

Using Eqs. (4), (5), and (15), let us calculate the differentiable operator

d̂Qdt=tQ+V̂iiQ+14dt2DF1DlklkQE16

where

Dlk=dlkid¯lkdlk=λ+lλ+kλlλk,d¯lk=λ+lλ+k+λlλkE17

Eq. (16) also allows us to define the covariant derivative in the bio-structure dynamics

d̂dt=t+V̂ii+14dt2DF1DlklkE18

Let us now consider the principle of scale covariance (the physics laws – bio-structure dynamics specific – are invariant with respect to scale transformations) and postulate that the passage from the classical (differentiable) physics to the fractal (non-differentiable) physics can be implemented by replacing the standard time derivative d/dtby the non-differentiable operator d̂/dt. In this way, this operator has the role of a scale covariant derivative. More precisely, it is used to write the bio-structure dynamics fundamental equations in the same form as in the classic (differentiable) case. In these conditions, applying the operator (18) to the complex velocity field (5), with no external constraint, the bio-structure geodesics take the form:

d̂V̂idt=tV̂i+V̂llV̂i+14dt2DF1DlklkV̂i=0E19

This means that the local acceleration tV̂i, the convection V̂llV̂iand the dissipation DlklkV̂i, make their balance at any point of the bio-structure fractal curve. Moreover, the presence of the complex coefficient of viscosity-type 41dt2DF1Dlkin the bio-structures dynamics specifies that it is a rheological medium. So, it has memory, as a datum, by its own structure.

If the fractalization is achieved by Markov type stochastic processes, which involve Lévy type movements [9, 13, 17] of the bio-structure structural units, then:

λ+iλ+l=λiλl=2λδilE20

where δilis the Kronecker’s pseudo-tensor.

Under these conditions, the equation of bio-structure geodesics takes the simple form

d̂V̂idt=tV̂i+V̂llV̂idt2DF1DlklkV̂i=0E21

or more, by separating the motions on differential and fractal scale resolutions,

d̂VDidt=tVDi+VDllVDiVFl+λdt2DF1llVFi=0d̂VFidt=tVFi+VDllVFi+VFl+λdt2DF1llVDi=0E22

Using the standard procedure from [18], the bio-structure dynamics at fractal scale resolution can be described by means of the following equations:

tVFi+VFllVFi=λdt2DF1llVFiE23
iVFi=0E24

Eq. (23) corresponds to the specific impulse conservation law at fractal scale resolution, while Eq. (24) corresponds to the states density conservation law at fractal scale resolution (we consider that the density of the bio-structure at fractal scale resolution is constant – incompressible bio-structure).

Since this equation system is non-linear, one could find relatively difficult finding the solutions for these equations [19, 20]. However, in the particular case of a stationary flow in a plane symmetry xy, there is an analytical solution of this system. Then, for VF=VxVy0, Eqs. (23) and (24) take the form:

VxVxx+VyVxx=λdt2DF12Vxy2E25
Vxx+Vyy=0E26

The boundary conditions of the flow are:

limy0Vyxy=0,limy0Vxy=0,limyVxxy=0E27

and the flux momentum per length unit is constant

Θ=ρ+Vx2dy=const.E28

Using the method from [18, 19, 20] for solving Eqs. (25) and (26), with the conditions (27) and (28), the following solutions result:

Vx=1.5Θ6ϱ23λdt2DF1x1/3sech20.5yΘ6ϱ13λdt2DF1x2/3E29
Vy=4.5Θ6ϱ233λdt2DF1x1/3[yΘ6ϱ13λdt2DF1x23sech20.5yΘ6ϱ13λdt2DF1x23tanh0.5yΘ6ϱ13λdt2DF1x23E30

Relations (29) and (30) suggest that the bio-structure velocity field is highly non-linear through topological solitons of kink type (tanh), various non-topological solitons of breather type (sech2), and through topological – non-topological soliton mixtures of kink-breather type (sech2-tanh). Given the structural complexity of the bio-structure (which is given by its various structural units, that retains their own velocity field) an accurate way of writing relations (29) and (30) will be the one in which we assign indexes for each component.

For y=0, we obtain in relation (29) the flow critical velocity of the bio-structure in the form

Vxxy=0=Vc=1.5Θ6ϱ23λdt2DF1x1/3E31

while taking into account (31), relation (28) becomes

Θ=ρ++Vx2xydy=dc+dcVc2x0dy,E32

so that the critical cross section of the strains lines tube of the bio-structure is given by:

dcxy=0=Θ2ρVc2=2.42λdt2DF1x23ρΘ1/3E33

Relations (29) and (30) can be strongly simplified if we introduce the normalized quantities

ζ=xx0,η=yy0,u=Vxw0,v=Vyw0,Ω=Θ6ϱ23w0λdt2DF1x01/3,ω=Θ6ϱ23y0λdt2DF1x02/3E34

where x0,y0,w0are specific lengths and the specific velocity, respectively, of the laminar flow of the bio-structure. It results that

uζη=1.5Ωζ13sech20.5Ωωηζ23,E35
vζη=4.52/3313Ωζ13ωηζ23sech20.5Ωωηζ23tanh0.5Ωωηζ23E36

We present in Figures 1a, b and 2a, b the dependence of the normalized velocity field u on the normalized spatial coordinates ξ, η for various nonlinearity degrees (ω = 0.3; 6). The results showcase that the velocity field on the bio-structure flow direction (ξ) is affected in a weak manner by the nonlinearity degree (the velocity always decreases on the flow axes regardless of the nonlinearity degree). Also, the bio-structure flow direction (η) is strongly affected. Bio-structure flow starts from constant values on the η axes and with the increase of ω, preferential bio-structure flow direction can be identified.

Figure 1.

The dependence of the normalized velocity field u on the normalized spatial coordinates ξ, η for the non-linearity degree ω = 0.3: (a) 3D representation; contour plot (b).

Figure 2.

The dependence of the normalized velocity field u on the normalized spatial coordinates ξ, η for the non-linearity degree ω = 6: (a) 3D representation; contour plot (b).

In Figures 3a, b and 4a, b the dependences of the normalized velocity field v on the normalized spatial coordinates ξ, η for various non-linearity degrees (ω = 0.3; 6) are represented. For small non-linearity degrees, the variations (increase/decrease) of the velocity field have similar behaviors on both directions (ξ, η), while for higher values of the non-linearity degree these variations are only focused on a single direction (ξ).

Figure 3.

The 3D representation (a) and the contour plot (b) of the normalized velocity field v on the normalized spatial coordinates ξ, η for the nonlinearity degree ω = 0.3.

Figure 4.

The 3D representation (a) and the contour plot (b) of the normalized velocity field v on the normalized spatial coordinates ξ, η for the non-linearity degree ω = 6.

Taking the above into account, the force that the bio-structure will exercise to the walls of the flow vessels is of great importance for the understanding of arterial occlusion and other circulatory system diseases.

In our case the normalized force is given by the relation:

fζη=ηuζv=1.5Ωsech20.5ωηζ23tanh0.5ωηζ23ωζ1ζ130.9Ω23ωηsech20.5ωηζ23ζ53+0.66ω2η2sech20.5ωηζ23tanh0.5ωηζ23ωζ73+0.331tanh20.5ωηζ23ωηζ53+0.3Ωωηsech20.5ωηζ23ζ23tanh0.5ωηζ2333/2ζ43E37

In Figures 5a, b and 6a, b the normalized force field evolution on the two-flow direction (ξ, η) for various non-linear degrees is represented. It results that with the increase of the non-linearity of the bio-structure the force toward the walls increases.

Figure 5.

The dependence of the normalized force field F of a bio-structure flow on the vessels, on the normalized spatial coordinates ξ, η for two resolution scales: 3D representation (a); contour plot (b) for the non-linearity degree ω = 0.3.

Figure 6.

The dependence of the normalized force field F of a bio-structure flow on the vessels, on the normalized spatial coordinates ξ, η for two resolution scales: 3D representation (a) and contour plot (b), for the non-linearity degree ω = 6.

3. Discussions

The theory proposed in this chapter explains from a fractal viewpoint the atherogenesis process [21], basically “molding” to the classical anatomical and histopathological descriptions and completely respecting the process they postulate. In consequence, the fractal physics model sustains already accumulated morpho-pathological information and research. There are plenty of electronic and optical microscopy images that describe the spatial-temporal hologram of the phenomenon; we can thus discuss about the non-fractal – fractal and microscopic – macroscopic translation through holographically reproducible auto-similarity [21]. In this way, we affirm that fractality is the mathematical and semantic quintessence for defining atherogenesis, a process that can be physically characterized by fractal physics. This physics becomes in this situation more of a component rather than an explanation for the complex biological system represented by the atheroma plaque [22, 23].

In what concerns the recovery of such biological diseases, there are a huge number of techniques. We recall that external electrical stimulation can cause changes in the bio-structure vessels. Although atherosclerosis cause vasodilatation in the affected area and bio-structure flow remains unchanged for an extended period of time, the vascular wall stiffness will increase the pulse pressure. The purpose of the study developed in [24, 25] was to measure the effects of electrical stimulation (ES) on bio-structure flow and bio-structure pressure. All subjects received electrical stimulation at intensity sufficient to produce torque equal to 15% of the predetermined maximal voluntary contraction of their right quadriceps femor is muscle. The conclusions were that the increase in bio-structure flow occurred within 5 min after the onset of ES and dropped to resting levels within 1 min after a 10-min period of ES [25]. Kinesiotherapy or Kinesitherapy or kinesiatrics, is the therapeutic treatment of disease by passive and active muscular movements (as by massage) and of exercise [26]. From the physiotherapeutic viewpoint, an efficient treatment is directed toward improving bio-structure flow and also decreasing the disparity between the demand for bio-structure and its supply [22]. An effective vascular rehabilitation training program for improving walking efficiency and vascular remodeling in patients with diabetic atherosclerosis suffering from intermittent claudication could be a supervised treadmill walking exercise combined with Allen-Burger exercises [23].

4. Conclusions

The present chapter proposes a fractal model for the dynamics analysis of bio-structure flows. The fractal hydrodynamic equations were obtained and applied for the laminar flow of bio-structure.

A second application was proposed for bio-structure flow, and of cholesterol deposition on the vessel walls. The results revealed the directional flow toward the walls. This could explain in our opinion the thickening effect which is one of the sources of arteriosclerosis. Moreover, our model imposes redefinition of “good” and “bad” cholesterol (which are traditionally associated with HDL and LDL respectively); instead they should be replaced by the following notions: specific cholesterol entities, associated with a certain non-differentiable curve, that have a major endothelial impact and specific cholesterol entities which have no or low endothelial impact.

There is currently a great number of works describing matter organization and behavior in all of its variations, from which we mention [27]. We consider that our bio-structure flow model could also be used to further development in the study of other complex systems dynamics (such as pulmonary and metabolic diseases or environmental systems). Moreover, possible therapeutic treatments can be developed, e.g., new drug release mechanisms.

Conflict of interest

The authors declare no conflict of interest.

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Emil Anton, Anna Gavrilut, Maricel Agop and Daniel Timofte (July 18th 2018). Non-Linear Behaviours in the Dynamics of Some Biostructures, Nonlinear Systems - Modeling, Estimation, and Stability, Mahmut Reyhanoglu, IntechOpen, DOI: 10.5772/intechopen.74435. Available from:

chapter statistics

273total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Soliton-Like Solutions in the Problems of Vibrations of Nonlinear Mechanical Systems: Survey

By Yury A. Rossikhin and Marina V. Shitikova

Related Book

First chapter

On Nonoscillatory Solutions of Two-Dimensional Nonlinear Dynamical Systems

By Elvan Akın and Özkan Öztürk

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us