Numerical solutions of generalized Burgers–Fisher and generalized Burgers–Huxley equations using collocation of cubic B-splines

In this work, we propose a numerical scheme to obtain approximate solutions of generalized Burgers–Fisher and Burgers–Huxley equations. The scheme is based on collocation of modified cubic B-spline functions and is applicable for a class of similar diffusion–convection–reaction equations. We use modified cubic B-spline functions for space variable and for its derivatives to obtain a system of first-order ordinary differential equations in time. We solve this system by using SSP-RK54 scheme. The stability of the method has been discussed and it is shown that the method is unconditionally stable. The approximate solutions have been computed without using any transformation or linearization. The proposed scheme needs less storage space and execution time. The test problems considered by the different researchers have been discussed to demonstrate the strength and utility of the proposed scheme. The computed numerical solutions are in good agreement with the exact solutions and competent with those available in the literature. The scheme is simple as well as computationally efficient. The scheme provides approximate solution not only at the grid points but also at any point in the solution range.     

A compact split-step finite difference method for solving the nonlinear Schrödinger equations with constant and variable coefficients

We propose a compact split-step finite difference method to solve the nonlinear Schrödinger equations with constant and variable coefficients. This method improves the accuracy of split-step finite difference method by introducing a compact scheme for discretization of space variable while this improvement does not reduce the stability range and does not increase the computational cost. This method also preserves some conservation laws. Numerical tests are presented to confirm the theoretical results for the new numerical method by using the cubic nonlinear Schrödinger equation with constant and variable coefficients and Gross-Pitaevskii equation.     

Fast and efficient numerical method for solving the Allen–Cahn equation on the cubic surface

In this study, we present a fast and efficient finite difference method (FDM) for solving the Allen–Cahn (AC) equation on the cubic surface. The proposed method applies appropriate boundary conditions in the two-dimensional (2D) space to calculate numerical solutions on cubic surfaces, which is relatively simpler than a direct computation in the three-dimensional (3D) space. To numerically solve the AC equation on the cubic surface, we first unfold the cubic surface domain in the 3D space into the 2D space, and then apply the FDM on the six planar sub-domains with appropriate boundary conditions. The proposed method solves the AC equation using an operator splitting method that splits the AC equation into the linear and nonlinear terms. To demonstrate that the proposed algorithm satisfies the properties of the AC equation on the cubic surface, we perform the numerical experiments such as convergence test, total energy decrease, and maximum principle.     

A conservative exponential time differencing method for the nonlinear cubic Schrödinger equation

A mass and energy conservative exponential time differencing scheme using the method of lines is proposed for the numerical solution of a certain family of first-order time-dependent PDEs. The resulting nonlinear system is solved with an unconditionally stable modified predictor–corrector method using a second-order explicit scheme. The efficiency of the method introduced is analyzed and discussed by applying it to the nonlinear cubic Schrödinger equation. The results arising from the experiments for the single, the double soliton waves and the system of two Schrödinger equations are compared with relevant known ones.     

Solving Index Form Equations in Fields of Degree 9 with Cubic Subfields

We describe an efficient algorithm for solving index form equations in number fields of degree 9 which are composites of cubic fields with coprime discriminants. We develop the algorithm in detail for the case of complex cubic fields, but the main steps of the procedure are also applicable for other cases. Our most important tool is the main theorem of a recent paper of Gaál (1998a). In view of this result the index form equation in the ninth degree field implies relative index form equations over the subfields. In our case these equations are cubic relative Thue equations over cubic fields. The main purpose of the paper is to show that this approach is much more efficient than the direct method, which consists of reducing the index form equation to unit equations over the normal closure of the original field. At the end of the paper we describe our computational experience. Many ideas of the paper can be applied to develop fast algorithms for solving index form equations in other types of higher degree fields which are composites of subfields.     

高精度三次参数样条曲线的构造

构造参数样条曲线的关键是选取节点，该文讨论了GC^2三次参数样条曲线需满足的连续性方程，提出了构造GC^2三次参数样条曲线的新方法，在讨论了平面有序五点确定一组三次多项式函数曲线，平面有序六点唯一确定一条三次多项式函数曲线的基础上，提出了计算相邻两区间上的节点的算法，构造的插值曲线具有三次多项式函数精,该文还以实例对新方法与其它方法构造的插值曲线的精度进行了比较。     

用MATLAB进行动力学实验数据处理

反应动力学实验数据误差较大且分布不均匀,一般的数据处理方法往往不适用。本文充分利用所得数据和化学反应方程组自身隐含的数学性质,首先采用自然三次样条进行数据拟合,然后利用拟合出的函数关系作为迭代计算反应速率常数的初值,给出了一种求解反应速率常数的方法,并采用数值软件MATLAB编程。给出了一个实例,计算表明本文方法对于处理误差较大且分布不均匀的动力学实验数据,取得了令人满意的效果。     

New numerical method for the solutions of the MCDF equations based on the CIP method

A new numerical method based on the constrained interpolation profile (CIP) method to solve the Multiconfiguration Dirac-Fock (MCDF) equations is presented. The radial wave functions are represented by the values and the spatial derivatives on an arbitrary grid system, and approximated by cubic polynomials. Owing to this representation, the values and the spatial derivatives of the effective charge distribution and inhomogeneous term are calculated using the previous cycle's wave functions. Then the homogeneous MCDF equations are integrated to obtain two linearly independent solutions, which are used to construct the Green function, by the adaptive stepsize controlled Runge-Kutta method controlling the truncation errors within a prescribed accuracy. The radial wave function is improved by taking the convolution of the Green function and the inhomogeneous term. The effectiveness of this numerical procedure is investigated after implementing it into the relativistic atomic structure code GRASP92.     

二维非结构网格Hamilton-Jacobi方程的一种简化的加权ENO格式

考虑标量Hamilton-Jacobi方程,对二维非结构网格给出了一种简化的三阶精度加权ENO格式.方法的主要思想是时间和空间分开处理,时间离散用三阶TVD Runge-Kutta 方法.对空间,在每一个三角形单元上构造一个三次多项式,该多项式是一些三次多项式的加权,并给出了加权因子的构造方法.最后用该格式对一些典型算例进行了数值试验,并分析了方法的精度,结果表明该格式是成功的.     

Non-polynomial cubic spline methods for the solution of parabolic equations

Second-order parabolic partial differential equations are solved by using a new three level method based on non-polynomial cubic spline in the space direction and finite difference in the time direction. Stability analysis of the method has been carried out and we have shown that our method is unconditionally stable. It has been shown that by suitably choosing the parameters most of the previous known methods for homogeneous and non-homogeneous cases can be obtained from our method. We also obtain a new high accuracy scheme of O(k ⁴+h ⁴). Numerical examples are given to illustrate the applicability and efficiency of the new method.     

求解对流扩散方程的混合三次B-样条配点法

通过对三次B-样条和三次三角B-样条基函数引入权因子[ω],给出了对流扩散方程的混合三次B-样条配点法。对对流扩散方程空间离散采用混合三次B-样条配点法和时间离散采用向前有限差分,引入参数[θ],建立差分格式。对差分格式的稳定性进行分析,得到稳定性条件。数值实验表明所构造方法的有效性,并且适当调整权因子[ω]和参数[θ]的值,可提高计算的精度。     

Numerical solution of nonlinear system of second-order boundary value problems using cubic B-spline scaling functions

A numerical technique is presented for the solution of nonlinear system of second-order boundary value problems. This method uses the cubic B-spline scaling functions. The method consists of expanding the required approximate solution as the elements of cubic B-spline scaling function. Using the operational matrix of derivative, we reduce the problem to a set of algebraic equations. Numerical examples are included to demonstrate the validity and applicability of the technique. The method is easy to implement and produces very accurate results.     

A new cubic convergent method for solving a system of nonlinear equations

We present a new cubic convergent method for solving a system of nonlinear equations. The new method can be viewed as a modified Chebyshev's method in which the difference of Jacobian matrixes replaces three order tensor. Therefore, the new method reduces the storage and computational cost. The new method possesses the local cubic convergence as well as Chebyshev's method. A rule is deduced to ensure the descent property of the search direction, and a nonmonotone line search technique is used to guarantee the global convergence. Numerical results indicate that the new method is competitive and efficient for some classical test problems.     

Numerical solution of Riccati equation using the cubic B-spline scaling functions and Chebyshev cardinal functions

Two numerical techniques are presented for solving the solution of Riccati differential equation. These methods use the cubic B-spline scaling functions and Chebyshev cardinal functions. The methods consist of expanding the required approximate solution as the elements of cubic B-spline scaling function or Chebyshev cardinal functions. Using the operational matrix of derivative, we reduce the problem to a set of algebraic equations. Some numerical examples are included to demonstrate the validity and applicability of the new techniques. The methods are easy to implement and produce very accurate results.     

The spline alternating group explicit (splage) method for two dimensional problems

A new group explicit iterative method based on cubic spline approximations is presented for the numerical solution of partial differential equations. The numerical results obtained confirm the viability of the method.     

Numerical solution of nonlinear system of Klein–Gordon equations by cubic B-spline collocation method

A technique to approximate the solutions of nonlinear Klein–Gordon equation and Klein–Gordon-Schrödinger equations is presented separately. The approach is based on collocation of cubic B-spline functions. The above-mentioned equations are decomposed into a system of partial differential equations, which are further converted to an amenable system of ODEs. The obtained system has been solved by SSP-RK54 scheme. Numerical solutions are presented for five examples, to show the accuracy and usefulness of proposed approach. The approximate solutions of both the equations are computed without using any transformation and linearization. The technique can be applied with ease to solve linear and nonlinear PDEs and also reduces the computational work.     

Multiple shooting and the calculation of some types of singularities in b.v.p.'s

We consider simple and multiple shooting methods for simple turning points, perturbed bifurcation points, cubic turning points and cusps in boundary value problems for ordinary differential equations. We discretize the original problem using the shooting technique to obtain a finite dimensional problem. A direct method is used for the characetrization and computation of simple turning points. Some suitable extension of this direct method is employed for the computation of bifurcation points and cubic turning points and cusps. Enough numerical examples are solved to demonstrate that the method is efficient.     

A B-spline collocation method for solving the incompressible Navier-Stokes equations using an ad hoc method: the Boundary Residual method

Collocation methods using piece-wise polynomials, including B-splines, have been developed to find approximate solutions to both ordinary and partial differential equations. Such methods are elegant in their simplicity and efficient in their application. The spline collocation method is typically more efficient than traditional Galerkin finite element methods, which are used to solve the equations of fluid dynamics. The collocation method avoids integration. Exact formulae are available to find derivatives on spline curves and surfaces. The primary objective of the present work is to determine the requirements for the successful application of B-spline collocation to solve the coupled, steady, 2D, incompressible Navier-Stokes and continuity equations for laminar flow. The successful application of B-spline collocation included the development of ad hoc method dubbed the Boundary Residual method to deal with the presence of the pressure terms in the Navier-Stokes equations. Historically, other ad hoc methods have been developed to solve the incompressible Navier-Stokes equations, including the artificial compressibility, pressure correction and penalty methods. Convergence studies show that the ad hoc Boundary Residual method is convergent toward an exact (manufactured) solution for the 2D, steady, incompressible Navier-Stokes and continuity equations. C¹ cubic and quartic B-spline schemes employing orthogonal collocation and C² cubic and C³ quartic B-spline schemes with collocation at the Greville points are investigated. The C³ quartic Greville scheme is shown to be the most efficient scheme for a given accuracy, even though the C¹ quartic orthogonal scheme is the most accurate for a given partition. Two solution approaches are employed, including a globally-convergent zero-finding Newton's method using an LU decomposition direct solver and the variable-metric minimization method using BFGS update.     

Direct and indirect boundary element methods for solving the heat conduction problem

The boundary element method is used to solve the stationary heat conduction problem as a Dirichlet, a Neumann or as a mixed boundary value problem. Using singularities which are interpreted physically, a number of Fredholm integral equations of the first or second kind is derived by the indirect method. With the aid of Green's third identity and Kupradze's functional equation further direct integral equations are obtained for the given problem. Finally a numerical method is described for solving the integral equations using Hermitian polynomials for the boundary elements and constant, linear, quadratic or cubic polynomials for the unknown functions.     

Cubic Spline Collocation for Volterra Integral Equations

In the standard step-by-step cubic spline collocation method for Volterra integral equations an initial condition is replaced by a not-a-knot boundary condition at the other end of the interval. Such a method is stable in the same region of collocation parameter as in the step-by-step implementation with linear splines. The results about stability and convergence are based on the uniform boundedness of corresponding cubic spline interpolation projections. The numerical tests given at the end completely support the theoretical analysis. Received: January 15, 2002; revised July 27, 2002 Published online: December 19, 2002