Theoretical investigations on a variablecoefficient generalized forced–perturbed Korteweg–de Vries–Burgers model for a dilated artery,blood vessel or circulatory system with experimental support

Xin-Yi Gao ,Yong-Jiang Guo and Wen-Rui Shan

1 State Key Laboratory of Information Photonics and Optical Communications,and School of Science,Beijing University of Posts and Telecommunications,Beijing 100876,China

2 College of Science,North China University of Technology,Beijing 100144,China

Abstract Recent theoretical physics efforts have been focused on the probes for nonlinear pulse waves in,for example,variable-radius arteries.With respect to the nonlinear waves in an artery full of blood with certain aneurysm,pulses in a blood vessel,or features in a circulatory system,this paper symbolically computes out an auto-B?cklund transformation via a noncharacteristic movable singular manifold,certain families of the solitonic solutions,as well as a family of the similarity reductions for a variable-coefficient generalized forced–perturbed Korteweg–de Vries–Burgers equation.Aiming,e.g.,at the dynamical radial displacement superimposed on the original static deformation from an arterial wall,our results rely on the axial stretch of the injured artery,blood as an incompressible Newtonian fluid,radius variation along the axial direction or aneurysmal geometry,viscosity of the fluid,thickness of the artery,mass density of the membrane material,mass density of the fluid,strain energy density of the artery,shear modulus,stretch ratio,etc.We also highlight that the shock-wave structures from our solutions agree well with those dusty-plasma-experimentally reported.

Keywords: dynamics in blood-filled artery or blood vessel,variable-coefficient generalized forced–perturbed Korteweg–de Vries–Burgers equation,solitons with experimental support,B?cklund transformation and similarity reductions,singular manifold

1.Introduction

1.1.Background around us and our consideration

Nonlinear waves are physically and currently interesting [1–3].Physical studies on the pulse waves in the human arteries have started from the ancient times until the present,while recent theoretical efforts have been focused on the probes for nonlinear pulse waves in the variable-radius arteries [4–16].

We hereby consider a variable-coefficient generalized forced–perturbed Korteweg–de Vries–Burgers equation:

where v(x,t),ρ4(x,t) and ρ6(x,t) are all the real differentiable functions of the variables x and t,the subscripts represent the partial derivatives,while ρ1(t) ≠0,ρ2(t),ρ3(t) ≠0 and ρ5(t)are all the real differentiable functions of t.As seen below,e.g.,for an artery full of blood with an aneurysm,v can be the dynamical radial displacement superimposed on the original static deformation from an arterial wall,t and x can be the stretched coordinates,respectively,related to the axial coordinate and to both the axial coordinate and time parameter,ρi's (i=1,...,6) can be linked with the axial stretch of the injured artery,blood as an incompressible Newtonian fluid,radius variation in the axial direction or aneurysmal geometry,viscosity of the fluid,thickness of the artery,mass density of the membrane material,mass density of the fluid,strain energy density of the artery,shear modulus,stretch ratio,etc.

The following special cases guide ρ1(t) ≠0 and ρ3(t) ≠0,while ρ2(t),ρ4(x,t),ρ5(t) and ρ6(x,t) are not restricted.

1.2.Special cases of equation (1)

There have been some special cases of equation (1) as follows:

· whenρ1(t)=,ρ2(t)=,ρ3(t)=,ρ4(x,t)=andρ6(x,t)=a forced–perturbed Korteweg–de Vries–Burgers equation with variable coefficients for the nonlinear waves in an artery full of blood with a local dilatation (representing a certain aneurysm) [15,16]3Reference [16] has discussed the situation with =0.

where v(x,t) is the wave amplitude,t and x are the scaled‘time’ and scaled ‘space’,ρ1(t),ρ3(t),ρ5(t) as well as the real functionsandrepresent the variable coefficients of the nonlinear,dispersive,perturbed,dissipative and external-force terms,respectively [17];4

· when ρ2(t)=0,a variable-coefficient Korteweg–de Vries equation for the pulse waves in a blood vessel,a circulatory system or a fluid-filled tube [18–23]

where v(x,t) is the wave amplitude,t and x are the scaled time coordinate and scaled space coordinate,ρ1(t),ρ3(t),ρ5(t),ρ4(x,t) and ρ6(x,t) represent the variable coefficients of the nonlinear,dispersive,perturbed,dissipative and external-force terms,respectively [18,19].5References [18,19] have also,with references therein,listed out other applications of equation (4) in the dusty plasmas,interactionless plasmas,two-layer liquids,atmospheric flows,shallow seas and deep oceans.

By the bye,more nonlinear evolution equations might be found,for instance,in [24–28].

1.3.Our work and its difference from the existing literature

However,to our knowledge,for equation (1),there has been no B?cklund-transformation work with solitons reported as yet.No experimental comparison,either.

Our objective: In this paper,for equation (1),linking ρi's,we will make use of symbolic computation [29–32] to erect a B?cklund transformation,address some solitons and present the relevant experimental support.In addition,we will employ symbolic computation to construct a family of the similarity reductions.

2.Auto-B?cklund transformation for equation (1)

From the view of a generalized Laurent series,6Similar to those in [33–36].we introduce a Painlevé expansion,

to equation (1),around a noncharacteristic movable singular manifold given by an analytic function φ=0,where J represents a positive integer,while vj's are all the analytic functions with v0≠0 and φx≠0.Equilibrating the powers of φ at the lowest orders in equation (1) leads to J=2,and cutting expansion (5) at the constant-level terms yields

With symbolic computation,we next introduce expression (6) to equation (1),make the coefficients of like powers of φ disappear and present the Painlevé-B?cklund equations:

with the ‘˙’ sign hereby denoting the derivation of t and v2meaning a seed solution for equation (1) [34,35].

For an artery full of blood with an aneurysm,e.g., v(x,t)can be the dynamical radial displacement superimposed on the original static deformation from an arterial wall,and equations (7)–(12) with expression (6) formulate an auto-B?cklund transformation,on account of the mutual consistency within equations (6)–(12), or,the explicit solvability within equations (6)–(12) with regard to φ,v0and v1,to be seen below.By the bye,more B?cklund transformations could be found in [37–40].

3.Three solitonic families for equation (1),their difference and experimental support7For simplicity,derivations are elided.

3.1.Three explicitly-solvable solitonic families for equation (1)We now choose that

where β1and β4are the real constants while β2(t) and β3(t)are the real differentiable functions with β1≠0 since φx≠0.

Symbolic computation on auto-B?cklund transformation(6)–(12) along with expressions (13) results in three explicitly-solvable solitonic families for equation (1):

under the variable-coefficient constraints

which can be looked on as a quasi-solitary-wave8The word ‘quasi’ implies that there still exist ρ1(t),ρ2(t) and ρ3(t),beyond the travelling-wave format.case of solutions (14),if β2(t)=β5t+β6,β3(t)=β7,β4=0,and variable-coefficient constraints (15)–(18) are reduced to

where β5,β6and β7are the real constants;

Those solutions indicate that,e.g.,for an artery full of blood with an aneurysm,v(x,t) is the solitonic radial displacement superimposed on the original static deformation from an arterial wall.

3.2.Difference among those solitonic families

The difference among solitonic solutions (14),(19) and (24) is ascribable to the respective variable-coefficient constraints among ρ1(t),ρ2(t),ρ3(t),ρ4(x,t),ρ5(t) and ρ6(x,t), while ρi's are linked with the axial stretch of the injured artery,blood as an incompressible Newtonian fluid,radius variation in the axial direction or aneurysmal geometry,viscosity of the fluid,thickness of the artery,mass density of the membrane material,mass density of the fluid,strain energy density of the artery,shear modulus,stretch ratio,etc.

3.3.Experimental support

Especially,we call the attention that the shock-wave structures from solutions (24) have been shown to agree well with those dusty-plasma-experimentally reported,as detailed in[41,42] and references among the rest.Graphs describing the dynamical behaviors of solutions (24),versus those experimental graphs,have been worked out and presented in[41,42].

We need to say that such a dusty-plasma-experimental agreement directly supports the correctness/validation of auto-B?cklund transformation (6)–(12) and solitonic solutions(24),which in fact supports the correctness/validation of our above analytic work towards the blood vessel or circulatory system.

By the way,other relevant solitonic issues might be found in [43–53].

4.Similarity reductions for equation (1)

Our assumption

which is similar to those in [54–59],can lead to certain similarity reductions for equation (1),with θ(x,t),κ(x,t) ≠0 and r(x,t) ≠0 as the real differentiable functions to be determined.

Making use of symbolic computation and inserting assumption (25) into equation (1) turn into

in which the ‘′’ sign hereby stands for the derivation with respect to r.

Taking into consideration that equation (26) can be designed to reduce to a single ordinary differential equation(ODE) as for q(r),one requires those ratios of derivatives and/or powers of q(r) to mean some functions with respect to r only,so that

with Γχ(r)'s (χ=1,…,6) as merely the real functions of r,of course to be determined.

Grounded on the 2nd freedom of remark 3 in [59],equation (27a) gives rise to

Because of the 1st freedom of remark 3 in [59],equation (27b) brings about

and then equation (27c) results in with ξ1indicating a real non-zero constant,while λ1(t) and λ2(t) implying two real non-zero differentiable functions with respect to t.

Because the 1st freedom of remark 3 in [59] helps us reduce equation (27d) to

equation (27e) turns to

and then equation (27f) develops into

with μ1and μ2as two real non-zero constants,while μ3and μ4as two real constants.

For an artery full of blood with an aneurysm,e.g., v(x,t)can be the dynamical radial displacement superimposed on the original static deformation from an arterial wall,and in general,under the variable-coefficient constraints

we build up the following family of the similarity reductions for equation (1):

ODE (35c) has been presented in [60] and thus can be considered as a known ODE.

5.Conclusions

With respect to the nonlinear waves in an artery full of blood with a certain aneurysm,pulses in a blood vessel,or features in a circulatory system,this paper has symbolically computed out auto-B?cklund transformation (6)–(12) and solitonic solutions (14),(19) and (24) for equation (1),i.e.,a variablecoefficient generalized forced–perturbed Korteweg–de Vries–Burgers equation.We have also built up similarity reductions (35), from equation (1) to a known ODE.Aiming,e.g.,at the dynamical radial displacement superimposed on the original static deformation from an arterial wall,our results rely on the axial stretch of the injured artery,blood as an incompressible Newtonian fluid,radius variation along the axial direction or aneurysmal geometry,viscosity of the fluid,thickness of the artery,mass density of the membrane material,mass density of the fluid,strain energy density of the artery,shear modulus,stretch ratio,etc.Relevant variablecoefficient constraints have also been given.Finally,we have highlighted that the shock-wave structures from our solutions agree well with those dusty-plasma-experimentally reported.

Acknowledgments

We express our sincere thanks to the Editors and Reviewers for their valuable comments.This work has been supported by the National Natural Science Foundation of China under Grant Nos.11871116 and 11772017,and by the Fundamental Research Funds for the Central Universities of China under Grant No.2019XD-A11.