PERTURBATION METHOD FOR REANALYSIS OF THE MATRIX SINGULAR VALUE DECOMPOSITION

The perturbation method for the reanalysis of the singular value decomposition(SVD)of general real matrices is presented in this paper.This is a simple but efficientreanalysis technique for the SVD,which is of great worth to enhance computationalefficiency of the iterative analysis problems that require matrix singular valuedecomposition repeatedly.The asymptotic estimate formulas for the singular values and thecorresponding left and right singular vectors up to second-order perturbation componentsare derived.At the end of the paper the way to extend the perturbation method to the case ofgeneral complex matrices is advanced.     

Direct perturbation method for reanalysis of matrix singular value decomposition

I.IntroductionMatrixsingularvaluedecomposition(SVD)isoneofthemostimportantandfundamentalcomputationalanalysistoolsformodernnumericallinearalgebra.ThematrixSVDhasverygoodnumericalstability,andsoitisthemostreliableandbeautifulnumericalanalysismethodinmanyth…     

Axisymmetric viscoplastic deformation by the boundary element method

This paper presents a boundary element formulation and numerical implementation of the problem of small axisymmetric deformation of viscoplastic bodies. While the extension from planar to axisymmetric problems can be carried out fairly simply for the finite element method (FEM), this is far from true for the boundary element method (BEM). The primary reason for this fact is that the axisymmetric kernels in the integral equations of the BEM contain elliptic functions which cannot be integrated analytically even over boundary elements and internal cells of simple shape. Thus, special methods have to be developed for the efficient and accurate numerical integration of these singular and sensitive kernels over discrete elements. The accurate determination of stress rates by differentiation of the displacement rates presents another formidable challenge.A successful numerical implementation of the boundary element method with elementwise (called the Mixed approach) or pointwise (called the pure BEM or BEM approach) determination of stress rates has been carried out. A computer program has been developed for the solution of general axisymmetric viscoplasticity problems. Comparisons of numerical results from the BEM and FEM, for several illustrative problems, are presented and discussed in the paper. It is possible to get direct solutions for the simpler class of problems for cylinders of uniform cross-section, and these solutions are also compared with the BEM and FEM results for such cases.     

A PIECEWISE–LINEARIZED METHOD FOR ORDINARY DIFFERENTIAL EQUATIONS: TWO–POINT BOUNDARY–VALUE PROBLEMS

Piecewise-linearized methods for the solution of two-point boundary value problems in ordinary differential equations are presented. These problems are approximated by piecewise linear ones which have analytical solutions and reduced to finding the slope of the solution at the left boundary so that the boundary conditions at the right end of the interval are satisfied. This results in a rather complex system of non-linear algebraic equations which may be reduced to a single non-linear equation whose unknown is the slope of the solution at the left boundary of the interval and whose solution may be obtained by means of the Newton–Raphson method. This is equivalent to solving the boundary value problem as an initial value one using the piecewise-linearized technique and a shooting method. It is shown that for problems characterized by a linear operator a technique based on the superposition principle and the piecewise-linearized method may be employed. For these problems the accuracy of piecewise-linearized methods is of second order. It is also shown that for linear problems the accuracy of the piecewise-linearized method is superior to that of fourth-order-accurate techniques. For the linear singular perturbation problems considered in this paper the accuracy of global piecewise linearizat ion is higher than that of finite difference and finite element methods. For non-linear problems the accuracy of piecewise-linearized methods is in most cases lower than that of fourth-order methods but comparable with that of second-order techniques owing to the linearization of the non-linear terms.     

Stable second-order accurate iterative solutions for second-order elliptic problems

The paper describes a numerical scheme for solving a convection–diffusion elliptic system with very small diffusion coefficients. This iterative numerical procedure is unconditionally stable and converges very rapidly. Although only linear equations are considered here, this technique can be easily extended to non-linear equations, while keeping its main features as for the linear case. The numerical experiments presented are quite general and confirm most of these features. These examples also show a good way of implementing this scheme.     

Singular perturbation of general boundary value problem for higher order quasilinear elliptic equation involving many small parameters

In this paper applying M. I. Visik’s and L. R. Lyuster-nik’s^[1] asymptotic method and principle of fixed point of functional analysis, we study the singular perturbation of general boundary value problem for higher order quasilinear elliptic equation in the case of boundary perturbation combined with equation perturbation. We prove the existence and uniqueness of solution for perturbed problem. We give its asymptotic approximation and estimation of related remainder term.     

The motion of organization center of scroll waves in excitable media with single diffusion

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎ１９７３,Ｗｉｎｆｒｅｆｉｒｓｔｄｉｓｃｏｖｅｒｅｄｔｈｅｔｈｒｅ＿ｄｉｍｅｎｓｉｏｎａｌｓｃｒｏｌｗａｖｅｉｎｔｈｅＢｅｌｏｕｓｏｖ＿Ｚｈａｂｏｔｉｎｓｋｙｒｅａｇｅｎｔ．Ｉｎ１９８３,Ｗｅｌｓｈｓｈｏｗｅｄｔｈｅｅｘｉ...     

The application of finite element analysis to the solution of Stokes wave diffraction problems

A new finite-element based method of calculating non-linear wave loads on offshore structures in extreme seas is presented in this paper. The diffraction wave field is modelled using Stokes wave theory developed to second order. Wave loads and free surface elevations are obtained for fixed surface-piercing structures by solving a boundary value problem for the second-order velocity potential. Special attention has been given to the radiation condition for the second-order diffraction field. Results are presented for three test examples, the vertical cylinders of Kim and Yue and of Chakrabarti, and an elliptic cylinder. These results demonstrate that early problems with the application of second-order theory arising from inadequate radiation conditions have been overcome.     

A uniform high-order method for a singular perturbation problem in conservative form

A uniform high order method is presented for the numerical solution of a singular perturbation problem in conservative form. We firest replace the original second-order problem (1.1) by two equivalent first-order problems (1.4), i.e., the solution of (1.1) is a linear combination of the solutions of (1.4). Then we derive a uniformly O(h~m+1)accurate scheme for the first-order problems (1.4), where m is an arbitrary nonnegative integer, so we can get a uniformly O(h~m+1) accurate solution of the original problem (1.1) by relation (1.3). Some illustrative numerical results are also given.     

Resonant non-linear waves—III. Elastic continuum with quadratic non-linearity

The paper derives the relevant non-linear integro-differential evolution equation by the method due to Collins expanding on a procedure by Keller. The quadratically non-linear case is not a trivial variation over the cubically non-linear case that was presented in preceding papers. As expected a different scaling and ordering of terms is needed and the first order perturbation solution provides no information on resonance. Nevertheless, although obtained by much longer calculations, the final equation for the present case is of identical form, with differences only in numerical coefficients, with the cubic case that was presented and solved earlier.     

Steady flows of non-Newtonian fluids past a porous plate with suction or injection

The problem of the steady flow of three classes of non-linear fluids of the differential type past a porous plate with uniform suction or injection is studied. The flow which is studied is the counterpart of the classical ‘asymptotic suction’ problem, within the context of the non-Newtonian fluid models. The non-linear differential equations resulting from the balance of momentum and mass, coupled with suitable boundary conditions, are solved numerically either by a finite difference method or by a collocation method with a B-spline function basis. The manner in which the various material parameters affect the structure of the boundary layer is delineated. The issue of paucity of boundary conditions for general non-linear fluids of the differential type, and a method for augmenting the boundary conditions for a certain class of flow problems, is illustrated. A comparison is made of the numerical solutions with the solutions from a regular perturbation approach, as well as a singular perturbation.     

An analytical solution for the laminar flow over a flat plate with similarity preserving suction

An exact analytical solution is presented for the laminar boundary-layer flow over a semi-infinite flat plate subjected to a type of similarity preserving suction. The solution is developed for the case of a plate immersed in either a uniform compressible stream with viscosity proportional to temperature or a uniform incompressible stream with constant viscosity. The problem is formulated in Crocco's variables. It is described by a second-order, non-linear, ordinary differential (and singular) boundary-value problem for the shear stress as a function of the velocity in the boundary layer. A unique solution is shown to exist and to possess a power series representation for all magnitudes of suction. The series is constructed explicitly and provides a transcendental equation for the shear stress at the plate (the important skin friction) which can be solved to any desired accuracy. Examples of upper and lower bounds for the wall shear are presented for several magnitudes of suction and confirm the reasonable accuracy of results obtained heretofore only by numerical solutions of the problem. In addition to the intrinsic value of the technique developed, it can be the basis of accurate checks for the numerical solution of more complex problems.     

The difference methods for the solution of singular-perturbation for the elliptic-parabolic differential equation

In this paper it is discussed the difference method for the solution of singular perturbation problems for the elliptic equations, involving small parameter in the higher derivatives. As ε= 0 the original equations are degenerated into the parabolic equations.Authors constructed special difference schemes by means of the boundary layer properties of the solutions of these problems as well as investigated the convergence of this scheme and asymptotic behaviour of the solutions. Finally, a numerical example is given.     

Periodic solutions of strongly non-linear Mathieu oscillators

The purpose of this paper is to continue our investigation into periodic solutions of strongly non-linear Mathieu oscillators. The modified version of the generalized averaging method which we developed recently is applied to derive highly accurate analytical expressions for these periodic solutions. These analytical results are used, together with the perturbation methods of multiple time scaling, to obtain second order expressions for the stability regions of these periodic solutions. The analytical research results are verified with numerical computations. Very good agreement is found, which shows the applicability of the modified version of the generalized averaging method to this kind of strongly non-linear oscillators. These oscillators may be used to model the beam-beam interaction in particle accelerators.     

Non-linear elastic analysis of thin rods subjected to bending with arbitrary kinetic conditions of their cross-sections

In this paper problems concerning the non-linear analysis of thin rods due to pure bending with constant initial curvatures and twist and with arbitrary kinetic conditions of their cross-sections are presented. Couples are not considered as being applied to the rods except at their ends. The solutions developed in this paper, which determine the curvature components and the twist of the rod after deformation, are exact in the form of elliptic integrals.     

An analytical approximate technique for a class of strongly non-linear oscillators

An analytical approximate technique for large amplitude oscillations of a class of conservative single degree-of-freedom systems with odd non-linearity is proposed. The method incorporates salient features of both Newton's method and the harmonic balance method. Unlike the classical harmonic balance method, accurate analytical approximate solutions are possible because linearization of the governing differential equation by Newton's method is conducted prior to harmonic balancing. The approach yields simple linear algebraic equations instead of non-linear algebraic equations without analytical solution. With carefully constructed iterations, only a few iterations can provide very accurate analytical approximate solutions for the whole range of oscillation amplitude beyond the domain of possible solution by the conventional perturbation methods or harmonic balance method. Three examples including cubic-quintic Duffing oscillators are presented to illustrate the usefulness and effectiveness of the proposed technique.     

Adaptive and anisotropic mesh strategy for thin shell problems. Case of inhibited parabolic shells

The aim of this paper is to show the reliability of an adaptive and anisotropic mesh procedure for thin shell problems. We consider singular perturbation problems only for parabolic shells whose behavior is described by the Koiter model. The corresponding system of equations, which depends on the relative thickness ε of the shell, is elliptic except at the limit for ε = 0 where it is parabolic. In a first part of this paper, we study theoretically the phenomena of internal layers appearing during the singular perturbation process, when the loading is somewhat singular. These layers have very different structures either they are along or across the asymptotic lines of the middle surface of the shell. In a second part, numerical computations are performed using a finite element software coupled with an adaptive anisotropic mesh generator. This technique enables to approach accurately the singularities and the layers predicted by the theory especially for very small values of the thickness. The efficiency of such a procedure in comparison with uniform meshes is put in a prominent position.     

A perturbation solution for non-Newtonian fluid two-phase flow through Porous media

In this paper, the mathematical problem of weak non-Newtonian fluid two-phase flow through porous media, including the effect of capillary pressure, is solved by singular perturbation method in combination with regular perturbation method. The asymptotic analytical solutions of the fractional flow function and the wetting phase saturation are obtained. The results are verified by numerical calculations and by classical solutions for corresponding Newtonian case. The influences of the non-Newtonian exponent and capillary pressure are discussed.     

An Elliptic Lindstedt--Poincaré Method for Certain Strongly Non-Linear Oscillators

An elliptic Lindstedt--Poincaré (L--P) method is presented for the steady-state analysis of strongly non-linear oscillators of the form , in which the Jacobian elliptic functions are employed instead of the usual circular functions in the classical L--P perturbation procedure. This method can be viewed as a generalization of the L--P method. As an application of this method, three types of the generalized Van der Pol equation with are studied in detail.