基于ALE有限元法的流固耦合强耦合数值模拟

针对不同流固耦合问题,提出一种基于任意拉格朗日--欧拉(ALE)有限元技术的分区强耦合算法. 运用半隐式特征线分裂算法求解ALE描述下的不可压缩黏性流体Navier-Stokes方程. 分别考虑一般平面运动刚体和几何非线性固体,采用复合隐式时间积分法推进结构运动方程,故可选用较大时间步长;进一步应用单元型光滑有限元法求解几何非线性固体大变形,获得更精确结构解且不影响计算效率. 运用子块移动技术结合正 交--半扭转弹簧近似法高效更新流体动网格;同时将一质量源项引入压力泊松方程满足几何守恒律,无需复杂构造网格速度差分格式. 采用简单高效的固定点法配合Aitken动态松弛技术实现各场耦合,可灵活选择先进单场求解技术,具备较好程序模块性. 运用本文算法分别模拟了H型桥梁截面颤振问题和均匀管道流内节气阀涡激振动问题. 研究表明,数值结果与已有文献数据吻合,计算精度和求解效率均令人满意.      

虚拟区域法在流固耦合问题中的应用

将多相流领域内的虚拟区域法引入到流固耦合问题的分析中,将固体视为应变率为零的虚拟流体,对流体和虚拟流体均以速度和压强作为基本变量,采用Navier-Stokes方程作为控制方程,同时求解流体域和虚拟流体域, 得到整个计算域的流场分布,应用分布式拉格朗日乘子法在虚拟流体域上施加刚体约束, 以保持虚拟流体的刚体外形和运动形式,最终建立一种流固耦合模型及其数值求解方法. 通过对粒子流问题和流固耦合问题进行数值模拟,验证了此模型的正确性和求解大变形/运动流固耦合问题的有效性.      

线性与非线性流固耦合动力学数值方法的进展及应用

本文综述了线性与非线性流固耦合问题数值方法的进展及工程应用。讨论了四种数值分析方法：(1)混合有限元–子结构–子区域数值模型,以求解有限域线性流固耦合问题,如流体晃动,声腔–结构耦合,流体中的压力波,化工容器的地震响应,坝水耦合等;(2)混合有限元–边界元数值模型,以求解涉及无限域的线性流固耦合问题,如大型浮体承受飞机降落冲击,船舰的炮击回应等;(3)混合有限元–有限差分(体积)数值模型,以求解不涉及破浪和两相分离的非线性流固耦合问题;(4)混合有限元–光滑粒子数值模型,以求解涉及破浪和两相分离的非线性流固耦合问题。文中推荐分区迭代求解过程,以便应用现有的固体及流体求解器,于毎一时间步长分别求解固体及流体的方程,通过耦合迭代收敛,向前推进以达问题求解。文中选用的工程应用例子包含气–液–壳三相耦合,液化天然气船水晃动,人体步行冲击引起的声腔–建筑结构耦合,大型浮体承受飞机降落冲击的瞬态动力回应,涉及破浪和两相分离的气–翼耦合及结构于水上降落的冲击。数值分析结果与可用的实验或计算结果作了比较,以说明所述方法的精度及工程应用价值。文中列出了基于流固耦合的波能采积装置模型,以应用线性系统的共振及非线性系统的周期解原理,有效地采积波能。本文列出了231篇参考文献,以便读者进一步研讨所感兴趣方法。     

基于ALE有限元法的流固耦合强耦合数值模拟

针对不同流固耦合问题,提出一种基于任意拉格朗日-欧拉(ALE)有限元技术的分区强耦合算法.运用半隐式特征线分裂算法求解ALE描述下的不可压缩黏性流体Navier-Stokes方程.分别考虑一般平面运动刚体和几何非线性固体,采用复合隐式时间积分法推进结构运动方程,故可选用较大时间步长;进一步应用单元型光滑有限元法求解几何非线性固体大变形,获得更精确结构解且不影响计算效率.运用子块移动技术结合正交-半扭转弹簧近似法高效更新流体动网格;同时将一质量源项引入压力泊松方程满足几何守恒律,无需复杂构造网格速度差分格式.采用简单高效的固定点法配合Aitken动态松弛技术实现各场耦合,可灵活选择先进单场求解技术,具备较好程序模块性.运用本文算法分别模拟了H型桥梁截面颤振问题和均匀管道流内节气阀涡激振动问题.研究表明,数值结果与已有文献数据吻合,计算精度和求解效率均令人满意.     

薄膜结构流固耦合的CFD数值模拟研究

基于弱耦合分区求解策略,在CompaqVisualFortran6．5环境下搭建了薄膜结构三维流固耦合效应的CFD数值模拟平台。程序采用模块化编程思想,主要包含几何建模、流体分析、结构分析和数据交换四个模块。其中几何建模模块采用自行编制的膜结构找形分析程序,流体分析模块采用经过二次开发的计算流体力学软件FLUENT6．0,结构分析模块采用自行编制的膜结构动力分析程序MDLFX;在数据交换模块中,编制了基于薄板样条法的插值计算程序,以实现流固交界面上不同区域网格间的数据传递问题,编制了基于代数法和迭代法的动网格变形程序,以实现流固耦合运算中的动网格更新。基于该软件平台,对单向柔性屋盖和鞍形膜结构屋盖进行了流固耦合数值模拟,验证了方法的有效性。     

非线性流体-刚体结构相互作用问题的一种数值模拟方法

给出了一种模拟非线性流体-刚体结构相互作用问题的数值方法．文中假定结构承受大的刚体运动,流体流动受非线性有粘或无粘的场方程支配并满足自由表面和两相耦合界面上的非线性边界条件,利用任意拉氏-欧氏（ALE）网格系统构造了数值模型．采用所探讨的多块数值格式,允许可动重造网格间有独立的相对运动,从而克服了流体网格与固体大运动匹配的困难．通过数值离散化,导出了描述非线性流固耦合动力学的数值方程并应用耦合迭代过程对其作了求解．通过算例,说明了所提出数值模型的应用．     

流固耦合分析的一种改进CBS有限元算法

针对流固耦合问题, 发展了一种基于任意拉格朗日-欧拉(ALE)描述有限元法的弱耦合分区算法. 运用半隐式特征线分裂算法求解Navier-Stokes方程, 在压力Poisson 方程中引入质量源项以满足几何守恒律; 运用子块移动技术更新动态网格, 并配以光滑处理防止网格质量下降; 采用Newmark-β 法求解结构运动方程. 为保持流体-结构界面处速度和动量守恒, 利用修正结合界面边界条件方法求解界面处速度通量和动量通量. 运用本方法分别模拟了不同雷诺数下单圆柱横向和两向流致振动、串列双圆柱两向流致振动. 计算表明, 本文方法计算效率高, 计算结果与已有实验和数值计算数据吻合.     

一种面向舰船结构毁伤的大变形流固耦合数值计算方法

水下爆炸导致舰船结构毁伤是一个复杂的非线性大变形流固耦合过程,高精度的流固耦合计算是获得高置信模拟结果的关键。基于浸没边界思想,本文提出一种面向大变形壳理论的流固耦合数值方法,可精确刻画流固耦合界面并高效求解流固界面约束方程。基于该方法,本文提出了完整的适用于水下爆炸舰船结构毁伤的大变形流固耦合数值计算方案,并基于大规模并行编程框架,研发形成适用于舰船结构毁伤的流固耦合大规模并行计算软件。与泰勒平板理论解和水下爆炸结构冲击响应实验数据等进行对比表明,本文方法可有效模拟大变形流固耦合工程问题,具备较高数值求解精度。在此基础上,完成了水下爆炸整船结构毁伤过程大规模数值模拟。该方法可有效应用于舰船毁伤等级评估,应用前景广阔。     

线性与非线性流固耦合动力学数值方法的进展及应用（英文）

本文综述了线性与非线性流固耦合问题数值方法的进展及工程应用.讨论了四种数值分析方法:(1)混合有限元–子结构–子区域数值模型,以求解有限域线性流固耦合问题,如流体晃动,声腔–结构耦合,流体中的压力波,化工容器的地震响应,坝水耦合等;(2)混合有限元–边界元数值模型,以求解涉及无限域的线性流固耦合问题,如大型浮体承受飞机降落冲击,船舰的炮击回应等;(3)混合有限元–有限差分(体积)数值模型,以求解不涉及破浪和两相分离的非线性流固耦合问题;(4)混合有限元–光滑粒子数值模型,以求解涉及破浪和两相分离的非线性流固耦合问题.文中推荐分区迭代求解过程,以便应用现有的固体及流体求解器,于毎一时间步长分别求解固体及流体的方程,通过耦合迭代收敛,向前推进以达问题求解.文中选用的工程应用例子包含气–液–壳三相耦合,液化天然气船水晃动,人体步行冲击引起的声腔–建筑结构耦合,大型浮体承受飞机降落冲击的瞬态动力回应,涉及破浪和两相分离的气–翼耦合及结构于水上降落的冲击.数值分析结果与可用的实验或计算结果作了比较,以说明所述方法的精度及工程应用价值.文中列出了基于流固耦合的波能采积装置模型,以应用线性系统的共振及非线性系统的周期解原理,有效地采积波能.本文列出了231篇参考文献,以便读者进一步研讨所感兴趣方法.     

考虑固体域几何非线性的固液耦合系统瞬态有限元分析

本文基于理想可压缩流体压力传输方程及增量型ＴｏｔａＩＬａｇｒａｎｇｉａｎ列式，建立了考虑固体域几何非线性的固液耦合系统有限元方程，并开展了非线性分析工作。本文的方法可用于流体域远大于固体域，从而可仅仅研究固体域非线性响应的情况。由于流体域每个节点仅有一个未知量（压力），故流体域的计算规模得以大大减缩。同时，因为本文方法仅需对固体域进行非线性分析，故显著地简化了固液耦合非线性问题的求解工作。本文的方法亦可用于三维分析。     

A new fully coupled solution of the Navier-Stokes equations

A fully coupled method for the solution of incompressible Navier-Stokes equations is investigated here. It uses a fully implicit time discretization of momentum equations, the standard linearization of convective terms, a cell-centred colocated grid approach and a block-nanodiagonal structure of the matrix of nodal unknowns. The Method is specific in the interpolation used for the flux reconstruction problem, in the basis iterative method for the fully coupled system and in the acceleration means that control the global efficiency of the procedure. The performance of the method is discussed using lid-driven cavity problems, both for two and three-dimensional geometries, for steady and unsteady flows.     

Flow simulation on moving boundary‐fitted grids and application to fluid–structure interaction problems

We present a method for the parallel numerical simulation of transient three‐dimensional fluid–structure interaction problems. Here, we consider the interaction of incompressible flow in the fluid domain and linear elastic deformation in the solid domain. The coupled problem is tackled by an approach based on the classical alternating Schwarz method with non‐overlapping subdomains, the subproblems are solved alternatingly and the coupling conditions are realized via the exchange of boundary conditions. The elasticity problem is solved by a standard linear finite element method. A main issue is that the flow solver has to be able to handle time‐dependent domains. To this end, we present a technique to solve the incompressible Navier–Stokes equation in three‐dimensional domains with moving boundaries. This numerical method is a generalization of a finite volume discretization using curvilinear coordinates to time‐dependent coordinate transformations. It corresponds to a discretization of the arbitrary Lagrangian–Eulerian formulation of the Navier–Stokes equations. Here the grid velocity is treated in such a way that the so‐called Geometric Conservation Law is implicitly satisfied. Altogether, our approach results in a scheme which is an extension of the well‐known MAC‐method to a staggered mesh in moving boundary‐fitted coordinates which uses grid‐dependent velocity components as the primary variables. To validate our method, we present some numerical results which show that second‐order convergence in space is obtained on moving grids. Finally, we give the results of a fully coupled fluid–structure interaction problem. It turns out that already a simple explicit coupling with one iteration of the Schwarz method, i.e. one solution of the fluid problem and one solution of the elasticity problem per time step, yields a convergent, simple, yet efficient overall method for fluid–structure interaction problems. Copyright © 2005 John Wiley & Sons, Ltd.     

基于混合网格Navier-Stokes方程的并行隐式计算方法研究

针对结构网格很难处理复杂外形和非结构网格无法计算具有边界层的粘性流动的缺点,发展了基于混合网格格点的隐式算法,成功地解决了在工程应用中难于处理的复杂外形粘性流场计算和效率问题。同时针对大规模的工程问题,发展了基于MPI通信技术的染色分层通讯并行计算方法。其中空间离散采用基于Roe格式发展的三阶迎风HLLEW(Harten-Lax-Van Leer-Einfeldt-Wada)或AUSM格式,湍流模型采用k??两方程湍流模型,时间推进考虑到LU-SGS并行等效较困难则采用基于DP-LUR(Data-Parallel Lower-Upper Relaxation)格式的隐式算法,计算CFL数可取到105量级,从2个到128个CPU的并行加速效率都保持在90%以上,大大提高了计算效率。算例对标模M6机翼模型流场进行计算,验证了方法的可靠性;然后对标模DLR-F6翼身组合体进行混合网格粘性与无粘计算结果进行比较,进一步验证混合网格方法;最后计算了DLR-WBNP外挂发动机翼身组合体模型,准确模拟了外挂和超临界机翼的相互干扰流动问题,采用4 CPU 16 CORE到24 CPU 96 CORE,2000步计算时间都不超...     

A parallel monolithic algorithm for the numerical simulation of large‐scale fluid structure interaction problems

A novel parallel monolithic algorithm has been developed for the numerical simulation of large‐scale fluid structure interaction problems. The governing incompressible Navier–Stokes equations for the fluid domain are discretized using the arbitrary Lagrangian–Eulerian formulation‐based side‐centered unstructured finite volume method. The deformation of the solid domain is governed by the constitutive laws for the nonlinear Saint Venant–Kirchhoff material, and the classical Galerkin finite element method is used to discretize the governing equations in a Lagrangian frame. A special attention is given to construct an algorithm with exact total fluid volume conservation while obeying both the global and the local discrete geometric conservation law. The resulting large‐scale algebraic nonlinear equations are multiplied with an upper triangular right preconditioner that results in a scaled discrete Laplacian instead of a zero block in the original system. Then, a one‐level restricted additive Schwarz preconditioner with a block‐incomplete factorization within each partitioned sub‐domains is utilized for the modified system. The accuracy and performance of the proposed algorithm are verified for the several benchmark problems including a pressure pulse in a flexible circular tube, a flag interacting with an incompressible viscous flow, and so on. John Wiley & Sons, Ltd.     

A fully implicit scheme for simulating ionized gas flows using the gas dynamics electrodynamics coupled system

The flow of ionized gases under the influence of electromagnetic fields is governed by the coupled system of the compressible flow equations and the Maxwell equations. In this system, coupling of the flow with the electromagnetic field is obtained through nonlinear and stiff source terms, which may cause difficulties with the numerical solution of the coupled system. The discontinuous Galerkin finite element method is used for the numerical solution of this system. For the magnetic field vector, discontinuous Galerkin discretization is performed using a divergence‐free vector base for the magnetic field to preserve zero divergence in the element and retain the implicit constraint of a divergence‐free magnetic field vector down to very low level both globally and locally. To circumvent difficulties resulting from the presence of the stiff source terms, implicit time marching is used for the fully coupled system to avoid wrong wave shapes and propagation speeds that are obtained when the coupling source terms are lagged in time or by using splitting iterative schemes. Numerical solutions for benchmark problems computed on collocated meshes for the flow and electromagnetic field variables with this fully coupled monolithic approach showed good agreement with other numerical solutions and exact results. Copyright © 2014 John Wiley & Sons, Ltd.     

A high order control volume finite element procedure for transient heat conduction analysis of functionally graded materials

A robust unstructured control-volume finite element method is presented for the solution of two-dimensional transient heat conduction in functionally graded materials (FGMs) with isotropic properties. The material properties at a point in the domain vary exponentially to spatial coordinates. A triangular mesh is chosen for spatial discretization and a fully implicit scheme is adapted for time discretization. Several problems are investigated and the results are successfully validated by using analytical and other numerical solutions available in the literature.     

基于非结构网格的高效求解方法研究

非结构网格的求解效率一直是计算流体力学工作者十分关注的问题。本文从一个新的角度分析了N-S(Euler/Navier-Stokes)方程求解效率的高低,表明计算效率不仅涉及时间离散的效率,空间离散和程序算法都与之息息相关。采用不同的计算状态,对目前非结构网格上广泛应用的LU-SGS、对称Gauss-Seidel和GMRES方法进行较详细地比较和分析,考查了空间离散的耗时对方程求解效率的影响。结果表明,LU-SGS方法的计算效率在所给的算例中均是最低的;在不考虑大量内存消耗时,GMRES算法求解Euler方程的效率较高,松耦合求解N-S方程时效率会有所降低;在大规模计算中,多次对称的Gauss-Seidel迭代方法应是较好的选择,特别是N-S方程的求解。     

一种基于块雅可比迭代的高阶FR格式隐式方法

最近, 基于非结构网格的高阶通量重构格式(flux reconstruction, FR)因其构造简单且通用性强而受到越来越多人的关注. 但将FR格式应用于大规模复杂流动的模拟时仍面临计算开销大、求解时间长等问题. 因此, 亟需发展与之相适应的高效隐式求解方法和并行计算技术. 本文提出了一种基于块Jacobi迭代的高阶FR格式求解定常二维欧拉方程的单GPU隐式时间推进方法. 由于直接求解FR格式空间和隐式时间离散后的全局线性方程组效率低下并且内存占用很大. 而通过块雅可比迭代的方式, 能够改变全局线性方程组左端矩阵的特征, 克服影响求解并行性的相邻单元依赖问题, 使得只需要存储和计算对角块矩阵. 最终将求解全局线性方程组转化为求解一系列局部单元线性方程组, 进而又可利用LU分解法在GPU上并行求解这些小型局部线性方程组. 通过二维无黏Bump流动和NACA0012无黏绕流两个数值实验表明, 该隐式方法计算收敛所用的迭代步数和计算时间均远小于使用多重网格加速的显式Runge-Kutta格式, 且在计算效率方面至少有一个量级的提升.      

A parallel implicit domain decomposition algorithm for the large eddy simulation of incompressible turbulent flows on 3D unstructured meshes

We present a parallel fully implicit algorithm for the large eddy simulation (LES) of incompressible turbulent flows on unstructured meshes in three dimensions. The LES governing equations are discretized by a stabilized Galerkin finite element method in space and an implicit second-order backward differentiation scheme in time. To efficiently solve the resulting large nonlinear systems, we present a highly parallel Newton-Krylov-Schwarz algorithm based on domain decomposition techniques. Analytic Jacobian is applied in order to obtain the best achievable performance. Two benchmark problems of lid-driven cavity and flow passing a square cylinder are employed to validate the proposed algorithm. We then apply the algorithm to the LES of turbulent flows passing a full-size high-speed train with realistic geometry and operating conditions. The numerical results show that the algorithm is both accurate and efficient and exhibits a good scalability and parallel efficiency with tens of millions of degrees of freedom on a computer with up to 4096 processors. To understand the numerical behavior of the proposed fully implicit scheme, we study several important issues, including the choices of linear solvers, the overlapping size of the subdomains, and, especially, the accuracy of the Jacobian matrix. The results show that an exact Jacobian is necessary for the efficiency and the robustness of the proposed LES solver.     

On the numerical solution to two fluid models via a cell centered finite volume method

A new method for the discretization of nonlinear systems of partial differential equations occurring in the numerical simulation of two phase flows is proposed. This method is based on a cell centered finite volume discretization on possibly unstructured meshes and aims to approximate three-dimensional stationary and evolution problems in arbitrary geometries. We are able to consider conservative and non-conservative systems of equations and the method belongs to the class of shock-capturing upwind ones. In the paper we put the emphasis on the treatment of terms involving first-order derivatives since we deal with the change of type (hyperbolic to non-hyperbolic). One of the features of the method is that it does not rely a priori on the hyperbolic character of the convection operator. The method is illustrated on a classical numerical benchmark and we refer to the bibliography concerning various and numerous applications in the context of two phase flows.