共查询到17条相似文献,搜索用时 140 毫秒
1.
计算平面运动激波和水柱群相互作用以及喷管流场.在Descartes网格中利用level-set方法分别追踪气/水和气/固界面,采用rGFM方法处理气/水和气/固界面边界条件.将喷管内壁简化为气/固界面并施加固壁边界条件,内壁型线数据拟合采用三次样条插值.采用5阶WENO格式分别求解Euler方程、level-set方程和界面重新初始化方程.给出激波和水柱群相互作用流场密度纹影图和指定点p-t曲线以及喷管流场压力、密度云图和速度场.改进界面法线确定方法可提高Riemann问题构造精度.可分辨运动激波和水柱群作用产生的复杂激波波系,表明激波在各水柱界面的透射和反射、在列和行水柱界面的多次反射和透射.水柱群下游区域的激波波后压力下降,表明激波加热水柱群附近气流和反向运动的反射激波造成了激波衰减.喷管流场数值解和理论解相符. 相似文献
2.
3.
多介质流动问题的求解一般是在结构网格上实现,而三角形网格对于复杂计算区域具有更好的适应性,本文结合rGFM方法,给出三角形网格上多介质流动问题界面处理方法.利用level-set方法跟踪界面,在界面处构造Riemann问题,得到界面处流体准确的流动状态.通过定义界面边界条件,将多介质流动问题转化为单介质流动问题,利用高精度RKDG方法求解.采用多个算例验证该方法的稳健性和有效性,结果表明该方法能准确捕捉界面和激波的位置,保持界面清晰. 相似文献
4.
5.
6.
7.
采用NGFM(New version of Ghost Fluid Method)处理复杂计算域的固壁边界,用RGFM(Real Ghost Fluid Method)求解气-水界面附近网格节点的状态参数,从而在直角坐标系下对复杂计算域的水下高压气泡膨胀问题进行数值模拟。流场控制方程选用Euler方程,用五阶WENO格式离散空间导数项,二阶Runge-Kutta法离散时间导数项;气-水界面追踪使用Level Set方法,对Level Set方程,用五阶HJ-WENO(Hamilton-Jacobi WENO)和三阶Runge-Kutta法求解。将计算结果与任意坐标系下的结果进行对比,验证了NGFM在笛卡尔网格中处理复杂形状固壁边界的可行性。得到了水下流场压力等值线图、高压气泡的演变过程以及特定点处的压力-时间曲线。计算结果表明,高压气泡在固壁反射激波的作用下,膨胀过程受到抑制;强激波在固壁的反射会导致固壁附近出现大范围的空化流动。 相似文献
8.
9.
10.
给出求解双向耦合可压缩气固两相流的间断有限元方法,对所得到的气固两相流方程组不需要采用分裂的方法离散,对气相、颗粒相方程及其对流部分和源项可以统一处理,两相都采用基于近似Riemann解的数值通量.数值模拟低压含尘激波管内的两相非平衡流动,并与平衡流、冻结流的结果进行比较.分析颗粒相的存在对气体运动的影响,及激波后松弛区域内两相间相互作用规律.发现颗粒质量比决定两相平衡后的最终状态,而颗粒直径决定两相流从非平衡到平衡的过渡过程,即不同尺寸颗粒对应的驰豫时间、松弛距离不同.结果表明:本文提出的计算方法对求解可压缩气固两相流是可行的,为研究复杂的气固两相流动问题奠定了基础. 相似文献
11.
This paper is devoted to developing a multi-material numerical scheme for non-linear elastic solids, with emphasis on the inclusion of interfacial boundary conditions. In particular for colliding solid objects it is desirable to allow large deformations and relative slide, whilst employing fixed grids and maintaining sharp interfaces. Existing schemes utilising interface tracking methods such as volume-of-fluid typically introduce erroneous transport of tangential momentum across material boundaries. Aside from combatting these difficulties one can also make improvements in a numerical scheme for multiple compressible solids by utilising governing models that facilitate application of high-order shock capturing methods developed for hydrodynamics. A numerical scheme that simultaneously allows for sliding boundaries and utilises such high-order shock capturing methods has not yet been demonstrated. A scheme is proposed here that directly addresses these challenges by extending a ghost cell method for gas-dynamics to solid mechanics, by using a first-order model for elastic materials in conservative form. Interface interactions are captured using the solution of a multi-material Riemann problem which is derived in detail. Several different boundary conditions are considered including solid/solid and solid/vacuum contact problems. Interfaces are tracked using level-set functions. The underlying single material numerical method includes a characteristic based Riemann solver and high-order WENO reconstruction. Numerical solutions of example multi-material problems are provided in comparison to exact solutions for the one-dimensional augmented system, and for a two-dimensional friction experiment. 相似文献
12.
通过在界面处构造Riemann问题,根据流体的法向速度和压力在界面(接触间断)处连续的特性,利用Riemann问题的解不仅定义了ghost流体的值,而且对真实流体中邻近界面的点值进行了更新,使得在界面处的流体的状态满足接触间断的性质,给出了更加精确的界面边界条件,守恒误差分析表明该方法在界面计算过程中引入较小的误差.数值试验表明该方法能准确地捕捉界面和激波的位置. 相似文献
13.
Euler多物质流体动力学数值方法中的界面处理算法 总被引:1,自引:1,他引:0
结合Euler型多物质流体动力学数值方法,将Youngs界面重构技术进行改进,改进后的算法中,混合网格周围网格物质的体积份额不但被用来计算物质界面的位置,还被用来确定混合网格中各物质的输运次序.将改进后的算法加入到自行开发的MMIC-2D通用多物质二维爆炸与冲击问题数值仿真程序中,对二维直角坐标系下圆环在平移流场中的运动过程进行模拟,以此对提出的改进界面处理算法进行数值考核.在此基础上,对聚能装药射流的形成过程进行数值模拟,模拟结果图像显示,其物质分界面清晰,并与实验结果吻合较好,从而验证了该方法的精度及有效性. 相似文献
14.
A systematic approach is presented for constructing higher-order immersed boundary and ghost fluid methods for CFD in general, and fluid–structure interaction problems in particular. Such methods are gaining popularity because they simplify a number of computational issues. These range from gridding the fluid domain, to designing and implementing Eulerian-based algorithms for challenging fluid–structure applications characterized by large structural motions and deformations or topological changes. However, because they typically operate on non body-fitted grids, immersed boundary and ghost fluid methods also complicate other issues such as the treatment of wall boundary conditions in general, and fluid–structure transmission conditions in particular. These methods also tend to be at best first-order space-accurate at the immersed interfaces. In some cases, they are also provably inconsistent at these locations. A methodology is presented in this paper for addressing this issue. It is developed for inviscid flows and prescribed structural motions. For the sake of clarity, but without any loss of generality, this methodology is described in one and two dimensions. However, its extensions to flow-induced structural motions and three dimensions are straightforward. The proposed methodology leads to a departure from the current practice of populating ghost fluid values independently from the chosen spatial discretization scheme. Instead, it accounts for the pattern and properties of a preferred higher-order discretization scheme, and attributes ghost values as to preserve the formal order of spatial accuracy of this scheme. It is illustrated in this paper by its application to various finite difference and finite volume methods. Its impact is also demonstrated by one- and two-dimensional numerical experiments that confirm its theoretically proven ability to preserve higher-order spatial accuracy, including in the vicinity of the immersed interfaces. 相似文献
15.
运动激波和气泡串相互作用的初步数值模拟 总被引:4,自引:0,他引:4
通过对激波和流体界面相互作用诱导的大变形界面演化的数值模拟,验证Level set方法精确模拟多个流体界面的有效性.采用2阶迎风TVD求解欧拉方程得到流场解,采用5阶WENO求解Level set方程追踪多流体界面,采用GFM方法处理流体内界面.利用文[1]的计算结果校核本文程序.在此基础上,对运动激波和气泡串相互作用过程进行了初步数值模拟,得到了不同时刻运动激波和圆管内的两个气泡作用后的演化图象,包括压力和密度等值线分布.计算结果表明:针对推广后的多界面Level set方程,该方法仍可高质量地捕捉多个流体界面. 相似文献
16.
Jeffrey W. Banks William D. Henshaw Donald W. Schwendeman 《Journal of computational physics》2012,231(9):3518-3547
We describe a mixed Eulerian–Lagrangian approach for solving fluid–structure interaction (FSI) problems. The technique, which uses deforming composite grids (DCG), is applied to FSI problems that couple high speed compressible flow with elastic solids. The fluid and solid domains are discretized with composite overlapping grids. Curvilinear grids are aligned with each interface and these grids deform as the interface evolves. The majority of grid points in the fluid domain generally belong to background Cartesian grids which do not move during a simulation. The FSI-DCG approach allows large displacements of the interfaces while retaining high quality grids. Efficiency is obtained through the use of structured grids and Cartesian grids. The governing equations in the fluid and solid domains are evolved in a partitioned approach. We solve the compressible Euler equations in the fluid domains using a high-order Godunov finite-volume scheme. We solve the linear elastodynamic equations in the solid domains using a second-order upwind scheme. We develop interface approximations based on the solution of a fluid–solid Riemann problem that results in a stable scheme even for the difficult case of light solids coupled to heavy fluids. The FSI-DCG approach is verified for three problems with known solutions, an elastic-piston problem, the superseismic shock problem and a deforming diffuser. In addition, a self convergence study is performed for an elastic shock hitting a fluid filled cavity. The overall FSI-DCG scheme is shown to be second-order accurate in the max-norm for smooth solutions, and robust and stable for problems with discontinuous solutions for a wide range of constitutive parameters. 相似文献