粘性不可压流体流动问题用直角坐标网格的贴体解法

研究一种新的全贴体的求解粘性不可压流体流动问题的非结构化直角坐标网格方法.该方法在于利用直角坐标网格但通过在边界附近保留不规则控制体,使得算法是完全贴体的.这有别于目前流行的各种非结构化直角坐标网格方法.通过对两个典型流动问题的计算对该数值方法进行验证.对比结果表明,本方法计算的结果与精确解和STAR-CD的结果在一定Re数和网格数时是很接近的,可以满足一定的精度要求,说明该数值计算方法是可行的.还对二维钝头体周围的流场进行了计算,计算的流场与STAR-CD的结果相当吻和,说明该算法还可计算较复杂的流动现象.     

基于黎曼解的一种SPH改进算法

将基于黎曼解的粒子间接触算法与基于Taylor展开的插值方法相结合,构造出一种改进的光滑粒子算法,给出了两个算例.数值计算表明,改进算法可以提高自由面计算精度、克服拉应力不稳定性,避免了CSPM在核函数的影响域包含较多粒子且流场梯度很大时出现的压缩失稳现象.     

二维管道充填过程的修正SPH模拟

通过施加一种密度初始化方法对传统光滑粒子动力学(SPH)方法进行修正,提出一种修正SPH方法.同时,为了提高边界上数值计算的准确性,提出一种新的固壁边界处理方法.通过修正SPH方法模拟液滴拉伸问题和溃坝问题,验证修正SPH方法的准确性和可靠性.随后,对研究很少的管道充填过程进行修正SPH模拟,并讨论Re对流场及涡的影响.数值结果表明,修正SPH方法能够准确模拟牛顿流体管道充填过程,且流动受Re的影响较大.     

基于光滑粒子动力学方法的液滴冲击固壁面问题数值模拟

采用改进的光滑粒子动力学(SPH)方法对液滴冲击固壁面问题进行了数值模拟. 为了提高传统SPH方法的计算精度和数值稳定性, 在传统的SPH方法的基础上对粒子方法中的密度和核梯度进行了修正, 采用了考虑黎曼解法的SPH流体控制方程, 构造了一种新型的粒子间相互作用力(IIF)模型来模拟表面张力的影响. 应用改进的SPH方法对液滴冲击固壁面问题进行了数值模拟. 计算结果表明：新型的IIF 模型能够较好地模拟表面张力的影响, 改进的SPH方法能够精细地描述液滴与固壁面相互作用过程中液滴的内部压力场演变和自由面形态变化, 液滴的铺展因子随初始韦伯数的增大而增大, 数值模拟结果与实验得到的结果基本一致. 关键词： 液滴 固壁面 光滑粒子动力学 表面张力     

K2_SPH方法及其对二维非线性水波的模拟

K2_SPH方法通过泰勒级数展开和联立求解积分方程组的办法得到具有二阶精度的核近似方法.随着核近似精度的提高,K2_SPH需要对-些关键数值技术进行改进才能成功模拟非线性水波问题,例如自由表面边界和固壁边界.通过与传统SPH方法计算结果比较,K2_SPH方法在非线性自由表面计算精度和整个粒子系统中有关变量分布都有显著提高.     

eXtended Pom-Pom黏弹性流体的改进光滑粒子动力学模拟

黏弹性流体广泛存在于自然界和工业生产中,对其复杂流变特性的研究具有重要的学术价值和应用意义.本文提出一种改进的光滑粒子动力学方法,对基于eXtended Pom-Pom模型的黏弹性流动进行了数值模拟.为了提高计算精度,采用一种不含核导数计算的核梯度修正离散格式.为了防止粒子穿透固壁,提出一种增强型的边界处理技术.为了消除张力不稳定性,将人工应力耦合到动量守恒方程中.运用改进光滑粒子动力学方法数值模拟了基于eXtended Pom-Pom模型的黏弹性Poiseuille流和黏弹性液滴撞击固壁问题,通过与解析解或有限差分方法解的比较以及对数值收敛性的评价,验证了改进光滑粒子动力学方法的有效性和优势,并在此基础上,深入分析了Reyonlds数、Weissenberg数、溶剂黏度比、各向异性参数、松弛时间比和分子链臂数等流变参数对流动过程的影响.     

液面冲击引起液滴飞溅问题的CIP方法数值模拟

利用自主研发的CIP-ZJU（Constrained Interpolation Profile Method in Zhejiang University）高精度数学模型研究强非线性自由表面流动问题.模型在直角坐标系统下建立,采用紧致插值曲线CIP方法作为流场的基本求解器,通过多相流的方式实现固-液-气耦合同步求解,采用平衡格式的VOF（VOF/WLIC：Volume of Fluid/Weighted LineInterface Calculation）自由面捕捉方法改进了原模型,利用浸入边界方法处理运动物体.利用改进的CIP模型开展了不同类型的物体冲击液面引起的液滴飞溅现象的数值模拟,重点分析液滴飞溅过程中的自由液面变形、作用荷载和物体位置等.通过数值结果与实验结果的比较验证模型的可靠性.结果表明：平衡格式的VOF自由面捕捉方法能更精确地重构自由面,本文的数学模型可精确预测强非线性自由表面流动问题.     

一种基于预估-修正策略的非定常流网格自适应加密/稀疏方法

给出一种非定常流动数值模拟的网格自适应处理方法.在"求解流动方程-自适应调整网格"的流程中,引入预估-修正步.根据自适应周期内每个时间步上的流场预估解,计算单元上的事后误差估算值.建立考虑解演变的网格自适应指示器,并进行多层次单元加密-稀疏的动态网格自适应处理.在自适应网格上重新计算流场.每个自适应周期中,流动演变区域的网格获得加密;而前一个周期中的特征现象已离开区域的网格被稀疏.应用边界非协调的当地DFD(Domain-Free Discretization)方法求解流动方程.为验证网格自适应处理方法,针对静止圆柱和自推进游鱼的流动进行了数值实验.     

一、二维弹塑性流体动力学数值计算的高精度新自由面格式

提出了适用于X射线或电子束辐照、爆炸和高速碰撞等问题的一、二维弹塑性流体动力学数值模拟计算的一种新自由面格式。对于X射线及电子束辐照材料与结构的热-力学效应的数值计算结果表明,该新自由面格式明显优于目前一直沿用的Richtmyer格式。从理论上证明了该新自由面格式具有三阶精度,而Richtmyer格式为一阶精度。     

高精度格式WCNS-E-5计算物面热流

采用加权紧致非线性格式WCNS-E-5,四阶精度的二阶导数差分近似以及四阶精度的边界格式构造了高精度算法,对高超声速粘性流动的物面热流进行数值研究.首先考察了壁面网格雷诺数对驻点热流的影响,然后开展了边界格式对热流计算结果影响的研究,最后对大攻角钝锥绕流进行了数值模拟.研究结果表明:WCNS-E-5能降低边界层内网格分辨率,全场高精度的WCNS-E-5计算得到的流场图像清晰、真实、分辨率高,热流值准确、可靠.     

On nonlinear periodic capillary-fluctuation wave flow in a thin liquid film on a solid substrate

Analytical calculation of a nonlinear periodic wave flow on the free surface of a charged layer of an ideal incompressible conducting liquid resting on a solid substrate is carried out for the case when fluctuation-induced forces (the dispersion component of the wedging pressure) have a decisive effect on the system. It is shown that wave flows emerge in the liquid in calculations of the second order of smallness in the wave amplitude, which is assumed to be small compared with the thickness of the liquid layer. These flows result from nonlinear interaction as nonlinear corrections to the waves set at the zero time. The field of fluctuation-induced forces displaces these flows toward the periphery of the area of influence of these forces. This effect takes place both in the presence of an external electric field near the free surface and in its absence. The sign and value of the nonlinear corrections depend on whether an electric field is present near the free surface of the liquid. In the presence of an electric field, the curvature of the crest of the nonlinear waves increases; in its absence, the curvature decreases.     

几种液体金属自由表面的MHD不稳定性实验研究

介绍了在核工业西南物理研究院的液态金属实验回路(LMEL)上获得的几种可供液态偏滤器-限制器系统选用的液体自由表面的磁流体动力学(MHD)效应不稳定性的实验结果。实验发现：自由表面射流在穿越梯度横磁场时射程变短、截面变扁平,但MHD效应稳定,调节射流与磁场的夹角可以控制射流的流动特性;自由表面膜流MHD效应存在三种现象,即层流、溪状流和湍流。层流是由多束射流打到固体表面产生的(简称“射-膜流”),从MHD效应角度考虑,“射-膜流”将是四种可选液态偏滤器-限制器系统的液体自由表面形式中最佳的选择。同时,探讨了从物理的角度来理解四种自由表面形式的MHD效应的现象。     

Boundary data immersion method for Cartesian-grid simulations of fluid-body interaction problems

A new robust and accurate Cartesian-grid treatment for the immersion of solid bodies within a fluid with general boundary conditions is described. The new approach, the Boundary Data Immersion Method (BDIM), is derived based on a general integration kernel formulation which allows the field equations of each domain and the interfacial conditions to be combined analytically. The resulting governing equation for the complete domain preserves the behavior of the original system in an efficient Cartesian-grid method, including stable and accurate pressure values on the solid boundary. The kernel formulation allows a detailed analysis of the method, and it is demonstrated that BDIM is consistent, obtains second-order convergence relative to the kernel width, and is robust with respect to the grid and boundary alignment. Formulation for no-slip and free slip boundary conditions are derived and numerical results are obtained for the flow past a cylinder and the impact of blunt bodies through a free surface. The BDIM predictions are compared to analytic, experimental and previous numerical results confirming the properties, efficiency and efficacy of this new boundary treatment for Cartesian grid methods.     

Һ̬�������ɱ���������������ѧЧӦ����ֵģ��

分析了磁流体力学效应对液态金属自由表面射流稳定性的影响。从射流的感应电势、电流、速度等方面,解释射流在磁场中稳定的原因。数值计算结果验证了理论分析。     

On energy transfer and generation of an electric current in the vicinity of an electrode dipped in an electrolyte and strongly heated by a current passing through it

Processes occurring when a metal electrode dipped in an electrolyte is heated by intense evaporation of the electrolyte are considered in terms of a physically rigorous model. Based on the Onsager principle of least energy dissipation rate in nonequilibrium processes, the fractions of thermal energy that are spent on heating and evaporating the electrolyte and on heating the vapor are found. The energy is released within the vapor-gas sheath when an electric current flows between the electrode and electrolyte surface. It is found that the electrolyte vapor temperature exceeds 1300 K. Analytical expressions are derived for the vapor-gas sheath thickness, the electrolyte vapor pressure, and the velocity of the vapor escaping the discharge zone. It is shown that field evaporation of thermally activated negative ions from the electrolyte surface cannot provide an electric current with densities found in experiments but is responsible for the generation of free electrons near the electrolyte surface. These electrons arise when the ions decay via collisions with excited molecules.     

Incompressible smoothed particle hydrodynamics for free-surface flows: A generalised diffusion-based algorithm for stability and validations for impulsive flows and propagating waves

The incompressible smoothed particle hydrodynamics (ISPH) method with projection-based pressure correction has been shown to be highly accurate and stable for internal flows and, importantly for many problems, the pressure field is virtually noise-free in contrast to the weakly compressible SPH approach (Xu et al., 2009 [31]). However for almost inviscid fluids instabilities at the free surface occur due to errors associated with the truncated kernels. A new algorithm is presented which remedies this issue, giving stable and accurate solutions to both internal and free-surface flows. Generalising the particle shifting approach of Xu et al. (2009) [31], the algorithm is based upon Fick’s law of diffusion and shifts particles in a manner that prevents highly anisotropic distributions and the onset of numerical instability. The algorithm is validated against analytical solutions for an internal flow at higher Reynolds numbers than previously, the flow due to an impulsively started plate and highly accurate solutions for wet bed dam break problems at zero and small times. The method is then validated for progressive regular waves with paddle motion defined by linear theory. The accurate predictions demonstrate the effectiveness of the algorithm in stabilising solutions and minimising the surface instabilities generated by the inevitable errors associated with truncated kernels. The test cases are thought to provide a more thorough quantitative validation than previously undertaken.     

Experimental-calculation simulation of the ejection of particles from a shock-loaded surface

The possibilities of measuring complex Pylenie, which is intended for studying the parameters of the particles ejected from the free surface of a shock-loaded material and was created at the Institute of Experimental Gas Dynamics and Detonation Physics RFNC-VNIIEF, are demonstrated. The operation of the complex is based on the following three methods, which are based on different physical principles and supplement each other: laser-optical method (macro- and microfilming), pulsed X-ray method, and piezoelectric pressure sensors. This complex is used to study the ejection of particles from the free surface of lead samples when a shock wave with a pressure of 7, 16, and 23 GPa at its front reaches this surface. The effect of the surface roughness and the pressure amplitude at the shock wave front on the quantitative characteristics of the process is shown. The calculation-theoretical simulation is performed by two-dimensional numerical calculations of gasdynamic flows, and the calculation results are used to estimate the “ejecta” characteristics in terms of the developed phenomenological model of the process.     

The influence of solid bodies on low Mach number vortex sound

The mechanisms of sound generation by unsteady, subsonic flows in the presence of solid boundaries are investigated. For this purpose an alternative integral representation for the radiated pressure field is applied which is different from the generally used integral representation introduced by Lighthill and Curle. The main advantage of the method is that there is a linear dependence of the integrand on the time derivative of the vorticity fluctuations in the hydrodynamic near field; instead of the ordinary Green function a “vector Green function” is used. This vector Green function can be chosen for a given flow field in such a way that surface integrals do not appear. Finally, the theory is illustrated by two- and three-dimensional model flows. Analytical solutions are determined by applying the method of matched asymptotic expansions.