共查询到20条相似文献,搜索用时 468 毫秒
1.
2.
流体系统在重力作用下,流体压力可与重力相平衡,达到静力平衡状态.常规的有限差分和有限体积方法无法在离散尺度上保持这种平衡态,常常产生虚假速度等非物理现象.通过将重力源项的体积分改写为网格界面值的加权平均,实现离散条件下压力与重力的精细平衡,得到一种Navier-Stokes (NS)方程的精细平衡气体动理学格式.该格式可以在机器精度上精确地保持等温条件下的静力平衡态以及捕捉平衡态附近的小扰动传播.同时,该格式还能求解自然界中更常见的非等温平衡态.此类平衡态除了要求静力平衡还要求热流平衡.利用提出的格式求解NS方程,可以得到静止(机器精度下的零速度)的非等温平衡态,并且密度和温度分布都具有二阶精度.通过多个算例验证格式的有效性,表明该格式可更好地模拟重力场下的温度、密度的分层流动. 相似文献
3.
4.
5.
提出一种数值模拟凝聚炸药爆轰问题的单元中心型Lagrange方法.利用有限体积离散爆轰反应流动方程组,基于双曲型偏微分方程组的特征理论获得离散网格节点的速度与压力,获得的网格节点速度与压力用于更新网格节点位置以及计算网格单元边的数值通量.以这种方式获得的网格节点解是一种"真正多维"的理论解,是一维Godunov格式在二维Riemann问题的推广.有限体积离散得到的爆轰反应流动的半离散系统使用一种显-隐Runge-Kutta格式来离散求解:显式格式处理对流项,隐式格式处理化学反应刚性源项.算例表明,提出的单元中心型Lagrange方法能够较好地模拟凝聚炸药的爆轰反应流动. 相似文献
6.
针对圆柱坐标系下原始变量Navier Stokes方程,在有限控制容积法和压力修正的基础上,引入多重交错网格算法及非线性方程的FAS全近似格式,并对封闭圆柱空腔内的旋转流动进行数值模拟. 相似文献
7.
一维非线性对流占优扩散方程的变网格特征差分方法 总被引:1,自引:0,他引:1
针对一维非线性对流占优扩散方程,提出了一类变网格特征差分格式,该格式能够根据解的梯度变化及时对计算网格进行调整.与均匀网格格式相比,给出的变网格特征差分格式对于对流占优扩散问题有着更好的计算效果. 相似文献
8.
9.
在二维柱坐标系下Lagrange流体力学的计算中,积分梯度法是动量方程的一种有效离散方法.积分梯度法中,IGT(Integral Gradient Total)格式不能保持柱几何下一维球对称性;IGA(Integral Gradient Average)格式可以保持一维球对称性,但当相邻网格质量相差比较大时,会得到远远脱离真实物理现象的加速度.深入研究IGA和IGT格式发现,当相邻网格边界压力取为质量加权时,即使相邻网格质量相差较大,对于一维平面和一维柱问题,IGT与IGA等价;在二维情形下,可以缩小IGT和IGA之间的差异.理论证明,IGA格式不能保持系统的动量守恒,IGT格式能保持系统的动量守恒性.数值模拟结果进一步显示了这两个格式的优缺点. 相似文献
10.
给出了柱坐标下Euler方程数值边界条件的一种处理方法.径向第一个网格点设在距离中心半点位置上.根据相应物理量的特性,在中心附近进行边界延拓,使得内点的高精度差分格式可以同样应用在网格中心附近,从而无需单侧差分格式,保持了一致的高阶精度.对于周向边界,也建立了一种周期延拓方法,使得在周向所有节点处都能够采用同样的高精度格式离散,并进行了数值试验. 相似文献
11.
求解Navier-Stokes方程组的组合紧致迎风格式 总被引:1,自引:0,他引:1
给出一种新的至少有四阶精度的组合紧致迎风(CCU)格式,该格式有较高的逼近解率,利用该组合迎风格式,提出一种新的适合于在交错网格系统下求解Navier-Stokes方程组的高精度紧致差分投影算法.用组合紧致迎风格式离散对流项,粘性项、压力梯度项以及压力Poisson方程均采用四阶对称型紧致差分格式逼近,算法的整体精度不低于四阶.通过对Taylor涡列、对流占优扩散问题和双周期双剪切层流动问题的计算表明,该算法适合于对复杂流体流动问题的数值模拟. 相似文献
12.
13.
高阶精度CE/SE算法及其应用 总被引:2,自引:0,他引:2
对时-空守恒元解元算法(CE/SE)的网格设置做较大改进,提出一种新的六面体解元和元定义;同时在解元中对物理量进行高阶Taylor展开,给出一种在时间和空间上均具有高阶精度CE/SE算法.在此基础上,把新型的高阶精度CE/SE算法推广应用于高速流动捕捉激波间断、气相化学反应流动、计及固体动态效应的流体-弹塑性流动和非稳态多相不可压缩粘性流动中.数值实践表明,提出的新型网格结构上的高阶精度CE/SE算法具有算法简单、计算精度高、计算效率和计算效果好的优点,并大大改进和拓展了CE/SE算法的应用范围. 相似文献
14.
We develop a new cell-centered control volume Lagrangian scheme for solving Euler equations of compressible gas dynamics in cylindrical coordinates. The scheme is designed to be able to preserve one-dimensional spherical symmetry in a two-dimensional cylindrical geometry when computed on an equal-angle-zoned initial grid. Unlike many previous area-weighted schemes that possess the spherical symmetry property, our scheme is discretized on the true volume and it can preserve the conservation property for all the conserved variables including density, momentum and total energy. Several two-dimensional numerical examples in cylindrical coordinates are presented to demonstrate the performance of the scheme in terms of symmetry, accuracy and non-oscillatory properties. 相似文献
15.
16.
Yih-Ferng Peng Rajat Mittal Amalendu Sau Robert R. Hwang 《Journal of computational physics》2010,229(19):7072-7101
In this work, the local grid refinement procedure is focused by using a nested Cartesian grid formulation. The method is developed for simulating unsteady viscous incompressible flows with complex immersed boundaries. A finite-volume formulation based on globally second-order accurate central-difference schemes is adopted here in conjunction with a two-step fractional-step procedure. The key aspects that needed to be considered in developing such a nested grid solver are proper imposition of interface conditions on the nested-block boundaries, and accurate discretization of the governing equations in cells that are with block-interface as a control-surface. The interpolation procedure adopted in the study allows systematic development of a discretization scheme that preserves global second-order spatial accuracy of the underlying solver, and as a result high efficiency/accuracy nested grid discretization method is developed. Herein the proposed nested grid method has been widely tested through effective simulation of four different classes of unsteady incompressible viscous flows, thereby demonstrating its performance in the solution of various complex flow–structure interactions. The numerical examples include a lid-driven cavity flow and Pearson vortex problems, flow past a circular cylinder symmetrically installed in a channel, flow past an elliptic cylinder at an angle of attack, and flow past two tandem circular cylinders of unequal diameters. For the numerical simulations of flows past bluff bodies an immersed boundary (IB) method has been implemented in which the solid object is represented by a distributed body force in the Navier–Stokes equations. The main advantages of the implemented immersed boundary method are that the simulations could be performed on a regular Cartesian grid and applied to multiple nested-block (Cartesian) structured grids without any difficulty. Through the numerical experiments the strength of the solver in effectively/accurately simulating various complex flows past different forms of immersed boundaries is extensively demonstrated, in which the nested Cartesian grid method was suitably combined together with the fractional-step algorithm to speed up the solution procedure. 相似文献
17.
In the present paper, we present some numerical methods to solve the equations of steady and unsteady flows, such as those in the microcirculatory bed and large blood vessels (arteries and veins), respectively. In the case of steady flows, the method does not need neither any boundary conditions on pressure nor any small parameter, and the main computation consists of solving some Poisson equations. In the case of unsteady flows, the scheme uses a consistent Neumann boundary condition for the pressure Poisson equation. At each time step, a Poisson and heat equation are solved for the pressure and each velocity component, respectively. The accuracy and efficiency of scheme are checked by a set of numerical tests. 相似文献
18.
19.
Basic function method is developed to treat the incompressible viscous flow. Artificial compressibility coefficient, the technique
of flux splitting method and the combination of central and upwind schemes are applied to construct the basic function scheme
of trigonometric function type for solving three-dimensional incompressible Navier-Stokes equations numerically. To prove
the method, flows in finite-length-pipe are calculated, the velocity and pressure distribution of which solved by our method
quite coincide with the exact solutions of Poiseuille flow except in the areas of entrance and exit. After the method is proved
elementary, the hemodynamics in two- and three-dimensional aneurysms is researched numerically by using the basic function
method of trigonometric function type and unstructured grids generation technique. The distributions of velocity, pressure
and shear force in steady flow of aneurysms are calculated, and the influence of the shape of the aneurysms on the hemodynamics
is studied.
Supported by the National Natural Foundation of China (Grant Nos. 40874077, 40504020, and 40536029) and the National Basic
Research Program of China (Grant No. 2006CB806304) 相似文献
20.
将求解偏微分方程的有限积分法应用于对流-扩散-反应问题,发现对于非对流占优的对流扩散问题,有限积分法的精度比QUICK法高一个数量级,比传统的有限体积法高两个数量级.处理对流占优的对流-扩散-反应问题时,对流项的离散时引进加权参数,通过调节该参数反映输运的方向性.结果表明这种改进的有限积分法的精度比传统的有限体积法至少高四个数量级,同时明显改进了原来的有限积分法的精度和稳定性.对于对流占优的对流-扩散-反应问题,即使采用粗网格,计算结果也未出现非物理振荡现象,表明改进的有限积分法具有很好的稳定性. 相似文献