An improved numerical simulator for multiphase flow in porous media

In basin modelling the thermodynamics of a multicomponent multiphase fluid flux are computationally too expensive when derived from an equation of state and the Gibbs equality constraints. In this article we present a novel implicit molar mass formulation technique using binary mixture thermodynamics. The two proposed solution methods, with and without cross derivative terms between components, are based on a preconditioned Newton‐GMRES scheme for each time‐step with analytical computation of the derivatives. These new algorithms reduce significantly the numerical effort for the computation of the molar masses, and we illustrate the behavior of these methods with numerical computations. Copyright © 2004 John Wiley & Sons Ltd.     

Assessment of two‐equation turbulence modelling for high Reynolds number hydrofoil flows

This paper presents an evaluation of the capability of turbulence models available in the commercial CFD code FLUENT 6.0 for their application to hydrofoil turbulent boundary layer separation flow at high Reynolds numbers. Four widely applied two‐equation RANS turbulence models were assessed through comparison with experimental data at Reynolds numbers of 8.284×106 and 1.657×107. They were the standard k–εmodel, the realizable k–εmodel, the standard k–ωmodel and the shear‐stress‐transport (SST) k–ωmodel. It has found that the realizable k–εturbulence model used with enhanced wall functions and near‐wall modelling techniques, consistently provides superior performance in predicting the flow characteristics around the hydrofoil. Copyright © 2004 John Wiley & Sons, Ltd.     

High‐order accurate numerical solutions of incompressible flows with the artificial compressibility method

A high‐order accurate, finite‐difference method for the numerical solution of incompressible flows is presented. This method is based on the artificial compressibility formulation of the incompressible Navier–Stokes equations. Fourth‐ or sixth‐order accurate discretizations of the metric terms and the convective fluxes are obtained using compact, centred schemes. The viscous terms are also discretized using fourth‐order accurate, centred finite differences. Implicit time marching is performed for both steady‐state and time‐accurate numerical solutions. High‐order, spectral‐type, low‐pass, compact filters are used to regularize the numerical solution and remove spurious modes arising from unresolved scales, non‐linearities, and inaccuracies in the application of boundary conditions. The accuracy and efficiency of the proposed method is demonstrated for test problems. Copyright © 2004 John Wiley & Sons, Ltd.     

A subdomain boundary element method for high‐Reynolds laminar flow using stream function‐vorticity formulation

The paper presents a new formulation of the integral boundary element method (BEM) using subdomain technique. A continuous approximation of the function and the function derivative in the direction normal to the boundary element (further ‘normal flux’) is introduced for solving the general form of a parabolic diffusion‐convective equation. Double nodes for normal flux approximation are used. The gradient continuity is required at the interior subdomain corners where compatibility and equilibrium interface conditions are prescribed. The obtained system matrix with more equations than unknowns is solved using the fast iterative linear least squares based solver. The robustness and stability of the developed formulation is shown on the cases of a backward‐facing step flow and a square‐driven cavity flow up to the Reynolds number value 50 000. Copyright © 2004 John Wiley & Sons, Ltd.     

Experimental investigation and theoretical modeling of liquid entrainment in a horizontal tee with a vertical-up branch

This article describes a comprehensive literature review of liquid entrainment in horizontal pipes with vertical-up branches. Deficiencies in the available data and correlations were identified. The Air–water Test Loop for Advanced Thermal–hydraulic Studies (ATLATS) was constructed and entrainment onset and steady-state data were collected for a wide range of flow conditions. Using new insights gained from experimental testing, the authors developed a model for predicting the onset of entrainment and steady-state entrainment rate. Previously published correlations, along with the new model, are compared against all available data. The new model shows very good agreement with the onset data, but is not very good at predicting branch quality at high liquid flow rates.     

Modeling Media with Oriented Structures

Heterogeneity is typically the result of space variability of soil parameters at different scales. Soil anisotropy may be defined as the spatial persistence in some direction only, across coarse-grid elements, of heterogeneous structures with different characteristic lengths in different directions. One can account for the effect of these structures by upscaling soil properties. Analyzing flow in a strongly anisotropic structured soil at different scales evidences how transverse dispersion reduces to a subscale process, leading to mixing within the conductive structures.     

On the hysteresis of vortex shedding in a periodically-varying flow

Efforts are made to explore the hysteresis characteristics of vortex shedding in a pipe flow, whose velocity varies periodically in time. Results obtained show that during acceleration of the flow, the vortex strength tends to be stronger, whereas during deceleration of the flow, the situation is reversed. As reconstructed from the velocity signals measured at a point in the flow field, the shed vortex arrays appear to possess uneven vortex strengths in response to periodically-varying incoming flows. Furthermore, in the hysteresis range, the streamwise spacings between the vortices appear to be unequal.     

Control of self‐organized low‐dimensional morphology in poly(styrene‐b‐4vinylpyridine)/polystyrene blend thin films

The surface morphologies of poly(styrene‐b‐4vinylpyridine) (PS‐b‐P4VP) diblock copolymer and homopolystyrene (hPS) binary blend thin films were investigated by atomic force microscopy as a function of total volume fraction of PS (?_PS) in the mixture. It was found that when hPS was added into symmetric PS‐b‐P4VP diblock copolymers, the surface morphology of this diblock copolymer was changed to a certain degree. With ?_PS increasing at first, hPS was solubilized into the corresponding domains of block copolymer and formed cylinders. Moreover, the more solubilized the hPS, the more cylinders exist. However, when the limit was reached, excessive hPS tended to separate from the domains independently instead of solubilizing into the corresponding domains any longer, that is, a macrophase separation occurred. A model describing transitions of these morphologies with an increase in ?_PS is proposed. The effect of composition on the phase morphology of blend films when graphite is used as a substrate is also investigated. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 3496–3504, 2004     

The exponential power law: partial wetting kinetics and dynamic contact angles

An equation for the kinetics of partial drop spreading is proposed. This equation was empirically derived from experimental data for the spreading kinetics of partially wetting liquids in terms of the wet area versus time. The equation has the form of an exponential power law (EPL), and transforms into the well-known power law for complete wetting, when the equilibrium contact angle approaches zero. The EPL fits very well available experimental data. To lend additional support to the validity of this generalized equation, it will be demonstrated that when it is transformed to present the dynamic contact angle (DCA), it fits very well DCA experimental data for other wetting processes, such as capillary flow and tape coating.     

密封容器中后处理制备Tl2Ba2CaCu2Oy高温超导薄膜及其性能研究

以铝酸镧（001）单晶为基片，采用两步法制备Tl2Ba2CaCu2Oy（Tl 2212）高 温超导薄膜.首先，利用脉冲激光沉积（PLD）工艺沉积Ba2CaCu2Ox非晶前驱体薄膜；然后，前驱体薄膜在高温（720—850℃）下密封钢容器里铊化结晶形成Tl 2212薄膜.XRD结果表明Tl2212 薄膜是沿c轴方向生长的，其相组成为Tl 2212，摇摆曲线（0012）的半高宽为0.72° ，SEM图像显示其表面光滑平整，其零电阻温度为106.2K. 关键词： Tl 2212超导薄膜 脉冲激光沉积