非结构动网格在三维可动边界问题中的应用

研究用于非结构动网格的弹簧近似方法，采用顶点弹簧描述，导出并讨论了弹簧倔强系数的取值。通过引入边界修正和扭转效应修正，对标准弹簧近似方法进行了改进，转动翼型算例的结果表明，改进后的方法大大提高了网格变形能力和网格质量，应用该动网格方法耦合求解基于（Arbitrary Lagrangian-Eulerian,ALE）描述的三维Euler方程，模拟了作俯抑振动的矩形机翼绕流，计算结果与实验数据及文献计算结果十分一致，作为多个自由刚与流体耦合运动问题的简单例证，耦合刚体动力学方程，模拟了激波与双立方体的相互作用，得到了非定常流场结构，研究表明，基于弹簧近似的非结构动网格与有限体积流式流场解算器相结合，是模拟包含运动边界的非定常流动问题的有效方法。     

非结构动网格在可动边界问题中的应用研究

本文对用于非结构动网格生成的弹簧近似方法进行了研究。通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进。改进后的方法可以大大提高网格变形能力和网格质量。应用本文发展的非结构动网格生成方法并通过耦合求解基于(ArbitraryLagrangian-Eulerian ALE)描述的Euler方程,模拟了谐和振动NACA0012翼型及M6机翼的跨音速绕流,计算结果与参考文献提供的结果及实验结果吻合良好。     

Dynamic stiffness method for free vibration of an axially moving beam with generalized boundary conditions

Axially moving beams are often discussed with several classic boundary conditions, such as simply-supported ends, fixed ends, and free ends. Here, axially moving beams with generalized boundary conditions are discussed for the first time. The beam is supported by torsional springs and vertical springs at both ends. By modifying the stiffness of the springs, generalized boundaries can replace those classical boundaries.Dynamic stiffness matrices are, respectively, established for axially moving Timoshenko beams and Euler-Bernoulli(EB) beams with generalized boundaries. In order to verify the applicability of the EB model, the natural frequencies of the axially moving Timoshenko beam and EB beam are compared. Furthermore, the effects of constrained spring stiffness on the vibration frequencies of the axially moving beam are studied. Interestingly, it can be found that the critical speed of the axially moving beam does not change with the vertical spring stiffness. In addition, both the moving speed and elastic boundaries make the Timoshenko beam theory more needed. The validity of the dynamic stiffness method is demonstrated by using numerical simulation.     

Error estimate and adaptive refinement for incompressible Navier‐Stokes equations using the discrete least squares meshless method

In this paper, an adaptive refinement strategy based on a node‐moving technique is proposed and used for the efficient solution of the steady‐state incompressible Navier–Stokes equations. The value of a least squares functional of the residual of the governing differential equation and its boundary conditions at nodal points is regarded as a measure of error and used to predict the areas of poor solutions. A node‐moving technique is then used to move the nodal points to the zones of higher numerical errors. The problem is then resolved on the refined distribution of nodes for higher accuracy. A spring analogy is used for the node‐moving methodology in which nodal points are connected to their neighbors by virtual springs. The stiffness of each spring is assumed to be proportional to the errors of its two end points and its initial length. The new positions of the nodal points are found such that the spring system attains its equilibrium state. Some numerical examples are used to illustrate the ability of the proposed scheme for the adaptive solution of the steady‐state incompressible Navier–Stokes equations. The results demonstrate a considerable improvement of the results with a reasonable computational effort by using the proposed adaptive strategy. Copyright © 2011 John Wiley & Sons, Ltd.     

A 3D mesh deformation technique for irregular in‐flight ice accretion

Aircraft holding around busy airports may be requested to sustain as much as 45 min of icing before landing or being diverted to another airport. In this paper, a three‐dimensional mesh deformation scheme, based on a structural frame analogy, is proposed for the numerical simulation of ice accretion during extended exposure to adverse weather conditions. The goal is to provide an approach that is robust and efficient enough to delay or altogether avoid re‐meshing while preserving (enforcing) nearly orthogonal elements at the highly distorted ice surface. Robustness is achieved by suitably modifying the axial and torsional stiffness components of the frame elements in order to handle large and irregular grid displacements typical of in‐flight icing. Computational efficiency is obtained by applying the mesh displacement to an automatically selected small subset of the entire computational domain. The methodology is validated first in the case of deformations typical of fluid‐structure interaction problems, including wing bending, a helicopter rotor in forward flight, and the twisting of a high‐lift wing configuration. The approach is then assessed for aero‐icing on two swept wings and compared against experimental measurements where available. Copyright © 2015 John Wiley & Sons, Ltd.     

THE STABILITY OF AN AXIALLY ACCELERATING BEAM ON SIMPLE SUPPORTS WITH TORSION SPRINGS

The axially moving beams on simple supports with torsion springs are studied. The general modal functions of the axially moving beam with constant speed have been obtained from the supporting conditions. The contribution of the spring stiffness to the natural frequencies has been numerically investigated. Transverse stability is also studied for axially moving beams on simple supports with torsion springs. The method of multiple scales is applied to the partialdifferential equation governing the transverse parametric vibration. The stability boundary is derived from the solvability condition. Instability occurs if the axial speed fluctuation frequency is close to the sum of any two natural frequencies or is two fold natural frequency of the unperturbed system. It can be concluded that the spring stiffness makes both the natural frequencies and the instability regions smaller in the axial speed fluctuation frequency-amplitude plane for given mean axial speed and bending stiffness of the beam.     

A moving mesh algorithm for aero‐thermo‐mechanical modelling in turbomachinery

This paper describes the development of a mesh deformation method used for aero‐thermo‐mechanical coupling of turbo‐engine components. The method is based on the nonlinear solution of an elastic medium analogy, solved using finite element discretisation and modified to let the boundary nodes be free to slide over the deflected surfaces. This sliding technique relies on a B‐spline reconstruction of the moving boundary and increases the robustness of the method in situations where the boundary deflection field presents significant gradients or large relative motion between two distinct boundaries. The performance of the method is illustrated with the application to an interstage cavity of a turbine assembly, subjected to the deformations computed by a coupled thermo‐mechanical analysis of the engine component. Copyright © 2011 John Wiley & Sons, Ltd.     

移动柔性梁作用下简支梁的动力响应

对移动结构作用下梁的响应问题进行了推广,采用柔性梁作为移动结构模型,在考虑结构柔性和悬挂连接的前提下对系统的耦合振动进行了分析.根据一般边界条件梁建立振动方程,通过量纲一参数以及模态叠加法处理系统动力学方程.以简支边界条件为例,得到了梁响应的数值结果,对系统主要参数即移动结构频率、移动速度及连接刚度对简支梁振动的影响进行了讨论.结果表明:考虑移动体的柔性频率对简支梁的振动会产生一定的影响.     

An extended finite element method for the simulation of particulate viscoelastic flows

We present an extended finite element method (XFEM) for the direct numerical simulation of the flow of viscoelastic fluids with suspended particles. For moving particle problems, we devise a temporary arbitrary Lagrangian–Eulerian (ALE) scheme which defines the mapping of field variables at previous time levels onto the computational mesh at the current time level. In this method, a regular mesh is used for the whole computational domain including both fluid and particles. A temporary ALE mesh is constructed separately and the computational mesh is kept unchanged throughout the whole computations. Particles are moving on a fixed Eulerian mesh without any need of re-meshing. For mesh refinements around the interface, we combine XFEM with the grid deformation method, in which nodal points are redistributed close to the interface while preserving the mesh topology. Our method is verified by comparing with the results of boundary fitted mesh problems combined with the conventional ALE scheme. The proposed method shows similar accuracy compared with boundary fitted mesh problems and superior accuracy compared with the fictitious domain method. If the grid deformation method is combined with XFEM, the required computational time is reduced significantly compared to uniform mesh refinements, while providing mesh convergent solutions. We apply the proposed method to the particle migration in rotating Couette flow of a Giesekus fluid. We investigate the effect of initial particle positions, the Weissenberg number, the mobility parameter of the Giesekus model and the particle size on the particle migration. We also show two-particle interactions in confined shear flow of a viscoelastic fluid. We find three different regimes of particle motions according to initial separations of particles.     

基于背景网格变形的动态网格移动方法

基于Delaunay背景网格插值技术的动态网格生成方法无需迭代计算,效率较高。但对复杂构形大幅运动的动边界问题,尤其当边界大幅转动时,背景网格极易交叉重叠。重新生成背景网格和重新定位网格节点信息不仅费时而且会导致网格质量的严重下降。本文提出改进的基于背景网格的动态网格变形方法,通过在初始Delaunay背景网格中添加辅助点,生成一层新的背景网格和新的映射关系;采用ball-vertex弹簧法驱动新背景网格的变形,进而牵动目标网格的变形。算例表明,本文提出的动态网格变形方法对所关心区域的网格具有良好保形性,边界可转动更大角度而不会出现网格交叉重叠问题,总体上提高了动态网格更新的效率和质量。     

基于点球弹簧修匀法的混合网格变形并行算法

针对二维/三维混合网格,提出基于点球弹簧修匀法的并行网格变形算法。按特定模板将混合网格中的非三角形/四面体单元分解成三角形/四面体单元。针对每个内部节点及其相邻节点建立相应的子弹簧系统,并通过增加Ball-Vertex弹簧避免弹簧系统的塌陷问题。由于点球弹簧法在计算中逐点对网格内部节点进行计算,在计算过程中具有良好的弱耦合性质,因此有利于算法并行化。在并行化时仅需对网格进行虚拟分区操作,不必进行复杂的几何分区操作,同时避免了混合网格不同单元之间的兼容性问题。该方法适用于具有复杂外形的大规模混合网格的变形问题,能够显著提高网格变形的效率,同时具有良好的适应性。     

A weakly compressible MPS method for modeling of open‐boundary free‐surface flow

A mesh‐free particle method, based on the moving particle semi‐implicit (MPS) interaction model, has been developed for the simulation of two‐dimensional open‐boundary free‐surface flows. The incompressibility model in the original MPS has been replaced with a weakly incompressible model. The effect of this replacement on the efficiency and accuracy of the model has been investigated. The new inflow–outflow boundary conditions along with the particle recycling strategy proposed in this study extend the application of the model to open‐boundary problems. The final model is able to simulate open‐boundary free surface flow in cases of large deformation and fragmentation of free surface. The models and proposed algorithms have been validated and applied to sample problems. The results confirm the model's efficiency and accuracy. Copyright © 2009 John Wiley & Sons, Ltd.     

外激励作用下弹性支撑浅拱的非线性动力学分析

研究了外激励下两端采用转动弹簧约束的铰支浅拱在发生1:1内共振时的非线性动力学行为。通过引入基本假定和无量纲化变量得到浅拱的动力学控制方程, 将阻尼项、外荷载项和非线性项去掉后,所得线性方程及对应边界条件即可确定考虑转动弹簧影响的频率和模态, 发现转动约束取不同刚度值时系统存在模态交叉与模态转向两种内共振形式。对动力方程进行Galerkin全离散, 并采用多尺度法对内共振进行了摄动分析, 得到了极坐标和直角坐标两种形式的平均方程, 其中平均方程系数与转动弹簧刚度一一对应。最低两阶模态之间1:1内共振的数值研究结果表明: 外激励能激发内共振模态的非线性相互作用, 参数处于某一范围时系统存在周期解、准周期解和混沌解窗口, 且通过(逆)倍周期分岔方式进入混沌。     

Numerical simulation of viscous flow interaction with an elastic membrane

A numerical fluid–structure interaction model is developed for the analysis of viscous flow over elastic membrane structures. The Navier–Stokes equations are discretized on a moving body‐fitted unstructured triangular grid using the finite volume method, taking into account grid non‐orthogonality, and implementing the SIMPLE algorithm for pressure solution, power law implicit differencing and Rhie–Chow explicit mass flux interpolations. The membrane is discretized as a set of links that coincide with a subset of the fluid mesh edges. A new model is introduced to distribute local and global elastic effects to aid stability of the structure model and damping effects are also included. A pseudo‐structural approach using a balance of mesh edge spring tensions and cell internal pressures controls the motion of fluid mesh nodes based on the displacements of the membrane. Following initial validation, the model is applied to the case of a two‐dimensional membrane pinned at both ends at an angle of attack of 4° to the oncoming flow, at a Reynolds number based on the chord length of 4 × 10³. A series of tests on membranes of different elastic stiffness investigates their unsteady movements over time. The membranes of higher elastic stiffness adopt a stable equilibrium shape, while the membrane of lowest elastic stiffness demonstrates unstable interactions between its inflated shape and the resulting unsteady wake. These unstable effects are shown to be significantly magnified by the flexible nature of the membrane compared with a rigid surface of the same average shape. Copyright © 2007 John Wiley & Sons, Ltd.     

A new modification of the immersed‐boundary method for simulating flows with complex moving boundaries

In this paper, a new immersed‐boundary method for simulating flows over complex immersed, moving boundaries is presented. The flow is computed on a fixed Cartesian mesh and the solid boundaries are allowed to move freely through the mesh. The present method is based on a finite‐difference approach on a staggered mesh together with a fractional‐step method. It must be noted that the immersed boundary is generally not coincident with the position of the solution variables on the grid, therefore, an appropriate strategy is needed to construct a relationship between the curved boundary and the grid points nearby. Furthermore, a momentum forcing is added on the body boundaries and also inside the body to satisfy the no‐slip boundary condition. The immersed boundary is represented by a series of interfacial markers, and the markers are also used as Lagrangian forcing points. A linear interpolation is then used to scale the Lagrangian forcing from the interfacial markers to the corresponding grid points nearby. This treatment of the immersed‐boundary is used to simulate several problems, which have been validated with previous experimental results in the open literature, verifying the accuracy of the present method. Copyright © 2006 John Wiley & Sons, Ltd.     

非结构混合网格在气动热标模模拟中的应用

非结构混合网格消除了结构网格节点的结构性限制,可以较好地处理边界,同时兼顾了粘性边界层模拟的需求,具有灵活性大、对复杂外形适应能力强和生成耗时短等优点,在飞行器气动特性模拟中得到广泛应用.本文针对非结构混合网格的特点,把前期针对非结构混合网格气动力高精度模拟发展改进的梯度计算方法和Roe格式熵修正方法推广应用到气动热流的数值模拟.以典型钝锥标模外形的高超声速绕流为研究对象,开展了不同网格形式和第一层网格不同间距的影响研究.结果 表明,热流计算时,头部物面网格最好采用四边形或四边形交叉剖分得到的三角形网格,物面法向的网格雷诺数取20左右,为热流计算时非结构混合网格的生成提供了指导,同时验证了计算方法的有效性和可靠性.     

运动车辆横向振动半车模型分析

通过半车模型, 数值研究平滑路面上运动车辆车体的前两阶横向振动频率. 将车体模型化为两端自由的Euler-Bernoulli梁, 半车模型的车轮模型化为两个弹性不等的弹簧. 建立半车模型的数学模型描述车体的横向振动. 以两端自由的静态梁的模态为试函数和权函数, 通过高阶Galerkin截断计算车体横向振动的频率, 并研究车辆运行速度、车体刚度、弹簧刚度等参数对车体振动频率的影响.     

Mesh Deformation Method Based on Mean Value Coordinates Interpolation

Mesh deformation technique is widely used in many application fields, and has received a lot of attentions in recent years. This paper focuses on the methodology and algorithm of algebraic type mesh deformation for unstructured mesh in numerical discretization. To preserve mesh quality effectively, an algebraic approach for two and three dimensional unstructured mesh is developed based on mean value coordinates interpolation combined with node visibility analysis.The proposed approach firstly performs node visibility analysis to find out the visible boundary for each grid point to be moved, then evaluates the mean value coordinates of each grid point with respect to all vertices on its visible boundary. Thus the displacements of grid points can be calculated by interpolating the boundary movement by the mean value coordinates. Compared with other methods, the proposed method has good deformation capability and predictable computational cost, with no need to select parameters or functions. Applications of mesh deformation in different fields are presented to demonstrate the effectiveness of the proposed approach. The results of numerical experiments exhibit not only superior deformation capability of the method in traditional applications of fluid dynamic grid, but also great potential in modeling for large deformation analysis and inverse design problems.     

A mesh patching method for finite volume modelling of shallow water flow

A new mesh‐patching model is presented for shallow water flow described by the 2D non‐linear shallow water (NLSW) equations. The mesh‐patching model is based on AMAZON, a high‐resolution NLSW engine with an improved HLLC approximate Riemann solver. A new patching algorithm has been developed, which not only provides improved spatial resolution of flow features in particular parts of the mesh, but also simplifies and speeds up the (structured) grid generation process for an area with complicated geometry. The new patching technique is also compatible with increasingly popular parallel computing and adaptive grid techniques. The patching algorithm has been tested with moving bores, and results of test problems are presented and compared to previous work. Copyright © 2005 John Wiley & Sons, Ltd.