An interpolating reproducing kernel particle method for two-dimensional scatter points

An interpolating reproducing kernel particle method for two-dimensional(2D) scatter points is introduced. It eliminates the dependency of gridding in numerical calculations. The interpolating shape function in the interpolating reproducing kernel particle method satisfies the property of the Kronecker delta function. This method offers a mathematics basis for recognition technology and simulation analysis, which can be expressed as simultaneous differential equations in science or project problems. Mathematical examples are given to show the validity of the interpolating reproducing kernel particle method.     

Improved reproducing kernel particle method for piezoelectric materials

In this paper, the normal derivative of the radial basis function(RBF) is introduced into the reproducing kernel particle method(RKPM), and the improved reproducing kernel particle method(IRKPM) is proposed. The method can decrease the errors on the boundary and improve the accuracy and stability of the algorithm. The proposed method is applied to the numerical simulation of piezoelectric materials and the corresponding governing equations are derived. The numerical results show that the IRKPM is more stable and accurate than the RKPM.     

The choosing of reproducing kernel particle shape function with mathematic proof

Many mechanical problems can be induced from differential equations with boundary conditions; there exist analytic and numerical methods for solving the differential equations. Usually it is not so easy to obtain analytic solutions. So it is necessary to give numerical solutions. The reproducing kernel particle (RKP) method is based on the Garlerkin Meshless method. According to the Sobolev space and Fourier transform, the RKP shape function is mathematically proved in this paper.     

A NOTE ON USING ADOMIAN DECOMPOSITION METHOD FOR SOLVING BOUNDARY VALUE PROBLEMS

In this paper we outline a reliable strategy to use Adomian decomposition method properly for solving nonlinear partial differential equations with boundary conditions. Our fundamental goal in this paper has two features: (i) it introduces an efficient way for using Adomian decomposition method for boundary value problems, and (ii) it also would present the framework in a general way so that it may be used in BVPs of the same type. A numerical example is included to dwell upon the importance of the analysis presented.     

Kernel method for nonlinear granger causality

Important information on the structure of complex systems can be obtained by measuring to what extent the individual components exchange information among each other. The linear Granger approach, to detect cause-effect relationships between time series, has emerged in recent years as a leading statistical technique to accomplish this task. Here we generalize Granger causality to the nonlinear case using the theory of reproducing kernel Hilbert spaces. Our method performs linear Granger causality in the feature space of suitable kernel functions, assuming arbitrary degree of nonlinearity. We develop a new strategy to cope with the problem of overfitting, based on the geometry of reproducing kernel Hilbert spaces. Applications to coupled chaotic maps and physiological data sets are presented.     

Dynamical systems in the theory of solitons in the presence of nonlocal interactions

The behavior of solitons in models which take into account complex dispersion or nonlocal interaction of nonlinear waves is examined. A method is proposed to reduce this problem to one involving special trajectories (homoclinic and heteroclinic) of the dynamic system. This method involves replacing the nonlinear integrodifferential equation with the differential equations which link the original nonlinear field with the auxiliary linear fields. The interaction of fields in such a model is a local interaction. The number of introduced linear fields is determined by the Laplace transform of the integral operator kernel of the basic integrodifferential equation. The problem involving topological solitons for the nonlocal generalization of the Klein-Gordon equation is considered. Nonlocal interactions are found to lead to a number of singularities (unrestricted increase in the slope of the topological soliton front, break in the solutions, and other singularities).     

瞬态热传导问题的复变量重构核粒子法

在重构核粒子法的基础上,引入复变量,讨论了复变量重构核粒子法.复变量重构核粒子法的优点是在构造形函数时采用一维基函数建立二维问题的修正函数.然后,将复变量重构核粒子法应用于瞬态热传导问题的求解,结合瞬态热传导问题的Galerkin积分弱形式,采用罚函数法引入本质边界条件,建立了瞬态热传导问题的复变量重构核粒子法,推导了相应的计算公式.与传统的重构核粒子法相比,复变量重构核粒子法具有计算量小、精度高的优点.最后通过数值算例证明了该方法的有效性. 关键词： 重构核粒子法 复变量重构核粒子法 修正函数 瞬态热传导问题     

A kernel-free boundary integral method for elliptic boundary value problems

This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green’s functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green’s functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong.     

Homotopy Continuation Approach to Electrostatic Boundary-Value Problems of Nonlinear Media

The homotopy continuation method is employed to solve electrostatic boundaryvalue problems of nonlinear media. The difficulty associated with matching the inherently nonlinear boundary conditions on the interface is overcome by the mode expansion method, by which the nonlinear partial differential equations of the original problem are transformed into an infinite set of nonlinear ordinary differential equations. In this regard, the homotopy method has to be modified to handle the nonlinear boundary conditions. As an illustration, we study two cases:(a) nonlinear inclusion in linear host and (b) linear inclusion-in nonlinear host, both in two dimensions. The homotopy method is validated by comparing the results with the exact solution of case (a) and the results derived by perturbation method in case (b).     

The variational iteration method for solving <Emphasis Type="Italic">n</Emphasis>-th order fuzzy differential equations

The variational iteration method (VIM) proposed by Ji-Huan He is a new analytical method for solving linear and nonlinear equations. In this paper, the variational iteration method has been applied in solving nth-order fuzzy linear differential equations with fuzzy initial conditions. This method is illustrated by solving several examples.     

弹性力学的重构核粒子边界无单元法

将重构核粒子法(RKPM)和边界积分方程方法结合，提出了一种新的边界积分方程无网格方法——重构核粒子边界无单元法(RKP-BEFM).对弹性力学问题，推导了其重构核粒子边界无单元法的公式，研究其数值积分方案，建立了重构核粒子边界无单元法离散化边界积分方程，并推导了重构核粒子边界无单元法的内点位移和应力积分公式.重构核粒子法形成的形函数具有重构核函数的光滑性，且能再现多项式在插值点的精确值，所以本方法具有更高的精度.最后给出了数值算例，验证了本方法的有效性和正确性. 关键词： 重构核粒子法 弹性力学 边界无单元法     

The complex variable reproducing kernel particle method for elasto-plasticity problems

On the basis of reproducing kernel particle method(RKPM),using complex variable theory,the complex variable reproducing kernel particle method(CVRKPM) is discussed in this paper.The advantage of the CVRKPM is that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is formed.Then the CVRKPM is applied to solve two-dimensional elasto-plasticity problems.The Galerkin weak form is employed to obtain the discretized system equation,the penal...     

Toward Analytic Solution of Nonlinear Differential Difference Equations via Extended Sensitivity Approach

In this paper an efficient computational method based on extending the sensitivity approach(SA) is proposed to find an analytic exact solution of nonlinear differential difference equations.In this manner we avoid solving the nonlinear problem directly.By extension of sensitivity approach for differential difference equations(DDEs),the nonlinear original problem is transformed into infinite linear differential difference equations,which should be solved in a recursive manner.Then the exact solution is determined in the form of infinite terms series and by intercepting series an approximate solution is obtained.Numerical examples are employed to show the effectiveness of the proposed approach.     

Introducing a Green–Volterra series formalism to solve weakly nonlinear boundary problems: Application to Kirchhoff's string

This paper introduces a formalism which extends that of “Green's function” and that of “the Volterra series”. These formalisms are typically used to solve, respectively, linear inhomogeneous space–time differential equations in physics and weakly nonlinear time-differential input-to-output systems in automatic control. While Green's function is a space–time integral kernel which fully characterizes a linear problem, the Volterra series expansions involve a sequence of multi-variate time integral kernels (of convolution type for time-invariant systems). The extension proposed here consists in combining the two approaches, by introducing a series expansion based on multi-variate space–time integral kernels. This series allows the representation of the space–time solution of weakly nonlinear boundary problems excited by an “input” which depends on space and time.     

Meshless analysis of three-dimensional steady-state heat conduction problems

Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attracted much attention and research recently. These problems are independent of time and involve only space coordinates, as in Poisson’s equation or the Laplace equation with Dirichlet, Neuman, or mixed conditions. When the problems are too complex, it is difficult to find an analytical solution, the only choice left is an approximate numerical solution. This paper deals with the numerical solution of three-dimensional steady-state heat conduction problems using the meshless reproducing kernel particle method (RKPM). A variational method is used to obtain the discrete equations. The essential boundary conditions are enforced by the penalty method. The effectiveness of RKPM for three-dimensional steady-state heat conduction problems is investigated by two numerical examples.     

瞬态热传导问题的一阶对称SPH方法模拟

为提高传统光滑粒子动力学(smoothed particle hydrodynamics, SPH)方法模拟瞬态热传导问题的精度和稳定性,本文提出了一种一阶对称光滑粒子动力学(first order symmetric SPH, FO-SSPH)方法.该方法将具有二阶热传导方程分解成两个一阶偏微分方程,然后基于梯度离散和Taylor级数展开思想,对一阶核梯度形式进行修正,并将得到的局部矩阵对称化.数值结果表明:与传统SPH方法相比,FO-SSPH方法精度高、数值稳定性好; 该方法能较准确地直接施加混合边值 关键词： 瞬态热传导 光滑粒子动力学 非线性     

Shallow rectangular TLD for structural control implementation

A simple and practical model for the application of shallow rectangular tuned liquid damper (TLD) in structural vibration control is presented in this paper. The dynamic properties of shallow liquid in rectangular containers subjected to forced horizontal oscillation are analysed directly from the continuity and momentum equations of fluids. Following some practical assumptions, the nonlinear partial differential equations describing the wave movement of shallow liquid in rectangular containers are established and a numerical procedure for the solutions of these equations is proposed based on the finite element method. The formula for determining the control force provided by a shallow rectangular TLD is presented. The advantage of the proposed approach for the modeling of shallow rectangular TLD is that it simplifies a three-dimensional problem into a one-dimensional problem and therefore reduces the computation efforts significantly. The whole process forms a solid foundation and provides a simplified procedure for the design and analysis of shallow rectangular TLD.     

基于小波核偏最小二乘回归方法的混沌系统建模研究

基于核学习的强大非线性映射能力,结合用于回归建模的线性偏最小二乘(PLS)算法,提出一种小波核偏最小二乘(WKPLS)回归方法. 该方法基于支持向量机使用的经典核函数技巧,将输入映射到高维非线性的特征空间,在特征空间中,构造线性的PLS回归模型. PLS方法利用输入与输出变量之间的协方差信息提取潜在特征,而可允许的小波核函数具有近似正交以及适用于信号局部分析的特性. 因此,结合它们优点的WKPLS方法显示了更好的非线性建模性能. 将WKPLS方法应用在非线性混沌动力系统建模上,并与基于高斯核的核偏最小二乘 关键词： 小波核 偏最小二乘回归 混沌系统 建模     

A kernel-free boundary integral method for elliptic boundary value problems

This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green’s functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green’s functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong.     

Using Numerical Solutions in the Fourier Method for the Fokker–Planck–Kolmogorov Equation in the Analysis of First-Order Stochastic Nonlinear Systems for Nonlinear Detectors

We present a method for analyzing the characteristics of nonlinear detectors using the algorithms of first-order nonlinear differential equations. This method is based on numerical solutions of the Fokker–Planck–Kolmogorov (FPK) equations in the form of series of functions over Hermite–Chebyshev polynomials for both nonlinear systems and their linear counterparts. The results of the solutions for the linear case are extended to nonlinear systems in a recurrent way.