Numerical computation of nonsimple turning points and cusps

Summary A method is developed for the numerical computation of a double turning point corresponding to a cusp catastrophe of a nonlinear operator equation depending on two parameters. An augmented system containing the original equation is introduced, for which the cusp point is an isolated solution. An efficient implementation of Newton's method in the finite-dimensional case is presented. Results are given for some chemical engineering problems and this direct method is compared with some other techniques to locate cusp points.     

Numerical computation of branch points in nonlinear equations

Summary The numerical computation of branch points in systems of nonlinear equations is considered. A direct method is presented which requires the solution of one equation only. The branch points are indicated by suitable testfunctions. Numerical results of three examples are given.     

椭圆外区域上Helmholtz问题的自然边界元法

本文研究椭圆外区域上Helmholtz方程边值问题的自然边界元法.利用自然边界归化原理,获得该问题的Poisson积分公式及自然积分方程,给出了自然积分方程的数值方法.由于计算的需要,我们详细地讨论了Mathieu函数的计算方法(当0相似文献   

广义非线性Sin-Gordon方程的整体解及数值计算

本文考察了一类广义非线性Sin-Gordon方程的周期初值问题，利用非线性Galerkin方法，证明了其整体解的存在性和唯一性，并给出了其有界吸引集的存在性．构造了全离散的Fourier拟谱显格式，利用有界延拓法证明了其格式的收敛性与稳定性，并给出了误差估计、算法分析及计算复杂度，最后，通过数值例子，检验了理论结果的可信性．为对此模型的数值分析提供了理论基础和一个有效的算法．     

Sobolev方程的全离散有限体积元格式及数值模拟

本文研究二维Sobolev方程的有限体积元方法, 给出一种全离散化有限体积元格式及其有限体积元解的误差估计,并用数值例子说明数值计算的结果与理论结果是相吻合的, 进一步说明了有限体积元方法比其他数值方法更优越.     

Componentwise inclusion and exclusion sets for solutions of quadratic equations in finite dimensional spaces

Summary In this paper we use interval arithmetic tools for the computation of componentwise inclusion and exclusion sets for solutions of quadratic equations in finite dimensional spaces. We define a mapping for which under certain assumptions we can construct an interval vector which is mapped into itself. Using Brouwer's fixed point theorem we conclude the existence of a solution of the original equation in this interval vector. Under different assumptions we can construct an interval vector such that the range of the mapping has no point in common with this interval vector. This implies that there is no solution in this interval vector. Furthermore we consider an iteration method which improves componentwise errorbounds for a solution of a quadratic. The theoretical results of this paper are demonstrated by some numerical examples using the algebraic eigenvalue problem which is probably the best known example of a quadratic equation.This paper contains the main results of a talk given by the author on the occasion of the 25th anniversary of the founding of Numerische Mathematik, March 19–21, 1984 at the Technische Universität of Munich, Germany     

Discretization of the stationary convection-diffusion-reaction equation

A finite volume method for the convection-diffusion-reaction equation is presented, which is a model equation in combustion theory. This method is combined with an exponential scheme for the computation of the fluxes. We prove that the numerical fluxes are second-order accurate, uniformly in the local Peclet numbers. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 607–625, 1998     

Versions of inexact Kleinman-Newton methods for Riccati equations

Optimal control problems involving PDEs often lead in practice to the numerical computation of feedback laws for an optimal control. This is achieved through the solution of a Riccati equation which can be large scale, since the discretized problems are large scale and require special attention in their numerical solution. The Kleinman-Newton method is a classical way to solve an algebraic Riccati equation. We look at two versions of an extension of this method to an inexact Newton method. It can be shown that these two implementable versions of Newton's method are identical in the exact case, but differ substantially for the inexact Newton method. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Wave generation in a computation domain

One of the possible methods to deal with the reflecting waves at the incident boundary in numerical modeling is to generate waves in the computation domain and absorb the outgoing waves at the incident boundary. A source function is introduced into the momentum equation of Boussinesq equations for generating wave in a computation domain in this paper. Typical numerical examples are given for the verification of the proposed method. Numerical examination for the wave diffraction through a breakwater gap shows that the proposed method is especially useful for multidirectional waves.     

Numerical proof of stability of roll waves in the small-amplitude limit for inclined thin film flow

We present a rigorous numerical proof based on interval arithmetic computations categorizing the linearized and nonlinear stability of periodic viscous roll waves of the KdV–KS equation modeling weakly unstable flow of a thin fluid film on an incline in the small-amplitude KdV limit. The argument proceeds by verification of a stability condition derived by Bar–Nepomnyashchy and Johnson–Noble–Rodrigues–Zumbrun involving inner products of various elliptic functions arising through the KdV equation. One key point in the analysis is a bootstrap argument balancing the extremely poor sup norm bounds for these functions against the extremely good convergence properties for analytic interpolation in order to obtain a feasible computation time. Another is the way of handling analytic interpolation in several variables by a two-step process carving up the parameter space into manageable pieces for rigorous evaluation. These and other general aspects of the analysis should serve as blueprints for more general analyses of spectral stability.     

A Numerical Method of Local Energy Decay for the Boundary Controllability of Time-Reversible Distributed Parameter Systems

This paper deals with the numerical computation of the boundary controls of linear, time-reversible, second-order evolution systems. Based on a method introduced by Russell ( Stud. Appl. Math. LII(3) (1973)) for the wave equation, a numerical algorithm is proposed for solving this type of problems. The convergence of the method is based on the local energy decay of the solution of a suitable Cauchy problem associated with the original control system. The method is illustrated with several numerical simulations for the Klein–Gordon and the Euler–Bernoulli equations in 1D, the wave equation on a rectangle, and the plate equation on a disk.     

An upwind center difference parallel method and numerical analysis for the displacement problem with moving boundary

A nonlinear coupled mathematical system of two‐phase seepage flow displacement is discussed in this paper including an elliptic equation for the pressure and a convection‐dominated diffusion equation for the saturation. In fact, the boundary of an underground region where the fluid flows through is nonstationary. So a moving boundary should be considered. The saturation equation is convection‐dominated, therefore the method of upwind finite difference is introduced for the accurate computation. The upwind approximation could eliminate numerical oscillation and strong stability is shown. Since the computational work of saturation is larger than the pressure, the authors apply a parallel method, decomposing the whole domain into several nonoverlapping subdomains, to simplify the computation. A domain decomposition method coupled with upwind differences is presented for the saturation. The pressure equation is discretized by a five‐point center finite difference method. By using a transformation and defining new inner products and norms, error estimates in l² norm is discussed. Finally, two experimental tests are given to illustrate the efficiency and accuracy of the parallel algorithm.     

Numerical Solution for the Helmholtz Equation with Mixed Boundary Condition

We consider the numerical solution for the Helmholtz equation in R~2 with mixed boundary conditions.The solvability of this mixed boundary value problem is estab- lished by the boundary integral equation method.Based on the Green formula,we express the solution in terms of the boundary data.The key to the numerical real- ization of this method is the computation of weakly singular integrals.Numerical performances show the validity and feasibility of our method.The numerical schemes proposed in this paper have been applied in the realization of probe method for inverse scattering problems.     

微分本构粘弹性轴向运动弦线横向振动分析的差分法

给出了微分本构粘弹性轴向运动弦线横向振动数值仿真的一种差分法．文中建立了具有微分本构的粘弹性运动弦线的横向振动模型;通过对系统的控制方程和本构方程在不同的分数节点离散,得到一种新的差分方法．利用这一方法,弦线振动方程的数值计算过程可以交替地显式进行,且有较小的截断误差和好的数值稳定性．与通用的方法比较,新的方法计算简单、方便．文中利用方程的不变量检验了数值结果的可靠性,并利用这一方法给出了一类弦线模型的参数振动分析．     

The computation of consistent initial values for nonlinear index-2 differential–algebraic equations

The computation of consistent initial values for differential–algebraic equations (DAEs) is essential for starting a numerical integration. Based on the tractability index concept a method is proposed to filter those equations of a system of index-2 DAEs, whose differentiation leads to an index reduction. The considered equation class covers Hessenberg-systems and the equations arising from the simulation of electrical networks by means of Modified Nodal Analysis (MNA). The index reduction provides a method for the computation of the consistent initial values. The realized algorithm is described and illustrated by examples.     

Wavelet method for a class of fractional convection-diffusion equation with variable coefficients

A wavelet method to the solution for a class of space–time fractional convection-diffusion equation with variable coefficients is proposed, by which combining Haar wavelet and operational matrix together and dispersing the coefficients efficaciously. The original problem is translated into Sylvester equation and computation became convenient. The numerical example shows that the method is effective.     

Compactly supported wavelet and the numerical solution of the vlasov equation

A new scheme for solving the Vlasov equation using a compactly supported wavelets basis is proposed. We use a numerical method which minimizes the numerical diffusion and conserves a reasonable time computing cost. So we introduce a representation in a compactly supported wavelet of the derivative operator. This method makes easy and simple the computation of the coefficients of the matrix representing the operator. This allows us to solve the two equations which result from the splitting technique of the main Vlasov equation. Some numerical results are exposed using different numbers of wavelets.     

Numerical analysis of the weighted particle method applied to the semiconductor Boltzmann equation

Summary The semiconductor Boltzmann equation involves an integral operator, the kernel of which is a measure supported by a surface. This feature introduces some singularities of the exact solution, which makes the numerical approximation of this equation difficult. This paper is devoted to the error analysis of the weighted particle method (introduced by Mas-Gallic and Raviart [14]) applied to the space homogeneous semiconductor Boltzmann equation. The results are commented in view of the practical use of the method. This paper is closely related to [12], where results of numerical simulations on both test and real problems are given.     

Non-separable splines and numerical computation of evolution equations by the Galerkin methods

We construct new non-separable splines and we apply the spline sampling approximation to the computation of numerical solutions of evolution equations. The non-separable splines are basis functions which give a fine sampling approximation which enables us to compute numerical solutions by means of the method of lines combined with the Galerkin method. To demonstrate our approach we compute numerical solutions of the Burgers equation and the Kadomtsev–Petviashvili equation.     

p-Adic Analogue of the Wave Equation

In this paper, a p-adic analogue of the wave equation with Lipschitz source is considered. Since it is hard to arrive the solution of the problem, we propose a regularized method to solve the problem from a modified p-adic integral equation. Moreover, we give an iterative scheme for numerical computation of the regularlized solution.