The sound velocity in ideal liquid mixtures from thermal volume fluctuations.

The selection of the correct mixing rule for sound velocity in ideal liquid mixtures determines the interpretation of the sound velocity in real mixtures. This is especially important for the determination of apparent properties of solutes, such as their apparent compressibility. There are different approaches reported in the literature, and this article presents a new derivation of the mixing rule based on statistical mechanics. It is shown that the correlation of volume fluctuations between adjacent components has a crucial influence on the ideal mixing rule.     

论可压缩平板边界层线性稳定性的分歧

Ｍａｃｋ和Ｗａｚｚａｎ关于可压缩平板边界层线性稳定性结论的主要分歧在于来流马赫数对粘性稳定性的影响．本文充分考虑了空气热力学参数的影响，并用配置点方法计算了绝热平板时间模式特征值问题．数值结果表明粘性对第一模式只起稳定作用，Ｍ＝３粘性第二模式总是稳定的，结论和Ｍａｃｋ一致     

高速液—固碰撞时固体弹性和液体可压缩性对撞击压力的影响

本文联立求解可压缩平面液体撞击弹性固体表面时液体激波面及液固相界面上的物理过程控制方程组,导出了激波速度、撞击压力和固体表面变形速度的计算公式.导出的公式中包括了表征液体可压缩性的撞击马赫数M_0和表征固体弹性的液固声阻抗比Γ.通过计算液体撞击PMMA,钙钠玻璃、ZnSe和1Cr13时的撞击压力,本文具体分析了M_0和Γ对撞击压力产生的综合影响。     

空化可压缩流动空穴溃灭激波特性研究

为深入研究空化可压缩流动中空泡/空泡团溃灭过程中激波产生、传播及其与空穴相互作用规律,本文采用数值模拟方法对空化可压缩流动空穴溃灭激波特性展开了研究.数值计算基于OpenFOAM开源程序,综合考虑蒸汽相和液相的压缩性,通过在原无相变两相可压缩求解器的控制方程中耦合模拟空化汽液相间质量交换的源项,实现了对空化流动的非定常可压缩计算.利用上述考虑汽/液相可压缩性的空化流动求解器,对周期性云状空化流动进行了数值模拟,并重点研究了空穴溃灭激波特性.结果表明:上述数值计算方法可以准确捕捉到空穴非定常演化过程及大尺度脱落空泡云团溃灭激波现象,大尺度脱落空泡云团溃灭过程分为3个阶段:(1) U型空泡团形成; (2) U型空泡团头部溃灭; (3) U型空泡团腿部溃灭.在U 型空泡团腿部溃灭瞬间,观察到激波产生,并向上游和下游传播,向上游传播的激波与空穴相互作用,导致水翼吸力面新生的附着型片状空穴回缩,直至完全溃灭.并且空穴溃灭激波存在回弹现象, 抑制了下一周期的空化发展.      

不同的分子间势对LJD状态方程的影响

本文分別采用6-exp,M-M和LJ(6,12)三种形式的分子间作用势,计算了LJD理论的压缩性因子,以检验由于势的不同对LJD状态方程造成的影响。     

Turbulence modelling and role of compressibility on oil spilling from a damaged double hull tank

The viscosity plays an important role, and a multiphase solver is necessary to numerically simulate the oil spilling from a damaged double hull tank (DHT). However, it is uncertain whether turbulence modelling is necessary, which turbulence model is suitable; and what the role of compressibility of the fluids is. This paper presents experimental and numerical investigations to address these issues for various cases representing different scenarios of the oil spilling, including grounding and collision. In the numerical investigations, various approaches to model the turbulence, including the large eddy simulation (LES), direct numerical simulation and the Reynolds average Navier–Stokes equation (RANS) with different turbulence models, are employed. Based on the investigations, it is suggested that the effective Reynolds numbers corresponding to both oil outflow and water inflow shall be considered when classifying the significance of the turbulence and selecting the appropriate turbulence models. This is confirmed by new lab tests considering the axial offset between the internal and the external holes on two hulls of the DHT. The investigations conclude for numerically simulating oil spilling from a damaged DHT that when the effective Re is smaller the RANS approaches should not be used and LES modelling should be employed; while when the effective Reynolds numbers is large, the RANS models may be used as they can give similar results to LES in terms of the height of the mixture in the ballast tank and discharge but costing much less CPU time. The investigation on the role of the compressibility of the fluid reveals that the compressibility of the fluid may be considerable in a small temporal‐spatial scale but plays an insignificant role on macroscopic process of the oil spilling. Copyright © 2016 John Wiley & Sons, Ltd.     

An improved progressive preconditioning method for steady non‐cavitating and sheet‐cavitating flows

An improved progressive preconditioning method for analyzing steady inviscid and laminar flows around fully wetted and sheet‐cavitating hydrofoils is presented. The preconditioning matrix is adapted automatically from the pressure and/or velocity flow‐field by a power‐law relation. The cavitating calculations are based on a single fluid approach. In this approach, the liquid/vapour mixture is treated as a homogeneous fluid whose density is controlled by a barotropic state law. This physical model is integrated with a numerical resolution derived from the cell‐centered Jameson's finite volume algorithm. The stabilization is achieved via the second‐and fourth‐order artificial dissipation scheme. Explicit four‐step Runge–Kutta time integration is applied to achieve the steady‐state condition. Results presented in the paper focus on the pressure distribution on hydrofoils wall, velocity profiles, lift and drag forces, length of sheet cavitation, and effect of the power‐law preconditioning method on convergence speed. The results show satisfactory agreement with numerical and experimental works of others. The scheme has a progressive effect on the convergence speed. The results indicate that using the power‐law preconditioner improves the convergence rate, significantly. Copyright © 2010 John Wiley & Sons, Ltd.     

Free (open) boundary condition: some experiences with viscous flow simulations

The free (or open) boundary condition (FBC, OBC) was proposed by Papanastasiou et al. (A new outflow boundary condition, International Journal for Numerical Methods in Fluids, 1992; 14:587–608) to handle truncated domains with synthetic boundaries where the outflow conditions are unknown. In the present work, implementation of the FBC has been tested in several benchmark problems of viscous flow in fluid mechanics. The FEM is used to provide numerical results for both cases of planar and axisymmetric domains under laminar, isothermal or non‐isothermal, steady‐state conditions, for Newtonian fluids. The effects of inertia, gravity, compressibility, pressure dependence of the viscosity, slip at the wall, and surface tension are all considered individually in the extrudate‐swell benchmark problem for a wide range of the relevant parameters. The present results extend previous ones regarding the applicability of the FBC and show cases where the FBC is inappropriate, namely in the extrudate‐swell problem with gravity or surface‐tension effects. Particular emphasis has been given to the pressure at the outflow, which is the most sensitive quantity of the computations. In all cases where FBC is appropriate, excellent agreement has been found in comparisons with results from very long domains. The formulation for Picard‐type iterations is given in some detail, and the differences with the Newton–Raphson formulation are highlighted regarding some computational aspects. Copyright © 2011 John Wiley & Sons, Ltd.     

The Hard Rod Fluid in a Uniform Gravitational Field

Abstract

The system of hard rods in an external, uniform gravitational field is studied and exact expressions obtained for the partition functions, and the one- and two-particle distributions. In principle all higher order distributions can be exactly expressed as finite sums and the Laplace transforms of the one-particle density and pressure for a semi-infinite system are reducible to single integrals. Limiting cases of weak field and small rod diameters are examined. In the former case, these results agree with Percus and Lebowitz's local density expansion. In the latter case, corrections to the barometric pressure and density laws are obtained. Finally some mathematical difficulties involved in the calculation of the virial expansion and distribution of roots of the grand partition function are mentioned.     

Ultrasonic study on binary liquid mixtures between some bromoalkanes and hydrocarbons

Densities and ultrasound velocities for the binary mixtures of 1-bromobutane+benzene and 1,4-dimethylbenzene and of 1-bromopentane+cyclohexane and benzene have been measured at 308.15 K. Adiabatic compressibilities (beta(ad)), and Wada's constants (W) have also been evaluated as a function of composition. The ultrasound velocities decrease, attains a minimum and then increase with increase in mole fractions of hydrocarbons in the binary mixtures except in the case of 1-bromopentane+benzene binary mixtures where the variation is just the reverse. Dependence of adiabatic compressibilities with mole fractions of hydrocarbons is sigmoid. The non-ideal behaviour of the systems studied is explained on the basis of dipole-induced dipole interactions.