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 high‐resolution scheme for low Mach number flows

A method for computing low Mach number flows using high‐resolution interpolation and difference formulas, within the framework of the Marker and Cell (MAC) scheme, is presented. This increases the range of wavenumbers that are properly resolved on a given grid so that a sufficiently accurate solution can be obtained without extensive grid refinement. Results using this scheme are presented for three problems. The first is the two‐dimensional Taylor–Green flow which has a closed form solution. The second is the evolution of perturbations to constant‐density, plane channel flow for which linear stability solutions are known. The third is the oscillatory instability of a variable density plane jet. In this case, unless the sharp density gradients are resolved, the calculations would breakdown. Under‐resolved calculations gave solutions containing vortices which grew in place rather than being convected out. With the present scheme, regular oscillations of this instability were obtained and vortices were convected out regularly. Stable computations were possible over a wider range of sensitive parameters such as density ratio and co‐flow velocity ratio. 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.     

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.     

Desorption behavior of a surfactant in a linear low‐density polyethylene blend at elevated temperatures

The desorption behavior of a surfactant in a linear low‐density polyethylene (LLDPE) blend at elevated temperatures of 50, 70, and 80 °C was studied with Fourier transform infrared spectroscopy. The composition of the LLDPE blend was 70:30 LLDPE/low‐density polyethylene. Three different specimens (II, III, and IV) were prepared with various compositions of a small molecular penetrant, sorbitan palmitate (SPAN‐40), and a migration controller, poly(ethylene acrylic acid) (EAA), in the LLDPE blend. The calculated diffusion coefficient (D) of SPAN‐40 in specimens II, III, and IV, between 50 and 80 °C, varied from 1.74 × 10^?11 to 6.79 × 10^?11 cm²/s, from 1.10 × 10^?11 to 5.75 × 10^?11 cm²/s, and from 0.58 × 10^?11 to 4.75 × 10^?11 cm²/s, respectively. In addition, the calculated activation energies (E_D) of specimens II, III, and IV, from the plotting of ln D versus 1/T between 50 and 80 °C, were 42.9, 52.7, and 65.6 kJ/mol, respectively. These values were different from those obtained between 25 and 50 °C and were believed to have been influenced by the interference of Tinuvin (a UV stabilizer) at elevated temperatures higher than 50 °C. Although the desorption rate of SPAN‐40 increased with the temperature and decreased with the EAA content, the observed spectral behavior did not depend on the temperature and time. For all specimens stored over 50 °C, the peak at 1739 cm^?1 decreased in a few days and subsequently increased with a peak shift toward 1730 cm^?1. This arose from the carbonyl stretching vibration of Tinuvin, possibly because of oxidation or degradation at elevated temperatures. In addition, the incorporation of EAA into the LLDPE blend suppressed the desorption rate of SPAN‐40 and retarded the appearance of the 1730 cm^?1 peak. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 1114–1126, 2004     

Effects of poly[styrene‐b‐(ethylene‐co‐butylene)‐b‐styrene] on the charge distributions of low‐density polyethylene/polystyrene blends

The effect of the triblock copolymer poly[styrene‐b‐(ethylene‐co‐butylene)‐b‐styrene] (SEBS) on the formation of the space charge of immiscible low‐density polyethylene (LDPE)/polystyrene (PS) blends was investigated. Blends of 70/30 (wt %) LDPE/PS were prepared through melt blending in an internal mixer at a blend temperature of 220 °C. The amount of charge that accumulated in the 70% LDPE/30% PS blends decreased when the SEBS content increased up to 10 wt %. For compatibilized and uncompatibilized blends, no significant change in the degree of crystallinity of LDPE in the blends was observed, and so the effect of crystallization on the space charge distribution could be excluded. Morphological observations showed that the addition of SEBS resulted in a domain size reduction of the dispersed PS phase and better interfacial adhesion between the LDPE and PS phases. The location of SEBS at a domain interface enabled charges to migrate from one phase to the other via the domain interface and, therefore, resulted in a significant decrease in the amount of space charge for the LDPE/PS blends with SEBS. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 2813–2820, 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.     

交通灯控制下城市主干道双车道多速元胞自动机交通流模型研究

提出一个城市主干道双车道多速元胞自动机交通流模型，此模型采用开放性 边界条件，在考虑诸多实际因素影响情况下，研究了主干道中车站的设置、交通灯绿信比对 车流量和车流速度的影响.计算机模拟所得到的基本图（流量 速度图）能较好地反映在交 通灯控制下城市主干道交通流真实状况. 关键词： 元胞自动机模型 交通流 主干道