SAMR网格上扩散方程有限体格式的逼近性与两层网格算法

针对结构自适应加密网格(SAMR)上扩散方程的求解,分析几种有限体格式的逼近性,同时设计和分析一种两层网格算法.首先,讨论一种常见的守恒型有限体格式,并给出网格加密区域和细化/粗化插值算子的条件;接着,通过在粗细界面附近引入辅助三角形单元,消除粗细界面处的非协调单元,设计了一种保对称有限体元(SFVE)格式,分析表明,该格式具有更好的逼近性,且对网格加密区域和插值算子的限制更弱;最后,为SFVE格式构造一种两层网格(TL)算法,理论分析和数值实验表明该算法的一致收敛性.     

考虑密度与温度耦合关系的非结构网格SIMPLEC算法

提出一种新的密度压力速度和温度耦合的算法,使用静焓方程构造温度修正量与压力修正量之间的关系,类似Rhie-Chow方法用动量方程插值,网格面上的温度压力修正量插值采用能量方程的形式.CGNS API(CFD General Notation System Application Programming Interface)作为非结构化网格求解器的前处理和后处理,基于网格的FVM(Finite Volume Method)作为偏微分方程求解方法.燃烧模拟采用Arrinius-Eddy Dissipation模型.最后给出两个算例作为说明.     

三维不可压N-S方程的多重网格求解

应用全近似存储(Full Approximation Storage,FAS)多重网格法和人工压缩性方法求解了三维不可压Navi-er-Stokes方程.在解粗网格差分方程时,对Neumann边界条件采用增量形式进行更新,离散方程用对角化形式的近似隐式因子分解格式求解,其中空间无粘项分别用MUSCL格式和对称TVD格式进行离散.对90°弯曲的方截面管道流动和4:1椭球体层流绕流的数值模拟表明,多重网格的计算时间比单重网格节省一半以上,且无限制函数的MUSCL格式比TVD格式对流动结构有更好的分辨能力.     

非定常粘性空化流动模型及其数值计算

基于气液两相当地均相介质模型,本文给出了一种模拟非定常粘性空化流动的计算模型。认为空化绕流流场中流动介质是一种当地均匀的气液混相物。控制方程采用了应用两相流模型的N-S方程。基于液相和气相的状态方程推导了混合介质密度的表达式。为了保证数值计算的稳定性,控制方程的数值求解采用了TVD-MacCormack格式。为了评价计算模型的可靠性,分别计算了绕台阶和管道水锤的空化流动,所得结果是合理的,说明该方法可以用于空化流动的数值计算。     

球坐标动量空间下相对论Vlasov方程的数值算法

针对相对论Vlasov方程动量区间跨度大、难以计算的困难,将相对论Vlasov方程在球坐标动量空间中进行数值求解.对相对论Vlasov方程球坐标动量空间构造4阶非分裂守恒型数值格式.数值模拟相对论Landau阻尼问题并与解析理论进行比较,验证数值模型和算法的有效性.对激光等离子体相互作用进行初步模拟分析,表明通过采用球坐标下的动量空间,可在相对较少动量网格情形下,获得与粒子模拟可相互验证的结果.     

交错网格下MHD相容守恒格式的发展

在低磁场雷诺数条件下,基于电势泊松方程,发展了交错网格下可以精确计算电流和洛伦兹力(电磁力)的相容守恒格式。采用压力为变量的原始变量法求解不可压缩Navier-Stokes方程,所计算的电流满足电荷守恒定律,所计算的电磁力满足动量守恒定律。对金属流体在Hartmann数50~5000范围内验证了格式的精确性。交错网格下相容守恒格式的发展为后续MHD稳定性分析、湍流的大涡模拟及直接数值模拟提供很好的选择。     

基于特征理论的二维无粘Lagrange流体力学有限体积法

提出一个求解二维无粘Lagrange流体力学方程的中心型有限体积方法.采用特征理论求解网格节点处的速度及压力,并利用这些物理量更新节点位置及计算网格界面通量.方法适用于结构网格与非结构网格.典型数值实验的结果表明,格式具有较好的收敛性、对称性、能量守恒性及鲁棒性,且能自然地求解多物质流动问题.     

移动网格与Level Set耦合方法在气-液两相流数值模拟中的应用

本文采用自适应移动网格与Level Set函数相耦合的方法来实现气-液两相流的数值模拟与计算.作为自适应网格方法的一种,移动网格方法主要是为了解决发展方程的计算问题而设计的方法.文中给出了移动网格的生成方程,并针对方程的非线性,给出了一种半隐式的离散方法用于进行求解.本文将移动网格方法与Level Set方法相耦合,将控制流体运动的Navier-Stokes方程以及追踪相界面的Level Set方程转换到曲线坐标下,应用一套曲线坐标方程组来同时描述气、液两相流的运动规律,成功实现了对气-液两相流问题的数值模拟.通过对顶盖驱动流的计算以及对液滴沉降现象的模拟计算,验证了本文方法的可靠性.本文对常重力与微重力下两气泡融合的发展规律进行了数值模拟,通过分析对比,得到了重力对两气泡融合变形的影响规律.     

大雷诺数Navier-Stokes方程的两水平亚格子模型稳定化方法

基于两重网格离散方法,提出三种求解大雷诺数定常Navier-Stokes方程的两水平亚格子模型稳定化有限元算法.其基本思想是首先在一粗网格上求解带有亚格子模型稳定项的Navier-Stokes方程,然后在细网格上分别用三种不同的校正格式求解一个亚格子模型稳定化的线性问题,以校正粗网格解.通过适当的稳定化参数和粗细网格尺寸的选取,这些算法能取得最优渐近收敛阶的有限元解.最后,用数值模拟验证三种算法的有效性.     

润滑油对风冷冷凝器管内换热和阻力影响的分析

文中抛弃传统的换热器两相区简单均相流模型 ,采用三相双流体模型 ,在合理假设的基础上 ,针对风冷冷凝器可能出现的流型 ,分别建立相应的控制方程。在保证质量、动量和能量守恒基础上 ,考虑相应气液界面关系 ,联立求解。并以 R1 2为矿物油、R1 34a与 P.A.G.油为例 ,分析不同含油率情况下对换热和流动的影响     

II. Recirculation flow induced by a bubble stream rising in a viscous liquid

The recirculation flow induced by the rising motion of a bubble stream in a viscous fluid within an open-top rectangular enclosure is studied. The three-dimensional volume averaged conservation equations are solved by a control-volume method using a hybrid finite differencing scheme to describe the liquid phase hydrodynamics. The momentum exhange between the bubbles and the liquid phase is modeled with a source term equals to the volumetric buoyancy force acting on the gas in the bubble stream. The volumetric buoyancy force accounts for in line interactions between bubbles through the average gas volume fraction in the gas liquid column which depends on the size and the rising velocity of bubbles. The fluid flow within an open-top rectangular enclosure is further investigated by particle image velocimetry for a bubble stream rising in a water-glycerol solution. The measured fluid velocities in a vertical plane are compared with the predictions of the numerical model over a wide range of fluid viscosity (43 mPa s-800 mPa s) and gas flow rates. Finally, the recirculation flows resulting from the interaction of two neighbouring vertical bubble streams are studied. Received: 23 July 1997 / Revised: 19 December 1997 / Accepted: 11 May 1998     

Numerical simulations of gas-liquid two-phase flows in a micro porous structure

The lattice Boltzmann method for two-phase fluid flows is applied to the simulations of gas-liquid two-phase flows in a micro porous structure for various capillary numbers at low Reynolds numbers. The behaviors of the gas-liquid interface and the velocities of the two-phase fluid in the structure are simulated, and the permeability of gas and liquid through the structure are estimated from the calculated results. By changing the void fraction, the contact angle of the interface on walls, and the surface tension, the effect of these properties on the behaviors and the permeability of the two-phase flows in the micro porous structure is investigated. It is found that the permeability of liquid flows depends on the contact angle and it increases for hydrophobic walls. It is also seen that liquid flows are choked in pores for large void fractions and low capillary numbers.     

三维非均匀不稳定渗流方程的自适应网格粗化算法

将渗透率自适应网格技术应用于三维非均匀不稳定渗流方程的网格粗化算法中,在渗透率或孔隙度变化异常区域自动采用精细网格,用直接解法求解渗透率或孔隙度变化异常区域的压强分布,在其它区域采用不均匀网格粗化的方法计算,即在流体流速大的区域采用精细网格.用该方法计算了三维非均匀不稳定渗流场的压降解,结果表明三维非均匀不稳定渗流方程的三维非均匀自适应网格粗化算法的解在渗透率或孔隙度异常区的压强分布规律与采用精细网格的解非常逼近,在其它区域压强分布规律与粗化算法的解非常逼近,计算速度比采用精细网格提高100多倍.     

A transient method for measuring the gas volume fraction in a mixed gas-liquid flow using acoustic resonance spectroscopy

In this paper, the feasibility of measuring the gas volume fraction in a mixed gas-liquid flow by using an acoustic resonant spectroscopy (ARS) method in a transient way is studied theoretically and experimentally. Firstly, the effects of sizes and locations of a single air bubble in a cylindrical cavity with two open ends on resonant frequencies are investigated numerically. Then, a transient measurement system for ARS is established, and the trends of the resonant frequencies (RFs) and resonant amplitudes (RAs) in the cylindrical cavity with gas flux inside are investigated experimentally. The measurement results by the proposed transient method are compared with those by steady-state ones and numerical ones. The numerical results show that the RFs of the cavity are highly sensitive to the volume of the single air bubble. A tiny bubble volume perturbation may cause a prominent RF shift even though the volume of the air bubble is smaller than 0.1% of that of the cavity. When the small air bubble moves, the RF shift will change and reach its maximum value as it is located at the middle of the cavity. As the gas volume fraction of the two-phase flow is low, both the RFs and RAs from the measurement results decrease dramatically with the increasing gas volume, and this decreasing trend gradually becomes even as the gas volume fraction increases further. These experimental results agree with the theoretical ones qualitatively. In addition, the transient method for ARS is more suitable for measuring the gas volume fraction with randomness and instantaneity than the steady-state one, because the latter could not reflect the random and instant characteristics of the mixed fluid due to the time consumption for frequency sweeping. This study will play a very important role in the quantitative measurement of the gas volume fraction of multiphase flows.     

黏弹性流体充模过程中凝固现象的数值模拟

建立了黏弹性流体在充模过程中带有相变的气一液两相模型,该模型分别由气、液两相的质量守恒方程、动量守恒方程、能量守恒方程描述,并通过引入Heaviside函数将气一液两相的方程组统一为一个方程组;建立了一个对型腔内熔体和气体都适用的修正的焓方法来描述充模过程中的相变,采用基于同位网格的有限体积方法对模型进行求解,水平集方法捕捉充模过程中的界面演化,模拟出了黏弹性流体在充模过程中的凝固现象,得出了充模过程中型腔内的温度、压力、第一法向应力差等随时间的变化;并讨论了型腔壁面温度、熔体温度、注射速度对充模过程中凝固现象的影响,研究结果表明:型腔壁面温度越高,凝固层越薄;熔体温度越高,凝固层越薄;注射速度越高,凝固层越薄,故提高型腔壁面温度、熔体温度、注射速度可以减少或消除型腔壁面附近的凝固层。     

颗粒增强湍流的新模型和气粒流动的湍流变动

用雷诺应力方程模型和极细的网格系对单个颗粒受湍流气体绕流进行了数值模拟,研究了改变颗粒直径和气体相对速度时颗粒增强气体湍流的规律.据此构造了颗粒尾涡增强气体湍流的新模型.将此子模型加入到两相流动模型中,对竖直和水平通道内气粒两相流动进行了数值模拟,和实验结果的对照表明,考虑颗粒尾涡增强气体湍流效应得到的气体湍流脉动速度的模拟结果比不考虑此效应的模拟结果好得多.     

An effective explicit pressure gradient scheme implemented in the two-level non-staggered grids for incompressible Navier–Stokes equations

In this paper, an improved two-level method is presented for effectively solving the incompressible Navier–Stokes equations. This proposed method solves a smaller system of nonlinear Navier–Stokes equations on the coarse mesh and needs to solve the Oseen-type linearized equations of motion only once on the fine mesh level. Within the proposed two-level framework, a prolongation operator, which is required to linearize the convective terms at the fine mesh level using the convergent Navier–Stokes solutions computed at the coarse mesh level, is rigorously derived to increase the prediction accuracy. This indispensable prolongation operator can properly communicate the flow velocities between the two mesh levels because it is locally analytic. Solution convergence can therefore be accelerated. For the sake of numerical accuracy, momentum equations are discretized by employing the general solution for the two-dimensional convection–diffusion–reaction model equation. The convective instability problem can be simultaneously eliminated thanks to the proper treatment of convective terms. The converged solution is, thus, very high in accuracy as well as in yielding a quadratic spatial rate of convergence. For the sake of programming simplicity and computational efficiency, pressure gradient terms are rigorously discretized within the explicit framework in the non-staggered grid system. The proposed analytical prolongation operator for the mapping of solutions from the coarse to fine meshes and the explicit pressure gradient discretization scheme, which accommodates the dispersion-relation-preserving property, have been both rigorously justified from the predicted Navier–Stokes solutions.     

集输管路上升管系统严重段塞流实验和理论模拟

严重段塞流的实验研究表明,在气泡进入上升管底部到运动至出口的过程中,上升管中气泡头部以下流型为弹状流型;当气泡头部流出上升管后,上升管中的流型可看作块状流型。根据实验结果,本文提出了采用漂移流模型简化计算上升管中两相流动、上游管道中气体膨胀满足质量守恒,同时考虑上升管内液体动量守恒的严重段塞流计算模型。计算值与测量值比较表明,模型可以正确预测出气体膨胀流动过程,气体流动时间不受入口气液流量的影响。模型可以准确计算出严重段塞流周期、液塞长度和倾斜管中液柱最大长度等参数。     

Calculation of shock wave propagation in water containing reactive gas bubbles

The entry of a shock wave from air into water containing reactive gas (stoichiometric acetylene–oxygen mixture) bubbles uniformly distributed over the volume of the liquid has been numerically investigated using equations describing two-phase compressible viscous reactive flow. It has been demonstrated that a steady-state supersonic self-sustaining reaction front with rapid and complete fuel burnout in the leading shock wave can propagate in this bubbly medium. This reaction front can be treated as a detonation-like front or “bubble detonation.” The calculated and measured velocities of the bubble detonation wave have been compared at initial gas volume fraction of 2 to 6%. The observed and calculated data are in satisfactory qualitative and quantitative agreement. The structure of the bubble detonation wave has been numerically studied. In this wave, the gas volume fraction behind the leading front is approximately 3–4 times higher than in the pressure wave that propagates in water with air bubbles when the other initial conditions are the same. The bubble detonation wave can form after the penetration of the shock wave to a small depth (~300 mm) into the column of the bubbly medium. The model suggested here can be used to find optimum conditions for maximizing the efficiency of momentum transfer from the pressure wave to the bubbly medium in promising hydrojet pulse detonation engines.