A numerical solution of problems in calculus of variation using direct method and nonclassical parameterization

In this paper a direct method for solving variational problems using nonclassical parameterization is presented. A nonclassical parameterization based on nonclassical orthogonal polynomials is first introduced to reduce a variational problem to a nonlinear mathematical programming problem. Then, using the Lagrange multiplier technique the problem is converted to that of solving a system of algebraic equations. Illustrative examples are given to demonstrate the validity and applicability of the technique.     

Finding the Periodic Solution of Differential Equation via Solving Optimization Problem

In this paper, we propose a new method to find the periodic solutions of differential equations. The key technique is to convert the problem of finding periodic solutions of differential equations into an optimization problem. Then by solving the corresponding optimization problem, we can find the periodic solutions of differential equations. Finally, some numerical results are presented to illustrate the utility of the technique.     

Construction of locally conservative fluxes for the SUPG method

We consider the construction of locally conservative fluxes by means of a simple postprocessing technique obtained from the finite element solutions of advection diffusion equations. It is known that a naive calculation of fluxes from these solutions yields nonconservative fluxes. We consider two finite element methods: the usual continuous Galerkin finite element method for solving nondominating advection diffusion equations and the streamline upwind/Petrov‐Galerkin method for solving advection dominated problems. We then describe the postprocessing technique for constructing conservative fluxes from the numerical solutions of the general variational formulation. The postprocessing technique requires solving an auxiliary Neumann boundary value problem on each element independently and it produces a locally conservative flux on a vertex centered dual mesh relative to the finite element mesh. We provide a convergence analysis for the postprocessing technique. Performance of the technique and the convergence behavior are demonstrated through numerical examples including a set of test problems for advection diffusion equations, advection dominated equations, and drift‐diffusion equations. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1971–1994, 2015     

Numerical treatment for solving fractional Riccati differential equation

This paper presents an accurate numerical method for solving fractional Riccati differential equation (FRDE). The proposed method so called fractional Chebyshev finite difference method (FCheb-FDM). In this technique, we approximate FRDE with a finite dimensional problem. The method is based on the combination of the useful properties of Chebyshev polynomials approximation and finite difference method. The Caputo fractional derivative is replaced by a difference quotient and the integral by a finite sum. By this method the given problem is reduced to a problem for solving a system of algebraic equations, and by solving this system, we obtain the solution of FRDE. Special attention is given to study the convergence analysis and estimate an error upper bound of the obtained approximate formula. Illustrative examples are included to demonstrate the validity and applicability of the proposed technique.     

A numerical solution of singular integrodifferential equations with variable coefficients

A method for the numerical solution of singular integrodifferential equations is presented where the integrals are discretized by using a convenient quadrature rule. Then the problem is reduced to a system of linear algebraic equations by applying the discretized functional equation to appropriately selected collocation points. This technique constitutes an extension of an analogous method convenient for solving singular integral equations which was proposed by the authors.     

A finite volume method based on radial basis functions for two-dimensional nonlinear diffusion equations

The finite volume method is the favoured numerical technique for solving (possibly coupled, nonlinear, anisotropic) diffusion equations. The method transforms the original problem into a system of nonlinear, algebraic equations through the process of discretisation. The accuracy of this discretisation determines to a large extent the accuracy of the final solution.     

Numerical Solution of Arbitrary-Order Ordinary Differential and Integro-Differential Equations with Separated Boundary Conditions Using Optimal Control Technique

In this paper, a new method for solving arbitrary order ordinary differential equations and integro-differential equations of Fredholm and Volterra kind is presented. In the proposed method, these equations with separated boundary conditions are converted to a parametric optimization problem subject to algebraic constraints. Finally, control and state variables will be approximated by a Chebychev series. In this method, a new idea has been used, which offers us the ability of applying the mentioned method for almost all kinds of ordinary differential and integro-differential equations with different types of boundary conditions. The accuracy and efficiency of the proposed numerical technique have been illustrated by solving some test problems.     

An Initial Value Technique for Singularly Perturbed Convection–Diffusion Problems with a Negative Shift

In this paper, a numerical method named as Initial Value Technique (IVT) is suggested to solve the singularly perturbed boundary value problem for the second order ordinary differential equations of convection–diffusion type with a delay (negative shift). In this technique, the original problem of solving the second order equation is reduced to solving two first order differential equations, one of which is singularly perturbed without delay and other one is regular with a delay term. The singularly perturbed problem is solved by the second order hybrid finite difference scheme, whereas the delay problem is solved by the fourth order Runge–Kutta method with Hermite interpolation. An error estimate is derived by using the supremum norm. Numerical results are provided to illustrate the theoretical results.     

New Approach to Solving a Two-Dimensional Problem for an Orthotropic Plate

A new approach to solving a two-dimensional problem for an orthotropic multiply connected plate is proposed in place of the standard approach which reduces the problem to integrating a fourth-order differential equation. The new method reduces the problem to integrating second-order differential equations that can be solved successively using a perturbation technique.     

General fractional variational problem depending on indefinite integrals

In this report, we consider two kind of general fractional variational problem depending on indefinite integrals include unconstrained problem and isoperimetric problem. These problems can have multiple dependent variables, multiorder fractional derivatives, multiorder integral derivatives and boundary conditions. For both problems, we obtain the Euler-Lagrange type necessary conditions which must be satisfied for the given functional to be extremum. Also, we apply the Rayleigh-Ritz method for solving the unconstrained general fractional variational problem depending on indefinite integrals. By this method, the given problem is reduced to the problem for solving a system of algebraic equations using shifted Legendre polynomials basis functions. An approximate solution for this problem is obtained by solving the system. We discuss the analytic convergence of this method and finally by some examples will be showing the accurately and applicability for this technique.     

Distributed optimal control of the viscous Burgers equation via a Legendre pseudo‐spectral approach

This paper presents a computational technique based on the pseudo‐spectral method for the solution of distributed optimal control problem for the viscous Burgers equation. By using pseudo‐spectral method, the problem is converted to a classical optimal control problem governed by a system of ordinary differential equations, which can be solved by well‐developed direct or indirect methods. For solving the resulting optimal control problem, we present an indirect method by deriving and numerically solving the first‐order optimality conditions. Numerical tests involving both unconstrained and constrained control problems are considered. Copyright © 2015 John Wiley & Sons, Ltd.     

Vieta–Fibonacci wavelets: Application in solving fractional pantograph equations

In this paper, the Vieta–Fibonacci wavelets as a new family of orthonormal wavelets are generated. An operational matrix concerning fractional integration of these wavelets is extracted. A numerical scheme is established based on these wavelets and their fractional integral matrix together with the collocation technique to solve fractional pantograph equations. The presented method reduces solving the problem under study into solving a system of algebraic equations. Several examples are provided to show the accuracy of the method.     

Factorization method for linear and quasilinear singularly perturbed boundary value problems for ordinary differential equations

For linear singularly perturbed boundary value problems, we come up with a method that reduces solving a differential problem to a discrete (difference) problem. Difference equations, which are an exact analog of differential equations, are constructed by the factorization method. Coefficients of difference equations are calculated by solving Cauchy problems for first-order differential equations. In this case nonlinear Ricatti equations with a small parameter are solved by asymptotic methods, and solving linear equations reduces to computing quadratures. A solution for quasilinear singularly perturbed equations is obtained by means of an implicit relaxation method. A solution to a linearized problem is calculated by analogy with a linear problem at each iterative step. The method is tested against solutions to the known Lagerstrom-Cole problem.     

Convergence of the forward-backward sweep method in optimal control

The Forward-Backward Sweep Method is a numerical technique for solving optimal control problems. The technique is one of the indirect methods in which the differential equations from the Maximum Principle are numerically solved. After the method is briefly reviewed, two convergence theorems are proved for a basic type of optimal control problem. The first shows that recursively solving the system of differential equations will produce a sequence of iterates converging to the solution of the system. The second theorem shows that a discretized implementation of the continuous system also converges as the iteration and number of subintervals increases. The hypotheses of the theorem are a combination of basic Lipschitz conditions and the length of the interval of integration. An example illustrates the performance of the method.     

利用随机模拟方法求解第二类积分方程

给出一种求解第二类Fredholm和Volterra积分方程的数值算法,算法在数值积分技术的基础上使用Monte Carlo随机模拟方法求积分方程的近似解.通过数值例子证明了该算法是有效的.     

A Generalization of Ritz-Variational Method for Solving a Class of Fractional Optimization Problems

This paper presents an approximate method for solving a class of fractional optimization problems with multiple dependent variables with multi-order fractional derivatives and a group of boundary conditions. The fractional derivatives are in the Caputo sense. In the presented method, first, the given optimization problem is transformed into an equivalent variational equality; then, by applying a special form of polynomial basis functions and approximations, the variational equality is reduced to a simple linear system of algebraic equations. It is demonstrated that the derived linear system has a unique solution. We get an approximate solution for the initial optimization problem by solving the final linear system of equations. The choice of polynomial basis functions provides a method with such flexibility that all initial and boundary conditions of the problem can be easily imposed. We extensively discuss the convergence of the method and, finally, present illustrative test examples to demonstrate the validity and applicability of the new technique.     

A numerical technique for solving a class of fractional variational problems

This paper presents a numerical method for solving a class of fractional variational problems (FVPs) with multiple dependent variables, multi order fractional derivatives and a group of boundary conditions. The fractional derivative in the problem is in the Caputo sense. In the presented method, the given optimization problem reduces to a system of algebraic equations using polynomial basis functions. An approximate solution for the FVP is achieved by solving the system. The choice of polynomial basis functions provides the method with such a flexibility that initial and boundary conditions can be easily imposed. We extensively discuss the convergence of the method and finally present illustrative examples to demonstrate validity and applicability of the new technique.     

求解非线性方程组的连续极小化方法

1.引言 求非线性方程组 F(x)=0 (1)(F:D?R~n→R~n)的各种方法中,牛顿法最为基本.但它只有局部收敛性和半局部收敛性,而且要求DF(x)~(-1)存在.为了扩大收敛范围及克服DF(x)奇异性带来的困难,用“连续化”的思想求方程(1)的解是一个有效的途径.这方面,已有许多工作,如[3—6].本文利用常微分方程几何理论,对连续化方法进行 些探讨,给出了沿积分曲线极小化求非线性方程组(1)的解的方法.考虑如下给定函数:     

Uzawa methods for a class of block three‐by‐three saddle‐point problems

In this work, we consider numerical methods for solving a class of block three‐by‐three saddle‐point problems, which arise from finite element methods for solving time‐dependent Maxwell equations and some other applications. The direct extension of the Uzawa method for solving this block three‐by‐three saddle‐point problem requires the exact solution of a symmetric indefinite system of linear equations at each step. To avoid heavy computations at each step, we propose an inexact Uzawa method, which solves the symmetric indefinite linear system in some inexact way. Under suitable assumptions, we show that the inexact Uzawa method converges to the unique solution of the saddle‐point problem within the approximation level. Two special algorithms are customized for the inexact Uzawa method combining the splitting iteration method and a preconditioning technique, respectively. Numerical experiments are presented, which demonstrated the usefulness of the inexact Uzawa method and the two customized algorithms.     

A numerical technique for solving fractional variational problems

This paper presents an accurate numerical method for solving a class of fractional variational problems (FVPs). The fractional derivative in these problems is in the Caputo sense. The proposed method is called fractional Chebyshev finite difference method. In this technique, we approximate FVPs and end up with a finite‐dimensional problem. The method is based on the combination of the useful properties of Chebyshev polynomials approximation and finite difference method. The Caputo fractional derivative is replaced by a difference quotient and the integral by a finite sum. The fractional derivative approximation using Clenshaw and Curtis formula introduced here, along with Clenshaw and Curtis procedure for the numerical integration of a non‐singular functions and the Rayleigh–Ritz method for the constrained extremum, is considered. By this method, the given problem is reduced to the problem for solving a system of algebraic equations, and by solving this system, we obtain the solution of FVPs. Special attention is given to study the convergence analysis and evaluate an error upper bound of the obtained approximate formula. Illustrative examples are included to demonstrate the validity and applicability of the proposed technique. A comparison with another method is given. Copyright © 2012 John Wiley & Sons, Ltd.