Research Article  Open Access
Approximation for Transient of Nonlinear Circuits Using RHPM and BPES Methods
Abstract
The microelectronics area constantly demands better and improved circuit simulation tools. Therefore, in this paper, rational homotopy perturbation method and Boubaker Polynomials Expansion Scheme are applied to a differential equation from a nonlinear circuit. Comparing the results obtained by both techniques revealed that they are effective and convenient.
1. Introduction
Industrial competition constantly pushes the area of electronic circuit design to the limits of technology. This has caused a rapid growth in the levels of integration for integrated circuits and the emergence of novel devices such as singleelectron transistors and memristors. Because of this, the development and improvement of mathematical and numerical tools, applied to circuit simulation for the transient domain, are important. In the dynamic domain (transient), the circuit analysis is carried out only numerically because the resulting differential equations are highly nonlinear. Nevertheless, several methods are focused to find approximate solutions to nonlinear differential equations like homotopy perturbation method (HPM) [1–12], rational homotopy perturbation method (RHPM) [5, 6], variational iteration method (VIM) [13–16], and Boubaker Polynomials Expansion Scheme (BPES) [17–36], among many others. Therefore, we propose the comparison between RHPM and BPES methods by solving the nonlinear differential equation that represents the dynamics of a nonlinear circuit. The results should be a meaningful supply for monitoring complex nonlinear circuits behaviours and responses. In fact, the used protocols try to embed boundary conditions instead of direct solving, as preceded in spectral or limitcycle bifurcations approaches.
This paper is arranged as follows. In Section 2, we present the differential equation of a nonlinear circuit. Sections 3 and 4 present the fundamentals of RHPM and BPES methods, respectively. The solutions obtained using both methods are explained in Section 5. Comparisons between the two methods and some other results presented in the recent literature have been illustrated in Section 6. Conclusions will be discussed in Section 7.
2. Nonlinear Circuit
The rapid increase in the number of transistors by integrated circuit and the increase of complexity for the models (a result of lowering the dimension of the components) results in a complex calculation for the transient. Furthermore, the task of tracing the transient for nonlinear circuits is a critical and difficult task. In fact, commercial circuit simulators do not provide any symbolic/analytic solution for the transient of any given circuit. Instead, the simulator provides only numerical data that allows circuit designers to explore a limited range of dynamics of nonlinear circuits.
Consider the analysis of the nonlinear circuit depicted in Figure 1 as a case study [37]. Let the branch relationship of the nonlinear capacitor to be where is the voltage, is the charge of the capacitor, and is a parameter of the capacitor.
By applying Kirchhoff Laws, we obtain the equation for the transient Therefore, If we consider the case for DC excitation, then, resulting in
3. Fundamentals of the Rational Homotopy Perturbation Method
The rational homotopy perturbation method RHPM [5, 6] can be considered as a combination of the classical perturbation technique [38, 39] and the homotopy (whose origin isin the topology) [40–42] but not restricted to a small parameter like traditional perturbation methods. For example, RHPM requires neither small parameter nor linearization, but only few iterations to obtain accurate solutions.
To figure out how RHPM method works, consider a general nonlinear equation in the form: with the following boundary conditions: where is a general differential operator, is a boundary operator, is a known analytical function, and is the domain boundary for. can be divided into two operators and , where is linear and nonlinear; from this last statement, (5) can be rewritten as Generally, a homotopy can be constructed in the form [1–3]: where is a homotopy parameter, whose values are within the range of 0 and is the first approximation for the solution of (6) that satisfies the boundary conditions.
When , (8) is reduced to where operator possesses trivial solution.
For, (8) is reduced to the original problem Assuming that the solution for (8) can be written as a power series of : where are unknown functions to be determined by the RHPM, and are known analytic functions of the independent variable.
Substituting (11) into (8) and equating identical powers of terms, it is possible to obtain values for the sequence .
Whenin (11), it yields in the approximate solution for (5) in the form: Convergence of RHPM method is studied in [5, 6].
4. Fundamentals of the Boubaker Polynomials Expansion Scheme BPES
The Boubaker Polynomials Expansion Scheme BPES [17–36] is a resolution protocol, which has been successfully applied to several appliedphysics and mathematical problems. The BPES protocol ensures the validity of the related boundary conditions regardless of main equation features. The Boubaker Polynomials Expansion Scheme BPES is based on the Boubaker polynomials first derivatives properties: with Several solutions have been proposed through the BPES in many fields like numerical analysis [17–20], theoretical physics [21–24], mathematical algorithms [25], heat transfer [26], homodynamic [27, 28], material characterization [29], fuzzy systems modelling [30–34], and biology [35, 36].
5. Application of RHPM and BPES
5.1. Solution Using RHPM Method
Using (8), we establish the following RHPM homotopy map: where the trial function.
Using (11), we propose the following rational solution: where and .
We substitute (16) into (15), regroup, and equate terms with identical powers of . In order to fulfil boundary condition of (16), it follows that for the homotopy map.
The results are recast in the following systems of differential equations:
Solving (17) yields Substituting (18) into (16) and calculating the limit when, we obtain the seventhorder approximation: If we consider , , and as reported in [37], it is possible to obtain the adjustment parameters using the procedure reported in [6, 12], resulting in and .
5.2. Solution Using the Boubaker Polynomials Expansion Scheme BPES
The Boubaker Polynomials Expansion Scheme BPES is applied to (4) using the setting expression: Using the properties provided by (13), boundary conditions are verified in advance of the resolution process. The system in (16) is reduced to Boundary conditions become redundant since they are already verified by the proposed expansion, consecutively, and thus, majoring and integrating along the given interval for the time variable transform the problem in a linear system with unknown real variables: . Calculations are reduced to approximately arithmetical operations. Solutions are obtained by using the Householder [39, 40] algorithm detailed elsewhere and are denoted by. The final solution is given as
6. Results and Discussion
From Table 1, we can observe that RHPM solution (19) and BPES solution (22) are in good agreement with the numerical results obtained using Fehlberg fourthfifthorder RungeKutta method with degree four interpolant (RKF45) [43, 44] (builtin function of Maple software). In order to guarantee a good numerical reference, RKF45 is configured using an absolute error of and a relative error of. The power of RHPM method is based on the capability of rational expressions containing a huge amount of information of dynamics from asymptotic problems.

Moreover, convergence of the BPES algorithm has been obtained for moderate values of (), since, as mentioned above, boundary conditions were verified in advance of the resolution process. Both methods generated analytical expressions useful for other analysis like circuit power consumption; such analytical expressions can provide more information about the nature and behaviour of circuits than numerical integration schemes with variable step size [41, 44–46]. Nonetheless, semianalytical techniques like RHPM and BPES may be combined with numerical methods [43–46] to improve the simulation tools of VLSI circuits.
7. Conclusion
In this paper, powerful analytical methods like rational homotopy perturbation method (RHPM) and Boubaker Polynomials Expansion Scheme (BPES) are presented to construct semianalytical solutions for the transient of a nonlinear circuit. The results exhibited that both techniques are powerful, obtaining highly accurate analytical expressions for the transient of a simple test circuit. While RHPM yielded accurate and reliable results, BPES exhibited the advantage of ensuring the validity of boundary conditions regardless of main equation features. This feature made the protocol yielding faster and provided more convergent solutions than many numerical integration schemes with variable step size. Further work is necessary to extend the use of both methods for larger circuits.
Acknowledgments
The authors gratefully acknowledge the financial support provided by the National Council for Science and Technology of Mexico (CONACyT) through Grant CB201001 no. 157024. The authors would like to express their gratitude to RogelioAlejandro CallejasMolina and Roberto RuizGomez for their contribution to this project.
References
 J. H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 34, pp. 257–262, 1999. View at: Publisher Site  Google Scholar
 J. H. He, “A coupling method of a homotopy technique and a perturbation technique for nonlinear problems,” International Journal of NonLinear Mechanics, vol. 35, no. 1, pp. 37–43, 2000. View at: Publisher Site  Google Scholar
 J. H. He, “Homotopy perturbation method: a new nonlinear analytical technique,” Applied Mathematics and Computation, vol. 135, no. 1, pp. 73–79, 2003. View at: Publisher Site  Google Scholar
 H. VazquezLeal, Y. Khan, G. FernandezAnaya et al., “A general solution for Troesch's problem,” Mathematical Problems in Engineering, vol. 2012, Article ID 208375, 14 pages, 2012. View at: Publisher Site  Google Scholar
 H. VazquezLeal, A. SarmientoReyes, Y. Khan, U. FilobelloNiño, and A. DiazSanchez, “Rational biparameter homotopy perturbation method and laplacepadé coupled version,” Journal of Applied Mathematics, vol. 2012, Article ID 923975, 21 pages, 2012. View at: Publisher Site  Google Scholar
 H. VázquezLeal, “Rational homotopy perturbation method,” Journal of Applied Mathematics, vol. 2012, Article ID 490342, 14 pages, 2012. View at: Publisher Site  Google Scholar
 Y. Khan, H. VázquezLeal, and L. HernandezMartínez, “Removal of noise oscillation term appearing in the nonlinear equation solution,” Journal of Applied Mathematics, vol. 2012, Article ID 387365, 9 pages, 2012. View at: Publisher Site  Google Scholar
 U. FilobelloNino, H. VazquezLeal, Y. Khan et al., “HPM applied to solve nonlinear circuits: a study case,” Applied Mathematical Sciences, vol. 6, no. 85–88, pp. 4331–4344, 2012. View at: Google Scholar
 H. VazquezLeal, R. CastanedaSheissa, A. Yıldırım et al., “Biparameter homotopybased direct current simulation of multistable circuits,” British Journal of Mathematics & Computer Science, vol. 2, no. 3, pp. 137–150, 2012. View at: Google Scholar
 U. FilobelloNino, Hector VazquezLeal, R. CastanedaSheissa et al., “An approximate solution of Blasius equation by using HPM method,” Asian Journal of Mathematics & Statistics, vol. 5, no. 2, Article ID 103923, pp. 50–59, 2012. View at: Publisher Site  Google Scholar
 H. VázquezLeal, U. FilobelloNiño, R. CastañedaSheissa, L. HernándezMartínez, and A. SarmientoReyes, “Modified HPMs inspired by homotopy continuation methods,” Mathematical Problems in Engineering, vol. 2012, Article ID 309123, 19 pages, 2012. View at: Publisher Site  Google Scholar
 H. VazquezLeal, R. CastanedaSheissa, U. FilobelloNino, A. SarmientoReyes, and J. S. Orea, “High accurate simple approximation of normal distribution integral,” Mathematical Problems in Engineering, vol. 2012, Article ID 124029, 22 pages, 2012. View at: Publisher Site  Google Scholar
 Y. Khan, H. VázquezLeal, L. HernandezMartínez, and N. Faraz, “Variational iteration algorithmII for solving linear and nonlinear ODEs,” International Journal of the Physical Sciences, vol. 7, no. 25, pp. 3099–4002, 2012. View at: Google Scholar
 J. H. He, “Variational iteration method—a kind of nonlinear analytical technique: some examples,” International Journal of NonLinear Mechanics, vol. 34, no. 4, pp. 699–708, 1999. View at: Google Scholar
 J. H. He, “Variational iteration method—some recent results and new interpretations,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 3–17, 2007. View at: Publisher Site  Google Scholar
 J. H. He, “Variational iteration method for autonomous ordinary differential systems,” Applied Mathematics and Computation, vol. 114, no. 23, pp. 115–123, 2000. View at: Google Scholar
 M. Agida and A. S. Kumar, “A Boubaker polynomials expansion scheme solution to random Love's equation in the case of a rational Kernel,” Electronic Journal of Theoretical Physics, vol. 7, no. 24, pp. 319–326, 2010. View at: Google Scholar
 A. Yildirim, S. T. MohyudDin, and D. H. Zhang, “Analytical solutions to the pulsed KleinGordon equation using Modified Variational Iteration Method (MVIM) and Boubaker Polynomials Expansion Scheme (BPES),” Computers and Mathematics with Applications, vol. 59, no. 8, pp. 2473–2477, 2010. View at: Publisher Site  Google Scholar
 J. Ghanouchi, H. Labiadh, and K. Boubaker, “An attempt to solve the heat transfer equation in a model of pyrolysis spray using 4qorder mboubaker polynomials,” International Journal of Heat and Technology, vol. 26, no. 1, pp. 49–53, 2008. View at: Google Scholar
 S. Slama, J. Bessrour, K. Boubaker, and M. Bouhafs, “A dynamical model for investigation of A3 point maximal spatial evolution during resistance spot welding using Boubaker polynomials,” The European Physical Journal Applied Physics, vol. 44, no. 3, pp. 317–322, 2008. View at: Publisher Site  Google Scholar
 S. Slama, M. Bouhafs, and K. B. Mahmoud, “A boubaker polynomials solution to heat equation for monitoring A3 point evolution during resistance spot welding,” International Journal of Heat and Technology, vol. 26, no. 2, pp. 141–146, 2008. View at: Google Scholar
 S. Lazzez, K. B. Ben Mahmoud, S. Abroug, F. Saadallah, and M. Amlouk, “A Boubaker polynomials expansion scheme (BPES)related protocol for measuring sprayed thin films thermal characteristics,” Current Applied Physics, vol. 9, no. 5, pp. 1129–1133, 2009. View at: Publisher Site  Google Scholar
 T. Ghrib, K. Boubaker, and M. Bouhafs, “Investigation of thermal diffusivitymicrohardness correlation extended to surfacenitrured steel using Boubaker polynomials expansion,” Modern Physics Letters B, vol. 22, no. 29, pp. 2893–2907, 2008. View at: Publisher Site  Google Scholar
 S. Fridjine, K. B. Ben Mahmoud, M. Amlouk, and M. Bouhafs, “A study of sulfur/selenium substitution effects on physical and mechanical properties of vacuumgrown ZnS1xSex compounds using Boubaker polynomials expansion scheme (BPES),” Journal of Alloys and Compounds, vol. 479, no. 12, pp. 457–461, 2009. View at: Publisher Site  Google Scholar
 C. Khélia, K. Boubaker, T. B. Nasrallah, M. Amlouk, and S. Belgacem, “Morphological and thermal properties of βSnS2 sprayed thin films using Boubaker polynomials expansion,” Journal of Alloys and Compounds, vol. 477, no. 12, pp. 461–467, 2009. View at: Publisher Site  Google Scholar
 K. B. Mahmoud and M. Amlouk, “The 3D AmloukBoubaker expansivityenergy gapVickers hardness abacus: a new tool for optimizing semiconductor thin film materials,” Materials Letters, vol. 63, no. 12, pp. 991–994, 2009. View at: Publisher Site  Google Scholar
 M. Dada, O. B. Awojoyogbe, and K. Boubaker, “Heat transfer spray model: an improved theoretical thermal timeresponse to uniform layers deposit using Bessel and Boubaker polynomials,” Current Applied Physics, vol. 9, no. 3, pp. 622–624, 2009. View at: Publisher Site  Google Scholar
 S. A. H. A. E. Tabatabaei, T. Zhao, O. B. Awojoyogbe, and F. O. Moses, “Cutoff cooling velocity profiling inside a keyhole model using the Boubaker polynomials expansion scheme,” Heat and Mass Transfer, vol. 45, no. 10, pp. 1247–1251, 2009. View at: Publisher Site  Google Scholar
 A. Belhadj, J. Bessrour, M. Bouhafs, and L. Barrallier, “Experimental and theoretical cooling velocity profile inside laser welded metals using keyhole approximation and Boubaker polynomials expansion,” Journal of Thermal Analysis and Calorimetry, vol. 97, no. 3, pp. 911–915, 2009. View at: Publisher Site  Google Scholar
 A. Belhadj, O. F. Onyango, and N. Rozibaeva, “Boubaker polynomials expansion schemerelated heat transfer investigation inside keyhole model,” Journal of Thermophysics and Heat Transfer, vol. 23, no. 3, pp. 639–640, 2009. View at: Publisher Site  Google Scholar
 P. Barry and A. Hennessy, “Meixnertype results for Riordan arrays and associated integer sequences,” Journal of Integer Sequences, vol. 13, no. 9, pp. 1–34, 2010. View at: Google Scholar
 A. S. Kumar, “An analytical solution to applied mathematicsrelated Love's equation using the Boubaker polynomials expansion scheme,” Journal of the Franklin Institute, vol. 347, no. 9, pp. 1755–1761, 2010. View at: Publisher Site  Google Scholar
 S. Fridjine and M. Amlouk, “A new parameter: an ABACUS for optimizig functional materials using the Boubaker polynomials expansion scheme,” Modern Physics Letters B, vol. 23, no. 17, pp. 2179–2191, 2009. View at: Publisher Site  Google Scholar
 M. Benhaliliba, C. E. Benouis, K. Boubaker, M. Amlouk, and A. Amlouk, “A new guide to thermally optimized doped oxides monolayer spraygrown solar cells: the amloukboubaker optothermal expansivity Ψab,” in Solar Cells—New Aspects and Solutions, L. A. Kosyachenko, Ed., pp. 27–41, InTech, 2011. View at: Google Scholar
 A. Milgram, “The stability of the Boubaker polynomials expansion scheme (BPES)based solution to LotkaVolterra problem,” Journal of Theoretical Biology, vol. 271, no. 1, pp. 157–158, 2011. View at: Publisher Site  Google Scholar
 H. Rahmanov, “A solution to the non lLinear kortewegdevries equation in the particular case dispersionadsorption problem in porous media using the spectral boubaker polynomials expansion scheme (BPES),” Studies in Nonlinear Sciences, vol. 2, no. 1, pp. 46–49, 2011. View at: Google Scholar
 M. Koksal and S. Herdem, “Analysis of nonlinear circuits by using differential Taylor transform,” Computers and Electrical Engineering, vol. 28, no. 6, pp. 513–525, 2002. View at: Publisher Site  Google Scholar
 U. FilobelloNino, H. VazquezLeal, Y. Khan et al., “Perturbation method and LaplacePadé approximation to solve nonlinear problems,” Miskolc Mathematical Notes. In press. View at: Google Scholar
 H. VazquezLeal, U. FilobelloNino, A. Yildirim et al., “Transient and DC approximate expressions for diode circuits,” IEICE Electronics Express, vol. 9, no. 6, pp. 522–530, 2012. View at: Publisher Site  Google Scholar
 H. VazquezLeal, L. HernandezMartinez, A. SarmientoReyes, R. CastañedaSheissa, and A. GallardoDelAngel, “Homotopy method with a formal stop criterion applied to circuit simulation,” IEICE Electronics Express, vol. 8, no. 21, pp. 1808–1815, 2011. View at: Publisher Site  Google Scholar
 H. VazquezLeal, L. HernandezMartinez, and A. SarmientoReyes, “Doublebounded homotopy for analysing nonlinear resistive circuits,” in Proceeding of the IEEE International Symposium on Circuits and Systems (ISCAS '05), pp. 3203–3206, Kobe, Japan, May 2005. View at: Publisher Site  Google Scholar
 H. VazquezLeal, L. HernandezMartinez, A. SarmientoReyes, and R. CastañedaSheissa, “Numerical continuation scheme for tracing the double bounded homotopy for analysing nonlinear circuits,” in Proceedings of the International Conference on Communications, Circuits and Systems, pp. 1122–1126, Hong Kong, China, May 2005. View at: Google Scholar
 W. H. Enright, K. R. Jackson, S. P. Norsett, and P. G. Thomsen, “Interpolants for rungekutta formulas,” ACM Transactions on Mathematical Software, vol. 12, no. 3, pp. 193–218, 1986. View at: Publisher Site  Google Scholar
 E. Fehlberg, “Klassischerungekuttaformelnvierter und niedrigererordnung mitschrittweitenkontrolle und ihreanwendung auf waermeleitungsprobleme,” Computing, vol. 6, no. 12, pp. 61–71, 1970. View at: Publisher Site  Google Scholar
 E. TleloCuautle, J. M. MuñozPacheco, and J. MartínezCarballido, “Frequency scaling simulation of Chua's circuit by automatic determination and control of stepsize,” Applied Mathematics and Computation, vol. 194, no. 2, pp. 486–491, 2007. View at: Publisher Site  Google Scholar
 L. Portero, A. Arrarás, and J. C. Jorge, “Variable stepsize fractional step RungeKutta methods for timedependent partial differential equations,” Applied Numerical Mathematics, vol. 62, no. 10, pp. 1463–1476, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 H. VazquezLeal et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.