Semi‐invariant solutions of the Navier‐Stokes equations

We present new exact solutions and reduced differential systems of the Navier‐Stokes equations of incompressible viscous fluid flow. We apply the method of semi‐invariant manifolds, introduced earlier as a modification of the Lie invariance method. We show that many known solutions of the Navier‐Stokes equations are, in fact, semi‐invariant and that the reduced differential systems we derive using semi‐invariant manifolds generalize previously obtained results that used ad hoc methods. Many of our semi‐invariant solutions solve decoupled systems in triangular form that are effectively linear. We also obtain several new reductions of Navier‐Stokes to a single nonlinear partial differential equation. In some cases, we can solve reduced systems and generate new analytic solutions of the Navier‐Stokes equations or find their approximations, and physical interpretation.     

Oscillation result for half‐linear dynamic equations on timescales and its consequences

We study oscillatory properties of half‐linear dynamic equations on timescales. Via the combination of the Riccati technique and an averaging method, we find the domain of oscillation for many equations. The presented main result is not the conversion of a known result from the theory of differential or difference equations, ie, we obtain new results for the timescales (for differential equations) and (for difference equations). Half‐linear equations generalize linear equations (in fact, they coincide with certain one‐dimensional PDEs with p‐Laplacian), but the main result is new also for linear differential and difference equations. The corresponding corollaries and examples are given as well.     

Computational methods of solving the boundary value problems for the loaded differential and Fredholm integro‐differential equations

The article presents a new general solution to a loaded differential equation and describes its properties. Solving a linear boundary value problem for loaded differential equation is reduced to the solving a system of linear algebraic equations with respect to the arbitrary vectors of general solution introduced. The system's coefficients and right sides are computed by solving the Cauchy problems for ordinary differential equations. Algorithms of constructing a new general solution and solving a linear boundary value problem for loaded differential equation are offered. Linear boundary value problem for the Fredholm integro‐differential equation is approximated by the linear boundary value problem for loaded differential equation. A mutual relationship between the qualitative properties of original and approximate problems is obtained, and the estimates for differences between their solutions are given. The paper proposes numerical and approximate methods of solving a linear boundary value problem for the Fredholm integro‐differential equation and examines their convergence, stability, and accuracy.     

Chebyshev polynomial approximation for high‐order partial differential equations with complicated conditions

In this article, a new method is presented for the solution of high‐order linear partial differential equations (PDEs) with variable coefficients under the most general conditions. The method is based on the approximation by the truncated double Chebyshev series. PDE and conditions are transformed into the matrix equations, which corresponds to a system of linear algebraic equations with the unknown Chebyshev coefficients, via Chebyshev collocation points. Combining these matrix equations and then solving the system yields the Chebyshev coefficients of the solution function. Some numerical results are included to demonstrate the validity and applicability of the method. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A Taylor polynomial approach for solving the most general linear Fredholm integro‐differential‐difference equations

In this study, a matrix method is developed to solve approximately the most general higher order linear Fredholm integro‐differential‐difference equations with variable coefficients under the mixed conditions in terms of Taylor polynomials. This technique reduces the problem into the linear algebraic system. The method is valid for any combination of differential, difference and integral equations. An initial value problem and a boundary value problem are also presented to illustrate the accuracy and efficiency of the method. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

Solution of high‐order linear Fredholm integro‐differential equations with piecewise intervals

In this study, a practical matrix method is presented to find an approximate solution for high‐order linear Fredholm integro‐differential equations with piecewise intervals under the initial boundary conditions in terms of Taylor polynomials. The method converts the integro differential equation to a matrix equation, which corresponds to a system of linear algebraic equations. Error analysis and illustrative examples are included to demonstrate the validity and applicability of the technique. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010 27: 1327–1339, 2011     

New Lyapunov‐type inequalities for a class of even‐order linear differential equations

In this paper, we obtain some new Lyapunov‐type inequalities for a class of even‐order linear differential equations, the results are new and generalize and improve some early results in this field.     

An iterative solver for the Oseen problem and numerical solution of incompressible Navier–Stokes equations

Incompressible unsteady Navier–Stokes equations in pressure–velocity variables are considered. By use of the implicit and semi‐implicit schemes presented the resulting system of linear equations can be solved by a robust and efficient iterative method. This iterative solver is constructed for the system of linearized Navier–Stokes equations. The Schur complement technique is used. We present a new approach of building a non‐symmetric preconditioner to solve a non‐symmetric problem of convection–diffusion and saddle‐point type. It is shown that handling the differential equations properly results in constructing efficient solvers for the corresponding finite linear algebra systems. The method has good performance for various ranges of viscosity and can be used both for 2D and 3D problems. The analysis of the method is still partly heuristic, however, the mathematically rigorous results are proved for certain cases. The proof is based on energy estimates and basic properties of the underlying partial differential equations. Numerical results are provided. Additionally, a multigrid method for the auxiliary convection–diffusion problem is briefly discussed. Copyright © 1999 John Wiley & Sons, Ltd.     

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.     

Amplitude–shape approximation as an extension of separation of variables

Separation of variables is a well‐known technique for solving differential equations. However, it is seldom used in practical applications since it is impossible to carry out a separation of variables in most cases. In this paper, we propose the amplitude–shape approximation (ASA) which may be considered as an extension of the separation of variables method for ordinary differential equations. The main idea of the ASA is to write the solution as a product of an amplitude function and a shape function, both depending on time, and may be viewed as an incomplete separation of variables. In fact, it will be seen that such a separation exists naturally when the method of lines is used to solve certain classes of coupled partial differential equations. We derive new conditions which may be used to solve the shape equations directly and present a numerical algorithm for solving the resulting system of ordinary differential equations for the amplitude functions. Alternatively, we propose a numerical method, similar to the well‐established exponential time differencing method, for solving the shape equations. We consider stability conditions for the specific case corresponding to the explicit Euler method. We also consider a generalization of the method for solving systems of coupled partial differential equations. Finally, we consider the simple reaction diffusion equation and a numerical example from chemical kinetics to demonstrate the effectiveness of the method. The ASA results in far superior numerical results when the relative errors are compared to the separation of variables method. Furthermore, the method leads to a reduction in CPU time as compared to using the Rosenbrock semi‐implicit method for solving a stiff system of ordinary differential equations resulting from a method of lines solution of a coupled pair of partial differential equations. The present amplitude–shape method is a simplified version of previous ones due to the use of a linear approximation to the time dependence of the shape function. Copyright © 2007 John Wiley & Sons, Ltd.     

Solution by quadratures of certain ordinary differential equations

This paper discusses methods that are applicable in the solution by quadratures of linear second‐order differential equations with variable coefficients. These same techniques, when applied to equations with constant coefficients, produce an extremely useful method in the teachingof ordinary differential equations.     

A shift‐adaptive meshfree method for solving a class of initial‐boundary value problems with moving boundaries in one‐dimensional domain

A new shift‐adaptive meshfree method for solving a class of time‐dependent partial differential equations (PDEs) in a bounded domain (one‐dimensional domain) with moving boundaries and nonhomogeneous boundary conditions is introduced. The radial basis function (RBF) collocation method is combined with the finite difference scheme, because, unlike with Kansa's method, nonlinear PDEs can be converted to a system of linear equations. The grid‐free property of the RBF method is exploited, and a new adaptive algorithm is used to choose the location of the collocation points in the first time step only. In fact, instead of applying the adaptive algorithm on the entire domain of the problem (like with other existing adaptive algorithms), the new adaptive algorithm can be applied only on time steps. Furthermore, because of the radial property of the RBFs, the new adaptive strategy is applied only on the first time step; in the other time steps, the adaptive nodes (obtained in the first time step) are shifted. Thus, only one small system of linear equations must be solved (by LU decomposition method) rather than a large linear or nonlinear system of equations as in Kansa's method (adaptive strategy applied to entire domain), or a large number of small linear systems of equations in the adaptive strategy on each time step. This saves a lot in time and memory usage. Also, Stability analysis is obtained for our scheme, using Von Neumann stability analysis method. Results show that the new method is capable of reducing the number of nodes in the grid without compromising the accuracy of the solution, and the adaptive grading scheme is effective in localizing oscillations due to sharp gradients or discontinuities in the solution. The efficiency and effectiveness of the proposed procedure is examined by adaptively solving two difficult benchmark problems, including a regularized long‐wave equation and a Korteweg‐de Vries problem. © 2016Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1622–1646, 2016     

Solvable systems of linear differential equations

The asymptotic iteration method (AIM) is an iterative technique used to find exact and approximate solutions to second-order linear differential equations. In this work, we employed AIM to solve systems of two first-order linear differential equations. The termination criteria of AIM will be re-examined and the whole theory is re-worked in order to fit this new application. As a result of our investigation, an interesting connection between the solution of linear systems and the solution of Riccati equations is established. Further, new classes of exactly solvable systems of linear differential equations with variable coefficients are obtained. The method discussed allow to construct many solvable classes through a simple procedure.     

A new method for exact solutions of variant types of time‐fractional Korteweg‐de Vries equations in shallow water waves

The current article devoted on the new method for finding the exact solutions of some time‐fractional Korteweg–de Vries (KdV) type equations appearing in shallow water waves. We employ the new method here for time‐fractional equations viz. time‐fractional KdV‐Burgers and KdV‐mKdV equations for finding the exact solutions. We use here the fractional complex transform accompanied by properties of local fractional calculus for reduction of fractional partial differential equations to ordinary differential equations. The obtained results are demonstrated by graphs for the new solutions. Copyright © 2016 John Wiley & Sons, Ltd.     

Positivity and stability of linear functional differential equations with infinite delay

General linear functional differential equations with infinite delay are considered. We first give an explicit criterion for positivity of the solution semigroup of linear functional differential equations with infinite delay and then a Perron‐Frobenius type theorem for positive equations. Next, a novel criterion for the exponential asymptotic stability of positive equations is presented. Furthermore, two sufficient conditions for the exponential asymptotic stability of positive equations subjected to structured perturbations and affine perturbations are provided. Finally, we applied the obtained results to problems of the exponential asymptotic stability of Volterra integrodifferential equations. To the best of our knowledge, most of the results of this paper are new.     

Stability of implicit-explicit linear multistep methods for ordinary and delay differential equations

Stability properties of implicit-explicit (IMEX) linear multistep methods for ordinary and delay differential equations are analyzed on the basis of stability regions defined by using scalar test equations. The analysis is closely related to the stability analysis of the standard linear multistep methods for delay differential equations. A new second-order IMEX method which has approximately the same stability region as that of the IMEX Euler method, the simplest IMEX method of order 1, is proposed. Some numerical results are also presented which show superiority of the new method.      

A new analytical technique to solve Volterra's integral equations

The purpose of the present paper is to show that the well‐known homotopy analysis method for solving ordinary and partial differential equations can be applied to solve linear and nonlinear integral equations of Volterra's type with high accuracy as well. Comparison of the present method with Adomian decomposition method (ADM), a well‐known method to solve integral equations, reveals that the ADM is only especial case of the present method. Furthermore, some illustrating examples such as linear, nonlinear and singular integral equations of Volterra's type are given to show high efficiency with reliable accuracy of the method. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical comparison of methods for solving linear differential equations of fractional order

In this article, we implement relatively new analytical techniques, the variational iteration method and the Adomian decomposition method, for solving linear differential equations of fractional order. The two methods in applied mathematics can be used as alternative methods for obtaining analytic and approximate solutions for different types of fractional differential equations. In these schemes, the solution takes the form of a convergent series with easily computable components. This paper will present a numerical comparison between the two methods and a conventional method such as the fractional difference method for solving linear differential equations of fractional order. The numerical results demonstrates that the new methods are quite accurate and readily implemented.