Full-Text PDF - Scientific Research Publishing

Journal of Modern Physics, 2014, 5, 394-401
Published Online April 2014 in SciRes. http://www.scirp.org/journal/jmp
Propagation of Dark Solitary Waves in the
Korteveg-Devries-Burgers Equation
Describing the Nonlinear RLC Transmission
Serge Doka Yamigno
Department of Physics, Higher Teachers’ Training College, University of Maroua, Maroua, Cameroon
Email: [email protected]
Received 2 January 2014; revised 1 February 2014; accepted 18 February 2014
Copyright © 2014 by author and Scientific Research Publishing Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
We investigate the propagation of dark solitons in a nonlinear dissipative electrical line. We show
that the dynamics of the line is reduced to an expanded Korteweg-de Vries-Burgers (KdVB) equation. By applying the perturbation theory to the KdVB equation, we obtain soliton-like pulse solutions. The numerical simulations of the discrete equation are carried out and show the possibility
of the founding solution to spread through the line. The effect of the dissipation through soliton is
also shown. A chaotic-like behavior can take place in the system during the propagation of dark
solitons through the line.
Nonlinear Electrical Lines, Dark Solitons
1. Introduction
In recent years, nonlinear discrete system has received considerable attention. There are two main reasons for
this interest: the development of experimental techniques making it possible to realize experiment in complicated periodic nonlinear structure and the potential for all-optical switching applications. The nonlinear electrical lines are very convenient tools for studying quantitatively the fascinating properties of wave propagation in
nonlinear dispersive media [1] [2]. In particular, they provide a useful way to check directly how the nonlinear
excitations behave inside the nonlinear medium by means of probes related to an oscilloscope. The first studies
of soliton on electrical lattices have been done by Hirota and Suzuki [2]. Qualitatively, the origin of soliton in
nonlinear electrical line is explained by the balance between the effect of dispersion (due to the periodic location
of capacitor in the nonlinear electrical lines) and nonlinearity (due to the voltage dependence of the capacitance).
How to cite this paper: Yamigno, S.D. (2014) Propagation of Dark Solitary Waves in the Korteveg-Devries-Burgers Equation
Describing the Nonlinear RLC Transmission. Journal of Modern Physics, 5, 394-401.
S. D. Yamigno
The nonlinear transmission lines (NLTLs) are discrete systems but approximate the continuum system quite well,
as will be seen below [1] [3]. The NLTLs provide a useful way to check how the nonlinear excitations behave
inside the nonlinear medium and to model the exotic properties of new systems [3]. Depending on device technology and design, NLTLs can easily be fabricated and integrated with resonant tunneling diodes (RTDs) using
InP technology [4] [5]. Let us also point out that, recently, NLTLs have proven to be of great practical use in extremely wideband (frequencies from dc to 100 GHz) focusing and shaping of signals [6] which is usually a hard
problem. Various investigators have discovered the existence of solitons in NLTLs, through both mathematical
models and physical experiments. Interestingly, in constructing proposed nonlinear RLC network, we were motivated with a very interesting example of nanobio electronics problem. Namely, Sekulić et al. [7] have recently
modeled ionic currents along microtubules in the eukaryotic cells, they have analytically analyzed the possible
use of microtubules as protein structure for building biomolecular nanoscale NLTLs in the context of the polyelectrolyte character of these cytoskeletal filaments [8]. Solitary wave propagation is one of the important phenomena associated with the NLTL. It has been analyzed theoretically that there could exist solitary waves with
both bright and dark soliton solutions in the NLTL with nonlinear element replacing the loaded shunt inductors
or the shunt capacitors [9]-[11]. Many soliton equations describing nonlinear systems are known and solitons
themselves have been observed (directly or indirectly) in various media. However, there are only a few systems
where solitons are easily and directly observed in controlled laboratory experiments. So, finding exact solutions
is important to understand the mechanism of the complicated nonlinear physical phenomena. In the recent decade, several methods for finding the exact solutions to nonlinear equations of mathematical physics have been
proposed, such as trigonometric function series method [12], the modified mapping method and the extended
mapping method [13] [14], homogeneous balance method [15], tanh function method [16], extended tanh function method [17], hyperbolic function method [18], rational expansion method [19], sine-cosine method [20],
Jacobi elliptic function method [21], F-expansion method [22], and so on. Our goal is to present a method to
computational study of solitary waves in nonlinear RLC transmission lines. A circuit model of a NLTL with N
identical cells is given in Figure 1. The series inductance is due to magnetic field effects, and the capacitance is
due to electric field coupling between the lines. The losses in the transmission media are depicted by the series
and the shunt resistors. These resistors represent the finite conductivity of the conductors and the dielectric insulator between the conductors, respectively. The resistors r1 and r2, accounting for the transmission line losses,
the linear inductance L, and the voltage-dependent capacitance C(Vn) are the circuit parameters. The resulting
circuit is referred to as a distributed model of a nonlinear RLC transmission lines.
In this work, after writing the discrete equation of the NLTL trough the Kirchhoff’s laws, we apply the perturbative method to reduce the previous equation to the well-known KdV-Burgers equation. This equation has
solitary wave solution namely the dark soliton. Then we numerically investigated the propagation of this dark
soliton through the line. Numerical investigations show also the impact of the dissipation on the soliton and the
possible chaotic behavior. The remainder of the paper is organized as follows: in Section 2 we present the model
of the NLTL understudy and we derive the equation governing modulated waves in the network. In Section 3,
numerical simulations confirm the propagation of solitary wave through the line. This paper is concluded in the
C (Un)
C (Un+1)
Figure 1. Schematic representation of the nonlinear electric line.
S. D. Yamigno
2. Model and Its Analytical Consideration
We consider a nonlinear transmission line which consists of a number of RLC blocks connected as illustrated in
Figure 1. The capacitance-voltage relationship can be approximated by a second-order polynomial [1] [7] [9]
[14] [17] [23]-[25]
C (U 0 + U n ) = C1 1 + 2αU n + 3β U n2 ,
where C1 is a constant corresponding to the capacitance of the nonlinear diode at the dcbias-voltage U0. The
nonlinear parameters α and β are assumed to be positive or negative constants. α and β are nonlinear
coefficients that determine electric charge Qn stored in the n-th capacitor in the line, and Un is the voltage across
the n-th capacitor. For the diode the sign of α is negative and the sign of β is positive and vice versa for the
MOS. By applying the Kirchhoff’s laws to transmission line we obtain
U n − αU n2 + r1C1
U n − αU n2 = U n −1 + U n +1 − 2U n + r2C1 (U n −1 + U n +1 − 2U n )
dt 2
U n2−1 + U n2+1 − 2U n2
+ r2C1
LC1 , we express Equation (2) in the following dimensionless form:
Next, measuring time in units of
r d
r d
U n − αU n2 + 1
U n − αU n2 = U n −1 + U n +1 − 2U n + 2
(U n−1 + U n+1 − 2U n )
Lu0 dt
Lu0 dt
r d
U n2−1 + U n2+1 − 2U n2
+ 2
Lu0 dt
Starting from the semi-discrete model given by Equation (3), we develop a continuum model in the standard
way, i.e. the right-hand side of (3) can be approximated with partial derivatives with respect to distance x = nδ ,
from the beginning of the line, assuming that the spacing between two adjacent sections is δ (This is spatial
sampling period equals the length of each element). Therefore, an approximate continuous partial differential
Equation (3) can be obtained by using the Taylor expansions of Un + 1(t) and Un − 1(t), both around the reference
value Un(t). We regard the following equation as a continuum model of the transmission line that retains the effect of discreteness in the fourth-order term:
r1 ∂
r2 ∂  2 ∂ 2U δ 4 ∂ 4U 
δ 4 ∂ 4U
2 ∂ U
Lu0 ∂t
∂t 2
∂x 2 12 ∂x 4 Lu0 ∂t 
∂x 2 12 ∂x 4 
r2 ∂  2  ∂U 
1 ∂ 4U 1 ∂U ∂ 3U  
4 1∂ U 
2δ 
 + δ   2  + U 4 +
Lu0 ∂t 
2 ∂x  6 ∂x
3 ∂x ∂x3  
 ∂x 
 
Setting dimensionless voltage F, such that F = 2αU , after some extensive manipulations [26], we obtain the
equation describing the developing solitonic pulses:
∂ 3W
∂ 2W
µ3W + µ2 2 ,
− 6W
+ µ1 3 =
12r2δ 2
To derive the previous equation, we used within the perturbative me2
L µ02
L µ02
thod [27] the Gardner-Morikawa transformation that guarantees the balance between the nonlinearity and dispersion in series expansion as follows:
ε 3 ( x − u0 t ) , τ =
and a further transformation of constants and variables such as
S. D. Yamigno
r1 → ε 3 r1 , r2 → ε 3 r2 , F →
where ε is an indicator of the slowness. Equation (5) governs the propagation of traveling on linear waves in
the network. This equation is well known as the perturbed KdV equation. In particular, in the absence of dissipative terms (r1 = r2 = 0) it can support the well-known pulse soliton solution of the KdV (Korteweg-de Vries)
equation. The KdV equation is the generic model for the study of weakly nonlinear long waves. It arises in
physical systems which involve a balance between nonlinearity and dispersion at leading-order. For example, it
describes surface waves of long wave length and small amplitude on shallow water and internal waves in a
shallow density-stratified fluid. Many other applications for the KdV equations so exist, such as plasma waves,
Ross waves and magma flow. Also, the KdV equation is integrable. This means that a collision between KdV
solitary waves is elastic; after the collision the solitons retain their original shape with the only memory of the
collision being a phase shift. The explicit solution for interacting KdV solitons was developed using the bilinear
transformation method, by Hirota [2]. Which has a soliton solution explicitly described as
W =−2Ω 2 sech 2 Ω ξ − 4Ω 2τ  .
In the case of finite r1 and r2 the exact solution of Equation (8) can no longer be obtained. However, in the
framework of the perturbation theory as applied to KdV equation on the basis of the inverse scattering method
[28], a soliton-like pulse
−2Ω 2 (τ ) sech 2 Ω (τ ) (ξ − η (τ ) )  ,
is obtained, where the quantity Ω (τ ) and η (τ ) obey the equations
1 +∞
= −
∫ dzP (W ) sech z,
4Ω −∞
1 +∞
4Ω 2 −
∫ dzP (W )  z + 2 sinh ( 2 z )  sech z,
4Ω3 −∞
z Ω (τ ) (ξ − η (τ ) ) and P (=
U ) µ3W + µ2
Ω2 =
∂ 2W
∂ξ 2
Ω 02 exp ( −τ τ 0 )
 τ 
1 + 2 δ 2 Ω 02 1 − exp  −  
5 r1
 τ0 
and Ω 0 corresponds to the unperturbed soliton’s amplitude. This shows that Ω decreas16r1C1
es as τ increases.
Another particular case, when we set r1 = 0 and r1 ≠ 0 Equation (8) is commonly called the Korteweg-de
Vries-Burgers equation. This equation plays an important role in applications. In certain cases, the coefficient
µ2 in equation (8) may be negative. In those cases, the perturbed KdV equation describes the development of
instability. For r1 = 0 and r1 ≠ 0 , perturbation describes a variable depth in shallow water theory; µ3 corresponds to a constant depth gradient. The KdV equation with µ3 = 0 finds a similar application in the description of nonlinear ion acoustic waves in an inhomogeneous plasma [28].
where τ 0 =
3. Numerical Investigations
In order to study the behavior of the soliton along the cascaded maps, we numerically solve Equation (3) using
the fourth-order Runge-Kutta scheme with a normalized integration time step ∇t =10−1 . The parameters of the
line are: L = 470 μH, C1 = 320 pF. At the input of the line, we launch as initial solution the dark solution obtained above. Figure 2, Figure 3 and Figure 4 show the propagation of these solutions through the line.
S. D. Yamigno
Figure 2. The pulse decay of an initial pulse V(n = 0) = ‒sech2(n − n0) at n0 =
25 cell and different times, with α = 0.21 V−1, r1 = 1 Ω; r2 = 0.1 Ω.
Figure 3. The pulse decay of an initial pulse V(n = 0) = −sech2(n−n0) at n0 =
25 cell and different times, with α = 0.21 V−1, r1 = 0.1 Ω; r2 = 1 Ω.
It can be seen that the soliton keep its shape when traveling down the line. Dark soliton appears as an intensity
dip in an infinitely extended constant background. It has been investigated in many theoretical and experimental
papers [29] [30]. The dark solitons have attracted much attention due to their potential applications [31]. For
example, using structures created during the propagation and interaction of the dark solitons, many types of
all-optical switches may be written [32]. Other applications of the dark solitons involve the optical logic devices
[33] and waveguide optics as dynamic switches and junctions [34]. They are also considered in signal processing
and communication because of their inherent stability [35]. Compared with the bright solitons, they have better
stability against various perturbations such as fiber loss, mutual interaction between neighboring pulses, the
Raman effect, and the superposition of noise emitted from optical amplifiers. These solutions are plotted for cell
25. Coherent structures and chaotic states are well known as two distinct states of nonlinear dissipative wave
S. D. Yamigno
Figure 4. The pulse decay of an initial pulse V(n = 0) = −sech2(n − n0) at n0 =
25 cell and different times, with α = 0.21 V−1, r1 = 0 Ω; r2 = 0.0 Ω.
Figure 5. The evolution of dark soliton solution of an initial pulse V(n = 0) = −sech (n − n0) at n = 25 cell with α = 0.21
V−1, r1 = 1Ω; r2 = 0.1 Ω.
Figure 6. The evolution of dark soliton solution of an initial pulse V(n = 0) = −sech (n − n0) at n = 25 cell with α = 0.21
V−1, r1 = 0.1 Ω; r2 = 1 Ω.
S. D. Yamigno
Figure 7. The evolution of dark soliton solution of an initial pulse V(n = 0) = −sech (n − n0) at n = 25 cell with α = 0.21
V−1, r1 = 0 Ω; r2 = 0 Ω.
systems. However, these states sometimes occur and propagate together in some systems. Figure 5, Figure 6,
Figure 7 depict the incoherent evolution of dark solitons through the line. This incoherent evolution of dark soliton can also be evidenced from the nonreproducibility of experiments devoted to their propagation in the nonlinear medium, as observed by Ablowitz et al. [36] in the context of fluid dynamics that is, considering modulated periodic Stokes waves in deep water. For two different experiments with initial identical signals generated by the wave maker, the resulting temporal evolutions of the surface displacement at a given position in the
tank are graphed against each other to produce a “phase plane” plot indicating the level of reproducibility.
4. Conclusion
In this paper, we have studied analytically and numerically dark soliton. These solutions have been also obtained
in microtubules protein structure. We have used the nonlinear transmission line analog of this structure to derive
a nonlinear lattice equation governing the voltage equation of the system. This lattice equation was then treated
analytically, by means of the continuum medium approximation. The latter is a variant of the perturbation method, which takes into regard the discreteness of the system by considering the carrier (dark) wave as a discrete
(continuum) object. This approach allowed us to derive, in the wavelength approximation, a nonlinear perturbed
KdV equation. The KdV model was then used to predict the formation of dark soliton. Numerical simulations
performed in the framework of the nonlinear lattice equation, with initial conditions borrowed from the exact
KdV solution, revealed that the dark soliton propagated through the line. Furthermore, it was shown how certain
physical parameters of the configuration may affect the results. We numerically found soliton profile and characteristics (amplitude and width) that are in very good agreement with the analytical predictions. Our study also
shows the presence of chaos in the system.
S.Y.D. is very grateful for the support of the Higher Teachers’ Training College of Maroua through his Director,
and to A. M and A. S.
Remoissenet, M. (1999) Waves Called Solitons. 3rd Edition, Springer-Verlag, Berlin.
Hirota, R. and Suzuki, K. (1970) Journal of the Physical Society of Japan, 28, 1366-1367.
Lonngren, K.E. and Scott, A.C. (1978) Solitons in Action. Academic, New York.
Reddy, M., Martin, S.C., Molnar, A.C., Muller, R.E., Smith, R.P., Siegel, P.H., Mondry, M.J., Rodwell, M.J.W.,
Kroemer, H. and Allen Jr., S.J. (1997) IEEE Electron Device Letters, 18, 218-221. http://dx.doi.org/10.1109/55.568771
S. D. Yamigno
Jaeger, I., Kalinowski, D., Stöhr, A., Stiens, J., Vounckx, R. and Jäger, D. (2007) Microwave and Optical Technology
Letters, 49, 2907-2909. http://dx.doi.org/10.1002/mop.22902
Kivshar, Y.S. and Agrawal, G.P. (2013) Optical Solitons: From Fibers to Photonic Crystals. Academic, New York.
Sekulić, D.L., Satarić, B.M., Tuszynski, J.A. and Satarić, M.V. (2011) European Physical E: Soft Matter, 34, 49.
Sekulić, D.L. and Sataric, M.V. (2012) 28th International Conference on Microelectronics-MIEL. in press.
Marquie, P., Bilbault, J.M. and Remoissenet, M. (1994) Physical Review E, 49, 828-835.
[10] Marquie, P., Bilbault, J.M. and Remoissenet, M. (1995) Physical Review E, 51, 6127-6133.
[11] Yemélé, D., Marquie, P. and Bilbault, J.M. (2003) Physical Review E, 68, 016605-016616.
[12] Zhang, Z.Y. (2013) Turkish Journal of Physics, 37, 259-267. http://dx.doi.org/10.3906/fiz-1205-13
[13] Zhang, Z.Y., Liu, Z.H., Miao, X.J. and Chen, Y.Z. (2010) Applied Mathematics and Computation, 216, 3064-3072.
[14] Abdoulkary, S., Ousmanou, Mohamadou, A. and Beda, T. (article submitted in CPB 2012).
[15] Wang, M.L., Zhou, Y.B. and Li, Z.B. (1996) Physics Letters A, 216, 67-75.
[16] Yan, Z.Y. and Zhang, H. (2001) Physics Letters A, 285, 355-362. http://dx.doi.org/10.1016/S0375-9601(01)00376-0
[17] Abdoulkary, S., Mohamadou, A. and Beda, T. (2011) Communications in Nonlinear Science and Numerical Simulation,
16, 3525-3532. http://dx.doi.org/10.1016/j.cnsns.2010.12.029
[18] Shang, Y., Huang, H. and Yuan, W.J. (2008) Applied Mathematics and Computation, 200, 110-122.
[19] Wang, Q., Chen, Y. and Zhan, H.Q. (2005) Chaos, Solitons and Fractals, 25, 1019-1028.
[20] Yan, C. (1996) Physics Letters A, 224, 77-84. http://dx.doi.org/10.1016/S0375-9601(96)00770-0
[21] Liu, S.K., Fu, Z.T., Liu, S.D. and Zhao, Q. (2001) Physics Letters A, 289, 69-74.
[22] Zhou, Y.B., Wang, M.L. and Wang, Y.M. (2003) Physics Letters A, 308, 31-36.
[23] Sekulic, D.L., Sataric, M.V., Zivanov, M.B. and Bajic, J.S. (2012) EEE-Kaunas, Technologija, 121, 53-58.
[24] Sekulic, D.L. and Zivanov, M.B. (2012) MIPRO.
[25] Peyrard, M. and Dauxois, T. (2004) Physique des Solitons. EDP Sciences/CNRS Editions, Les Ulis.
[26] Kivshar, Y.S. and Malomed, B.A. (1989) Reviews of Modern Physics, 61, 763-915.
[27] Essimbi, B.Z. and Jager, D. (2010) Physica Scripta, 81, Article ID: 035801.
[28] Kivshar, Y.S. and Luther-Davies, B. (1998) Physics Reports, 298, 81-197.
[29] Weiner, A.M. (1992) Optical Solitons: Theory and Experiment. Cambridge University, Cambridge.
[30] Zhang, X.F., Yang, Q., Zhang, J.F., Chen, X.Z. and Liu, W.M. (2008) Physical Review A, 77, 023613-023621.
[31] Luther-Davies, B. and Yang, X.P. (1992) Optics Letters, 17, 1755-1757. http://dx.doi.org/10.1364/OL.17.001755
[32] Swartzlander Jr., G.A. (1992) Optics Letters, 17, 493-495. http://dx.doi.org/10.1364/OL.17.000493
[33] Luther-Davies, B. and Yang, X.P. (1992) Optics Letters, 17, 496-498. http://dx.doi.org/10.1364/OL.17.000496
[34] Krökel, D., Halas, N.J., Giuliani, G. and Grischkowsky, D. (1988) Physical Review Letters, 60, 29-32.
[35] Kivshar, Y.S., Haelterman, M., Emplit, P. and Hamaide, J.P. (1994) Optics Letters, 19, 19-21.
[36] Ablowitz, M.J., Hammack, J., Henderson, D. and Schober, C.M. (2001) Physica D: Nonlinear Phenomena, 152-153,
416-433. http://dx.doi.org/10.1016/S0167-2789(01)00183-X