arxiv:nlin/ v1 [nlin.ps] 12 Jul 2001

Similar documents
Second Order Lax Pairs of Nonlinear Partial Differential Equations with Schwarzian Forms

arxiv:nlin/ v2 [nlin.si] 14 Sep 2001

A Note on Nonclassical Symmetries of a Class of Nonlinear Partial Differential Equations and Compatibility

HOMOTOPY ANALYSIS METHOD FOR SOLVING KDV EQUATIONS

Exact Solutions of Supersymmetric KdV-a System via Bosonization Approach

A NEW VARIABLE-COEFFICIENT BERNOULLI EQUATION-BASED SUB-EQUATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS

New explicit solitary wave solutions for (2 + 1)-dimensional Boussinesq equation and (3 + 1)-dimensional KP equation

Generalized bilinear differential equations

Periodic, hyperbolic and rational function solutions of nonlinear wave equations

Bäcklund transformation and soliton solutions in terms of the Wronskian for the Kadomtsev Petviashvili-based system in fluid dynamics

Breaking soliton equations and negative-order breaking soliton equations of typical and higher orders

Three types of generalized Kadomtsev Petviashvili equations arising from baroclinic potential vorticity equation

No. 11 A series of new double periodic solutions metry constraint. For more details about the results of this system, the reader can find the

Biao Li a,b, Yong Chen a,b, and Hongqing Zhang a,b. 1. Introduction

Exact Solutions for a Fifth-Order Two-Mode KdV Equation with Variable Coefficients

Bäcklund transformation and special solutions for Drinfeld Sokolov Satsuma Hirota system of coupled equations

Building Generalized Lax Integrable KdV and MKdV Equations with Spatiotemporally Varying Coefficients

Soliton solutions of Hirota equation and Hirota-Maccari system

EXACT BREATHER-TYPE SOLUTIONS AND RESONANCE-TYPE SOLUTIONS OF THE (2+1)-DIMENSIONAL POTENTIAL BURGERS SYSTEM

Prolongation structure for nonlinear integrable couplings of a KdV soliton hierarchy

Similarity Reductions of (2+1)-Dimensional Multi-component Broer Kaup System

A New Integrable Couplings of Classical-Boussinesq Hierarchy with Self-Consistent Sources

A MULTI-COMPONENT LAX INTEGRABLE HIERARCHY WITH HAMILTONIAN STRUCTURE

Exact Solutions of Kuramoto-Sivashinsky Equation

Homotopy perturbation method for the Wu-Zhang equation in fluid dynamics

Periodic and Soliton Solutions for a Generalized Two-Mode KdV-Burger s Type Equation

Modified Simple Equation Method and its Applications for some Nonlinear Evolution Equations in Mathematical Physics

Painlevé Analysis and Darboux Transformation for a Variable-Coefficient Boussinesq System in Fluid Dynamics with Symbolic Computation

Yong Chen a,b,c,qiwang c,d, and Biao Li c,d

Numerical Simulation of the Generalized Hirota-Satsuma Coupled KdV Equations by Variational Iteration Method

Invariant Sets and Exact Solutions to Higher-Dimensional Wave Equations

Auto-Bäcklund transformation and exact solutions for compound KdV-type and compound KdV Burgers-type equations with nonlinear terms of any order

Lump solutions to dimensionally reduced p-gkp and p-gbkp equations

New Exact Solutions for the (2+1)-dimensional Boiti-Leon-Manna-Pempinelli Equation

Boundary value problems for integrable equations compatible with the symmetry algebra

A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning

Rational Form Solitary Wave Solutions and Doubly Periodic Wave Solutions to (1+1)-Dimensional Dispersive Long Wave Equation

arxiv: v1 [math-ph] 25 Jul Preliminaries

Differential characteristic set algorithm for PDEs symmetry computation, classification, decision and extension

MACSYMA PROGRAM FOR THE PAINLEVÉ TEST FOR NONLINEAR ORDINARY AND PARTIAL DIFFERENTIAL EQUATIONS

New Analytical Solutions For (3+1) Dimensional Kaup-Kupershmidt Equation

Elsayed M. E. Zayed 1 + (Received April 4, 2012, accepted December 2, 2012)

Symmetry Reductions of (2+1) dimensional Equal Width. Wave Equation

New Application of the (G /G)-Expansion Method to Excite Soliton Structures for Nonlinear Equation

New Exact Solutions of the Modified Benjamin-Bona-Mahony Equation Yun-jie YANG and Li YAO

Exact Periodic Solitary Wave and Double Periodic Wave Solutions for the (2+1)-Dimensional Korteweg-de Vries Equation*

A multiple Riccati equations rational expansion method and novel solutions of the Broer Kaup Kupershmidt system

On bosonic limits of two recent supersymmetric extensions of the Harry Dym hierarchy

Symbolic Computation and New Soliton-Like Solutions of the 1+2D Calogero-Bogoyavlenskii-Schif Equation

The (2+1) and (3+1)-Dimensional CBS Equations: Multiple Soliton Solutions and Multiple Singular Soliton Solutions

New Exact Travelling Wave Solutions for Regularized Long-wave, Phi-Four and Drinfeld-Sokolov Equations

Grammian and Pfaffian solutions as well as Pfaffianization for a (3+1)-dimensional generalized shallow water equation

arxiv:nlin/ v1 [nlin.si] 7 Sep 2005

EXACT SOLITARY WAVE AND PERIODIC WAVE SOLUTIONS OF THE KAUP-KUPERSCHMIDT EQUATION

Multiple-Soliton Solutions for Extended Shallow Water Wave Equations

Periodic and Solitary Wave Solutions of the Davey-Stewartson Equation

A Note On Solitary Wave Solutions of the Compound Burgers-Korteweg-de Vries Equation

Lax Representations and Zero Curvature Representations by Kronecker Product

Application of Laplace Adomian Decomposition Method for the soliton solutions of Boussinesq-Burger equations

Freedom in the Expansion and Obstacles to Integrability in Multiple-Soliton Solutions of the Perturbed KdV Equation

A note on the G /G - expansion method

The (G'/G) - Expansion Method for Finding Traveling Wave Solutions of Some Nonlinear Pdes in Mathematical Physics

arxiv:nlin/ v2 [nlin.si] 9 Oct 2002

New families of non-travelling wave solutions to a new (3+1)-dimensional potential-ytsf equation

Exact Travelling Wave Solutions of the Coupled Klein-Gordon Equation by the Infinite Series Method

New approach for tanh and extended-tanh methods with applications on Hirota-Satsuma equations

A Generalized Extended F -Expansion Method and Its Application in (2+1)-Dimensional Dispersive Long Wave Equation

The General Form of Linearized Exact Solution for the KdV Equation by the Simplest Equation Method

(Received 05 August 2013, accepted 15 July 2014)

Homotopy Analysis Method for Nonlinear Differential Equations with Fractional Orders

Anti-synchronization of a new hyperchaotic system via small-gain theorem

Application of the trial equation method for solving some nonlinear evolution equations arising in mathematical physics

Alexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011

Extended tanh-function method and its applications to nonlinear equations

The Traveling Wave Solutions for Nonlinear Partial Differential Equations Using the ( G. )-expansion Method

EXACT SOLUTION TO TIME FRACTIONAL FIFTH-ORDER KORTEWEG-DE VRIES EQUATION BY USING (G /G)-EXPANSION METHOD. A. Neamaty, B. Agheli, R.

Multi-Soliton Solutions to Nonlinear Hirota-Ramani Equation

New Integrable Decomposition of Super AKNS Equation

Lump-type solutions to nonlinear differential equations derived from generalized bilinear equations

Math 575-Lecture 26. KdV equation. Derivation of KdV

New Homoclinic and Heteroclinic Solutions for Zakharov System

The Riccati equation with variable coefficients expansion algorithm to find more exact solutions of nonlinear differential equations

The first three (of infinitely many) conservation laws for (1) are (3) (4) D t (u) =D x (3u 2 + u 2x ); D t (u 2 )=D x (4u 3 u 2 x +2uu 2x ); D t (u 3

NEW PERIODIC WAVE SOLUTIONS OF (3+1)-DIMENSIONAL SOLITON EQUATION

A Further Improved Tanh Function Method Exactly Solving The (2+1)-Dimensional Dispersive Long Wave Equations

Continuous limits and integrability for a semidiscrete system Zuo-nong Zhu Department of Mathematics, Shanghai Jiao Tong University, P R China

Extended Jacobi Elliptic Function Expansion Method for Nonlinear Benjamin-Bona-Mahony Equations

Hamiltonian partial differential equations and Painlevé transcendents

Exact Solutions of The Regularized Long-Wave Equation: The Hirota Direct Method Approach to Partially Integrable Equations

Exact traveling wave solutions of nonlinear variable coefficients evolution equations with forced terms using the generalized.

An Improved F-Expansion Method and Its Application to Coupled Drinfel d Sokolov Wilson Equation

arxiv: v1 [nlin.si] 10 Oct 2011

arxiv:patt-sol/ v1 25 Sep 1995

The Modified (G /G)-Expansion Method for Nonlinear Evolution Equations

Approximate Similarity Reduction for Perturbed Kaup Kupershmidt Equation via Lie Symmetry Method and Direct Method

Linearization of Mirror Systems

HOMOTOPY PERTURBATION METHOD TO FRACTIONAL BIOLOGICAL POPULATION EQUATION. 1. Introduction

Travelling Wave Solutions for the Gilson-Pickering Equation by Using the Simplified G /G-expansion Method

Symbolic computation of conservation laws for nonlinear partial differential equations in multiple space dimensions

New Approach of ( Ǵ/G ) Expansion Method. Applications to KdV Equation

Transcription:

Higher dimensional Lax pairs of lower dimensional chaos and turbulence systems arxiv:nlin/0107028v1 [nlin.ps] 12 Jul 2001 Sen-yue Lou CCAST (World Laboratory), PO Box 8730, Beijing 100080, P. R. China Physics Department of Shanghai Jiao Tong University, Shanghai 200030, P. R. China Physics Department of Ningbo University, Ningbo 315211, P. R. China Abstract In this letter, a definition of the higher dimensional Lax pair for a lower dimensional system which may be a chaotic system is given. A special concrete (2+1)-dimensional Lax pair for a general (1+1)-dimensional three order autonomous partial differential equation is studied. The result shows that any (1+1)-dimensional three order semi-linear autonomous system (no matter it is integrable or not) possesses infinitely many (2+1)-dimensional Lax pairs. Especially, every solution of the KdV equation and the Harry-Dym equation with their space variable being replaced by the field variable can be used to obtain a (2+1)-dimensional Lax pair of any three order (1+1)-dimensional semi-linear equation. In the nonlinear mathematical physics, if an (n+1)-dimensional nonlinear system can be considered as a consistent condition of an(n+1)-dimensional linear system, then many types of interesting properties of the nonlinear system can be obtained by analyzing the linear system. The linear system is called as the Lax pair of the original nonlinear system[1]. Inthisletter, weusetheconceptofthelaxpairunderamoregeneralmeaning. Ifthecompatible condition [L 1, L 2 ]ψ (L 1 L 2 L 2 L 1 )ψ = 0, (1) of a pair of linear equations L 1 (x 1,x 2,...,x n,t,u(x 1,x 2,...,x m,t))ψ(x 1,x 2,...,x n,t) L 1 ψ = 0, (2) L 2 (x 1,x 2,...,x n,t,u(x 1,x 2,...,x m,t))ψ(x 1,x 2,...,x n,t) L 2 ψ = 0, m < n, (3) Email: sylou@mail.sjtu.edu.cn 1

is only a ψ-independent nonlinear equation of u(x 1,x 2,...,x m,t) u F(u) = 0, (4) then we call the system (2) and (3) the (n+1)-dimensional Lax pair of the (m+1)-dimensional nonlinear equation (4) (or equivalently (1)). To find some nontrivial examples, we restrict ourselves to discussion the case for L 1 2 x y F y(t, y, u,u x ), u = u(x,t) (5) L 2 t + 3 x 3 G(t, y, u, u x, u xx ) x, (6) where F F(t, y, u,u x ), G G(t, y, u, u x,,u xx ) should be selected appropriately such that Eq. (1) is only a (1+1)-dimensional partial differential equation of u. For the notation simplicity later, the undetermined function in the equation (5) is written as a partial derivative form. In the next step we try to find some autonomous evolution models u t = H 0 (u, u x, u xx, u xxx, u xxxx,...,) (7) which are given by the compatibility conditions of (2) and (3) with (5) and (6) or equivalently (1) with (5) and (6). ψ xyt = ψ txy (8) Using the relations (5), (6) and (7), the compatibility condition (8) has the following form W 1 (t,y,u,u x,...)ψ +W 2 (t,y,u,u x,...)ψ x +(3F ux u xx G+3F u u x ) y ψ xx = 0, (9) where W 1 (t,y,u,u x,...) W 1 and W 2 (t,y,u,u x,...) W 2 are two complicated functions of the indicated variables and F ux and F u denotes the partial derivatives of F with respect to u x and u. The ψ-independent condition of (8) requires that three terms of (9) should all be zero. Vanishing the last term of (9) yields G = 3F ux u xx +3F u u x +G 0 (t,u,u x,u xx ), (10) whereg 0 (t,u,u x,u xx ) is a y-independentarbitraryfunction of the indicated variables. Substituting (10) into (9), we find that W 2 becomes zero identically while W 1 = 0 is simplified to W 11 (t,y,u,u x,...,u k,x )+F yux H 0uk,x u k+1,x = 0, k > 3, (11) where u k,x denotes the kth order derivative of u with respect to x, k is the highest derivative order of u included in H 0 (u, u x, u xx, u xxx, u xxxx,...,) H 0 and W 11 is u K,x (K k +1) independent 2

complicated function. Because (11) should be true for the general solution u of (7), two terms of (11) should all be zero. For F yux 0, vanishing second term of (11), we have the result that H is u k,x independent for k 4. For k = 3, (8) with (10) has the form W 11 (t,y,u,u x,u xx,u xxx )+F yux (H 0uxxx +1)u xxxx = 0. (12) Vanishing the second term of (12) yields the result that the only possible form of H 0 can be written as for F yux 0. Substituting (13) back to (12) yields H 0 = u xxx +H(u, u x, u xx ) (13) W 12 (t,y,u,u x,u xx )+(F yux H uxx 3F y F ux +3F yuux u x +3F yuxu x u xx F y G 0uxx )u xxx = 0. (14) Because W 12 (t,y,u,u x,u xx ) is u xxx independent, the second term of (14) should be also zero. Vanishing the coefficients of u xxx of (14) and partial integrating once with respect to u xx we obtain 3 2 F yu xu x u 2 xx +3(F yuu x u x F y F ux )u xx F y G 0 +F yux H +F 1 (t,y,u,u x ) = 0, (15) where F 1 (t,y,u,u x ) F 1 is an integrating function. Because F and F 1 in (15) are y-dependent while all other functions in (15) are all y-independent, we can conclude that the only possible solution of (15) has the form for general H with G 0 = 3 2 Γ 2(t,u,u x )u 2 xx +Γ 3(t,u,u x )u xx +Γ 1 (t,u,u x )H +Γ 4 (t,u,u x ), (16) F yux = F y Γ 1 (t,u,u x ), (17) F yuxu x = F y Γ 2 (t,u,u x ), (18) 3(F yuux u x F y F ux ) = F y Γ 3 (t,u,u x ), (19) F 1 (t,y,u,u x ) = F y Γ 4 (t,u,u x ), (20) and Γ 1 (t,u,u x ), Γ 2 (t,u,u x ), Γ 3 (t,u,u x ) and Γ 4 (t,u,u x ) being four unknown functions which will be determined later. According to the compatibility condition F uux = F uxu and the fact that F is the only y- dependent function, we can obtain the unique possible solution of (17)-(19): Γ 1 (t,u,u x ) = u 1 x, (21) 3

Γ 2 (t,u,u x ) = 0, (22) Γ 3 (t,u,u x ) = 3G 1ux (t,u,u x ), (23) F = G 1 (t,u,u x )+u x F 0 (t,y,u), (24) F 0u (t,y,u) = 1 2 F 0(t,y,u) 2 +Γ 5 (t,u), (25) where G 1 (t,u,u x ) is an arbitrary function of {t,u,u x } and Γ 5 (t,u) is an arbitrary function of {t,u}. Using the relations (16), (21)-(25), Eq. (12) is simplified to F 0ty (t,y,u) F 0y (t,y,u) [2u3 x Γ 5(t,u)+Γ 4 (t,u,u x )u x +3G 1u (t,u,u x )u 2 x ]F 0(t,y,u) u xx [3u x (2Γ 5 (t,u)+g 1uux (t,u,u x ))+Γ 4ux (t,u,u x )+6G 1u (t,u,u x )+u 1 x Γ 4 (t,u,u x )] u x Γ 4u (t,u,u x ) 2u 3 x Γ 5u 3u 2 x G 1uu(t,u,u x ) = 0. (26) Because F 0 (t,y,u) is y-dependent, and other functions in (26) are y and u xx independent, (26) is true only for F 0ty (t,y,u) = F 0y (t,y,u)(γ 6 (t,u)f 0 (t,y,u)+γ 7 (t,u)), (27) 3u x (2Γ 5 (t,u)+g 1uux (t,u,u x ))+Γ 4ux (t,u,u x )+6G 1u (t,u,u x )+u 1 x Γ 4 (t,u,u x ) = 0, (28) Γ 7 (t,u) u x Γ 4u (t,u,u x ) 2u 3 x Γ 5u 3u 2 x G 1uu(t,u,u x ) = 0, (29) 2u 3 x Γ 5(t,u)+Γ 4 (t,u,u x )u x +3G 1u (t,u,u x )u 2 x Γ 6(t,u) = 0, (30) where Γ 6 (t,u) and Γ 7 (t,u) are arbitrary functions of {t,u}. The general solution of (28)-(30) reads Γ 4 (t,u,u x ) = 3u x G 1u (t,u,u x ) 2u 2 xγ 5 (t,u)+u 1 x Γ 6 (t,u), (31) Γ 7 (t,u) = Γ 6u (t,u). (32) Integrating (27) once with respect to y leads to F 0t (t,y,u) = 1 2 Γ 6(t,u)F 0 (t,y,u) 2 +Γ 6u (t,u)f 0 (t,y,u)+γ 8 (t,u), (33) where Γ 8 (t,u) is a further integrating function. Because one function, F 0 (t,y,u), should satisfy two equations (25) and (33), the compatibility condition F tu = F ut yields the constraints Γ 8 (t,u) = Γ 6 (t,u)γ 5 (t,u)+γ 6uu (t,u), (34) and Γ 5t (t,u) = Γ 6uuu (t,u)+γ 6 (t,u)γ 5u (t,u)+2γ 5 (t,u)γ 6u (t,u). (35) Finally, we write down the obtained results here. A (1+1)-dimensional three order equation u t = u xxx +H(u,u x,u xx ) (36) 4

with an arbitrary nonlinear interaction term H(u,u x,u xx ) possesses a (2+1)-dimensional Lax pair where ψ xy u x F 0y (t,y,u)ψ = 0, (37) ψ t +ψ xxx [3F 0 (t,y,u)u xx +3u 2 xf 0u (t,y,u)+h 1 (t,u,u x,u xx )]ψ x = 0, (38) H 1 (t,u,u x,u xx ) = u 1 x (H(u,u x,u xx )+Γ 6 (t,u)) 2u 2 xγ 5 (t,u). (39) In the Lax pair equations (37) and (38), the function F 0 (t,y,u) F 0 is determined by a pair of the consistent Riccati equations F 0t (t,y,u) = 1 2 Γ 6(t,u)F 0 (t,y,u) 2 +Γ 6u (t,u)f 0 (t,y,u)+γ 6 (t,u)γ 5 (t,u)+γ 6uu, (40) F 0u (t,y,u) = 1 2 F 0(t,y,u) 2 +Γ 5 (t,u) (41) while the consistent condition of (40) and (41) gives the constraint equation (35) for the functions Γ 6 (t,u) and Γ 5 (t,u). It is known that for the usual Lax pair of a given nonlinear model, there may be infinitely many Lax pairs[2]. The similar situation appears for the extended higher dimensional Lax pairs. For instance, for the given evolution system (36), the arbitrariness in the selections of the functions Γ 5 (t,u) and Γ 6 (t,u) and thesolutions of F 0 (t,y,u) means that the evolution equation (36) possesses infinitely many different kinds of Lax pairs. To give out some more concrete results, we may fix the functions Γ 6 (t,u) and Γ 5 (t,u). The first interesting selection is Γ 6 (t,u) = Γ 5 (t,u). (42) Under the selection (42), (36) becomes the well known Korteweg de-vries (KdV) equation Γ 5t (t,u) = Γ 5uuu (t,u)+3γ 5 (t,u)γ 5u (t,u). (43) with the independent variables {t,u} and the linearized system of (40) and (41) by using the Cole- Hopf transformation is just the Lax pair of the KdV equation (after neglecting the variable y in (40) and (41)). Now substituting every known special solution of the KdV equation and the related solution of the psudopotential F 0 of (40) and (41) into (37) and (38), we get a concrete Lax pair of the system (36). The simplest trivial solution of the KdV equation is and the related solution of (40) and (41) reads Γ 5 (t,u) = Γ 6 (t,u) = 0, (44) 2 F 0 (t,y,u) = u+q(y), (45) 5

where q(y) is an arbitrary function of y. The second interesting selection is Γ 6 (t,u) = Γ 1/2 5 (t,u). (46) In this case, (36) becomes the well known Harry-Dym (HD) equation Γ 5t (t,u) = [Γ 1/2 5 (t,u)] uuu. (47) with the independent variables {t,u} and the linearized system of (40) and (41) is the Lax pair of the HD equation. That means substituting every known special solution of the HD equation and the related solution of the psudopotential F 0 of (40) and (41) into (37) and (38) will yield a concrete Lax pair of the semi-linear system (36). Generally, the evolution equation (36) is nonintegrable under the traditional meanings. In Ref. [3], the authors had claimed that there exist only six nonequivalent three order semi-linear integrable models under the usual meanings. In other words, the evolution equation (36) is chaotic or turbulent for most of the selections of H. The results of this letter show us that no matter of the model (36) is integrable or not, it may have infinitely many (2+1)-dimensional Lax pairs. A simple special nonintegrable example of (36) is the so-called KdV-Burgers (KdVB) equation[4], u t +uu x νu xx +u xxx = 0. (48) The KdVB is one of the possible candidate to describe the turbulence phenomena in fluid physics and plasma physics[5, 4]. Though one has not yet find any (1+1)-dimensional Lax pair of the turbulence system KdVB equation, we can do find many types of (2+1)-dimensional Lax pairs. Especially, every solution of the KdV equation and the HD equation can be used to form a (2+1)- dimensional Lax pair of the KdVB equation. In summary, a lower dimensional nonlinear system may have some (perhaps infinitely many) types of higher dimensional Lax pairs. In terms of a special example we show that any three order semi-linear equation (no mater it is integrable or not) possesses infinitely many Lax pairs. The next interesting important problem should be studied in the future work may be: What kinds of information about lower dimensional systems especially the lower dimensional turbulent and chaotic systems can be obtained from some types of special higher dimensional Lax pairs? The work was supported by the Outstanding Youth Foundation and the National Natural Science Foundation of China (Grant. No. 19925522), the Research Fund for the Doctoral Program of Higher Education of China (Grant. No. 2000024832) and the Natural Science Foundation of Zhejiang Province, China. The author is in debt to thanks the helpful discussions with the 6

professors Q. P. Liu, G-x Huang and C-p Sun and the Drs. X-y Tang, S-l Zhang, C-l Chen and B. Wu. 7

References [1] M. J. Ablowitz and P. A. Clarkson, Solitons Nonlinear Evolution Equations and Inverse Scattering, London Mathematical Society Lecture Note Series 149 (1991 Cambridge University Press). [2] S-Lou and X-b Hu, J. Math. Phys. 38(1997)6401. [3] S. I. Svinolupov, V. V. Sokolov, and R. I. Yamilov, Sov. Math. Dokl., 28 (1983) 165; S. I. Svinolupov and V. V. Sokolov, Func. Anal. Appl. 16 (1982) 317. [4] Y. Nakamura, H. Bailung and P. K. Shukla, Phys. Rev. Lett. 83 (1999) 1602; E. P. Raposo and D. Bazeia, Phys. Lett. A253 (1999) 151; G. Karch, Nonlinear Analysis: Theory Methods & Applications, 35 (1999)199; [5] S-d Liu and S-k Liu, Sincia Sinica, A 9 (1991) 938 (in Chinese); N. Antar, and H. Demiray, Int. J. Engineering Sci. 37 (1999) 1859. 8