首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 93 毫秒
1.
Lagrange方法中,当流场发生大变形时,跟踪流体运动的Lagrange网格发生扭曲,使计算无法进行下去,此时必须重分网格,把网格修复成较好的形状。另外,网格自适应技术中的重构、合并与加密,以及同一问题不同程序相继计算的连接,并行计算中相邻块边界区域的数据传递等,这些情况都需要利用旧网格上的物理量来确定新网格上的物理量,是一个物理量重映过程。质点重映方法是基于物理上守恒规律的一种离散的物理量守恒映射方法,既可实现分片常数分布的一阶精度重映计算,又可实现分片线性分布的二阶精度重映计算。这种方法可严格保证守恒量的守恒性,且可以实现任意多边形网格以及节点上物理量的守恒重映。但是,基于分片线性分布的二阶精度重映方法,如果新网格的守恒量没有进行保界调整,那么相应的强度量有可能在其局部的限制范围之外,破坏了原网格物理量的单调性。因而,对二阶精度的质点重映方法进行了进一步研究。在分片线性分布的基础上,将基于结构网格的保界算法扩展到非结构网格上,给出了二阶保界的质点守恒重映方法。  相似文献   

2.
引进一种守恒的分片抛物线对流重映方法,通过交替扫描平均法提高对流重映方法的对称性,使用分片抛物线分布函数提高对流重映方法的精度.给出一维算例和二维算例检验分片抛物线对流重映方法的精度和对称性.  相似文献   

3.
提出基于细分和数值积分思想的一种离散的守恒重映方法——质点重映方法,密度分布可采用一阶精度的分片常数分布,或二阶精度的分片线性分布,分片线性密度分布函数采用面平均方法构造,重映过程中,借助四边形辅助网格,实现了交错网格节点量的重映.质点重映方法既适用于结构网格,也适用于非结构网格,且不要求新旧网格之间一一对应。数值结果表明,一阶精度重映算法健壮性好,但会产生较大的扩散效应;二阶精度重映算法可较好地保持密度分布的特性,但存在单调性问题,为改善二阶精度重映方法单调性,将结构网格质量守恒调整算法推广到非结构网格上,以限制新网格的质量密度,给出了一些重映的例子,并进行了误差分析。  相似文献   

4.
为保证重映过程的高守恒精度和单调性,并且在间断处具有极高的分辨率,基于径向基函数(RBF)插值方法构造了一类适用于任意网格的RBF守恒重映算法,通过计算守恒误差测试重映算法的守恒精度.将该方法用于光滑函数和含有间断的函数,并与其它守恒重映方法比较,表明该方法数值结果较好.  相似文献   

5.
基于MOF界面重构的多物质ALE方法   总被引:2,自引:0,他引:2  
提出一种基于MOF(Moment-of-Fluid)界面重构的多物质ALE(Arbitrary Lagrangian-Eulerian)方法.流体力学方程组采用相容有限元方法进行空间离散.提出一种新的二维子网格力学模型,用来计算混合网格中的物理量经过一个拉氏步后发生的变化,混合网格内的界面重构采用MOF方法.提出一种精确积分守恒重映方法.给出数值算例,如空气和水的Riemann问题,Dukowicz问题,水中强激波与空气泡相互作用问题等.结果表明,方法具有较高的精度,能够处理物质界面和网格的大变形问题.  相似文献   

6.
发展了一种基于MOF(Moment of Fluid)界面重构的二维中心型MMALE(Multi-Material Arbitrary Lagrangian-Eulerian)方法.其中,流体力学方程组采用中心型拉氏方法进行离散求解.混合网格的热力学封闭采用Tipton压力松弛模型.混合网格内的界面重构采用MOF方法,并对MOF方法作了简化和改进.重映步采用一种基于多边形剪裁算法的精确积分守恒重映方法.计算了若干数值例子,包括二维漩涡发展问题、Sedov问题、激波与氦气泡相互作用问题、水中强激波与空气泡相互作用问题、二维RT不稳定性问题等.数值算例表明,该方法具有二阶精度,能够计算界面两侧密度比和压力比很大的问题,并且其健壮性优于交错型MMALE方法,适合计算多介质复杂流体动力学问题.  相似文献   

7.
在计算流体力学(CFD)领域,几乎所有的方法都离不开网格,网格是各种数值方法的基础。网格质量的好坏直接影响数值结果的精度,甚至影响到数值计算的成败。为此CFD工作者发展了许多方法。如迭合网格、贴体网格和非结构网格,为了更好地数值模拟大变形问题,又进一步发展了结构/非结构混合网格的技术,尤其是发展了网格跟随流场智能化调整的网格自适应技术。这些网格技术的发展,几乎都涉及网格的变动,只要改动网格就涉及物理量的重映。重映方法一般被分为两种类型:插植重映和积分重映。所谓插植重映方法就是在计算区域D上,用已知网格上的物理量分布,通过插值理论把它插值到新网格或任意定义的规则网格上的一个过程,通过这个过程给出新网格的物理量的分布。所谓积分重映方法就是用积分的形式把旧网格上的守恒量重新映射到新网格上。如对某种体积密度分布q,简单的积分形式为  相似文献   

8.
赵丰祥  潘亮  王双虎 《计算物理》2018,35(5):525-534
基于非结构四边形网格发展求解双曲守恒律的三阶加权基本无振荡(WENO)格式.针对任意非结构四边形网格选取重构模板,并给出基于线性多项式的三阶线性重构.但对于一般的非结构四边形网格,会出现非常大的线性权和负权,使得非线性重构的WENO格式对光滑问题也不稳定.本文给出一个处理非常大的线性权的优化重构方法,对优化后得到的负线性权采用分裂方法进行处理.对于非线性权,提出一种考虑局部网格和物理量间断的新光滑度量因子.采用优化重构方法和新的非线性权,当前的三阶WENO格式在质量很差的网格上也具有很好的稳定性.理论的三阶精度在数值精度测试算例中得到验证,同时一范数和无穷范数的误差绝对值不依赖于网格质量;具有强间断的数值结果证明了当前格式的有效性.  相似文献   

9.
在ENO(Essentially Non-osciuatory)守恒插值方法的基础上,分析和研究现今流体力学计算中涉及的几类网格技术:重叠网格技术、自适应加密技术和运动网格技术.基于ENO插值多项式构造的重映方法具有良好的守恒性,可以有效保证数据传递中物理量的总体守恒.提出该类守恒插值方法在以上几种网格技术中的一些应用前景,并给出一些数值算例.  相似文献   

10.
介绍了二维非结构网格上的守恒重映算法,重点是基于SFB/DC思想的通量重映算法。用统一的公式表示不同的单元量重映算法,包括原始的贡献网格法、Barth—Jespersen方法、最小二乘法,不同算法间的区别体现为梯度求法的差异。对于交错网格上速度的重映,介绍了SALE和HIS算法。此外,为保证重映算法的有界性,引入了修补方法。  相似文献   

11.
A new, second-order accurate, volume conservative, material-order-independent interface reconstruction method for multi-material flow simulations is presented. First, materials are located in multi-material computational cells using a piecewise linear reconstruction of the volume fraction function. These material locator points are then used as generators to reconstruct the interface with a weighted Voronoi diagram that matches the volume fractions. The interfaces are then improved by minimizing an objective function that smoothes interface normals while enforcing convexity and volume constraints for the pure material subcells. Convergence tests are shown demonstrating second-order accuracy. Static and dynamic examples are shown illustrating the superior performance of the method over existing material-order-dependent methods.  相似文献   

12.
This paper presents a new variant of the volume-of-fluid (VOF) color function C advection algorithm based on the piecewise linear interface construction (PLIC) method suitable for use on general moving grids. From several existing methods for reconstructing the linear interface we adopted the least squares volume-of-fluid interface reconstruction algorithm (LVIRA) which can be easily implemented on general grids. The distinguishing step in the advection algorithm that takes into account the grid movement is the construction of the donating region containing the fluid passing through corresponding cell-faces in a single time-step. The donating regions are constructed utilizing fluid velocity in cell corners relative to grid (corner) velocities. The method is conservative as it complies with the space conservation law (SCL) and requires a proper definition of the grid velocities and fluxes due to the grid movement. The accuracy of the presented advection algorithm is assessed with standard test cases. It is comparable with other PLIC based algorithms on fixed grids, while the applicability on adaptive moving grids enables a considerable reduction in the number of grid cells.  相似文献   

13.
 WENO有限差分格式有较高的分辨精度,适合复杂流场的计算,在国际上被广泛采用。本文利用WENO有限差分格式求解2维守恒型欧拉方程,实现了对无粘流体中Kelvin-Helmholtz不稳定性的数值模拟。速度剪切方向采用周期边界条件;扰动增长方向采用嵌边出流边界条件,一个不稳定波长分布64个网格。数值模拟给出的扰动幅值线性增长率与线性稳定性分析给出的结果很好符合,显示了该格式的有效性和精度。数值模拟给出了清晰的密度等值线,表明该方法还具有较好的界面变形捕捉能力。  相似文献   

14.
A novel numerical method for two-fluid flow computations is presented, which combines the space–time discontinuous Galerkin finite element discretization with the level set method and cut-cell based interface tracking. The space–time discontinuous Galerkin (STDG) finite element method offers high accuracy, an inherent ability to handle discontinuities and a very local stencil, making it relatively easy to combine with local hp-refinement. The front tracking is incorporated via cut-cell mesh refinement to ensure a sharp interface between the fluids. To compute the interface dynamics the level set method (LSM) is used because of its ability to deal with merging and breakup. Also, the LSM is easy to extend to higher dimensions. Small cells arising from the cut-cell refinement are merged to improve the stability and performance. The interface conditions are incorporated in the numerical flux at the interface and the STDG discretization ensures that the scheme is conservative as long as the numerical fluxes are conservative. The numerical method is applied to one and two dimensional two-fluid test problems using the Euler equations.  相似文献   

15.
在自适应网格上,采用VOF方法捕捉界面,相容守恒格式计算电流及电磁力,发展了金属流体自由界面MHD数值方法。通过数值模拟磁场作用下不同Hartmann数的气泡在导电溶液中的运动和变形,分析磁场对气泡以及流场的影响,同时给出诱导电场和电流的分布。为进一步深入研究冶金及热核聚变相关的金属流体在强磁场作用下的自由界面流打下基础。  相似文献   

16.
热传导方程的一类无网格方法   总被引:1,自引:0,他引:1  
李寿佛  张瑗  刘玉珍 《计算物理》2007,24(5):573-580
构造求解热传导方程的一类无网格方法,只要选择好每个节点的适当的邻点集合,便可利用节点信息顺利进行计算.作为特殊情形,也可在各种结构或非结构网格上进行计算.在矩形或均匀平行四边形网格上进行计算时具有二阶精度,当在任意的不规则四边形或三角形网格上计算时仍然是守恒的和相容的,且至少具有一阶精度.作为数值试验,将该方法用于在不规则四边形网格上及四边形与三角形混合网格上求解二维非线性抛物型方程,并在不规则四边形网格上求解二维三温辐射热传导方程,均获得了较为理想的数值结果.  相似文献   

17.
Diffuse interface methods have been recently proposed and successfully used for accurate compressible multi-fluid computations Abgrall [1]; Kapila et al. [20]; Saurel et al. [30]. These methods deal with extended systems of hyperbolic equations involving a non-conservative volume fraction equation and relaxation terms. Following the same theoretical frame, we derive here an Eulerian diffuse interface model for elastic solid-compressible fluid interactions in situations involving extreme deformations. Elastic effects are included following the Eulerian conservative formulation proposed in Godunov [16], Miller and Colella [23], Godunov and Romenskii [17], Plohr and Plohr [27] and Gavrilyuk et al. [14]. We apply first the Hamilton principle of stationary action to derive the conservative part of the model. The relaxation terms are then added which are compatible with the entropy inequality. In the limit of vanishing volume fractions the Euler equations of compressible fluids and a conservative hyperelastic model are recovered. It is solved by a unique hyperbolic solver valid at each mesh point (pure fluid, pure solid and mixture cell). Capabilities of the model and methods are illustrated on various tests of impacts of solids moving in an ambient compressible fluid.  相似文献   

18.
In this paper, we present the development of a sharp numerical scheme for multiphase electrohydrodynamic (EHD) flows for a high electric Reynolds number regime. The electric potential Poisson equation contains EHD interface boundary conditions, which are implemented using the ghost fluid method (GFM). The GFM is also used to solve the pressure Poisson equation. The methods detailed here are integrated with state-of-the-art interface transport techniques and coupled to a robust, high order fully conservative finite difference Navier–Stokes solver. Test cases with exact or approximate analytic solutions are used to assess the robustness and accuracy of the EHD numerical scheme. The method is then applied to simulate a charged liquid kerosene jet.  相似文献   

19.
In this paper we present a method to treat interface jump conditions for constant coefficients Poisson problems that allows the use of standard “black box” solvers, without compromising accuracy. The basic idea of the new approach is similar to the Ghost Fluid Method (GFM). The GFM relies on corrections applied on nodes located across the interface for discretization stencils that straddle the interface. If the corrections are solution-independent, they can be moved to the right-hand-side (RHS) of the equations, producing a problem with the same linear system as if there were no jumps, only with a different RHS. However, achieving high accuracy is very hard (if not impossible) with the “standard” approaches used to compute the GFM correction terms.  相似文献   

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

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