Propagation of surface (seismic) waves: Ordinary differential equations with strongly nonlinear damping

Second-order ordinary differential equations (ODEs) with strongly nonlinear damping (cubic nonlinearities) govern surface wave motions that entail nonlinear surface seismic motions. They apply to dynamic crack propagation and nonlinear oscillation problems in physics and nonlinear mechanics. It is shown that the nonlinear surface seismic wave equation (Rayleigh equation) admits several functional transformations and it is possible to reduce it to an equivalent first-order Abel ODE of the second kind in normal form. Based on a recently developed methodology concerning the construction of exact analytic solutions for the type of Abel equations under consideration, exact solutions are obtained for the nonlinear seismic wave (NLSW) equation for initial conditions of the physical problem. The method employed is general and can be applied to a large class of relevant ODEs in mathematical physics and nonlinear mechanics.     

ADAPTIVE INTERVAL WAVELET PRECISE INTEGRATION METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS

IntroductionThepreciseintegrationmethod(PIM) [1],whichwasproposedforsolvingstructuraldynamicequations.Thismethodissimplerandpossesseshigherprecision .Forlinearsteadystructuraldynamicsystems,itsnumericalresultsattheintegrationpointsarealmostequaltothatoftheexactsolutioninmachineaccuracy .InthepreciseintegrationmethodforsolvingPDEs,theequationsshouldbediscretizedinthephysicalspaceforobtainingthesystemofODEsintime ,whichisoftenexecutedbythefinitedifferencemethodorthefiniteelementmethod .Inrec…     

Asymptotic investigation of a weakly nonlinear multifrequency oscillating system

We consider a weakly nonlinear multifrequency autonomous system of differential equations with a small parameter on the right-hand side under the condition that the unperturbed system has a quasiperiodic general solution. The system is reduced to a simpler form by averaging and separation. We establish sufficient conditions for the preservation of an invariant torus under a small perturbation.     

Effect of Viscous Dissipation on the Boundary Layer Flow and Heat Transfer Past a Permeable Stretching Surface Embedded in a Porous Medium with a Second-Order Slip Using Chebyshev Finite Difference Method

This article investigates a theoretical and numerical study for the effect of viscous dissipation on the steady flow with heat transfer of Newtonian fluid toward a permeable stretching surface embedded in a porous medium with a second-order slip and thermal slip. The governing nonlinear partial differential equations are converted into nonlinear ordinary differential equations (ODEs) using similarity variables. The resulting ODEs are successfully solved numerically with the help of Chebyshev finite difference method. Graphically results are shown for non-dimensional velocities and temperature. The effects of the porous parameter, the suction (injection) parameter, Eckert number, first- and second-order velocity slip parameter, the thermal slip parameter and the Prandtl number on the flow and temperature profiles are presented. Moreover, the local skin-friction and Nusselt numbers are presented. A comparison of numerical results is made with the earlier published results under limiting cases.     

Nonlinear theory of hydrodynamic stability and bifurcations of solutions of the Navier-Stokes equations

A new approach to the investigation of the nonlinear development of disturbances, based on the theory of invariant manifolds, is outlined. A method of obtaining projections of the Navier-Stokes equations on finite-dimensional invariant manifolds is proposed. The behavior of the disturbances is finally described by a system of ordinary differential equations with right sides in the form of power series in the amplitudes. It is an important property of the method that these series are convergent. A two-dimensional invariant projection is calculated numerically for plane Poiseuille flow. As a result, it is possible to clarify the nature of the change from subcritical to supercritical bifurcation and investigate new bifurcations of periodic flow regimes.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 9–15, January–February, 1990.The author is grateful to A. Zharilkasinov for carrying out the numerical calculations.     

Group classification and exact solutions of a radially symmetric porous-medium equation

In this paper we perform a group classification for the generalized radial porous-medium equation. We also classify symmetry reductions of the equation to first- or second-order ordinary differential equations (ODEs) and hence construct invariant solutions in a systematic manner. We show that the reduced second-order equations are invariant under either a two-parameter or one-parameter Lie groups. In the first case, they are completely integrated by a pair of quadratures. In the latter, they are often reduced to first-order ODEs of Abel type.     

Invariant manifolds of one class of systems of impulsive differential equations

We consider general problems related to the existence of invariant toroidal sets for linear and weakly nonlinear systems of impulsive differential equations defined in the direct product of an m-dimensional torus and an n-dimensional Euclidean space. We investigate classes of problems for which the conditions for the existence of invariant toroidal manifolds are satisfied.     

Existence of an invariant torus for one class of systems of differential equations

We consider the problem of the existence of an asymptotically stable toroidal set for a system of linear differential equations defined on an m-dimensional torus. We establish conditions under which a nonlinear system of differential equations has an invariant toroidal manifold. Translated from Neliniini Kolyvannya, Vol. 11, No. 4, pp. 520–529, October–December, 2008.     

Van der Pol type self-excited micro-cantilever probe of atomic force microscopy

A technique for dimensional reduction of nonlinear delay differential equations (DDEs) with time-periodic coefficients is presented. The DDEs considered here have a canonical form with at most cubic nonlinearities and periodic coefficients. The nonlinear terms are multiplied by a perturbation parameter. Perturbation expansion converts the nonlinear response problem into solutions of a series of nonhomogeneous linear ordinary differential equations (ODEs) with time-periodic coefficients. One set of linear nonhomogeneous ODEs is solved for each power of the perturbation parameter. Each ODE is solved by a Chebyshev spectral collocation method. Thus we compute a finite approximation to the nonlinear infinite-dimensional map for the DDE. The linear part of the map is the monodromy operator whose eigenvalues characterize stability. Dimensional reduction on the map is then carried out. In the case of critical eigenvalues, this corresponds to center manifold reduction, while for the noncritical case resonance conditions are derived. The accuracy of the nonlinear Chebyshev collocation map is demonstrated by finding the solution of a nonlinear delayed Mathieu equation and then a milling model via the method of steps. Center manifold reduction is illustrated via a single inverted pendulum including both a periodic retarded follower force and a nonlinear restoring force. In this example, the amplitude of the limit cycle associated with a flip bifurcation is found analytically and compared to that obtained from direct numerical simulation. The method of this paper is shown by example to be applicable to systems with strong parametric excitations.     

Observer design and output feedback stabilization for nonlinear multivariable systems with diffusion PDE-governed sensor dynamics

This paper addresses the problems of observer design and output feedback stabilization for a class of nonlinear multivariable systems, where the nonlinear system dynamics are described by ordinary differential equations (ODEs), and the sensor dynamics are governed by diffusion partial differential equations (PDEs). Based on the Luenberger observer theory, a Luenberger-type PDE-ODE cascaded observer is derived to estimate the state variables of the system. Then, an observer-based output feedback stabilizing controller is developed. The exponential stability of both the observer error system and closed-loop control system is proven via the Lyapunov direct method. Finally, numerical examples are provided to illustrate the effectiveness of the proposed design methods.     

ARC-length method for differential equations

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｏｒｄｉｎａｒｙａｎｄｐａｒｔｉａｌｄｉｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｓｏｆｃｏｎｔｉｎｕｕｍｐｒｏｂｌｅｍａｒｅｏｆｔｅｎｗｉｔｈｃｅｒｔａｉｎｔｙｐｅｓｏｆｓｉｎｇｕｌａｒｉｔｙａｓｓｔｉｆｐｒｏｐｅｒｔｙ,ｏｒ...     

插值矩阵法分析双材料平面V形切口奇异阶

对二维V形切口问题提出奇异阶分析的一个新方法.首先,以V形切口尖端附近位移场沿其径向渐近展开为基础,将其线弹性理论控制方程转换成切口尖端附近关于周向变量的常微分方程组特征值问题,然后将数值求解两点边值问题的插值矩阵法进一步拓展为求解一般常微分方程组特征值问题,插值矩阵法是在离散节点上采用微分方程中待求函数的最高阶导数作为基本未知量.由此,V形切口的应力奇性阶问题通过插值矩阵法获得,同时相应的切口附近位移场和应力场特征向量一并求出.     

DYNAMICAL APPROACH STUDY OF SPURIOUS STEADY-STATE NUMERICAL SOLUTIONS OF NONLINEAR DIFFERENTIAL EQUATIONS II. GLOBAL ASYMPTOTIC BEHAVIOR OF TIME DISCRETIZATIONS ∗

SUMMARY

The global asymptotic nonlinear behavior of 11 explicit and implicit time discretizations for four 2 × 2 systems of first-order autonomous nonlinear ordinary differential equations (ODEs) is analyzed. The objectives are to gain a basic understanding of the difference in the dynamics of numerics between the scalars and systems of nonlinear autonomous ODEs and to set a baseline global asymptotic solution behavior of these schemes for practical computations in computational fluid dynamics. We show how “numerical” basins of attraction can complement the bifurcation diagrams in gaining more detailed global asymptotic behavior of time discretizations for nonlinear differential equations (DCs). We show how in the presence of spurious asymptotes the basins of the true stable steady states can be segmented by the basins of the spurious stable and unstable asymptotes. One major consequence of this phenomenon which is not commonly known is that this spurious behavior can result in a dramatic distortion and, in most cases, a dramatic shrinkage and segmentation of the basin of attraction of the true solution for finite time steps. Such distortion, shrinkage and segmentation of the numerical basins of attraction will occur regardless ofthe stability ofthe spurious asymptotes, and will occur for unconditionally stable implicit linear multistep methods. In other words, for the same (common) steady-state solution the associated basin of attraction of the DE might be very different from the discretized counterparts and the numerical basin of attraction can be very different from numerical method to numerical method. The results can be used as an explanation for possible causes of error, and slow convergence and nonconvergence of steady-state numerical solutions when using the time-dependent approach for nonlinear hyperbolic or parabolic PDEs.     

METHOD–OF–LINES SOLUTION OF TIME–DEPENDENT TWO–DIMENSIONAL NAVIER–STOKES EQUATIONS

A novel approach to the development of a code for the solution of the time-dependent two-dimensional Navier–Stokes equations is described. The code involves coupling between the method of lines (MOL) for the solution of partial differential equations and a parabolic algorithm which removes the necessity of iterative solution on pressure and solution of a Poisson-type equation for the pressure. The code is applied to a test problem involving the solution of transient laminar flow in a short pipe for an incompressible Newtonian fluid. Comparisons show that the MOL solutions are in good agreement with the previously reported values. The proposed method described in this paper demonstrates the ease with which the Navier–Stokes equations can be solved in an accurate manner using sophisticated numerical algorithms for the solution of ordinary differential equations (ODEs).     

On a computational approach for the approximate dynamics of averaged variables in nonlinear ODE systems: Toward the derivation of constitutive laws of the rate type

A non-perturbative approach to the time-averaging of nonlinear, autonomous ordinary differential equations is developed based on invariant manifold methodology. The method is implemented computationally and applied to model problems arising in the mechanics of solids.     

Heat transfer characteristics of thin power-law liquid films over horizontal stretching sheet with internal heating and variable thermal coefficient

The effect of internal heating source on the film momentum and thermal transport characteristic of thin finite power-law liquids over an accelerating unsteady horizontal stretched interface is studied. Unlike most classical works in this field, a general surface temperature distribution of the liquid film and the generalized Fourier’s law for varying thermal conductivity are taken into consideration. Appropriate similarity transformations are used to convert the strongly nonlinear governing partial differential equations (PDEs) into a boundary value problem with a group of two-point ordinary differential equations (ODEs). The correspondence between the liquid film thickness and the unsteadiness parameter is derived with the BVP4C program in MATLAB. Numerical solutions to the self-similarity ODEs are obtained using the shooting technique combined with a Runge-Kutta iteration program and Newton’s scheme. The effects of the involved physical parameters on the fluid’s horizontal velocity and temperature distribution are presented and discussed.     

Simultaneous impacts of MHD and variable wall temperature on transient mixed Casson nanofluid flow in the stagnation point of rotating sphere

A numerical analysis is provided to scrutinize time-dependent magnetohydrodynamics(MHD) free and forced convection of an electrically conducting non-Newtonian Casson nanofluid flow in the forward stagnation point region of an impulsively rotating sphere with variable wall temperature. A single-phase flow of nanofluid model is reflected with a number of experimental formulae for both effective viscosity and thermal conductivity of nanofluid. Exceedingly nonlinear governing partial differential equations(PDEs)subject to their compatible boundary conditions are mutated into a system of nonlinear ordinary differential equations(ODEs). The derived nonlinear system is solved numerically with implementation of an implicit finite difference procedure merging with a technique of quasi-linearization. The controlled parameter impacts are clarified by a parametric study of the entire flow regime. It is depicted that from all the exhibited nanoparticles,Cu possesses the best convection. The surface heat transfer and surface shear stresses in the x-and z-directions are boosted with maximizing the values of nanoparticle solid volume fraction ? and rotation λ. Besides, as both the surface temperature exponent n and the Casson parameter γ upgrade, an enhancement of the Nusselt number is given.     

Nonlinear nonplanar dynamics of a parametrically excited inextensional elastic beam

The nonlinear dynamics of a clamped-clamped/sliding inextensional elastic beam subject to a harmonic axial load is investigated. The Galerkin method is used on the coupled bending-bending-torsional nonlinear equations with inertial and geometric nonlinearities and the resulting two second order ordinary differential equations are studied by the method of multiple time seales and by direct numerical integration. The amplitude equations are analyzed for steady and Hopf bifurcations. Depending on the amplitude of excitation, the damping and the ratio of principal flexural rigidities, various qualitatively distinct frequency response diagrams are uncovered and limit cycles and chaotic motions are found. In the truncated two-degree-of-freedom system the transition from periodic to chaotic amplitude-modulated motions is via the process of torus doubling and subsequent destruction of the torus.     

Constructing invariant tori for two weakly coupled van der Pol oscillators

An algorithm is developed for the construction of an invariant torus of a weakly coupled autonomous oscillator. The system is put into angular standard form. The determining equations are found by averaging and are solved for the approximate amplitudes of the torus. A perturbation series is then constructed about the approximate amplitudes with unknown coefficients as periodic functions of the angular variables. A sequence of solvable partial differential equations is developed for determining the coefficients. The algorithm is applied to a system of nonlinearly coupled van der Pol equations and the first order coefficients are generated in a straightforward manner. The approximation shows both good numerical accuracy and reproducibility of the periodicities of the van der Pol system. A comparitive analysis of integrating the van der Pol system with integrating the phase equations from the angular standard form on the approximate torus shows numerical errors of the order of the perturbation parameter =0.05 for integrations of up to 10,000 steps. Applying FFT to the numerical periodicities generated by integrating the van der Pol system near the tours reveals the same predominant frequencies found in the perturbation coefficients. Finally an expected rotation number is found by integrating the phase equations on the approximate torus.Contribution of the National Institute of Standards and Technology, a Federal agency.     

Postbuckling and mode jumping analysis of composite laminates using an asymptotically correct, geometrically non-linear theory

An analytic method is presented in this paper to study the postbuckling and mode jumping behavior of bi-axially compressed composite laminates. The governing partial differential equations (PDEs) are derived rigorously from an asymptotically correct, geometrically non-linear theory. A novel and relatively simpler solution approach is developed to solve the two coupled fourth-order PDEs, namely, the compatibility equation and the dynamic governing equation. The generalized Galerkin method is used to solve boundary value problems corresponding to antisymmetric angle-ply and cross-ply composite plates, respectively. The variety of possible modal interactions is expressed in an explicit and concise form by transforming the coupled non-linear governing equations into a system of non-linear ordinary differential equations (ODEs).

The comparison between the present method and the finite element analysis (FEA) shows a pretty good match in their numerical results in the primary postbuckling region. While the FEA may lose its convergence when solution comes close to the secondary bifurcation point, the analytic approach has the capability of exploring deeply into the post-secondary buckling realm and capture the mode jumping phenomenon for various combinations of plate configurations and in-plane boundary conditions. Free vibration along the stable primary postbuckling and the jumped equilibrium paths are also studied.