Meshless solution of two-dimensional incompressible flow problems using the radial basis integral equation method

The two-dimensional incompressible fluid flow problems governed by the velocity–vorticity formulation of the Navier–Stokes equations were solved using the radial basis integral (RBIE) equation method. The RBIE is a meshless method based on the multi-domain boundary element method with overlapping subdomains. It solves at each node for the potential and its spatial derivatives. This feature of the RBIE is advantageous in solving the velocity–vorticity formulation of the Navier–Stokes equations since the calculated velocity gradients can be used to compute the vorticity that is prescribed as a boundary condition to the vorticity transport equation. The accuracy of the numerical solution was examined by solving the test problem with known analytical solution. Two benchmark problems, i.e. the lid driven cavity flow and the thermally driven cavity flow were also solved. The numerical results obtained using the RBIE showed very good agreement with the benchmark solutions.     

基于三角小波的弹性薄板高精度分析方法

采用同时具有三角函数良好逼近特性和小波多分辨率与局部特性的Hermite插值型三角小波,基于二维张量积三角小波,推导了求解各种不同边界条件下的矩形弹性薄板的弯曲、振动和屈曲问题的统一列式,同时给出了两种提高计算精度的方法一升阶法和多分辨率法。数值算例表明,三角小波法求解弹性薄板的弯曲、振动和屈曲问题时,能方便地处理各类边界条件,计算效果良好;自振特性分析更具优势,升阶法和多分辨率法能有效地提高分析精度。     

A finite element formulation satisfying the discrete geometric conservation law based on averaged Jacobians

In this article, a new methodology for developing discrete geometric conservation law (DGCL) compliant formulations is presented. It is carried out in the context of the finite element method for general advective–diffusive systems on moving domains using an ALE scheme. There is an extensive literature about the impact of DGCL compliance on the stability and precision of time integration methods. In those articles, it has been proved that satisfying the DGCL is a necessary and sufficient condition for any ALE scheme to maintain on moving grids the nonlinear stability properties of its fixed‐grid counterpart. However, only a few works proposed a methodology for obtaining a compliant scheme. In this work, a DGCL compliant scheme based on an averaged ALE Jacobians formulation is obtained. This new formulation is applied to the θ family of time integration methods. In addition, an extension to the three‐point backward difference formula is given. With the aim to validate the averaged ALE Jacobians formulation, a set of numerical tests are performed. These tests include 2D and 3D diffusion problems with different mesh movements and the 2D compressible Navier–Stokes equations. Copyright © 2011 John Wiley & Sons, Ltd.     

弹塑性摩擦接触问题形状设计灵敏度分析

提出了一种计算二维有限变形弹塑性摩擦接触问题形状设计灵敏度的算法. 采用主动集策略和mortar方法处理接触边线上的约束条件. 在mortar接触边线的切线和法线方向上采用相同的名义罚函数,提出基于名义罚函数的移动摩擦锥算法来正则化接触约束条件,发展了一种新的二维多体有限变形摩擦接触算法. 在此基础上, 通过将离散形式的摩擦接触问题控制方程对形状设计变量微分,得到了该路径相关问题的直接微分法解析设计灵敏度计算格式, 其节点位移灵敏度方程在每个增量步不用迭代、直接求解. 与国际上现有的二维多体有限变形摩擦接触问题的解析设计灵敏度算法相比,本算法不需分解为法向和切向推导,表达式较简洁,便于编程实现. 数值算例验证了算法的精度和有效性.      

Dynamic analysis and regulation of the flexible pipe conveying fluid with a hard-magnetic soft segment

Applied Mathematics and Mechanics - The recently developed hard-magnetic soft (HMS) materials can play a significant role in the actuation and control of medical devices, soft robots, flexible...     

On the numerical solution of tracked vehicle dynamic equations

In this investigation, the solution of the nonlinear dynamic equations of the multibody tracked vehicle systems are obtained using different procedures. In the first technique, which is based on the augmented formulation that employes the absolute Cartesian coordinates and Lagrange multipliers, the generalized coordinate partitioning of the constraint Jacobian matrix is used to determine the independent coordinates and the associated independent differential equations. An iterative Newton-Raphson algorithm is used to solve the nonlinear constraint equations for the dependent variables. The numerical problems encountered when one set of independent coordinates is used during the simulation of large scale tracked vehicle systems are demonstrated and their relationship to the track dynamics is discussed. The second approach employed in this investigation is the velocity transformation technique. One of the versions of this technique is discussed in this paper and the numerical problems that arise from the use of inconsistent system of kinematic equations are reported. In the velocity transformation technique, the tracked vehicle system is assumed to consist of two kinematically decoupled subsystems; the first subsystem consists of the chassis, the rollers, the sprocket and the idler, while the second subsystem consists of the track which is represented as a closed kinematic chain that consists of rigid links connected by revolute joints. It is demonstrated that the use of one set of recursive equations leads to numerical difficulties because of the change in the track configuration. Singular configurations can be avoided by repeated changes in the recursive equations. The sensitivity of the predictor-corrector multistep numerical integration schemes to the method of formulating the state equations is demonstrated. The numerical results presented in this investigation are obtained using a planner tracked vehicle model that consists of fifty four rigid bodies.     

Interfacial roughness and angle dependence of giant magnetoresistance in magnetic superlattices

Using the two-point conductivity formula, we numerically evaluate the giant magnetoresistance (GMR) in magnetic superlattices with currents in the plane of the layers (CIP), from which the effect of the interfacial roughness and magnetization configuration on the GMR is studied. With increasing interfacial roughness, the maximal GMR ratio is found to first increase and then decrease, exhibiting a peak at an optimum strength of interfacial roughness. For systems composed of relatively thick layers, the GMR is approximately proportional to ,where is the angle between the magnetizations in two successive ferromagnetic layers, but noticeable departures from this dependence are found when the layers become sufficiently thin. Received 21 September 1998 and Received in final form 22 December 1998     

Three-dimensional incompressible Navier–Stokes equations on non-orthogonal staggered grids using the velocity–vorticity formulation

This paper is concerned with the numerical resolution of the incompressible Navier–Stokes equations in the velocity–vorticity form on non-orthogonal structured grids. The discretization is performed in such a way, that the discrete operators mimic the properties of the continuous ones. This allows the discrete equivalence between the primitive and velocity–vorticity formulations to be proved. This last formulation can thus be seen as a particular technique for solving the primitive equations. The difficulty associated with non-simply connected computational domains and with the implementation of the boundary conditions are discussed. One of the main drawback of the velocity–vorticity formulation, relative to the additional computational work required for solving the additional unknowns, is alleviated. Two- and three-dimensional numerical test cases validate the proposed method. © 1998 John Wiley & Sons, Ltd.     

Mixed boundary elements for laminar flows

This paper presents a mixed boundary element formulation of the boundary domain integral method (BDIM) for solving diffusion–convective transport problems. The basic idea of mixed elements is the use of a continuous interpolation polynomial for conservative field function approximation and a discontinuous interpolation polynomial for its normal derivative along the boundary element. In this way, the advantages of continuous field function approximation are retained and its conservation is preserved while the normal flux values are approximated by interpolation nodal points with a uniquely defined normal direction. Due to the use of mixed boundary elements, the final discretized matrix system is overdetermined and a special solver based on the least squares method is applied. Driven cavity, natural and forced convection in a closed cavity are studied. Driven cavity results at Re=100, 400 and 1000 agree better with the benchmark solution than Finite Element Method or Finite Volume Method results for the same grid density with 21×21 degrees of freedom. The average Nusselt number values for natural convection 10³≤Ra≤10⁶ agree better than 0.1% with benchmark solutions for maximal calculated grid densities 61×61 degrees of freedom. Copyright © 1999 John Wiley & Sons, Ltd.