共查询到20条相似文献,搜索用时 216 毫秒
1.
基于非结构化同位网格的SIMPLE算法 总被引:4,自引:1,他引:4
通过基于非结构化网格的有限体积法对二维稳态Navier—Stokes方程进行了数值求解。其中对流项采用延迟修正的二阶格式进行离散;扩散项的离散采用二阶中心差分格式;对于压力-速度耦合利用SIMPLE算法进行处理;计算节点的布置采用同位网格技术,界面流速通过动量插值确定。本文对方腔驱动流、倾斜腔驱动流和圆柱外部绕流问题进行了计算,讨论了非结构化同位网格有限体积法在实现SIMPLE算法时,迭代次数与欠松弛系数的关系、不同网格情况的收敛性、同结构化网格的对比以及流场尾迹结构。通过和以往结果比较可知,本文的方法是准确和可信的。 相似文献
2.
针对不同流固耦合问题,提出一种基于任意拉格朗日--欧拉(ALE)有限元技术的分区强耦合算法. 运用半隐式特征线分裂算法求解ALE描述下的不可压缩黏性流体Navier-Stokes方程. 分别考虑一般平面运动刚体和几何非线性固体,采用复合隐式时间积分法推进结构运动方程,故可选用较大时间步长;进一步应用单元型光滑有限元法求解几何非线性固体大变形,获得更精确结构解且不影响计算效率. 运用子块移动技术结合正 交--半扭转弹簧近似法高效更新流体动网格;同时将一质量源项引入压力泊松方程满足几何守恒律,无需复杂构造网格速度差分格式. 采用简单高效的固定点法配合Aitken动态松弛技术实现各场耦合,可灵活选择先进单场求解技术,具备较好程序模块性. 运用本文算法分别模拟了H型桥梁截面颤振问题和均匀管道流内节气阀涡激振动问题. 研究表明,数值结果与已有文献数据吻合,计算精度和求解效率均令人满意. 相似文献
3.
针对不同流固耦合问题,提出一种基于任意拉格朗日-欧拉(ALE)有限元技术的分区强耦合算法.运用半隐式特征线分裂算法求解ALE描述下的不可压缩黏性流体Navier-Stokes方程.分别考虑一般平面运动刚体和几何非线性固体,采用复合隐式时间积分法推进结构运动方程,故可选用较大时间步长;进一步应用单元型光滑有限元法求解几何非线性固体大变形,获得更精确结构解且不影响计算效率.运用子块移动技术结合正交-半扭转弹簧近似法高效更新流体动网格;同时将一质量源项引入压力泊松方程满足几何守恒律,无需复杂构造网格速度差分格式.采用简单高效的固定点法配合Aitken动态松弛技术实现各场耦合,可灵活选择先进单场求解技术,具备较好程序模块性.运用本文算法分别模拟了H型桥梁截面颤振问题和均匀管道流内节气阀涡激振动问题.研究表明,数值结果与已有文献数据吻合,计算精度和求解效率均令人满意. 相似文献
4.
5.
基于同位网格下求解N-S方程的快速算法 总被引:1,自引:1,他引:0
在有限容积法基础上建立了基于同位网格的SIMPLEM算法。此算法使初始压力场与速度场耦合,让压力场和速度场同时更好地满足动量方程和连续性方程,且兼顾考虑扩散对流项对计算节点速度修正值的影响及源项与速度场之间的同步性,详细给出了算法的推导过程且对方腔顶盖驱动流进行了数值模拟。计算节点的布置采用同位网格技术,界面流速通过动量插值确定,在不同条件下讨论了迭代次数与残差的关系和不同算法的收敛性,同时验证了算法及程序是准确和可信的。 相似文献
6.
气体动理学统一算法的隐式方法研究 总被引:1,自引:0,他引:1
目前的气体动理学统一算法(unified gas kinetic scheme, 简称UGKS) 在求解高速流动问题时的计算效率,难以满足求解复杂工程问题的需求. 为了提高该算法的计算效率, 本文对模型方程的对流项和碰撞项进行了隐式处理, 并针对UGKS 界面通量与演化时间相关的特点, 引入了演化时间平均界面通量, 通过对控制方程矩阵进行近似LU 分解(lower-upper decomposition), 实现了隐式UGKS. 不同来流马赫数的圆柱绕流算例测试表明, 只要演化时间选取得当, 隐式方法可以得到与显式方法完全相同的结果, 且计算效率可以提高1~2 个量级. 相似文献
7.
串列双圆柱绕流下游圆柱两自由度涡致振动研究 总被引:1,自引:2,他引:1
数值研究了串列双圆柱绕流下游圆柱两自由度涡致振动问题,研究发现:(1) 双自由度的圆柱振幅峰值及出现振峰的频率比都比单自由度的大;(2) 尾流圆柱中的升力远大于均匀来流的,而阻力却相反;(3) 下游圆柱的位移响应对于频率比的变化没有均匀来流中的"敏感";(4) 尾流中,在频率比1.16和0.87之间,出现了明显的"拍"现象,即圆柱的振幅响应包含不同的频率,而在均匀来流中,并无明显的"拍"现象. 采用ALE方法,计算网格采用H-O非交错网格系统,结合分块耦合方法. N-S方程的对流项和扩散项分别采用三阶迎风紧致格式和四阶中心紧致格式离散. 圆柱振动采用弹簧柱体阻尼器模型,柱体的振动方程采用龙格-库塔法求解. 通过模拟柱体和流体之间的非线性耦合作用,成功地捕捉到了"拍"和"相位开关"等现象. 相似文献
8.
9.
10.
11.
三维流固耦合问题的非结构网格数值算法在很多工程领域都有重要应用, 目前现有的数值方法主要基于分区算法, 即流体和固体区域分别进行求解, 因此存在收敛速度较慢以及附加质量导致的稳定性问题, 此外, 该类算法的并行可扩展性不高, 在大规模应用计算方面也受到一定限制.本文针对三维非定常流固耦合问题, 提出一种基于区域分解的全隐全耦合可扩展并行算法.首先基于任意拉格朗日?欧拉框架建立流固耦合控制方程, 然后时间方向采用二阶向后差分隐式格式、空间方向采用非结构稳定化有限元方法进行离散.对于大规模非线性离散系统, 构造一种结合非精确Newton法、Krylov子空间迭代法与区域分解Schwarz预条件子的Newton-Krylov-Schwarz (NKS) 并行求解算法, 实现流体、固体和动网格方程的一次性整体求解.采用弹性障碍物绕流的标准测试算例对数值方法的准确性进行了验证, 数值性能测试结果显示本文构造的全隐全耦合算法具有良好的稳定性, 在不同的物理参数下具有良好的鲁棒性, 在“天河二号”超级计算机上, 当并行规模从192增加到3072个处理器核时获得了91%的并行效率.性能测试结果表明本文构造的NKS算法有望应用于复杂区域流固耦合问题的大规模数值模拟研究中. 相似文献
12.
Jong Sun Park Seung Mo Kim Min Soo Kim 《International Journal of Computational Fluid Dynamics》2013,27(3):263-276
A numerical study was performed on flow and heat transfer involving moving free surfaces that occurs in mold filling processes such as casting and injection molding. In these problems, the calculation domain changes continuously and the numerical treatment of the moving interface tends to cause artificial diffusion. Among the solution algorithms based on the Eulerian method, the volume-of-fluid (VOF) method was used because the method is simple and efficient in handling the complex flow patterns inside the cavity. To solve the transport equation of free surface without artificial smearing of the interface the baby-cell method was employed in the geometric reconstruction of the free surface. Furthermore, a predictor–corrector method was adopted in the time integration of volume-of-fluid (VOF) transport equation to increase the accuracy. The proposed scheme was verified through several benchmark problems. In order to show the capability of the proposed method, several three-dimensional mold filling processes were solved. The current algorithm was applied to the floating body problem. Three-dimensional floating body problems were tested. 相似文献
13.
作为一种混合拉格朗日欧拉法,物质点法在流固耦合问题中具有重要的应用前景。对于自由液面的流动问题,基于物质点法框架已建立了弱可压物质点法和完全不可压物质点法,但在处理流固耦合问题时遇到了困难。弱可压物质点法由于采用可压缩状态方程,导致求解时间步长过小,压力振荡严重,产生了非物理的飞溅现象;完全不可压物质点法基于投影算法和不可压条件,消除了弱可压物质点法的压力振荡,提高了时间步长,但难以处理移动边界问题。基于变分形式的投影算法提出了一种新型流固耦合不可压物质点法,得到了体积加权的压力泊松方程PPE(Pressure Poisson Equation),解决了完全不可压物质点法无法处理不规则边界和移动边界的问题。采用流固耦合不可压物质点法研究了运动刚体容器中的液体晃动问题,并与已有实验和数值结果进行对比,验证了算法的正确性和精度。 相似文献
14.
给出了一种新的适用于流体力学问题的并行自适应有限元算法。首先,基于初始稀网格上获得的事后误差估算值,应用反复谱对剖分方法对初网格进行划分,使各子域上总体误差近似相等,从而解决并行自适应计算中的负载平衡问题。然后在各处理器上独立地求解整体问题,并进行指定子域上的网格自适应处理。最后将各子域上的自适应网格组合成一个整体网格,应用基于粘接元技术的区域分裂法在该网格上获得最终解。文末给出了数值实验结果。 相似文献
15.
The new capability has been added as the numerical method for modeling volumeless and thin rigid bodies to the direct forcing immersed boundary (DFIB) method. The DFIB approach is based on adding a virtual force to the Navier–Stokes equations of incompressible flow to account for the interaction between the fluid and structures. The volume of a solid function (VOS) identifies the stationary or moving solid structures in a given fluid domain. A new VOS-based algorithm was developed to identify thin, rigid structure boundary points in fluid flow and ensure that the fluid cannot cross through the boundary of a thin rigid structure while moving or stationary. The DFIB method was first validated in a three-dimensional (3D) turbulent flow over a circular cylinder. The large-eddy simulation simulated the turbulent flow scales. The proposed algorithm was tested using a 3D turbulent flow past a stationary and rotating Savonius wind turbine that functions as a thin, rigid body. The validation results showed that the selected DFIB approach, combined with the novel algorithm, could simulate a thin, volumeless, rigid structure that is stationary and rotating in incompressible turbulent flows. The current method is also applicable for two-way fluid-structure interaction problems. 相似文献
16.
17.
《Journal of Fluids and Structures》2007,23(6):839-857
In this paper, we present a detailed derivation of the numerical method, Immersed Finite Element Method (IFEM), for the solution of fluid-structure interaction problems. This method is developed based on the Immersed Boundary (IB) method that was initiated by Peskin, with additional capabilities in handling nonuniform and independent meshes and applying arbitrary boundary conditions on both fluid and solid domains. A higher order interpolation function is adopted from one of the mesh-free methods, the Reproducing Kernel Particle Method (RKPM), which relieves the uniformity constraint of fluid meshes. Two 2-D example problems are presented to illustrate the capabilities of the algorithm. The accuracy in the numerical analysis demonstrates that the IFEM algorithm is a reliable and robust numerical approach to solve fluid and deformable solid interactions. 相似文献
18.
19.
An analytical interface reconstruction algorithm in the PLIC‐VOF method for 2D polygonal unstructured meshes 下载免费PDF全文
Piecewise linear interface calculation (PLIC) schemes have been extensively employed in the volume‐of‐fluid (VOF) method for interface capturing in numerical simulations of multiphase flows. Polygonal unstructured meshes are often adopted because of their geometric flexibility and superiority in gradient calculation. An analytical interface reconstruction algorithm in the PLIC‐VOF method for arbitrary convex polygonal cells has been proposed in this study. The line interface at a given orientation within a polygonal cell is located by an analytical technique. It has been tested successfully for four different geometric shapes that are common in polygonal meshes. The computational efficiency of the present algorithm has been compared with several published schemes in the literature. The proposed algorithm has been shown to yield higher accuracy with reduction in computational complexity. A numerical simulation of a dam‐breaking problem has been performed using the proposed analytical PLIC technique on polygonal meshes. The results are in good agreement with experimental data available in the literature, which serves as a demonstration of its performance in a real multiphase flow. 相似文献
20.
A modified front‐tracking method was proposed for the simulation of fluid‐flexible body interactions with large deformations. A large deformable body was modeled by restructuring the body using a grid adaptation. Discontinuities in the viscosity at the fluid‐structure interface were incorporated by distributing the viscosity across the interface using an indicator function. A viscosity gradient field was created near the interface, and a smooth transition occurred between the structure and the fluid. The fluid motion was defined on the Eulerian domain and was solved using the fractional step method on a staggered Cartesian grid system. The solid motion was described by Lagrangian variables and was solved by the finite element method on an unstructured triangular mesh. The fluid motion and the structure motion were independently solved, and their interaction force was calculated using a feedback law. The interaction force was the restoring force of a stiff spring with damping, and spread from the Lagrangian coordinates to the Eulerian grid by a smoothed approximation of the Dirac delta function. In the numerical simulations, we validated the effect of the grid adaptation on the solid solver using a vibrating circular ring. The effects of the viscosity gradient field were verified by solving the deformation of a circular disk in a linear shear flow, including an elastic ring moving through a channel with constriction, deformation of a suspended catenary, and a swimming jellyfish. A comparison of the numerical results with the theoretical solutions was presented. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献