不同交通流状态下桥梁气动导数的风洞试验研究

为了考虑实际运营车辆对桥梁气动导数的影响,根据车辆密度模拟了三种交通流状态,基于强迫振动装置,分别对每个交通流和无车状态下的桥梁气动导数进行风洞试验研究,讨论了不同攻角下不同车流的车辆对桥梁气动导数的影响,探究了车辆对气动导数影响的百分比以及气动导数变化量的变化规律。研究结果表明:不同攻角下不同车流的车辆均对直接导数A*2、H*4和交叉导数A*4、H*2影响显著,A*2、A*3变化量随着折减风速有一定的变化规律。虽然不同攻角下不同车流的车辆对气动导数的影响程度及影响规律不同,并且车流的繁忙程度对大多数气动导数的影响规律不明显,但是车辆对桥梁气动导数的影响不容忽略。     

Existence of supersonic zones in three-dimensional separation flows

Up till now the region of three-dimensional separation flows which occur with supersonic flow past obstacles has received insufficient study. Supersonic flow with a Mach number of 2.5 past a cylinder mounted on a plate was studied in [1]. A local zone with supersonic velocities was found in the reverse subsonic flow region ahead of the cylinder. Its presence is explained by the three-dimensional nature of the flow. Similar supersonic zones are not observed in the case of supersonic flow over plane and axisymmetric steps.The present paper presents the results of experimental studies whose objective was refinement of the flow pattern ahead of a cylinder on a plate and the study of the local supersonic zones.The experiments were performed in a supersonic wind tunnel with a freestream Mach number M₁=3.11. The 24-mm-diameter cylinder with pressure taps along the generating line was mounted perpendicular to the surface of a sharpened plate. The distance from the plate leading edge to the cylinder axis wasl ₀=140 mm. The plate was pressure tapped along the flow symmetry axis. The Reynolds number was Rl ₀=u₀ l ₀/v ₁, Rl ₀=1.87.10⁷, where u₁ andv ₁ are the freestream velocity and the kinematic viscosity, respectively. The pressures were measured using a Pilot probe with internal and external diameters of 0.15 and 0.9 mm, respectively.The probe was displaced in the flow symmetry plane at a distance of 1.6 mm from the plate surface and at a distance of 1.1 mm along the leading generator of the cylinder. The flow on the surface of the plate and cylinder was studied with the aid of a visualization composition and the flow past the model was photographed with a schlieren instrument. Typical patterns of the visualization composition distribution and the pressure distribution curves over the plate surface, and also photographs of the flow past the model, are shown in [1].     

Tailoring turbulence with an active grid

Using an active grid in a wind tunnel, we generate homogeneous shear turbulence and initiate turbulent boundary layers with adjustable properties. Homogeneous shear turbulence is characterized by a constant gradient of the mean velocity and a constant turbulence intensity. It is the simplest anisotropic turbulent flow thinkable, and it is generated traditionally by equipping a wind tunnel with screens which have a varying transparency and flow straighteners. This is not done easily, and the reachable turbulence levels are modest. We describe a new technique for generating homogeneous shear turbulence using an active grid only. Our active grid consists of a grid of rods with attached vanes which can be rotated by servo motors. We control the grid by prescribing the time-dependent angle of each axis. We tune the vertical transparency profile of the grid by setting appropriate angles of each rod such as to generate a uniform velocity gradient, and set the rods in flapping motion around these angles to tailor the turbulence intensity. The Taylor Reynolds number reached was R _λ = 870, the shear rate S = ∂U/∂y = 9.2 s⁻¹, the nondimensional shear parameter S ^*≡ Sq ²/ε = 12 and u = 1.4 ms⁻¹. As a further application of this idea we demonstrate the generation of a simulated atmospheric boundary layer in a wind tunnel which has tunable properties. This method offers a great advantage over the traditional one, in which vortex-generating structures need to be placed in the wind tunnel to initiate a fat boundary layer.     

Air flow aerodynamic on a wall-mounted hemisphere for various turbulent boundary layers

Air flow field around a surface-mounted hemisphere of a fixed height for two different turbulent boundary layers (thin and thick) are investigated experimentally and numerically. Flow measurements are performed in a wind tunnel using hot-wire anemometer and streamwise component of velocity fluctuation are calculated using a special developed program of the hardware system. Mean surface pressure coefficients and velocity field for the same hemisphere are determined by the numerical simulation. Turbulent flow field and intensity are measured for two types of boundary layers and compared at various sections in both streamwise and spanwise directions. Numerical scheme based on finite volume and SIMPLE algorithm is used to treat pressure and velocity coupling. Studies are performed for Reynolds number, Re_H = 32,000. Based on the numerical simulation using RNG k–ε turbulence model, flow pathlines, separation region and recirculation area are determined for the two types of turbulent boundary layer flows and complex flow field and recirculation regions are identified and presented graphically.     

Study on the instability in separating–reattaching flow over a surface-mounted rib

A numerical study on the flow structure and instability in the separated–reattached flow over a surface-mounted rib at Re = 1000 is performed using large eddy simulation. It is found that the phenomenon of vortex pairing, which has been extensively observed in similar flows, exists in the separation zone. Based on the spectral analysis, the Kelvin–Helmholtz (K-H) instability of shear layer at St ≈ 0.361 (St ≡ fh/U₀) and its subharmonic at St ≈ 0.18 are found. It is assumed that the K-H instability reduces to its subharmonic through the vortex pairing. This process is further confirmed by the flow visualisation. The two-dimensional (2D) structures are subjected to sinusoidal undulation along the spanwise and observed to undergo helical pairing process, which is attributed to the transformation of 2D structures into 3D. However, the low frequency due to flapping of the shear layer is not found.     

An improved ghost‐cell immersed boundary method for compressible flow simulations

This study presents an improved ghost‐cell immersed boundary approach to represent a solid body in compressible flow simulations. In contrast to the commonly used approaches, in the present work, ghost cells are mirrored through the boundary described using a level‐set method to farther image points, incorporating a higher‐order extra/interpolation scheme for the ghost‐cell values. A sensor is introduced to deal with image points near the discontinuities in the flow field. Adaptive mesh refinement is used to improve the representation of the geometry efficiently in the Cartesian grid system. The improved ghost‐cell method is validated against four test cases: (a) double Mach reflections on a ramp, (b) smooth Prandtl–Meyer expansion flows, (c) supersonic flows in a wind tunnel with a forward‐facing step, and (d) supersonic flows over a circular cylinder. It is demonstrated that the improved ghost‐cell method can reach the accuracy of second order in L¹ norm and higher than first order in L^∞ norm. Direct comparisons against the cut‐cell method demonstrate that the improved ghost‐cell method is almost equally accurate with better efficiency for boundary representation in high‐fidelity compressible flow simulations. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical solutions of the shallow water equations with discontinuous bed topography

A simple scheme is developed for treatment of vertical bed topography in shallow water flows. The effect of the vertical step on flows is modelled with the shallow water equations including local energy loss terms. The bed elevation is denoted with z_b^‐ for the left and z_b⁺ for the right values at each grid point, hence exactly representing a discontinuity in the bed topography. The surface gradient method (SGM) is generalized to reconstruct water depths at cell interfaces involving a vertical step so that the fluxes at the cell interfaces can accurately be calculated with a Riemann solver. The scheme is verified by predicting a surge crossing a step, a tidal flow over a step and dam‐break flows on wet/dry beds. The results have shown good agreements compared with analytical solutions and available experimental data. The scheme is efficient, robust, and may be used for practical flow calculations. Copyright © 2002 John Wiley & Sons, Ltd.     

Adaptive wall functions for the v2‐f turbulence model

Adaptive wall functions for the v²‐f turbulence model have been derived for the flow over a flat plate at zero pressure gradient. These wall functions were implemented via tables for the turbulence quantities and the friction velocity u_τ. A special treatment for the ε and f boundary conditions is proposed. On fine grids (y⁺<1) this approach yields results consistent with the wall integration solution. Detailed numerical results are presented for a zero pressure gradient boundary layer and separated flow over a ramp. Copyright © 2006 John Wiley & Sons, Ltd.     

Simulation of Wind Flow Around a Building with a k–ε Model

The three-dimensional numerical simulation of airflow around a building using a k–ε two-equation turbulence model is presented in this paper. Several cases of numerical simulation of airflow around a building are carried out to estimate the influence of mesh spacing on simulated results. The accuracy of simulations is examined by comparing the predicted results with wind-tunnel experiments. It is confirmed that numerical simulations by means of the k–ε model reproduce the velocity fields well when using fine mesh resolution. In the latter part of the paper, the simulation method is applied to predict the flow field around a building with different width-to-height ratios, under light wind conditions. Received 16 June 1999 and accepted 20 July 2000     

The wake flow control behind a circular cylinder using ion wind

Active and passive flow control methods have been studied for decades, but there have been only a few studies of flow control methods using ion wind, which is the bulk motion of neutral molecules driven by locally ionized air from a corona discharge. This paper describes an experimental study of ion wind wake control behind a circular cylinder. The experimental conditions consisted of a range of electrohydrodynamic numbers—the ratio of an electrical body force to a fluid inertial force—from 0 to 2 and a range of Reynolds numbers from 4×10³ to 8×10³. Pressure distributions over the cylinder surface were measured and flow visualizations were carried out using a smoke-wire method. The flow visualizations confirmed that ion wind significantly affects the wake structure behind a circular cylinder, and that the pressure drag can be dramatically reduced by superimposing ion wind.List of symbols BR blockage ratio - C _d coefficient of the pressure drag - C _p coefficient of the surface pressure, 2(p–p ₀)/(U ₀ ²) - C _pb coefficient of the base surface pressure, 2(p _b–p ₀)/(U ₀ ²) - D diameter of the cylinder - D _P pressure drag - d _p diameter of particle - E the electric field - F _e Coulombian force (qE) - F _v viscous force - H wire-to-cylinder spacing - I total electric current (A) - L the axial length of cylinder (m) - N _EHD electrohydrodynamic number - p _b base pressure of cylinder at =180° - p ₀ reference static pressure at 10D upstream - q the charge on the particle - R radius of the cylinder - V applied voltage (kV) - U ₀ mean flow velocity (m/s) - ion mobility in air (m²/(s V)) - ₀ permittivity of free space - viscosity of fluid (kg/ms) - density of fluid (kg/m³) - installation angle of a wire electrode (°)     

不同展长比旗帜在均匀来流中摆动的实验研究

为研究展长比对旗帜摆动特性的影响,在低速风洞内开展了质量比M*为0.6,展长比H*为0.2～1.5的柔性旗帜在不同来流速度下的失稳实验.利用高速摄影机和图像处理技术分析了无量纲速度U*和展长比H*对旗帜运动学特征(频率和振幅)的影响;利用自行研制的天平测量了旗帜受到的阻力,首次给出了旗帜的动力学特征.结果发现,保持展长...     

Visualization of Supersonic Low-Noise Flow Over Surface-Mounted Two-Dimensional Prisms

An experimental study of supersonic flow over two-dimensional surface-mounted prisms is carried out in a Mach 3 low-noise wind tunnel. The noise level of this supersonic wind tunnel, defined as the root mean-square Pitot pressure fluctuation normalized by the mean Pitot pressure, can be reduced to about 0.37%. The nanotracer planar laser scattering (NPLS) technique is used to analyze the influence of the prism geometry and the oncoming flow conditions on the typical flow structures including separation and reattachment shocks. With increase in the prism height the induced shocks move upstream. At a constant streamwise length L of a prism the timeaveraged NPLS images show that the length of the downstream recirculation region increases from 0.8L to 1.2L, when the prism height H changes from 3 to 5 mm. As compared with the flow structures occurring downstream of the prisms, the upstream flow structures are more susceptible to the oncoming boundary layer and are considerably different in laminar and turbulent flows. The separation shock wave is clearly visible in turbulent flow even for the 1-mm prism, whereas in the case of laminar flow there is no a distinct shock wave upstream of this prism. At the same time, the location of the flow reattachment and the angle of the reattachment shock wave in the downstream flow remain almost the same in both two flow regimes.     

Boundary layer structure of oscillatory open-channel shallow flows over smooth and rough beds

The boundary layer structure of oscillatory shallow open channel flows has been studied in a wide flume. Fluorescence solution was released at a porous rough bed through a diffuser covered by gravel of 0.5 cm grain size. A planar laser-induced fluorescence (PLIF) system was used to visualise the dye plumes in both vertical and horizontal planes for a qualitative understanding of the roles of large-scale flow structures in mass transport. A variety of tests were conducted for a range of oscillatory periods (30–240 s), water depths (3–16 cm) and velocity amplitudes (0.027–0.325 m/s), which cover a wide range of oscillatory flows with Reynolds numbers Re _a varied from 0.3 × 10⁴ (laminar) to 2.1 × 10⁶ (fully turbulent). For quantitative investigation, a novel technique, namely combined laser-induced fluorescence (LIF) and 2D laser Doppler velocimetry (LDV) (LIF/LDV), was developed and used to measure the velocity and solute concentration simultaneously in a vertical plane over 50 cycles. From the dye plumes revealed by the PLIF in transitional flows, there are different patterns of flow structure and solute transport with three representative stages of acceleration, deceleration and flow reversal. In the acceleration stage, turbulence was suppressed with dye layers adhering to the surface with little vertical mass transport. In the deceleration stage, flame-like turbulent structures occurred when turbulence generation was prominent. This was investigated quantitatively by recording the percentage occurrence of the adhered smooth layers per cycle. For those smooth bed cases with Re _a < 1.8 × 10⁵, the adhered smooth dye layers type of boundary layer occupied 100% of the oscillation period. Over a sufficiently high Re _a , a rough bed can generate fully turbulent oscillatory flows without the appearance of adhering dye layers. Between these two extremes, a transitional flow regime occurs in a wide range of flow conditions: Re _a > 2.7 × 10⁴ over the rough bed and Re _a > 8.3 × 10⁶ over a smooth bed.     

Accuracy of out-of-plane vorticity measurements derived from in-plane velocity field data

 A study of the errors in out-of-plane vorticity (ω _z ) calculated using a local χ² fitting of the measured velocity field and analytic differentiation has been carried out. The primary factors of spatial velocity sampling separation and random velocity measurement error have been investigated. In principle the ω _z error can be decomposed into a bias error contribution and a random error contribution. Theoretical expressions for the transmission of the random velocity error into the random vorticity error have been derived. The velocity and vorticity field of the Oseen vortex has been used as a typical vortex structure in this study. Data of different quality, ranging from exact velocity vectors of analytically defined flow fields (Oseen vortex flow) sampled at discrete locations to computer generated digital image frames analysed using cross-correlation DPIV, have been investigated in this study. This data has been used to provide support for the theoretical random error results, to isolate the different sources of error and to determine their effect on ω _z measurements. A method for estimating in-situ the velocity random error is presented. This estimate coupled with the theoretically derived random error transmission results for the χ² vorticity calculation method can be used a priori to estimate the magnitude of the random error in ω _z . This random error is independent of a particular flow field. The velocity sampling separation is found to have a profound effect on the precise determination of ω _z by introducing a bias error. This bias error results in an underestimation of the peak vorticity. Simple equations, which are based on a local model of the Oseen vortex around the peak vorticity region, allowing the prediction of the ω _z bias error for the χ² vorticity calculation method, are presented. An important conclusion of this study is that the random error transmission factor and the bias error cannot be minimised simultaneously. Both depend on the velocity sampling separation, but with opposing effects. The application of the random and bias vorticity error predictions are illustrated by application to experimental velocity data determined using cross-correlation DPIV (CCDPIV) analysis of digital images of a laminar vortex ring. Received: 31 October 1997/Accepted: 6 February 1998     

Natural Convection in a Cavity Filled with Porous Layers on the Top and Bottom Walls

Natural convection in a partially filled porous square cavity is numerically investigated using SIMPLEC method. The Brinkman-Forchheimer extended model was used to govern the flow in the porous medium region. At the porous-fluid interface, the flow boundary condition imposed is a shear stress jump, which includes both the viscous and inertial effects, together with a continuity of normal stress. The thermal boundary condition is continuity of temperature and heat flux. The results are presented with flow configurations and isotherms, local and average Nusselt number along the cold wall for different Darcy numbers from 10⁻¹ to 10⁻⁶, porosity values from 0.2 to 0.8, Rayleigh numbers from 10³ to 10⁷, and the ratio of porous layer thickness to cavity height from 0 to 0.50. The flow pattern inside the cavity is affected with these parameters and hence the local and global heat transfer. A modified Darcy–Rayleigh number is proposed for the heat convection intensity in porous/fluid filled domains. When its value is less than unit, global heat transfer keeps unchanged. The interfacial stress jump coefficients β ₁ and β ₂ were varied from  −1 to +1, and their effects on the local and average Nusselt numbers, velocity and temperature profiles in the mid-width of the cavity are investigated.     

Computation of conjugate heat transfer in the turbulent mixed convection regime in a vertical channel with multiple heat sources

This paper presents the results of a comprehensive numerical study to analyze conjugate, turbulent mixed convection heat transfer from a vertical channel with four heat sources, uniformly flush-mounted to one of the channel walls. The results are presented to study the effect of various parameters like thermal conductivity of wall material (k _s), thermal conductivity of flush-mounted discrete heat source (k _c), Reynolds number of fluid flow (Re _s), modified Richardson number (Ri ⁺) and aspect ratio (AR) of the channel. The standard k-ε turbulence model, modified by including buoyancy effects with physical boundary conditions, i.e. without wall functions, has been used for the analysis. Semi-staggered, non-uniform grids are used to discretise the two dimensional governing equations, using finite volume method. A correlation, encompassing a wide range of parameters, is developed for the non-dimensional maximum temperature (T ^*) using the asymptotic computational fluid dynamics (ACFD) technique.     

Drag reduction of motor vehicles by active flow control using the Coanda effect

A test facility has been constructed to realistically simulate the flow around a two dimensional car shaped body in a wind tunnel. A moving belt simulator has been employed to generate the relative motion between model and ground. In a first step, the aerodynamic coefficients c _L and c _D of the model are determined using static pressure and force measurements. LDA-measurements behind the model show the large vortex and turbulence structures of the near and far wake. In a second step, the ambient flow around the model is modified by way of an active flow control which uses the Coanda effect, whereby the base-pressure increases by nearly 50% and the total drag can be reduced by 10%. The recirculating region is completely eliminated. The current work reveals the fundamental physical phenomena of the new method by observing the pressure forces on the model surface as well as the time averaged velocities and turbulence distributions for the near and far wake. A theory resting on this empirical information is developed and provides information about the effectiveness of the blowing method. For this, momentum and energy equations were applied to the flow around the vehicle to enable a validation of the theoretical results using experimental values. Received: 9 June 1998 / Accepted: 20 July 1999     

Numerical simulation of the flow around a porous covering square cylinder in a channel via lattice Boltzmann method

In this paper, a detailed investigation on the flow past a porous covering cylinder is presented through the lattice Boltzmann method. The Brinkman‐Forchheimer‐extended Darcy model is adopted for the entire flow field with the solid, fluid, and porous medium. The effects of several parameters, such as porous layer thickness, Darcy number, porosity, and Reynolds number on flow field are discussed. Compared with the case of a solid cylinder, the present work shows that the porous layer may play an important role on the flow, the lift and drag force exerted on the cylinder. The numerical results indicate that the maximal drag coefficient C_d and maximal amplitude of lift coefficient C_l exist at certain Darcy number which is in the range of 10^?6–10^?2. Copyright © 2010 John Wiley & Sons, Ltd.     

加权最小二乘法求解Rayleigh阻尼系数的讨论

合理确定Rayleigh阻尼矩阵比例系数对于准确计算场地地震响应有重要影响。本文提出以H2iH′2iA2i为权重系数,采用加权最小二乘法求解Rayleigh阻尼矩阵比例系数,以提高土层地震反应时程分析的计算精度。以上海某处地质剖面建立深覆盖土层有限元模型,计算其在28条各具代表性的基岩场地地震波激励下的动力反应,与其他研究成果对比分析,验证本文方法的精度及适用性。     

Model representations of the interaction between the solar wind and the supersonic interstellar medium flow. prediction and interpretation of experimental data

The planning and conducting of physical experiments requires the development of theoretical models capable either of predicting possible experimental data or explaining those already obtained. The processes taking place in the physical world can be understood only in terms of the close interaction between theory and experiment. Developing any quantitative or qualitative model of a physical phenomenon requires a mathematical apparatus, on the basis of which such models can be constructed. The branch of theoretical science using the methods of magnetohydrodynamics and hydroaeromechanics for studying space physics problems is usually called cosmic gasdynamics; it is mostly used in developing models of physical phenomena occurring under space conditions. In order to emphasize the importance of cosmic gasdynamics in the development of astrophysics and space research, we will present several examples of models constructed by aerodynamicists. These models not only played an important role in qualitative predictions but are still being developed due to the need for the quantitative interpretation of the experimental data. The solar corona was long thought to be a formation in a state of gravitational equilibrium (Chapman model). However, it turned out that the pressure at infinity obtained on the basis of this equilibrium solution is considerably greater than the estimated pressure in the interstellar gas surrounding the solar corona. In [1] it was concluded that in this case the solar corona gas must expand and a solution describing this expansion was obtained by invoking the steady-state hydrodynamics equations in the spherically-symmetric approximation. The solution of these equations led to the theoretical prediction of the solar wind, a radial flow of fully ionized hydrogen plasma issuing from the solar corona at a low subsonic velocity but already hypersonic at the Earth’s orbit. Subsonic-to-supersonic transition is ensured by solar gravitation which in this case plays the role of a convergent-divergent nozzle. Within a year, the theoretical prediction of the solar wind [1] was confirmed by its experimental detection [2] onboard the Soviet spacecraft Luna-2. It turned out that at the Earth’s orbit the mean velocity of the solar wind V _E ≈ 450 km·s^?1, the mean proton temperature T _E ≈ 6 · 10⁴ K (the electron temperature is somewhat higher), and the mean concentration of protons (and electrons) n _E ≈ 10 cm^?3. The first hydrodynamic model of the supersonic solar-wind flow past the Earth’s magnetosphere [3] was only qualitative, since it considered a flow past a plane magnetic dipole in the approximation of a thin layer between the bow shock and an “obstacle” embedded in the flow. However, it was constructed before the actual discovery of the solar wind and provided further important impetus to the development of models of the supersonic solar wind flow past planets with a detached shock. One more example is furnished by the gasdynamicmodel of the solar wind flow past cometary atmospheres, first suggested in In this work, a model of the interaction between the supersonic solar wind and the supersonic flow of the local, i.e., surrounding the Sun, interstellar medium is considered; it was first suggested in [6] in a much simplified formulation. This model has been actively developed in connection with the flights of the spacecraft Voyager 1 and 2, Ulysses, Hubble Space Telescope, SOHO, and others, exploring the outer regions of the solar system.