Extended homotopy perturbation method and the axisymmetric flow past a porous stretching sheet

The extended homotopy perturbation method, which is an extension of the celebrated homotopy perturbation method (HPM), is applied to obtain a solution to the problem of the steady, laminar, axisymmetric flow of a viscous, incompressible fluid past a porous stretching sheet. The solution so obtained is totally analytical and is expressible in terms of the cross‐flow velocity of the fluid past the stretching sheet. Its hallmark is that it does not depend upon computation of any auxiliary parameter for enlarging the convergence region of the solution. Rather, it calculates the solution automatically adjusting the scaling factor of the independent similarity variable normal to the sheet. The results obtained by the extended HPM are in excellent agreement with the exact numerical solution. Also, an asymptotic solution valid for large suction parameter is developed, which matches well with the exact solution even for moderate values of the suction parameter. Copyright © 2011 John Wiley & Sons, Ltd.     

Laminar flow through a channel with uniformly porous walls of different permeability

Summary The steady laminar flow of a viscous incompressible fluid through a two-dimensional channel, having fluid sucked or injected with different velocities through its uniformly porous parallel walls is considered. A solution for small suction Reynolds number has been given by the authors in a previous paper. The purpose of this paper is to present a solution valid for large Reynolds numbers for the cases of (i) suction at both walls, and (ii) suction at one wall and injection at the other. A technique of matching outer and inner expansions is used to obtain an asymptotic solution for both of these cases. Further a perturbation solution for the case of suction at one wall and injection at the other is obtained by choosing the difference between two wall velocities as the perturbation parameter. Both asymptotic and perturbation solutions are confirmed by exact numerical solutions. As expected, the resulting solutions show the presence of the usual suction boundary layers in both types of flow considered in this paper.     

Stability of nonlinear time-varying perturbed differential equations

The goal of this paper is twofold. The first part presents a converse Lyapunov theorem for the notion of uniform practical exponential stability of nonlinear differential equations in presence of small perturbation. This class of nonlinear differential equations can be viewed as parametric differential equations. The second part provides the classical perturbation method of seeking an approximate solution as a finite Taylor expansion of the exact solution. The practical asymptotic validity on the approximate is established on infinite-time interval. Finally, we give a numerical example to prove the validity of our methods.     

均布荷载作用下功能梯度悬臂梁弯曲问题的解析解

采用弹性力学半逆解法,假设所有材料常数沿梁厚度方向按同一函数规律变化,求得了功能梯度悬臂梁在均布载荷作用下的解析解.该解退化到各向同性均匀弹性情况时与已有的理论解相一致.对弹性模量按指数函数梯度变化的算例进行了分析.所得到的解对任意梯度函数均成立,可作为数值解以及简化理论的检验依据.     

High order symplectic conservative perturbation method for time-varying Hamiltonian system

This paper presents a high order symplectic conservative perturbation method for linear time-varying Hamiltonian system.Firstly,the dynamic equation of Hamiltonian system is gradually changed into a high order perturbation equation,which is solved approximately by resolving the Hamiltonian coefficient matrix into a "major component" and a "high order small quantity" and using perturbation transformation technique,then the solution to the original equation of Hamiltonian system is determined through a series of inverse transform.Because the transfer matrix determined by the method in this paper is the product of a series of exponential matrixes,the transfer matrix is a symplectic matrix;furthermore,the exponential matrices can be calculated accurately by the precise time integration method,so the method presented in this paper has fine accuracy,efficiency and stability.The examples show that the proposed method can also give good results even though a large time step is selected,and with the increase of the perturbation order,the perturbation solutions tend to exact solutions rapidly.     

Radiocarbon dating in trondheim

Summary Naturally occurring radiocarbon is registered in a proportional counter filled with pure CO₂. THe CO₂ gas is purified in almost the same manner as has been developed by De Vries and Barendsen. Checking of the gas purity is described. A grid proportional counter with an effective volume of about 2.5 liter has been used with CO₂-pressures from 1 to 4 atmospheres. At a pressure of 4 atmospheres the net count from contemporary wood is 58 counts/min above a background of 15 eounts/min. Dating may at the present stage be extended to 35 000 years by 24 hours counting with 4 atmospheres CO₂. A method for checking the discriminating level of the counting apparatus is described.     

多自由度参激系统稳定性的数值分析

提出多自由度周期参激系统稳定性的数值直接法。通过将扰动方程表示成状态方程形式,再根据Flo-quet理论将扰动解表示成指数特征分量与周期分量之积,并将其周期分量与系统周期系数展成Fourier级数,导出一系列代数方程,建立矩阵特征值问题,从而由数值求解特征值可直接确定参激系统的稳定性。该方法可用于一般周期参激阻尼系统,特征值矩阵不含逆子阵。应用于斜拉索在支座周期运动激励下的参激振动不稳定性分析,数值结果表明该方法的有效性。     

Application of a modified Lindstedt–Poincaré method in coupled TDOF systems with quadratic nonlinearity and a constant external excitation

An analytical approach is developed for the nonlinear oscillation of a conservative, two-degree-of-freedom (TDOF) mass-spring system with serial combined linear–nonlinear stiffness excited by a constant external force. The main idea of the proposed approach lies in two categories, the first one is the transformation of two nonlinear differential equations of a two-mass system using suitable intermediate variables into a single nonlinear differential equation. Another is the treatment a quadratic nonlinear oscillator (QNO) by the modified Lindstedt–Poincaré (L-P) method presented recently by the authors. The first-order and second-order analytical approximations for the modified L-P method are established for the QNOs with satisfactory results. After solving the nonlinear differential equation, the displacements of two-mass system can be obtained directly from the governing linear second-order differential equation. Unlike the common perturbation method, the modified L-P method is valid for weak as well as strong nonlinear oscillation systems. On the other hand, the new approach yields simple approximate analytical expressions valid for small as well as large amplitudes of oscillation. In short, this new approach yields extended scope of applicability, simplicity, flexibility in application, and avoidance of complicated numerical integration as compared to the previous approaches such as the perturbation and classical harmonic balance methods. Two examples of nonlinear TDOF mass-spring systems excited by a constant external force are selected and the approximate solutions are verified with the exact solutions derived from the Jacobi elliptic function and also the numerical fourth-order Runge–Kutta solutions.     

Numerical modeling of perturbation propagation in a supersonic boundary layer

The growth of two-dimensional disturbances generated in a supersonic (M = 6) boundary layer on a flat plate by a periodic perturbation of the injection/suction type is investigated on the basis of a numerical solution of the Navier-Stokes equations. For small initial perturbation amplitudes, the second-mode growth rate obtained from the numerical modeling coincides with the growth rate calculated using linear theory with account for the non-parallelism of the main flow. Calculations performed for large initial perturbation amplitudes reveal the nonlinear dynamics of the perturbation growth downstream, with rapid growth of the higher multiple harmonics.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, 2004, pp. 33–44. Original Russian Text Copyright © 2004 by Egorov, Sudakov, Fedorov.     

Nonlinear vibration of a two-mass system with nonlinear stiffnesses

An analytical approach is developed for nonlinear free vibration of a conservative, two-degree-of-freedom mass–spring system having linear and nonlinear stiffnesses. The main contribution of the proposed approach is twofold. First, it introduces the transformation of two nonlinear differential equations of a two-mass system using suitable intermediate variables into a single nonlinear differential equation and, more significantly, the treatment a nonlinear differential system by linearization coupled with Newton’s method and harmonic balance method. New and accurate higher-order analytical approximate solutions for the nonlinear system are established. After solving the nonlinear differential equation, the displacement of two-mass system can be obtained directly from the governing linear second-order differential equation. Unlike the common perturbation method, this higher-order Newton–harmonic balance (NHB) method is valid for weak as well as strong nonlinear oscillation systems. On the other hand, the new approach yields simple approximate analytical expressions valid for small as well as large amplitudes of oscillation unlike the classical harmonic balance method which results in complicated algebraic equations requiring further numerical analysis. In short, this new approach yields extended scope of applicability, simplicity, flexibility in application, and avoidance of complicated numerical integration as compared to the previous approaches such as the perturbation and the classical harmonic balance methods. Two examples of nonlinear two-degree-of-freedom mass–spring system are analyzed and verified with published result, exact solutions and numerical integration data.     

Nonlinear oscillations with parametric excitation solved by homotopy analysis method

An analytical technique, namely the homotopy analysis method （HAM）, is used to solve problems of nonlinear oscillations with parametric excitation. Unlike perturbation methods, HAM is not dependent on any small physical parameters at all, and thus valid for both weakly and strongly nonlinear problems. In addition, HAM is different from all other analytic techniques in providing a simple way to adjust and control convergence region of the series solution by means of an auxiliary parameter h. In the present paper, a periodic analytic approximations for nonlinear oscillations with parametric excitation are obtained by using HAM, and the results are validated by numerical simulations.     

Perturbation methods for nonlinear autonomous discrete-time dynamical systems

Two perturbation methods for nonlinear autonomous discrete-time dynamical systems are presented. They generalize the classical Lindstedt-Poincaré and multiple scale perturbation methods that are valid for continuous-time systems. The Lindstedt-Poincaré method allows determination of the periodic or almost-periodic orbits of the nonlinear system (limit cycles), while the multiple scale method also permits analysis of the transient state and the stability of the limit cycles. An application to the discrete Van der Pol equation is also presented, for which the asymptotic solution is shown to be in excellent agreement with the exact (numerical) solution. It is demonstrated that, when the sampling step tends to zero the asymptotic transient and steady-state discrete-time solutions correctly tend to the asymptotic continuous-time solutions.     

Influence of the velocity gradient on the stagnation point heating in hypersonic flow

In a number of experimental and numerical publications a deviation has been found between the measured or computed stagnation point heat flux and that given by the theory of Fay and Riddell. Since the formula of Fay and Riddell is used in many applications to yield a reference heat flux for experiments performed in wind tunnels, for flight testing and numerical simulations, it is important that this reference heat flux is as accurate as possible. There are some shortcomings in experiments and numerical simulations which are responsible in some part for the deviations observed. But, as will be shown in the present paper, there is also a shortcoming on the theoretical side which plays a major role in the deviation between the theoretical and experimental/numerical stagnation point heat fluxes. This is caused by the method used so far to determine the tangential velocity gradient at the stagnation point. This value is important for the stagnation point heat flux, which so far has been determined by a simple Newtonian flow model. In the present paper a new expression for the tangential velocity gradient is derived, which is based on a more realistic flow model. An integral method is used to solve the conservation equations and, for the stagnation point, yields an explicit solution of the tangential velocity gradient. The solution achieved is also valid for high temperature flows with real gas effects. A comparison of numerical and experimental results shows good agreement with the stagnation point heat flux according to the theory of Fay and Riddell, if the tangential velocity gradient is determined by the new theory presented in this paper.This article was processed by the author using theLAT_EX style filepljour2 from Springer-Verlag.     

Exponential dichotomies of nonlinear discrete systems and its application to numerical analysis and computation

In this pepar we consider the upwind difference scheme of a kind of boundary value problems for nonlinear, second order, ordinary differential equations. Singular perturbation method is applied to construct the asymptotic approximation of the solution to the upwind difference equation. Using the theory of exponential dichotomies we show that the solution of an order-reduced equation is a good approximation of the solution to the upwind difference equation except near boundaries. We construct correctors which yield asymptotic approximations by adding them to the solution of the order-reduced equation. Finally, some numerical examples are illustrated.     

Finite element solution of the Orr–Sommerfeld equation using high precision Hermite elements: plane Poiseuille flow

This paper presents a comprehensive review of the numerical techniques used during the past half century and their accuracy in hydrodynamic stability analysis of plane parallel flows. The paper also describes a finite element solution of the Orr–Sommerfeld equation using high precision Hermite elements. A stability analysis technique is performed by imposing an infinitesimal perturbation to the laminar base flow to determine the thresholds of neutral instabilities or the growth rate of the perturbation for any Reynolds and wave numbers. Validation of the present numerical technique is performed for plane Poiseuille flow. The numerical results, obtained with uniform and nonuniform meshes, show excellent agreement with the most accurate results available in the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

非牛顿流体两相渗流问题的摄动解

本文用奇异摄动法结合正则摄动法求解了考虑毛管力因素时多孔介质中弱非牛顿流体的两相驱替问题,得到了分流函数和湿相饱和度的渐近解析解。所得结果同数值解和经典的牛顿流体两相渗流结果进行了比较,并着重讨论了非牛顿因素的影响。     

多尺度法二义性的一种解释

多尺度法是为解决含小参数系统发展起来的应用最广泛的摄动法之一. 在求解高阶近 似方程时，多尺度法一般只求特解. 用多尺度法求解van der Pol 方程的三阶解时 将出现矛盾. 以van der Pol方程为例，证明了忽略一阶修正量中的一阶谐波 项使得混合偏导数不能交换顺序，从而导致了多尺度法的二义性和另一个数学矛盾. 在求解一阶修正量时采用含有一阶谐波项的全解，消除了二义性和该矛盾. 该 方法所求得的近似解与数值解进行了比较，结果非常吻合，验证了其合理性.     

Singular perturbation for the weakly nonlinear reaction diffusion equation with boundary perturbation

In this paper, a class of nonlinear singularly perturbed initial boundary value problems for reaction diffusion equations with boundary perturbation are considered under suitable conditions. Firstly, by dint of the regular perturbation method, the outer solution of the original problem is obtained. Secondly, by using the stretched variable and the expansion theory of power series the initial layer of the solution is constructed. And then, by using the theory of differential inequalities, the asymptotic behavior of the solution for the initial boundary value problems is studied. Finally, using some relational inequalities the existence and uniqueness of solution for the original problem and the uniformly valid asymptotic estimation are discussed.     

A bending quasi-front generated by an instantaneous impulse loading at the edge of a semi-infinite pre-stressed incompressible elastic plate

A refined membrane-like theory is used to describe bending of a semi-infinite pre-stressed incompressible elastic plate subjected to an instantaneous impulse loading at the edge. A far-field solution for the quasi-front is obtained by using the method of matched asymptotic expansions. A leading-order hyperbolic membrane equation is used for an outer problem, whereas a refined (singularly perturbed) membrane equation of an inner problem describes a boundary layer, which smoothes a discontinuity predicted by the outer problem at the wave front. The inner problem is then reduced to one-dimensional by an appropriate choice of inner coordinates, motivated by the wave front geometry. Using the inherent symmetry of the outer problem, a solution for the quasi-front is derived that is valid in a vicinity of the tip of the wave front. Pre-stress is shown to affect geometry and type of the generated quasi-front; in addition to a classical receding quasi-front the pre-stressed plate can support propagation of an advancing quasi-front. Possible responses may even feature both types of quasi-front at the same time, which is illustrated by numerical examples. The case of a so-called narrow quasi-front, associated with a possible degeneration of contribution of singular perturbation terms to the governing equation, is studied qualitatively.     

FD-method for a nonlinear eigenvalue problem with discontinuous eigenfunctions

An algorithm for the solution of a nonlinear eigenvalue problem with discontinuous eigenfunctions is developed. The numerical technique is based on a perturbation of the coefficients of a differential equation combined with the Adomian decomposition method for the nonlinear term of the equation. The proposed approach provides an exponential convergence rate dependent on the index of the trial eigenvalue and on the transmission coefficient. Numerical examples support the theory. Published in Neliniini Kolyvannya, Vol. 10, No. 1, pp. 126–143, January–March, 2007.