Transformation of an axialsymmetric disk problem for the helmholtz equation into an ordinary differential equation

The problem of solving the three-dimensional Helmholtz equation in the exterior of a circular disk is considered where radially symmetric Dirichlet data on the disk are assumed to be prescribed. This problem for example arises in the scattering of plane (sound) waves at an infinite plane screen with a circular aperture if the direction of the incident wave is normal to the screen, as well as in the process of diffusion through a circular hole. By applying the factorization technique developed in [N. GORENFLO, M. WERNER,Solution of a finite convolution equation with a Hankel kernel by matrix factorization, SIAM J. Math. Anal., 28 (1997), pp. 434–451] to the disk problem an equivalent ordinary differential equation is derived, whose solution leads directly to the solution of the disk problem. This differential equation belongs to a class of ordinary differential equations which are of higher complexity than the standard ordinary differential equations of mathematical physics. The examination of this new class of differential equations therefore is motivated.     

On double reductions from symmetries and conservation laws for a damped Boussinesq equation

In this work, we study a Boussinesq equation with a strong damping term from the point of view of the Lie theory. We derive the classical Lie symmetries admitted by the equation as well as the reduced ordinary differential equations. Some nontrivial conservation laws are derived by using the multipliers method. Taking into account the relationship between symmetries and conservation laws and applying the double reduction method, we obtain a direct reduction of order of the ordinary differential equations and in particular a kink solution.     

Be careful with the Exp-function method – Additional remarks

It is shown that the solution produced by the Exp-function method may not hold for all initial conditions. Riccati and Maccari nonlinear differential equations are used to illustrate that fact. Conditions of existence for the produced solution in the space of initial conditions and in the space of system’s parameters are derived using the operator method based on the generalized operator of differentiation. The concept of the expansion of an ordinary differential equation is introduced and it is shown that the algebraic–analytical solution of Maccari equation can be produced by solving Riccati equation.     

有限挠度下Timoshenko梁中的非线性弯曲波及其混沌行为

以Timoshenko梁理论为基础,引入了有限挠度和轴向惯性,建立了支配梁运动的非线性偏微分方程组,采用行波法求解,通过某些积分技巧,将其转化为一个非线性常微分方程．常微分方程的定性分析表明,在一定条件下,系统存在异宿轨道,预示着有冲击波解存在．借助Jacobi椭圆函数展开求解,得到了非线性波动方程的准确周期解及其当模数m→1退化情况下的冲击波解．进而考虑阻尼和外加横向载荷对系统的摄动,利用Melnikov函数给出了横截异宿点出现的阈值条件,从而表明系统具有Smale马蹄意义下的混沌性质．     

The Sinc-collocation method for solving the Thomas–Fermi equation

A numerical technique for solving nonlinear ordinary differential equations on a semi-infinite interval is presented. We solve the Thomas–Fermi equation by the Sinc-Collocation method that converges to the solution at an exponential rate. This method is utilized to reduce the nonlinear ordinary differential equation to some algebraic equations. This method is easy to implement and yields very accurate results.     

Similarity solutions to nonlinear heat conduction and Burgers/Korteweg–deVries fractional equations

We analyze self-similar solutions to a nonlinear fractional diffusion equation and fractional Burgers/Korteweg–deVries equation in one spatial variable. By using Lie-group scaling transformation, we determined the similarity solutions. After the introduction of the similarity variables, both problems are reduced to ordinary nonlinear fractional differential equations. In two special cases exact solutions to the ordinary fractional differential equation, which is derived from the diffusion equation, are presented. In several other cases the ordinary fractional differential equations are solved numerically, for several values of governing parameters. In formulating the numerical procedure, we use special representation of a fractional derivative that is recently obtained.     

二阶线性常微分方程的两点边值问题的泰勒展开式解法

本文用泰勒展开公式求解二阶线性常微分方程的两点边值问题.首先将两点边值问题化为一个F redho lm积分方程,进一步通过泰勒展开公式化F redho lm积分方程为线性方程组,利用G ramm er法则可求得问题的近似解.     

Laplace transform for the solution of higher order deformation equations arising in the homotopy analysis method

A new analytic approach for solving nonlinear ordinary differential equations with initial conditions is proposed. First, the homotopy analysis method is used to transform a nonlinear differential equation into a system of linear differential equations; then, the Laplace transform method is applied to solve the resulting linear initial value problems; finally, the solutions to the linear initial value problems are employed to form a convergent series solution to the given problem. The main advantage of the new approach is that it provides an effective way to solve the higher order deformation equations arising in the homotopy analysis method.     

Singularity analysis of a variant of the Painlevé–Ince equation

We examine by singularity analysis an equation derived by reduction using Lie point symmetries from the Euler–Bernoulli Beam equation which is the Painlevé–Ince Equation with additional terms. The equation possesses the same leading-order behaviour and resonances as the Painlevé–Ince Equation and has a Right Painlevé Series. However, it has no Left Painlevé Series. A conjecture for the existence of Left Painlevé Series for ordinary differential equations is given.     

An indirect variable transformation approach and Jacobi elliptic solutions to Korteweg de Vries equation

Based on a variable change and the variable separated ODE method, an indirect variable transformation approach is proposed to search exact solutions to special types of partial differential equations (PDEs). The new method provides a more systematical and convenient handling of the solution process for the nonlinear equations. Its key point is to reduce the given PDEs to variable-coefficient ordinary differential equations, then we look for solutions to the resulting equations by some methods. As an application, exact solutions for the KdV equation are formally derived.     

Complex Lie symmetries for scalar second-order ordinary differential equations

A scalar complex ordinary differential equation can be considered as two coupled real partial differential equations, along with the constraint of the Cauchy–Riemann equations, which constitute a system of four equations for two unknown real functions of two real variables. It is shown that the resulting system possesses those real Lie symmetries that are obtained by splitting each complex Lie symmetry of the given complex ordinary differential equation. Further, if we restrict the complex function to be of a single real variable, then the complex ordinary differential equation yields a coupled system of two ordinary differential equations and their invariance can be obtained in a non-trivial way from the invariance of the restricted complex differential equation. Also, the use of a complex Lie symmetry reduces the order of the complex ordinary differential equation (restricted complex ordinary differential equation) by one, which in turn yields a reduction in the order by one of the system of partial differential equations (system of ordinary differential equations). In this paper, for simplicity, we investigate the case of scalar second-order ordinary differential equations. As a consequence, we obtain an extension of the Lie table for second-order equations with two symmetries.     

A Note on the Numerical Solution of Quasi-linear Parabolic Partial Differential Equations with Error Estimates

A method is considered for the numerical solution of quasi-linearpartial differential equations. The partial differential equationis reduced to a set of ordinary differential equations usinga Chebyshev series expansion. The exact solution of this setof ordinary differential equations is shown to be the solutionof a perturbed form of the original equation. This enables errorestimates to be found for linear and mildly non-linear problems.     

Efficient hybrid method for solving special type of nonlinear partial differential equations

In this article, an efficient hybrid method has been developed for solving some special type of nonlinear partial differential equations. Hybrid method is based on tanh–coth method, quasilinearization technique and Haar wavelet method. Nonlinear partial differential equations have been converted into a nonlinear ordinary differential equation by choosing some suitable variable transformations. Quasilinearization technique is used to linearize the nonlinear ordinary differential equation and then the Haar wavelet method is applied to linearized ordinary differential equation. A tanh–coth method has been used to obtain the exact solutions of nonlinear ordinary differential equations. It is easier to handle nonlinear ordinary differential equations in comparison to nonlinear partial differential equations. A distinct feature of the proposed method is their simple applicability in a variety of two‐ and three‐dimensional nonlinear partial differential equations. Numerical examples show better accuracy of the proposed method as compared with the methods described in past. Error analysis and stability of the proposed method have been discussed.     

Numerical solutions for convection–diffusion equation using El-Gendi method

The method of El-Gendi [El-Gendi SE. Chebyshev solution of differential integral and integro-differential equations. J Comput 1969;12:282–7; Mihaila B, Mihaila I. Numerical approximation using Chebyshev polynomial expansions: El-gendi’s method revisited. J Phys A Math Gen 2002;35:731–46] is presented with interface points to deal with linear and non-linear convection–diffusion equations.The linear problem is reduced to two systems of ordinary differential equations. And, then, each system is solved using three-level time scheme.The non-linear problem is reduced to three systems of ordinary differential. Each one of these systems is, then, solved using three-level time scheme. Numerical results for Burgers’ equation and modified Burgers’ equation are shown and compared with other methods. The numerical results are found to be in good agreement with the exact solutions.     

A method of solving a nonlinear boundary value problem with a parameter for a loaded differential equation

A nonlinear loaded differential equation with a parameter on a finite interval is studied. The interval is partitioned by the load points, at which the values of the solution to the equation are set as additional parameters. A nonlinear boundary value problem for the considered equation is reduced to a nonlinear multipoint boundary value problem for the system of nonlinear ordinary differential equations with parameters. For fixed parameters, we obtain the Cauchy problems for ordinary differential equations on the subintervals. Substituting the values of the solutions to these problems into the boundary condition and continuity conditions at the partition points, we compose a system of nonlinear algebraic equations in parameters. A method of solving the boundary value problem with a parameter is proposed. The method is based on finding the solution to the system of nonlinear algebraic equations composed.     

A remark on connection between conjugacy of half-linear differential equation and equation with mixed nonlinearities

In this paper new criteria for conjugacy of half-linear ordinary differential equations are derived by using a Riccati transformation. These criteria are used to derive nonexistence and oscillation results for an equation with mixed nonlinearities, which is viewed as a perturbation of a half-linear equation.     

Symmetry Analysis of Differential Equations Using MATHLIE

A general procedure for solving ordinary differential equations of arbitrary order is discussed. The method used is based on symmetries of a differential equation. The known symmetries allow the derivation of first integrals of the equation. The knowledge of at least r symmetries of an ordinary differential equation of order n with r n is the basis for deriving the solution. Our aim is to show that Lie's theory is a useful tool for solving ordinary differential equations of higher orders. Bibliography: 12 titles.     

Boundary element solution of the two-dimensional groundwater flow equation with stochastic free surface boundary condition

An innovative approach to the approximate solution of stochastic partial differential equations in groundwater flow is presented. The method uses a formulation of the Ito's lemma in Hilbert spaces to derive partial differential equations satisfying the moments of the solution process. Since the moments equations are deterministic, they could be solved by any analytical or numerical method existing in the literature. This permits the analysis and solution of stochastic partial differential equations occurring in two-dimensional or three-dimensional domains of any geometrical shape. The method is tested for the first time in the present paper through a practical application in a sandy phreatic aquifer at the Chalk River Nuclear Laboratories, Ontario, Canada. The equation solved is the two-dimensional LaPlace equation with a dynamic, randomly perturbed, free surface boundary condition. The moments equations are derived and solved by using the boundary integral equation method. A comparison is made with a previous analytical solution obtained by applying the randomly forced one-dimensional Boussinesq equation, and some observations on modeling procedures are given.     

Error bounds for numerical solutions of ordinary differential equations

Summary An a posteriori error bound, for an approximate solution of a system of ordinary differential equations, is derived as the solution of a Riccati equation. The coefficients of the Riccati equation depend on an eigenvalue of a matrix related to a Jacobian matrix, on a Lipschitz constant for the Jacobian matrix, and on the approximation defect. An upper bound is computable as the formal solution of a sequence of Riccati equations with constant coefficients. This upper bound may sometimes be used to control step length in a numerical method.     

Three dimensional elastodynamics of 2D quasicrystals: The derivation of the time-dependent fundamental solution

The time-dependent differential equations of elasticity for 2D quasicrystals with general structure of anisotropy (dodecagonal, octagonal, decagonal, pentagonal, hexagonal, triclinic) are considered in the paper. These equations are written in the form of a vector partial differential equation of the second order with symmetric matrix coefficients. The fundamental solution (matrix) is defined for this vector partial differential equation. A new method of the numerical computation of values of the fundamental solution is suggested. This method consists of the following: the Fourier transform with respect to space variables is applied to vector equation for the fundamental solution. The obtained vector ordinary differential equation has matrix coefficients depending on Fourier parameters. Using the matrix computations a solution of the vector ordinary differential equation is numerically computed. Finally, applying the inverse Fourier transform numerically we find the values of the fundamental solution. Computational examples confirm the robustness of the suggested method for 2D quasicrystals with arbitrary type of anisotropy.