不同网格扭曲率下压力修正全隐算法——IDEAL求解性能研究

2008年,本文作者和陶文铨等提出了一种用于速度和压力耦合求解的高效稳定压力修正全隐算法IDEAL,该算法通过在每个迭代层次上对压力方程进行两次内迭代计算,完全克服了SIMPLE算法的两个假设,充分满足了速度和压力之间的耦合,从而大大提高了计算的收敛性和健壮性.为了进一步实现IDEAL算法的推广应用,本文基于三维倾斜方腔顶盖驱动流动,研究了IDEAL算法在不同网格扭曲率下的求解特性.研究发现,在不同网格扭曲率下,IDEAL算法的健壮性和收敛性均优于SIMPLE算法,特别在高网格扭曲率情况下,IDEAL算法求解性能更加优于SIMPLE算法.在不同网格扭曲率下,IDEAL算法健壮性保持不变,几乎可以在任意速度亚松弛因子下获得收敛的解,同时IDEAL算法最短计算耗时较SIMPLE算法减少了56%~89%,验证了IDEAL算法的优越性.     

不同网格扭曲率下压力修正全隐算法——IDEAL求解性能研究

2008年,本文作者和陶文铨等提出了一种用于速度和压力耦合求解的高效稳定压力修正全隐算法IDEAL,该算法通过在每个迭代层次上对压力方程进行两次内迭代计算,完全克服了SIMPLE算法的两个假设,充分满足了速度和压力之间的耦合,从而大大提高了计算的收敛性和健壮性。为了进一步实现IDEAL算法的推广应用,本文基于三维倾斜方腔顶盖驱动流动,研究了IDEAL算法在不同网格扭曲率下的求解特性。研究发现,在不同网格扭曲率下,IDEAL算法的健壮性和收敛性均优于SIMPLE算法,特别在高网格扭曲率情况下,IDEAL算法求解性能更加优于SIMPLE算法。在不同网格扭曲率下,IDEAL算法健壮性保持不变,几乎可以在任意速度亚松弛因子下获得收敛的解,同时IDEAL算法最短计算耗时较SIMPLE算法减少了56%~89%,验证了IDEAL算法的优越性。     

基于半光滑牛顿法的润滑液膜有限元空化算法

针对润滑液膜中空化问题,引入Fischer-Burmeister函数,提出一种求解满足质量守恒雷诺方程的半光滑牛顿迭代算法.该算法将空化问题的非线性互补关系转化为等式约束方程,避免了迭代计算中的不等式约束识别问题.算法可将空化约束方程与雷诺方程、力平衡方程、变形方程等同时纳入牛顿迭代方程组,有效解决了传统松弛迭代算法需要多重嵌套循环带来的效率低下问题及压力与膜厚的强耦合性带来的收敛困难问题.计算实例表明,该算法计算效率高、收敛性好,且易应用于弹流润滑分析中,在滑动轴承和机械端面密封等多种物理模型下均有良好的适用性.     

求解接触问题的一种新的实验误差法

提出了一种带松弛因子的UZAW算法求解实验误差法中给定状态下的位移和接触力满足的等式方程，并证明了该算法是R超线性收敛的。整个区域被划分为多个子区域，不同子区域位移场的求解是独立的。还提出了一种带参数的以不完全因子分解为基础的预条件子共轭梯度法求解不同子区域位移场，该算法在块体规模较大时更加有效。     

基于贴体网格的VOF方法数模流场研究

提出了一种基于VOF方法的模拟具有复杂边界形状结构物附近流场的新算法,BFC—SIMPLE—VOF算法。采用坐标变换方法实现了任意复杂区域的结构化网格划分,在贴体网格下对二维不可压缩粘性流体的控制方程进行了离散。提出了基于交错网格的修正SIMPLE算法来迭代求解压力一速度场,修正了贴体坐标下的界面跟踪方法（VOF方法）...     

基于非结构化同位网格的SIMPLE算法

通过基于非结构化网格的有限体积法对二维稳态Navier—Stokes方程进行了数值求解。其中对流项采用延迟修正的二阶格式进行离散；扩散项的离散采用二阶中心差分格式；对于压力-速度耦合利用SIMPLE算法进行处理；计算节点的布置采用同位网格技术，界面流速通过动量插值确定。本文对方腔驱动流、倾斜腔驱动流和圆柱外部绕流问题进行了计算，讨论了非结构化同位网格有限体积法在实现SIMPLE算法时，迭代次数与欠松弛系数的关系、不同网格情况的收敛性、同结构化网格的对比以及流场尾迹结构。通过和以往结果比较可知，本文的方法是准确和可信的。     

结构可靠度分析FORM迭代算法的混沌控制

利用混沌控制原理对FORM收敛失败进行控制. 理清了全局性和局部性两类混沌反馈 控制各种方法的内在联系，说明稳定转换法和自适应调节法属于全局混沌反馈控制 方法，自适应调节法可视为稳定转换法的特例. 参 数调节混合法不过是松弛牛顿法的另一种表达形式，它们都属于局部混沌反馈控制方法. 阐 明了混沌反馈控制表达式与工程力学收敛控制迭代算法的对应关系. 也揭示了这些迭代算法 收敛控制措施的功效和局限性. 提出了一个以稳定转换法为主联合松弛牛顿法的混 沌反馈控制方法，对可靠度分析FORM迭代算法实现了周期振荡、分岔和混沌控制.     

自适应SVD-UKF算法及在组合导航的应用

提出一种新的自适应奇异值分解Unscented卡尔曼滤波(UKF)算法.该算法利用自适应因子平衡动力学模型信息与观测信息的权比,控制动力学模型误差对导航参数解的影响.用奇异值分解阵(SVD)的迭代计算代替协方差矩阵的迭代变换,提高了协方差矩阵的数值稳定性.将新算法应用于组合导航系统进行计算仿真,结果证明,新算法具有良好的鲁棒性,能有效改善滤波性能,提高组合导航系统的精度.     

基于多项式迭代的等效旋转矢量微分方程精确数值算法

提出了一种基于多项式迭代的求解等效旋转矢量微分方程(Bortz方程)的新算法,以角速度多项式作为输入,利用泰勒级数将Bortz方程中的余切函数展开成多项式形式,将不可交换误差补偿中的叉乘和求模运算变换为多项式的卷积运算,通过迭代求取等效旋转矢量的多项式精确解,有效地解决了捷联惯导系统在大角度机动环境下的姿态精确确定问题。此外,还对新算法的数值运算量和迭代收敛精度进行了分析。最后,进行了半锥角90°、频率1 Hz、时长1 s的大幅值圆锥运动姿态更新仿真实验,与传统圆锥优化算法的5'姿态误差相比,新算法的误差仅为10~(-7)',具有明显的精度优势。     

基于噪声和残缺投影图像的重建算法研究

滤波反投影法和迭代法是SR-CT重建算法中的两大类算法,松弛迭代法是迭代法中的一种;本文主要研究滤波反投影方法和松弛迭代法.为了研究两种方法的影响因素,讨论了滤波反投影法中滤波函数的选取,分析了松弛迭代法中重要参数,如松弛因子、迭代次数对重建效果的影响,给出了重建效果和这些参数的影响关系;对比了两种方法在处理噪声和残缺...     

On the improvement of the SIMPLE-like method for flows with complex geometry

For strongly nonorthogonal grids, the convergence difficulty is often encountered in the SIMPLE method unless very low relaxation factors are adopted. This is because that the cross derivatives of pressure corrections are usually neglected in deriving the pressure-correction equation. In this study, an approximate treatment is proposed to deal with the cross derivatives in the pressure-correction equation. Computations are performed for various laminar flows in a nonorthogonal cavity and the turbulent flows through the compressor cascades. It is found that the convergence difficulty encountered in the SIMPLE method for strongly nonorthogonal grids can be overcome when the SIMPLEC method is employed. In addition, when the present approximate treatment for the cross derivatives is adopted in the SIMPLEC method, the convergence rate can be speeded up in comparison with the original SIMPLEC method. In some cases, the CPU time can be decreased by more than 60%. The computational results are also compared with the numerical results of others and the available experimental data. Received on 30 November 1998     

Nonlinear adaptive task space control for a 2-DOF redundantly actuated parallel manipulator

A nonlinear adaptive (NA) controller in the task space is developed for the trajectory tracking of a 2-DOF redundantly actuated parallel manipulator. The dynamic model with nonlinear friction is established in the task space for the parallel manipulator, and the linear parameterization expression of the dynamic model is formulated. Based on the dynamic model, a new control law including adaptive dynamics compensation, adaptive friction compensation and error elimination items is designed. After defining a quadratic performance index, the parameter update law is derived with the gradient descent algorithm. The stability of the parallel manipulator system is proved by the Lyapunov theorem, and the convergence of the tracking error and the error rate is proved by the Barbalat’s lemma. The NA controller is implemented in the trajectory tracking experiments of an actual 2-DOF redundantly actuated parallel manipulator, and the experiment results are compared with the APD controller.     

A novel implicit numerical treatment for non‐linear turbulence models using high and low Reynolds number formulations

Non‐linear turbulence models can be seen as an improvement of the classical eddy‐viscosity concept due to their better capacity to simulate characteristics of important flows. However, application of non‐linear models demand robustness of the numerical method applied, requiring a stable discretization scheme for convergence of all variables involved. Usually, non‐linear terms are handled in an explicit manner leading to possible numerical instabilities. Thus, the present work shows the steps taken to adapt a general non‐linear constitutive equation using a new semi‐implicit numerical treatment for the non‐linear diffusion terms. The objective is to increase the degree of implicitness of the solution algorithm to enhance convergence characteristics. Flow over a backward‐facing step was computed using the control volume method applied to a boundary‐fitted coordinate system. The SIMPLE algorithm was used to relax the algebraic equations. Classical wall function and a low Reynolds number model were employed to describe the flow near the wall. The results showed that for certain combination of relaxation parameters, the semi‐implicit treatment proposed here was the sole successful treatment in order to achieve solution convergence. Also, application of the implicit method described here shows that the stability of the solution either increases (high Reynolds with non‐orthogonal mesh) or preserves the same (low Reynolds number applications). Additional advantages of the procedure proposed here lie in the possibility of testing different non‐linear expressions if one considers the enhanced robustness and stability obtained for the entire numerical algorithm. Copyright © 2010 John Wiley & Sons, Ltd.     

一种快速逼近Fourier级数和的实用算子

采用逐次平均方法,构造一种快速逼近Ｆｏｕｒｉｅｒ级数和实用算子（逐次平均算子）,用这种算子来取代Ｆｏｕｒｉｅｒ级数无穷项求和,具有误差小,收敛速度快的功能。志Ｆｅｊｅｒ算子比,它的逼近性要好得多,与Ｋｏｒｂｉｎ算子比,它更易于实用,理论证明其一般性,计算著名的Ｓａｆｆｍａｎ－Ｔａｙｌｏｒ粘性指进问题时应用逐次平均算子,既不产生大的截断误差,又能避免高频振荡而使计算过程稳定,得到很好的结果。     

Accelerated convergence of the numerical simulation of incompressible flow in general curvilinear co‐ordinates by discretizations on the double‐staggered grids

The convergence rate of a methodology for solving incompressible flow in general curvilinear co‐ordinates is analyzed. Double‐staggered grids (DSGs), each defined by the same boundaries as the physical domain, are used for discretization. Both grids are MAC quadrilateral meshes with scalar variables (pressure, temperature, etc.) arranged at the center and the Cartesian velocity components at the middle of the sides of the mesh cells. The problem was checked against benchmark solutions of natural convection in a squeezed cavity, heat transfer in concentric horizontal cylindrical annuli, and a hot cylinder in a duct. Poisson's pressure‐correction equations that arise from the SIMPLE‐like procedure are solved by several methods: successive overrelaxation, symmetric overrelaxation, modified incomplete factorization preconditioner, conjugate gradient (CG), and CG with preconditioner. A genetic algorithm was developed to solve problems of numerical optimization of SIMPLE‐like calculation time in a space of iteration numbers and relaxation parameters. The application provides a means of making an unbiased comparison between the DSGs method and the widely used interpolation method. Furthermore, the convergence rate was demonstrated by application to the calculation of natural convection heat transfer in concentric horizontal cylindrical annuli. Calculation times when DSGs were used were 2–10 times shorter than those achieved by interpolation. With the DSGs method, calculation time increases slightly with increasing non‐orthogonality of the grids, whereas an interpolation method calls for very small iteration parameters that lead to unacceptable calculation times. Copyright © 2007 John Wiley & Sons, Ltd.     

A new method for accelerating the rate of convergence of the simple-like algorithm

A pressure correction formula is proposed for the SIMPLE-like algorithm in order to improve the rate of the convergence when solving laminar Navier–Stokes equations when there is rapidly varying pressure. Based on global mass conservation, a line average pressure correction is derived by integration of the momentum equation for approximate one-dimensional flow. The use of this formula with the SIMPLE-like algorithm can rapidly build up the pressure distribution in the region where the pressure undergoes a very large change, which normally causes the rate of convergence of the SIMPLE or the SIMPLEC schemes to be slow. In order to illustrate the technique, the performances of SIMPLE and of SIMPLEC with the average pressur correction are investigated for axisymmetric flow past and through a sampler. A comparison of these two techniques shows that the average pressure correction proposed in this paper significantly accelerates the rate of convergence.     

Dynamic model based nonlinear tracking control of a planar parallel manipulator

Based on the dynamic model, a novel nonlinear tracking controller is developed to overcome the nonlinear dynamics and friction of a planar parallel manipulator. The dynamic model is formulated in the active joint space, and the active joint friction is described with the Coulomb + viscous friction model. A nonlinear tracking controller is designed to eliminate the tracking error by using the power function. The nonlinear tracking controller is proven to guarantee asymptotic convergence to zero of both the tracking error and error rate with the Barbalat’s lemma. The trajectory tracking experiment of the proposed controller is implemented on an actual five-bar planar parallel manipulator both at the low-speed and high-speed motion. Moreover, the control performances of the proposed controller are compared with the results of the augmented PD (APD) controller.     

A numerical method for accelerating the rate of convergence of the SIMPLE-like algorithm for flow through a thin filter

The control volume, finite difference method and the k-? tubulence model are employed in a numerical simulation of the turbulent fluid flow both outside and inside a blunt cylindrical sampler which houses a paper filter in its chamber. The presence of a paper filter, which has a very large resistance, results in a large pressure drop across the filter and this causes difficulties in making the SIMPLE or the SIMPLEC scheme converge. In order to improve the rate of convergence of the SIMPLE-like algorithm when the resistance of the filter is very large, an average pressure correction formula is proposed. Based on global mass conservation, a line average pressure correction for the paper filter is derived using a modified Darcy law for a porous medium. A combination of this formula and the SIMPLE-like algorithm can rapidly build up the pressure drop across the filter and hence dramatically improve the rate of convergence of the iterative scheme. Comparisons of the convergence histories and the numerical results for the fluid flow when using SIMPLE and SIMPLEC with the average pressure correction method show that the average pressure correction method for dealing with the paper filter significantly accelerates the rate of convergence of the iterative scheme.     

Vibrational self-alignment of a rigid object exploiting friction

In this paper, one-dimensional self-alignment of a rigid object via stick-slip vibrations is studied. The object is situated on a table, which has a prescribed periodic motion. Friction is exploited as the mechanism to move the object in a desired direction and to stop and self-align the mass at a desired end position with the smallest possible positioning error. In the modeling and analysis of the system, theory of discontinuous dynamical systems is used. Analytic solutions can be derived for a model based on Coulomb friction and an intuitively chosen table acceleration profile, which allows for a classification of different possible types of motion. Local stability and convergence is proven for the solutions of the system, if a constant Coulomb friction coefficient is used. Next, near the desired end position, the Coulomb friction coefficient is increased (e.g. by changing the roughness of the table surface) in order to stop the object. In the transition region from low friction to high friction coefficient, it is shown that, under certain conditions, accumulation of the object to a unique end position occurs. This behavior can be studied analytically and a mapping is given for subsequent stick positions.     

Effects of under‐relaxation factors on turbulent flow simulations

Under‐relaxation factors are significant parameters affecting the convergence of a numerical scheme. Some earlier work has been done to optimize these parameters, but this was restricted to special flow domains, and the range of changes for under‐relaxation factors and convective algorithms are limited. In this paper, the effects of changing under‐relaxation factors for different variables, different convective schemes and grid sizes on the convergence of the numerical solution of three 2D turbulent flow situations are studied. These three flows are duct flow, trench flow and inclined free falling jet flow. Copyright © 2003 John Wiley & Sons, Ltd.