A stochastic Lagrangian representation of the three‐dimensional incompressible Navier‐Stokes equations

In this paper we derive a probabilistic representation of the deterministic three‐dimensional Navier‐Stokes equations based on stochastic Lagrangian paths. The particle trajectories obey SDEs driven by a uniform Wiener process; the inviscid Weber formula for the Euler equations of ideal fluids is used to recover the velocity field. This method admits a self‐contained proof of local existence for the nonlinear stochastic system and can be extended to formulate stochastic representations of related hydrodynamic‐type equations, including viscous Burgers equations and Lagrangian‐averaged Navier‐Stokes alpha models. © 2007 Wiley Periodicals, Inc.     

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.     

Nonautonomous first integrals and general solutions of the KdV‐Burgers and mKdV‐Burgers equations with the source

The method for constructing first integrals and general solutions of nonlinear ordinary differential equations is presented. The method is based on index accounting of the Fuchs indices, which appeared during the Painlevé test of a nonlinear differential equation. The Fuchs indices indicate us the leading members of the first integrals for the origin differential equation. Taking into account the values of the Fuchs indices, we can construct the auxiliary equation, which allows to look for the first integrals of nonlinear differential equations. The method is used to obtain the first integrals and general solutions of the KdV‐Burgers and the mKdV‐Burgers equations with a source. The nonautonomous first integrals in the polynomials form are found. The general solutions of these nonlinear differential equations under at some additional conditions on the parameters of differential equations are also obtained. Illustrations of some solutions of the KdV‐Burgers and the mKdV‐Burgers are given.     

Interpolants for sixth‐order Numerov‐type methods

The classical four‐stage family of explicit sixth‐order Numerov‐type method is considered. We provide two kinds of interpolants: (a) a three‐step interpolation based on all available data at mesh points and (b) a local interpolant (ie, two steps) that is constructed after solving scaled equations of condition. These latter equations are explained and provided here. Applying these interpolants in a set of tests, we conclude that they produce global errors of the same magnitude with the underlying method.     

Mixed two‐grid finite difference methods for solving one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations

The aim of this paper is to propose mixed two‐grid finite difference methods to obtain the numerical solution of the one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations. The finite difference equations at all interior grid points form a large‐sparse linear system, which needs to be solved efficiently. The solution cost of this sparse linear system usually dominates the total cost of solving the discretized partial differential equation. The proposed method is based on applying a family of finite difference methods for discretizing the spatial and time derivatives. The obtained system has been solved by two‐grid method, where the two‐grid method is used for solving the large‐sparse linear systems. Also, in the proposed method, the spectral radius with local Fourier analysis is calculated for different values of h and Δt. The numerical examples show the efficiency of this algorithm for solving the one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations. Copyright © 2016 John Wiley & Sons, Ltd.     

Strong convergence of the split‐step θ‐method for stochastic age‐dependent population equations

In this paper, we constructed the split‐step θ (SSθ)‐method for stochastic age‐dependent population equations. The main aim of this paper is to investigate the convergence of the SS θ‐method for stochastic age‐dependent population equations. It is proved that the proposed method is convergent with strong order 1/2 under given conditions. Finally, an example is simulated to verify the results obtained from the theory, and comparative analysis with Euler method is given, the results show the higher accuracy of the SS θ‐method. Copyright © 2014 John Wiley & Sons, Ltd.     

Bianchi Permutability for the Anti‐Self‐Dual Yang‐Mills Equations

The anti‐self‐dual Yang‐Mills equations are known to have reductions to many integrable differential equations. A general Bäcklund transformation (BT) for the anti‐self‐dual Yang‐Mills (ASDYM) equations generated by a Darboux matrix with an affine dependence on the spectral parameter is obtained, together with its Bianchi permutability equation. We give examples in which we obtain BTs of symmetry reductions of the ASDYM equations by reducing this ASDYM BT. Some discrete integrable systems are obtained directly from reductions of the ASDYM Bianchi system.     

WEB‐spline based mesh‐free finite element analysis for the heat equation and the time‐dependent Navier–Stokes equation: A survey

In this article, we give some numerical techniques and error estimates using web‐spline based mesh‐free finite element method for the heat equation and the time‐dependent Navier–Stokes equations on bounded domains. The web‐spline method uses weighted extended B‐splines on a regular grid as basis functions and does not require any grid generation. We demonstrate the method by providing numerical results for the Poisson's and stationary Stokes equation. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Numerical solutions of the space‐time fractional advection‐dispersion equation

Fractional advection‐dispersion equations are used in groundwater hydrologhy to model the transport of passive tracers carried by fluid flow in a porous medium. In this paper we present two reliable algorithms, the Adomian decomposition method and variational iteration method, to construct numerical solutions of the space‐time fractional advection‐dispersion equation in the form of a rabidly convergent series with easily computable components. The fractional derivatives are described in the Caputo sense. Some examples are given. Numerical results show that the two approaches are easy to implement and accurate when applied to space‐time fractional advection‐dispersion equations. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

Second‐order Galerkin‐Lagrange method for the Navier‐Stokes equations (retracted article)

It has come to the attention of the editors and publisher that an article published in Numerical Methods and Partial Differential Equations, “Second‐order Galerkin‐Lagrange method for the Navier‐Stokes equations,” by Mohamed Bensaada, Driss Esselaoui, and Pierre Saramito, Numer Methods Partial Differential Eq 21(6) (2005), 1099–1121 included large portions that were copied from the following paper without proper citation: “Convergence and nonlinear stability of the Lagrange‐Galerkin method for the Navier‐Stokes equations,” Endre Suli, Numerische Mathematik, Vol. 53, No. 4, pp. 459–486 (July, 1988). We have retracted the paper and apologize to Dr. Suli Numer Methods Partial Differential Eq (2007)23(1)211 .     

Application of the modified operational matrices in multiterm variable‐order time‐fractional partial differential equations

In this paper, we present a novel discrete scheme based on Genocchi polynomials and fractional Laguerre functions to solve multiterm variable‐order time‐fractional partial differential equations (M‐V‐TFPDEs) in the large interval. In this purpose, the accurate modified operational matrices are constructed to reduce the problems into a system of algebraic equations. Also, the computational algorithm based on the method and modified operational matrices in the large interval is easily implemented. Furthermore, we discuss the error estimation of the proposed method. Ultimately, to confirm our theoretical analysis and accuracy of numerical approach, several examples are presented.     

Boundary‐value problems for the third‐order loaded equation with noncharacteristic type‐change boundaries

The purpose of this paper is to establish unique solvability for a certain generalized boundary‐value problem for a loaded third‐order integro‐differential equation with variable coefficients. Moreover, the method of integral equations is applied to obtain an equation related to the Riemann‐Liouville operators.     

A uniform estimate for the MMOC for two‐dimensional advection‐diffusion equations

We prove an optimal‐order error estimate in a weighted energy norm for the modified method of characteristics (MMOC) and the modified method of characteristics with adjusted advection (MMOCAA) for two‐dimensional time‐dependent advection‐diffusion equations, in the sense that the generic constants in the estimates depend on certain Sobolev norms of the true solution but not on the scaling diffusion parameter ε. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

A preconditioned low‐rank CG method for parameter‐dependent Lyapunov matrix equations

This paper is concerned with the numerical solution of symmetric large‐scale Lyapunov equations with low‐rank right‐hand sides and coefficient matrices depending on a parameter. Specifically, we consider the situation when the parameter dependence is sufficiently smooth, and the aim is to compute solutions for many different parameter samples. On the basis of existing results for Lyapunov equations and parameter‐dependent linear systems, we prove that the tensor containing all solution samples typically allows for an excellent low multilinear rank approximation. Stacking all sampled equations into one huge linear system, this fact can be exploited by combining the preconditioned CG method with low‐rank truncation. Our approach is flexible enough to allow for a variety of preconditioners based, for example, on the sign function iteration or the alternating direction implicit method. Copyright © 2013 John Wiley & Sons, Ltd.     

Negative norm least‐squares methods for the velocity‐vorticity‐pressure Navier–Stokes equations

We develop and analyze a least‐squares finite element method for the steady state, incompressible Navier–Stokes equations, written as a first‐order system involving vorticity as new dependent variable. In contrast to standard L² least‐squares methods for this system, our approach utilizes discrete negative norms in the least‐squares functional. This allows us to devise efficient preconditioners for the discrete equations, and to establish optimal error estimates under relaxed regularity assumptions. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 237–256, 1999     

Numerical solutions of nonlinear 3‐dimensional Volterra integral‐differential equations with 3D‐block‐pulse functions

This paper studies nonlinear 3‐dimensional Volterra integral‐differential equations, by implementing 3‐dimensional block‐pulse functions. First, we prove a theorem and corollary about sufficient condition for the minimum of mean square error under the block pulse coefficients and uniqueness of solution of the nonlinear Volterra integral‐differential equations. Then, we convert the main problem to a nonlinear system to the 3‐dimensional block‐pulse functions. In addition, illustrative examples are included to demonstrate the validity and applicability of the presented method.     

Finite element method for two‐dimensional time‐fractional tricomi‐type equations

In this article, we consider the finite element method (FEM) for two‐dimensional linear time‐fractional Tricomi‐type equations, which is obtained from the standard two‐dimensional linear Tricomi‐type equation by replacing the first‐order time derivative with a fractional derivative (of order α, with 1 <α< 2 ). The method is based on finite element method for space and finite difference method for time. We prove that the method is unconditionally stable, and the error estimate is presented. The comparison of the FEM results with the exact solutions is made, and numerical experiments reveal that the FEM is very effective. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Analysis of a particle method for the one‐dimensional Vlasov‐Maxwell system

In this article we present a particle method for solving numerically the one‐dimensional Vlasov‐Maxwell equations. This method is based on the formulation by characteristics. We perform the error analysis and we investigate the properties of this scheme. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A two‐level finite element method for the stationary Navier‐Stokes equations based on a stabilized local projection

In this article, we propose a two‐level finite element method to analyze the approximate solutions of the stationary Navier‐Stokes equations based on a stabilized local projection. The local projection allows to circumvent the Babuska‐Brezzi condition by using equal‐order finite element pairs. The local projection can be used to stabilize high equal‐order finite element pairs. The proposed method combines the local projection stabilization method and the two‐level method under the assumption of the uniqueness condition. The two‐level method consists of solving a nonlinear equation on the coarse mesh and solving a linear equation on fine mesh. The nonlinear equation is solved by the one‐step Newtonian iteration method. In the rest of this article, we show the error analysis of the lowest equal‐order finite element pair and provide convergence rate of approximate solutions. Furthermore, the numerical illustrations coincide with the theoretical analysis expectations. From the view of computational time, the results show that the two‐level method is effective to solve the stationary Navier‐Stokes equations. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011