Implicit scheme for meshless compressible Euler solver

In this paper, an implicit scheme is presented for a meshless compressible Euler solver based on the Least Square Kinetic Upwind Method (LSKUM). The Jameson and Yoon's split flux Jacobians formulation is very popular in finite volume methodology, which leads to a scalar diagonal dominant matrix for an efficient implicit procedure (Jameson & Yoon, 1987). However, this approach leads to a block diagonal matrix when applied to the LSKUM meshless method. The above split flux Jacobian formulation, along with a matrix-free approach, has been adopted to obtain a diagonally dominant, robust and cheap implicit time integration scheme. The efficacy of the scheme is demonstrated by computing 2D flow past a NACA 0012 airfoil under subsonic, transonic and supersonic flow conditions. The results obtained are compared with available experiments and other reliable computational fluid dynamics (CFD) results. The present implicit formulation shows good convergence acceleration over the RK4 explicit procedure. Further, the accuracy and robustness of the scheme in 3D is demonstrated by computing the flow past an ONERA M6 wing and a clipped delta wing with aileron deflection. The computed results show good agreement with wind tunnel experiments and other CFD computations.

Meshless Local Petrov-Galerkin Method for Rotating Euler-Bernoulli Beam

Free vibration problem of a rotating Euler-Bernoulli beam is solved with a truly meshless local Petrov-Galerkin method. Radial basis function and summation of two radial basis functions are used for interpolation. Radial basis function satisfies the Kronecker delta property and makes it simpler to apply the essential boundary conditions. Interpolation with summation of two radial basis functions increases the node carrying capacity within the sub-domain of the trial function and higher natural frequencies can be computed by selecting the complete domain as a sub-domain of the trial function. The mass and stiffness matrices are derived and numerical results for frequencies are obtained for a fixed-free beam and hinged-free beam simulating hingeless and articulated helicopter blades. Stiffness and mass distribution suitable for wind turbine blades are also considered. Results show an accurate match with existing literature.

求解Euler方程的区域分解方法与并行算法

将复杂形状区域划分成多块子区域，研究发展了一种多块区域之间迎风守恒型的内边界耦合方法，实现相邻子区域解的光滑过渡，使多区耦合得到总体流场的数值解。对二维翼型跨音速流动和圆板形隆起物超音速流动等进行了分区数值计算，并将计算结果与单区计算结果和实验结果作了比较。并行分区计算引入“先进先出”的同步控制等待机制，实现了高效率并行计算，还分析了影响并行效率的主要因素。

A Lagrangian lattice Boltzmann method for Euler equations

A Lagrangian lattice Boltzmann method for solving Euler equations is proposed. The key step in formulating this method is the introduction of the displacement distribution function. The equilibrium distribution function consists of macroscopic Lagrangian variables at time steps n and n + 1. It is different from the standard lattice Boltzmann method. In this method the element, instead of each particle, is required to satisfy the basic law. The element is considered as one large particle, which results in simpler version than the corresponding Eulerian one, because the advection term disappears here. Our numerical examples successfully reproduce the classical results.

一种易于并行求解Euler方程组的分区技术

在具有复杂边界的计算区域内,求解偏微分方程组时,经常需要分区和并行计算,分区方法直接关系到数值计算的并行化程度,本文在应用时间算子分裂方法求解Euler方程组的过程中,提出了一种非常容易实现并行化计算的分区技术.

二维Euler方程的非结构网格再生成自适应方法

给出了提高二维Euler方程定常解质量的非结构网格自适应方法和模拟结果。计算了无粘激波在固壁上的反射、NAC0012翼型跨声速绕流和马赫数为3的前台阶绕流，自适应效果较好。

Euler方程的双分布函数格子Boltzmann Godunov方法

应用双分布函数系统，通过Ｇｏｄｕｎｏｖ分解，构造了一维Ｅｕｌｅｒ方程的格子Ｂｏｌｔｚｍａｎｎ算法。解决了传统格子气固有的ＧＣ问题与能量方程之间的矛盾，实现了分布函数与宏观物理量之间的一一对应。

基于CE/SE方法的二维Euler型多物质流体弹塑性问题计算

将CE/SE方法推广到二维固体流体弹塑性问题的数值计算,同时结合杂交粒子水平集方法追踪物质界面和合适的边界条件,提出一套完整的二维Euler型流体弹塑性计算方案.通过长钨杆侵彻装甲钢实验的数值模拟,对方法的精度和有效性进行验证.

全机绕流Euler方程多重网格分区计算方法

全机三维复杂形状绕流数值求解只能采用分区求解的方法，本文采用可压缩Euler方程有限体积方法以及多重网格分区方法对流场进行分区计算。数值方法采用改进的van Leer迎风型矢通量分裂格式和MUSCL方法，基于有限体积方法和迎风型矢通量分裂方法，建立一套处理子区域内分界面的耦合条件。各个子区域之间采用显式耦合条件，区域内部采用隐式格式和局部时间步长等，以加快收敛速度。计算结果飞机表面压力分布等气动力特性与实验值进行了比较，二者基本吻合。计算结果表明采用分析“V”型多重网格方法，能提高计算效率，加快收敛速度达到接近一个量级。根据全机数值计算结果和可视化结果讨论了流场背风区域旋涡的形成过程。