首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 125 毫秒
1.
为满足亚声速和跨声速飞机概念设计中快速气动计算的需求,研究和发展一种基于自适应直角网格的非线性全速势方程有限体积解法。要点如下。(1)在几何自适应直角网格的基础上,使用结合单元融合的网格切割算法处理物面边界,提出一种修正非贴体切割网格的方法。(2)采用隐式格式结合GMRES算法求解该非线性位流方程,针对流场的自适应来捕捉激波。(3)采用镜像法处理物面边界处的无穿透条件,并提出解析的方法来修正镜像单元的值。(4)针对直角网格的特点,提出在库塔线上插入库塔单元的方法施加库塔条件。NACA0012翼型绕流的算例结果表明,该方法用于亚声速和跨声速气动计算能得到令人满意的结果,且自动化程度高、收敛速度快。  相似文献   

2.
采用自适应直角网格计算三维增升装置绕流   总被引:2,自引:0,他引:2  
针对三维增升装置绕流,对存在剪刀叉的不连续外形,基于自适应直角网格,提出并介绍了分区和面搭接技术,采用变长宽比网格,进行了直角网格生成和流场Euler方程数值计算. 根据几何外形的特点,在直角网格生成过程中,以外形不连续面作为分区边界,对初始``根'网格实施分区处理,降低了整个网格的生成难度. 通过基于外形的自适应网格加密,详细描述了剪刀叉外形和缝道,提高了网格质量. 在分区边界面上,基于面搭接技术,构造重叠面积切割算法,实现边界两侧网格间的流场信息传递,保证流场计算中的通量守恒. 采用中心有限体积方法,结合双时间推进算法,完成了两段机翼、带增升襟翼翼身组合体绕流流场的Euler方程数值模拟,对计算结果与实验数据进行了对比,验证了所提方法、算法的合理性和实用性.  相似文献   

3.
网格自适应技术在复杂外形流场模拟中的应用   总被引:2,自引:0,他引:2  
建立了一套适用于非结构混合网格自适应方法,针对激波和涡的不同特征采用不同加密探测器,各向异性加密棱柱单元并沿物面法向方向剖分所有棱柱层,各向异性剖分四面体单元,并保证四面体与棱柱交界面上网格协调。构造Hermit插值近似投影物面新加网格点和基于Laplacian光滑方法对空间网格进行优化。通过网格自适应加密,使用Roe格式计算高超声速球头绕流的红玉现象得到明显减轻。F16飞机含激波和脱体涡的流场自适应计算表明,网格加密集中在激波面和涡核附近区域,激波和涡计算更准确。  相似文献   

4.
非结构混合网格消除了结构网格节点的结构性限制,可以较好地处理边界,同时兼顾了粘性边界层模拟的需求,具有灵活性大、对复杂外形适应能力强和生成耗时短等优点,在飞行器气动特性模拟中得到广泛应用.本文针对非结构混合网格的特点,把前期针对非结构混合网格气动力高精度模拟发展改进的梯度计算方法和Roe格式熵修正方法推广应用到气动热流的数值模拟.以典型钝锥标模外形的高超声速绕流为研究对象,开展了不同网格形式和第一层网格不同间距的影响研究.结果 表明,热流计算时,头部物面网格最好采用四边形或四边形交叉剖分得到的三角形网格,物面法向的网格雷诺数取20左右,为热流计算时非结构混合网格的生成提供了指导,同时验证了计算方法的有效性和可靠性.  相似文献   

5.
基于贴体网格的VOF方法数模流场研究   总被引:1,自引:0,他引:1  
提出了一种基于VOF方法的模拟具有复杂边界形状结构物附近流场的新算法,BFC—SIMPLE—VOF算法。采用坐标变换方法实现了任意复杂区域的结构化网格划分,在贴体网格下对二维不可压缩粘性流体的控制方程进行了离散。提出了基于交错网格的修正SIMPLE算法来迭代求解压力一速度场,修正了贴体坐标下的界面跟踪方法(VOF方法)...  相似文献   

6.
提出一种新的网格自适应方法:在需要加密的网格单元中心加入新结点,并对加密后的相邻三角形网格单元进行公共边变换,构成新的网格单元.与传统的在网格边界中点加入新节点的自适应方法相比,新方法可以更加灵活地控制网格密度,加密后的网格继承原先的网格质量不发生畸变,并且算法编程简便,容易实现.将自适应网格生成方法和基于特征线方程的分离算法相结合,对空腔内不可压缩黏性流动进行了计算.在特征线方向上进行时间步离散,动量方程求解过程中采用非增量型分离算法.计算中,把求解变量梯度值作为判定准则,在变化剧烈的区域进行网格局部加密.计算结粜表明该组合算法有很好的计算精度,并有效减少了计算时间和存储量.  相似文献   

7.
在激波区使用自适应壁对跨音速翼型的激波/边界层的相互作用(干扰)进行控制,可改变机翼的气动性能,这种被动控制可通过在翼型的激波区开一凹腔,其上覆盖一弹性橡胶膜柔壁来,本文给出用Navier-Stoker方程数值模拟这一自适应控制翼型的跨音速粘性绕流,提出了一个适应于本特殊情况(物面边界局部地区在求解过程中有变化)的处理办法。并探讨了自适应柔壁对当代跨音速翼绕流的影响。  相似文献   

8.
介绍一种基于Delaunay算法的四面体自适应网格的自动划分方法。该方法用单元尺度场控制生成网格的疏密分布,在不满足尺度场要求的单元面形心处插入新节点,同时计算新节点单元尺寸参数,实现三维实体的Delaunay四面体自动划分。此方法具有几个特点:一是表面网格与体内网格同步划分,无需区分两者;二是结点与单元同时生成;三是生成网格自适应性好,疏密分布任意。另外,还介绍了三维网格划分中两个相关算法:一个是约束面恢复算法,该算法基于约束面不允许有单元边与之相交的性质而提出的;另一个是将二维射线法推广至三维空间,判断一个点是否在一多面体内,实现了凹多面体的划分。最后通过算例对单元质量进行了评价。本文所述方法是一种有效的四面体自适应单元生成算法。  相似文献   

9.
利用同位非结构化网格上的压力加权修正算法 ,对翼型湍流绕流进行了数值分析。详细地给出了一孤立翼型在不同攻角下的分离流结构及翼型表面压力分布 ,为了显示非结构化网格方法在求解多连通流动区域的优越性 ,对双翼型绕流进行了数值计算。在数值分析中 ,对阵面推进法进行改进来生成三角形网格 ,采用有限控制体方法直接在物理空间中的非结构化网格单元上离散 Navier- Stokes方程及 k- ε方程 ,形成的代数方程组通过预条件矩阵共轭梯度平方法求解。计算结果表明 :当流动为附着流时 ,计算结果与实验值吻合程度令人相当满意 ;而在分离区内 ,计算结果与实验值存在一定的误差 ,需对分离区内的湍流模型做进一步的改进。  相似文献   

10.
无网格算法在多段翼型流动计算中的应用   总被引:6,自引:1,他引:5  
研究了一种求解欧拉方程的无网格算法,发展出了一套布点及点云自动生成的方法;在点云离散的基础上,采用最小二乘法求解矛盾方程的方法来求取空间导数,进而获得数值通量;采用四步龙格-库塔方法进行时间推进,并引入当地时间步长和残值光顺等加速收敛措施。通过对NA-CA0012翼型的跨音速流动和多段翼型复杂绕流的数值模拟,验证了上述无网格算法的正确性和实用性。  相似文献   

11.
The second of a two‐paper series, this paper details a solver for the characteristics‐bias system from the acoustics–convection upstream resolution algorithm for the Euler and Navier–Stokes equations. An integral formulation leads to several surface integrals that allow effective enforcement of boundary conditions. Also presented is a new multi‐dimensional procedure to enforce a pressure boundary condition at a subsonic outlet, a procedure that remains accurate and stable. A classical finite element Galerkin discretization of the integral formulation on any prescribed grid directly yields an optimal discretely conservative upstream approximation for the Euler and Navier–Stokes equations, an approximation that remains multi‐dimensional independently of the orientation of the reference axes and computational cells. The time‐dependent discrete equations are then integrated in time via an implicit Runge–Kutta procedure that in this paper is proven to remain absolutely non‐linearly stable for the spatially‐discrete Euler and Navier–Stokes equations and shown to converge rapidly to steady states, with maximum Courant number exceeding 100 for the linearized version. Even on relatively coarse grids, the acoustics–convection upstream resolution algorithm generates essentially non‐oscillatory solutions for subsonic, transonic and supersonic flows, encompassing oblique‐ and interacting‐shock fields that converge within 40 time steps and reflect reference exact solutions. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

12.
We establish the existence and stability of multidimensional steady transonic flows with transonic shocks through an infinite nozzle of arbitrary cross-sections, including a slowly varying de Laval nozzle. The transonic flow is governed by the inviscid potential flow equation with supersonic upstream flow at the entrance, uniform subsonic downstream flow at the exit at infinity, and the slip boundary condition on the nozzle boundary. Our results indicate that, if the supersonic upstream flow at the entrance is sufficiently close to a uniform flow, there exists a solution that consists of a C 1,α subsonic flow in the unbounded downstream region, converging to a uniform velocity state at infinity, and a C 1,α multidimensional transonic shock separating the subsonic flow from the supersonic upstream flow; the uniform velocity state at the exit at infinity in the downstream direction is uniquely determined by the supersonic upstream flow; and the shock is orthogonal to the nozzle boundary at every point of their intersection. In order to construct such a transonic flow, we reformulate the multidimensional transonic nozzle problem into a free boundary problem for the subsonic phase, in which the equation is elliptic and the free boundary is a transonic shock. The free boundary conditions are determined by the Rankine–Hugoniot conditions along the shock. We further develop a nonlinear iteration approach and employ its advantages to deal with such a free boundary problem in the unbounded domain. We also prove that the transonic flow with a transonic shock is unique and stable with respect to the nozzle boundary and the smooth supersonic upstream flow at the entrance.  相似文献   

13.
An improved formulation of the inverse integral equation method proposed in Reference 1 is presented which allows, in particular, a well-posed problem to be ensured. The corresponding computation code is tested in an exhaustive manner for axial and radial compressor and turbine cascades. The agreement between the velocity field obtained with the inverse method and that resulting from a direct calculation is examined for subsonic, transonic and supersonic flows. Accuracy and reliability of the solution to the boundary condition problem are excellent for the subsonic and transonic flows. However, for the supersonic flow, the application of the method seems to be limited by the use of elementary solutions of the Laplace operator.  相似文献   

14.
A kinetic flux-vector-splitting method has been used to solve the Euler equations for inviscid, compressible flow on unstructured grids. This method is derived from the Boltzmann equation and is an upwind, cell-centered, finite volume scheme with an explicit time-stepping procedure. The Delaunay triangulation has been used to generate the grids. The approach is demonstrated for three flow field simulations, namely the subsonic flow over a two-component high-lift aerofoil, the transonic flow over an aerofoil and the supersonic flow in a channel.  相似文献   

15.
计算流体力学(computational fluid dynamics,CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度. 其次,通过NASATurbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶). 最后,通过NASA CommonResearchModel标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.   相似文献   

16.
复杂无粘流场数值模拟的矩形/三角形混合网格技术   总被引:5,自引:0,他引:5  
张来平  张涵信 《力学学报》1998,30(1):104-108
建立了一套模拟复杂无粘流场的矩形/三角形混合网格技术,其中三角形仅限于物面附近,发挥非结构网格的几何灵活性,以少量的网格模拟复杂外型;同时在以外的区域采用矩形结构网格,发挥矩形网格计算简单快速的优势,有效地克服全非结构网格计算方法需要较大内存量和较长CPU时间的不足.混合网格系统由修正的四分树方法生成.将NND有限差分与NND有限体积格式有机地融合于混合网格计算,消除了全矩形网格模拟曲边界的台阶效应,同时保证了网格间的通量守恒.数值实验表明本方法在模拟复杂无粘流场方面的灵活性和高效性.  相似文献   

17.
Dynamic fluid–solid interactions are widely found in chemical engineering, such as in particle-laden flows, which usually contain complex moving boundaries. The immersed boundary method (IBM) is a convenient approach to handle fluid–solid interactions with complex geometries. In this work, Uhlmann's direct-forcing IBM is improved and implemented on a supercomputer with CPU–GPU hybrid architecture. The direct-forcing IBM is modified as follows: the Poisson's equation for pressure is solved before evaluation of the body force, and the force is only distributed to the Cartesian grids inside the immersed boundary. A multidirect forcing scheme is used to evaluate the body force. These modifications result in a divergence-free flow field in the fluid domain and the no-slip boundary condition at the immersed boundary simultaneously. This method is implemented in an explicit finite-difference fractional-step scheme, and validated by 2D simulations of lid-driven cavity flow, Couette flow between two concentric cylinders and flow over a circular cylinder. Finally, the method is used to simulate the sedimentation of two circular particles in a channel. The results agree very well with previous experimental and numerical data, and are more accurate than the conventional direct-forcing method, especially in the vicinity of a moving boundary.  相似文献   

18.
A novel nonreflecting boundary condition, which converges to the specified time‐dependent boundary condition within any degree of accuracy, is introduced for the numerical simulation of hyperbolic systems and validated against the solution of two fundamental boundary value problems in fluids. First, transonic nozzle flow with backward acoustic disturbance is considered. Using high‐order aeroacoustic numerical schemes, the proposed nonreflecting boundary condition yields results that are in excellent agreement with those obtained using conventional nonreflecting boundary conditions based on the method of characteristics as well as with the results of the exact solution. The novel nonreflecting boundary condition, implemented into a semi‐analytical solution algorithm of unsteady bubbly cavitating nozzle flows, is also validated against results obtained using a Lagrangian finite volume scheme. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

19.
A nodally exact convection–diffusion–reaction scheme developed in Cartesian grids is applied to solve the flow equations in irregular domains within the framework of immersed boundary (IB) method. The artificial momentum forcing term applied at certain points in the flow and inside the body of any shape allows the imposition of no‐slip velocity condition to account for the body of complex boundary. Development of an interpolation scheme that can accurately lead to no‐slip velocity condition along the IB is essential since Cartesian grid lines generally do not coincide with the IB. The results simulated from the proposed IB method agree well with other numerical and experimental results for several chosen benchmark problems. The accuracy and fidelity of the IB flow solver to predict flows with irregular IBs are therefore demonstrated. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

20.
The inverse problem of the theory of the Laval nozzle is considered, which leads to the Cauchy problem for the gasdynamic equations; the streamlines and the flow parameters are found from the known velocity distribution on the axis of symmetry.The inverse problem of Laval nozzle theory was considered in 1908 by Meyer [1], who expanded the velocity potential into a series in powers of the Cartesian coordinates and constructed the subsonic and supersonic solutions in the vicinity of the center of the nozzle. Taylor [2] used a similar method to construct a flowfield which is subsonic but has local supersonic zones in the vicinity of the minimal section. Frankl [3] and Fal'kovich [4] studied the flow in the vicinity of the nozzle center in the hodograph plane. Their solution, just as the Meyer solution, made it possible to obtain an idea of the structure of the transonic flow in the vicinity of the center of the nozzle.A large number of studies on transonic flow in the vicinity of the center of the nozzle have been made using the method of small perturbations. The approximate equation for the transonic velocity potential in the physical plane, obtained in [3–6], has been studied in detail for the plane and axisymmetric cases. In [7] Ryzhov used this equation to study the question of the formation of shock waves in the vicinity of the center of the nozzle, and conditions were formulated for the plane and axisymmetric cases under which the flow will not contain shock waves. However, none of the solutions listed above for the inverse problem of Laval nozzle theory makes it possible to calculate the flow in the subsonic and transonic parts of the nozzles with large gradients of the gasdynamic parameters along the normal to the axis of symmetry.Among the studies devoted to the numerical calculation of the flow in the subsonic portion of the Laval nozzle we should note the study of Alikhashkin et al., and the work of Favorskii [9], in which the method of integral relations was used to solve the direct problem for the plane and axisymmetric cases.The present paper provides a numerical solution of the inverse problem of Laval nozzle theory. A stable difference scheme is presented which permits analysis with a high degree of accuracy of the subsonic, transonic, and supersonic flow regions. The result of the calculations is a series of nozzles with rectilinear and curvilinear transition surfaces in which the flow is significantly different from the one-dimensional flow. The flowfield in the subsonic and transonic portions of the nozzles is studied. Several asymptotic solutions are obtained and a comparison is made of these solutions with the numerical solution.The author wishes to thank G. D. Vladimirov for compiling the large number of programs and carrying out the calculations on the M-20 computer.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号