近似保角的自适应坐标变换方法

自适应坐标变换方法是为解决多介质和大变形问题而提出的一类网格生成方法，该方法中的一种为近似保持网格夹角不变，保持物质界面为拉氏描述，并要求网格速度在最小二乘意义下尽量靠近流体运动速度。这里所讨论的坐标变换的自适应性，指的是新坐标系自动适应流体流场的一些重要特性（接近流体速度）以及保持网格的几何特性（保角）。为了处理多介质情况，网格方程应在子区域的所有边界上给出边界条件。     

强爆炸火球辐射流体自适应网格高精度数值模拟

提出一种自适应网格方法,应用于基于Euler方法的强爆炸辐射流体高精度数值求解。通过与Zinn数值结果对比,验证该方法的正确性。研究自适应网格对冲击波和光辐射输出模拟精度的影响,对比不同网格尺度下的计算耗时。在相同的条件下,使用自适应网格与均匀网格加密3倍得到的冲击波超压分布、光辐射输出演化接近,计算效率提升约8.5倍。由此可知,提出的自适应网格方法可用于强爆炸问题的高精度数值模拟,显著地提升模拟精度和模拟效率。     

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

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

矩形网格上VOF运动界面重构的流体体积分数保持法

针对处理运动界面问题的流体体积函数(VOF)法,提出一种新的流体体积分数保持界面重构算法.该方法在单个网格内用斜线段近似运动界面,并要求相邻网格内的斜线段在公共边上的交点重合,通过保持网格内的流体体积分数不变,在边界网格上建立非线性方程组,通过求解非线性方程组,确定运动界面边界线与网格线的交点位置,最终重构出运动界面的形状.应用算例模拟结果表明,该方法能够高精度,较少尖点地重构运动界面,证实了该方法的有效性和可行性.     

振动NACA0012翼型的移动网格数值模拟

利用流体力学方程的积分形式给出非结构移动网格上离散格式,利用自适应移动网格方法移动网格,进而得到网格速度.对振动Naca0012翼型问题,分三种类型确定网格速度,再结合Riemann问题的解法器构造数值通量,得到移动网格单元上新的物理量.数值实验表明这种格式同时具有高效、高分辨率的特点.     

基于WENO重构的高阶移动网格动理学格式

提出一种基于WENO重构的高阶(至少三阶)移动网格动理学格式.利用流体力学方程的积分形式得到移动网格上离散格式,再利用自适应移动网格方法移动网格,进而得到网格速度,利用WENO重构得到高阶插值多项式,最后使用时间方向上精确的动理学数值方法构造数值通量,得到移动网格单元上新的物理量.数值实验表明这种格式同时具有高精度、高分辨率的特点.     

纳米流体多相流动的多尺度模拟方法

针对纳米流体多相流动的微观特征,提出一种基于格子Boltzmann的多尺度耦合方法,在速度和悬浮纳米粒子分布变化比较剧烈的区域采用细网格多相模型,在其它区域视纳米流体为均匀混合的单相流体,使用粗网格单相模型.为保证不同尺度区域之间的物理信息(参数)的准确传递,运用质量和动量守恒原理,建立跨区域的边界耦合模型.几个算例表明,该方法既可以反映纳米流体流动的微观特征,又能提高计算效率,与单纯使用多相模型相比,节省大量时间.     

自适应间断有限元方法求解双曲守恒律方程

研究自适应Runge-Kutta间断Galerkin (RKDG)方法求解双曲守恒律方程组,并提出两种生成相容三角形网格的自适应算法.第一种算法适用于规则网格,实现简单、计算速度快.第二种算法基于非结构网格,设计一类基于间断界面的自适应网格加密策略,方法灵活高效.两种方法都具有令人满意的计算效果,而且降低了RKDG的计算量.     

用自适应流线网格计算定常超声速流场

提出一个用自适应流线网格计算定常超声速流场的数值方法。Friedrichs-Lax格式推广到任意四边形网格,在迭代过程中得到速度场后,从流线方程计算流线的位置,最后算出由本方法得到高分辨率结果     

二维自适应结构网格的变分生成方法

提出-一种二维结构网格的自适应生成方法,即采用变分方程控制网格的疏密程度以及光滑性、正交性和正规性,并且所有的控制方程具有-致的量纲.几个数值算例表明,给出的网格自适应生成方法能够较好地收敛到方程的唯-解,同时产生性质良好的二维自适应结构网格.     

入射速度对长杆弹垂直侵彻行为的影响规律

 以长杆弹垂直侵彻半无限厚靶板为研究对象,分析了弹体最大侵彻深度与入射速度的关系,研究了弹体入射速度对侵彻最大深度的影响规律。研究表明：靶板的强度和界面效应使弹体在侵彻过程中存在一个临界速度,当入射速度大于临界速度时,弹体的侵彻才能通过开坑阶段进入准稳定阶段,它是造成当入射速度较小时侵彻深度随入射速度的提高而几乎不变或缓慢增加的主要原因;准稳定侵彻过程中弹体速度和侵彻速度基本不变,并且两者存在线性关系,这种关系只与弹体和靶板的材料性能有关,是造成当入射速度较大时侵彻深度随入射速度的提高呈快速线性增大的主要原因。     

Interaction of two parallel plane jets of different velocities

The interaction of two parallel plane jets of different velocities is studied by flow visualization and PIV measurement to examine the influence of velocity ratio on the development ofjets in the initial region. It is found that the parallel plane jets develop toward the high velocity side and the jet width is reduced with a decrease in the jet velocity ratio. Corresponding to the variation of mean velocity field to the velocity ratio, the magnitudes of turbulence intensities, Reynolds stress and static pressure are weakened in the merging region of the jets and their peak locations of the properties are shifted to the high velocity side. These results indicate that the interaction of two parallel jets is weakened with a decrease in the velocity ratio of the jets.     

组合时延估计及人工时延方法用于脉冲相干多普勒测速去模糊

提出了一种改进的脉冲相干多普勒测速的速度模糊消除方法,将待测速度空间以两倍最大可测速度的间隔划分为一系列速度区间,使得对任意速度的测量转化为对速度区间的确定和区间内速度偏移量的测量,通过"相关-相位"组合法进行精确时延估计,实现无模糊速度测量;提出"人工时延"的方法,消除在速度区间边界值附近可能存在的速度跳变;与传统的速度解模糊方法相比,该方法仅需要对多次人工时延的方案进行计算,无需多次测量,且理论上不受流速范围的限制;通过仿真和实验对该方法的有效性进行了验证。      

Molecular dynamics simulations of the nano-droplet impact process on hydrophobic surfaces

Large-scale molecular dynamics simulations are used to study the dynamic processes of a nano-droplet impacting on hydrophobic surfaces at a microscopic level. Both the impact phenomena and the velocity distributions are recorded and analyzed. According to the simulation results, similar phenomena are obtained to those in macro-experiments. Impact velocity affects the spread process to a greater degree than at a level of contact angle when the velocity is relatively high. The velocity distribution along the X axis during spread is wave-like, either W- or M-shaped, and the velocity at each point is oscillatory; while the edges have the highest spread velocity and there are crests in the distribution curve which shift toward the edges over time. The distribution along the Y axis is 〈- or 〉-shaped, and the segments above the middle have the lowest decrease rate in the spreading process and the highest increase rate in the retraction process.     

On the dynamics of the clavichord: from tangent motion to sound

An experimental study of variations in the sound of clavichord notes at different dynamic levels is described. Radiated acoustic signal, tangent velocity and two tangent-string contact signals are synchronously measured for all 51 notes of an unfretted instrument. More than ten repeated measures are recorded in order to obtain as much variation in dynamic level as possible. The tangent motion, expressed in terms of velocity, is studied in the time and frequency domains. A model of the tangent-string contact point velocity is proposed. Then, three aspects of the sounded tones are analyzed: SPL and its relationship to tangent velocity, spectral slope, and pitch variations. These results indicate a linear relationship between sound pressure level and tangent peak log velocity. Spectral slope seems almost constant independent of tangent velocity and dynamic level. Both tangent velocity and finger pressure are shown to influence the fundamental frequency. In conclusion, controlling both finger velocity and finger pressure may prove challenging for the player, and this may explain why the sound quality of the clavichord depends so much on the players ability.     

甲烷/氧气层流反扩散火焰形态及滞后特性研究

对空气气氛中甲烷/氧气反扩散火焰的形态和推举滞后特性进行了实验研究. 实验中通过改变气体流量考察了气速变化对火焰形态演变及滞后特性的影响, 并利用紫外相机系统研究了气速对不同形态火焰中OH^*分布的影响. 研究结果表明: 甲烷气速、氧气气速和火焰的历史状态是决定火焰形态的三个重要参数, 并以此对实验范围内的火焰形态进行了分区; 氧气气速对不同形态反扩散火焰轴线上的OH^*分布有相似的影响, 当氧气缺乏时, 反扩散反应区较短, 当氧气富余时, 反扩散反应区在轴向分布较广; 同轴甲烷的气速对反扩散火焰的滞后特性影响显著, 随着甲烷气速的增加, 反扩散火焰的推举速度和再附着速度呈线性减小, 部分预混火焰向反扩散火焰转变的速度呈线性增加.     

激光多普勒自主测速仪天线安装误差对测速精度的影响

为了减小激光多普勒自主测速仪天线安装误差引起的测速误差,以多普勒频移原理为基础,通过旋转矩阵的复合变换推导出测速误差与安装误差间的函数关系,分析了三个安装误差角对三维测速精度的影响,并对旋转矩阵采用Bursa模型线性化后在整体最小二乘准则下给出安装误差角的最优估计,从而实现测速误差的补偿。仿真结果表明:测速误差随速度的增大而增大,影响纵向、横向测速精度的主要因素是偏航角,影响垂向测速精度的主要因素是纵倾角,补偿后三维测速误差显著减小。     

Numerical simulation of transit-time ultrasonic flowmeters: uncertainties due to flow profile and fluid turbulence

Flowmeter measurement using the ultrasonic transit-time method is based on the apparent difference of the sound velocity in the flow direction and in the opposite direction. This method gives a flow velocity averaged along a particular acoustical path. To convert this path velocity to a velocity averaged over the entire cross-section of the flowing medium, the knowledge of the flow velocity profile is essential. However, the acoustical paths joining the two transducers are supposed to be straight and fluid turbulence phenomena are neglected. In this paper, we describe a numerical procedure to estimate the uncertainties due to these approximations in the case of fully developed turbulence. The ultrasonic propagation is modelled in 2-D moving inhomogeneous media via a ray tracing algorithm. Influence of mean profiles of temperature and velocity is studied on simple examples. Fluid temperature fluctuations and fluid velocity turbulence are considered in the stochastic framework to obtain average uncertainties on the measurements of the liquid flow rate.     

A new method of ultrasound color flow mapping

Conventional ultrasound color flow mapping systems estimate and visualize only the axial velocity component. To obtain the transverse velocity component a modification of a multiple-beam method is proposed. The new two-dimensional color flow mapping system has a small size and consists of three transducers. The central transducer is an appodized and focused phased array. The other transducers are unfocused probes. Three transducers act as receivers and the central transducer operates as a transmitter. All receivers acquire rf scan lines that are then processed to estimate three axial velocity components using an autocorrelation method. These estimates are then combined to estimate the transverse velocity component, taking into account the geometric relationships among three transducers. Two algorithms for transverse velocity estimation are proposed. The first uses the Doppler angle estimate for calculation of the transverse velocity component. The other algorithm calculates the transverse velocity component directly from the axial components. The accuracy of the flow velocity estimators is estimated by simulations. Analysis of accuracy allows choosing the more effective algorithm for two-dimensional velocity estimation, which is insensitive to variations of the Doppler angle.