Finite element solutions for laminar and turbulent compressible flow

The time-split finite element method is extended to compute laminar and turbulent flows with and without separation. The examples considered are the flows past trailing edges of a flat plate and a backward-facing step. Eddy viscosity models are used to represent effects of turbulence. It is found that the time-split method produces results in agreement with previous experimental and computational results. The eddy viscosity models employed are found to give accurate predictions in all regions of flow except downstream of reattachment.     

Experimental study on the confined 3D laminar flow past a square prism with a high blockage ratio

An experimental PIV study is presented that addresses the confined 3D laminar flow behaviour past a square prism. The Reynolds number (Re), based on prism cross-section height varies between 100 and 256. The channel aspect ratio is 1/1 and the blockage ratio is 1/2.5. This geometry is representative of a passive method to enhance mixing in otherwise highly ordered laminar channel flow. It is found that the lateral walls exert a strong effect on the flow behaviour with two main consequences: (a) the onset of vortex shedding is delayed to a Re in the vicinity of 170, as opposed to the unconfined case where the critical Re is reported to be between 50 and 60 and (b) transition from the steady closed recirculation bubble regime to the vortex shedding regime is not abrupt. In particular, there is a range of Re for which the closed recirculation bubble pulsates with increasing amplitude prior to the onset of the Karman street regime. The experimental results are supported by numerical computations, and the relation of the results with the practical design of engineering systems is also discussed.     

Experimental study of turbulent round jet flow impinging on a square cylinder laid on a flat plate

A turbulent axisymmetric air jet impinging on a square cylinder mounted on a flat plate has been studied experimentally. Turbulence statistics and flow’s topology were investigated. When the surface was heated through uniform heat flux, local heat transfer coefficient was measured. The jet from a long round pipe, 75 pipe diameters (D) in length, at Reynolds number of 23,000, impinged vertically on the square cylinder (3D × 3D × 43D). Measurements were performed using particle image velocimetry, flow visualization using fluorescent dye and infrared thermography. The flow’s topology demonstrated a three-dimensional recirculation after separating from the square cylinder and a presence of foci between the bottom corner and the recirculation’s detachment line. The distribution of heat transfer coefficient was explained by the influence of these flow’s structures and the advection of kinetic energy. On the impingement wall of the square cylinder, a secondary peak in heat transfer coefficient was observed. Its origin can be attributed to very pronounced shear production coupled with the external turbulence coming from the free jet.     

Finite volume methods for laminar and turbulent flows using a penalty function approach

A penalty function, finite volume method is described for two-dimensional laminar and turbulent flows. Turbulence is modelled using the k-? model. The governing equations are discretized and the resulting algebraic equations are solved using both sequential and coupled methods. The performance of these methods is gauged with reference to a tuned SIMPLE-C algorithm. Flows considered are a square cavity with a sliding top, a plane channel flow, a plane jet impingement and a plane channel with a sudden expansion. A sequential method is employed, which uses a variety of dicretization practices, but is found to be extremely slow to converge; a coupled method, evaluated using a variety of matrix solvers, converges rapidly but, relative to the sequential approach, requires larger memory.     

Experimental and numerical study on a laminar fluid-structure interaction reference test case

With the rapid development of numerical codes for fluid-structure interaction computations, the demand for validation test cases increases. In this paper we present a comparison between numerical and experimental results for such a fluid-structure interaction reference test case. The investigated structural model consists of an aluminum front cylinder with an attached thin metal plate and a rear mass at the trailing edge. All the structure is free to rotate around the axle mounted in the center of the front cylinder. The model's geometry and mechanical properties are chosen in such a way as to attain a self-exciting periodical swiveling movement when exposed to a uniform laminar flow. Reproducibility of the coupled fluid-structure motion is the key criterion for the selection of the model in order to permit an accurate reconstruction of the results in the time-phase space. The Reynolds number of the tests varies up to 270 and within that range the structure undergoes large deformations and shows a strong nonlinear behavior. It also presents two different self-excitation mechanisms depending on the flow velocity. Hence, challenging tasks arise for both the numerical solution algorithm and the experimental measurements. To account for the two different excitation mechanisms observed on increasing the speed of the flow, results for two different velocities are considered: the first at 1.07 m/s (Re=140) and the second at 1.45 m/s (Re=195). The comparisons presented in this paper are carried out on the basis of the time trace of the front body angle, trailing edge coordinates, structure deformation and the time-phase resolved flow velocity field. They reveal very good agreement in some of the fluid-structure interaction modes whereas in others deficiencies are observed that need to be analyzed in more detail.     

The turbulent character and pressure loss produced by periodic symmetric ribs in a circular duct

The structural character and steady-state statistics of the turbulence inside a rib-wall circular duct is investigated by the large-eddy simulation (LES) methodology. The impetus of this study is to gain an understanding of the principle physics attributing to minimizing the pressure recovered (or maximizing the pressure loss) within the core flow. For a rib periodicity with height (h) to pitch (p) ratio p/h=5, the computational results show that the majority of turbulence produced due to the rib’s presence is concentrated near the rib crest leading edge. Pairs of counter-rotating streamwise vortices form soon after the leading edge that are quickly convected radially toward the core flow. The turbulent activity within the duct trough region is negligible compared to the turbulence levels of the core flow. At this rib periodicity, the separated shear layers from the trailing edge of each rib nearly reattach to the trough floor before reaching the next rib. The resultant irrecoverable pressure loss in the form a centerline frictional coefficient is verified by an ‘at-sea’ test on board a US Navy submarine. Based on the duct diameter, their Reynolds numbers are Re_D^LES=8×10³ and (Re_D^exp)_avg=4×10⁶, respectively.     

超临界CO2直旋混合射流破岩特性的实验研究

超临界二氧化碳(CO₂)射流破岩既能降低岩石门限压力又能有效保护储层,直旋混合射流兼具直射流和旋转射流特点可提高破岩效率,基于此提出了超临界CO₂直旋混合射流的破岩方法。为了揭示超临界CO₂直旋混合射流破岩特性,设计加工出叶轮式直旋混合射流喷嘴,通过岩石定点冲击破碎实验对比了该射流与常规水射流的破岩效果,并研究了叶轮长度、叶轮中心孔直径、混合腔长度、喷射距离、射流压力等重要参数对超临界CO₂直旋混合射流破岩效果的影响。结果表明:相同实验条件下,该射流方法的平均破岩能力比常规水射流提高了42.9%;超临界CO₂直旋混合射流破岩易出现较大体积岩屑崩落现象;随着叶轮长度、混合腔长度、喷射距离的增大破岩效果均先增强后减弱,实验条件下上述参数存在最优范围值;叶轮中心孔直径的增大会导致岩石破碎孔深度增加、直径减小;随着射流压力的升高,超临界CO₂直旋混合射流破岩效果有着较为明显的提升。研究结果可为超临界CO₂直旋混合射流破岩方法的进一步研究提供实验依据。     

Numerical study on the stability and mixing of vertical round buoyant jet in shallow water

IntroductionAsimplewaytodisposeofthelargequantitiesofwasteheatresultingfromsteamelectricpowergeneration ,istodischargetheheatedcondenserwaterthroughasubmergedroundoutfalllocatedatthebottomofthereceivingwater.Inrecentyears,suchkindsofhighvelocityinjectionsystemshaveincreasinglybeenusedaseffectivemixingdevices,frequently ,thesedevicesarelocatedinwatersofonlyafewportdiametersdeep ,anddischargingareaisrelativelysmall.Thusanunderstandingoftheinducedexcesstemperaturedistribution ,inrelationshiptothe…     

Experimental study on the wake behind tapered circular cylinders

The flow around tapered cylinders can act as basic models for numerous bluff body flows with a spanwise variation of either the body shape or the inflow conditions. The well-known vortex street is influenced by strong three-dimensional effects from the spanwise variation of the shedding frequency, namely oblique vortex shedding and vortex dislocations. Stereo-PIV was chosen to study these phenomena, since it allows analyzing planes with the full three-component, instantaneous velocity fields and local, time-dependent variations in the same setting. Hence, detailed aspects of the vortex dislocation phenomenon are presented. Single vortex dislocation events are presented through the local variation of the three measured velocity components u, v and w. Longer time-series reveal both period and location of these dislocation events, as well as quantity and sizes of the cells of constant shedding velocity in between them. The influence of the Reynolds number and the cylinder aspect ratio on the vortex cells could be shown. The analysis of the vortex shedding behavior shows good agreement with previously published results. At the same time, the applied PIV technique provides more spatial information than point-based measurements and offers insight into a Reynolds number range that is currently out of reach of Direct Numerical Simulations.     

An experimental study on the efficiency of transpiration cooling in laminar and turbulent hypersonic flows

An experimental study on the efficiency of transpiration cooling in hypersonic laminar and turbulent flow regimes is carried out in the Hypersonic Windtunnel Cologne with a focus on the aerothermal problems downstream of the cooled model part. The model is made of a material of low thermal conductivity (PEEK) with an integrated probe of a porous material. The experimental setup allows the direct comparison of the thermal behavior of transpiration cooling to a well-defined and radiatively cooled reference surface. Experiments are performed at Mach number of 6 and two different Reynolds numbers. Air, argon and helium are used as coolants at various flow rates, in order to identify the influence of coolant medium on cooling efficiency. The cooling efficiency of air and argon is comparable. Helium provides significantly higher cooling efficiency at the same blowing ratio, i.e. same coolant mass flow rate. The experimental data shows that the efficiency of the transpiration cooling in turbulent flows is much lower than in laminar flow.     

Effect of an inhibitor on high-speed turbulent flames and the transition to detonation

The influence of an inhibitor (CF₃Br or Halon 1301) on the propagation of high-speed turbulent flames, quasi-detonations and the transition to detonation has been investigated for methane-air, propane-air and acetylene-air mixtures. The experiments are carried out in a 13 m tube (15 cm diameter) filled with regularly spaced orifice plates (blockage ratio of 0.39) to ensure rapid flame acceleration. In all cases, the addition of the inhibitor reduces the turbulent flame velocity and extinguishes the flame with sufficient inhibitor concentration (2.7% and 7.5% for methane-air and propane-air, respectively). For acetylene-air mixtures, the quasi-detonation speed is progressively reduced with increasing inhibitor concentration and eventually causes the failure of the quasi-detonation and transition back to a fast turbulent flame. The inhibitor also narrows the propagation limits in all cases. To elucidate the inhibition mechanism, detailed modelling of both the turbulent flame structure as well as the chemical kinetics are required.     

Effect of span length and temperature on the 3-D confined flow around a vortex promoter

This article presents a numerical study on the influence of span length and wall temperature on the 3-D flow pattern around a square section vortex promoter located inside a micro-channel in the low Reynolds number regime. The first objective of the work is to quantify the critical Reynolds number that defines the onset of vortex shedding and to identify the different regimes that appear as a function of the channel aspect ratio (span to height ratio). We found that the critical Reynolds number for the onset of the Karman street regime increases as the aspect ratio decreases. In particular, for the aspect ratio of 1/2 the critical Reynolds number is nearly six times the critical Reynolds number of the 2-D problem. An intermediate oscillating regime between the steady and the Karman street solutions was also found to exist within a rather wide range of Reynolds numbers for small channel aspect ratios. The second objective was to investigate the influence of the vortex promoter wall temperature on both vortex shedding and flow pattern. This has practical engineering implications because the working fluid considered in the article is water that has a viscosity that depends significantly on temperature and promotes a strong coupling between the momentum and energy equations that influences the system behaviour. Results indicate that high surface temperature on the prism promotes the onset of the Karman street, suggesting design guidelines for micro-channel based heat sinks that make use of vortex promoters.     

Experimental study on the flow regimes past a confined prism undergoing self-sustained oscillations

An experimental study based on Particle Image Velocimetry (PIV) is presented with the objective of studying the flow regimes that appear in the flow past a confined prism undergoing self-sustained oscillations at low Reynolds numbers (Re). The square-section prism, placed inside a 3D square cross-section vertical channel with a confinement ratio of 1/2.5, was tethered to the channel walls and, therefore, it was allowed to move freely transverse to the incoming flow. Re (based on the prism cross-section height) was varied in the range from 100 to 700. Three different prism to fluid density ratios (m¹) were considered: 0.56, 0.70, and 0.91. These two parameters, Re and m¹, were used to map the results obtained. In particular, it was found that five different regimes appear: (1) steady prism with steady recirculation bubble, (2) steady prism with unsteady vortex shedding wake, (3) large amplitude low frequency oscillating prism with unsteady vortex shedding wake, (4) small amplitude high frequency oscillating prism with unsteady vortex shedding wake, and (5) irregular/chaotic motion of both the prism and the wake. The PIV results and associated numerical simulations were used to analyze the different prism and wake states.     

Al2O3陶瓷复合靶抗长杆弹侵彻性能和机理实验研究

用DOP(depthofpenetration)实验研究了质量分数为 90 %的Al2 O3 陶瓷复合靶抗长杆弹侵彻性能。实验表明 :随陶瓷厚度的增加 ,以差分效益系数和靶平均阻力表征的陶瓷复合靶的抗长杆弹侵彻性能降低 ;盖板能增强陶瓷复合靶抗侵彻性能。对回收陶瓷和盖板的破坏特征分析表明 :长杆弹对陶瓷复合靶侵彻过程大致可分为初期侵彻和稳定侵彻两个阶段 ,但对薄靶和厚靶 ,两个阶段在整个侵彻过程中所占比重不同。     

Experimental study on turbulent features in the negative transport region of asymmetric plane channel flow

Turbulent features of streamwise and vertical components of velocity in the negative transport region of asymmetric plane channel flow have been studied experimentally in details. Experiments show that turbulent fluctuations in negative transport region are suppressed, and their probability distributions are far from Gaussian. Besides, the skewness factors attain their negative maxima at the position of the maximum mean velocity, whereas the flatness factors attain their positive maxima at the same position. The project supported by the National Natural Science Foundation of China (19872043)     

An investigation of thermally stratified turbulent channel flow with temperature oscillation on the bottom wall by large eddy simulation

Thermally stratified shear turbulent channel flow with temperature oscillation on the bottom wall of the channel is calculated to investigate the behavior of turbulent flow and heat transfer by use of large eddy simulation (LES) approach coupled with dynamic subgrid-scale (SGS) models. The objective of this study is to deal with the effect of the temperature oscillation on turbulent behavior of thermally stratified turbulent channel flow and to examine the effectiveness of the LES technique for predicting statistically unsteady turbulent flow driven by time-varying buoyancy force. To validate the present calculation, thermally stratified shear turbulent channel flow is computed and compared with available data obtained by direct numerical simulation (DNS), which confirm that the present approach can be used to predict thermally stratified turbulent channel flow satisfactorily. Further, to illustrate the effect of the temperature oscillation with different Richardson numbers and periods of the oscillation on turbulence characteristics, the phase-averaged mean value and fluctuation of the resolved velocities and temperature, and instantaneous velocity fluctuation structures are analyzed.