Commun Nonlinear Sci Numer Simulat
|
|
- Mark Greene
- 5 years ago
- Views:
Transcription
1 Commun Nonlinear Sci Numer Simulat 6 (2) Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: The scaled boundary FEM for nonlinear problems Zhiliang Lin, Shijun Liao * State Key Laboratory of Ocean Engineering, School of Naval Architecture, Ocean and Civil Engineering, Shanghai JiaoTong University, Shanghai 224, China article info abstract Article history: Received 5 February 2 Received in revised form 4 March 2 Accepted 7 March 2 Available online 5 March 2 Keywords: Scaled boundary finite-element Homotopy analysis method Nonlinear problem The traditional scaled boundary finite-element method (SBFEM) is a rather efficient semianalytical technique widely applied in engineering, which is however valid mostly for linear differential equations. In this paper, the traditional SBFEM is combined with the homotopy analysis method (HAM), an analytic technique for strongly nonlinear problems: a nonlinear equation is first transformed into a series of linear equations by means of the HAM, and then solved by the traditional SBFEM. In this way, the traditional SBFEM is extended to nonlinear differential equations. A nonlinear heat transfer problem is used as an example to show the validity and computational efficiency of this new SBFEM. Ó 2 Elsevier B.V. All rights reserved.. Introduction The scaled boundary finite-element method (SBFEM), developed recently by Wolf and Song [ 3], is a novel semi-analytical method to solve partial differential equations. In the frame of the SBFEM, a scaled boundary coordinate system is introduced. The weighted residual approximation of finite elements is applied in the circumferential direction, and then the governing partial differential equations are transformed to ordinary differential equations in the radial direction. Those ordinary differential equations can be solved analytically in the radial direction. Like the boundary element method, only the boundary of the problem domain is discretized, but no fundamental solution is required. Thus, this semi-analytical method combines the advantages of the finite element and boundary element methods, and besides, it presents appealing features of its own. The SBFEM has been used for the problems of elasto-statics and elasto-dynamics, especially for soil-structure interaction problems in unbounded domains [3,4]. Recently, the SBFEM has been extended to the fluid flow problems [5 8]. Meanwhile, researchers reported successful algorithms for the coupling of SBFEM with other numerical techniques [9 ] to analyze geomechanics problems. By coupling the finite element method and the SBFEM, Doherty and Deeks [9] captured the nonlinearity of problems in the near field. However, the SBFEM only accurately modeled the linear elastic far field response. It must be emphasized that, up to now, all problems solved by SBFEM are governed by linear partial differential equations. The homotopy analysis method (HAM), which is an analytic approach to get series solutions of strongly nonlinear equations, was first proposed by Liao [2] in 992. Thereafter, the HAM has been improved step by step [3 9]. Briefly speaking, by means of the HAM, one constructs a continuous mapping of an initial guess approximation to the exact solution of considered equations. An auxiliary linear operator is chosen to construct such kind of continuous mapping, and a convergencecontrol parameter is used to ensure the convergence of solution series. More importantly, this method enjoys great freedom in choosing initial guesses and auxiliary linear operators. By means of this kind of freedom, a complicated nonlinear problem can be transferred into an infinite number of simple, linear sub-problems, as shown by Liao and Tan [7]. Thus, if the SBFEM is successfully combined with the HAM, it could be extended to many engineering problems governed by nonlinear partial differential equations. This is the motivation of our present work. * Corresponding author. addresses: linzhiliang@sjtu.edu.cn (Z. Lin), sjliao@sjtu.edu.cn (S. Liao) /$ - see front matter Ó 2 Elsevier B.V. All rights reserved. doi:.6/j.cnsns.2.3.5
2 64 Z. Lin, S. Liao / Commun Nonlinear Sci Numer Simulat 6 (2) The aim of this paper is to combine the SBFEM with the HAM and then apply it to solve nonlinear problems. Shortly speaking, we first transfer an original nonlinear problem into an infinite number of linear sub-problems using the HAM, and then solve all those sub-problems by the SBFEM. To show its validity, 2D nonlinear Poisson-type problems are solved, and numerical results demonstrate the efficiency and accuracy of this approach. 2. Theoretical description 2.. Governing equations Consider a 2D Poisson equation with a varying conductivity kðx; yþ, which can be expressed as r ½kðx; yþruðx; yþš ¼ f ðx; yþ; ðx; yþ 2X; ðþ under boundary conditions u ¼ g ðx; yþ; ðx; yþ 2C ¼ g 2ðx; yþ; ðx; yþ 2C 2 ; ð3þ where r is the Hamilton operator, uðx; yþ is the dependent variable, the real functions f ðx; yþ; g ðx; yþ and g 2 ðx; yþ are prescribed functions or given values, the x and y are spatial coordinates, n denotes the normal to the boundary, X the domain and C ¼ C [ C 2 the boundary. Obviously, Eq. () can be rewritten as kðx; yþr 2 uðx; yþþrkðx; yþruðx; yþ ¼f ðx; yþ; ðx; yþ 2X: ð4þ Assuming that kðx; yþ ¼ þ uðx; yþ, i.e., the conductivity is dependent on the unknown uðx; yþ. Then Eq. (4) becomes r 2 uðx; yþþuðx; yþr 2 uðx; yþþruðx; yþruðx; yþ ¼f ðx; yþ; ðx; yþ 2X: ð5þ Thus, the 2D nonlinear Poisson-type problems which we will solve are governed by Eq. (5) with the boundary conditions (2) and (3) Homotopy transformation According to Eq. (5), a nonlinear operator N is defined as follows: N½UŠ ¼r 2 U þ Ur 2 U þ ru ru f ðx; yþ: ð6þ Then Eq. (5) can be written as N½uðx; yþš ¼ ; ðx; yþ 2 X: ð7þ In the frame of the HAM, we first construct the zeroth-order deformation equation ð qþl½uðx; y; qþ u ðx; yþš ¼ c qn½uðx; y; qþš; q 2½; Š; c ; ðx; yþ 2X ð8þ with the boundary conditions Uðx; y; qþ ¼g ðx; yþ; ðx; yþ 2C ; y; ¼ g 2 ðx; yþ; ðx; yþ 2C 2 ; ðþ where c is the convergence-control parameter, u ðx; yþ is an initial guess of uðx; yþ; L is the auxiliary linear operator, and q 2 ½; Š is called homotopy-parameter, respectively. As q increases from to, the solution Uðx; y; qþ deforms continuously from the initial guess u ðx; yþ to the solution uðx; yþ of Eqs. (5), (2) and (3). The homotopy-series solution is given by uðx; yþ ¼u ðx; yþþ Xþ u m ðx; yþ; ðx; yþ 2X: ðþ m¼ u m ðx; yþ ðm P Þ are governed by the so-called mth-order deformation equations L½u m ðx; yþ v m u m ðx; yþš ¼ c d m ðx; yþ; ðx; yþ 2X; ð2þ with boundary conditions u m ðx; yþ ¼; ðx; yþ 2C ; m ðx; ¼ ; ðx; yþ 2C 2 ; ð4þ
3 Z. Lin, S. Liao / Commun Nonlinear Sci Numer Simulat 6 (2) where d m ðx; yþ m N½Uðx; y; qþš ðm m q¼ ð5þ and ; m 6 ; v m ¼ ; otherwise: ð6þ According to Eqs. (6) and (5), we can obtain d m ðx; yþ ¼r 2 u m ðx; yþþ Xm u j ðx; yþr 2 u m j ðx; yþþ Xm ru j ðx; yþru m j ðx; yþ ð v m Þf ðx; yþ: j¼ j¼ ð7þ Therefore, the original nonlinear problem governed by Eqs. (5), (2) and (3) has been transferred into an infinite number of linear sub-problems governed by Eqs. (2) (4) Implementation of SBFEM It should be emphasized that, using the HAM, one has great freedom to choose the auxiliary linear operators and an initial guess under some basic rules [5,7]. Here, we select r 2 as the auxiliary linear operator. Then Eq. (2) can be rewritten as where r 2 u m ðx; yþ ¼R m ðx; yþ; ðx; yþ 2X; ð8þ R m ðx; yþ ¼v m r 2 u m ðx; yþþc d m ðx; yþ: ð9þ That is, the mth-order deformation Eq. (2) becomes a 2D linear Poisson equation and all the linear sub-problems become 2D linear Poisson-type problems which can be solved by the scaled boundary finite-element method. To describe the implementation of the SBFEM, the 2D linear Poisson equation is written as r 2 uðx; yþ ¼Rðx; yþ in X ð2þ with boundary conditions u ¼ u on C u ¼ u n on C m ; ð22þ where Rðx; yþ is the so-called non-homogeneous term and the overbar denotes a prescribed value. The weak form of the weighted-residual statement is formulated as Z Z Z r T wrudx þ wrdx wu n dc ¼ ; ð23þ X X C m where w is the weighting function, and it is the starting point of the derivation of the SBFEM. In the frame of the SBFEM, a scaled boundary coordinate system is introduced at first. As shown in Fig., the circumferential coordinate s is anticlockwise along the boundary and the dimensionless radial coordinate n is a scaling factor, defined as at the scaling centre O and at the boundary S. The two straight sections s ¼ s and s ¼ s are termed side-faces. If the boundary is closed, the side-faces will coincide. Thus, the whole solution domain X is in the range of n 6 n 6 n and s 6 s 6 s. Bounded domains containing the scaling centre are modeled by taking n ¼ and n ¼ ; whereas, unbounded domains are modeled by taking n ¼ and n ¼. The mapping between the Cartesian coordinate system and this scaled boundary coordinate system is expressed by the scaling equations x ¼ x þ n x s ðsþ; y ¼ y þ n y s ðsþ; ð24þ where ðx ; y Þ is the scaling center and x s ðsþ and y s ðsþ are C continuous piecewise smooth functions describing the boundary S in Cartesian coordinates relative to the scaling center. Using conventional techniques the operator r can be mapped to the scaled boundary coordinate system as [3] r ¼ þ n ; ð25þ where b ðsþ and b 2 ðsþ are dependent only on the definition of the boundary S:
4 66 Z. Lin, S. Liao / Commun Nonlinear Sci Numer Simulat 6 (2) b ðsþ ¼ jjj ( ) y s ðsþ ;s x s ðsþ ;s ; b 2 ðsþ ¼ jjj Fig.. The scaled boundary coordinate definition. y s ðsþ x s ðsþ ; ð26þ and jjj is the Jacobian at the boundary jjj ¼x s ðsþy s ðsþ ;s y s ðsþx s ðsþ ;s : In the SBFEM, the unknown variables u at any point ðn; sþ in the domain are sought in the form uðn; sþ ¼NðsÞaðnÞ; ð27þ ð28þ where NðsÞ is the shape function, the vector aðnþ is analogous to the nodal values in the standard finite-element method. The weighting function w can be chosen as the same shape function as (28) by Galerkin approach wðn; sþ ¼NðsÞwðnÞ ¼wðnÞ T NðsÞ T : ð29þ Substituting Eqs. (25), (28) and (29) into Eq. (23) and integrating the domain integrals containing wðnþ ;n by parts with respect to n using Green s theorem, and noting that dx ¼jJjndnds, yields Z qðn Þ¼ NðsÞ T u n ðn ; sþn ds; ð3þ SZ qðn Þ¼ NðsÞ T u n ðn ; sþn ds; ð3þ S E n 2 aðnþ ;nn þðe þ E T E ÞnaðnÞ ;n E 2 aðnþþfðnþ ¼; ð32þ where and qðnþ ¼E naðnþ ;n þ E T aðnþ; FðnÞ ¼nF s ðnþ n 2 F b ðnþ; B ðsþ ¼b ðsþnðsþ; B 2 ðsþ ¼b 2 ðsþnðsþ ;s ; Z E ¼ B ðsþ T B ðsþjjjds; ZS E ¼ B 2 ðsþ T B ðsþjjjds; ZS E 2 ¼ B 2 ðsþ T B 2 ðsþjjjds; S Z F b ðnþ ¼ NðsÞ T RjJjds; S F s ðnþ ¼Nðs Þ T u n ðn; s ÞjJðs Þj þ Nðs Þ T u n ðn; s ÞjJðs Þj: ð33þ ð34þ ð35þ ð36þ ð37þ ð38þ ð39þ ð4þ ð4þ Eq. (32) is the so-called scaled boundary finite-element equation. By introducing the shape function, the Poisson equation has been weakened in the circumferential direction in a finite element manner, but remains strong in the radial direction. Thus, the governing partial differential equation is transformed to an ordinary matrix differential equation in the radial direction. The rank of matrices E ; E ; E 2 and vectors aðnþ; FðnÞ; F b ðnþ; F s ðnþ is N (where N is the number of nodes in the boundary S). The solution technique of the ordinary matrix differential equation is described in Appendix A.
5 Z. Lin, S. Liao / Commun Nonlinear Sci Numer Simulat 6 (2) Thus, when a proper initial guess u is chosen, we can gradually solve the mth-order deformation equations by means of SBFEM as m increases. It must be noted that the right-hand side term d m of the mth-order deformation Eq. (2) becomes more and more complicated as m increases. It means that the non-homogeneous term FðnÞ will become more and more complicated. This will result in a significant reduction in computational efficiency. To overcome this drawback, Chebyshev interpolation is employed to obtain a polynomial approximation P n ðnþ of degree n to FðnÞ, where n is the highest order of Chebyshev polynomials. Replacing FðnÞ with P n ðnþ, we can solve the scaled boundary finite element equations with high efficiency. For details of Chebyshev interpolation, please refer to [2] Optimal c The Mth-order approximation of uðx; yþ is given by ~u M ðx; yþ ¼u ðx; yþþ XM m¼ u m ðx; yþ: Note that, ~u M ðx; yþ is not only dependent upon x and y, but also upon the parameter c. Recently, some researches [9,2] have been done to find out the optimal convergence-control parameter c to get a faster convergent series solution. At the Mth-order approximation, the exact square residual error is defined as follow: Z Z n Z s Z n Z s D M ¼ ðn½~u M ðx; yþšþ 2 dx ¼ ðn½~u M ðn; sþšþ 2 njjjdsdn ¼ Ef ðn; sþdsdn; ð43þ s where X n s n ð42þ Ef ðn; sþ ¼ðN½~u M ðn; sþšþ 2 njjj: ð44þ However, a large amount of CPU time is needed to calculate the exact residual error D M, especially for high order of approximation. Thus, to decrease the CPU time, the so-called averaged residual error is defined by E M ¼ X I I J i¼ X J j¼ Ef ðn þ idn; s þ jdsþ; ð45þ where Dn ¼ðn n Þ=I and Ds ¼ðs s Þ=J. E M is the function of the unknown convergence-control parameter c. Obviously, the more quickly E M decreases to zero, the faster the corresponding homotopy-series solution converges. For a given M, the corresponding optimal values of the convergence-control parameter c are given by the minimum of E M, corresponding to an algebraic equation de M dc ¼ : ð46þ 3. Numerical examples To demonstrate the validity of the homotopy-based scaled boundary FEM, assume that the solution of Eq. (5) is exactly equal to a prescribed function uðx; yþ. And then the non-homogeneous term f ðx; yþ can be determined by substituting the exact solution into Eq. (5). Meanwhile, the functions g ðx; yþ and g 2 ðx; yþ can be set according to uðx; yþ. In this way, we can compare our numerical results with the known exact solution. All of our calculations were done on a laptop PC with a 2.49 GHz CPU and 2 GB RAM. 3.. For a bounded domain First of all, we consider the nonlinear 2D Poisson equation (5) defined in a unit bounded domain X ¼fðx; yþ : x 2 þ y 2 6 g with the Dirichlet-type boundary condition uðx; yþ ¼gðx; yþ on its boundary C. Let uðx; yþ ¼x 2 þ y 2 be the solution of Eq. (5), then we have f ðx; yþ ¼4 þ 8ðx 2 þ y 2 Þ and gðx; yþ ¼. Due to the symmetry, only one quarter of this problem is modeled, as illustrated in Fig. 2. The boundary condition remains the same as uðx; yþ ¼g ðx; yþ ¼onC, but the boundary condition u ;n ¼ g 2 ðx; yþ ¼ is enforced on C 2 due to the symmetry, where the comma in the subscript p designates partial derivative with respect to the variable following the comma. The initial guess is chosen as u ðx; yþ ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi x 2 þ y 2, which satisfies the boundary condition automatically. The scaling centre is chosen at the original of the Cartesian coordinates, and then the boundary C 2 is the so-called side-faces which need not to be discretized. The boundary C is discretized with three-noded quadratic elements. For a quarter unite circle, we have x s ðsþ ¼cosðsÞ; y s ðsþ ¼sinðsÞ; ð47þ where s is in the range of ½; p=2š. One kind of mesh is constructed by binary subdivision of the boundary C, i.e. the mesh consists of two scaled boundary elements, as shown in Fig. 2. To calculate the averaged residual error E M (45), we use I ¼ 2 and J ¼ 2.
6 68 Z. Lin, S. Liao / Commun Nonlinear Sci Numer Simulat 6 (2) Fig. 2. Bounded domain. pffiffiffiffiffiffi Firstly, let n, the highest order of Chebyshev polynomials, equal to. The curves of logð E MÞ versus c at different order of approximation M ¼ 4; 8; 2; 6 and 2 are shown in Fig. 3, which indicate that the optimal value of c is about -.6. In case of n ¼ 2, the optimal value of c is close to.63, as shown in Table. For different values of n, the optimal value of c changes very little. As long as the choice of c is near the optimal value, the homotopy-series solution could converge faster. For the sake of comparison, c is set as :6. The variations of the square root of E m and the CPU time with the order of approximation M at different n ¼ ; 5; 2; 25 are plotted in Figs. 4 and 5, respectively. Fig. 4 illustrates that, for a given n, the value of E M decreases as M increases, but when M increases to a certain value M p ; E M is almost the same. The lowest value of E M is denoted as E M, which depends on n. The larger value of n, the smaller the value of E M.InFig. 5, for a given order of approximation M, more CPU time is needed as n increases. It means that the highest order of Chebyshev polynomials affect not only the accuracy but also the computational efficiency. p As shown in Figs. 4 and 5, in the case of n ¼ 25 and M ¼, the CPU time is about 9236 seconds with ffiffiffiffiffiffi E M ¼ 3:6 7 ; pffiffiffiffiffiffi in the case of n ¼ and M ¼ 2; E M has decreased to : 5, and the CPU time is about 37 s, less than 3 min. In the frame of the SBFEM, properties of the standard finite-element method are maintained in the circumferential direction, but the analytical solution is provided in the radial direction. As long as the high precision is ensured in the radial direction and the number of surface elements is large enough in the circumferential direction, we could guarantee the accuracy of the numerical solution in the whole domain. In this example, since the boundary S is a quarter of circle and the exact solution on any path from the scaling centre to the boundary is the same, the impact from the discretization of the boundary can be ignored. So we only need to compare the numerical approximations with the exact solution in a given radial path. As shown in Figs. 6 and 7, the higher the order of approximation, the less the absolute error between the numerical approximation and the exact solution n= ) log ( EM M=4,8,2,6, C pffiffiffiffiffiffi Fig. 3. The logarithm of E M versus c in the case of n ¼ for the bounded domain.
7 Z. Lin, S. Liao / Commun Nonlinear Sci Numer Simulat 6 (2) Table Optimal value of c in the case of n ¼ 2 for the bounded domain. Order of approximation M Optimal value of c pffiffiffiffiffiffi Minimum value of E M E E E E E E E E E 6-2 C = n= EM -8 - n= n=2-6 n= M Fig. 4. The square root of E M versus M in the case of c ¼ :6 for the bounded domain. 4 CPU time (sec.) 3 2 n=,5,2,25 C = M Fig. 5. The CPU time versus M in the case of c ¼ :6 for the bounded domain.
8 7 Z. Lin, S. Liao / Commun Nonlinear Sci Numer Simulat 6 (2) n=2 c = -.6 u M ~.8.6 M= (initialguess) M= M=2 M=6 exactsolution Fig. 6. The approximation ~u M versus n in the case of n ¼ 2 and c ¼ :6 for the bounded domain. ξ.7 7E-6.6 n = 2 c = -.6 6E-6 n=2 c = E-6 Err.4 Err 4E-6.3 M=4 3E-6 M=6.2. M=8 M= ξ 2E-6 E-6 M=2 M= ξ Fig. 7. The error (j~u M u exact j) versus n in the case of n ¼ 2 and c ¼ :6 for the bounded domain For an unbounded domain Secondly, the nonlinear 2D Poisson equation (5) defined in an unbounded domain X ¼fðx; yþ : x 2 þ y 2 P g with the Dirichlet-type boundary is considered. Let the exact solution uðx; yþ equal to =ðx 2 þ y 2 Þ, then we have f ðx; yþ ¼4ðx 2 þ y 2 þ 2Þ=ðx 2 þ y 2 Þ 3 and gðx; yþ ¼gðx; yþ ¼. Similarly, due to the symmetry, only one quarter of this problem is modeled, as illustrated in Fig. 8. The boundary condition remains the same as uðx; yþ ¼g ðx; yþ ¼onC, and the condition u ;n ¼ g 2 ðx; yþ ¼ is enforced on C 2 due to the sym- Fig. 8. Unbounded domain.
9 Z. Lin, S. Liao / Commun Nonlinear Sci Numer Simulat 6 (2) p metry. The initial guess is chosen as u ðx; yþ ¼= ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi x 2 þ y 2. Since the boundary S is a quarter of unit circle, the scaled boundary transformation is the same as the bounded domain, but the scaling factor n varies from to þ. In order to calculate E M, we use Dn ¼ :, I ¼ and J ¼ 2 in Eq. (45). So the optimal convergence-control parameter c can be calculated by Eq. (46). As shown in Table 2, in the case of n ¼, the optimal value of c is about.59. Fig. 9 shows that if c is chosen near the optimal value, E M decreases faster as the order of approximation M increases. The value of the highest order of Chebyshev polynomials has a great influence on the computational efficiency and accuracy, as shown in Figs. and. In the case of n ¼ 2 and c ¼ :6, as M increases to 2, the square root of the averaged residual error E M decreases to 8:6 7, but the total CPU time is only about 42 s. Similarly, let n equal to 2 and the convergence-control parameter c equal to.6, then we compare the numerical approximations with the exact solution in a given radial path. The exact solution uðn; sþ is equal to =n 2, which tends to zero as n!. Figs. 2 and 3 illustrate that the absolute error between the numerical approximation and the exact solution decreases quickly as M increases. 4. Conclusions A new approach has been proposed for 2D nonlinear Poisson-type problems by combining the scaled boundary finite-element method with the homotopy analysis method. In this approach, the original nonlinear problem is transferred into an infinite number of linear Poisson-type sub-problems by the HAM, and then all those linear Poisson-type sub-problems are solved by means of the SBFEM. Numerical examples governed by a 2D nonlinear Poisson equation, including a bounded Table 2 Optimal value of c in the case of n ¼ for the unbounded domain. Order of approximation M Optimal value of c pffiffiffiffiffiffi Minimum value of E M E E E E E E E E E n= log ( EM ) M=4,8,2,6, C pffiffiffiffiffiffi Fig. 9. The logarithm of E M versus c in the case of n ¼ for the unbounded domain.
10 72 Z. Lin, S. Liao / Commun Nonlinear Sci Numer Simulat 6 (2) C = n= E M - n= n=2-7 n= M Fig.. The square root of E M versus M in the case of c ¼ :6 for the unbounded domain. 4 3 n=,5,2,25 CPU time (sec.) 2 C = M Fig.. The CPU time versus M in the case of c ¼ :6 for the unbounded domain. domain and an unbounded domain, are solved by this new approach, respectively. All the numerical results demonstrate the high efficiency and accuracy of our approach. Note that the HAM enjoys great freedom in choosing initial guess u. If the convergence-control parameter c is selected properly, the Mth-order approximation ~u M is better than the free selected initial approximation. So, it can provide us with a family of iterative formulae. Clearly, this new approach has great potential to the nonlinear problems. In this paper, we have successfully extended the SBFEM to solve nonlinear problems. Meanwhile, the idea of the HAM has been further applied to numerical calculation. Though more research effort is needed to develop and improve this new approach, it may not be far away to apply our approach to solve more nonlinear engineering problems with increased complexity.
11 Z. Lin, S. Liao / Commun Nonlinear Sci Numer Simulat 6 (2) n =2 c = M= (initialguess) M= M=2 M=6 exactsolution ~ u M ξ Fig. 2. The approximation ~u M versus n in the case of n ¼ 2 and c ¼ :6 for the unbounded domain..7 7E-6.6 n= 2 c =-.6 6E-6 n=2 c = E-6 Err.4 M=4 Err 4E-6 M=6.3 3E-6.2 2E-6 M=8. M=8 M=6 2 3 ξ E-6 M=2 2 3 ξ Fig. 3. The error (j~u M u exact j) versus n in the case of n ¼ 2 and c ¼ :6 for the unbounded domain. Acknowledgements The work is partly supported by National Natural Science Foundation of China (Approval No ) and State Key Lab of Ocean Engineering (Approval No. GKZD2). Thanks to Prof. Tao (Newcastle University) for the suggestive discussions and valuable comments. Appendix A. Solution technique of SBFEM By introducing the variable XðnÞ ¼ aðnþ ; ða:þ qðnþ Eq. (32) is transformed into the first-order ordinary differential equations nxðnþ ;n ¼ ZXðnÞ FðnÞ ða:2þ with a Hamiltonian matrix " # Z ¼ E ET E : ða:3þ E 2 þ E E ET E E
12 74 Z. Lin, S. Liao / Commun Nonlinear Sci Numer Simulat 6 (2) The eigenvalue problem is formulated as ZU ¼ UK ða:4þ with the eigenvalues of the Hamiltonian matrix K ¼ K ; K ¼dk i c ða:5þ K and the eigenvector matrix partitioned conformably U ¼ U U 2 : ða:6þ U 2 U 22 The real parts of all elements k i are non-positive. The solution of Eq. (A.2) equals XðnÞ ¼ U " # U 2 n K C ðnþ ; ða:7þ U 2 U 22 n K C 2 ðnþ where Z n C ðnþ ¼c C 2 ðnþ ¼c 2 Z n s s K A 2 FðsÞds; s s K A 22 FðsÞds; c ; c are integration constants and A ¼ A A 2 ¼ U : A 2 A 22 ða:8þ ða:9þ ða:þ Substituting (A.7) into (A.) yields aðnþ ¼U n K C ðnþþu 2 n K C 2 ðnþ qðnþ ¼U 2 n K C ðnþþu 22 n K C 2 ðnþ ða:þ ða:2þ For a bounded domain, aðnþ at n ¼ must remain finite, leading to C 2 ðnþj n¼ ¼, then c 2 ¼ Z s s K A 22 FðsÞds: ða:3þ Substituting (A.3) into (A.9) yields Z n C 2 ðnþ ¼ s s K A 22 FðsÞds: And then substituting Eqs. (A.8) and (A.4) into Eqs. (A.) and (A.2) results in Z Z n aðnþ ¼U n K c þ U 2 n K n s s K A 2 FðsÞds s s K A 22 FðsÞds ða:4þ ða:5þ and qðnþ ¼U 2 n K c þ Z s s K A 2 FðsÞds U 22 n K Z n n s s K A 22 FðsÞds: ða:6þ Note that c is determined by the boundary condition aðnþj n¼ ¼ U c U 2 Z s s K A 22 FðsÞds ða:7þ or qðnþj n¼ ¼ U 2 c U 22 Z s s K A 22 FðsÞds: ða:8þ For an unbounded domain, aðnþ as n!must remain finite, leading to C ðnþj n! ¼. The solution procedure is similar, please refer to [3].
13 Z. Lin, S. Liao / Commun Nonlinear Sci Numer Simulat 6 (2) References [] Wolf JP, Song C. The scaled boundary finite-element method a primer: derivation. Comput Struct 2;78:9 2. [2] Wolf JP, Song C. The scaled boundary finite-element method a fundamental solution-less boundary-element method. Comput Meth Appl Mech Eng 2;9: [3] Wolf JP. The scaled boundary finite-element method. Chichester (UK): Wiley; 23. [4] Song C, Wolf JP. The scaled boundary finite-element method alias consistent infinitesimal finite-element cell method for elastodynamics. Comput Meth Appl Mech Eng 997;47: [5] Deeks AJ, Cheng L. Potential flow around obstacles using the scaled boundary finite-element method. Int J Numer Meth Fluids 23;4:72 4. [6] Tao L, Song H, Chakrabarti S. Scaled boundary FEM solution of short-crested wave diffraction by a vertical cylinder. Comput Meth Appl Mech Eng 27;97( 4): [7] Song H, Tao L, Chakrabarti S. Modelling of water wave interaction with multiple cylinders of arbitrary shape. J Comput Phys 2;229(5): [8] Song H, Tao L. An efficient scaled boundary FEM model for wave interaction with a nonuniform porous cylinder. Int J Numer Meth Fluids. doi:.2/ fld.28. [9] Doherty JP, Deeks AJ. Adaptive coupling of the finite-element and scaled boundary finite-element methods for non-linear analysis of unbounded media. Comput Geotech 25;32: [] Deeks AJ, Augarde CE. A hybrid meshless Petrov Galerkin method for unbounded domains. Comput Meth Appl Mech Eng 27;96: [] Chidgzey SR, Trevelyan J, Deeks AJ. Coupling of the boundary element method and the scaled boundary finite element method for computations in fracture mechanics. Comput Struct 28;86: [2] Liao SJ. The proposed homotopy analysis technique for the solution of nonlinear problems. Ph.D. Thesis, Shanghai Jiao Tong University; 992. [3] Liao SJ. A kind of approximate solution technique which does not depend upon small parameters (II): an application in fluid mechanics. Int J Nonlinear Mech 997;32: [4] Liao SJ. An explicit, totally analytic approximation of Blasius viscous flow problems. Int J Nonlinear Mech 999;34(4): [5] Liao SJ. Beyond perturbation: introduction to the homotopy analysis method. Boca Raton (FL): Chapman & Hall; 23. [6] Liao SJ. On the homotopy analysis method for nonlinear problems. Appl Math Comput 24;47: [7] Liao SJ, Tan Y. A general approach to obtain series solutions of nonlinear differential equations. Stud Appl Math 27;9: [8] Liao SJ. Notes on the homotopy analysis method: some definitions and theorems. Commun Nonlinear Sci Numer Simulat 29;4: [9] Liao SJ. An optimal homotopy-analysis approach for strongly nonlinear differential equations. Commun Nonlinear Sci Numer Simulat. doi:.6/ j.cnsns [2] Mason JC, Handscomb DC. Chebyshev polynomials. Boca Raton (FL): Chapman & Hall; 23. [2] Niu Z, Wang Ch. A one-step optimal homotopy analysis method for nonlinear differential equations. Commun Nonlinear Sci Numer Simulat. doi:.6/j.cnsns
Commun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat xxx (2009) xxx xxx Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns A one-step optimal homotopy
More informationNewton-homotopy analysis method for nonlinear equations
Applied Mathematics and Computation 188 (2007) 1794 1800 www.elsevier.com/locate/amc Newton-homotopy analysis method for nonlinear equations S. Abbasbandy a, *, Y. Tan b, S.J. Liao b a Department of Mathematics,
More informationSemi-analytical solution of Poisson's equation in bounded domain
Semi-analytical solution of Poisson's equation in bounded domain Author Song, Hao, Tao, L. Published 21 Journal Title ANZIAM Journal Copyright Statement 21 Australian Mathematical Society. The attached
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 16 (2011) 2730 2736 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns Homotopy analysis method
More informationarxiv: v2 [math.na] 10 Nov 2015
On the method of directly defining inverse mapping for nonlinear differential equations arxiv:1507.03755v2 [math.na] 10 Nov 2015 Shijun Liao a,b,c, Yinlong Zhao a a State Key Laboratory of Ocean Engineering
More informationAn analytic approach to solve multiple solutions of a strongly nonlinear problem
Applied Mathematics and Computation 169 (2005) 854 865 www.elsevier.com/locate/amc An analytic approach to solve multiple solutions of a strongly nonlinear problem Shuicai Li, Shi-Jun Liao * School of
More informationComparison of homotopy analysis method and homotopy perturbation method through an evolution equation
Comparison of homotopy analysis method and homotopy perturbation method through an evolution equation Songxin Liang, David J. Jeffrey Department of Applied Mathematics, University of Western Ontario, London,
More informationA Numerical Study of One-Dimensional Hyperbolic Telegraph Equation
Journal of Mathematics and System Science 7 (2017) 62-72 doi: 10.17265/2159-5291/2017.02.003 D DAVID PUBLISHING A Numerical Study of One-Dimensional Hyperbolic Telegraph Equation Shaheed N. Huseen Thi-Qar
More informationSchur decomposition in the scaled boundary finite element method in elastostatics
IOP Conference Series: Materials Science and Engineering Schur decomposition in the scaled boundary finite element method in elastostatics o cite this article: M Li et al 010 IOP Conf. Ser.: Mater. Sci.
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 15 (2010) 2003 2016 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns An optimal homotopy-analysis
More informationA MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC PROBLEMS IN SEMI-INFINITE MEDIA
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 A MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC
More informationSeries solutions of non-linear Riccati differential equations with fractional order
Available online at www.sciencedirect.com Chaos, Solitons and Fractals 40 (2009) 1 9 www.elsevier.com/locate/chaos Series solutions of non-linear Riccati differential equations with fractional order Jie
More informationAnalytical Solution of BVPs for Fourth-order Integro-differential Equations by Using Homotopy Analysis Method
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.9(21) No.4,pp.414-421 Analytical Solution of BVPs for Fourth-order Integro-differential Equations by Using Homotopy
More informationImproving homotopy analysis method for system of nonlinear algebraic equations
Journal of Advanced Research in Applied Mathematics Vol., Issue. 4, 010, pp. -30 Online ISSN: 194-9649 Improving homotopy analysis method for system of nonlinear algebraic equations M.M. Hosseini, S.M.
More informationImproving convergence of incremental harmonic balance method using homotopy analysis method
Acta Mech Sin (2009) 25:707 712 DOI 10.1007/s10409-009-0256-4 RESEARCH PAPER Improving convergence of incremental harmonic balance method using homotopy analysis method Yanmao Chen Jike Liu Received: 10
More informationHomotopy Analysis Method to Water Quality. Model in a Uniform Channel
Applied Mathematical Sciences, Vol. 7, 201, no. 22, 1057-1066 HIKARI Ltd, www.m-hikari.com Homotopy Analysis Method to Water Quality Model in a Uniform Channel S. Padma Department of Mathematics School
More informationShijun Liao. Homotopy Analysis Method in Nonlinear Differential Equations
Shijun Liao Homotopy Analysis Method in Nonlinear Differential Equations Shijun Liao Homotopy Analysis Method in Nonlinear Differential Equations With 127 figures Author Shijun Liao Shanghai Jiao Tong
More informationA new method for homoclinic solutions of ordinary differential equations
Chaos, Solitons and Fractals xxx (2007) xxx xxx www.elsevier.com/locate/chaos A new method for homoclinic solutions of ordinary differential equations F. Talay Akyildiz a, K. Vajravelu b, *, S.-J. Liao
More informationJournal of Engineering Science and Technology Review 2 (1) (2009) Research Article
Journal of Engineering Science and Technology Review 2 (1) (2009) 118-122 Research Article JOURNAL OF Engineering Science and Technology Review www.jestr.org Thin film flow of non-newtonian fluids on a
More informationApproximate Analytical Solutions of Two. Dimensional Transient Heat Conduction Equations
Applied Mathematical Sciences Vol. 6 2012 no. 71 3507-3518 Approximate Analytical Solutions of Two Dimensional Transient Heat Conduction Equations M. Mahalakshmi Department of Mathematics School of Humanities
More informationPeriodic wave solution of a second order nonlinear ordinary differential equation by Homotopy analysis method
Periodic wave solution of a second order nonlinear ordinary differential equation by Homotopy analysis method Author Song, Hao, Tao, L. Published 2010 Journal Title ANZIAM Journal Copyright Statement 2010
More informationAN AUTOMATIC SCHEME ON THE HOMOTOPY ANALYSIS METHOD FOR SOLVING NONLINEAR ALGEBRAIC EQUATIONS. Safwan Al-Shara
italian journal of pure and applied mathematics n. 37 2017 (5 14) 5 AN AUTOMATIC SCHEME ON THE HOMOTOPY ANALYSIS METHOD FOR SOLVING NONLINEAR ALGEBRAIC EQUATIONS Safwan Al-Shara Department of Mathematics
More informationSecond-order wave diffraction by a circular cylinder using scaled boundary finite element method
IOP Conference Series: Materials Science and Engineering Second-order wave diffraction by a circular cylinder using scaled boundary finite element method To cite this article: H Song and L Tao 010 IOP
More informationBasic Ideas and Brief History of the Homotopy Analysis Method
1 Basic Ideas and Brief History of the Homotopy Analysis Method 1 Introduction Nonlinear equations are much more difficult to solve than linear ones, especially by means of analytic methods. In general,
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 15 (21) 3965 3973 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns A similarity solution in
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat xxx (9) xxx xxx Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns Short communication Simple
More informationarxiv: v1 [math.ap] 21 Apr 2016
arxiv:164.6711v1 [math.ap] 21 Apr 216 Analytic Solutions of Von Kármán Plate under Arbitrary Uniform Pressure Part (II): Equations in Integral Form Xiaoxu Zhong 3, Shijun Liao 1,2,3 1 State Key Laboratory
More informationAPPROXIMATING THE FORTH ORDER STRUM-LIOUVILLE EIGENVALUE PROBLEMS BY HOMOTOPY ANALYSIS METHOD
APPROXIMATING THE FORTH ORDER STRUM-LIOUVILLE EIGENVALUE PROBLEMS BY HOMOTOPY ANALYSIS METHOD * Nader Rafatimaleki Department of Mathematics, College of Science, Islamic Azad University, Tabriz Branch,
More informationExample 2: a system of coupled ODEs with algebraic property at infinity
Example 2: a system of coupled ODEs with algebraic property at infinity Consider a set of two coupled nonlinear differential equations 5 subject to f (η) + θ(η) f 2 = 0, (10) θ (η) = 3σf (η)θ(η), (11)
More informationA new scaled boundary finite element method using Fourier shape. functions
APCOM & ISCM -4 th December, 3, Singapore A new caled boundary finite element method uing Fourier hape Abtract function Yiqian He¹, Haitian Yang¹, and *Andrew J. Deek Department of Engineering Mechanic,
More informationChaos, Solitons and Fractals
Chaos, Solitons and Fractals 42 (2009) 3047 3052 Contents lists available at ScienceDirect Chaos, Solitons and Fractals journal homepage: www.elsevier.com/locate/chaos Solutions of the SIR models of epidemics
More informationarxiv: v2 [nlin.si] 27 Sep 2012
1 On peaked solitary waves of the Degasperis - Procesi equation Shijun Liao arxiv:1205.5396v2 [nlin.si] 27 Sep 2012 State Key Laboratory of Ocean Engineering, Dept. of Mathematics School of Naval Architecture,
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat xxx ) xxx xxx Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: wwwelseviercom/locate/cnsns A note on the use of Adomian
More informationMULTISTAGE HOMOTOPY ANALYSIS METHOD FOR SOLVING NON- LINEAR RICCATI DIFFERENTIAL EQUATIONS
MULTISTAGE HOMOTOPY ANALYSIS METHOD FOR SOLVING NON- LINEAR RICCATI DIFFERENTIAL EQUATIONS Hossein Jafari & M. A. Firoozjaee Young Researchers club, Islamic Azad University, Jouybar Branch, Jouybar, Iran
More informationAn explicit solution of the large deformation of a cantilever beam under point load at the free tip
Journal of Computational and Applied Mathematics 212 (2008) 320 330 www.elsevier.com/locate/cam An explicit solution of the large deformation of a cantilever beam under point load at the free tip Ji Wang
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 14 (2009) 3833 3843 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns Homotopy analysis solution
More informationApproximate Analytical Solution to Time-Fractional Damped Burger and Cahn-Allen Equations
Appl. Math. Inf. Sci. 7, No. 5, 1951-1956 (013) 1951 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.1785/amis/070533 Approximate Analytical Solution to Time-Fractional
More informationInternational Journal of Modern Mathematical Sciences, 2012, 3(2): International Journal of Modern Mathematical Sciences
Article International Journal of Modern Mathematical Sciences 2012 3(2): 63-76 International Journal of Modern Mathematical Sciences Journal homepage:wwwmodernscientificpresscom/journals/ijmmsaspx On Goursat
More informationHOMOTOPY ANALYSIS METHOD FOR SOLVING KDV EQUATIONS
Surveys in Mathematics and its Applications ISSN 1842-6298 (electronic), 1843-7265 (print) Volume 5 (21), 89 98 HOMOTOPY ANALYSIS METHOD FOR SOLVING KDV EQUATIONS Hossein Jafari and M. A. Firoozjaee Abstract.
More informationHomotopy Analysis Method for Nonlinear Jaulent-Miodek Equation
ISSN 746-7659, England, UK Journal of Inforation and Coputing Science Vol. 5, No.,, pp. 8-88 Hootopy Analysis Method for Nonlinear Jaulent-Miodek Equation J. Biazar, M. Eslai Departent of Matheatics, Faculty
More informationNumerical solution of hyperbolic heat conduction in thin surface layers
International Journal of Heat and Mass Transfer 50 (007) 9 www.elsevier.com/locate/ijhmt Numerical solution of hyperbolic heat conduction in thin surface layers Tzer-Ming Chen * Department of Vehicle Engineering,
More information2233. An improved homotopy analysis method with accelerated convergence for nonlinear problems
2233. An improved homotopy analysis method with accelerated convergence for nonlinear problems Hong-wei Li 1, Jun Wang 2, Li-xin Lu 3, Zhi-wei Wang 4 1, 2, 3 Department of Packaging Engineering, Jiangnan
More informationImproved Lebesgue constants on the triangle
Journal of Computational Physics 207 (2005) 625 638 www.elsevier.com/locate/jcp Improved Lebesgue constants on the triangle Wilhelm Heinrichs Universität Duisburg-Essen, Ingenieurmathematik (FB 10), Universitätsstrasse
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat (9) Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns On the use of homotopy analysis method
More informationA multiple Riccati equations rational expansion method and novel solutions of the Broer Kaup Kupershmidt system
Chaos, Solitons and Fractals 30 (006) 197 03 www.elsevier.com/locate/chaos A multiple Riccati equations rational expansion method and novel solutions of the Broer Kaup Kupershmidt system Qi Wang a,c, *,
More informationAnalytical solution for nonlinear Gas Dynamic equation by Homotopy Analysis Method
Available at http://pvau.edu/aa Appl. Appl. Math. ISSN: 932-9466 Vol. 4, Issue (June 29) pp. 49 54 (Previously, Vol. 4, No. ) Applications and Applied Matheatics: An International Journal (AAM) Analytical
More informationV. G. Gupta 1, Pramod Kumar 2. (Received 2 April 2012, accepted 10 March 2013)
ISSN 749-3889 (print, 749-3897 (online International Journal of Nonlinear Science Vol.9(205 No.2,pp.3-20 Approimate Solutions of Fractional Linear and Nonlinear Differential Equations Using Laplace Homotopy
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 7 (0) 3499 3507 Contents lists available at SciVerse ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns Application of the
More informationResearch Article On a New Reliable Algorithm
Hindawi Publishing Corporation International Journal of Differential Equations Volume 2009, Article ID 710250, 13 pages doi:10.1155/2009/710250 Research Article On a New Reliable Algorithm A. K. Alomari,
More informationAnalytical solution for determination the control parameter in the inverse parabolic equation using HAM
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 12, Issue 2 (December 2017, pp. 1072 1087 Applications and Applied Mathematics: An International Journal (AAM Analytical solution
More informationNumerical solution of Maxwell equations using local weak form meshless techniques
Journal of mathematics and computer science 13 2014), 168-185 Numerical solution of Maxwell equations using local weak form meshless techniques S. Sarabadan 1, M. Shahrezaee 1, J.A. Rad 2, K. Parand 2,*
More informationANALYTICAL SOLUTION FOR VIBRATION OF BUCKLED BEAMS
IJRRAS August ANALYTICAL SOLUTION FOR VIBRATION OF BUCKLED BEAMS A. Fereidoon, D.D. Ganji, H.D. Kaliji & M. Ghadimi,* Department of Mechanical Engineering, Faculty of Engineering, Semnan University, Iran
More informationAn Analytical Scheme for Multi-order Fractional Differential Equations
Tamsui Oxford Journal of Mathematical Sciences 26(3) (2010) 305-320 Aletheia University An Analytical Scheme for Multi-order Fractional Differential Equations H. M. Jaradat Al Al Bayt University, Jordan
More informationSolving the One Dimensional Advection Diffusion Equation Using Mixed Discrete Least Squares Meshless Method
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 2015 Paper No. 292 Solving the One Dimensional Advection Diffusion Equation
More informationA new approach to solve fuzzy system of linear equations by Homotopy perturbation method
Journal of Linear and Topological Algebra Vol. 02, No. 02, 2013, 105-115 A new approach to solve fuzzy system of linear equations by Homotopy perturbation method M. Paripour a,, J. Saeidian b and A. Sadeghi
More informationSOLUTION TO BERMAN S MODEL OF VISCOUS FLOW IN POROUS CHANNEL BY OPTIMAL HOMOTOPY ASYMPTOTIC METHOD
SOLUTION TO BERMAN S MODEL OF VISCOUS FLOW IN POROUS CHANNEL BY OPTIMAL HOMOTOPY ASYMPTOTIC METHOD Murad Ullah Khan 1*, S. Zuhra 2, M. Alam 3, R. Nawaz 4 ABSTRACT Berman developed the fourth-order nonlinear
More informationApparent thermal conductivity of periodic two-dimensional composites
Computational Materials Science 25 (2002) 329 338 www.elsevier.com/locate/commatsci Apparent thermal conductivity of periodic two-dimensional composites M. Jiang a, I. Jasiuk b, *, M. Ostoja-Starzewski
More informationDynamic Analysis of SSI Systems via a Coupled Finite-element/Scaled Boundary. Finite-element Model
APCOM & ISCM -4 th December, 23, Singapore Dynamic Analysis of SSI Systems via a Coupled Finite-element/Scaled Boundary Finite-element Model *D. Chen, C. Birk 2, and S. Dai 3 College of Civil Engineering
More informationIdentification of model parameters from elastic/elasto-plastic spherical indentation
Thomas Niederkofler a, Andreas Jäger a, Roman Lackner b a Institute for Mechanics of Materials and Structures (IMWS), Department of Civil Engineering, Vienna University of Technology, Vienna, Austria b
More informationHomotopy Analysis Transform Method for Time-fractional Schrödinger Equations
International Journal of Modern Mathematical Sciences, 2013, 7(1): 26-40 International Journal of Modern Mathematical Sciences Journal homepage:wwwmodernscientificpresscom/journals/ijmmsaspx ISSN:2166-286X
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationA new modified Halley method without second derivatives for nonlinear equation
Applied Mathematics and Computation 189 (2007) 1268 1273 www.elsevier.com/locate/amc A new modified Halley method without second derivatives for nonlinear equation Muhammad Aslam Noor *, Waseem Asghar
More informationOn the stability of the three classes of Newtonian three-body planar periodic orbits
. Article. SCIENCE CHINA Physics, Mechanics & Astronomy November 2014 Vol. 57 No. 11: 2121 2126 doi: 10.1007/s11433-014-5563-5 On the stability of the three classes of Newtonian three-body planar periodic
More informationOn the convergence of the homotopy analysis method to solve the system of partial differential equations
Journal of Linear and Topological Algebra Vol. 04, No. 0, 015, 87-100 On the convergence of the homotopy analysis method to solve the system of partial differential equations A. Fallahzadeh a, M. A. Fariborzi
More informationApplied Thermal Engineering
Applied Thermal Engineering 29 (2009) 372 379 Contents lists available at ScienceDirect Applied Thermal Engineering journal homepage: www.elsevier.com/locate/apthermeng Numerical simulations of the equilibrium
More information(Received 1 February 2012, accepted 29 June 2012) address: kamyar (K. Hosseini)
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.14(2012) No.2,pp.201-210 Homotopy Analysis Method for a Fin with Temperature Dependent Internal Heat Generation
More informationCONTINUED-FRACTION ABSORBING BOUNDARY CONDITIONS FOR THE WAVE EQUATION
Journal of Computational Acoustics, Vol. 8, No. 1 (2) 139 156 c IMACS CONTINUED-FRACTION ABSORBING BOUNDARY CONDITIONS FOR THE WAVE EQUATION MURTHY N. GUDDATI Department of Civil Engineering, North Carolina
More informationDirect method for deriving equilibrium equations in solid continuous systems
Engineering Solid Mechanics 2 (2014) 321-330 Contents lists available at GrowingScience Engineering Solid Mechanics homepage: www.growingscience.com/esm Direct method for deriving equilibrium equations
More informationRadial Growth of a Micro-Void in a Class of. Compressible Hyperelastic Cylinder. Under an Axial Pre-Strain *
dv. Theor. ppl. Mech., Vol. 5, 2012, no. 6, 257-262 Radial Growth of a Micro-Void in a Class of Compressible Hyperelastic Cylinder Under an xial Pre-Strain * Yuxia Song, Datian Niu and Xuegang Yuan College
More informationInternational Communications in Heat and Mass Transfer
International Communications in Heat and Mass Transfer 39 (12) 82 86 Contents lists available at SciVerse ScienceDirect International Communications in Heat and Mass Transfer journal homepage: www.elsevier.com/locate/ichmt
More informationAnalytic solution of fractional integro-differential equations
Annals of the University of Craiova, Mathematics and Computer Science Series Volume 38(1), 211, Pages 1 1 ISSN: 1223-6934 Analytic solution of fractional integro-differential equations Fadi Awawdeh, E.A.
More informationComparison of Optimal Homotopy Asymptotic Method with Homotopy Perturbation Method of Twelfth Order Boundary Value Problems
Abstract Comparison of Optimal Homotopy Asymptotic Method with Homotopy Perturbation Method of Twelfth Order Boundary Value Problems MukeshGrover grover.mukesh@yahoo.com Department of Mathematics G.Z.S.C.E.T
More informationAdomain Decomposition Method for Solving Non Linear Partial Differential Equations
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 10, Issue 5 Ver. V (Sep-Oct. 2014), PP 60-66 Adomain Decomposition Method for Solving Non Linear Partial Differential
More informationExplicit Analytic Solution for an. Axisymmetric Stagnation Flow and. Heat Transfer on a Moving Plate
Int. J. Contep. Math. Sciences, Vol. 5,, no. 5, 699-7 Explicit Analytic Solution for an Axisyetric Stagnation Flow and Heat Transfer on a Moving Plate Haed Shahohaadi Mechanical Engineering Departent,
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationANALYTICAL APPROXIMATE SOLUTIONS OF THE ZAKHAROV-KUZNETSOV EQUATIONS
(c) Romanian RRP 66(No. Reports in 2) Physics, 296 306 Vol. 2014 66, No. 2, P. 296 306, 2014 ANALYTICAL APPROXIMATE SOLUTIONS OF THE ZAKHAROV-KUZNETSOV EQUATIONS A. JAFARIAN 1, P. GHADERI 2, ALIREZA K.
More informationComputation Time Assessment of a Galerkin Finite Volume Method (GFVM) for Solving Time Solid Mechanics Problems under Dynamic Loads
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 215 Paper o. 31 Computation Time Assessment of a Galerkin Finite Volume
More informationScaled boundary FEM solution of short-crested wave diffraction by a vertical cylinder
Scaled boundary FEM solution of short-crested wave diffraction by a vertical cylinder Author Tao, Longbin, Song, Hao, Chakrabarti, Subrata Published 27 Journal Title Computer Methods in Applied Mechanics
More informationfield using second order edge elements in 3D
The current issue and full text archive of this journal is available at http://www.emerald-library.com using second order edge elements in 3D Z. Ren Laboratoire de GeÂnie Electrique de Paris, UniversiteÂs
More informationOn the steady-state resonant acoustic gravity waves
J. Fluid Mech. (2018), vol. 849, pp. 111 135. c Cambridge University Press 2018 doi:10.1017/jfm.2018.422 111 On the steady-state resonant acoustic gravity waves Xiaoyan Yang 3, Frederic Dias 4 and Shijun
More informationJournal of Computational Physics
Journal of Computational Physics 23 (211) 8573 862 Contents lists available at SciVerse ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp An eigen-based high-order
More informationResearch Article Series Solution of the Multispecies Lotka-Volterra Equations by Means of the Homotopy Analysis Method
Hindawi Publishing Corporation Differential Equations and Nonlinear Mechanics Volume 28, Article ID 816787, 14 pages doi:1.1155/28/816787 Research Article Series Solution of the Multispecies Lotka-Volterra
More informationBack Analysis of Measured Displacements of Tunnels
Rock Mechanics and Rock Engineering 16, 173--180 (1983) Rock Mechanics and Rock Engineering 9 by Springer-Verlag 1983 Back Analysis of Measured Displacements of Tunnels By S. Sakurai and K. Takeuchi Kobe
More informationInitial inverse problem in heat equation with Bessel operator
International Journal of Heat and Mass Transfer 45 (22) 2959 2965 www.elsevier.com/locate/ijhmt Initial inverse problem in heat equation with Bessel operator Khalid Masood *, Salim Messaoudi, F.D. Zaman
More informationGalerkin method for the numerical solution of the RLW equation using quintic B-splines
Journal of Computational and Applied Mathematics 19 (26) 532 547 www.elsevier.com/locate/cam Galerkin method for the numerical solution of the RLW equation using quintic B-splines İdris Dağ a,, Bülent
More informationHaar wavelet method for nonlinear integro-differential equations
Applied Mathematics and Computation 176 (6) 34 333 www.elsevier.com/locate/amc Haar wavelet method for nonlinear integro-differential equations Ülo Lepik Institute of Applied Mathematics, University of
More informationDynamics of four coupled phase-only oscillators
Communications in Nonlinear Science and Numerical Simulation 13 (2008) 501 507 www.elsevier.com/locate/cnsns Short communication Dynamics of four coupled phase-only oscillators Richard Rand *, Jeffrey
More informationApplication of Homotopy Analysis Method for Linear Integro-Differential Equations
International Mathematical Forum, 5, 21, no. 5, 237-249 Application of Homotopy Analysis Method for Linear Integro-Differential Equations Zulkifly Abbas a, Saeed Vahdati a,1, Fudziah Ismail a,b and A.
More informationThe Series Solution of Problems in the Calculus of Variations via the Homotopy Analysis Method
The Series Solution of Problems in the Calculus of Variations via the Homotopy Analysis Method Saeid Abbasbandy and Ahmand Shirzadi Department of Mathematics, Imam Khomeini International University, Ghazvin,
More informationJOURNAL OF INTERNATIONAL ACADEMIC RESEARCH FOR MULTIDISCIPLINARY Impact Factor 1.393, ISSN: , Volume 2, Issue 7, August 2014
HOMOTOPY ANALYSIS TO THERMAL RADIATION EFFECTS ON HEAT TRANSFER OF WALTERS LIQUID-B FLOW OVER A STRETCHING SHEET FOR LARGE PRANDTL NUMBERS HYMAVATHI TALLA* P.VIJAY KUMAR** V.MALLIPRIYA*** *Dept. of Mathematics,
More informationQuintic B-Spline Galerkin Method for Numerical Solutions of the Burgers Equation
5 1 July 24, Antalya, Turkey Dynamical Systems and Applications, Proceedings, pp. 295 39 Quintic B-Spline Galerkin Method for Numerical Solutions of the Burgers Equation İdris Dağ 1,BülentSaka 2 and Ahmet
More informationAn a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element
Calcolo manuscript No. (will be inserted by the editor) An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Dietrich Braess Faculty of Mathematics, Ruhr-University
More informationRational Energy Balance Method to Nonlinear Oscillators with Cubic Term
From the SelectedWorks of Hassan Askari 2013 Rational Energy Balance Method to Nonlinear Oscillators with Cubic Term Hassan Askari Available at: https://works.bepress.com/hassan_askari/4/ Asian-European
More informationCOMPUTATIONAL ELASTICITY
COMPUTATIONAL ELASTICITY Theory of Elasticity and Finite and Boundary Element Methods Mohammed Ameen Alpha Science International Ltd. Harrow, U.K. Contents Preface Notation vii xi PART A: THEORETICAL ELASTICITY
More informationEigenvalue inverse formulation for optimising vibratory behaviour of truss and continuous structures
Computers and Structures 8 () 397 43 www.elsevier.com/locate/compstruc Eigenvalue inverse formulation for optimising vibratory behaviour of truss and continuous structures H. Bahai *, K. Farahani, M.S.
More information1 Introduction IPICSE-2016
(06) DOI: 0.05/ matecconf/06860006 IPICSE-06 Numerical algorithm for solving of nonlinear problems of structural mechanics based on the continuation method in combination with the dynamic relaxation method
More informationSelection of Rayleigh Damping Coefficients for Seismic Response Analysis of Soil Layers
Selection of Rayleigh Damping Coefficients for Seismic Response Analysis of Soil Layers Huai-Feng Wang, Meng-Lin Lou, Ru-Lin Zhang Abstract One good analysis method in seismic response analysis is direct
More informationComparison of Homotopy-Perturbation Method and variational iteration Method to the Estimation of Electric Potential in 2D Plate With Infinite Length
Australian Journal of Basic and Applied Sciences, 4(6): 173-181, 1 ISSN 1991-8178 Comparison of Homotopy-Perturbation Method and variational iteration Method to the Estimation of Electric Potential in
More informationComputers and Mathematics with Applications. A modified variational iteration method for solving Riccati differential equations
Computers and Mathematics with Applications 6 (21) 1868 1872 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa A modified
More informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More information