Haar wavelet solutions of nonlinear oscillator equations

In this paper, we present a numerical scheme using uniform Haar wavelet approximation and quasilinearization process for solving some nonlinear oscillator equations. In our proposed work, quasilinearization technique is first applied through Haar wavelets to convert a nonlinear differential equation into a set of linear algebraic equations. Finally, to demonstrate the validity of the proposed method, it has been applied on three type of nonlinear oscillators namely Duffing, Van der Pol, and Duffing–van der Pol. The obtained responses are presented graphically and compared with available numerical and analytical solutions found in the literature. The main advantage of uniform Haar wavelet series with quasilinearization process is that it captures the behavior of the nonlinear oscillators without any iteration. The numerical problems are considered with force and without force to check the efficiency and simple applicability of method on nonlinear oscillator problems.     

An inexact Newton method for nonlinear two-point boundary-value problems

The method of quasilinearization for nonlinear two-point boundary-value problems is an application of Newton's method to a nonlinear differential operator equation. Since the linear boundary-value problem to be solved at each iteration must be discretized, it is natural to consider quasilinearization in the framework of an inexact Newton method. More importantly, each linear problem is only a local model of the nonlinear problem, and so it is inefficient to try to solve the linear problems to full accuracy. Conditions on size of the relative residual of the linear differential equation can then be specified to guarantee rapid local convergence to the solution of the nonlinear continuous problem. If initial-value techniques are used to solve the linear boundary-value problems, then an integration step selection scheme is proposed so that the residual criteria are satisfied by the approximate solutions. Numerical results are presented that demonstrate substantial computational savings by this type of economizing on the intermediate problems.This work was supported in part by DOE Contract DE-AS05-82-ER13016 and NSF Grant RII-89-17691 and was part of the author's doctoral thesis at Rice University. It is a pleasure to thank the author's thesis advisors, Professor R. A. Tapia and Professor J. E. Dennis, Jr.     

B-Spline Collocation Method for Nonlinear Singularly-Perturbed Two-Point Boundary-Value Problems

A B-spline collocation method is presented for nonlinear singularly-perturbed boundary-value problems with mixed boundary conditions. The quasilinearization technique is used to linearize the original nonlinear singular perturbation problem into a sequence of linear singular perturbation problems. The B-spline collocation method on piecewise uniform mesh is derived for the linear case and is used to solve each linear singular perturbation problem obtained through quasilinearization. The fitted mesh technique is employed to generate a piecewise uniform mesh, condensed in the neighborhood of the boundary layers. The convergence analysis is given and the method is shown to have second-order uniform convergence. The stability of the B-spline collocation system is discussed. Numerical experiments are conducted to demonstrate the efficiency of the method.     

Numerical solution of the Falkner-Skan equation based on quasilinearization

We present two iterative methods for solving the Falkner-Skan equation based on the quasilinearization method. We formulate the original problem as a new free boundary value problem. The truncated boundary depending on a small parameter is an unknown free boundary and has to be determined as part of solution. Using a change of variables, the free boundary value problem is transformed to a problem defined on [0, 1]. We apply the quasilinearization method to solve the resulting nonlinear problem. Then we propose two different iterative algorithms by means of a cubic spline solver. Numerical results for various instances are compared with those reported previously in the literature. The comparisons show the accuracy, robustness and efficiency of the presented methodology.     

Nonlinear asymptotic analysis of a system of two free coupled oscillators with cubic nonlinearities

In this paper a two degrees of freedom undamped nonlinear system of two unforced coupled oscillators with cubic nonlinearities is analyzed. Through a decoupling procedure and using admissible functional transformations it is proved that this system can be reduced to an intermediate second order nonlinear ordinary differential equation (ODE) connecting both displacements to each other. By nonlinear asymptotic approximations the above equation can be further reduced to new nonlinear ODE that can be analytically solved. The solutions in the physical plane are extracted in parametric form. As generalization, the model of a damped system of two masses connected with stiffness with linear and nonlinear coefficient of rigidities respectively is analyzed and exact analytical solutions are extracted. Finally an application has been given in the case of a two mass system with cubic strong non-linearity.     

Numerical treatment of singularly perturbed two point boundary value problems using initial-value method

In this paper, we describe an initial-value method for linear and nonlinear singularly perturbed boundary value problems in the interval [p,q]. For linear problems, the required approximate solution is obtained by solving the reduced problem and one initial-value problems directly deduced from the given problem. For nonlinear problems the original second-order nonlinear problem is linearized by using quasilinearization method. Then this linear problem is solved as previous method. The present method has been implemented on several linear and non-linear examples which approximate the exact solution. We also present the approximate and exact solutions graphically.     

Initial-value methods for second-order singularly perturbed boundary-value problems

In this paper, we present a numerical method for solving linear and nonlinear second-order singularly perturbed boundary-value-problems. For linear problems, the method comes from the well-known WKB method. The required approximate solution is obtained by solving the reduced problem and one or two suitable initial-value problems, directly deduced from the given problem. For nonlinear problems, the quasilinearization method is applied. Numerical results are given showing the accuracy and feasibility of the proposed method.This work was supported in part by the Consiglio Nazionale delle Ricerche (Contract No. 86.02108.01 and Progetto Finalizzatto Sistemi Informatia e Calcolo Paralello, Sottoprogetto 1), and in part by the Ministero della Pubblica Istruzione, Rome, Italy.     

Gegenbauer wavelet collocation method for the extended Fisher-Kolmogorov equation in two dimensions

Gegenbauer wavelets operational matrices play an important role in the numeric solution of differential equations. In this study, operational matrices of rth integration of Gegenbauer wavelets are derived and used to obtain an approximate solution of the nonlinear extended Fisher-Kolmogorov (EFK) equation in two-space dimension. Nonlinear EFK equation is converted into the linear partial differential equation by quasilinearization technique. Numerical examples have shown that present method is convergent even in the case of a small number of grid points. The results of the presented method are in a good agreement with the results in literature.     

A numerical study of boundary layer flow of viscous incompressible fluid past an inclined stretching sheet and heat transfer using nonpolynomial spline method

In the present paper, we study the boundary layer flow of viscous incompressible fluid over an inclined stretching sheet with body force and heat transfer. Considering the stream function, we convert the boundary layer equation into nonlinear third-order ordinary differential equation together with appropriate boundary conditions in an infinite domain. The nonlinear boundary value problem has been linearized by using the quasilinearization technique. Then, we develop a nonpolynomial spline method, which is used to solve the flow problem. The convergence analysis of the method is also discussed. We study the velocity function for different angles of inclination and Froude number with the help of various graphs and tables. Then using these in heat convection flow, we obtain the expression for temperature field. Skin friction is also calculated. The various results have been given in tables. At last, we calculated the Nusselt number.     

The effect of a reinforcing ring on the stress-strain state of flexible cylindrical shells

We investigate the effect of a reinforcing ring on the stress-strain state of a cylindrical shell in the geometrically nonlinear problem with a nonaxisymmetric load on the edge. The nonlinear boundary-value problem is reduced to a sequence of linear problems by the quasilinearization method. The linear problems are solved by the discrete orthogonalization method. The results obtained using linear and nonlinear theory are compared.Translated from Vychislitel'naya i Prikladnaya Matematika, No. 55, pp. 92–96, 1985.     

Galerkin's procedure,quasilinearization, and nonlinear boundary-value problems

Among the popular and successful techniques for solving boundary-value problems for nonlinear, ordinary differential equations (ODE) are quasilinearization and the Galerkin procedure. In this note, it is demonstrated that utilizing the Galerkin criterion followed by the Newton-Raphson scheme results in the same iteration process as that obtained by applying quasilinearization to the nonlinear ODE and then the Galerkin criterion to each linear ODE in the resulting sequence. This equivalence holds for only the Galerkin procedure in the broad class of weighted-residual methods.This work was supported in part by the National Science Foundation, Grant No. GJ-1075.     

Approximate solution to the time–space fractional cubic nonlinear Schrodinger equation

By introducing the fractional derivatives in the sense of Caputo, we use the adomian decomposition method to construct the approximate solutions for the cubic nonlinear fractional Schordinger equation with time and space fractional derivatives. The exact solution of the cubic nonlinear Schrodinger equation is given as a special case of our approximate solution. This method is efficient and powerful in solving wide classes of nonlinear evolution fractional order equation.     

A model trust-region modification of Newton's method for nonlinear two-point boundary-value problems

The method of quasilinearization for nonlinear two-point boundary-value problems is Newton's method for a nonlinear differential operator equation. A model trust-region approach to globalizing the quasilinearization algorithm is presented. A double-dogleg implementation yields a globally convergent algorithm that is robust in solving difficult problems.This work was supported in part by DOE Contract DE-AS05-82-ER13016 and NSF Grant RII-89-17691 and was part of the author's doctoral thesis at Rice University. It is a pleasure to thank the author's thesis advisors, Professor J. E. Dennis, Jr., and Professor R. A. Tapia.     

Optimal input system identification for nonlinear dynamic systems

The estimation accuracy for nonlinear dynamic system identification is known to be maximized by the use of optimal inputs. Few examples of the design of optimal inputs for nonlinear dynamic systems are given in the literature, however. The performance criterion is selected such that the sensitivity of the measured state variables to the unknown parameters is maximized. The application of Pontryagin's maximum principle yields a nonlinear two-point boundary-value problem. In this paper, the boundary-value problem for a simple nonlinear example is solved using two different methods, the method of quasilinearization and the Newton-Raphson method. The estimation accuracy is discussed in terms of the Cramer-Rao lower bound.     

Third-order iterative methods free from second derivatives for nonlinear equation

In this paper, we suggest and analyze a new two-step iterative method for solving nonlinear equations, which is called the modified Householder method without second derivatives for nonlinear equation. We also prove that the modified method has cubic convergence. Several examples are given to illustrate the efficiency and the performance of the new method. New method can be considered as an alternative to the present cubic convergent methods for solving nonlinear equations.     

Some modifications of the quasilinearization method with higher-order convergence for solving nonlinear BVPs

In this paper, modifications of the quasilinearization method with higher-order convergence for solving nonlinear differential equations are constructed. A general technique for systematically obtaining iteration schemes of order m (?>?2) for finding solutions of highly nonlinear differential equations is developed. The proposed iterative schemes have convergence rates of cubic, quartic and quintic orders. These schemes were further applied to bifurcation problems and to obtain critical parameter values for the existence and uniqueness of solutions. The accuracy and validity of the new schemes is tested by finding accurate solutions of the one-dimensional Bratu and Frank-Kamenetzkii equations.     

Exact solutions for logistic reaction–diffusion equations in biology

Reaction–diffusion equations with a nonlinear source have been widely used to model various systems, with particular application to biology. Here, we provide a solution technique for these types of equations in N-dimensions. The nonclassical symmetry method leads to a single relationship between the nonlinear diffusion coefficient and the nonlinear reaction term; the subsequent solutions for the Kirchhoff variable are exponential in time (either growth or decay) and satisfy the linear Helmholtz equation in space. Example solutions are given in two dimensions for particular parameter sets for both quadratic and cubic reaction terms.     

Existence and approximation of solutions of forced duffing type integro-differential equations with three-point nonlinear boundary conditions

A generalized quasilinearization technique is developed to obtain an analytic approximation of the solutions of the forced Duffing type integro-differential equation with nonlinear three-point boundary conditions. Monotone sequences of approximate solutions converging uniformly and quadratically to a unique solution of the problem are presented.     

Solving the 2-D elliptic Monge-Ampère equation by a Kansa’s method

In this paper, a Kansa’s method is designed to solve numerically the Monge-Ampère equation. The primitive Kansa’s method is a meshfree method which applying the combination of some radial basis functions (such as Hardy’s MQ) to approximate the solution of the linear parabolic, hyperbolic and elliptic problems. But this method is deteriorated when is used to solve nonlinear partial differential equations. We approximate the solution in some local triangular subdomains by using the combination of some cubic polynomials. Then the given problems can be computed in each subdomains independently. We prove the stability and convergence of the new method for the elliptic Monge-Ampère equation. Finally, some numerical experiments are presented to demonstrate the theoretical results.     

A quasilinearization method for a class of integro-differential equations with mixed nonlinearities

We develop an extended form of the quasilinearization method for an initial value problem involving a nonlinear integro-differential equation with mixed nonlinearities, occurring in distributed-infective and distributed-contact models. In fact, a monotone sequence of iterates converging uniformly and quadratically to a solution of the problem is obtained. Some interesting observations are presented.