Wavelet method for solving the unsteady porous-medium flow problem with discontinuous coefficients

A method is proposed for constructing stable approximate wavelet decompositions of weak solutions to boundary value problems for the unsteady porous-medium flow equation with discontinuous coefficients and inexact data. The method is based on the general scheme for finite-dimensional approximation in Tikhonov regularization and on multiresolution analysis with basis functions defined as the product of one-dimensional Daubechies wavelets.     

Transition to chaos in lid-driven square cavity flow

To date, there are very few studies on the transition beyond second Hopf bifurcation in a lid-driven square cavity, due to the difficulties in theoretical analysis and numerical simulations. In this paper, we study the characteristics of the third Hopf bifurcation in a driven square cavity by applying a consistent fourth-order compact finite difference scheme rectently developed by us. We numerically identify the critical Reynolds number of the third Hopf bifurcation located in the interval of (13944.7021,13946.5333) by the method of bisection. Through Fourier analysis, it is discovered that the flow becomes chaotic with a characteristic of period-doubling bifurcation when the Reynolds number is beyond the third bifurcation critical interval. Nonlinear time series analysis further ascertains the flow chaotic behaviors via the phase diagram, Kolmogorov entropy and maximal Lyapunov exponent. The phase diagram changes interestingly from a closed curve with self-intersection to an unclosed curve and the attractor eventually becomes strange when the flow becomes chaotic.     

明渠非定常水流的垂向结构实验研究

本文报告了明渠非定常水流垂向结构的实验研究,通过 LDV 激光测速,我们给出了非定常水流的速度削面,湍流度,Eulerian质量输运速度的垂向分布。实验结果表明:在非定常周期很长时非定常流动的速度剖面,湍流强度及其粗糙度的影响基本上与定常流动相一致,当非定常流动为纯往复流动时,存在一个逆时针方向的环流,这种环流有利于底部的质量输运。     

Measurements in the vaneless diffuser of a radial flow compressor

Results from an experimental study of flow behaviour at the inlet of a vaneless diffuser of a centrifugal compressor are presented. Measurements from a crossed hot-wire probe are given for operating points having inlet flow coefficients ranging from 0.006 to 0.019 at different Reynolds numbers. Instantaneous, time-averaged, and phase-averaged absolute velocity and flow angle at the diffuser inlet are deduced from the hot-wire signals after correction for mean density variations. These results show how flow behaviour varies in stable, rotating stall and surge regimes of compressor operation     

Characteristic Regimes of Transitional Separation Bubbles in Unsteady Flow

This experimental investigation deals with transition phenomena of a separated boundary layer under unsteady inlet flow conditions. The main purpose of this investigation is to understand the influence of the rotor-stator interaction in turbomachinery on the subsequent, highly loaded boundary layer. The research project is divided into two phases. In the first phase, which has been completed recently, only the variation of mean velocity caused by upstream blades was simulated in the experiments while the free-stream turbulence intensity was retained at a constant low level. The experiments are carried out in an Eifel-type wind tunnel to investigate the laminar separated boundary layer of a flat plate under oscillating inlet conditions. The adverse pressure gradient, similar to that of turbomachines, is generated by the contoured upper wall. The unsteadiness is produced by a rotating flap located downstream of the test section. The reduced frequency, the amplitude and the mean Reynolds number are varied to simulate the conditions prevailing in turbomachines. In addition to the Kelvin–Helmholtz instability of the separated shear layer, a lower frequency instability was observed. This is frequently referred to as `free shear layer flapping' and results in two distinctly different ways of re-attachment, depending primarily on the Reynolds number. For low momentum thickness Reynolds numbers at the separation point, large-scale vortices locked to the frequency of the unsteady main flow are identified. They originate nearly at the top of the separation bubble and are ejected downstream. A fully turbulent boundary layer develops after these vortices mix out. For higher Reynolds numbers, transition is completed within a short length of the free shear layer and there-attachment region. The characteristic momentum thickness Reynolds number separating these two regimes in unsteady flow is about 125. The Strouhal number (reduced frequency) does not appear to have any significant effect. Based on the experimental results, this behaviour is discussed in some detail. This revised version was published online in July 2006 with corrections to the Cover Date.     

Maxwell流体管内起动流的研究

对Maxwell流体流体管内起动流的振荡特性进行研究，得到了描述振荡特性的解析解。研究了黏弹性参数对各时刻速度剖面的影响，获得了轴心速度，平均速度和壁面摩擦力随时间的变化规律以及它们的频率特征。结果表明振荡的基频成分决定了流动的主要特性，给出了并分析了基频频率与振幅和黏弹性参数之间的关系。     

动失速型非定常分离涡结构的控制

本文对动失速型非定常分离涡结构的控制方法，在低速风洞中应用相平均测压技术进行了实验研究。在二元平板模型中部安装一作俯仰振荡的扰流板产生动失速型分离涡，在其上游安装另一用作控制的小扰流板。实验结果表明，应用前置的振荡小扰流板可影响并改变动失速分离涡的强度和对流特性。在最有利的控制相位下，涡吸力峰可降低４８％，涡对流时间可以推迟０．１９周期。对于间歇式振荡扰流板，采用相位提前控制方式比相位滞后控制方式更有效。     

跨音轴流压气机转子叶尖喷气扩稳机理分析

对跨音压气机Rotor35进行了多通道全三维定常/非定常叶顶喷气数值模拟。数值计算所获得实壁机匣总性能与试验结果符合良好。计算表明使用3.6%转子堵塞流量的叶顶喷气量可以获得21.4%的扩稳效果。定常计算结果显示叶顶喷气重点影响0.9叶展以上区域,使得该区域进气攻角和扩散因子减小,从而降低叶顶载荷,减小了由激波和泄漏涡相互作用形成的通道堵塞。非定常计算结果显示,叶尖喷气的扩稳效果来自两方面:一是对某一叶片叶顶的卸载作用;二是对激波/泄漏涡干扰形成的低能区重新注入轴向动量。后者对通道流通的改善作用大于前者。非常高的喷射频率使得叶顶喷气能够抑制每个通道中低速区的进一步增长,从而实现了对压气机的扩稳。     

Unsteady CFD computations using vertex‐centered finite volumes for unstructured grids on Graphics Processing Units

This paper presents a Navier–Stokes solver for steady and unsteady turbulent flows on unstructured/hybrid grids, with triangular and quadrilateral elements, which was implemented to run on Graphics Processing Units (GPUs). The paper focuses on programming issues for efficiently porting the CPU code to the GPU, using the CUDA language. Compared with cell‐centered schemes, the use of a vertex‐centered finite volume scheme on unstructured grids increases the programming complexity since the number of nodes connected by edge to any other node might vary a lot. Thus, delicate GPU memory handling is absolutely necessary in order to maximize the speed‐up of the GPU implementation with respect to the Fortran code running on a single CPU core. The developed GPU‐enabled code is used to numerically study steady and unsteady flows around the supercritical airfoil OAT15A, by laying emphasis on the transonic buffet phenomenon. The computations were carried out on NVIDIA's Ge‐Force GTX 285 graphics cards and speed‐ups up to ～46 × (on a single GPU, with double precision arithmetic) are reported. Copyright © 2010 John Wiley & Sons, Ltd.     

Cell face velocity alternatives in a structured colocated grid for the unsteady Navier–Stokes equations

The use of a colocated variable arrangement for the numerical solution of fluid flow is becoming more and more popular due to its coding simplicity. The inherent decoupling of the pressure and velocity fields in this arrangement can be handled via a special interpolation procedure for the calculation of the cell face velocity named pressure‐weighted interpolation method (PWIM) (AIAA J. 1983; 21 (11):1525–1532). In this paper a discussion on the alternatives to extend PWIM to unsteady flows is presented along with a very simple criterion to ascertain if a given interpolation practice will produce steady results that are relaxation dependent or time step dependent. Following this criterion it will be shown that some prior schemes presented as time step independent are actually not, although by using special interpolations can be readily adapted to be. A systematic way of deriving different cell face velocity expressions will be presented and new formulae free of Δt dependence will be derived. Several computational exercises will accompany the theoretical discussion to support our claims. Copyright © 2009 John Wiley & Sons, Ltd.