Optimal control of cylinder wake flow by electro-magnetic force based on adjoint flow field

A nonlinear adjoint-based optimal control approach of cylinder wake flow by using electro-magnetic forcing has been proposed and investigated numerically in the paper. A cost functional representing the balance between the enstrophy Ω² and the interaction parameter N has been developed, and its corresponding adjoint equations have been derived. The sensitivity of the cost functional is found to be a simple function of the adjoint stream function in the adjoint field. Under the action of optimal force, N(t), the flow separation is suppressed successfully, the oscillations of drag and lift disappear and the total drag coefficient decreases dramatically. Furthermore, the global enstrophy is proved to be equal to the total squared strain during the control, which validates the conservation rule of Okubo-Weiss function.     

Two-phase damping for internal flow: Physical mechanism and effect of excitation parameters

Two-phase flow induced-vibration is a major concern for the nuclear industry. This paper provides experimental data on two-phase damping that is crucial to predict vibration effects in steam generators. An original test section consisting of a tube subjected to internal two-phase flow was built. The tube is supported by linear bearings and compression springs allowing it to slide in the direction transverse to the flow. An excitation system provides external sinusoidal force. The frequency and magnitude of the force are controlled through extension springs. Damping is extracted from the frequency response function of the system. It is found that two-phase damping depends on flow pattern and is fairly proportional to volumetric fraction for bubbly flow. Measurements are completed by the processing of high-speed videos which allow to characterize the transverse relative motion of the gas phase with respect to the tube for bubbly flow. It is shown that the bubble drag forces play a significant role in the dissipation mechanism of two-phase damping.     

Analysis of the flow in inhomogeneous particle beds using the spatially averaged two-fluid equations

The drag force term appearing in two-fluid models for fluid–particle flows is commonly closed by expressing it as a function of the local quantities, such as the local particle volume fraction, the local slip velocity between the particle and fluid phases, and the local mean-squared fluctuating velocity of the particles. The adequacy of such closures for inhomogeneous suspensions has been debated in the literature and some researchers have suggested the need for additional terms involving spatial gradients in these quantities. To test this proposition, simulations of flow in inhomogeneous steady beds of particles have been performed using the lattice-Boltzmann method. The particle beds consisted of disordered assemblies with a density profile on a scale much larger than the particle radius. Inhomogeneous beds with a controlled density profile were generated in three different ways, (i) by inhomogeneous stretching of the particle bed in one direction, (ii) by applying an inhomogeneous force to the particle phase during random motion of the particles, and (iii) by taking snapshots of a direct simulation of a traveling wave in a fluidization simulation. The global structure of the three beds was comparable, while assessment of the radial distribution functions showed that the three beds exhibited clearly different microscopic structures.     

Fluidelastic instability study in a rotated triangular tube array subject to two-phase cross-flow. Part I: Fluid force measurements and time delay extraction

Fluidelastic instability is a key issue in steam generator tube bundles subjected to cross-flow. The extension to two-phase flow of the existing theoretical models, developed and tested mostly for single phase flow, is investigated in this paper. The time delay is one of the key parameter for modeling fluidelastic instability, especially the damping controlled mechanism. The direct measurement of the time delay between the tube motion and the fluid force faces certain difficulties in two-phase flow since the high turbulence due to the interaction of the two components of the flow may increase the randomness of the measured force. To overcome this difficulty, an innovative method for extracting the time delay inherent to the quasi-steady model for fluidelastic instability is proposed in this study.Firstly, experimental measurements of unsteady and quasi-static fluid forces (in the lift direction) acting on a tube subjected to air–water two-phase flow were conducted. The unsteady fluid forces were measured by exciting the tube using a linear motor. These forces were measured for a wide range of void fractions, flow velocities and excitation frequencies. The experimental results showed that the unsteady fluid forces could be represented as single valued function of the reduced flow velocity. It was also found that for a given frequency, the unsteady fluid force phase was weakly dependent on the void fraction for the range of flow velocities considered.The time delay was determined by equating the unsteady fluid forces with the quasi-steady forces. The results given by this innovative method of measuring the time delay in two-phase flow were consistent with theoretical expectations. The time delay could be expressed as a linear function of the convection time and the time delay parameter was determined for void fractions ranging from 60% to 90%.     

Eulerian spectrum of finite-time Lyapunov exponents in compound channels

Fluid flows reveal a wealth of structures, such as vortices and barriers to transport. Usually, either an Eulerian or a Lagrangian frame of reference is employed in order to detect such features of the flow. However, the two frameworks detect structures that have different properties. Indeed, common Eulerian diagnostics (Hua-Klein and Okubo-Weiss criterion) employed in order to detect vortices do not always agree with Lagrangian diagnostics such as finite-time Lyapunov exponents. Besides, the former are Galilean-invariant whereas the latter is objective. However, both the Lagrangian and the Eulerian approaches to coherent structure detection must show some links under any inertial-frame. Compound channels flows have been accurately studied in the past, both from a Lagrangian and an Eulerian point of view. The features detected do not superimpose: Eulerian vortices do not coincide with barriers to transport. The missing link between the two approaches is here recovered thanks to a spectral analysis.

     

Drying performance and surface forces of a pair of impinging radial jets

The drying performance and the forces generated by a pair of radial jets impinging on a plane surface are studied experimentally and compared with a pair of in-line (axial) jets. The reattachment rings of the radial jets are visualized by liquid-crystal-thermography. The results of the drying experiments are presented in terms of Sherwoodnumbers. The results of the surface force measurements are presented in terms of force coefficients. The effects of flow exit angle, exit velocity, horizontal nozzle spacing and vertical nozzle to surface distance on mass transfer and net surface force are discussed. Results show that drying performance of a pair of radial jets with an exit angle towards the surface of 10° is 18% higher than that of a pair of in-line jets. The surface force of impinging radial jets depends on the exit angle. A global pressure force on the surface is typical for jets directed towards the surface, while a global lifting force is found when the jet is initially directed away from the surface. Radial jets have a high potential in drying applications for force sensitive products and when high drying rates and/or small suction forces are of interest.     

Experimental analysis of flow structures and forces on a 3D-bluff-body in constant cross-wind

The effect of cross-wind on a bluff-body in ground proximity is analysed experimentally. A recently introduced wind-tunnel model is used to analyse the effect of flow structures on the global forces and global moments. This study aims at creating a large database of the flow produced in the proximity of this bluff-body and giving some new insight into the physical mechanisms responsible for the increase in force in cross-wind situations. The importance of the model’s lee side to the global yawing moment is put forward. Two vortices are identified and changes in the slope of the yawing moment versus yaw angle are attributed to their presence.     

翼型绕流的电磁力控制

将表面包覆电磁激活板的翼型,按一定的攻角,置于流动的弱电介质溶液中,电磁激活板可产生作用于流体的切向电磁力(Lorentz力),从而改变流体边界层的结构. 在转动水槽中,对翼型绕流及电磁力控制下的绕流形态进行了实验研究. 结果表明,未加电磁力时,前缘涡的脱落点是不确定的,与流场具体条件有关,而后缘涡仅在尖角处脱落. 前缘涡与后缘涡相互影响,并周期性的脱体,在尾部形成涡街. 施加电磁力后,当力的方向与流动方向相同时,可以在一定程度上抑制分离,消除涡街,其效果与减小攻角类似. 加反向电磁力时,则相当于加大攻角,在翼型体的背风面形成涡街.      

Control of the vortex flow around a cone with a spark discharge

A subsonic flow around a cone mounted at an angle of attack is considered. A separated flow is formed on the cone with formation of a pair of large-scale streamwise vortices arranged symmetrically or asymmetrically, depending on flow parameters. The possibility of controlling the vortex flow by means of an electric discharge is studied in experiments. It is demonstrated that the use of an electric discharge makes it possible both to transform an asymmetric flow to symmetric and to generate directed asymmetry. Controlling the flow asymmetry, in turn, allows the side force direction to be controlled.     

Instability of cylinder wake under open-loop active control

Instability of a wake controlled by a streamwise Lorentz force is investigated through a Floquet stability analysis. The streamwise Lorentz force, which is a two-dimensional control input created by an electromagnetic actuator located on the cylinder surface,adjusts the base flow to affect the three-dimensional wake instability and achieve wake stabilization and transition delay. The instability mode at a Reynolds number Re = 300 can be transformed from B to A with N = 1.0, where N is an interaction number representing the strength of the Lorentz force relative to the inertial force in the fluid. The wake flow is Floquet stable when N increases to 1.3. The spanwise perturbation wavelengths are 3.926 D and 0.822 D in the modes A and B, respectively, where D is the cylinder diameter. In addition, the oscillating amplitudes of drag and lift are reduced with the increase in the interaction number. Particle tracing is used to explore the essential physical mechanism for mode transformation. The path lines show that suppression of flow separation hinders the fluid deformation and rotation, leading to the decrease in elliptic and hyperbolic instability regions, which is the material cause of mode transformation.All of the results indicate that wake stabilization and transition delay can be achieved under open-loop active control via the streamwise Lorentz force.     

考虑等效曲率的超二次曲面单元非线性接触模型

基于连续函数包络的超二次曲面单元可有效地描述自然界和工业生产中的非球体颗粒形态, 并通过非线性迭代方法精确计算单元间的接触力. 对于具有复杂几何形态的超二次曲面单元, 线性接触模型不能准确地计算不同接触模式下的作用力. 考虑超二次曲面单元相互作用时不同颗粒形状及表面曲率的影响, 本文发展了相应的非线性黏弹性接触模型. 该模型将不同接触模式下的法向刚度和黏滞力统一表述为单元间局部接触点处等效曲率半径的函数; 切向接触作用则借鉴基于Mohr-Coulomb摩擦定律的球体单元非线性接触模型的计算方法. 为检验超二次曲面单元接触模型的可靠性, 对球形颗粒间的法向碰撞、椭球体颗粒间的斜冲击过程、圆柱体的静态堆积和椭球体的动态卸料过程进行离散元模拟, 并与有限元数值结果及试验结果进行对比验证. 计算表明, 考虑接触点处等效曲率半径的超二次曲面非线性接触模型可准确地计算单元间的接触碰撞作用, 并合理地反映非球形颗粒体系的运动规律. 在此基础上进一步分析了不同长宽比和表面尖锐度对卸料过程中颗粒流动特性的影响, 为非球形颗粒材料的流动特性分析提供了一种有效的离散元方法.      

FBRE系统偏心力下的功率流特性及主动控制

利用机械阻抗方法推导了偏心力作用下，机器－隔振器－柔性板基础所组成的FBRE隔振系统的动态传递方程。以传递到柔性板基础的功率流为目标函数，通过计算机仿真计算，研究了在对称系统中偏心力作用下的功率流特性，分析了和对称力、力矩作用效果的不同和联系，并探讨了主动控制策略的效果。     

预处理方法在低速弹箭流场计算中的应用

基于Weiss-Smith预处理矩阵和全局截断预处理参数,采用有限体积方法对雷诺平均Navier-Stokes方程进行离散。对流项离散采用二阶线性重构和AUSM ＋-up格式,时间推进方法采用多重网格下的LU-SGS方法。结合M PI消息传递方法,建立了一套计算低速流动的并行数值方法。计算了低速椭球体的流场和气动力,压力系数和切应力系数计算结果与文献实验结果对比吻合度较好。生成了末敏弹的流场计算网格,对绕末敏弹流场进行了数值模拟。对多重网格下多进程的加速比和并行效率进行了测试,显示了程序良好的并行效率。计算的气动力结果与实验结果吻合。综合结果表明：本文的数值方法能够用于低速弹箭流场和气动力计算,为新型弹箭的设计和定型提供保证。     

双曲冷却塔塔筒水平地震响应分析

为明确冷却塔在水平地震下的内力环向分布特征及内在原因,同时探究不同地震波时程与规范反应谱所得内力差异的原因,以某大型双曲冷却塔为例,在动力特性分析的基础上,通过反应谱方法和时程方法的水平地震响应计算及对比分析,对上述两个问题进行了研究。研究表明:由于仅侧弯振型对水平地震有贡献,而塔筒的侧弯振型和实际响应均表现为整体侧倾并伴随微弱的截面“流动”变形,这也使塔筒各内力的环向分布分别呈现正弦、余弦分布特征;其整体侧倾可类比于悬臂杆结构,塔筒子午向轴力FY、子午向弯矩MY、剪力FXY和扭矩MXY的环向分布可借助悬臂杆侧倾时截面正应力和剪应力分布来解释;而截面“流动”变形则决定了环向轴力FX和环向弯矩MX的环向分布;整体侧移显著而截面变形极小也使FY和FXY的幅值在塔筒中下部明显大于FX;由于冷却塔第1阶侧弯振型在水平地震响应中往往起绝对主导作用,因此可先对所选地震波计算得到其反应谱,对比第1阶侧弯振型周期对应的水平地震影响系数α值,即可初步推断不同时程及规范反应谱方法所得结果的大小关系。     

?????????????????????????

基于计算结构力学理论与递推公式相结合的方法，推导出随机刚架结构在退化成机构 前各阶段的内力与结构处于初始状态时的内力关系. 因此，在基于随机有限元的可靠性分 析过程中，不必反复修改刚阵和进行反向节点力的增加，在极限状态函数对设计变量求导 过程中，导数形式只用到了原结构总刚度矩阵和原结构外载向量，这对编写电算程序是有 益的. 应用所提的方法与其它基于随机有限元的可靠性分析方法所求得的结果是一致 的，且大大提高了计算效率.     

包覆电磁场激活板的圆柱尾迹的数值研究

利用电磁场作用于电介质溶液的Lorentz力可以控制溶液的流动.对置于弱电介质溶液中,包覆电磁场激活板的圆柱周围的绕流进行数值研究,讨论了电磁场激活板的安装位置,激活方式等对圆柱尾迹的影响.     

Interactions of a propeller with a stator-induced circumferentially varying flow

The interactions of a circumferentially varying stator cascade and a downstream fixed pitch propeller were investigated experimentally. The global performance of the components and the coupled system were systematically investigated through force and moment measurements on the propulsor model in a water tunnel. In addition, the wake of the cyclic stator cascade with and without the propeller was investigated downstream from a propulsor model using the Stereoscopic PIV technique. A cyclic distribution of the stators’ deflections resulted in non-axisymmetric distributions of the flow field downstream of the stator array. The stator distribution alone produced a significant side force that increased linearly with stator pitch amplitude. When a propeller was incorporated downstream from the cyclic cascade, the side force from the stator cascade was reduced, but a small normal force and pitching moment were created. The generation of these secondary forces and moments can be related to the redistribution of the tangential flow from the cyclic cascade into the axial direction by the retreating and advancing blade states of the fixed pitch propeller.     

Turbulent boundary layer separation control and loss evaluation of low profile vortex generators

The present paper analyses the results of a detailed experimental study on low profile vortex generators used to control the turbulent boundary layer separation on a large-scale flat plate with a prescribed adverse pressure gradient, typical of aggressive turbine intermediate ducts. This activity is part of a joint European research program on Aggressive Intermediate Duct Aerodynamics (AIDA). Laser Doppler Velocimetry and a Kiel total pressure probe have been employed to perform measurements in the test section symmetry plane and in cross-stream planes to investigate the turbulent boundary layer, with and without control device application.Velocity fields, Reynolds stresses, and total pressure distributions are analysed and compared for the controlled and non controlled flow conditions to characterize the mean flow behaviour. The detail and the accuracy of the measurements allow the evaluation of the deformation works of the mean motion in the test section symmetry plane. Normal and shear contributions of viscous and turbulent deformation works have been analysed and employed to explain the distribution of the total pressure loss. For the controlled flow the discussion of the flow field is extended considering the effects of the vortex generated in the cross-stream planes. The experimental data allow the evaluation of the global amount of losses, considering a balance of total pressure fluxes in the different measuring planes.     

Global modes in a confined impinging jet: application to heat transfer and control

We investigate the stability and control of a plane, laminar jet impinging on a flat plate in a channel, a geometry used to cool down a hot wall with a cold air jet in many industrial configurations. The global stability analysis indicates that, even for a strong confinement, the two-dimensional (2-D) steady flow is unstable to three-dimensional (3-D), steady perturbations. In the simplest limit case where dilatation effects are neglected, we show that the development of the instability induces a significant spanwise modulation of the heat flux at the impacted wall. To control the leading global mode, we propose adjoint-based 3-D harmonic and 2-D steady forcing in the bulk or at the wall. We show for instance that the unstable mode is controllable using a spanwise uniform blowing at the upper wall, in a specific domain corresponding to the footprint of the upper recirculating bubble. These techniques are applied to a novel open-loop control, in which we introduce into the flow a small airfoil, modelled by the lift force it exerts on the flow.     

A material force method for inelastic fracture mechanics

A material force method is proposed for evaluating the energy release rate and work rate of dissipation for fracture in inelastic materials. The inelastic material response is characterized by an internal variable model with an explicitly defined free energy density and dissipation potential. Expressions for the global material and dissipation forces are obtained from a global balance of energy-momentum that incorporates dissipation from inelastic material behavior. It is shown that in the special case of steady-state growth, the global dissipation force equals the work rate of dissipation, and the global material force and J-integral methods are equivalent. For implementation in finite element computations, an equivalent domain expression of the global material force is developed from the weak form of the energy-momentum balance. The method is applied to model problems of cohesive fracture in a remote K-field for viscoelasticity and elastoplasticity. The viscoelastic problem is used to compare various element discretizations in combination with different schemes for computing strain gradients. For the elastoplastic problem, the effects of cohesive and bulk properties on the plastic dissipation are examined using calculations of the global dissipation force.