Journal of Modern Physics, 2014, 5, 394-401 Published Online April 2014 in SciRes. http://www.scirp.org/journal/jmp http://dx.doi.org/10.4236/jmp.2014.56051 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). http://creativecommons.org/licenses/by/4.0/ Abstract 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. Keywords 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. http://dx.doi.org/10.4236/jmp.2014.56051 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 last. L L r1 L r1 Vn Vn+1 In-1 r1 In+1 C (Un) Un C (Un+1) r2 r2 Figure 1. Schematic representation of the nonlinear electric line. 395 Un+1 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 , (1) 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 LC1 ( ) ( ) d2 d d 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 dt dt 2 d U n2−1 + U n2+1 − 2U n2 + r2C1 dt ( LC1 , we express Equation (2) in the following dimensionless form: Next, measuring time in units of ( ) ) (2) ( ) r d r d d2 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 ) 2 Lu0 dt Lu0 dt dt ( r d U n2−1 + U n2+1 − 2U n2 + 2 Lu0 dt ) (3) 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: 2 r1 ∂ r2 ∂ 2 ∂ 2U δ 4 ∂ 4U δ 4 ∂ 4U ∂2 2 2 2 ∂ U U U U U α α δ + + + = + + + δ Lu0 ∂t ∂t 2 ∂x 2 12 ∂x 4 Lu0 ∂t ∂x 2 12 ∂x 4 2 2 2 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 ( ) ( ) (4) Setting dimensionless voltage F, such that F = 2αU , after some extensive manipulations [26], we obtain the equation describing the developing solitonic pulses: ∂W ∂W ∂ 3W ∂ 2W µ3W + µ2 2 , − 6W + µ1 3 = ∂τ ∂ξ ∂ξ ∂ξ (5) 12r1 12r2δ 2 δ2 = µ = µ , , To derive the previous equation, we used within the perturbative me2 1 L µ02 L µ02 µ02 thod [27] the Gardner-Morikawa transformation that guarantees the balance between the nonlinearity and dispersion in series expansion as follows: = µ3 with, 1 3 ξ= ε 3 ( x − u0 t ) , τ = ε3 and a further transformation of constants and variables such as 396 u0 t, 24 (6) S. D. Yamigno 3 1 r1 → ε 3 r1 , r2 → ε 3 r2 , F → W 2 (7) 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τ . (8) 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 Ω (τ ) (ξ − η (τ ) ) , W= (9) is obtained, where the quantity Ω (τ ) and η (τ ) obey the equations dΩ 1 +∞ 2 = − ∫ dzP (W ) sech z, dτ 4Ω −∞ (10) dη 1 +∞ 1 2 4Ω 2 − = ∫ dzP (W ) z + 2 sinh ( 2 z ) sech z, dτ 4Ω3 −∞ (11) = z Ω (τ ) (ξ − η (τ ) ) and P (= where U ) µ3W + µ2 Ω2 = ∂ 2W ∂ξ 2 Ω 02 exp ( −τ τ 0 ) τ 4r 1 + 2 δ 2 Ω 02 1 − exp − 5 r1 τ0 , (12) 1 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. 397 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 398 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 Ω. (a) (b) 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 Ω. 2 (a) (b) 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 Ω. 2 399 S. D. Yamigno (a) (b) 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 Ω. 2 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. Acknowledgements 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. References [1] Remoissenet, M. (1999) Waves Called Solitons. 3rd Edition, Springer-Verlag, Berlin. http://dx.doi.org/10.1007/978-3-662-03790-4 [2] Hirota, R. and Suzuki, K. (1970) Journal of the Physical Society of Japan, 28, 1366-1367. http://dx.doi.org/10.1143/JPSJ.28.1366 [3] Lonngren, K.E. and Scott, A.C. (1978) Solitons in Action. Academic, New York. [4] 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 400 [5] 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 [6] Kivshar, Y.S. and Agrawal, G.P. (2013) Optical Solitons: From Fibers to Photonic Crystals. Academic, New York. [7] Sekulić, D.L., Satarić, B.M., Tuszynski, J.A. and Satarić, M.V. (2011) European Physical E: Soft Matter, 34, 49. http://dx.doi.org/10.1140/epje/i2011-11049-0 [8] Sekulić, D.L. and Sataric, M.V. (2012) 28th International Conference on Microelectronics-MIEL. in press. [9] Marquie, P., Bilbault, J.M. and Remoissenet, M. (1994) Physical Review E, 49, 828-835. http://dx.doi.org/10.1103/PhysRevE.49.828 [10] Marquie, P., Bilbault, J.M. and Remoissenet, M. (1995) Physical Review E, 51, 6127-6133. http://dx.doi.org/10.1103/PhysRevE.51.6127 [11] Yemélé, D., Marquie, P. and Bilbault, J.M. (2003) Physical Review E, 68, 016605-016616. http://dx.doi.org/10.1103/PhysRevE.68.016605 [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. http://dx.doi.org/10.1016/j.amc.2010.04.026 [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. http://dx.doi.org/10.1016/0375-9601(96)00283-6 [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. http://dx.doi.org/10.1016/j.amc.2007.10.059 [19] Wang, Q., Chen, Y. and Zhan, H.Q. (2005) Chaos, Solitons and Fractals, 25, 1019-1028. http://dx.doi.org/10.1016/j.chaos.2005.01.039 [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. http://dx.doi.org/10.1016/S0375-9601(01)00580-1 [22] Zhou, Y.B., Wang, M.L. and Wang, Y.M. (2003) Physics Letters A, 308, 31-36. http://dx.doi.org/10.1016/S0375-9601(02)01775-9 [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. http://mipro-proceedings.com/sites/mipro-proceedings.com/files/upload/meet/meet_046.pdf [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. http://dx.doi.org/10.1088/0031-8949/81/03/035801 [28] Kivshar, Y.S. and Luther-Davies, B. (1998) Physics Reports, 298, 81-197. http://dx.doi.org/10.1016/S0370-1573(97)00073-2 [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. http://dx.doi.org/10.1103/PhysRevA.77.023613 [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. http://dx.doi.org/10.1103/PhysRevLett.60.29 [35] Kivshar, Y.S., Haelterman, M., Emplit, P. and Hamaide, J.P. (1994) Optics Letters, 19, 19-21. http://dx.doi.org/10.1364/OL.19.000019 [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 401

© Copyright 2017 ExploreDoc