Adomian修正分解法的后处理方法比较

Adomian修正分解法在求解非线性微分方程中得到广泛应用。Adomian修正分解法的主要特点在于计算简单快速,并且不需要进行线性化或离散化。但是Adomian修正分解法的计算精度取决于其收敛域。为了扩大Adomian修正分解法的收敛域,需要对所得解进行后处理,目前常见的后处理方法包括Padé近似、LaplacePadé近似和多步迭代方法。本文首先简要回顾了Adomian修正分解法,然后讨论了这三种后处理方法,最后通过Duffing振子为例对这些后处理方法的优缺点进行讨论和分析。数值计算结果表明,多步迭代方法能够加速Adomian修正分解法解的收敛,并扩大其收敛域。     

Modified asymptotic Adomian decomposition method for solving Boussinesq equation of groundwater flow

The Adomian decomposition method （ADM） is an approximate analytic method for solving nonlinear equations. Generally, an approximate solution can be ob- tained by using only a few terms. However, in applications, we need to use it flexibly according to the real problem. In this paper, based on the ADM, we give a modified asymptotic Adomian decomposition method and use it to solve the nonlinear Boussinesq equation describing groundwater flows. The example shows effectiveness of the modified asymptotic Adomian decomposition method.     

通过Adomian分解法求解二维Helmholtz方程

提出基于Adomian分解法求解二维Helmholtz方程。通过Adomian分解法可以把Helmholtz微分方程和边界条件分别转换成递归代数公式和适用符号计算的简单代数公式。利用边界条件可以很容易得到方程的解析解表达式。Adomian分解法的主要特点在于计算简单快速,并且不需要进行线性化或离散化。最后给出数值实例以验证Adomian分解法求解二维Helmholtz方程的有效性。通过数值计算可以发现,基于Adomian分解法的计算结果非常接近精确解,并且该方法具有良好的收敛性。这表明Adomian分解法能够快速有效求解Helmholtz方程。     

Analytical solution of fractionally damped beam by Adomian decomposition method

The analytical solution of a viscoelastic continuous beam whose damping characteristics are described in terms of a fractional derivative of arbitrary order was derived by means of the Adomian decomposition method.The solution contains arbitrary initial conditions and zero input.For specific analysis,the initial conditions were assumed homogeneous,and the input force was treated as a special process with a particular beam. Two simple cases,step and impulse function responses,were considered respectively. Subsequently,some figures were plotted to show the displacement of the beam under different sets of parameters including different orders of the fractional derivatives.     

Analytical investigation of Jeffery-Hamel flow with high magnetic field and nanoparticle by Adomian decomposition method

In this study, the effects of magnetic field and nanoparticle on the Jeffery-Hamel flow are studied using a powerful analytical method called the Adomian decomposition method (ADM). The traditional Navier-Stokes equation of fluid mechanics and Maxwell’s electromagnetism governing equations are reduced to nonlinear ordinary differential equations to model the problem. The obtained results are well agreed with that of the Runge-Kutta method. The present plots confirm that the method has high accuracy for different α, Ha, and Re numbers. The flow field inside the divergent channel is studied for various values of Hartmann number and angle of channel. The effect of nanoparticle volume fraction in the absence of magnetic field is investigated.     

Blow-up solutions of the generalized Boussinesq equation obtained by variational iteration method

This paper deals with obtaining explicit solutions of a generalized non-linear Boussinesq equation using He’s variational iteration method. Both finite and blow-up solutions can be obtained.     

Generalized variational principles of the viscoelastic body with voids and their applications

From the Boltzmann‘ s constitutive law of viscoelastic materials and the linear theory of elastic materials with voids, a constitutive model of generalized force fields for viscoelastic solids with voids was given. By using the variational integral method, the convolution-type functional was given and the corresponding generalized variational principles and potential energy principle of viscoelastic solids with voids were presented. It can be shown that the variational principles correspond to the differential equations and theinitial and boundary conditions of viscoelastic body with voids. As an application, a generalized variational principle of viscoelastic Timoshenko beams with damage was obtained which corresponds to the differential equations of generalized motion and the initial and boundary conditions of beams. The variational principles provide a way for solving problems of viscoelastic solids with voids.     

降维变分方法

依据能量原理建立降维变分法,该方法将物体离散为条,在一定约束下,一条独立求解能量极值,可有效简化数值计算,本文研究分条离散方法,极值约束,边界条件处理并给出证明.更重要的,该方法揭示了新的物体形变运动规律.在此也做了分析.     

AN ITERATIVE METHOD FOR THE DISCRETE PROBLEMS OF A CLASS OF ELLIPTICAL VARIATIONAL INEQUALITIES

ＡＮＩＴＥＲＡＴＩＶＥＭＥＴＨＯＤＦＯＲＴＨＥＤＩＳＣＲＥＴＥＰＲＯＢＬＥＭＳＯＦＡＣＬＡＳＳＯＦＥＬＬＩＰＴＩＣＡＬＶＡＲＩＡＴＩＯＮＡＬＩＮＥＱＵＡＬＩＴＩＥＳＺｈｅｎｇＴｉｅ－ｓｈｅｎｓ（郑铁生）ＬｉＬｉ（李立）ＸｕＱｉｎｇ－ｙｕ（许庆余）（Ｄ...     

The variational principles and application of nonlinear numerical manifold method

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｆｉｎｉｔｅｃｏｖｅｒｔｅｃｈｎｉｑｕｅｓ,ｓｉｍｉｌａｒｔｏｔｈｅｃｏｎｃｅｐｔｉｏｎｏｆｆｉｎｉｔｅｃｏｖｅｒｕｓｅｄｉｎｍａｎｉｆｏｌｄａｎａｌｙｓｉｓｏｆｍｏｄｅｒｎｍａｔｈｅｍａｔｉｃｓ,ａｒｅｉｎｔｒｏｄｕｃｅｄｉｎｔｏＮｕｍｅｒｉｃａｌＭａｎｉｆｏｌｄｍｅｔｈｏｄ (ＮＭＭ)ａｎｄｔｈｅｆｉｎｉｔｅｃｏｖｅｒｓｃｏｎｓｉｓｔｏｆｍａｔｈｅｍａｔｉｃａｌｃｏｖｅｒｓａｎｄｐｈｙｓｉｃａｌｃｏｖｅｒｓｗｈｉｃｈｃａｎｂｅｓｅｐａｒａｔｅｄ .ＴｈｅＮＭＭｉｓｃｏ…     

Application of the modified iteration method to nonlinear postbuckling analysis of thin circular plates

In this paper,the modified iteration method is further generalized to the study ofaxisymmetrical postbuckling of thin circular plates and hereby a new approximate analyticsolution of the problem is obtained.Further utilizations of this method to postbucklinganalyses of plates of more complicated structure are expected.     

Numerical solution of stiff systems from nonlinear dynamics using single-term Haar wavelet series

The paper presents single-term Haar wavelet series (STHWS) approach to the solution of nonlinear stiff differential equations arising in nonlinear dynamics. The properties of STHWS are given. The method of implementation is discussed. Numerical solutions of some model equations are investigated for their stiffness and stability and solutions are obtained to demonstrate the suitability and applicability of the method. The results in the form of block-pulse and discrete solutions are given for typical nonlinear stiff systems. As compared with the TR BDF2 method of Shampine and Gill’s method, the STHWS turns out to be more effective in its ability to solve systems ranging from mildly to highly stiff equations and is free from stability constraints.     

Study of size effects in the Dugdale model through the case of a crack in a semi-infinite plane under anti-plane shear loading

The main objective of this work is to prove that, with the Dugdale model, the small size defects, comparatively to the material characteristic length, are practically without influence on the limit load of structures. For that, we treat the case of a crack in a semi-infinite plane under anti-plane shear loading. Using integral transforms, the elasticity equations are converted analytically into a singular integral equation. The singular integral equation is solved numerically using Chebychev polynomials. Special care is needed to take into account the presence of jump discontinuities in the loading distribution along the crack lips.      

Frictional passage of fluid through inhomogeneous porous system: a variational principle

A Fermat-like principle of minimum time is formulated for nonlinear steady paths of fluid flow in inhomogeneous isotropic porous media where fluid streamlines are curved by a location dependent hydraulic conductivity. The principle describes an optimal nature of nonlinear paths in steady Darcy’s flows of fluids. An expression for the total path resistance leads to a basic analytical formula for an optimal shape of a steady trajectory. In the physical space an optimal curved path ensures the maximum flux or shortest transition time of the fluid through the porous medium. A sort of “law of bending” holds for the frictional fluid flux in Lagrange coordinates. This law shows that—by minimizing the total resistance—a ray spanned between two given points takes the shape assuring that a relatively large part of it resides in the region of lower flow resistance (a ‘rarer’ region of the medium).     

Asymptotic approximation method for elliptic variational inequality of first kind简

An asymptotic approximation method is proposed to solve a particular elliptic variational inequality of first kind associated with unilateral obstacle problems. In this method, the free boundary is first captured, and then the method of the fundamental solution (MFS) is used to find the solution of the Dirichlet problem for Laplace’s equation in the non-coincidence set. Numerical examples are given to show the efficiency of the method.     

One-dimensional unsteady inertial flow in phreatic aquifers induced by a sudden change of the boundary head

We are examining the classical problem of unsteady flow in a phreatic semi-infinite aquifer, induced by sudden rise or drawdown of the boundary head, by taking into account the influence of the inertial effects. We demonstrate that for short times the inertial effects are dominant and the equation system describing the flow behavior can be reduced to a single ordinary differential equation. This equation is solved both numerically by the Runge-Kutta method and analytically by the Adomian’s decomposition approach and an adequate polynomial-exponential approximation as well. The influence of the viscous term, occurring for longer times, is also taken into account by solving the full Forchheimer equation by a finite difference approach. It is also demonstrated that as for the Darcian flow, for the case of small fluctuations of the water table, the computation procedure can be simplified by using a linearized form of the mass balance equation. Compact analytical expressions for the computation of the water stored or extracted from an aquifer, including viscous corrections are also developed.     

Implicit solution of harmonic balance equation system using the LU-SGS method and one-step Jacobi/Gauss-Seidel iteration

ABSTRACT

This paper presents efficient alternative numerical methods for an implicit solution of the harmonic balance equation system for analysing temporal periodic unsteady flows. The proposed method employs approximate factorisation to decouple the common residual term and the time spectral source term of a harmonic balance equation system when it is discretised implicitly. With this approximate factorisation, the complexity of implicit solution of the discrete system is greatly reduced. The common residual term can be dealt with using a lower-upper symmetric-Gauss-Seidel (LU-SGS) method and the time spectral source term is integrated using a Jacobi iteration (JI) or one step Gauss-Seidel (GS) iteration, leading to the LU-SGS/JI method or LU-SGS/GS method. The NASA stage 35 compressor and the 1.5 stage Aachen turbine were used to demonstrate the effectiveness of the proposed methods in stabilising solution and its advantages in comparison with the existing lower-upper symmetric-Gauss-Seidel/block Jacobi (LU-SGS/BJ) method. The LU-SGS/GS method and the LU-SGS/JI method are more robust than the LU-SGS/BJ method in stabilising solution. The LU-SGS/GS method also has faster and tighter convergence and lower memory consumption in comparison with the LU-SGS/BJ method.     

A proper orthogonal decomposition variational multiscale meshless interpolating element-free Galerkin method for incompressible magnetohydrodynamics flow

In the recent decade, the meshless methods have been handled for solving most of PDEs due to easiness of the meshless methods. One of the popular meshless methods is the element-free Galerkin (EFG) method that was first proposed for solving some problems in the solid mechanics. The test and trial functions of the EFG are based on the special basis. Recently, some modifications have been developed to improve the EFG method. One of these improvements is the variational multiscale EFG procedure. In the current article, the shape functions of interpolation moving least squares approximation have been applied to the variational multiscale EFG technique for solving the incompressible magnetohydrodynamics flow. In order to reduce the elapsed CPU time of simulation, we employ a reduced-order model based on the proper orthogonal decomposition technique. The current combination can be referred to as the reduced-order variational multiscale EFG technique. To illustrate the reduction in CPU time used as well as the efficiency of the proposed method, we applied it for the two-dimensional cases.     

PARAMETRIC VARIATIONAL PRINCIPLE BASED ELASTIC-PLASTIC ANALYSIS OF HETEROGENEOUS MATERIALS WITH VORONOI FINITE ELEMENT METHOD

The Voronoi cell finite element method (VCFEM) is adopted to overcome the limitations of the classic displacement based finite element method in the numerical simulation of heterogeneous materials. The parametric variational principle and quadratic programming method are developed for elastic-plastic Voronoi finite element analysis of two-dimensional problems. Finite element formulations are derived and a standard quadratic programming model is deduced from the elastic-plastic equations. Influence of microscopic heterogeneities on the overall mechanical response of heterogeneous materials is studied in detail. The overall properties of heterogeneous materials depend mostly on the size, shape and distribution of the material phases of the microstructure. Numerical examples are presented to demonstrate the validity and effectiveness of the method developed.