Novel solution for acceleration motion of a vertically falling spherical particle by HPM–Padé approximant

In this paper the acceleration motion of a vertically falling spherical particle in incompressible Newtonian media is investigated. The velocity is evaluated by using homotopy perturbation method (HPM) and Padé approximant which is an analytical solution technique. The current results are then compared with those derived from HPM and the established fourth order Runge–Kutta method in order to verify the accuracy of the proposed method. It is found that this method can achieve more suitable results in comparison to HPM.     

Weighting parameters for time‐step integration algorithms with predetermined coefficients

In this paper, the effect of using the predetermined coefficients in constructing time‐step integration algorithms is investigated. Both first‐ and second‐order equations are considered. The approximate solution is assumed to be in a form of polynomial in the time domain. It can be related to the truncated Taylor's series expansion of the exact solution. Therefore, some of the coefficients can be predetermined from the known initial conditions. If there are m predetermined coefficients and r unknown coefficients in the approximate solution, the unknowns can be solved by the weighted residual method. The weighting parameter method is used to investigate the resultant algorithm characteristics. It is shown that the formulation is consistent with a minimum order of accuracy m+r. The maximum order of accuracy achievable is m+2r. Unconditionally stable algorithms exist if m⩽r for first‐order equations and m+1⩽r for second‐order equations. Hence, the Dahlquist's theorem is generalized. Algorithms equivalent to the Padé approximations and unconditionally stable algorithms equivalent to the generalized Padé approximations are constructed. The corresponding weighting parameters and weighting functions for the Padé and generalized Padé approximations are given explicitly. Copyright © 2000 John Wiley & Sons, Ltd.     

Weighting parameters for unconditionally stable higher‐order accurate time step integration algorithms. Part 1—first‐order equations

In this paper, unconditionally stable higher‐order accurate time step integration algorithms suitable for linear first‐order differential equations based on the weighted residual method are presented. Instead of specifying the weighting functions, the weighting parameters are used to control the algorithm characteristics. If the numerical solution is approximated by a polynomial of degree n, the approximation is at least nth‐order accurate. By choosing the weighting parameters carefully, the order of accuracy can be improved. The generalized Padé approximations with polynomials of degree n as the numerator and denominator are considered. The weighting parameters are chosen to reproduce the generalized Padé approximations. Once the weighting parameters are known, any set of linearly independent basic functions can be used to construct the corresponding weighting functions. The stabilizing weighting factions for the weighted residual method are then found explicitly. The accuracy of the particular solution due to excitation is also considered. It is shown that additional weighting parameters may be required to maintain the overall accuracy. The corresponding equations are listed and the additional weighting parameters are solved explicitly. However, it is found that some weighting functions could satisfy the listed equations automatically. Copyright © 1999 John Wiley & Sons, Ltd.     

A boundary condition in Padé series for frequency‐domain solution of wave propagation in unbounded domains

A boundary condition satisfying the radiation condition at infinity is frequently required in the numerical simulation of wave propagation in an unbounded domain. In a frequency domain analysis using finite elements, this boundary condition can be represented by the dynamic stiffness matrix of the unbounded domain defined on its boundary. A method for determining a Padé series of the dynamic stiffness matrix is proposed in this paper. This method starts from the scaled boundary finite‐element equation, which is a system of ordinary differential equations obtained by discretizing the boundary only. The coefficients of the Padé series are obtained directly from the ordinary differential equations, which are not actually solved for the dynamic stiffness matrix. The high rate of convergence of the Padé series with increasing order is demonstrated numerically. This technique is applicable to scalar waves and elastic vector waves propagating in anisotropic unbounded domains of irregular geometry. It can be combined seamlessly with standard finite elements. Copyright © 2006 John Wiley & Sons, Ltd.     

Padé approximants and the modal connection: Towards increased robustness for fast parametric sweeps

To increase the robustness of a Padé‐based approximation of parametric solutions to finite element problems, an a priori estimate of the poles is proposed. The resulting original approach is shown to allow for a straightforward, efficient, subsequent Padé‐based expansion of the solution vector components, overcoming some of the current convergence and robustness limitations. In particular, this enables for the intervals of approximation to be chosen a priori in direct connection with a given choice of Padé approximants. The choice of these approximants, as shown in the present work, is theoretically supported by the Montessus de Ballore theorem, concerning the convergence of a series of approximants with fixed denominator degrees. Key features and originality of the proposed approach are (1) a component‐wise expansion which allows to specifically target subsets of the solution field and (2) the a priori, simultaneous choice of the Padé approximants and their associated interval of convergence for an effective and more robust approximation. An academic acoustic case study, a structural‐acoustic application, and a larger acoustic problem are presented to demonstrate the potential of the approach proposed.     

Matrix‐Padé via Lanczos solutions for vibrations of fluid–structure interaction

For multiple‐frequency full‐field solutions of the boundary value problem describing small fluid–structure interaction vibration superimposed on a nominal state with prestress, we propose an efficient reduced order method by constructing the full‐field matrix‐Padé approximant of its finite element matrix function. Exploiting the matrix‐Padé via Lanczos connection, the Padé coefficients are computed in a stable and efficient way via an unsymmetric, banded Lanczos process. The full‐field Padé‐type approximant is the result of one‐sided projection onto Krylov subspace, we established its order of accuracy, which is not maximal. The superiority of this method in terms of various problem dimensions and parameters is established by complexity analysis via flop counts. Numerical examples obtained by using a model problem verified the accuracy of this full‐field matrix‐Padé approximant. Copyright © 2010 John Wiley & Sons, Ltd.     

Adaptive frequency windowing for multifrequency solutions in structural acoustics based on the matrix Padé‐via‐Lanczos algorithm

For many problems in structural acoustics, it is desired to obtain solutions at many frequencies over a large range in the frequency domain. A reduced‐order multifrequency algorithm based on matrix Padé approximation, using the matrix Padé‐via‐Lanczos (MPVL) connection, has been previously used to solve both exterior and interior acoustic problems. However, the method is not guaranteed to give the correct solution across the entire frequency region of interest, but only locally around a reference frequency. An adaptive frequency windowing scheme is introduced to address this shortcoming for practical application of this method. The application of this algorithm to tightly coupled problems in interior structural acoustics is presented. Copyright © 2007 John Wiley & Sons, Ltd.     

A closed-form solution for wave propagation in optical waveguides of longitudinally invariant structures without using beam propagation algorithm

In this paper, we propose a different method to study wave propagation in longitudinally invariant waveguides with arbitrary index profile. In our method, both the electric field and the refractive index profile are expanded into two Fourier cosine series. With these series substituted into the wave equation, a differential matrix equation can then be obtained. We show here that such a matrix equation can be solved and an explicit expression for the wave field at any longitudinal position along an optical waveguide can be obtained. The solution proposed in this method can thus exclude the use of the beam propagation algorithm. This study demonstrates that our approach yields the same results as those obtained by using commercial softwares in which a beam propagation method with the Padé approximation is used.     

The solution of higher order integration formulae for dynamic response equations by the conjugate gradient method

Higher order implicit integration techniques for solving dynamic response equations are derived utilizing Padé approximations. In an effort to minimize the disadvantages of using these higher order formulae to obtain solutions to systems with large numbers of degrees-of-freedom, the conjugate gradient method is employed to solve for the displacements. The accuracy and efficiency of the techniques are evaluated by making comparisons between known analytical and calculated results.     

Padé extensions to ray series methods and analysis of transients in inhomogeneous viscoelastic media of hereditary type

In this paper we demonstrate that wavefront expansions for the analysis of transient phenomena are far from adequate when numerical information back of the wavefront is required. However, by employing Padé approximants together with ray series methods, we can obtain directly and greatly extend the range of validity of these expansions. The procedure, which is straightforward and requires very little computing time, is here applied to a non-trivial problem involving impact-generated shear transients in inhomogenoues viscoelastic media whose stress-strain laws are given in integral form. For a special combination of the material parameters an exact solution is recovered and used to check the validity of our approximate Padé-extended wavefront solution. We also compare our results from the extended wavefront solution with numerical solutions obtained using Bellman's approximate inversion scheme for Laplace transforms. A further advantage of our approach is that, unlike transform techniques, it does not depend upon being able to find tabulated special function equations for the transformed dependent variables. All numerical results are presented graphically for ease of comparison.     

Numerical accuracy of a Padé‐type non‐reflecting boundary condition for the finite element solution of acoustic scattering problems at high‐frequency

The present text deals with the numerical solution of two‐dimensional high‐frequency acoustic scattering problems using a new high‐order and asymptotic Padé‐type artificial boundary condition. The Padé‐type condition is easy‐to‐implement in a Galerkin least‐squares (iterative) finite element solver for arbitrarily convex‐shaped boundaries. The method accuracy is investigated for different model problems and for the scattering problem by a submarine‐shaped scatterer. As a result, relatively small computational domains, optimized according to the shape of the scatterer, can be considered while yielding accurate computations for high‐frequencies. Copyright © 2005 John Wiley & Sons, Ltd.     

Three‐dimensional magnetostatic problem

The paper is devoted to an approximation of the solution of Maxwell's equations in three‐dimensional space. We present two methods which couple a finite element method inside the magnetic materials with a boundary integral method which uses Poincaré–Steklov's operator to describe the exterior domain. A computer code has been implemented for each method and a number of numerical experiments have been performed to validate each proposed methodology. Namely, we present numerical results concerning a non‐linear magnetostatic problem in ℝ³. Copyright © 1999 John Wiley & Sons, Ltd.     

A methodology for the generation of low‐cost higher‐order methods for linear dynamics

This work presents a methodology which generates efficient higher‐order methods for linear dynamics by improving the accuracy properties of Nørsett methods towards those of Padé methods. The methodology is based on a simple and low‐cost iterative procedure which is used to implement a set of higher‐order methods with controllable dissipation. A sequence of improved solutions is obtained which correspond to algorithms offering an effective compromise between the efficiency of Nørsett methods and the accuracy of Padé methods. Moreover, a direct control over high‐frequency dissipation is possible by means of an algorithmic parameter. Numerical tests are reported which confirm that this set of algorithms is really attractive for linear dynamic analysis. Copyright © 2003 John Wiley & Sons, Ltd.     

Polynomial chaos‐based extended Padé expansion in structural dynamics

The response of a random dynamical system is totally characterized by its probability density function (pdf). However, determining a pdf by a direct approach requires a high numerical cost; similarly, surrogate models such as direct polynomial chaos expansions are not generally efficient, especially around the eigenfrequencies of the dynamical system. In the present study, a new approach based on Padé approximants to obtain moments and pdf of the dynamic response in the frequency domain is proposed. A key difference between the direct polynomial chaos representation and the Padé representation is that the Padé approach has polynomials in both numerator and denominator. For frequency response functions, the denominator plays a vital role as it contains the information related to resonance frequencies, which are uncertain. A Galerkin approach in conjunction with polynomial chaos is proposed for the Padé approximation. Another physics‐based approach, utilizing polynomial chaos expansions of the random eigenmodes, is proposed and compared with the proposed Padé approach. It is shown that both methods give accurate results even if a very low degree of the polynomial expansion is used. The methods are demonstrated for two degree‐of‐freedom system with one and two uncertain parameters. Copyright © 2016 John Wiley & Sons, Ltd.     

High‐order predictor–corrector algorithms

New predictor–corrector algorithms are presented for the computation of solution paths of non‐linear partial differential equations. The predictors and the correctors are based on perturbation techniques and Padé approximants. This extends the Asymptotic Numerical Method (ANM), which is an efficient high‐order continuation technique without corrector. The efficiency and the reliability of the new technique are assessed by several examples within thin shell theory and Navier–Stokes equations. Many variants have been tested to establish an optimal algorithm. Copyright © 2002 John Wiley & Sons, Ltd.     

Bifurcation points and bifurcated branches by an asymptotic numerical method and Padé approximants

The aim of this work is to develop a reliable and fast algorithm to compute bifurcation points and bifurcated branches. It is based upon the asymptotic numerical method (ANM) and Padé approximants. The bifurcation point is detected by analysing the poles of Padé approximants or by evaluating, along the computed solution branch, a bifurcation indicator well adapted to ANM. Several examples are presented to assess the effectiveness of the proposed method, that emanate from buckling problems of thin elastic shells. Especially problems involving large rotations are discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

Padé approximants for the numerical solution of singular integral equations

Summary A method for accelerating the convergence of the numerical solution of a singular integral equation, based on Padé Approximants, is given in this paper. At first the general form of the Padé Table and of the “epsilon” algorithm are presented. Taking into consideration the classical quadrature method, based on the Gauss-Jacobi quadrature rule, an approximate formula is derived for the unknown density function of the Cauchy-type singular integral equation or of the equivalent Fredholm integral equation. In this formula applying the “epsilon” algorithm to the solution for the stress intensity factors, the convergence is achieved after a few operations. The number of numerical operations required for the determination of stress intensity factors is considerable reduced, when compared to the number of operations required for a classical type of solution. Illustrative examples are given, indicating the efficiency of the method.     

High-order hierarchical A- and L-stable integration methods

Numerical integration of stiff first-order systems of differential equations is considered. It is shown how a C⁰-continuous polynomial time discretization, in conjunction with a weighted residual method, can be used to derive methods corresponding to the diagonal and first subdiagonal Padé approximants. In this manner A -stable schemes of order 2k and L-stable schemes of order 2k - l are obtained, at least for k ? 4. The methods are hierarchical in the sense that a scheme with a given accuracy embraces the equations of all lower-order methods that have the same stability type. We show how this feature may be utilized to perform partial refinements in the integration process, and how an error estimation by an embedding approach naturally follows. An adaptive algorithm based on the error estimator is suggested. Some numerical experiments that illustrate the ideas are included.     

An application of theN-strip method of integral relations for analyzing the flow around a circular cone

Summary The solution of sets of non-linear partial differential equations using the method of integral relations is considered. Emphasis is laid on the derivation of a generalN-strip approximation algorithm. In order to check the applicability of this algorithm a program has been written to obtain the solution of the flow field around a circular cone at incidence in supersonic flow. Using the method of Stone, the angle of attack has been taken into account up to the second order. Thus a comparison can be made with the results given by Kopal.The results show that theN-strip algorithm in the case studied is a very attractive method which leads straight-forward to results of high accuracy.Formerly National Aerospace Laboratory (NLR).     

An adaptive strategy for the bivariate solution of finite element problems using multivariate nested Padé approximants

Most engineering applications involving solutions by numerical methods are dependent on several parameters, whose impact on the solution may significantly vary from one to the other. At times an evaluation of these multivariate solutions may be required at the expense of a prohibitively high computational cost. In the present paper, an adaptive approach is proposed as a way to estimate the solution of such multivariate finite element problems. It is based upon the integration of so‐called nested Padé approximants within the finite element procedure. This procedure includes an effective control of the approximation error, which enables adaptive refinements of the converged intervals upon reconstruction of the solution. The main advantages lie in a potential reduction of the computational effort and the fact that the level of a priori knowledge required about the solution in order to have an accurate, efficient, and well‐sampled estimate of the solution is low. The approach is introduced for bivariate problems, for which it is validated on elasto‐poro‐acoustic problems of both academic and more industrial scale. It is argued that the methodology in general holds for more than two variables, and a discussion is opened about the truncation refinements required in order to generalize the results accordingly. Copyright © 2014 John Wiley & Sons, Ltd.