一种简单快速的多块网格生成方法

利用有限元中的形函数进行网格插值映射分块生成网格,构造了一种块网格拓扑结构,使得块与相邻块之间自动组合,并利用该方法对二维、三维复杂区域进行分块网格生成,说明该方法生成的网格质量较好,适合于复杂区域分块求解。     

基于多块嵌套结构动网格的拦截弹防热保护罩的二维分离模拟

在末制导阶段,拦截弹上的红外热保护罩要进行分离,分离引起的气动力变化将直接影响到导弹的制导精度.采用多块嵌套结构动网格技术数值模拟了拦截导弹的热保护罩的二维分离过程.给出了在整个分离过程中,两片保护罩运动的轨迹,并且给出了各个时间步对应的压力云图.预测了由于前激波而产生的作用在导弹上气动力的振荡.跟踪了在各个时刻流场对于导弹上红外窗口的影响.     

多网格防阻块护栏系统的冲击响应

防阻块薄壁结构是护栏系统中重要的吸能元件,通过防阻块的优化设计提升护栏防护能力,最大限度地减小对乘员的伤害具有重要的意义。应用有限元软件ABAQUS,建立了单跨多网格防阻块波形梁护栏系统受质量块冲击的有限元分析模型,计算了壁厚分别为2、3和4 mm的多网格防阻块护栏系统的冲击响应过程,分析了3种不同壁厚的防阻块在给定冲击速度下,内能的变化情况。通过对比分析发现:冲击过程中,壁厚2 mm的双网格防阻块最终内能质量比明显高于其他情形,能量吸收能力最强;并且壁厚2 mm双网格防阻块情形的质量块加速度峰值比壁厚3 mm单网格防阻块情形小约47.6%,更有利于乘员安全。此外,根据WSTC曲线,模拟工况下,冲击引起的人体头部减速度均在安全忍受范围内。     

格子玻尔兹曼法多块网格的交界结构优化研究

基于格子玻尔兹曼方法LBM(Lattice Boltzmann Method)对多块网格方法(Multi-Block)的粗细网格交界结构进行了研究,提出了一种新的优化处理方案。解决了原有网格交界结构存在的三个问题,即两套插值运算造成的程序结构复杂的问题,存储前几个时间步的节点流场数据以备插值运算造成内存浪费的问题和基于时间插值结果进行空间插值计算造成插值误差积累的问题。用一次多点二维空间插值的方式,将原方法的空间和时间双插值,简并成一次空间插值。通过对经典的非定常的圆柱绕流算例和定常的标准顶盖方腔驱动流算例的仿真模拟,验证了交界面处质量、动量及应力的连续性以及网格交界面数据过渡的流畅度,最终验证了改进方法的正确性。数值模拟结果表明,改进后多块算法可实现局部网格细化,进一步推动LBM方法在实际工程问题中的应用。     

在动网格上的多介质界面数值处理方法研究

运用动网格上的ALE方法对一维可压缩多介质Riemann问题进行求解,在处理介质界面上的数值通量时提出了3种不同的数值方法：Lagrange方法、GFM和HLLC方法,并且对这3种方法的数值结果进行了比较,认为GFM方法和HLLC方法在介质界面上出现大压力梯度时能够明显消除界面上密度的非物理震荡,提高介质界面上数值解的精度.     

结构动响应的样条插值加权残量算法

本文以样条函数近似表达结构的瞬态响应,然后利用加权残量法建立起多步和单步形式的逐步积分算法。着重分析了应用三次B样条函数的算法SW32,并推荐了一个具有合适特性的p方案,     

一种基于双前沿推进的边界层网格生成算法

高雷诺数粘性流动模拟对边界层内的网格正交性有特殊要求.对于复杂外形,这类问题的网格自动化生成十分困难.面向该问题,提出一种双前沿推进思想,并形成一种面向复杂几何外形的边界层网格全自动生成算法.结合多种网格技术处理局部几何特征以保证边界层网格的质量.双前沿推进思想同时适用于多块结构网格和混合网格的边界层网格生成.多个模型...     

一种新的准结构网格生成方法

在有限元分析中,高质量的结构网格可以有效地提高有限元分析的精度,但结构网格的几何适应性差,针对复杂边界的二维计算模型,现有的方法很难自动生成高质量的结构网格;而非结构网格几何适应性很好,但存在计算效率低和精度差等问题。提出了一种新的准结构网格生成方法,能够实现复杂区域的网格自动生成并且具有高网格质量。该方法首先对计算区域运用Delaunay三角剖分技术生成粗背景网格;然后利用背景网格,使用优化的Voronoi图生成过渡的蜂巢网格;最后,通过中心圆方法对蜂巢网格单元进行结构网格剖分。分析NACA0012翼型数值模拟结果表明,提出的新准结构网格生成方法能够对边界复杂的模型自动生成高质量的网格,并且通过三种不同拓扑类型网格计算结果相互对比及与实验结果对比,证明准结构网格具有高计算精度。     

基于二重网格的动网格松弛法改进

相比传统的弹簧法等方法,基于球松弛算法的动网格松弛法在复杂边界大变形条件下可以得到质量更高的边界网格以及更大的极限变形量,但该方法在时间效率上还有提升的空间.引入二重网格,采用动网格松弛法进行稀疏网格的网格变形,将边界位移传递到整个网格计算域;再利用二重网格映射,将稀疏网格位移映射到原有计算网格的节点上.算例表明,改进...     

时刻追踪多介质界面运动的动网格方法

在对可压缩多介质流动的数值模拟中,定义介质界面为一种内部边界,由网格的边组成,界面边两侧对应两种不同介质中的网格。通过求解Riemann问题追踪介质界面上网格节点的运动,同时采用局部重构的动网格技术处理介质界面的大变形问题,并将介质界面定义为网格变形边界,以防止该边界上网格体积为负。运用HLLC格式求解ALE方程组得到整个多介质流场的数值解。最后从几个多介质流模型的计算结果可以看出,本文的动网格方法是可行的,而且可以时刻追踪介质界面的运动状态。     

用鲁棒Riemann求解器和运动重叠网格计算旋翼粘性绕流

发展了一种基于鲁棒Riemann求解器和运动重叠网格技术计算直升机悬停旋翼流场的方法。基于惯性坐标系,悬停旋翼流场是非定常流场,控制方程为可压缩Reynolds平均Navier-Stoke方程,其对流项采用Roe近似Reimann求解器离散,使用改进的五阶加权基本无振荡格式进行高阶重构,非定常时间推进采用含牛顿型LUSGS子迭代的全隐式双时间步方法。为实施旋转运动和便于捕捉尾迹,计算采用运动重叠网格技术。计算得到的桨叶表面压力分布及桨尖涡涡核位置都与实验结果吻合较好。数值结果表明:所发展方法对桨尖涡具有较高的分辨率,对激波具有较好的捕捉能力,该方法可进一步推广到前飞旋翼粘性绕流的计算。     

A moving mesh BGK scheme for multi‐material flows

In this paper, a moving mesh BGK scheme (MMBGK) for multi‐material flow computations is proposed. The basic idea of constructing the MMBGK is to couple the Lagrangian method, which tracks material interfaces and keeps the interfaces sharp, with a remapping‐free ALE‐type kinetic method within each single material region, where the kinetic method is based on the BGK (Bhatnagar–Gross–Krook) model. Within each single material region, a numerical flux formulation is developed on moving meshes from motion of microscope particles, and the mesh velocity is determined by requiring both mesh adaptation for accuracy and robustness, such that the grids are moving towards to the regions with high flow gradients in a way of diffusive mechanism (velocity) to adjust the distances between neighboring cells, thus increasing the numerical accuracy. To keep the sharpness of material interfaces, the Lagrangian velocity and flux are constructed at the interfaces only. Consequently, a BGK‐scheme‐based ALE‐type method (i.e., the MMBGK scheme) for multi‐material flows is constructed. Numerical examples in one and two dimensions are presented to illustrate the accuracy and robustness of the MMBGK scheme. Copyright © 2011 John Wiley & Sons, Ltd.     

Accurate and robust adaptive mesh refinement for aerodynamic simulation with multi‐block structured curvilinear mesh

A multi‐block curvilinear mesh‐based adaptive mesh refinement (AMR) method is developed to satisfy the competing objectives of improving accuracy and reducing cost. Body‐fitted curvilinear mesh‐based AMR is used to capture flow details of various length scales. A series of efforts are made to guarantee the accuracy and robustness of the AMR system. A physics‐based refinement function is proposed, which is proved to be able to detect both shock wave and vortical flow. The curvilinear mesh is refined with cubic interpolation, which guarantees the aspect ratio and smoothness. Furthermore, to enable its application in complex configurations, a sub‐block‐based refinement strategy is developed to avoid generating invalid mesh, which is the consequence of non‐smooth mesh lines or singular geometry features. A newfound problem of smaller wall distance, which negatively affects the stability and is never reported in the literature, is also discussed in detail, and an improved strategy is proposed. Together with the high‐accuracy numerical scheme, a multi‐block curvilinear mesh‐based AMR system is developed. With a series of test cases, the current method is verified to be accurate and robust and be able to automatically capture the flow details at great cost saving compared with the global refinement. Copyright © 2015 John Wiley & Sons, Ltd.     

Moving mesh methods for Boussinesq equation

The Boussinesq equation is a challenging problem both analytically and numerically. Owing to the complex dynamic development of small scales and the rapid loss of solution regularity, the Boussinesq equation pushes any numerical strategy to the limit. With uniform meshes, the amount of computational time is too large to enable us to obtain useful numerical approximations. Therefore, developing effective and robust moving mesh methods for these problems becomes necessary. In this work, we develop an efficient moving mesh algorithm for solving the two‐dimensional Boussinesq equation. Our moving mesh algorithm is an extension of Tang and Tang (SIAM J. Numer. Anal. 2003; 41 :487–515) for hyperbolic conservation laws and Zhang and Tang (Commun. Pure Appl. Anal. 2002; 1 :57–73) for convection‐dominated equations. Several numerical fluxes (Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction (2nd edn). Springer: Berlin, 1999; WASCOM 99”: 10th Conference on Waves and Stability in Continuous Media, Porto Ercole, Italy, 1999; 257–264; High‐order Methods for Computational Physics. Springer: Berlin, 1999; 439–582; J. Sci. Comput. 1990; 5 :127–149; SIAM J. Numer. Anal. 2003; 41 :487–515; Commun. Pure Appl. Anal. 2002; 1 :57–73) are also discussed. Numerical results demonstrate the advantage of our moving mesh method in resolving the small structures. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

A staggered control volume scheme for unstructured triangular grids

The purpose of this work is to introduce and validate a new staggered control volume method for the simulation of 2D/axisymmetric incompressible flows. The present study introduces a numerical procedure for solving the Navier–Stokes equations using the primitive variable formulation. The proposed method is an extension of the staggered grid methodology to unstructured triangular meshes for a control volume approach which features ease of handling of irregularly shaped domains. Two alternative elements are studied: transported scalars are stored either at the sides of an element or at its vertices, while the pressure is always stored at the centre of an element. Two interpolation functions were investigated for the integration of the momentum equations: a skewed mass-weighted upwind function and a flow-oriented exponential shape function. The momentum equations are solved over the covolume of a side or of a vertex and the pressure–velocity coupling makes use of a localized linear reconstruction of the discontinuous pressure field surrounding an element in order to obtain the pressure gradient terms. The pressure equation is obtained through a discretization of the continuity equation which uses the triangular element itself as the control volume. The method is applied to the simulation of the following test cases: backward-facing step flow, flow over a two-dimensional obstacle and flow in a pipe with sudden contraction of cross-sectional area. All numerical investigations are compared with experimental data from the literature. A grid convergence and error analysis study is also carried out for flow in a driven cavity. Results compared favourably with experimental data and so the new control volume scheme is deemed well suited for the prediction of incompressible flows in complex geometries. © 1997 John Wiley & Sons, Ltd.     

A finite volume Navier-Stokes algorithm for adaptive grids

A compact, finite volume, time-marching scheme for the two-dimensional Navier-Stokes equations of viscous fluid flow is presented. The scheme is designed for unstructured (locally refined) quadrilateral meshes. An earlier inviscid equation (Euler) scheme is employed for the convective terms and the emphasis is on treatment of the viscous terms. An essential feature of the algorithm is that all necessary operations are restricted to within each cell, which is very important when dealing with unstructured grids. Numerical issues which have to be addressed when developing a Navier-Stokes scheme are investigated. These issues are not limited to the particular Navier-Stokes scheme developed in the present work but are general problems. Specifically, the extent of the numerical molecule, which is related to the compactness of the scheme and to its suitability for unstructured grids, is examined. An approach which considers suppression of odd-even mode decoupling of the solution when designing a scheme is presented. In addition, accuracy issues related to grid stretching as well as boundary layer solution contamination due to artificial dissipation are addressed. Although the above issues are investigated with respect to the specific scheme presented, the conclusions are valid for an entire class of finite volume algorithms. The Navier-Stokes solver is validated through test cases which involve comparisons with analytical, numerical and experimental results. The solver is coupled to an adaptive algorithm for high-Reynolds-number aerofoil flow computations.     

求解带自由面渗流问题移动网格法的改进及其工程应用

采用有限元固定网格和移动网格相结合的方法,求解带自由面的坝体稳定渗流问题。本方法是针对一般坝体的结构,在有限元网格生成之后,又自动生成了网格移动所需的信息,实现了网格在迭代计算中随自由面的变化自动均匀移动;讨论了网格奇异的不同类型,提出了避免单元形状怪异的移动网格方法;在欧拉空间中分区记录坝体不同的渗透性能,在网格移动到新位置后,自动识别网格所在的物性空间,从而确定单元在新位置上的渗流系数。本文将该方法应用到吉林丰满混凝土重力坝防渗加固的实际工程当中,在考虑大坝原有防渗措施的前提下,利用数值分析,研究了防渗芯墙对大坝渗流性能的影响。计算结果表明,加芯墙后坝体自由面位置明显下移,湿水部分大为减少,坝基扬压力没有明显提升,下游出水点位置没有明显下降。计算过程表明,本文方法及措施具有很好的收敛性与稳定性,且收敛性与自由面初始位置的选择无关。     

A new algorithm for surface tension model in moving particle methods

A new algorithm for the surface tension model was developed for moving particle methods. The algorithm is based on the link‐list search algorithm and the continuum surface tension (CST) model. The developed algorithm with the CST model was implemented to a kind of moving particle approach, the finite volume particle (FVP) method. The FVP method with the new algorithm was tested by oscillatory behaviour of a two‐dimensional droplet. The oscillatory period agrees well with analytical one, and the transient shape of the droplet is also in good agreement with that obtained by other numerical methods. The droplet impact on a liquid surface was also studied using the new algorithm. The deposition and splashing phenomena were clearly reproduced. Simulated spread radius of the splashing phenomena was consistent with a power law. Copyright © 2007 John Wiley & Sons, Ltd.