Wall-modeling for large-eddy simulation of flows around an axisymmetric body using the diffuse-interface immersed boundary method

A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.     

Motion of a slender body made of magnetizable composite in a traveling magnetic field

The motion of a slender body made of magnetizable composite in a channel, along which coils producing a heterogeneous “traveling” magnetic field are mounted, is investigated. The coil axes are vertical and lie in the same plane. A mathematical model of a slender body made of viscoelastic magnetizable material is proposed. The magnetic force is calculated from a formula used in ferrohydrodynamics of magnetic fluids with equilibrium magnetization. The problem of the motion of this body in a channel in a vertical plane under the action of the magnetic field produced in an experimental setup is numerically solved. The dependence of the body velocity on the coil switching frequency is calculated and the effect of different problem parameters on the form of this dependence is studied. The theoretical results are compared with the experimental data.     

柔性双尾鳍效率的测量和优化

建立了仿生拍动推进推力和效率的测量平台,通过电机的转动扭矩、转动角速度及推进力计算得到仿生推进系统的推进效率.对本课题组研制的组合仿生无人水下航行器中柔性双尾鳍模型进行了测量与分析,结果表明柔性尾鳍推进时小摆幅摆动运动效率要高于大摆幅运动;双尾鳍推进方式在小摆幅模式下无论在推进力还是在效率方面都要优于单尾鳍.最后根据模型实验数据,优化了双尾鳍推进的运动参数,并用于优化UUV性能.     

An investigation on the lift force of a wing pitching in dynamic stall for a comfort control vessel

On the basis of a comfort control system for ocean vessels, the control forces and moments in the form of lift forces from active wings are of important interest. In an ocean vessel comfort control system, active wings or fins are commonly used and constantly adjust their angles of attack to produce optimal sea-keeping conditions. The unsteady nature of the flow field around a wing, and the behaviour of the generated lift force must be understood in order to optimize the comfort control system. This paper presents experimental data on the flow past a pitching wing, paying particular attention to the lagging effects between the fluid dynamic lift force and the motion of the wing at large angles of attack as a function of peak angle of attack and reduced frequency of oscillation. The range of motion investigated has been chosen according to the applicability of a comfort control wing surface. Numerical data is also included to aid explanation on some of the witnessed phenomena.     

近空间飞行器的气动复合控制原理及研究进展

简单回顾与讨论了在飞行器穿越大气层时,使用气动复合控制方案的必要性与可能性.其中最复杂情况为发动机喷流推力与舵面空气动力共同使用所形成的复合控制.喷流与外流相撞引起的强干扰形成了十分复杂的干扰流场,文中介绍了复杂流动形成的原因、流场结构的特点以及干扰引起的流场改变影响了飞行器性能的预估.基于3种研究途径:理论建模与数值模拟技术、地面试验模拟技术、飞行试验技术的研究,以及它们的发展及互相验证, 用来预估飞行器的性能.为了保证地面模拟与真实飞行之间存在相似关系, 研究相似准则的作用,及分析目前的模拟能力, 涉及到许多空气动力学界至今尚未解决的难题,为了解决这些困难对今后的研究及应用提出了多方面的需求.      

Simulation of dynamic fluid–solid interactions with an improved direct-forcing immersed boundary method

Dynamic fluid–solid interactions are widely found in chemical engineering, such as in particle-laden flows, which usually contain complex moving boundaries. The immersed boundary method (IBM) is a convenient approach to handle fluid–solid interactions with complex geometries. In this work, Uhlmann's direct-forcing IBM is improved and implemented on a supercomputer with CPU–GPU hybrid architecture. The direct-forcing IBM is modified as follows: the Poisson's equation for pressure is solved before evaluation of the body force, and the force is only distributed to the Cartesian grids inside the immersed boundary. A multidirect forcing scheme is used to evaluate the body force. These modifications result in a divergence-free flow field in the fluid domain and the no-slip boundary condition at the immersed boundary simultaneously. This method is implemented in an explicit finite-difference fractional-step scheme, and validated by 2D simulations of lid-driven cavity flow, Couette flow between two concentric cylinders and flow over a circular cylinder. Finally, the method is used to simulate the sedimentation of two circular particles in a channel. The results agree very well with previous experimental and numerical data, and are more accurate than the conventional direct-forcing method, especially in the vicinity of a moving boundary.     

跨音速大迎角绕流的Euler方程数值模拟

采用保角变换与代数方法相结合，生成全场统一的贴体、正交Ｏ－Ｈ型网格．采用有限体积法求解Ｅｕｌｅｒ方程，模拟具有歼击机外形的全机及翼身组合体大迎角跨音速绕流．计算表明，法向力系数、气动中心位置及压力分布的计算结果与实验结果吻合良好     

Application of Differing Forcing Function Models on Simulated Flow Past an Oscillating Cylinder in a Uniform Low Reynolds Number Flow

A Computational Fluid Dynamics (CFD) model is presented for the uniform viscous two dimensional flow past an oscillating cylinder at low Reynolds number. Numerical simulations are made to study the effect of differing forced induced oscillation mechanisms with a large range of cylinder forcing frequencies. In the first case sinusoidal velocity slip boundary conditions are adopted for the cylinder surface to simulate cylinder oscillation. The implication suggests that no modification or additional term need to be added to the Navier-Stokes equations. In the second case this time extra body force terms which are assumed to account for velocity effects due to cylinder movement are included in the Navier-Stokes equations with the imposition of same boundary conditions. Drag and lift coefficients are extracted from present numerical results and other detailed computations of these coefficients are made at a Reynolds number of 80 and an amplitude-to diameter ratio 0.14. The results are found to be in agreement with each other at low force driving frequencies below and near lock-in. However, differences are found at higher frequencies above lock-in. Agreement are also found with experimental results at some frequency ranges.     

Fish locomotion: kinematics and hydrodynamics of flexible foil-like fins

The fins of fishes are remarkable propulsive devices that appear at the origin of fishes about 500 million years ago and have been a key feature of fish evolutionary diversification. Most fish species possess both median (midline) dorsal, anal, and caudal fins as well as paired pectoral and pelvic fins. Fish fins are supported by jointed skeletal elements, fin rays, that in turn support a thin collagenous membrane. Muscles at the base of the fin attach to and actuate each fin ray, and fish fins thus generate their own hydrodynamic wake during locomotion, in addition to fluid motion induced by undulation of the body. In bony fishes, the jointed fin rays can be actively deformed and the fin surface can thus actively resist hydrodynamic loading. Fish fins are highly flexible, exhibit considerable deformation during locomotion, and can interact hydrodynamically during both propulsion and maneuvering. For example, the dorsal and anal fins shed a vortex wake that greatly modifies the flow environment experienced by the tail fin. New experimental kinematic and hydrodynamic data are presented for pectoral fin function in bluegill sunfish. The highly flexible sunfish pectoral fin moves in a complex manner with two leading edges, a spanwise wave of bending, and substantial changes in area through the fin beat cycle. Data from scanning particle image velocimetry (PIV) and time-resolved stereo PIV show that the pectoral fin generates thrust throughout the fin beat cycle, and that there is no time of net drag. Continuous thrust production is due to fin flexibility which enables some part of the fin to generate thrust at all times and to smooth out oscillations that might arise at the transition from outstroke to instroke during the movement cycle. Computational fluid dynamic analyses of sunfish pectoral fin function corroborate this conclusion. Future research on fish fin function will benefit considerably from close integration with studies of robotic model fins.     

Forces on oscillating bodies in viscous fluid

This paper describes a method for determining the fluid forces on oscillating bodies in viscous fluid when the corresponding flow problem has been solved using the finite element method. These forces are characterized by the concept of added mass, added damping and added force. Numerical results are obtained for several example body shapes. Comparison is made with exact analytical results and other finite element results for the limiting cases of Stoke's flow and inviscid flow, and good agreement is obtained. The results for finite values of the body amplitude parameter β show the appearance of added force from the steady streaming component of the flow for asymmetric bodies. Results are also obtained for the associated flow where the fluid remote from a fixed body is oscillating.     

Influence of channel walls on the formation of vortices behind a wedge

A study is made of the influence of channel walls on the occurrence of periodic vortex formation in the wake of a wedge in a turbulent flow of incompressible fluid. The results of numerical simulation of the flow are confirmed by experimental data. The influence of the blocking of the flow by the body on the existence of the regime of vortex formation is explained by means of a linear analysis of the stability of the non-parallel flow in the near wake.     

Experimental investigation of side force control on cone-cylinder slender bodies with flexible micro balloon actuators

Side forces on slender bodies of revolution at medium to high angles of attack (AOA > 30°) has been known from a large number of investigations. Asymmetric vortex pairs over a slender body are believed to be the principle cause of the side forces. Under some flight conditions, this side force may be as large as the normal force acting on the slender body. In this paper, experimental results are presented for side force control on a cone-cylinder slender body by using microfabricated balloon actuators. The micro balloon actuators are made of polydimethylsiloxane (PDMS) elastomer by using micro molding techniques. They can be packaged on curve surfaces of a cone-cylinder slender body. As actuator is actuated, the micro balloon actuator inflates about 1.2 mm vertically, which is about 2.4% of the cylinder diameter D (=50 mm) of the cone-cylinder slender body. Micro balloon actuators are actuated at different roll angles of a cone-cylinder slender body. Aerodynamic force measurement results indicate the effects of micro balloon actuators vary at different actuation locations on the cone-cylinder slender body. The side forces can be significantly reduced if the actuators are actuated in the weak vortex side (the side corresponding to the asymmetric vortex which is far from the surface) and actuation angles are located at about 50–60° (the actuation angle here is measured from stagnation line of the incidence plane toward weak vortex side direction). Significant changes are noticed from the surface pressure, as well as leeside vortex flow field, measurement. Micro balloon actuators change nose shapes of the slender body which decide adverse-pressure-gradient values and directly influence the origin of the separation lines and characteristics of the separated vortices over the leeside surface.     

平流层飞艇空气动力估算

本文采用有限基本解方法与工程估算方法相结合的气动力工程计算方法,用以计算平流层飞艇的气动力。将飞艇所受的气动力分成飞艇艇身和尾翼所受气动力两部分,每一部分的气动力按照无粘性流产生的线性气动力和粘性引起的非线性气动力分别进行计算。根据势流理论对飞艇艇身线性气动力进行分析计算,由于飞艇艇体是旋成体,故根据Allen的横流阻力理论对其所受的非线性气动力进行计算;尾翼的线性气动力采用有限基本解方法进行计算,非线性气动力用Polhamus-Lamar吸力比拟方法估算。该方法中考虑了由于尾翼安装在体上后,处于艇体产生的上洗流场中,尾翼气动力的变化和尾翼对艇身气动力的干扰作用。通过算例的计算与实验结果比较得出该方法可以快速、准确的计算飞艇所受的气动力。     

Fins and turbulence promoters for heat transfer enhancement in latent heat storage systems

One of the serious problems associated with the operation of PCM storage system is the heat transfer in and out of the element containing the PCM. This paper presents the results of an experimental investigation of the effects of radial fins and turbulence promoters on the enhancement of phase change heat transfer external to a horizontal tube submersed in the PCM with the working fluid flowing through it. The experimental measurements were realized on a bare cupper tube and an identical cupper tube fitted with radial fins. The fins investigated are 40, 60, 120 and 180 mm diameters. A turbulence promoter made of stainless steel wire of 1.0 mm diameter coiled in a helical form with a pitch of 25.0 mm was inserted into the cupper tubes. The tests were realized on bare tubes, finned tubes and finned tubes with the turbulence promoter inserted into the finned tubes. The measurements were realized for the working fluid temperatures in the range of −10 °C, to −25 °C and six values of the mass flow rate ranging from 0.013 to 0.031 kg/s. The position of the phase interface was photographed by a high resolution digital camera and scanned to determine the real interface position by comparison with a precision measuring scale. The results of the phase interface position, velocity of the interface, solidified mass fraction and the time for complete solidification are presented in function of the working fluid temperature, the working fluid mass and the tube arrangements. The results are presented and discussed.     

Aerodynamic force measurement using 3-component accelerometer force balance system in a hypersonic shock tunnel

A new three-component accelerometer force balance has been designed, calibrated and tested in hypersonic shock tunnel (HST2) of Indian Institute of Science. The newly designed balance is able to measure aerodynamic forces (within test time of one millisecond) on test models at angles of attack from 0 to 12°. Two models, a blunt cone with after body and a blunt cone with after body and frustum are used to establish the accuracy of the force balance. The tests were conducted for the above two configurations with a constant Mach number of 8 and total enthalpy of 2.0 MJ/kg. The effectiveness of the balance is demonstrated by comparing the forces and moments of measured data with AGARD models. The flow fields around the test model are simulated using a 3D axisymmetric Navier–Stokes solver and the simulated results were compared with the measured values. Measured and computed force data are matched within ±10% for two different models tested here. The accuracy of the force balance is also estimated with the Newtonian theory and the values are approximately ±10% for the axial component and ±8% for the normal and pitching moment components.      

Heat transfer characteristics from an array of thin strips pin fins due to their exposures to a single downward jet impingement

This paper investigates the heat transfer characteristics from thin strips pin fins due to their exposure to a single circular downward air jet impingement. Five aluminum specimens were considered; each one has a rectangular base of 84 mm × 78 mm and it has an array of about 300 thin strips pin fins. A test rig consists mainly of air compressor; nozzle and protractor mechanism was setup. Experiments were done to find out the effects of attack angle, Reynolds number, nozzle-to-target spacing, lateral pitch and parallel pitch among the fins on the heat transfer characteristics. Empirical correlations were deduced to describe the experimental data. A CFD-numerical model was introduced to monitor the flow characteristics on a scale of more details than that possible in the experimental work. The comparison among the results of the present work and those by the literature shows about 50% improvement in heat transfer characteristics rather than the single jet impingement onto flat plates, cylindrical surfaces, ribbed walls and multiple jets impingement onto flat plates.     

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

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

Measurements in the flow around a marine propeller at the stern of an axisymmetric body

An experimental study of the flow around and behind an axisymmetric body driven by a marine propeller is reported. Experiments were performed in a wind tunnel to document this complex, unsteady, three-dimensional, turbulent shear flow. Measurements were made in the boundary layer and wake of the bare body with a fixed dummy hub for the propeller, with the dummy hub rotating, and finally, with the propeller in operation. A five-hole yaw probe was employed for the mean-flow measurements, and two- and threesensor hotwires were used to obtained the mean and turbulent velocity fields. Part 1 of this two-part paper describes the experimental arrangement and circumferentially-averaged results which clarify certain overall aspects of the flow when it is viewed as a rotationally-symmetric flow. These are of special interest in marine hydrodynamics. In Part 2, the triple-sensor hotwire data are analyzed using phase-averaging techniques to reconstruct the instantaneous velocity and Reynolds-stress fields downstream of the propeller to show the evolution of the wakes of individual blades, blade-tip vortices, and the complex flow associated with vortices generated at hub-blade junctions.     

Shear rate dependence of the normal force of a magnetorheological suspension

An experimental study was made of the normal force in a magnetorheological suspension under shear flow with the magnetic field applied normally to the shear planes. Under no deformation, the normal force was found to increase with magnetic flux density, and acted to push apart the rheometer plates. At the start-up of shearing under a constant shear rate, the normal force decreased with strain and reached a plateau value which became smaller as shear rate was raised. This behaviour can be interpreted in terms of the aggregate model of magnetorheological suspensions: the normal force arises from the elongated aggregates spanning the gap between the plates, and these are broken up into smaller particle clusters under shearing.