势问题的径向基函数——局部边界积分方程方法

将基于径向基函数构造的具有插值特性的近似函数和局部边界积分方程方法相结合，建立了求解势问题的径向基函数——局部边界积分方程方法，推导了相应离散方程.与其他边界积分方程的无网格方法相比，本文方法具有数值实现过程简单、计算量小、精度高的优点，并可直接施加边界条件.最后通过算例说明了该方法的有效性. 关键词： 径向基函数 无网格方法 局部边界积分方程 势问题     

Shape and topology optimization for electrothermomechanical microactuators using level set methods

In this short note, a shape and topology optimization method is presented for multiphysics actuators including geometrically nonlinear modeling based on an implicit free boundary parameterization method. A level set model is established to describe structural design boundary by embedding it into the zero level set of a higher-dimensional level set function. The compactly supported radial basis functions (CSRBF) are introduced to parameterize the implicit level set surface with a high level of accuracy and smoothness. The original more difficult shape and topology optimization driven by the Hamilton–Jacobi partial differential equation (PDE) is transferred into a relatively easier parametric (size) optimization, to which many well-founded optimization algorithms can be applied. Thus the structural optimization is transformed to a numerical process that describes the design as a sequence of motions of the design boundaries by updating the expansion coefficients of the size optimization. Two widely studied examples are chosen to demonstrate the effectiveness of the proposed method.     

基于基函数展开方法处理具有陡峭界面的Abel反演问题

在激光聚变研究中,存在大量具有陡峭不连续界面特性的靶丸等离子体诊断问题,需要借助Abel反演方法对结果进行分析。 提出了一种基于基函数展开处理不连续径向分布的Abel反演方法,采用带约束的Landweber迭代法求得优化结果。利用该算法可以对不连续阶跃分布进行比较准确的反演,获得阶跃结构特征。在理想情况下,针对不连续阶跃分布的反演方差好于10-4,即使在存在一定噪声的情况下也仍然可以获得较好的结果。该方法可用于对激光聚变研究中压缩过程中球壳的背光照相、自发光成像、光谱测量等诊断结果的数据分析。     

Adiabatic representation for a hydrogen atom photoionization in a uniform magnetic field

A new effective method of calculating wave functions of discrete and continuous spectra of a hydrogen atom in a strong magnetic field is developed on the basis of the adiabatic approach to parametric eigenvalue problems in spherical coordinates. The two-dimensional spectral problem for the Schrödinger equation at a fixed magnetic quantum number and parity is reduced to a spectral parametric problem for a one-dimensional angular equation and a finite set of ordinary second-order radial differential equations. The results are in good agreement with the photoionization calculations by other authors and have a true threshold behavior.     

Semiclassical approximation for the twodimensional Fisher–Kolmogorov–Petrovskii– Piskunov equation with nonlocal nonlinearity in polar coordinates

The two-dimensional Kolmogorov–Petrovskii–Piskunov–Fisher equation with nonlocal nonlinearity and axially symmetric coefficients in polar coordinates is considered. The method of separation of variables in polar coordinates and the nonlinear superposition principle proposed by the authors are used to construct the asymptotic solution of a Cauchy problem in a special class of smooth functions. The functions of this class arbitrarily depend on the angular variable and are semiclassically concentrated in the radial variable. The angular dependence of the function has been exactly taken into account in the solution. For the radial equation, the formalism of semiclassical asymptotics has been developed for the class of functions which singularly depend on an asymptotic small parameter, whose part is played by the diffusion coefficient. A dynamic system of Einstein–Ehrenfest equations (a system of equations in mean and central moments) has been derived. The evolution operator for the class of functions under consideration has been constructed in explicit form.     

Shape and topology optimization of compliant mechanisms using a parameterization level set method

In this paper, a parameterization level set method is presented to simultaneously perform shape and topology optimization of compliant mechanisms. The structural shape boundary is implicitly embedded into a higher-dimensional scalar function as its zero level set, resultantly, establishing the level set model. By applying the compactly supported radial basis function with favorable smoothness and accuracy to interpolate the level set function, the temporal and spatial Hamilton–Jacobi equation from the conventional level set method is then discretized into a series of algebraic equations. Accordingly, the original shape and topology optimization is now fully transformed into a parameterization problem, namely, size optimization with the expansion coefficients of interpolants as a limited number of design variables.Design of compliant mechanisms is mathematically formulated as a general optimization problem with a nonconvex objective function and two additionally specified constraints. The structural shape boundary is then advanced as a process of renewing the level set function by iteratively finding the expansion coefficients of the size optimization with a sequential convex programming method. It is highlighted that the present method can not only inherit the merits of the implicit boundary representation, but also avoid some unfavorable features of the conventional discrete level set method, such as the CFL condition restriction, the re-initialization procedure and the velocity extension algorithm. Finally, an extensively investigated example is presented to demonstrate the benefits and advantages of the present method, especially, its capability of creating new holes inside the design domain.     

Determination of an Extremal in Two-Dimensional Variational Problems Based on the RBF Collocation Method

This paper introduces a direct method derived from the global radial basis function (RBF) interpolation over arbitrary collocation nodes occurring in variational problems involving functionals that depend on functions of a number of independent variables. This technique parameterizes solutions with an arbitrary RBF and transforms the two-dimensional variational problem (2DVP) into a constrained optimization problem via arbitrary collocation nodes. The advantage of this method lies in its flexibility in selecting between different RBFs for the interpolation and parameterizing a wide range of arbitrary nodal points. Arbitrary collocation points for the center of the RBFs are applied in order to reduce the constrained variation problem into one of a constrained optimization. The Lagrange multiplier technique is used to transform the optimization problem into an algebraic equation system. Three numerical examples indicate the high efficiency and accuracy of the proposed technique.     

耦合径向基函数与多项式基函数的无网格方法

耦合径向基函数和多项式基函数,形成一种新的近似函数.该近似函数对散乱分布的离散数据点进行逼近时,只需节点信息,不需要划分网格.详细描述了耦合近似函数的建立、属性、插值行为及其形函数和形函数导数的性质.最后引入修正变分原理和单位分解积分技术求解边值问题,并给出了计算实例,表明耦合径向基函数和多项式基函数是一种非常有效的方法.     

Shape and topology optimization based on the phase field method and sensitivity analysis

This paper discusses a structural optimization method that optimizes shape and topology based on the phase field method. The proposed method has the same functional capabilities as a structural optimization method based on the level set method incorporating perimeter control functions. The advantage of the method is the simplicity of computation, since extra operations such as re-initialization of functions are not required. Structural shapes are represented by the phase field function defined in the design domain, and optimization of this function is performed by solving a time-dependent reaction diffusion equation. The artificial double well potential function used in the equation is derived from sensitivity analysis. The proposed method is applied to two-dimensional linear elastic and vibration optimization problems such as the minimum compliance problem, a compliant mechanism design problem and the eigenfrequency maximization problem. The numerical examples provided illustrate the convergence of the various objective functions and the effect that perimeter control has on the optimal configurations.     

A Particle-Mesh Method for the Shallow Water Equations Near Geostrophic Balance

In this paper we outline a new particle-mesh method for rapidly rotating shallow water flows based on a set of regularized equations of motion. The time-stepping method uses an operator splitting of the equations into an Eulerian gravity wave part and a Lagrangian advection part. An essential ingredient is the advection of absolute vorticity by means of translated radial basis functions. We show that this implies exact conservation of enstrophy. The method is tested on two model problems based on the qualitative features of the solutions obtained (i.e., dispersion or smoothness of potential vorticity contours) as well as on the increase in mean divergence level.     

截断展开方法和广义变系数KdV方程新的精确类孤子解

利用特殊的截断展开方法求出了广义变系数KdV方程新的类孤子解.这种方法的基本思想是假定形式解具有截断展开形式,以致可把广义变系数KdV方程转化为一组待定函数的代数方程组,进而给出待定函数容易积分的常微分方程.利用例子证明了这种方法是十分有效的. 关键词： 截断展开法 变系数 KdV方程 孤波解     

Ultrasonic flaw detection using radial basis function networks (RBFNs)

Ultrasonic flaw detection has been studied many times in the literature. Schemes based on thresholding after a previous matched filter use to be the best solution, but results obtained with this method are only satisfactory when scattering and attenuation are not considered. In this paper, we propose an alternative solution to thresholding detection method. We deal with the usage of different flaw detection methods comparing them with the proposed one. The experiment tries to determinate whether a given ultrasonic signal contains a flaw echo or not. Starting with a set of 24,000 patterns with 750 samples each one, two subsets are defined for the experiments. The first one, the training set, is used to obtain the detection parameters of the different methods, and the second one is used to test the performance of them. The proposed method is based on radial basis functions networks, one of the most powerful neural network techniques. This signal processing technique tries to find the optimal decision criterion. Comparing this method with thresholding based ones, an improvement over 25-30% is obtained, depending on the probability of false alarm. So our new method is a good alternative to flaw detection problem.     

Adaptive meshless centres and RBF stencils for Poisson equation

We consider adaptive meshless discretisation of the Dirichlet problem for Poisson equation based on numerical differentiation stencils obtained with the help of radial basis functions. New meshless stencil selection and adaptive refinement algorithms are proposed in 2D. Numerical experiments show that the accuracy of the solution is comparable with, and often better than that achieved by the mesh-based adaptive finite element method.     

一种基于压缩感知的三维导体目标电磁散射问题的快速求解方法

提出了一种将压缩感知和特征基函数结合的方法来计算三维导体目标的雷达散射截面.利用压缩感知理论,将随机选择的矩量法阻抗矩阵作为测量矩阵,将激励电压视为测量值,然后再用恢复算法可实现二维或二维半目标感应电流的求解.对于三维导体目标,使用Rao-Wilton-Glisson基函数表示的感应电流在常用的离散余弦变换基、小波基等稀疏基上不稀疏.为此,本文将计算出的目标特征基函数作为稀疏基,用广义正交匹配追踪算法作为恢复算法来加速恢复过程,并应用到三维导体目标的雷达散射截面计算中.数值结果证明了本文方法的准确性与高效性.     

Modified level set method with Canny operator for image noise removal

The level set method is commonly used to address image noise removal. Existing studies concentrate mainly on determining the speed function of the evolution equation. Based on the idea of a Canny operator, this letter introduces a new method of controlling the level set evolution, in which the edge strength is taken into account in choosing curvature flows for the speed function and the normal to edge direction is used to orient the diffusion of the moving interface. The addition of an energy term to penalize the irregularity allows for better preservation of local edge information. In contrast with previous Canny-based level set methods that usually adopt a two-stage framework, the proposed algorithm can execute all the above operations in one process during noise removal.     

Cutoff frequencies in planar optical waveguides with arbitrary index profiles: An efficient numerical method

A numerical method for obtaining mode cutoffs for planar waveguides with arbitrary index profiles is developed. The method is based on the Galerkin method in which the wave equation for modes at cutoff is converted to a matrix eigenvalue equation using a set of orthogonal basis functions. Due to different boundary conditions, we have identified two separate cases; one, in which the cover and the substrate indices are identical leading to same behavior of the field at cutoff in these two regions and, the other, in which the two indices are different and hence, the field behaves differently. We have accordingly chosen different basis functions for the two cases. The method results in a generalized matrix eigenvalue problem which has been converted to a standard symmetric matrix eigenvalue analytically. The method has been used to obtain mode cutoffs for waveguides with a variety of index profiles. Comparisons with available exact results show that very good accuracies can be obtained with moderate matrix sizes.     

A level set method for the semiclassical limit of the Schrödinger equation with discontinuous potentials

We propose a level set method for the semiclassical limit of the Schrödinger equation with discontinuous potentials. The discontinuities in the potential corresponds to potential barriers, at which incoming waves can be partially transmitted and reflected. Previously such a problem was handled by Jin and Wen using the Liouville equation – which arises as the semiclassical limit of the Schrödinger equation – with an interface condition to account for partial transmissions and reflections (S. Jin, X. Wen, SIAM J. Num. Anal. 44 (2006) 1801–1828). However, the initial data are Dirac-delta functions which are difficult to approximate numerically with a high accuracy. In this paper, we extend the level set method introduced in (S. Jin, H. Liu, S. Osher, R. Tsai, J. Comp. Phys. 210 (2005) 497–518) for this problem. Instead of directly discretizing the Delta functions, our proposed method decomposes the initial data into finite sums of smooth functions that remain smooth in finite time along the phase flow, and hence can be solved much more easily using conventional high order discretization schemes.     

Ein neues Verfahren zur Berechnung von Einelektronenzuständen in Kristallen

A new method is presented for the calculation of wave-functions and energy-bands of electrons moving in a periodic potential in the one-electron approximation, which is based on different representations of the wave-functions. The wave-functions are expanded in spherical harmonies and radial solutions of the wave equation within a sphere inscribed in the atomic polyhedron. By suitable choice of the radial functions it is possible to obtain wave-functions which are well orthogonal to the functions of the core electrons. Outside the sphere the wave-functions are expanded in plane waves, having the correct translation symmetry. In contrast to other methods using similar expansions, the wave-functions and their first derivatives in our method are continuous in the whole space. These trial functions are introduced into a variational principle equivalent to the Schrödinger equation in which the conditions of continuity of the wave-functions and their first derivatives are taken into account by the method of undetermined multipliers. By an additional variation of the radial functions inside the sphere one may obtain functions which approximate the real functions leaving only an uniform error in all space.     

The use of PDE centres in the local RBF Hermitian method for 3D convective-diffusion problems

In this work an extension is proposed to the Local Hermitian Interpolation (LHI) method; a meshless numerical method based on interpolation with small and heavily overlapping radial basis function (RBF) systems. This extension to the LHI method uses interpolation functions which themselves satisfy the partial differential equation (PDE) to be solved. In this way, a much improved reconstruction of partial derivatives can be obtained, resulting in significantly improved accuracy in many cases.     

Variationally optimized numerical orbitals for molecular calculations

The problem of variational optimization of the atomic orbitals used in molecular calculations is investigated. It is shown that the variational principle leads to an equation similar to the radial Schrodinger equation but containing an inhomogeneous term. As an example, the equations are solved for the minimum basis set orbitals for the methane molecule. The results show a substantial improvement over those of a previous calculation optimizing in a minimum basis of Slater orbitals.