Open access peer-reviewed chapter

Complex Systems with Self‐Elimination of Dissipation with Implication in Bio‐Structural Behavior Via Nondifferentiability

By Maricel Agop, Decebal Vasincu, Daniel Timofte, Elena Simona Bacaita, Andrei Agop and Stefan Andrei Irimiciuc

Submitted: October 13th 2016Reviewed: February 17th 2017Published: July 26th 2017

DOI: 10.5772/67939

Downloaded: 408

Abstract

In the present chapter, we show that the use of the nondifferentiable mathematical procedures, developed in the Scale Relativity Theory with constant arbitrary fractal dimension, simplifies very much the dynamics analyses in the case of complex systems. By applying such a procedure to various complex systems dynamics (biological structures, ablation or discharge plasmas, etc.), we are able to observe that it starts from a steady (oscillating state) and as the external factor is varied the system undergoes significant changes. The systems evolve asymptotically through various transition, toward a chaotic regime (like bifurcations or intermittencies), but never reaching it. Another important reveal from the study of the system’s dynamics was the presence of various steady states depending on the resolution scale at which the theoretical investigations are performed.

Keywords

  • complex systems
  • fractal model
  • chaos
  • biological systems

1. Introduction

Complex systems, encountered in human societies, neural networks, the Internet, ecosystems, biological evolution, stock markets, economies, and many others, represent interdisciplinary research topics that have been studied by means of some fundamental theories, especially from physics and computer simulation. These systems are composed of a high number of individual components and their evolution cannot be predicted simply by analyzing the individual behavior of their elements or by adding their behavior. Their global evolution is decisively influenced by the manner in which the elements relate to each other, leading to emergence, self‐organization, and adaptability [13]. Correspondingly, the theoretical models corresponding to the complex systems dynamics become sophisticated [2, 3].

The models can be strongly simplified by taking into account that the complexity of interaction process imposes various temporal resolution scales, and the pattern evolution imposes different degrees of freedom [4]. Thus, new theoretical models can be developed by admitting that the complex systems that display chaotic behavior are recognized to acquire self‐similarity in association with strong fluctuations at all possible space‐time scales [5, 6]. Then, for large temporal scales compared to the inverse of the highest Lyapunov exponent [710], instead of deterministic trajectories we will work with a collection of potential trajectories and instead of definite positions with probability density. One interesting example is the case of collisions in complex system, where the dynamics of the particles can be described by nondifferentiable curves.

Therefore, in order to build a theoretical model for the evolution of complex systems, the fundaments of a nondifferentiable physics must be defined. For this, the system complexity is substituted with the nondifferentiable character of physical quantities, idea that lies at the base of Scale Relativity Theory (SRT) [11, 12] and of the nonstandard Scale Relativity Theory (NSRT), a particular case of the first one for a fractal dimension arbitrary constant [13]. In the framework of the above two theories, assuming that the individual components of the complex system move on continuous, but nondifferentiable curves, named fractal curves, the evolution of complex systems can be described by physical quantities depending both on the well‐known space‐time coordinates and, besides that, on the space‐time scales resolution; thus, it can be considered fractal functions [11, 13]. Moreover, the complex system will act as a “fluid” without interactions because its individual components are reduced to and identified with their own trajectories, named geodesics.

Further, we show that the dynamics “control” of various complex systems (biological structures, ablation, discharge plasmas, etc.) can be realized by means of fractality.

2. Complex systems: mathematical model

Taking into consideration the above, the following consequences emerge regarding the complex systems dynamics [1113]:

  1. any trajectory of the individual components (fractal fluid line) is explicitly dependent on scale resolution δt;

  2. applying the substitution principle, δtdt, the scale resolution can be considered as an independent variable. We reserve the notation dt for the usual time as in the Hamiltonian complex system dynamics;

  3. the variables that appear in the dynamics of complex system are fractal, i.e., are functions that depends both on the space‐time coordinates and on the scale resolution. Then, in any point of the fractal curve, two derivatives of the variable field M(t,dt)as explicit functions of the two variables t and dt, can be defined:

    d+M(t,dt)dt=limΔt0+M(t+Δt,Δt)M(t,Δt)ΔtdM(t,dt)dt=limΔt0M(t,Δt)M(tΔt,Δt)ΔtE1

    The “+” sign indicates forward processes and the “–“sign to the backwards ones;

  4. the differential of the spatial coordinate field dXi(t,dt)is expressed as:

    d±Xi(t,dt)=d±xi(t)+d±ξi(t,dt)E2

    where d±xi(t)is the “classical part” (differentiable, scale resolution independent) and d±ξi(t,dt)is the “fractal part” (nondifferentiable, scale resolution dependent);

  5. the nondifferentiable part of the spatial coordinate field, by means of which we can describe the complex system dynamics, satisfies the fractal equation [5]:

    d±ξi=λ±i(dt)1/DFE3

    where λ±iare constant coefficients through which the fractalization type describing the complex system dynamics is specified and DFdefines the fractal dimension of the motion non‐differentiable curve.

    In our opinion, the complex system processes imply a dynamics on geodesics with various fractal dimensions. The variety of these fractal dimensions of the complex system geodesics comes as a result of complex system structure. Precisely, for DF=2, quantum type processes are generated in the complex system, for DF<2correlative type processes are induced and for DF>2noncorrelative type ones can be found (for details, see [8, 14]);

  6. the differential time reflection invariance of any dynamical variable of the complex system is recovered by combining the derivatives d+/dtand d/dtin the nondifferentiable operator

    ^t=12(d++ddt)i2(d+ddt)E4

    This is a natural result of the complex prolongation procedure applied to complex system dynamics [13, 15]. Applying now the nondifferentiable operator to the spatial coordinate field, by means of which we can describe the fractal fluid dynamics, yields the complex velocity field of the fractal fluid:

    Vi^=d^Xidt=VDiVFiE5

    with:

    VDi=12(d+Xi+dXidt)VFi=12(d+XidXidt)E6

    The real part, VDi, represents the “classical velocity” (differentiable, i.e., scale resolution independent), while the imaginary part, VFi, is the “fractal velocity” (nondifferentiable, i.e., scale resolution dependent), induced by fractality;

  7. any external constraint is equivalent with a selection of fractal fluid geodesics, whose nondifferentiability, together with the two values of the derivative imply a generalized statistical fluid‐like description. Then, the average of d±Xiis:

    d±Xi=d±xiE7

    with

    d±ξi=0E8

    The previous relation (8) means that the average of the fractal fluctuations is null;

  8. in order to determine the scale covariant derivative, which describes the fractal fluid dynamics, we will consider that the individual components of the complex fluid move in a three‐dimensional space and that Xiare the spatial coordinates of a point on their nondifferentiable curves. The second‐order Taylor expansion for a variable field M(Xi,t)is:

    d±M(Xi,t)=tMdt+iMd±Xi+12lkMd±Xld±XkE9

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

d±M=tMdt+iMd±Xi+12lkMd±Xld±XkE10

Assuming that the differentials d±Xi, dtare independent and the mediated values of all variables and their derivatives coincide with themselves, relation (10) becomes:

d±M=tMdt+iMd±Xi+12lkMd±Xld±XkE11

Using Eq. (3), for ijwe can write:

d±ξld±ξk=±λ±lλ±k(dt)(2/DF)1dtE12

where the signs “+” and “−” have the same significance as in relation (1).

Then, Eq. (11) takes the form:

d±M=tMdt+iMd±Xi+12lkMd±xld±xk+12lkM[λ±lλ±k(dt)(2/DF)1dt]E13

Using the method from [1113], which involves dividing by dt and neglecting the terms with differential factors, we obtain:

d±Mdt=tM+v±iiM+12λ±lλ±k(dt)(2/DF)1lkME14

with

v±i=d±xidt

These relations also allow us to define the operator:

d±dt=t+v±ii+12λ±lλ±k(dt)(2/DF)1lkE15

Considering Eqs. (4), (5), and (15), it results:

dM^dt=tM+Vi^iM+14(dt)(2/DF)1DlklkME16

where

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

Using Eq. (16), the scale covariant derivative in the fractal fluid dynamics is defined as:

d^dt=t+Vi^i+14(dt)(2/DF)1DlklkE18

3. Fractal fluid geodesics

Based on the scale covariance principle (physics laws are invariant in respect to all scale transformations), we will consider that the transition from classical to nondifferentiable physics is equivalent with the replacement of standard derivative d/dtwith the nondifferentiable operator d^/dt. In consequence, it will act as scale covariant derivative, applied for writing the equations of fractal fluid dynamics in similar forms to those from the classical (differentiable) physics. Thus, the operator (18) applied to the complex velocity (5), with no external constraint, implies the geodesics in the following form:

d^Vi^dt=tVi^+Vl^lVi^+14(dt)(2/DF)1DlklkVi^=0E19

This means that the local acceleration tVi^, the convection Vl^lVi^, and the dissipation DlklkVi^make their balance in any point of the nondifferentiable curve. Moreover, the presence of the complex coefficient, 4(dt)1(2/DF)1Dlk, indicates that the fractal fluid has rheological properties, such as memory, determined by its internal structure.

Considering that the fractalization is induced by Markov‐type stochastic processes, with individual components displacements of Lévy type [16, 17], then:

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

where δilis the Kronecker’s pseudo‐tensor.

In these conditions, Eq. (19) becomes:

d^Vi^dt=tVi^+Vl^lVi^iλ(dt)(2/DF)1llVi^=0E21

Separating the individual components movements on differential and fractal scale resolutions, Eq. (21) becomes equivalent with the equations set:

d^VDidt=tVDi+VDllVDi[VFlλ(dt)(2/DF)1l]lVFi=0d^VFidt=tVFi+VDllVFi[VFlλ(dt)(2/DF)1l]lVDi=0E22

4. Geodesics in the Schrödinger‐type representation

For irrotational motions, the complex velocity satisfies the condition:

εiklkVl^=0E23

where εiklis the Levi‐Civita pseudo‐tensor. Then, the complex velocity field becomes:

Vi^=2iλ(dt)(2/DF)1ilnΨE24

where for the moment, lnΨdefines the scalar potential of the complex velocity field. Substituting Eq. (24) in Eq. (22), it results:

d^Vi^dt=2iλ(dt)(2/DF)1{tilnΨi[2λ(dt)(2/DF)1(llnΨl)ilnΨ+λ(dt)(2/DF)1llilnΨ]}=0E25

or using the identities:

lllnΨ+ilnΨilnΨ=llΨΨi(llΨΨ)=2(llnΨl)ilnΨ+llilnΨE26

the equation:

d^Vi^dt=2iλ(dt)(2/DF)1i[tlnΨ2iλ(dt)(2/DF)1lllnΨΨ]=0E27

By integration up to an arbitrary factor, set to zero by a suitable choice of Ψphase, from Eq. (27) it results:

λ2(dt)(4/DF)2llΨ+iλ(dt)(2/DF)1tΨ=0E28

This is an equation of Schrödinger type, describing geodesics in the Schrödinger representation. It is reduced to the usual Schrödinger equation for motions on Peano‐type curves, DF=2[1214], at Compton scale, λ=ħ/2m0, where ħis the reduced Planck’s constant and m0the rest mass of the particle. In the presence of an external constrain, defined by a scalar potential U, Eq. (28) can be written in the form:

λ2(dt)(4/DF)2llΨ+iλ(dt)(2/DF)1tΨU2Ψ=0E29

The relation (29) corresponds to the equation of motion for the “one body problem” in the Schrödinger‐type representation of the nondifferentiable model.

The standard equation of motion for the “one body problem” in the Scale Relativity Theory [1214]:

D2llΨ+iDtΨU2Ψ=0E30

results from Eq. (29) for movements of the complex system particles on Peano‐type curves, i.e., DF=2and the correspondence λD, where Dis the coefficient of the fractal‐nonfractal transition from Scale Relativity Theory [1214].

From such a view, in order to obtain the equations of motion for the “two‐body problem” in the Schrödinger‐type representation of the nondifferentiable model, we will apply a procedure which involves the next steps:

  1. the equations of motion for the “two‐body problem” are written in the complex momentum representation:

    d^P1l^dt=1lΦ,d^P2l^dt=2lΦE31

    where

    P1l^=m1V1l^,P2l^=m2V2l^E32

    are the complex momenta, V1l^and V2l^are the complex velocities, m1and m2are the rest masses of the “bodies” and Φis the scalar potential of the interaction forces.

  2. let’s consider that the motions are irrotational, i.e.,

    V1l^=2iλ1(dt)(2/DF)11llnΨ1,V2l^=2iλ2(dt)(2/DF)12llnΨ2E33

    where lnΨ1and lnΨ2are the scalar potentials of the complex velocities associated to the physical objects.

  3. substituting Eq. (33) in Eq. (31) and following the method described above, the equations of motions of the “bodies” in the absence of direct contact are obtained in the form:

    2m1λ1(dt)(2/DF)1tΨ1=[2m1λ12(dt)(4/DF)21l1l+Φ]Ψ12m2λ2(dt)(2/DF)1tΨ2=[2m2λ22(dt)(4/DF)22l2lΦ]Ψ2E34

In the presence of the direct contact, these equations become:

2m1λ1(dt)(2/DF)1tΨ1=[2m1λ12(dt)(4/DF)21l1l+Φ]Ψ1+Γ12Ψ22m2λ2(dt)(2/DF)1tΨ2=[2m2λ22(dt)(4/DF)22l2lΦ]Ψ2+Γ21Ψ1E35

where Γ12and Γ21are the coupling coefficients of the “bodies.”

Since Ψ1and Ψ2have a direct physical significations only through as |Ψ1|2and |Ψ2|2probability densities, the deterministic trajectories are replaced by a collection of “potential routes.” In its turn, the concept of “definite position” is replaced by that of probability density. Moreover, the complex fluid particles may be reduced to and identified with their own trajectories (i.e., their geodesics) so that complex fluid should behave as a special “fluid” free of interactions—a fractal fluid. In such a conjecture, quantum‐type effects (tunneling effect, entanglement effect, etc.) can be extended to the macroscopic complex fluid.

5. Geodesics in the hydrodynamic‐type representation

If Ψ=ρeiS, with ρthe amplitude and Sthe phase of Ψ, the complex velocity (24) takes the forms:

VDi^=2iλ(dt)(2/DF)1iSVFi^=2iλ(dt)(2/DF)1ilnρE36

Substituting Eq. (36) into Eq. (21) in the presence of an external potential and separating the real and imaginary parts, up to an arbitrary factor, set to zero by a suitable choice of Ψphase, we obtain:

tVDi+(VDl^l)VDi^=l(Q+U)E37
tρ+i(ρVDl^)=0E38

with Q, the specific nondifferentiable potential:

Q=2λ2(dt)(4/DF)2llρρ=VFlVFl2+λ(dt)(2/DF)1iVFiE39

Eq. (37) represents the specific momentum conservation law, while Eq. (38) represents the states density conservation law. Eqs. (37)(39) define the fractal hydrodynamic model and imply the following:

  1. any individual component is in permanent interaction with the surrounding fractal medium;

  2. the fractal medium can be identified with a nondifferentiable fluid, described by Eqs. (37)(39);

  3. even if the fractal velocity VFi^does not reflect a real motion, it influence the specific momentum and energy transfer (confirmed by the its absence from the conservation law of states density); and

  4. most part of the energy is the form of kinetic and potential energy, other parts in other forms, but the overall energy is constant; the energy and the specific momentum conservation is the one that ensures reversibility and the presence of eigenstates, but does not confirm a Levy type motion in an external field.

6. Fractal behaviors

Since, in our model, the fractality plays an essential role in the dynamics of the “complex fluid,” next we will expand the model, admitting the following equations:

(VFiλ(dt)(2/DF)1l)lVFi=0E40
lVFl=0E41

The first equation presents the fact that “the fractal force” is null, while the second one presents the property of incompressibility of the fractal fluid.

6.1. Fractal laminar flow

To find the solutions for these equations can be a relatively difficult, due to the fact that this equation system is a nonlinear one [1820]. However, there is an analytical solution of this system, in the particular case of a stationary flow in a plane symmetry (x,y). In these circumstances, Eqs. (40) and (41), with VFi=Vi, take the form:

VxVxx+VyVxx=λ(dt)(2/DF)12Vxy2E42
Vxx+Vyy=0E43

where Vx=Vx(x,y)is the velocity along axis Ox, Vy=Vy(x,y)is the velocity along axis Oy. The boundary conditions of the flow are:

limy0Vy(x,y)=0,limy0Vxy=0,limyVx(x,y)=0E44

and the flux momentum per length unit is constant:

Θ=ρ+Vx2dy=const.E45

Using the method from [1820] for resolving Eqs. (42) and (43), with the limit conditions (44) and (45), the following solutions result:

Vx=[1.5(Θ6ρ)23][λ(dt)(2/DF)1x]13sech2[(0.5y)(Θ6ρ)13][λ(dt)(2/DF)1x]23E46
Vy=[4.5(Θ6ρ)23][3λ(dt)(2/DF)1x]13[[y(Θ6ρ)13][λ(dt)(2/DF)1x]23sech2[(0.5y)(Θ6ρ)13][λ(dt)(2/DF)1x]23tanh[(0.5y)(Θ6ρ)13][λ(dt)(2/DF)1x]23]E47

Relations (46) and (47) suggest that the fractal fluid velocity field is highly nonlinear by means of soliton and soliton‐kink‐type solutions. Given the structural complexity of the fluid (given by its various structural units, that retains their own velocity field), an accurate way of writing relations (46) and (47) will be the one in which we assign indexes for each component.

For y=0, we obtain in relation (46) the critical velocity of the flow in the form:

Vx(x,y=0)=Vc=[1.5(Θ6ρ)23][λ(dt)(2/DF)1x]13E48

while relation (45), taking into account Eq. (48), becomes:

Θ=ρ+Vx2(x,y)dy=dc+dcVc2(x,0)dyE49

so that the critical cross section of the strain line tube is given by:

dc(x,y=0)=Θ2ρVc2=2.42[[λ(dt)(2/DF)1x]23](ρΘ)13E50

Relations (46) and (47) can be strongly simplified if we introduce the normalized quantities:

ξ=xx0,η=yy0,u=Vxw0,v=Vyw0,Ω=(Θ6ρ)23w0[λ(dt)(2/DF)1x0]13,ω=(Θ6ρ)13y0[λ(dt)(2/DF)1x0]23E51

where x0,y0,w0are the specific lengths and the specific velocity, respectively, of the laminar flow of the fractal fluid. It results that

u(ξ,η)=1.5Ωξ1/3sech2(0.5Ωωηξ2/3)E52
v(ξ,η)=4.52/331/3Ωξ1/3[ωηξ2/3sech2(0.5Ωωηξ2/3)tanh(0.5Ωωηξ2/3)]E53

We present in Figures 1(a, b)4(a, b) the dependence of the normalized velocity field uon the normalized spatial coordinates ξ,ηfor various nonlinearity degrees (ω=0.1,0.2,0.5,5). The results show that the velocity field on the 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). On the other hand, the flow direction (η) is strongly affected. The flow starts from constant values on the ηaxis, and with the increase of ω, preferential flow direction can be identified.

Figure 1.

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

Figure 2.

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

Figure 3.

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

Figure 4.

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

Figures 5(a, b)8(a, b) represent the dependence of the normalized velocity field von the normalized spatial coordinates ξ,ηfor various nonlinearity degrees (ω=0.1,0.2,0.5,5). For small nonlinearity degrees, the variations (increase/decrease) of the velocity field have similar behaviors on both directions (ξ,η), while for higher values of the nonlinearity degree, these variations are only focused on a single direction (ξ).

Figure 5.

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

Figure 6.

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

Figure 7.

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

Figure 8.

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

6.2. Dynamics of laser ablation plasma assimilated as a fractal fluid

Due to the complexity of the interactions between the particles inside the plume as well as to those between the expanding plume and the background gas, it remains difficult to connect the laser‐matter interactions with the plume expansion as well with the interactions between the plume and the laser beam, and the buffer gas. In the case of short laser wavelength, the ablation process can be generally divided in several stages [21]. In the first stages, laser-target interactions, such as laser absorption by the target, target heating, occur. If the laser fluency surpasses the ablation threshold of the target, the laser beam evaporates and ionizes material, generating a plasma plume above the material surface. First, the ejected particles undertake strong collisions in a high‐density area near the target creating the Knudsen layer [22]. This leads to a directional evolution, orthogonal to the target surface. After the pulse ends, no more particles are ejected from the target surface. An adiabatic expansion of the plasma happens where the temperatures can be related to the dimensions of the plasma; the thermal energy is then converted into kinetic energy, with the plasma reaching very high expansion velocities [23, 24]. During the initial stages of plasma expansion, defined by high densities, the mean‐free path of the particles is short and the plasma dynamic can be considered as one of a continuum fluid. As the plasma expands, the thermal energy decreases at a quick rate; however, the drop diminishes over time since the energy is recovered in the processes of ions recombination. Once the ambient pressure increases, the plume evolution is determined by the interactions of the laser plasma and ambient gas [25].

Given the complexity of the elementary processes involved in the particle removal and plasma expansion presented above, the dynamic of the laser ablation plasmas can be considered a particular case for the fractal fluid previous presented. An exhaustive model is presented in [26, 27]. In such a conjecture, the evolution of the plasma plume can be seen as series of successive dynamics at various resolutions scales, described by a constant particle density. Figures 9(a, b)10(a, b) present the flow of a fractal fluid at two simultaneously different resolution scales. We observed that the field velocity udepends on both expansion direction and its values decrease exponentially. In Figure 9(a, b), it is observed a single maximum, generally attributed to specific plasma structure. The field velocity vpresents a similar evolution on the ξaxis while on the ηdirection we observe a split, described by two maxima, attributed to the splitting process often observed for laser‐produced plasmas [26, 27].

Figure 9.

The dependence of the normalized velocity field u on the normalized spatial coordinates ξ,η for two simultaneously resolution scales: 3D representation (a), the contour plot (b).

Figure 10.

The dependence of the normalized velocity field v on the normalized spatial coordinates ξ,η for two simultaneously resolution scales: 3D representation (a), the contour plot (b).

6.3. Blood flow dynamics through fractal fluid approach

Blood is a physical fluid that carries both substances needed by cells for body functioning and, also, metabolic waste products from cells. It consists of 45% blood cells (red and white globules, platelets) and 55% plasma (92% water and various substances in dissipated forms, i.e., protein, glucose, ions, minerals, hormones, carbon dioxide, and cholesterol). Among blood cells, most numerous are red globules containing hemoglobin, a protein with iron in its structure, which facilitates the transport of oxygen by creating ionic bonds with it. In contrast, carbon dioxide is mostly transported extracellular as bicarbonate ion transported in plasma.

Blood is kept in continuous motion in the blood vessels by the heart, through its muscle contraction. Oxygen is transported by the arterial blood to all body cells, from where the waste product, carbon dioxide, is taken by the venous blood and carried to the lungs to be exhaled.

From a physical‐chemical point of view, the blood is a suspension, a mixture of liquids, gases, and solids (cells).

In terms of anatomy and histology, blood is considered a specialized form of connective tissue, given its origin in the bones and the presence of potential molecular fibers in the form of fibrinogen.

Since the circulatory system has a fractal structure, it is expected like its functionality to be also fractal. This allows us to assimilate the dynamics of the blood flow with the one of the fractal fluid. In this context, although the velocity fields will remain the same as the one presented in Figures 18, it is of great importance for the understanding of arterial occlusion and other circulatory system diseases, the force that the fluid will exercise to the walls of the flow vessels.

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

F=ηuξυ=1.5Ωsech2(0.5ωηξ23)tanh(0.5ωηξ23)ξ1ξ13(0.9Ω(23ωηsech2(0.5ωηξ23)ξ53+0.66ω2η2sech2(0.5ωηξ23)tanh(0.5ωηξ23)ξ73+0.33[(1tanh2(0.5ωηξ23))ωη]ξ53)332)+0.3Ω(ωηsech2(0.5ωηξ23)ξ23tanh(0.5ωηξ23))332ξ43E54

In Figures 11(a, b)14(a, b), it is represented the normalized force field evolution on the two flow direction (ξ,η) for various nonlinear degrees. It results that with the increase of the nonlinearity of the fluid the force toward the walls increases. This can be a starting point for understanding the complexity of the mechanisms involved in the arterial occlusion.

Figure 11.

The dependence of the normalized force field F on the normalized spatial coordinates ξ,η for the nonlinearity degree ω=0.1: 3D representation (a), the contour plot (b).

Figure 12.

The dependence of the normalized force field F on the normalized spatial coordinates ξ,η for the nonlinearity degree ω=0.2: 3D representation (a), the contour plot (b).

Figure 13.

The dependence of the normalized force field F on the normalized spatial coordinates ξ,η for the nonlinearity degree ω=0.5: 3D representation (a), the contour plot (b).

Figure 14.

The dependence of the normalized force field F on the normalized spatial coordinates ξ,η for the nonlinearity degree ω=5: 3D representation (a), the contour plot (b).

The proposed theory has the advantage that it explains from a fractal point of view the atherogenesis process [28], basically “molding” to the classical anatomical and histopathological descriptions, completely respecting the process postulated by them. Thus, the fractal physics model represents a dynamic and novel argument for sustaining already accumulated morphopathological information and research. The electronic and optical microscopy images (Figure 15) describe the spatiotemporal hologram of the phenomenon; we can thus talk about the nonfractal‐fractal and microscopic‐macroscopic translation through holographically reproducible autosimilarity [28]. We can thus say that fractality represents the mathematical and semantic quintessence for defining atherogenesis, a process that can be physically characterized perfectly by fractal physics, physics becoming in this situation more of a component rather than an explanation for the complex biological system represented by the atheroma plaque [28].

Figure 15.

Endoluminal view—significant atherosclerotic plaque with thrombus (optical microscopy).

Regarding the recovery of such biological diseases, there are a vast number of techniques. We remind that external electrical stimulation can cause changes in the blood vessels. Although atherosclerosis cause vasodilatation in the affected area and blood flow remains unchanged for an extended period of time, the vascular wall stiffness will increase the pulse pressure [29]. Tracy et al. developed, in 1950, a study whose purpose was to measure the effects of electrical stimulation (ES) on blood flow and blood 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 femoral muscle. The conclusions were that the increase in blood 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 [30].

Kinesiotherapy or kinesitherapy or kinesiatrics is the therapeutic treatment of disease by passive and active muscular movements (as by massage and by exercise) [31].

From the physiotherapeutic point of view, the treatment is directed toward improving blood flow and toward decreasing the disparity between the demand for blood and its supply [32].

An effective vascular rehabilitation training program for improving walking efficiency and vascular remodeling in patients with diabetic atherosclerosis suffering from intermittent claudicating could be a supervised treadmill walking exercise combined with Allen‐Burger exercises [33].

7. Complex system with the self‐elimination of dissipation

If the fractalization is achieved by Markov‐type stochastic processes, the scale covariant derivative (18) takes the form:

d^dt=t+Vl^l±iλ(dt)(2/DF)1llE55

Postulating now the scale covariance principle the states density conservation laws become:

dρ^dt=tρ+Vl^lρ±iλ(dt)(2/DF)1llρE56

or more, separating the movement on the scale resolution:

tρ+Vllρ=0E57

for the differentiable scale resolution and:

VFllρ±λ(dt)(2/DF)1llρ=0E58

for the nondifferentiable scale resolution. From such a perspective, the fractal‐nonfractal dynamic transition of the state density can be obtained by summing Eqs. (57) and (58), taking the form:

tρ+(VDlVFl)lρ±λ(dt)(2/DF)1llρ=0E59

From here, by identifying the movements at the two resolution scales VDl=VFl, the fractal‐type diffusion equations become:

tρ±λ(dt)(2/DF)1llρ=0E60

Let us now use Eq. (60) to analyze the dynamic of an electron beam accelerated in a strong electric field which impinges onto a neutral medium. As a result of these interactions, ionizations are produced by: (i) the primary electrons (from the fascicle), which are accelerated by an external electrical field, αj, where αis the primary ionization coefficient and jis the fascicle current density and by (ii) secondary electrons which results from the direct ionization process, βjρe, with βthe secondary ionization coefficient and ρethe electron density. In this circumstance, we will further focus on the study of the dynamics induced only by the electronic branch through Eq. (60) written in the following form:

tρe+λ(dt)(2/DF)1llρe=αjβjρeE61

Since the previous dynamics implies a dimensional symmetry, Eq. (61) by means of substitutions:

τ=xνt,M=Λβjν2λ(dt)(2/DF)1,2R=Λβj,K=Λ,αj+βjρe=ΛqE62

becomes a damped oscillator‐type equation:

Mq¨+2Rq˙+Kq=0E63

Rewritten as:

p˙=2RMpKMq,q=pE64

Eq. (63) induces a two‐dimensional manifold of phase space type (p,q)in which pwould correspond to a momentum‐type variable and qto a “position” type. Then, the parameters M, R, and Kcan have the following significance:

  1. Mcharacterizes the “inertial” type effect through the connection with ionization processes (global ionization described by αj+βjand the local ones described by βj) as well as in respect with the fractal diffusion (λ(dt)(2/DF)1). All these are done with respect to a traveling wave‐type movement based on the self‐similar dynamic solutions (τ=xνt);

  2. Rrepresents the “dissipative” type effects through the connection with the ionization processes, the same as above;

  3. Krepresents the structural type efforts in connection with the ionization processes, consisting, in this case, only from global ionization effects αj+βj.

The second equation from Eq. (64) corresponds to the momentum definition. Eq. (64) is not a Hamiltonian system, because the matrix traces are not null, and, as consequence, the associated matrix is not an involution. The observation becomes more obvious if we write it in matrix form:

(p˙q˙)=(2RMKM10)(pq)E65

With M, R, Kconstant values, this matrix equation, written in an equivalent form, evidences the position of the energy, i.e., of the Hamiltonian. From Eq. (64) we can write:

12M(pq˙qp˙)=12(Mp2+2Rpq+Kq2)E66

which proves that the energy in its quadratic form, i.e., the right‐hand side of Eq. (66) is the variation rate of the physical action represented by the elementary area from the phase space. We would like to show here that the energy does not have to satisfy the conservation laws in order to act like a variation rate for the action. Another issue is the form of the conservation law, if it exists. For this, Eq. (66) will be written as a Riccati‐type differential equation:

w˙+w2+2μw+ω02=0,w=pq,μ=RM,ω02=KME67

The Riccati‐type equation (67) represents always a Hamiltonian system that describes a dynamic of harmonic oscillating type:

(p˙q˙)=(RMKM1RM)(pq)E68

This is a general characteristic of the Riccati‐type equations and of the Hamiltonian’s dynamic [34]. The 1‐differential form for the elementary area from the phase space, obtained from Eq. (68), is identical with Eq. (60). Considering that the energy does not conserve anymore, but another, more complicated, dynamics variable is conserved, Eq. (65), by integration becomes [35]:

12(Mp2+2Rpq+Kq2)exp{2RMKR2arctan(Mp+RqqMKR2)}=constE69

The energy conservation, in the classical meaning, occurs if either Ris null, either the movement in the phase space takes place on a line though origin, with the slope R/M.

Moreover, by comparison with the case of the thermal radiation, regarding the distribution function of a preestablished ensemble [36]:

P(r,w)=11+2rw+w2exp[2r(1r2)1/2arctanw(1r2)1/21rw]E70

with rthe correlation coefficient and w2=ε0/uthe ratio between the quanta of thermal energy and the reference energy u, relation (69) becomes:

Kq22=const1+2rw+w2exp[2r(1r2)1/2arctanw(1r2)1/21rw],w2=Mp2kq2,r2=RkE71

From here we can emphasize the statistic character of energy: the potential energy is constructed as a functional of a specific statistical variable, given by the ratio between the kinetic energy and the potential one of a local oscillator.

Through the correlation of all the oscillator‐type ensembles, a “quantization” procedure can be established, by using the condition:

P(r,1)=eε0kT=12(1+r)exp[2r(1r2)1/2arctan(1+r1r)1/2]

where kis the Boltzmann constant and Tthe characteristic temperature of the thermal radiation (for details, see [37]). Figure 16(a, b) presents the quantization” procedure through correlation of all statistical ensembles associated with “local oscillators.”

Figure 16.

The probability function: in (r,w) space (a), and in the contour plot representation (b).

It shows explicitly the connection between the “quanta” and the statistical correlation of the process represented by the thermal radiation. Moreover, it also specifies the expression of the “quanta” in the small correlation limit, i.e., for r0, as ε0kTln2.

Thus, in such limit, the quanta and, implicitly, the frequency υ(through ε0=hυ, where his the Plank constant) are proportional with the “color” temperature.

Since we are focused on identifying these dissipative forces, we will try and present a physical significance for the Riccati equation (67) and of the associated Hamiltonian (68). For this, we observe that Eq. (63) is the expression of a variation principle:

δt0t1Ldt=0E72

regarding the Lagrangian:

L(q,q˙,t)=12(Mq˙2Kq2)e2RMtE73

This represents the Lagrangian form of a harmonic oscillator with explicit time‐dependent parameters. The Lagrangian integral defined on a finite interval [t0,t1]is the physical action of an oscillator on that specific time interval, describing the difference between the kinetic energy and the potential one. In order to obtain Eq. (63), it is necessary to consider the variation of this action under explicit conditions in such way that the variance of the coordinate at the interval extremes is null:

δq|t0=δq|t1=0E74

In order to obtain a closed trajectory, a supplementary condition must be imposed; for example, that the coordinates values at the interval extremes are identical:

q(t0)=q(t1)E75

Moreover, if this trajectory is closed in the phase space, the same condition results also for velocities.

Let us focus now on the movement principle. The Lagrangian is defined until an additive function which needs to be derivative in respect with the time of another function. The procedure is largely used in theoretical physics by defining the gauge transformation. In our case we define a gauge transformation in which the Lagrangian is a perfect square. This is known and explored in the control theory [34]. The procedure consists in adding to the Lagrangian the following term:

12ddt(we2RMtq2)E76

where wis a continuous function in time, so that the Lagrangian is a perfect square. The function variation is null, only due to the conditions presented in Eq. (75) and, therefore, the motion equation does not change. The new Lagrangian, written in relevant coordinates, takes the form:

L(q,q˙,t)=12Me2RMt(q˙+wMq)2E77

with the condition that wsatisfies the following Riccati‐type equation:

w˙1Mw2+2RMwK=0E78

The Lagrangian from Eq. (77) can be considered here as representing the whole system energy. As before, there is a relationship between the Riccati equation (78) and the Hamiltonian dynamic. Henceforth, we will find a similar relation to the one presented in Eq. (68):

(η˙ξ˙)=(RMKM1RM)(ηξ),w=ηξE79

In this case, our system is obviously a Hamiltonian one. Thus, we can identify the wfactors with the phase space coordinates. Eq. (78) specifies the fact that wis a dissipation factor (a mass variation rate, for the variable mass case). It is important to find the most general solution of this equation. Carena and Ramos presented a modern approach to integrate a Ricccati equation [38]. Applied to our case, it is enough to note the complex numbers as:

w0R+iMΩ,w0*RiMΩ,Ω2=KM(RM)2E80

The roots of the quadratic polynomial from the right‐hand part of Eq. (78) are two constant solution of the equation. Being constant, their derivative is null; thus, the polynomial is also null. In order to avoid this situation, we first perform the homographic transformation:

z=ww0ww0*E81

In these circumstances, it results that through direct determination, zis a solution of a linear and homogeneous first‐order equation:

z˙=2iΩzz(t)=z(0)e2iΩtE82

Hence, if we express the initial condition z(0)in a right manner, we can obtain the general solution of Eq. (78) by applying an inverse transformation to Eq. (81). We find:

w=w0+re2iΩ(ttr)w0*1+re2iΩ(ttr)E83

where rand trare two real constants which are characterizing the solution. Using Eq. (78) we can write the same solution in real terms:

z=R+MΩ(2rsin[2Ω(ttr)]1+r2+2rcos[2Ω(ttr)]+i1r21+r2+2rcos[2Ω(ttr)])E84

This relation shows frequency modulation through a Stoler transformation [39], which leads to the complex representation of this parameter.

The theoretical model proposed here analyzes the dynamic of charged particles in a plasma discharge where there is a strong flux of electrons from one plasma structure to another. Basically, the electrons dynamics are described using a forced damped oscillating system, with the aim to investigate the response of the global discharge current to different changes in resolution scale, oscillation frequency and dapping coefficient. Since our mathematical approach is sensitive to the changes in the resolution scales, we plotted in Figure 17(a–f) the 3D maps and the corresponding contour plot representations of the discharge current as functions of time and oscillation frequency, for a fixed value of the damping constant: 0,1. We observe that, for small resolutions, the current is described by a simple oscillatory regime, and, as the frequency resolution scale increases, we notice the appearance of some patterns. The patterns become denser and are foreshadowing the presence of some modulation in the oscillating frequency.

Figure 17.

3D representation of the discharge current for different resolution scales: 0–1 (a), 0–10 (b), 0–15 (c), 0–20 (d), 0–25 (e), 0–30 (f) with respect to the oscillation frequency and the corresponding contour plot representations.

The damping of the oscillatory state describes the losses through dissipative mechanisms. In order to study the effects induced by these mechanisms on the global current, we pointed two different oscillation frequencies, and observe the temporal response at different values of the damping coefficient. The results can be seen in Figure 18(a, b).

Figure 18.

(a and b) Discharge current temporal traces for various values of the damping coefficient extracted for two different frequencies scales.

We can identify competing oscillatory behaviors, described by two oscillation frequencies with comparable amplitudes. As the damping increases, the ratio between the two changes oscillation frequencies, and in the end the system oscillates on a single frequency.

The effect of the forced oscillations, which can be attributed to one of the plasma structures, is presented in Figure 19, where the discharge current for a fixed value of the damping coefficient and various values of the forced oscillations is represented. The systems seemingly start from a state described by period doubling, and goes through frequency modulation as the values increase. The important aspect is that the oscillation frequencies found for the current are not the ones induced through forced oscillatory system. This means that the system although forced to get on a specific state, it will define its own dynamic influenced, but not determined, by the external parameters.

Figure 19.

Discharge current temporal traces for different values of the oscillation frequency at a constant damping coefficient (r=0.1).

Finally, we investigated the evolution of our system with the increase of the control parameter ω. We investigate for a large range of values the evolution of the system and we observe (Figure 20a) that the system starts from double period state (see Figure 19), which arises periodically as the values of the control parameter increases. Also, the amplitude of such oscillations increases of about 10 times. In Figure 20(b), the bifurcation map is presented, where we observe that our system starts from steady state (double period state) and passes toward a chaotic one (ω=7,22,44), but never reaches it.

Figure 20.

The evolution of the double‐period state (a) and the characteristic bifurcation map (b) with relation to the control parameter ω.

Thus, the evolution of complex system, in our case a plasma, not only intends to reach a chaotic state, but each transition back to the steady state leads to an increase in the system energy, which can be seen in the increase of the amplitude.

8. Conclusions

This chapter proposes a new approach for the analysis of dynamics in complex systems based on fractal models, based on the assumption that the “individual components” motion takes place on continuous, but nondifferentiable curves. Consequently, the standard and the “exotic” properties of complex system flow at meso and nano scale can be highlighted.

In this framework, we analyze possible implications of complex fluids in blood flow dynamics. We can affirm that the fractal physics can explain the initiation phase of arteriosclerosis, overcoming the common expression found in classical medicine, “at a certain moment,” which lacks a precise temporal‐spatial definition. Moreover, our model imposes redefinition of “good” and “bad” cholesterol, traditionally associated with HDL, respectively LDL. Instead, they should be replaced by the following notions: specific cholesterol entities, associated with a certain nondifferentiable curve, that have a major endothelial impact, i.e., HDL, and specific cholesterol particles which have no or low endothelial impact, i.e., LDL.

We consider that our complex system flow model could also be used for the further development of the study of other complex systems dynamics, such as pulmonary and metabolic diseases or environmental systems.

The evolution of a complex system (plasma discharge or biological systems) was discussed as nonconservative and presented in a simpler way (i.e., through a forced damped oscillating system). Here, the driving force is used to simulate the response of the immune system, while the damping covers the energy loses during the evolution of the tumor cells. Once such a system is chosen, by simply investigating its dynamics, we were able to observe that the systems start from a steady (oscillating state), and as the external factor is varied, the system undergoes significant changes. The system evolves asymptotically through various transition, toward a chaotic regime (like bifurcations or intermittencies), but never reaching it. Another important reveal from the study of the system’s dynamics, was the presence of various steady states depending on the resolution scale at which the theoretical investigations are performed.

Our study not only shows the chaotic evolution tendencies of the complex system, but also reveals that there is a high probability for the presence of steady states. These states arise during the system’s evolution toward chaos, and could be related to dormant stages of the disease and the transitions toward chaos correspond to the uncontrollable growth. The fact that the system would still reach a healthy (steady) state followed by a chaotic evolution couple with an increase in energy would connect to some rebound of cancer in some cases.

© 2017 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

Maricel Agop, Decebal Vasincu, Daniel Timofte, Elena Simona Bacaita, Andrei Agop and Stefan Andrei Irimiciuc (July 26th 2017). Complex Systems with Self‐Elimination of Dissipation with Implication in Bio‐Structural Behavior Via Nondifferentiability, Fractal Analysis - Applications in Health Sciences and Social Sciences, Fernando Brambila, IntechOpen, DOI: 10.5772/67939. Available from:

chapter statistics

408total 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

The Fractal Analysis of the Images and Signals in Medical Diagnostics

By Tayurskii Dmitrii Albertovich and Rusanova Inna Aleksandrovna

Related Book

First chapter

Applications of Radial Basis Function Schemes to Fractional Partial Differential Equations

By Carlos Alberto Torres Martínez and Carlos Fuentes

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