The effects of flow split ratio and flow rate in manifolds

This paper reports a combined experimental and numerical investigation of three-dimensional steady turbulent flows in inlet manifolds of square cross-section. Predictions and measurements of the flows were carried out using computational fluid dynamics and laser Doppler anemometry techniques respectively. The flow structure was characterized in detail and the effects of flow split ratio and inlet flow rate were studied. These were found to cause significant variations in the size and shape of recirculation regions in the branches, and in the turbulence levels. It was then found that there is a significant difference between the flow rates through different branches. The performance of the code was assessed through a comparison between predictions and measurements. The comparison demonstrates that all important features of the flow are well represented by the predictions.     

Decay of confined,two-dimensional,spatially periodic arrays of vortices: A numerical investigation

The disarrangement of a perturbed lattice of vortices was studied numerically. The basic state is an exponentially decaying, exact solution of the Navier-Stokes equations. Square arrays of vortices with even numbers of vortex cells along each side were perturbed and their evolution was investigated. Whether the energy in the perturbation grows somewhat before it decays or decays monotonically depends on the initial strength of the vortices of the basic state, the extent of lateral confinement and the structure of the perturbation. The critical condition for temporally local instability, i.e. the critical amplitude of the basic state that must be exceeded to allow energy transfer from the basic state to the perturbation, is discussed. In the strongly confined case of a square lattice of four vortices the appearance of enchancement of global rotation is the result of energy transfer from the basic state to a temporally local unstable mode. Energy is transferred from the basic state to larger-scaled structures (inverse cascade) only if the scales of the larger structures are inherently contained in the initial structure of the perturbation. The initial structure of the double array of vortices is not maintained except for a very special form of perturbation. The facts that large scales decay more slowly than small scales and that, when non-linearities are sufficiently strong, energy is transferred from one scale to another explain the differences in the disarrangement process for different initial strengths of the vortices of the basic state. The stronger vortices, i.e. the vortices perturbed in a manner that increases their strength, tend to dominate the weaker vortices. The pairing and subsequent merging (or capture) of vortices of like sense into larger-scale vortices are described in terms of peaks in the evolution of the square root of the palinstrophy divided by the enstrophy.     

潜热型热功能流体强化换热分析

本文建立了分析带有相变微胶囊的潜热型热功能流体的流动和换热过程数理模型,应用有限差分法和移动热源法进行联合求解。计算结果表明,相变微胶囊的加入:较好地提高了流体的换热性能。获得了相变颗粒大小和体积分数对强化换热影响等结果。计算结果对该类流体的设计和应用提供了理论依据。     

The modified characteristic finite difference fractional steps method for the coupled system of fluid dynamics in porous media and its analysis

For the coupled system of multilayer fluid dynamics in porous media, the modified characteristic finite difference fractional steps method applicable to parallel arithmetic is put forward and two‐dimensional and three‐dimensional schemes are used to form a complete set. Some techniques, such as calculus of variations, energy method, piecewise biquadratic interpolation, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates in L² norm are derived to determine the error in the approximate solution. This method has already been applied to the numerical simulation of multilayer fluid dynamics in porous media. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 665–681, 2003.     

A new VOF‐based numerical scheme for the simulation of fluid flow with free surface. Part II: application to the cavity filling and sloshing problems

Finite element analysis of fluid flow with moving free surface has been performed in 2‐D and 3‐D. The new VOF‐based numerical algorithm that has been proposed by the present authors (Int. J. Numer. Meth. Fluids, submitted) was applied to several 2‐D and 3‐D free surface flow problems. The proposed free surface tracking scheme is based on two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby‐cell to determine the fluid volume flux at each cell boundary. The proposed numerical algorithm has been applied to 2‐D and 3‐D cavity filling and sloshing problems in order to demonstrate the versatility and effectiveness of the scheme. The proposed numerical algorithm resolved successfully the free surfaces interacting with each other. The simulated results demonstrated applicability of the proposed numerical algorithm to the practical problems of large free surface motion. It has been also demonstrated that the proposed free surface tracking scheme can be easily implemented in any irregular non‐uniform grid systems and can be extended to 3‐D free surface flow problems without additional efforts. Copyright © 2003 John Wiley & Sons, Ltd.     

The dynamics of a two-component fluid in the presence of capillary forces

In the present paper we study the qualitative behavior ast→∞ of the solution of the Cauchy problem for a system of equations describing a dynamics of a two-component viscous fluid. The model under consideration takes into account the mutual diffusion of the fluid components as well as their capillary interaction. We describe the ω-limit set of trajectories of the dynamical system generated by the problem. It is proved that the stationary solution of the problem, is a homogeneous stationary distribution of one of the components, is asymptotically stable. Any other stationary solution is not asymptotically stable and is even unstable if there are no close stationary solutions corresponding to a smaller energy level. Translated fromMatematicheskie Zametki, Vol. 62, No. 2, pp. 293–305, August, 1997. Translated by A. M. Chebotarev     

Bianchi type-V model with a perfect fluid and A-term

A self-consistent system of gravitational field with a binary mixture of perfect fluid and dark energy given by a cosmological constant has been considered in Bianchi Type-V universe. The perfect fluid is chosen to be obeying either the equation of state p=γρ with γ ε |0,1| or a van der Waals equation of state. The role of A-term in the evolution of the Bianchi Type-V universe has been studied.     

存在滑移时两圆球间的幂律流体挤压流动

基于Reynolds润滑理论分析了壁面滑移对任意圆球颗粒间幂律流体的挤压流动的影响。研究表明有壁面滑移时挤压流动的粘性力可通过引进本文定义的滑移修正系数分离出无滑移解。推导出的挤压力滑移修正系数是一积分表达式，依赖于滑移参数、幂律指数、球间隙和积分上限。一般地壁面滑移导致粘性力减小，粘性力的减小量随幂律指数的增大而增大，表明壁面滑移对剪切增稠流变材料有更大的影响；粘性力的减小量还随着滑移参数的增大而增加，而这恰与假设相符合；粘性力的减小量又随球间隙减小或积分上限的增大（从液桥情况到完全浸渍）而减小直到趋于常数，这一特性在离散元模拟时可以有效地减少计算量。     

One case of vortex motion in fluid

This text presents an English translation of the significant paper [6] on vortex dynamics published by the outstanding Russian scientist S. A. Chaplygin, which seems to have almost escaped the attention of later investigators in this field. Although it was published more than a century ago, in our opinion it is still interesting and valuable. __________ Originally published in: Trudy otdeleniya fizicheskikh nauk imperatorskogo Moskovskogo obshchestva lyubitelei estestvoznaniya, antropologii i etnografii (Transactions of the Physical Section of Moscow Society of Friends of Natural Sciences, Anthropology and Ethnography), 1903, Vol. 11, No. 2, pp. 11–14. Translated from Russian by S. Ramodanov; edited by D. Blackmore; commented by V. V. Meleshko (Department of Theoretical and Applied Mechanics, Faculty of Machanics and Mathematics, Kiev National Taras Shevchenko University, 01030 Kiev, Ukraine. E-mail: meleshko@univ.kiev.ua) and G.J.F. van Heijst (Fluid Dynamics Laboratory, Faculty of Applied Physics, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands. E-mail: g.j.f.v.heijst@tue.nl).     

Measurement and prediction of solubilities and diffusion coefficients of carbon dioxide in starch–water mixtures at elevated pressures

The solubility and diffusion coefficient of carbon dioxide in intermediate‐moisture starch–water mixtures were determined both experimentally and theoretically at elevated pressures up to 16 MPa at 50 °C. A high‐pressure decay sorption system was assembled to measure the equilibrium CO₂ mass uptake by the starch–water system. The experimentally measured solubilities accounted for the estimated swollen volume by Sanchez–Lacombe equation of state (S‐L EOS) were found to increase almost linearly with pressure, yielding 4.0 g CO₂/g starch–water system at 16 MPa. Moreover, CO₂ solubilities above 5 MPa displayed a solubility increase, which was not contributed by the water fraction in the starch–water mixture. The solubilities, however, showed no dependence on the degree of gelatinization (DG) of starch. The diffusion coefficient of CO₂ was found to increase with concentration of dissolved CO₂, which is pressure‐dependent, and decrease with increasing DG in the range of 50–100%. A free‐volume‐based diffusion model proposed by Areerat was employed to predict the CO₂ diffusivity in terms of pressure, temperature, and the concentration of dissolved CO₂. S‐L EOS was once more used to determine the specific free volume of the mixture system. The predicted diffusion coefficients showed to correlate well with the measured values for all starch–water mixtures. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 607–621, 2006