一类求解非线性代数方程组的并行多分裂AOR算法

本文提出了求解大型非线性代数方程组Aф(x)+Bψ(x)=b的并行多分裂AOR(Accelerated Overrela Xation)算法。在一定的条件下,证明了非线性代数方程组解的存在唯一性,并建立了新算法的全局收敛性理论。     

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.     

Numerical solutions of systems of nonlinear integro-differential equations by Homotopy-perturbation method

In this study, the numerical solutions of a system of two nonlinear integro-differential equations, which describes biological species living together, are derived employing the well-known Homotopy-perturbation method. The approximate solutions are in excellent agreement with those obtained by the Adomian decomposition method. Furthermore, we present an analytical approximate solution for a more general form of the system of nonlinear integro-differential equations. The numerical result indicates that the proposed method is straightforward to implement, efficient and accurate for solving nonlinear integro-differential equations.     

A Positivity-preserving Finite Element Method for Quantum Drift-diffusion Model

In this paper, we propose a positivity-preserving finite element method for solving the three-dimensional quantum drift-diffusion model. The model consists of five nonlinear elliptic equations, and two of them describe quantum corrections for quasi-Fermi levels. We propose an interpolated-exponential finite element (IEFE) method for solving the two quantum-correction equations. The IEFE method always yields positive carrier densities and preserves the positivity of second-order differential operators in the Newton linearization of quantum-correction equations. Moreover, we solve the two continuity equations with the edge-averaged finite element (EAFE) method to reduce numerical oscillations of quasi-Fermi levels. The Poisson equation of electrical potential is solved with standard Lagrangian finite elements. We prove the existence of solution to the nonlinear discrete problem by using a fixed-point iteration and solving the minimum problem of a new discrete functional. A Newton method is proposed to solve the nonlinear discrete problem. Numerical experiments for a three-dimensional nano-scale FinFET device show that the Newton method is robust for source-to-gate bias voltages up to 9V and source-to-drain bias voltages up to 10V.     

A spectral solution of nonlinear mean field dynamo equations: Without inertia

In this paper we describe in detail spectral solution method for solving nonlinear mean field dynamo equations without inertial effects in a rapidly rotating spherical shell filled with electrically conducting fluid taken into account no-slip velocity boundary conditions in the core for a finitely conducting inner core and an insulating mantle. Sample results are provided for nonlinear alpha squared dynamos. The method is suitable for solving many geophysical, astrophysical and engineering problems.     

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.     

Constructive proof for existence of nonlinear two-point boundary value problems

Owing to the importance of differential equations in physics, the existence of solutions for differential equations has been paid much attention. In this paper, the existence of solution are obtained for the nonlinear second order two-point boundary value problem in the reproducing kernel space. Under certain assumptions on right-hand side, we propose constructive proof for the existence result, and a method is presented to obtain the exact solution expressed by the form of series. This paper is a extension of previous paper [Wei Jiang, Minggen Cui, The exact solution and stability analysis for integral equation of third or first kind with singular kernel, Appl. Math. Comput. 202 (2) (2008) 666-674], which extends a method of solving linear problems to present method for solving nonlinear problems.     

Analytical approximations for a population growth model with fractional order

In this paper, we apply the homotopy analysis method (HAM) to solve the fractional Volterra’s model for population growth of a species in a closed system. This technique is extended to give solutions for nonlinear fractional integro–differential equations. The whole HAM solution procedure for nonlinear fractional differential equations is established. Further, the accurate analytical approximations are obtained for the first time, which are valid and convergent for all time t. This indicates the validity and great potential of the homotopy analysis method for solving nonlinear fractional integro–differential equations.     

Efficient analytic method for solving nonlinear fractional differential equations

In this paper, we propose a new approach of the generalized differential transform method (GDTM) for solving nonlinear fractional differential equations. In GDTM, it is a key to derive a recurrence relation of generalized differential transform (GDT) associated with the solution in the given fractional equation. However, the recurrence relations of complex nonlinear functions such as exponential, logarithmic and trigonometry functions have not been derived before in GDTM. We propose new algorithms to construct the recurrence relations of complex nonlinear functions and apply the GDTM with the proposed algorithms to solve nonlinear fractional differential equations. Several illustrative examples are demonstrated to show the effectiveness of the proposed method. It is shown that the proposed technique is robust and accurate for solving fractional differential equations.     

Numerical solution of Fredholm integral equations of the second kind by using integral mean value theorem II. High dimensional problems

In this work, we generalize the numerical method discussed in [Z. Avazzadeh, M. Heydari, G.B. Loghmani, Numerical solution of Fredholm integral equations of the second kind by using integral mean value theorem, Appl. math. modelling, 35 (2011) 2374–2383] for solving linear and nonlinear Fredholm integral and integro-differential equations of the second kind. The presented method can be used for solving integral equations in high dimensions. In this work, we describe the integral mean value method (IMVM) as the technical algorithm for solving high dimensional integral equations. The main idea in this method is applying the integral mean value theorem. However the mean value theorem is valid for multiple integrals, we apply one dimensional integral mean value theorem directly to fulfill required linearly independent equations. We solve some examples to investigate the applicability and simplicity of the method. The numerical results confirm that the method is efficient and simple.     

The solution of high-order nonlinear ordinary differential equations by Chebyshev Series

By the use of the Chebyshev series, a direct computational method for solving the higher order nonlinear differential equations has been developed in this paper. This method transforms the nonlinear differential equation into the matrix equation, which corresponds to a system of nonlinear algebraic equations with unknown Chebyshev coefficients, via Chebyshev collocation points. The solution of this system yields the Chebyshev coefficients of the solution function. An algorithm for this nonlinear system is also proposed in this paper. The method is valid for both initial-value and boundary-value problems. Several examples are presented to illustrate the accuracy and effectiveness of the method.     

Finite-difference methods for solving the reaction-diffusion equations of a simple isothermal chemical system

Numerical methods are proposed for the numerical solution of a system of reaction-diffusion equations, which model chemical wave propagation. The reaction terms in this system of partial differential equations contain nonlinear expressions. Nevertheless, it is seen that the numerical solution is obtained by solving a linear algebraic system at each time step, as opposed to solving a nonlinear algebraic system, which is often required when integrating nonlinear partial differential equations. The development of each numerical method is made in the light of experience gained in solving the system of ordinary differential equations, which model the well-stirred analogue of the chemical system. The first-order numerical methods proposed for the solution of this initialvalue problem are characterized to be implicit. However, in each case it is seen that the numerical solution is obtained explicitly. In a series of numerical experiments, in which the ordinary differential equations are solved first of all, it is seen that the proposed methods have superior stability properties to those of the well-known, first-order, Euler method to which they are compared. Incorporating the proposed methods into the numerical solution of the partial differential equations is seen to lead to two economical and reliable methods, one sequential and one parallel, for solving the travelling-wave problem. © 1994 John Wiley & Sons, Inc.     

Stein implicit Runge-Kutta methods with high stage order for large-scale ordinary differential equations

The Runge-Kutta method is one of the most popular implicit methods for the solution of stiff ordinary differential equations. For large problems, the main drawback of such methods is the cost required at each integration step for computing the solution of a nonlinear system of equations. In this paper, we propose to reduce the cost of the computation by transforming the linear systems arising in the application of Newton's method to Stein matrix equations. We propose an iterative projection method onto block Krylov subspaces for solving numerically such Stein matrix equations. Numerical examples are given to illustrate the performance of our proposed method.     

A Reproducing Kernel Hilbert Space Method for Solving Integro-Differential Equations of Fractional Order

In this article, we implement a relatively new analytical technique, the reproducing kernel Hilbert space method (RKHSM), for solving integro-differential equations of fractional order. The solution obtained by using the method takes the form of a convergent series with easily computable components. Two numerical examples are studied to demonstrate the accuracy of the present method. The present work shows the validity and great potential of the reproducing kernel Hilbert space method for solving linear and nonlinear integro-differential equations of fractional order.     

Inexact implicit methods for monotone general variational inequalities

Solving a variational inequality problem is equivalent to finding a solution of a system of nonsmooth equations. Recently, we proposed an implicit method, which solves monotone variational inequality problem via solving a series of systems of nonlinear smooth (whenever the operator is smooth) equations. It can exploit the facilities of the classical Newton–like methods for smooth equations. In this paper, we extend the method to solve a class of general variational inequality problems Moreover, we improve the implicit method to allow inexact solutions of the systems of nonlinear equations at each iteration. The method is shown to preserve the same convergence properties as the original implicit method. Received July 31, 1995 / Revised version received January 15, 1999? Published online May 28, 1999     

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.     

在再生核空间中方程AuBu+Cu＝f的求解方法

$ 1 引言 本文研究下面一类非线性算子方程求解问题 AμBμ Cμ=f, (1.1)其中f,μ∈W(Ω),μ(O)=1,||f ||=1,A,B,C∈(W(Ω)→W(Ω)),(W(Ω)→W(Ω))是W(Ω)到W(Ω)的连续线性算子空间,W(Ω)是定义在Ω域上的(Ω是实数域R的有界域)再生核空间。 本文是在再生核空间上,通过将一维非线性算子方程(1.1)转化为二维线性算子方     

Chebyshev operational matrix method for solving multi-order fractional ordinary differential equations

The aim of this article is to present an analytical approximation solution for linear and nonlinear multi-order fractional differential equations (FDEs) by extending the application of the shifted Chebyshev operational matrix. For this purpose, we convert FDE into a counterpart system and then using proposed method to solve the resultant system. Our results in solving four different linear and nonlinear FDE, confirm the accuracy of proposed method.     

Approximations of the nonlinear Volterra's population model by an efficient numerical method

In this paper, we apply the new homotopy perturbation method to solve the Volterra's model for population growth of a species in a closed system. This technique is extended to give solution for nonlinear integro‐differential equation in which the integral term represents the total metabolism accumulated fromtime zero. The approximate analytical procedure only depends on two components. The newhomotopy perturbationmethodwas applied to nonlinear integro‐differential equations directly and by converting the problem into nonlinear ordinary differential equation. We also compare this method with some other numerical results and show that the present approach is less computational and is applicable for solving nonlinear integro‐differential equations and ordinary differential equations as well. Copyright © 2011 John Wiley & Sons, Ltd.     

FINITE ELEMENT NONLINEAR GALERKIN COUPLING METHOD FOR THE EXTERIOR STEADY NAVIER-STOKES PROBLEM

1.IntroductionNonlinearGalerkinmethodsaremultilevelschemesforthedissipativeevolutionpartialdifferentialequations.Theycorrespondtothesplittingsoftheunknownu:u=y z)wherethecomponentsareofdifferentorderofmagnitudewithrespecttoaparameterrelatedtothespati...