Spline R-Function and Applications in FEM

R-function is a widely used tool when considering objects obtained through the Boolean operations start from simple base primitives. However, there is square root operation in the representation. Considering that the use of splines will facilitate the calculations within the CAD system, in this paper, we propose a system of R-functions represented in spline form called Spline R-function (SR). After transforming the function ranges of two base primitives to a new coordinate system, a series of sign constraints following a specific Boolean operation are derived and the spline R-function can be formulated as a piecewise function. Representation of SR in both Bézier form and B-spline form have been given. Among which the Bézier ordinates are determined with the help of the B-net method through setting up a series of relations according to the sign constraints and properties of R-functions. The construction processes for both Boolean intersection and union operations with different smoothness are discussed in detail. Numerical experiments are conducted to show the potential of the proposed spline R-function.     

平行电磁场中的Rydberg锂原子吸收谱的模型势计算

用B样条基组展开方法结合模型势计算了Rydberg锂原子在平行电磁场下的振子强度谱. 径向和角向均采用高阶Ｂ样条基组.计算结果与已有的实验结果符合得很好.利用分波分析法,对部分谱线的振子强度的强弱进行了分析. 本文方法简单有效，易于推广到交叉电磁场中Rydberg原子的精确谱的计算. 关键词： B样条 能谱 振子强度 平行电磁场 Rydberg原子     

Analytical and Approximate Solutions for Complex Nonlinear Schrödinger Equation via Generalized Auxiliary Equation and Numerical Schemes

This article studies the performance of analytical, semi-analytical and numerical scheme on the complex nonlinear Schr¨odinger(NLS) equation. The generalized auxiliary equation method is surveyed to get the explicit wave solutions that are used to examine the semi-analytical and numerical solutions that are obtained by the Adomian decomposition method, and B-spline schemes(cubic, quantic, and septic). The complex NLS equation relates to many physical phenomena in different branches of science like a quantum state, fiber optics, and water waves. It describes the evolution of slowly varying packets of quasi-monochromatic waves, wave propagation, and the envelope of modulated wave groups, respectively. Moreover, it relates to Bose-Einstein condensates which is a state of matter of a dilute gas of bosons cooled to temperatures very close to absolute zero. Some of the obtained solutions are studied under specific conditions on the parameters to constitute and study the dynamical behavior of this model in two and three-dimensional.     

B-spline solver for one-electron Schrödinger equation

A numerical algorithm for solving the one-electron Schrödinger equation is presented. The algorithm is based on the Finite Element method, and the basis functions are tensor products of univariate B-splines. The application of cubic or higher order B-splines guarantees that the searched solution belongs to a continuous and one time differentiable function space, which is a desirable property in the Kohn–Sham equation context from the Density Functional Theory with pseudopotential approximation. The theoretical background of the numerical algorithm is presented, and additionally, the implementation on parallel computers with distributed memory is described. The current implementation of the algorithm uses the MPI, HYPRE and ParMETIS libraries to distribute matrices on processing units. Additionally, the LOBPCG algorithm from HYPRE library is used to solve the algebraic generalized eigenvalue problem. The proposed algorithm works for any smooth interaction potential, where the domain of the problem is a finite subspace of the ?³ space. The accuracy of the algorithm is demonstrated for a selected interaction potential. In the current stage, the algorithm can be applied to solve the linearized Kohn–Sham equation for molecular systems.     

Least squares approximation by splines with free knots

Suppose we are given noisy data which are considered to be perturbed values of a smooth, univariate function. In order to approximate these data in the least squares sense, a linear combination of B-splines is used where the tradeoff between smoothness and closeness of the fit is controlled by a smoothing term which regularizes the least squares problem and guarantees unique solvability independent of the position of knots. Moreover, a subset of the knot sequence which defines the B-splines, the so-calledfree knots, is included in the optimization process.The resulting constrained least squares problem which is linear in the spline coefficients but nonlinear in the free knots is reduced to a problem that has only the free knots as variables. The reduced problem is solved by a generalized Gauss-Newton method. The method developed can be combined with a knot removal strategy in order to obtain an approximating spline with as few parameters as possible.Dedicated to Professor Dr.-Ing. habil. Dr. h.c. Helmut Heinrich on the occasion of his 90th birthdayResearch of the second author was partly supported by Deutsche Forschungsgemeinschaft under grant Schm 968/2-1.     

A unified approach to B-spline recursions and knot insertion, with application to new recursion formulas

We present a unified approach to and a generalization of almost all known recursion schemes concerning B-spline functions. This includes formulas for the computation of a B-spline's values, its derivatives (ordinary and partial), and for a knot insertion method for B-spline curves. Furthermore, our generalization allows us to derive also some new relations for these purposes.     

A geometric nonuniform fast Fourier transform

An efficient algorithm is presented for the computation of Fourier coefficients of piecewise-polynomial densities on flat geometric objects in arbitrary dimension and codimension. Applications range from standard nonuniform FFTs of scattered point data, through line and surface potentials in two and three dimensions, to volumetric transforms in three dimensions. Input densities are smoothed with a B-spline kernel, sampled on a uniform grid, and transformed by a standard FFT, and the resulting coefficients are unsmoothed by division. Any specified accuracy can be achieved, and numerical experiments demonstrate the efficiency of the algorithm for a gallery of realistic examples.     

Numerical solution of Hammerstein integral equations by using combination of spline-collocation method and Lagrange interpolation

The numerical solutions to the nonlinear integral equations of the Hammerstein-type:

with using B-splines functions are investigated. An interpolation method based on B-splines functions combined with a new collocation method is presented. Finally, for showing efficiency of the method we give some numerical examples.     

Exponential B-spline collocation method for self-adjoint singularly perturbed boundary value problems

We present an exponential B-spline collocation method for self-adjoint singularly perturbed boundary value problem. The convergence analysis is given and the method is shown to have second order uniform convergence. Numerical experiments are conducted to demonstrate the efficiency of the method.     

软X光Dante谱仪的平滑解谱算法及其在Z箍缩诊断中的应用

针对自行研制的8通道Dante谱仪,提出了一种使用B样条曲线进行平滑解谱的算法,并用黑体谱对其进行了验证。计算结果表明：新算法与使用矩形函数解谱的算法相比,对能谱轮廓与总辐射强度的反演精度均有明显提高。给出了“阳”加速器氩气喷气Z箍缩实验的测量波形,由8通道Dante谱仪的测量结果解谱得到的软X光峰值辐射功率和总能量分别为40.5 GW和851 J,与软X光功率计经修正后的测量结果相一致。