Numerical Solution of the Diffusion Equation with Restrictive Pade Approximation

The problem of solving the linear diffusion equation by a method related to the Restrictive Pade Approximation (RPA) is considered. The advantage is that it has the exact value at certain r. This method will exhibit several advantages for example highly accurate, fast and with good results, etc. The absolutely error is still very small. The obtained results are compared with the exact solution and the other methods. The numerical results are in agreement with the exact solution.

PADE APPROXIMANTS AS LIMITS OF RATIONAL FUNCTIONS OF BEST APPROXIMATION IN ORLICZ SPACE

In this paper, we prove that the best rational approximation of a given analytic function in Orlicz space L~*(G), where G = {|z|≤∈}, converges to the Pade approximants of the function as the measure of G approaches zero.

Normal moveout for long offset in isotropic media using the Padé approximation

The nomaal moveout correction is important to long-offset observations, especially deep layers. For isotropic media, the conventional two-term approximation of the normal moveout function assumes a small offset-to-depth ratio and thus fails at large offset-to-depth ratios. We approximate the long-offset moveout using the Pade approximation. This method is superior to typical methods and flattens the seismic gathers over a wide range of offsets in multilayered media. For a four-layer model, traditional methods show traveltime errors of about 5 ms for offset-to-depth ratio of 2 and greater than 10 ms for offset-to-depth ratio of 3; in contrast, the maximum traveltime error for the [3, 3]-order Pade approximation is no more than 5 ms at offset-to-depth ratio of 3. For the Cooper Basin model, the maximum oft＇set-to-depth ratio for the [3, 3]-order Pade approximation is typically double of those in typical methods. The [7, 7]-order Pade approximation performs better than the [3.3]-order Pade armroximation.

Magnetic interface forward and inversion method based on Padé approximation

The magnetic interface forward and inversion method is realized using the Taylor series expansion to linearize the Fourier transform of the exponential function. With a large expansion step and unbounded neighborhood, the Taylor series is not convergent, and therefore, this paper presents the magnetic interface forward and inversion method based on Pade approximation instead of the Taylor series expansion. Compared with the Taylor series, Pade＇s expansion＇s convergence is more stable and its approximation more accurate. Model tests show the validity of the magnetic forward modeling and inversion of Pade approximation proposed in the paper, and when this inversion method is applied to the measured data of the Matagami area in Canada, a stable and reasonable distribution of underground interface is obtained.

分子阻障内旋转配分函数的二点pade’逼近

一、分子阻障内旋转配分函数的计算在理论化学中一直为人们所研究,早在五十年代Pitzer就用数值方法进行了计算,但求得这一配分函数的近似解析表示式仍有重要意义,在文献中作者曾用Wigner-Kirkwood微扰展开法求得阻障内旋转配分函数的近似解析表示式,根据此式所作的数值计算结果表明:在温度较高,内旋转动惯量较大时和Pitzer的数值计算结果符合得甚好,但随着温度的降低偏离也就增大,计算Wigner-Kirkwood展开式的的高次项固然可以改善所得之配分函数,但此时计算烦杂,困难较大,本文则采用二点Pade’逼近法将Wigner-Kirkwood微扰展开式（包括一级量子改正项）作为高温展开式,并选用配分函数Q=Σe-tEi的前三项作为其低温展开式,然后用二点Pade’逼近进行内插,求得阻障内旋转配分函数的[5/5]Pade’逼近式,它在整个温度范围内均能给出良好的数值结果。

CONVERGENCE IN MEASURE OF MULTIVARIATE PADE APPROXIMANTS

Convergence conclusions of Pade approximants in the univariate case can be found in various papers. However,resuhs in the multivariate case are few.A.Cuyt seems to be the only one who discusses convergence for multivariate Pade approximants,she gives in[2]a de Montessus de Bollore type theorem.In this paper,we will discuss the zero set of a real multivariate polynomial,and present a convergence theorem in measure of multivariate Pade approximant.The proof technique used in this paper is quite different from that used in the univariate case.

Artificial boundary condition for Klein-Gordon equation by constructing mechanics structure

An innovative local artificial boundary condition is proposed to numerically solve the Cauchy problem of the Klein-Gordon equation in an unbounded domain.Initially,the equation is considered as the axial wave prop-agation in a bar supported on a spring foundation.The numerical model is then truncated by replacing the half-infinitely long bar with an equivalent mechanical structure.The effective frequency-dependent stiffness of the half-infinitely long bar is expressed as the sum of rational terms using Pade approximation.For each term,a corresponding substructure composed of dampers and masses is constructed.Finally,the equivalent mechan-ical structure is obtained by parallelly connecting these substructures.The proposed approach can be easily implemented within a standard finite element framework by incorporating additional mass points and damper elements.Numerical examples show that with just a few extra degrees of freedom,the proposed approach effec-tively suppresses artificial reflections at the truncation boundary and exhibits first-order convergence.

Delay-dependent stability and added damping of SDOF real-time dynamic hybrid testing

It is well-recognized that a transfer system response delay that reduces the test stability inevitably exists in real-time dynamic hybrid testing （RTDHT）. This paper focuses on the delay-dependent stability and added damping of SDOF systems in RTDHT. The exponential delay term is transferred into a rational fraction by the Pad6 approximation, and the delay-dependent stability conditions and instability mechanism of SDOF RTDHT systems are investigated by the root locus technique. First, the stability conditions are discussed separately for the cases of stiffness, mass, and damping experimental substructure. The use of root locus plots shows that the added damping effect and instability mechanism for mass are different from those for stiffness. For the stiffness experimental substructure case, the instability results from the inherent mode because of an obvious negative damping effect of the delay. For the mass case, the delay introduces an equivalent positive damping into the inherent mode, and instability occurs at an added high frequency mode. Then, the compound stability condition is investigated for a general case and the results show that the mass ratio may have both upper and lower limits to remain stable. Finally, a high-emulational virtual shaking table model is built to validate the stability conclusions.

A New Form of Boussineq Equation with Improved Linear Dispersion Characteristics

It has been shown that Boussinesq type equations, which include the lowest order effects of nonlinearity and frequency dispersion, can provide an accurate description of wave evolution in coastal regions. But different linear dispersion characteristics of the equation can be obtained by different integrating method. In this paper, a new form of the Boussinesq equation is derived by use of two different layer horizontal velocity variables instead of the commonly used depth-averaged velocity or an arbitrary layer velocity. This significantly improves the linear dispersion properties of the Boussinesq equation and enables it to be applied to a wider range of water depth.

Solving multi-object radar cross section based on wide-angle parabolic equation method

Based on a Pade approximation, a wide-angle parabolic equation method is introduced for computing the multiobject radar cross section （RCS） for the first time. The method is a paraxial version of the scalar wave equation, which solves the field by marching them along the paraxial direction. Numerical results show that a single wide-angle parabofic equation run can compute multi-object RCS efficiently for angles up to 45 ° . The method provides anew and efficient numerical method for computation electromagnetics.

Numerical Simulation of Modified Kortweg-de Vries Equation by Linearized Implicit Schemes

In this paper, we are going to derive four numerical methods for solving the Modified Kortweg-de Vries (MKdV) equation using fourth Pade approximation for space direction and Crank Nicolson in the time direction. Two nonlinear schemes and two linearized schemes are presented. All resulting schemes will be analyzed for accuracy and stability. The exact solution and the conserved quantities are used to highlight the efficiency and the robustness of the proposed schemes. Interaction of two and three solitons will be also conducted. The numerical results show that the interaction behavior is elastic and the conserved quantities are conserved exactly, and this is a good indication of the reliability of the schemes which we derived. A comparison with some existing is presented as well.

Robust Control of Time Delay Systems

This study presents the design of both H∞ Loop Shaping Control （HLSC） and Internal Model Control （IMC） strategies for linear time delay systems. For first order time delay system, a systematic approach for weight selection based on the sensitivity function was proposed, then compared to the internal model control strategy. For both methods, the synthesis was based on the Pade approximation. Two cases are considered for time delay： upper or lower than system time constant. Simulation results for the proposed approaches are acceptable ever in presence of disturbances and model mismatches.

Mean-field and high temperature series expansion calculations of some magnetic properties of Ising and XY antiferromagnetic thin-films

In this work, the magnetic properties of Ising and XY antiferromagnetic thin-films are investigated each as a function of Neel temperature and thickness for layers （n = 2, 3, 4, 5, 6, and bulk （∞） by means of a mean-field and high temperature series expansion （HTSE） combined with Pade approximant calculations. The scaling law of magnetic susceptibility and magnetization is used to determine the critical exponent γ, veff （mean）, ratio of the critical exponents γ/v, and magnetic properties of Ising and XY antiferromagnetic thin-films for different thickness layers n = 2, 3, 4, 5, 6, and bulk （∞）.

The construction of homoclinic and heteroclinic orbitals in asymmetric strongly nonlinear systems based on the Pad'e approximant

In this paper, the extended Pade approximant is used to construct the homoclinic and the heteroclinic trajectories in nonlinear dynamical systems that are asymmetric at origin. Meanwhile, the conservative system, the autonomous system, and the nonautonomous system equations with quadratic and cubic nonlinearities are considered. The disturbance parameter ~ is not limited to being small. The ranges of the values of the linear and the nonlinear term parameters, which are variables, can be determined when the boundary values are satisfied. New conditions for the potentiality and the convergence are posed to make it possible to solve the boundary-value problems formulated for the orbitals and to evaluate the initial amplitude values.

Adomian decomposition method and Padé approximants for solving the Blaszak-Marciniak lattice

The Adomian decomposition method （ADM） and Pade approximants are combined to solve the well-known Blaszak-Marciniak lattice, which has rich mathematical structures and many important applications in physics and mathematics. In some cases, the truncated series solution of ADM is adequate only in a small region when the exact solution is not reached. To overcome the drawback, the Pade approximants, which have the advantage in turning the polynomials approximation into a rational function, are applied to the series solution to improve the accuracy and enlarge the convergence domain. By using the ADM-Pade technique, the soliton solutions of the Blaszak-Marciniak lattice are constructed with better accuracy and better convergence than by using the ADM alone. Numerical and figurative illustrations show that it is a promising tool for solving nonlinear problems.

Adomian Decomposition Method and Padé Approximants for Nonlinear Differential-Difference Equations

Combining Adomian decomposition method （ADM） with Pade approximants, we solve two differentiaidifference equations （DDEs）： the relativistic Toda lattice equation and the modified Volterra lattice equation. With the help of symbolic computation Maple, the results obtained by ADM-Pade technique are compared with those obtained by using ADM alone. The numerical results demonstrate that ADM-Pade technique give the approximate solution with faster convergence rate and higher accuracy and relative in larger domain of convergence than using ADM.

Brightness Temperature Model of Sea Foam Layer at L-band

Permittivity of a sea foam layer is very important in investigating ocean brightness temperature model. At microwave frequency, the Rayleigh method is developed to estimate the effective permittivity of the sea foam layer. To simplify the tedious calculation of sea foam effective permittivity at L band （1.4GHz）, Pade＇ approximation is adopted to fit the sea foam effective permittivity computed by the Rayleigh method. With this fitting formula, a new brightness temperature model of sea foam layer defined by certain geophysical parameters, such as air volume fraction （AVF）, sea surface temperature （SST）, sea surface salinity （SSS） and thickness of foam layer d, is given. Furthermore, the sensitivities of the brightness temperature model to SST, SSS, d and AVF of a sea foam layer at L band are discussed. The sensitivities are ranked from most to least in the order： （1） d; （2） AVF; （3） SSS; （4） SST. This result indicates that the measurement errors old and AVF have significant impacts on the retrievals of SSS and SST. With the experimental brightness temperature data, the SSS and AFV are retrieved by cost function.

Free Convection of a Viscous Electrically Conducting Fluid Past a Stretching Surface

Free convection of a viscous electrically conducting liquid past a vertical stretching surface is investigated in the presence of a transverse magnetic field.Natural convection is driven by both thermal and solutal buoyancy.The original partial differential equations governing the problem are turned into a set of ordinary differential equations through a similar variables transformation.This alternate set of equations is solved through a Differential Transform Method(DTM)and the Pade approximation.The response of the considered physical system to the non-dimensional parameters accounting for the relative importance of different effects is assessed considering different situations.

Numerical Treatment of Initial Value Problems of Nonlinear Ordinary Differential Equations by Duan-Rach-Wazwaz Modified Adomian Decomposition Method

We employ the Duan-Rach-Wazwaz modified Adomian decomposition method for solving initial value problems for the systems of nonlinear ordinary differential equations numerically. In order to confirm practicality, robustness and reliability of the method, we compare the results from the modified Adomian decomposition method with those from the MATHEMATICA solutions and also from the fourth-order Runge Kutta method solutions in some cases. Furthermore, we apply Padé approximants technique to improve the solutions of the modified decomposition method whenever the exact solutions exist.

COMPUTATION FORMULAS OF GENERALIZED INVERSE PAD APPROXIMANTS USING FOR SOLUTION OF INTEGRAL EQUATIONS

For the generalized inverse function-valued Pade approximants, its intact computation formulas are given. The explicit determinantal formulas for the denominator scalar polynomials and the numerator function-valued polynomials are first established. A useful existence condition is given by means of determinant form.