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.     

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.     

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.     

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

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

多孔介质BISQ模型中的慢纵波

着重研究了多孔介质BISQ模型中慢纵波的基本特性.给出了BISQ模型下慢纵波速度 和衰减的低频近似公式.与Biot理论对比，BISQ模型中慢纵波的衰减随频率降低急剧增大， 且随喷射流长度的减小而增加；相速度随喷射流长度的减小而增加，其低频极限值不是零； 孔隙流体位移与固相骨架位移之比的幅值随喷射流长度的增加而减小，其相位特点与Biot模 型预测的不同；在流体与孔隙介质的边界上可产生更大的渗流.为对比，同时也给出快纵波 的行为.依据BISQ模型可推断：非黏滞流体饱和孔隙介质中不存在喷射流机理；BISQ模型中 关键词： 多孔介质 喷射流 慢纵波 动力协调     

Hydromagnetic rotating flow of a fourth‐order fluid past a porous plate

The rotating flow in the presence of a magnetic field is a problem belonging to hydromagnetics and deserves to be more widely studied than it has been to date. In the non‐linear regime the literature is scarce. We develop the governing equations for the unsteady hydromagnetic rotating flow of a fourth‐order fluid past a porous plate. The steady flow is governed by a boundary value problem in which the order of differential equations is more than the number of available boundary conditions. It is shown that by augmenting the boundary conditions based on asymptotic structures at infinity it is possible to obtain numerical solutions of the nonlinear hydromagnetic equations. Effects of uniform suction or blowing past the porous plate, exerted magnetic field and rotation on the flow phenomena, especially on the boundary layer structure near the plate, are numerically analysed and discussed. The flow behaviours of the Newtonian fluid and second‐, third‐ and fourth‐order non‐Newtonian fluids are compared for the special flow problem, respectively. Copyright © 2004 John Wiley & Sons, Ltd.     

The long‐time behaviour of the thermoconvective flow in a porous medium

For the Boussinesq approximation of the equations of coupled heat and fluid flow in a porous medium we show that the corresponding system of partial differential equations possesses a global attractor. We give lower and upper bounds of the Hausdorff dimension of the attractor depending on a physical parameter of the system, namely the Rayleigh number of the flow. Numerical experiments confirm the theoretical findings and raise new questions on the structure of the solutions of the system. Copyright © 2004 John Wiley & Sons, Ltd.     

Some non‐standard problems in viscous flow

Energy bounds are derived for Dirichlet type boundary value problems for the Navier–Stokes and Stokes equations when a combination of the solution values initially and at a later time is prescribed. The bounds are obtained by means of a differential inequality and imply uniqueness and continuous data dependence of the solutions for a range of values of the parameter in the non‐standard auxiliary condition. Copyright © 2004 John Wiley & Sons, Ltd.