Eigenvalue analysis for a crack in a power-law material

A nonlinear eigenvalue problem related to determining the stress and strain fields near the tip of a transverse crack in a power-law material is studied. The eigenvalues are found by a perturbation method based on representations of an eigenvalue, the corresponding eigenfunction, and the material nonlinearity parameter in the form of series expansions in powers of a small parameter equal to the difference between the eigenvalues in the linear and nonlinear problems. The resulting eigenvalues are compared with the accurate numerical solution of the nonlinear eigenvalue problem.     

Analytical solution of time‐fractional Navier–Stokes equation in polar coordinate by homotopy perturbation method

In this letter, we implement a relatively new analytical technique, the homotopy perturbation method (HPM), for solving linear partial differential equations of fractional order arising in fluid mechanics. The fractional derivatives are described in Caputo derivatives. This method can be used as an alternative to obtain analytic and approximate solutions of different types of fractional differential equations applied in engineering mathematics. The corresponding solutions of the integer order equations are found to follow as special cases of those of fractional order equations. Some numerical examples are presented to illustrate the efficiency and reliability of HPM. He's HPM, which does not need small parameter is implemented for solving the differential equations. In this method, a homotopy is introduced to be constructed for the equation. The initial approximations can be freely chosen with possible unknown constants that can be determined by imposing the boundary and initial conditions. It is predicted that HPM can be found widely applicable in engineering. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

Singular Perturbation Analysis of Boundary Value Problems for Differential-Difference Equations. IV. A Nonlinear Example with Layer Behavior

A study is made of boundary value problems for a class of singularly perturbed nonlinear, second-order, differential-difference equations, i.e., where the highest-order derivative is multiplied by a small parameter. Depending on the region of parameter space, solutions of the nonlinear problem may not be unique, can exhibit extreme sensitivity to the values of the parameters, or may not exist. Typically, solutions exhibit layer behavior and/or exponentially large amplitudes. Approximate solutions of these boundary value problems are obtained by using singular perturbation methods and numerical computations and are then compared. Numerical computations of representative solutions illustrate the wide variety of possible behaviors.     

粘性阻尼线性振动系统的复模态特征值问题的一种新的矩阵摄动分析解法

本文提出了对粘性阻尼线性振动系统的复模态二次广义特征值问题进行高效近似求解的一种新的矩阵摄动分析方法,即先将阻尼矩阵分解为比例阻尼部分和非比例阻尼部分之和,并求得系统的比例阻尼实模态特征解;然后以此为初始值,将阻尼矩阵的非比例部分作为对其比例部分的小量修改,利用摄动分析方法简捷地得到系统的复模态特征值问题的近似解.这一新方法适用于振系阻尼分布不十分偏离比例阻尼情况的问题,因此对大阻尼(非过阻尼)振动系统也有效.这是它优于以前提出的基于无阻尼实模态特征解的类似摄动分析方法的重要特点.文中建立了复模态特征值和特征向量的二阶摄动解式,并通过算例证实了其有效性.此外还讨论了利用比例阻尼假定估计阻尼系统固有振动的复特征值的可行性.     

Spectral corrections for Sturm–Liouville problems

The numerical solution of the Sturm–Liouville problem can be achieved using shooting to obtain an eigenvalue approximation as a solution of a suitable nonlinear equation and then computing the corresponding eigenfunction. In this paper we use the shooting method both for eigenvalues and eigenfunctions. In integrating the corresponding initial value problems we resort to the boundary value method. The technique proposed seems to be well suited to supplying a general formula for the global discretization error of the eigenfunctions depending on the discretization errors arising from the numerical integration of the initial value problems. A technique to estimate the eigenvalue errors is also suggested, and seems to be particularly effective for the higher-index eigenvalues. Numerical experiments on some classical Sturm–Liouville problems are presented.     

Perturbation method for solving the nonlinear eigenvalue problem arising from fatigue crack growth problem in a damaged medium

An analytical solution of the nonlinear eigenvalue problem arising from the fatigue crack growth problem in a damaged medium in coupled formulation is obtained. The perturbation technique for solving the nonlinear eigenvalue problem is used. The method allows to find the analytical formula expressing the eigenvalue as the function of parameters of the damage evolution law. It is shown that the eigenvalues of the nonlinear eigenvalue problem are fully determined by the exponents of the damage evolution law. In the paper the third-order (four-term) asymptotic expansions of the angular functions determining the stress and continuity fields in the neighborhood of the crack tip are given. The asymptotic expansions of the angular functions permit to find the closed-form solution for the problem considered.     

Preconditioned iterative methods for a class of nonlinear eigenvalue problems

This paper proposes new iterative methods for the efficient computation of the smallest eigenvalue of symmetric nonlinear matrix eigenvalue problems of large order with a monotone dependence on the spectral parameter. Monotone nonlinear eigenvalue problems for differential equations have important applications in mechanics and physics. The discretization of these eigenvalue problems leads to nonlinear eigenvalue problems with very large sparse ill-conditioned matrices monotonically depending on the spectral parameter. To compute the smallest eigenvalue of large-scale matrix nonlinear eigenvalue problems, we suggest preconditioned iterative methods: preconditioned simple iteration method, preconditioned steepest descent method, and preconditioned conjugate gradient method. These methods use only matrix-vector multiplications, preconditioner-vector multiplications, linear operations with vectors, and inner products of vectors. We investigate the convergence and derive grid-independent error estimates for these methods. Numerical experiments demonstrate the practical effectiveness of the proposed methods for a model problem.     

Analytical solutions of the boundary layer flow of power-law fluid over a power-law stretching surface

This article discusses analytical solutions for a nonlinear problem arising in the boundary layer flow of power-law fluid over a power-law stretching surface. Using perturbation method analytical solution is presented for linear stretching surface. This solution covers large range of shear thinning and shear thickening fluids and matches excellently with the numerical solution. Furthermore, some new exact solutions are found for particular combination of m (power-law stretching index) and n (power-law fluid index). This leads to generalize the case of linear stretching to nonlinear stretching surface. The effects of fluid index n on the boundary layer thickness and the skin friction for nonlinear stretching surface is analyzed and discussed. It is observed that the boundary layer thickness and the skin friction coefficient increase as non-linear parameter n decreases. This study gives a new dimension to obtain analytical solutions asymptotically for highly nonlinear problems which to the best of our knowledge has not been examined so far.     

Linear stability of solitary waves of the Green‐Naghdi equations

We investigate the eigenvalue problem obtained from linearizing the Green‐Naghdi equations about solitary wave solutions. Unlike weakly nonlinear water wave models, the physical system considered here has nonlinearity in its highest derivative term. This results in more detailed asymptotic analysis of the eigenvalue problem in the presence of a large parameter. Combining the technique of singular perturbation with the Evans function, we show that for solitary waves of small amplitude, the problem has no eigenvalues of positive real part and the Evans function is nonvanishing everywhere except the origin. This fact then leads to the linear stability of these solitary waves. © 2001 John Wiley & Sons, Inc.     

几类微分-代数方程的神经网络求解法

在非线性科学中，寻求微分方程的近似解析解一直是重要的研究课题和研究热点.利用人工神经网络原理，结合最优化方法，研究了几类微分-代数方程的近似解析解，包括指标1，2，3型Hessenberg方程及指标3型Euler-Lagrange方程，得到了方程近似解析解的表达式.通过与精确解或Runge-Kutta(龙格-库塔)数值计算结果对比，表明神经网络方法的结果有很高的精度.     

The numerical treatment of perturbed bifurcation problems in ordinary differential equations

A new numerical method is presented to analyze perturbations of bifurcations of the solutions of nonlinear boundary value problems. The perturbations may result from imperfections or other inhomogeneities in the corresponding scientific problem. The nonisolated solutions are calculated in dependence of the perturbation parameter. Therefore, it is possible to determine the singular solution as well as a solution branch through this nonisolated solution simultaneously. Standard procedures of numerical analysis are applicable.     

在材料非线性问题中的摄动有限元法

摄动法是解决非线性连续介质力学问题的一种有效方法.这种方法是建立在该问题的线性解析解的基础上的,因此,若得不到一个简单的解析解,应用这种方法去解决一些复杂的非线性问题将遇到困难.有限元法对解非线性问题也是一种十分有用的工具,然而一般来说,它需要相当长的计算时间. 本文介绍摄动有限元法.这种方法吸取上述两种方法的优点,能够解决更复杂的非线性问题,而且也能大量节省计算机的计算时间. 本文讨论了比例加载下的弹塑性力学问题,并提出一个带孔拉板的数值解.     

Comparison of the Adomian decomposition method with homotopy perturbation method for the solutions of seventh order boundary value problems

The aim of this paper is to compare the Adomian decomposition method and the homotopy perturbation method for solving the linear and nonlinear seventh order boundary value problems. The approximate solutions of the problems obtained with a small amount of computation in both methods. Two numerical examples have been considered to illustrate the accuracy and implementation of the methods.     

源于膨胀波边界层理论中的一类奇异边值问题

对源于膨胀波边界层理论中的一类奇异边值问题进行了研究,利用单调逼近方法得到了满足物理意义上正解的存在性和唯一性的充分条件,同时给出了用速度比例参数表示的壁摩擦的估计公式, 利用打靶法技巧得到了问题的数值解.数值结果证明了估计公式的可靠性和有效性.     

Application of He's variational iteration method to nonlinear Jaulent–Miodek equations and comparing it with ADM

Instead of finding a small parameter for solving nonlinear problems through perturbation method, a new analytical method called He's variational iteration method (VIM) is introduced to be applied to solve nonlinear Jaulent–Miodek, coupled KdV and coupled MKdV equations in this article. In this method, general Lagrange multipliers are introduced to construct correction functionals for the problems. The multipliers can be identified optimally via the variational theory. The results are compared with exact solutions.     

Multiscaling analysis of a slowly varying anaerobic digestion model

We construct the slowly varying limiting state solutions to a nonlinear dynamical system for anaerobic digestion with Monod-based kinetics involving slowly varying model parameters arising from slow environmental variation. The advantage of these approximate solutions over numerical solutions is their applicability to a wide range of parameter values. We use these limiting state solutions to develop analytic approximations to the full nonlinear system by applying a multiscaling technique. The approximate solutions are shown to compare favorably with numerical solutions.     

椭圆-抛物偏微分方程奇异摄动问题的差分解法

本文讨论了含有小参数在高阶导数项的椭圆型方程奇异摄动问题的差分解法.当ε=0时椭圆型方程退化为抛物型方程.作者根据此问题解的边界层性质,构造了特殊的差分格式:研究了它的收敛性和解的渐近性态.最后给出一个数值例题.     

Finding Eigenvalues of Holomorphic Fredholm Operator Pencils Using Boundary Value Problems and Contour Integrals

Investigating the stability of nonlinear waves often leads to linear or nonlinear eigenvalue problems for differential operators on unbounded domains. In this paper we propose to detect and approximate the point spectra of such operators (and the associated eigenfunctions) via contour integrals of solutions to resolvent equations. The approach is based on Keldysh’ theorem and extends a recent method for matrices depending analytically on the eigenvalue parameter. We show that errors are well-controlled under very general assumptions when the resolvent equations are solved via boundary value problems on finite domains. Two applications are presented: an analytical study of Schrödinger operators on the real line as well as on bounded intervals and a numerical study of the FitzHugh–Nagumo system. We also relate the contour method to the well-known Evans function and show that our approach provides an alternative to evaluating and computing its zeros.     

A Jacobi–Davidson method for two‐real‐parameter nonlinear eigenvalue problems arising from delay‐differential equations

The critical delays of a delay‐differential equation can be computed by solving a nonlinear two‐parameter eigenvalue problem. The solution of this two‐parameter problem can be translated to solving a quadratic eigenvalue problem of squared dimension. We present a structure preserving QR‐type method for solving such quadratic eigenvalue problem that only computes real‐valued critical delays; that is, complex critical delays, which have no physical meaning, are discarded. For large‐scale problems, we propose new correction equations for a Newton‐type or Jacobi–Davidson style method, which also forces real‐valued critical delays. We present three different equations: one real‐valued equation using a direct linear system solver, one complex valued equation using a direct linear system solver, and one Jacobi–Davidson style correction equation that is suitable for an iterative linear system solver. We show numerical examples for large‐scale problems arising from PDEs. Copyright © 2012 John Wiley & Sons, Ltd.     

Boundary layer analysis for the eigenvalues of hanging rod

The motion of a naturally straight inextensible flexible elastic hanging rod is formulated and then linearized about the straight solution. To solve this equation by separation of variables, an eigenvalue problem is derived. When the stiffness of the rod is small, the eigenvalue equation is a singular perturbation problem. This paper is devoted to solving this eigenvalue problem by boundary layer analysis when the stiffness is suitably small, especially on the analytic approximate solutions of the first several eigenvalues and eigenfunctions. The first three eigenvalues are also compared with the numerical results computed by a finite difference method. The excellent agreement shows the efficiency of the boundary layer analysis.