适应复杂外形的三维粘性混合网格生成算法

提出一类适应复杂外形的粘性混合网格生成算法。表面网格由前沿推进三角形曲面网格程序获得,边界层布置各向异性的三棱柱体网格,远物面区域采用Delaunay方法生成四面体网格。针对模型的复杂几何特征,综合采用了各种网格处理技术,以保证边界层网格的质量,并避免算法失效问题。网格实例及计算结果表明了本文算法的实用性及和效性。     

Stabilized finite element formulation for incompressible flow on distorted meshes

Flow computations frequently require unfavourably meshes, as for example highly stretched elements in regions of boundary layers or distorted elements in deforming arbitrary Lagrangian Eulerian meshes. Thus, the performance of a flow solver on such meshes is of great interest. The behaviour of finite elements with residual‐based stabilization for incompressible Newtonian flow on distorted meshes is considered here. We investigate the influence of the stabilization terms on the results obtained on distorted meshes by a number of numerical studies. The effect of different element length definitions within the elemental stabilization parameter is considered. Further, different variants of residual‐based stabilization are compared indicating that dropping the second derivatives from the stabilization operator, i.e. using a streamline upwind Petrov–Galerkin type of formulation yields better results in a variety of cases. A comparison of the performance of linear and quadratic elements reveals further that the inconsistency of linear elements equipped with residual‐based stabilization introduces significant errors on distorted meshes, while quadratic elements are almost unaffected by moderate mesh distortion. Copyright © 2008 John Wiley & Sons, Ltd.     

多尺度六面体单元网格连接方法研究

提出多结点六面体单元的结点形函数构建方法,解决了数字化分析时两个多尺度有限元网格在接触界面由于单元结点不匹配而导致的结点属性不能连续传递问题.首先将全局坐标下的六面体单元及其表面上多结点通过等参逆变换转成局部坐标下的规则六面体单元及其表面上多结点;在规则单元中,以每个结点为基点,分别沿三个正交的局部坐标方向在单元内寻找...     

High-order discontinuous Galerkin solver on hybrid anisotropic meshes for laminar and turbulent simulations

Efficient and robust solution strategies are developed for discontinuous Galerkin （DG） discretization of the Navier-Stokes （NS） and Reynolds-averaged NS （RANS） equations on structured/unstructured hybrid meshes. A novel line-implicit scheme is devised and implemented to reduce the memory gain and improve the computational eificiency for highly anisotropic meshes. A simple and effective technique to use the mod- ified Baldwin-Lomax （BL） model on the unstructured meshes for the DC methods is proposed. The compact Hermite weighted essentially non-oscillatory （HWENO） limiters are also investigated for the hybrid meshes to treat solution discontinuities. A variety of compressible viscous flows are performed to examine the capability of the present high- order DG solver. Numerical results indicate that the designed line-implicit algorithms exhibit weak dependence on the cell aspect-ratio as well as the discretization order. The accuracy and robustness of the proposed approaches are demonstrated by capturing com- plex flow structures and giving reliable predictions of benchmark turbulent problems.     

Stabilized finite‐element computation of compressible flow with linear and quadratic interpolation functions

The unsteady compressible flow equations are solved using a stabilized finite‐element formulation with C⁰ elements. In 2D, the performance of three‐noded linear and six‐noded quadratic triangular elements is compared. In 3D, the relative performance is evaluated for 6‐noded linear and 18‐noded quadratic wedge elements. Results are compared for the solutions to Euler, laminar, and turbulent flows at different Mach numbers for several flow problems. The finite‐element meshes considered for comparison have same location of nodes for the linear and quadratic interpolations. For the turbulent flow, the Spalart–Allmaras model is used for closure. It is found that the quadratic elements yield better performance than the linear elements. This is attributed to accurate representation of the stabilization terms that involve second‐order derivatives in the formulation. When these terms are dropped from the formulation with quadratic interpolation, the numerical results are similar to those obtained with linear interpolation. The absence of these terms result in added numerical diffusion that seems to give the effect of a relatively reduced Reynolds number. For the same location of nodes, the computations with the linear triangular and wedge elements are approximately 20% and 100% faster than those with quadratic triangular and wedge elements, respectively. However, if the same quadrature rule for numerical integration is used for both interpolations, the computations with quadratic elements are approximately 20% and 45% faster in 2D and 3D, respectively. Copyright © 2014 John Wiley & Sons, Ltd.     

Wave envelope and infinite elements for acoustical radiation

Finite element models are presented for the calculation of near and far field acoustical radiation. These models are applied to the specific problem of fan noise radiation from axisymmetric turbofan inlets. In all cases conventional acoustic finite elements are used within an inner region close to the inlet. The far field is represented by infinite elements or wave envelope elements. Theory and results are presented for the case with zero mean flow. Comparisons of computed data with analytic solutions and measured values establish the utility of both the infinite element and wave envelope element schemes in determining the near field values of acoustical pressure. The wave envelope scheme is shown to be effective also in the far field. Both schemes use meshes an order of magnitude more sparse that would be required in conventional numerical discretizations, and may consequently be applied at modest computational cost.     

一种中高频激励下动力系统响应的级数展开改进算法

针对频率响应函数的级数展开法在中高频激励时计算发散的问题,提出一种新的级数展开改进算法.将系统的结构模态划分为低阶和截断的高阶模态,在模态叠加分析的基础上,将频率响应函数进行泰勒级数展开.根据高低阶模态对质量矩阵和刚度矩阵的耦合特性,用低阶模态及系统矩阵表达高阶模态对响应的影响.研究结果表明,该算法将频率响应函数的级数展开法扩展到高频激励和中频激励范围阶段,在非完备模态条件下提高了频率响应函数的计算精度,数值计算检验了该方法准确可靠并有很好的收敛性.     

An Assessment of Detached-Eddy Simulations of Unsteady Crosswind Aerodynamics of Road Vehicles

In this paper, the availability of Detached-Eddy Simulations (DES) for simulating unsteady crosswind aerodynamics in industrial applications is explored and reported. The flow around a simple car geometry under headwind, steady crosswind and time-dependent wind gust respectively, is analysed for two different types of mesh using a commercial software, STAR-CD. The typical Reynolds number of the flow cases studied is 2.0 · 10⁶ based on the vehicle length. The performances of polyhedral meshes together with a comparable hexahedral mesh are evaluated for capturing the time development of the flow structures during a gust. From the results presented, it can be concluded that the aerodynamic coefficients can be reliably calculated using both polyhedral and hexahedral meshes, although the fine turbulent structures in the wake are better captured with the hexahedral mesh. Using turbulence methods like DES in order to resolve the flow scales provides a significant insight for designing a ground vehicle and, due to the reasonable computational times involved, can hence be foreseen to be incorporated in a design process in a near future.     

二维正交各向异性位势问题的高阶单元快速多极边界元法

研制了一种适用于二维正交各向异性位势问题的高阶单元(线性单元和二次单元)快速多极边界元法. 在快速多极边界元法中, 源点对于远场区域的积分采用快速多极展开式计算, 而对于近场区域的积分则直接进行计算. 高阶单元的使用使得近场积分, 尤其是奇异积分和几乎奇异积分的计算更加复杂. 通过引入复数表达对其进行简化, 若边界采用线性单元插值, 近场积分可直接解析计算; 若采用二次单元插值, 则给出一个半解析算法计算近场积分. 高阶单元奇异积分和几乎奇异积分计算难题的解决, 使得高阶单元快速多极边界元法不仅能够计算一般结构, 也能被应用于超薄体结构, 拓宽了高阶单元快速多极边界元法的适用范围. 数值算例表明, 若计算精度一定, 高阶单元快速多极边界元法较常值单元快速多极边界元法使用的单元数量显著减少, 且高阶单元快速多极边界元法计算时间与自由度数量成线性关系, 其计算效率仍处于$O(N)$量级, 因此高阶单元快速多极边界元法可更加高效求解大规模问题.      

Stabilized mixed three‐field formulation for a generalized incompressible Oldroyd‐B model

This paper presents a generalization of the incompressible Oldroyd‐B model based on a thermodynamic framework within which the fluid can be viewed to exist in multiple natural configurations. The response of the fluid is viewed as a combination of an elastic component and a dissipative component. The dissipative component leads to the evolution of the underlying natural configurations, while the response from the natural configuration to the current configuration is considered elastic and therefore non‐dissipative. For an incompressible fluid, it is necessary that both the elastic behavior as well as the dissipative behavior is isochoric. This is achieved by ensuring that the determinant of the stretch tensor associated with the elastic response meets the constraint that its determinant is unity. A new stabilized mixed method is developed for the velocity, pressure and the kinematic tensor fields. Analytical models for fine scale fields are derived via the solution of the fine‐scale equations facilitated by the Variational Multiscale framework that are then variationally embedded in the coarse‐scale variational equations. The resulting method inherits the attributes of the classical SUPG and GLS methods, while a significant new contribution is that the form of the stabilization tensors is explicitly derived. A family of linear and quadratic tetrahedral and hexahedral elements is developed with equal‐order interpolations for the various fields. Numerical tests are presented that validate the incompressibility of the elastic stretch tensor, show optimal L₂ convergence for the conformation tensor field, and present stable response for high Weissenberg number flows. Copyright © 2016 John Wiley & Sons, Ltd.     

A quadratic-element formulation of the complex variable boundary element method

The complex variable boundary element method (CVBEM) for simply connected domains is extended to include the use of quadratic elements and interpolating functions. The derivation follows the format for linear elements given in the literature, with second-degree Lagrange polynomials taken as the interpolating functions. The quadratic-element CVBEM nodal- and interior-point equations are given in detail, and the resulting formulation is successfully tested by solving example problems with available analytical solutions. Comparisons of computational efficiency and accuracy are made between the solutions obtained using linear and quadratic elements. Additional comparisons are made using published results from other boundary element methods.     

圆筒内横向流体与弹性管耦合的分域求解算法

海洋丛式井组隔水管和换热器管束等在流体作用下,会诱导管束振动及碰撞,导致管束断裂失效。将弹性管离散成梁单元,采用非线性结构动力学方程描述;将圆筒流体域离散成实体单元,采用计算流体动力学方程描述。在流固耦合界面处,推导了界面位移、速度和载荷计算公式及收敛判断准则,建立了圆筒内横向流体与弹性管耦合的分域求解算法。算例表明,其分域与全域求解计算结果吻合较好,本文算法为复杂流体域内多根管束的振动及碰撞问题求解提供行之有效的计算方法。     

黏性边界层网格自动生成

高雷诺数黏性流动在壁面附近存在边界层,在计算模拟中自动生成可靠且有效的计算网格仍然是计算流体力学存在的瓶颈.三棱柱/四面体混合网格技术在一定程度上缓解了这个困难.然而,对于复杂外形的情况,在边界层内自动高效生成高质量的三棱柱单元仍然十分困难.常用的层推进法在凹凸区域及角点处生成的边界层网格单元质量较差,边界层网格最外层尺寸不均匀.针对这些问题,发展了一种黏性边界层网格自动生成方法,通过顶点周围边的二面角识别物面网格特征确定多生长方向,预估并调整生长高度处理相交情况.同时提出一种三维前沿尺寸调节方式,提高了边界层网格单元的正交性,保证了边界层网格与远场网格尺寸的光滑过渡.通过ONERA M6翼型以及带发动机短舱的DLR-F6翼身组合体等外形的网格生成实例及绕流数值模拟,将计算值与标准实验值进行对比,结果表明:该方法能够自动高效地生成满足数值计算需求的混合网格.     

Discretization of incompressible vorticity–velocity equations on triangular meshes

This paper describes a new approach to discretizing first- and second-order partial differential equations. It combines the advantages of finite elements and finite differences in having both unstructured (triangular/tetrahedral) meshes and low-order physically intuitive schemes. In this ‘co-volume’ framework, the discretized gradient, divergence, curl, (scalar) Laplacian, and vector Laplacian operators satisfy relationships found in standard vector field theory, such as a Helmholtz decomposition. This article focuses on the vorticity–velocity formulation for planar incompressible flows. The algorithm is described and some supporting numerical evidence is provided.     

Nonuniform residual-stress measurement by hole-drilling method

Considering the general stress field as the summation of two terms of a power series, a method for the measurement of nonuniform stress fields in thin plates by the hole-drilling method has been devised. The relieved-strain equations for the uniform and linear terms of the assumed power series have been calculated and the related constants of these equations for a range of hole diameter have been plotted.From the relieved-strain equations, it is shown that for a linear approximation of a field, a rosette gage with at least five grid elements is needed. A special rosette is proposed for the linear approximation of the residual-stress fields. In addition the equations used to determine the uniform parts, the direction, and the slopes are given. An example of the linear approximation is presented. It is shown that for some residual-stress fields, the conventional equations based on a uniform stress field produce erroneous results. The improved equations, however, provide the correct solution.     

IMPOSING DISPLACEMENT BOUNDARY CONDITIONS WITH NITSCHE＇S METHOD IN ISOGEOMETRIC ANALYSIS

等几何分析使用 NURBS 基函数统一表示几何和分析模型, 消除了传统有限元的网格离散误差, 容易构造高阶连续的协调单元. 对于结构分析, 选择合适的几何参数可以得到光滑的应力解, 避免了后置处理的应力磨平. 但是由于 NURBS 基函数不具备插值性, 难以直接施加位移边界条件. 针对这一问题, 提出一种基于 Nitsche 变分原理的边界位移条件“弱”处理方法, 它具有一致稳定的弱形式, 不增加自由度, 方程组对称正定和不会产生病态矩阵等优点. 同时给出方法的稳定性条件, 并通过求解广义特征值问题计算稳定性系数. 最后, 数值算例表明 Nitsche 方法在h细化策略下能获得最优收敛率, 其结果要明显优于在控制顶点处直接施加位移约束.}     

Hybrid mesh generation with embedded surfaces using a multiple marching direction approach

This paper proposes a method for the creation of hybrid meshes with embedded surfaces for viscous flow simulations as an extension of the multiple marching direction approach (AIAA J. 2007; 45 (1):162–167). The multiple marching direction approach enables to place semi‐structured elements around singular points, where valid semi‐structured elements cannot be placed using conventional hybrid mesh generation methods. This feature is discussed first with a couple of examples. Elements sometimes need to be clustered inside a computational domain to obtain more accurate results. For example, solution features, such as shocks, vortex cores and wake regions, can be extracted during the process of adaptive mesh generation. These features can be represented as surface meshes embedded in a computational domain. Semi‐structured elements can be placed around the embedded surface meshes using the multiple marching direction approach with a pretreatment method. Tetrahedral elements can be placed easily instead. A couple of results are presented to demonstrate the capability of the mesh generation method. Copyright © 2008 John Wiley & Sons, Ltd.     

等几何分析中采用Nitsche法施加位移边界条件

等几何分析使用NURBS基函数统一表示几何和分析模型,消除了传统有限元的网格离散误差,容易构造高阶连续的协调单元.对于结构分析,选择合适的几何参数可以得到光滑的应力解,避免了后置处理的应力磨平.但是由于NURBS基函数不具备插值性,难以直接施加位移边界条件.针对这一问题,提出一种基于Nitsche变分原理的边界位移条件"弱"处理方法,它具有一致稳定的弱形式,不增加自由度,方程组对称正定和不会产生病态矩阵等优点.同时给出方法的稳定性条件,并通过求解广义特征值问题计算稳定性系数.最后,数值算例表明Nitsche方法在h细化策略下能获得最优收敛率,其结果要明显优于在控制顶点处直接施加位移约束.     

An efficient and accurate algebraic interface capturing method for unstructured grids in 2 and 3 dimensions: The THINC method with quadratic surface representation

We present in this paper an efficient and accurate volume of fluid (VOF) type scheme to compute moving interfaces on unstructured grids with arbitrary quadrilateral mesh elements in 2D and hexahedral elements in 3D. Being an extension of the multi‐dimensional tangent of hyperbola interface capturing (THINC) reconstruction proposed by the authors in Cartesian grid, an algebraic VOF scheme is devised for arbitrary quadrilateral and hexahedral elements. The interface is cell‐wisely approximated by a quadratic surface, which substantially improves the numerical accuracy. The same as the other THINC type schemes, the present method does not require the explicit geometric representation of the interface when computing numerical fluxes and thus is very computationally efficient and straightforward in implementation. The proposed scheme has been verified by benchmark tests, which reveal that this scheme is able to produce high‐quality numerical solutions of moving interfaces in unstructured grids and thus a practical method for interfacial multi‐phase flow simulations. Copyright © 2014 John Wiley & Sons, Ltd.     

平面广义四节点等参元GQ4及其性能探讨

广义节点有限元是将传统有限元方法中的节点广义化,在不增加节点个数的前提下,仅通过提高广义节点的插值函数的阶次,从而达到提高有限元解精度的目的.与现有的p型和hp型有限元不同,在这种新的有限元中,节点自由度全部定义在节点处,在理论与程序实现上与传统有限元方法具有很好的相容性,传统有限元方法是这种新方法的广义节点退化为0阶时的特殊情形.文中主要讨论了这一新方法的四节点等参元（记为GQ4）的形式.对GQ4进行的各种数值试验表明,所发展的广义四节点等参单元具有精度高且无剪切自锁与体积自锁等的特点.