Analysis and computational method based on quadratic B‐spline FEM for the Rosenau–Burgers equation

L_∞‐error estimates for B‐spline Galerkin finite element solution of the Rosenau–Burgers equation are considered. The semidiscrete B‐spline Galerkin scheme is studied using appropriate projections. For fully discrete B‐spline Galerkin scheme, we consider the Crank–Nicolson method and analyze the corresponding error estimates in time. Numerical experiments are given to demonstrate validity and order of accuracy of the proposed method. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 877–895, 2016     

B‐spline collocation algorithm for numerical solution of the generalized Burger's‐Huxley equation

The cubic B‐spline collocation scheme is implemented to find numerical solution of the generalized Burger's–Huxley equation. The scheme is based on the finite‐difference formulation for time integration and cubic B‐spline functions for space integration. Convergence of the scheme is discussed through standard convergence analysis. The proposed scheme is of second‐order convergent. The accuracy of the proposed method is demonstrated by four test problems. The numerical results are found to be in good agreement with the exact solutions. Results are compared with other results given in literature. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

广义非参数回归的B样本贝叶斯估计

本文主要研究广义非参数模型B样条Bayes估计 .将回归函数按照B样条基展开 ,我们不具体选择节点的个数 ,而是节点个数取均匀的无信息先验 ,样条函数系数取正态先验 ,用B样条函数的后验均值估计回归函数 .并给出了回归函数B样条Bayes估计的MCMC的模拟计算方法 .通过对Logistic非参数回归的模拟研究 ,表明B样条Bayes估计得到了很好的估计效果     

Sommerfeld Diffraction Problems with Oblique Derivatives

A class of diffraction problems for the two-dimensional Helmholtz equation is studied where the boundary conditions on a canonical screen contain oblique derivatives. The problem in weak formulation is reduced to a system of pseudodifferential equations on a half-line. The Fredholm property is characterized in dependence of the parameters, index formulas are given, and, if necessary, a normalization of the problem is implemented. © 1997 by B.G. Teubner Stuttgart–John Wiley & Sons Ltd.     

Four techniques based on the B‐spline expansion and the collocation approach for the numerical solution of the Lane–Emden equation

Four numerical techniques based on the linear B‐spline functions are presented for the numerical solution of the Lane–Emden equation. Some properties of the B‐spline functions are presented and are utilized to reduce the solution of the Lane–Emden equation to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the new techniques. Copyright © 2013 John Wiley & Sons, Ltd.     

Exact Fourier expansion in cylindrical coordinates for the three-dimensional Helmholtz Green function

A new method is presented for Fourier decomposition of the Helmholtz Green function in cylindrical coordinates, which is equivalent to obtaining the solution of the Helmholtz equation for a general ring source. The Fourier coefficients of the Green function are split into their half advanced + half retarded and half advanced–half retarded components, and closed form solutions for these components are then obtained in terms of a Horn function and a Kampé de Fériet function respectively. Series solutions for the Fourier coefficients are given in terms of associated Legendre functions, Bessel and Hankel functions and a hypergeometric function. These series are derived either from the closed form 2-dimensional hypergeometric solutions or from an integral representation, or from both. A simple closed form far-field solution for the general Fourier coefficient is derived from the Hankel series. Numerical calculations comparing different methods of calculating the Fourier coefficients are presented. Fourth order ordinary differential equations for the Fourier coefficients are also given and discussed briefly.     

The uniqueness and existence of solutions for the 3D Helmholtz equation in a step‐index waveguide with unbounded perturbation

In this paper, we study the 3D Helmholtz equation in a step‐index waveguide with unbounded perturbation, allowing the presence of guided waves. Our assumptions on the perturbed and source terms are too few. On the basis of the Green's function for the 3D homogeneous Helmholtz equation in a step‐index waveguide without perturbation, we introduce a generalized (out‐going) Sommerfeld–Rellich radiation condition, and then we prove the uniqueness and existence of solutions for the studied 3D Helmholtz equation satisfying our radiation condition. Copyright © 2012 John Wiley & Sons, Ltd.     

A multigrid solver to the Helmholtz equation with a point source based on travel time and amplitude

The Helmholtz equation arises when modeling wave propagation in the frequency domain. The equation is discretized as an indefinite linear system, which is difficult to solve at high wave numbers. In many applications, the solution of the Helmholtz equation is required for a point source. In this case, it is possible to reformulate the equation as two separate equations: one for the travel time of the wave and one for its amplitude. The travel time is obtained by a solution of the factored eikonal equation, and the amplitude is obtained by solving a complex‐valued advection–diffusion–reaction equation. The reformulated equation is equivalent to the original Helmholtz equation, and the differences between the numerical solutions of these equations arise only from discretization errors. We develop an efficient multigrid solver for obtaining the amplitude given the travel time, which can be efficiently computed. This approach is advantageous because the amplitude is typically smooth in this case and, hence, more suitable for multigrid solvers than the standard Helmholtz discretization. We demonstrate that our second‐order advection–diffusion–reaction discretization is more accurate than the standard second‐order discretization at high wave numbers, as long as there are no reflections or caustics. Moreover, we show that using our approach, the problem can be solved more efficiently than using the common shifted Laplacian multigrid approach.     

Local modification of Pythagorean-hodograph quintic spline curves using the B-spline form

The problems of determining the B–spline form of a C ² Pythagorean–hodograph (PH) quintic spline curve interpolating given points, and of using this form to make local modifications, are addressed. To achieve the correct order of continuity, a quintic B–spline basis constructed on a knot sequence in which each (interior) knot is of multiplicity 3 is required. C ² quintic bases on uniform triple knots are constructed for both open and closed C ² curves, and are used to derive simple explicit formulae for the B–spline control points of C ² PH quintic spline curves. These B-spline control points are verified, and generalized to the case of non–uniform knots, by applying a knot removal scheme to the Bézier control points of the individual PH quintic spline segments, associated with a set of six–fold knots. Based on the B–spline form, a scheme for the local modification of planar PH quintic splines, in response to a control point displacement, is proposed. Only two contiguous spline segments are modified, but to preserve the PH nature of the modified segments, the continuity between modified and unmodified segments must be relaxed from C ² to C ¹. A number of computed examples are presented, to compare the shape quality of PH quintic and “ordinary” cubic splines subject to control point modifications.     

Collocation with WEB–Splines

We describe a collocation method with weighted extended B–splines (WEB–splines) for arbitrary bounded multidimensional domains, considering Poisson’s equation as a typical model problem. By slightly modifying the B–spline classification for the WEB–basis, the centers of the supports of inner B–splines can be used as collocation points. This resolves the mismatch between the number of basis functions and interpolation conditions, already present in classical univariate schemes, in a simple fashion. Collocation with WEB–splines is particularly easy to implement when the domain boundary can be represented as zero set of a weight function; sample programs are provided on the website http://www.web-spline.de. In contrast to standard finite element methods, no mesh generation and numerical integration is required, regardless of the geometric shape of the domain. As a consequence, the system equations can be compiled very efficiently. Moreover, numerical tests confirm that increasing the B–spline degree yields highly accurate approximations already on relatively coarse grids. Compared with Ritz-Galerkin methods, the observed convergence rates are decreased by 1 or 2 when using splines of odd or even order, respectively. This drawback, however, is outweighed by a substantially smaller bandwidth of collocation matrices.     

Helmholtz方程周期Green函数及其偏导数截断误差收敛阶的分析

在应用边界元方法求解Helmholtz方程周期边值问题时,需要构造以周期Green函数或其偏导数为核函数的积分算子形式的解.由于Helmholtz方程的周期Green函数G~P是一个函数项级数,该级数的通项是Hankel函数,在数值求解中,需要对其进行截断,从而很有必要研究其截断误差.本文根据Hankel函数在变量趋于无穷大时的渐近展开式,并结合Abel不等式,证明了G~P及其一阶偏导和二阶混合偏导一致收敛,且其截断误差收敛阶均为O(1/p~(1/2)).最后,通过数值实验验证了理论证明的正确性.本文的证明方法也可被用于证明其它一些方程周期Green函数的收敛性问题.     

B‐spline collocation algorithms for numerical solution of the RLW equation

Both sextic and septic B‐spline collocation algorithms are presented for the numerical solutions of the RLW equation. Numerical results resolve the fine structure of the single solitary wave propagation, two and three solitary waves interaction, and evolution of solitary waves. Comparison of the numerical results is done by the results of some earlier schemes mentioned in the article. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 581–607, 2011     

Optimal Design of Periodic Diffractive Structures

The problem of designing a periodic interface between two different materials, which gives rise to a specified far-field diffraction pattern for a given incoming plane wave, is considered. The time harmonic waves are assumed to be TM (transverse magnetic) polarized. The diffraction problem is modeled by a generalized Helmholtz equation with transparent boundary conditions. In this paper the design problem is relaxed to include highly oscillatory profiles. Existence of an optimal design is established. The principal method is based on the theory of homogenization for the model equation. Accepted 31 May 2000. Online publication 26 February 2001.     

Impedance eduction from energy dissipation measurements with flow

An optimization method for the calculation of the wall impedance of a perforate from measurements under grazing flow conditions is presented. It is in particular characterized by the simultaneous eduction of the parameter set defining the impedance as function of frequency, which allows interpolation and extrapolation. The Extended Helmholtz Resonator is applied as impedance model and the simultaneous eduction is performed using a time–domain CAA method. Furthermore, the method uses processed energy data from experiments to define the objective function. The results give an insight in the effect of bias flow on the wall impedance under grazing flow conditions. It proves to be an important acoustic design parameter. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Error estimate and extrapolation of a quadrature formula derived from a quartic spline quasi-interpolant

In this paper we analyze a quadrature rule based on integrating a C ³ quartic spline quasi-interpolant on a bounded interval which has been introduced in Sablonnière (Rend. Semin. Mat. Univ. Pol. Torino 63(3):107–118, 2005). By studying the sign structure of its associated Peano kernel we derive an explicit formula of the quadrature error with an approximation order O(h ⁶). A comparison of this rule with the composite Boole’s and the three-point Gauss-Legendre rules is given. We also compare the Nyström methods associated with the above quadrature formulae for solving the linear Fredholm integral equation of the second kind. Then, by combining the proposed rule with composite Boole’s rule, we construct a new quadrature rule of order O(h ⁷). All the obtained results are illustrated by several numerical tests.     

Quartic B‐spline collocation algorithms for numerical solution of the RLW equation

The collocation method based on quartic B‐spline interpolation is studied for numerical solution of the regularized long wave (RLW) equation. The time‐split RLW equation is also solved with the quartic B‐spline collocation method. Numerical accuracy is tested by obtaining the single solitary wave solution. Then, interaction, undulation and evolution of solitary waves are studied. Solutions are compared with available results. Conservation quantities are computed for all test experiments. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

Nonlocal Boundary Conditions for Higher–Order Parabolic Equations

This work deals with the efficient numerical solution of the two–dimensional one–way Helmholtz equation posed on an unbounded domain. In this case one has to introduce artificial boundary conditions to confine the computational domain. Here we construct with the Z –transformation so–called discrete transparent boundary conditions for higher–order parabolic equations schemes. These methods are Padé “Parabolic” approximations of the one–way Helmholtz equation and frequently used in integrated optics and (underwater) acoustics. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

B‐spline solution of fractional integro partial differential equation with a weakly singular kernel

The main objective of the paper is to find the approximate solution of fractional integro partial differential equation with a weakly singular kernel. Integro partial differential equation (IPDE) appears in the study of viscoelastic phenomena. Cubic B‐spline collocation method is employed for fractional IPDE. The developed scheme for finding the solution of the considered problem is based on finite difference method and collocation method. Caputo fractional derivative is used for time fractional derivative of order α, . The given problem is discretized in both time and space directions. Backward Euler formula is used for temporal discretization. Collocation method is used for spatial discretization. The developed scheme is proved to be stable and convergent with respect to time. Approximate solutions are examined to check the precision and effectiveness of the presented method.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1565–1581, 2017     

On constructing the solution to the exterior Dirichlet problem for the Helmholtz equation

In this paper a method is given for constructing the solution to the exterior Dirichlet problem for the Helmholtz equation in three dimensions. This method is modeled after the procedure of Colton and Kleinman (Proc. Roy. Soc. Edin. 86A(1980), 29–42) for solving the corresponding two-dimensional problem. The scattering problem is reformulated as an integral equation and it is shown that its solution can be represented as a convergent Neumann series for small values of the wave number. Comparisons are made between the present method and known results. Examples are given which illustrate the method.     

Structural optimization of an acoustic horn

The objective of the present article is to find an optimal design of an acoustic horn in the case that the magnitude of the reflection wave integrated over the inflow boundary is to be minimized meanwhile the Helmholtz equation models the wave propagation. In contrast to the current approaches such as gradient-based optimization algorithms, we employ here a non-iterative method based on measure theory which dose not require any information of gradients and the differentiability of objective function in the optimization problem is not as a rule. Implementation of our fast convergence approach shows that the resulting horns, not only for single frequency optimization but also for a band of frequencies, are very efficient.