Multi-modal method for chatter stability prediction and control in milling of thin-walled workpiece

Chatter stability in milling can be predicted by analytical methods or numerical methods. The system should be considered as multi-modal in milling of thin-walled workpiece. This paper proposes a numerical difference method based on Adams-Bashforth scheme. Moreover, multi-modal scheme of numerical methods is proposed. Analytical methods and numerical methods are verified by performing a series of milling trials. The experimental results are consistent with the predicted critical stability boundaries. Moreover, a new method for analyzing the computational time of analytical methods and numerical methods, which is based on the time complexity modeling, is presented. Computational time can be expressed as exact mathematical expression. By using the expression, the rate of increase of computational time can be derived.     

求解流函数涡度方程的特征线—差分方法

本文讨论了用特征线方法与有限差分方法相结合的数值方法(特征线—差分方法)求解流函数涡度形式的Navier-Stokes方程的问题.证明了该方法的收敛性,给出了数值例子.     

基于人工鱼群算法求任意函数的数值积分

将不等距离分割方法与人工鱼群算法相结合,提出一种基于人工鱼群算法求任意函数数值积分的方法,该方法除能计算通常意义下任意函数的定积分外,还能计算奇异函数积分、振荡函数积分以及原函数不易求得的被积函数的积分.最后给出几个数值积分算例,并与传统数值积分方法作了比较,仿真结果分析表明,该算法十分有效,能够快速有效地获得任意函数的数值积分值.     

用有限差分—Chebyshev Tau方法求Poisson方程的准确解

提出了一种求解二维Poisson方程的新方法－有限差分－Chebyshev Tau方法，并给出了一些有关的数值结果，结果表明，这一方法是令人满意的，且与其它方法相容。     

一类变形的牛顿法求解矩阵平方根

提出一些改进的方法来计算矩阵A的平方根,也就是应用一些牛顿法的变形来解决二次矩阵方程．研究表明,改进的方法比牛顿算法和一些已有的牛顿算法的变形效果要好．通过迭代方法,举出一些数值例子说明改进的方法的性能．     

Variable-Elliptic-Vortex Method for Incompressible Flow Simulation

A variable-elliptic-vortex method, which is a generalization of the elliptic-vortex method proposed by the author in [1], is presented for the numerical simulation of incompressible flows. The most attractive feature of the new method is that the numerical vortex blobs used in this model like actual vortex blobs can be translated, rotated and deformed in elliptic shape. The new method provides a more reasonable and more accurate approach for flow simulation than the fixed-vortex methods. Numerical examples are presented to demonstrate the performance of the new method.     

多体系统动力学微分／代数方程组的一类新的数值分析方法

本文讨论了多体系统动力学微分／代数混合方程组的数值离散问题．首先把参数t并入广义坐标讨论，简化了方程组及其隐含条件的结构，并将其化为指标1的方程组．然后利用方程组的特殊结构，引入一种局部离散技巧并构造了相应的算法．算法结构紧凑，易于编程，具有较高的计算效率和良好的数值性态，且其形式适合于各种数值积分方法的的实施．文末给出了具体算例．     

Contributions to the development of differential systems exactly solved by multistep finite-difference schemes

The motivation underlying this contribution has been to complete some of the topics concerning exact schemes for numerically solving ordinary differential equations. A procedure for obtaining differential systems exactly solved by a given finite-difference method is described. Examples illustrating the application of the procedure for obtaining first-order, second-order and systems of differential equations exactly solved by different numerical methods are given. Among the numerical methods considered there are the Trapezoidal Rule, the two-step Adams-Bashforth method and the Numerov method. Some numerical examples are presented to provide evidence that the procedure works properly.     

Numerical Methods for the Maxnear Criterion of Multiple-Sets Canonical Analysis

This paper deals with numerical methods for the Maxnear criterion of multiple-sets canonical analysis. Optimality conditions are derived. Upper and lower bounds of the optimal objective function value are presented. Two iterative methods are proposed. One is an alternating variable method, and the other called Gauss-Seidel method is an inexact version of the alternating variable method. Convergence of these methods are analyzed. A starting point strategy is suggested for both methods. Numerical results are presented to demonstrate the efficiency of these methods and the starting point strategy.     

Practical finite‐analytic method for solving differential equations by compact numerical schemes

Methodology for development of compact numerical schemes by the practical finite‐analytic method (PFAM) is presented for spatial and/or temporal solution of differential equations. The advantage and accuracy of this approach over the conventional numerical methods are demonstrated. In contrast to the tedious discretization schemes resulting from the original finite‐analytic solution methods, such as based on the separation of variables and Laplace transformation, the practical finite‐analytical method is proven to yield simple and convenient discretization schemes. This is accomplished by a special universal determinant construction procedure using the general multi‐variate power series solutions obtained directly from differential equations. This method allows for direct incorporation of the boundary conditions into the numerical discretization scheme in a consistent manner without requiring the use of artificial fixing methods and fictitious points, and yields effective numerical schemes which are operationally similar to the finite‐difference schemes. Consequently, the methods developed for numerical solution of the algebraic equations resulting from the finite‐difference schemes can be readily facilitated. Several applications are presented demonstrating the effect of the computational molecule, grid spacing, and boundary condition treatment on the numerical accuracy. The quality of the numerical solutions generated by the PFAM is shown to approach to the exact analytical solution at optimum grid spacing. It is concluded that the PFAM offers great potential for development of robust numerical schemes. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A scheme for stable numerical differentiation

A method for stable numerical differentiation of noisy data is proposed. The method requires solving a Volterra integral equation of the second kind. This equation is solved analytically. In the examples considered its solution is computed analytically. Some numerical results of its application are presented. These examples show that the proposed method for stable numerical differentiation is numerically more efficient than some other methods, in particular, than variational regularization.     

约束优化问题的修正共轭梯度投影算法

对闭凸集约束的非线性规划问题构造了一个修正共轭梯度投影下降算法,在去掉迭代点列有界的条件下,分析了算法的全局收敛性.新算法与共轭梯度参数结合,给出了三类结合共轭梯度参数的修正共轭梯度投影算法.数值例子表明算法是有效的.     

Zonal modeling of radiative heat transfer in industrial furnaces using simplified model for exchange area calculation

The zonal analysis of industrial furnaces is considered with three-dimensional radiative heat transfer, incorporated with the mathematical zone method. In this method exchange areas are determined by simplified numerical integration in three dimensions for surface-surface, surface-gas and gas–gas zones for absorbing and emitting media. By focusing on new strategies to overcome the drawbacks of evaluating direct exchange areas, it is shown that the zone method is an effective numerical method for modeling three-dimensional thermal performance of gas-filled enclosures. Also the developed method for evaluating of exchange area is presented and compared with other methods in both sides of CPU time and accuracy. The method can decrease about 70% in error of calculation of some exchange areas as compared with the other numerical methods.     

Implicit–explicit Runge–Kutta methods for financial derivatives pricing models

Implicit–explicit Runge–Kutta methods are investigated for application to financial derivatives pricing models in the partial differential equations approach. The methods are showed to be an alternative to other existing procedures for the numerical valuation of American type contracts. We follow the method of lines in order to have a numerical method that can be used with a variety of state variable discretizations including finite elements, finite differences and finite volume methods. Some numerical experiments are presented.     

一类块隐式混合单步并行计算方法

本文讨论了一类解常微分方程初值问题的块隐式混合单步并行算法，这种算法的块数为K，精度阶为2d＋2，可在S台处理机上进行并行计算，其中K＝S·d．本文讨论了方法的一般性质，给出了方法的稳定性定理，最后给出了一个数值例子．     

Discrete-analytic schemes for solving an inverse coefficient heat conduction problem in a layered medium with gradient methods

A method for constructing numerical schemes for an inverse coefficient heat conduction problem with boundary measurement data and piecewise-constant coefficients is considered. Some numerical schemes for a gradient optimization algorithm to solve the inverse problem are presented. The method is based on locally-adjoint problems in combination with approximation methods in Hilbert spaces.     

A method of local convex majorants for solving variational-like inequalities

A numerical method based on convex approximations that locally majorize a gap function is proposed for solving a variational-like inequality. The algorithm is theoretically validated and the results of comparison of its numerical efficiency to that of the conventional methods are presented.     

Two-level defect-correction Oseen iterative stabilized finite element methods for the stationary Navier–Stokes equations

Two-level defect-correction Oseen iterative stabilized finite element methods for the stationary Navier–Stokes equations based on local Gauss integration are considered in this paper. The methods combine the defect-correction method and the two-level strategy with the locally stabilized method. Moreover, the stability and convergence of the presented methods are deduced. Finally, numerical tests confirm the theoretical results of the presented methods.     

解带线性或非线性约束最优化问题的三项记忆梯度Rosen投影算法

利用Rosen投影矩阵，建立求解带线性或非线性不等式约束优化问题的三项记忆梯度Rosen投影下降算法，并证明了算法的收敛性．同时给出了结合FR，PR，HS共轭梯度参数的三项记忆梯度Rosen投影算法，从而将经典的共轭梯度法推广用于求解约束规划问题．数值例子表明算法是有效的。     

A numerical approach for solving an extended Fisher-Kolomogrov-Petrovskii-Piskunov equation

In the present paper a numerical method, based on finite differences and spline collocation, is presented for the numerical solution of a generalized Fisher integro-differential equation. A composite weighted trapezoidal rule is manipulated to handle the numerical integrations which results in a closed-form difference scheme. A number of test examples are solved to assess the accuracy of the method. The numerical solutions obtained, indicate that the approach is reliable and yields results compatible with the exact solutions and consistent with other existing numerical methods. Convergence and stability of the scheme have also been discussed.