Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been developed to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries

Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by \delt T_{Lp} and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of \delt T_{Lp} on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading

中国科学院力学研究所非线性力学国家重点实验室;

湍流的小尺度混合特性对湍流燃烧过程至关重要。在大涡模拟中,由于采用空间滤波操作,湍流的亚格子尺度运动被过滤掉,亚格子尺度的湍流混合特性需要进一步模化。假设的标量在亚格子尺度中的概率密度分布函数是一种主要的模化方法。作者之前的研究显示,使用标量的一阶和二阶滤波矩作为参数的假设双高斯滤波密度函数模型可以很好的描述标量在亚格子中的混合状态,用于预测三阶和四阶滤波矩时也得到了很好的结果。这预示着,标量的滤波密度函数可能是一阶和二阶标量滤波矩的某种函数。因此,以一阶和二阶标量滤波矩作为条件变量,从解析的角度推导出湍流亚格子混合的条件滤波密度函数表达式,并与湍流标量混合的DNS结果互相验证。更进一步地,将..

Self-propulsion of flapping bodies in viscous fluids: Recent advances and perspectives

Flapping-powered propulsion is used by many animals to locomote through air or water. Here we review recent experimental and numerical studies on self-propelled mechanical systems powered by a flapping motion. These studies improve our understanding of the mutual interaction between actively flapping bodies and surrounding fluids. The results obtained in these works provide not only new insights into biolocomotion but also useful information for the biomimeticdesign of artificial flyers and swimmers.</div

How flexibility affects the wake symmetry properties of a self-propelled plunging foil

The wake symmetry properties of a flapping-foil system are closely associated with its propulsive performance. In the present work, the effect of the foil flexibility on the wake symmetry properties of a self-propelled plunging foil is studied numerically. We compare the wakes of a flexible foil and a rigid foil at a low flapping Reynolds number of 200. The two foils are of the same dimensions, flapping frequency, leading-edge amplitude and cruising velocity but different bending rigidities. The results indicate that flexibility can either inhibit or trigger the symmetry breaking of the wake. We find that there exists a threshold value of vortex circulation above which symmetry breaking occurs. The modification of vortex circulation is found to be the pivotal factor in the influence of the foil flexibility on the wake symmetry properties. An increase in flexibility can result in a reduction in the vorticity production at the leading edge because ofthe decrease in the effective angle of attack, but it also enhances vorticity production at the trailing edge because of the increase in the trailing-edge flapping velocity. The competition between these two opposing effects eventually determines the strength of vortex circulation, which, in turn, governs the wake symmetry properties. Further investigation indicates that the former effect is related to the streamlined shape of the deformed foil while the latter effect is associated with structural resonance. The results of this work provide new insights into the functional role of passive flexibility inflapping-based biolocomotion

On applicability of the Kutta-Joukowski theorem to low-Reynolds-number unsteady flows

The limitations of the Kutta-Joukowski (K-J) theorem in prediction of the time-averaged and instantaneous lift of an airfoil and a wing in low-Reynolds-number unsteady flows are examined. A general lift formula for a rectangular control volume is given in a very simple form in the framework of viscous flow theory, which provides a rational foundation for a direct comparison with the K-J theorem considered as a reduced case. Direct numerical simulations on the stationary and flapping flat plate and rectangular wing are conducted to assess the accuracy of both the K-J theorem and the general lift formula. In particular, the Lamb vector integral for the vortex force and the acceleration term of fluid for the unsteady inertial effect are evaluated as the main contributions to the unsteady lift generation of a flapping wing

近海生态系统碳汇过程、调控机制及增汇模式

海洋是地球上最大的碳库,发挥着全球气候变化缓冲器的作用.蓝色碳汇,简称蓝碳,即由海洋生态系统捕获的碳(主要是有机碳),是海洋储碳的重要机制之一.; 蓝碳最初认识的形式是可见的海岸带植物固碳.其实之前没有得到足够重视的、看不见的微型生物(浮游植物、细菌、古菌、病毒、原生动物)占海洋生物量90%; 以上,是蓝碳的主要贡献者.中国陆架边缘海占国土总面积的1/3,碳汇潜力巨大,亟待研发.本文以近海生态系统碳汇过程、调控机制及增汇模式为主线,论述; 了近海生态系统结构与碳循环功能特征、碳汇形成过程与机理,并结合近海碳汇在沉积记录中的地史过程演变探讨了自然过程和人类活动对碳汇的可能影响,展望了; 碳汇工程在增加近海海洋储碳能力方面的应用前景.国家重大科学研究计划项目; 国家重点研发计划项目; 国家自然科学基金项目; 国家海洋局全球变化与海气相互作用专项项

Mapping Closure Approximation for Reactive Scalars in Random Flows

The Mapping Closure Approximation (MCA) approach is developed to describe the statistics of both conserved and reactive scalars in random flows. The statistics include Probability Density Function (PDF), Conditional Dissipation Rate (CDR) and Conditional Laplacian (CL). The statistical quantities are calculated using the MCA and compared with the results of the Direct Numerical Simulation (DNS). The results obtained from the MCA are in agreement with those from the DNS. It is shown that the MCA approach can predict the statistics of reactive scalars in random flows

Self-propelled swimming of a flexible plunging foil near a solid wall

Numerical simulations are conducted to investigate the influences of a solid wall on the self-propelled swimming of a flexible plunging foil. It is found that the presence of a solid wall enhances the cruising speed, with the cost of increasing input power. Rigid foil can achieve high percentage increase in cruising speed when swimming near a solid wall, but the propulsive efficiency may be reduced. Foils with some flexibility can enjoy the enhancements in both cruising speed and propulsive efficiency. Another advantage of the flexible foils in near-wall swimming is that smaller averaged lateral forces are produced. The effects of wall confinement on the wake structure and the vortex dynamics are also studied in this paper. The results obtained in this study shed some light on the unsteady wall effect experienced by aquatic animals and also inform the design of bio-mimetic underwater vehicles which are capable of exploiting the wall effect

1978～2008年中国湿地类型变化

分别基于美国陆地卫星(Landsat MSS/TM/ETM+)和中巴资源卫星(CBERS-02B)影像数据,以人工目视解译为主,完成了中国1978～2008年4期(基准年分别为1978,1990,2000和2008年)湿地遥感制图,并进行了大量的室内外验证.在此基础上,对我国湿地现状及近30年来湿地变化进行了初步分析,得到以下主要结论:(ⅰ)截止2008年,中国湿地面积约为324097km2,其中以内陆沼泽(35%)和湖泊湿地(26%)为主.(ⅱ)1978～2008年,中国湿地面积减少了约33%,而人工湿地增加了约122%.过去30年里湿地减少的速度大幅降低,由最初5523km2/a(1978～1990年)降为831km2/a(2000～2008年).(ⅲ)减少的自然湿地(包括滨海湿地和内陆湿地),其类型变化由湿地向非湿地转化的比例逐渐降低.初期(1978～1990年)几乎全部(98%)转换为非湿地;在1990～2000年间减少的自然湿地约有86%转化为非湿地,而在2000～2008年,这一比例下降为77%.(ⅳ)气候变化和农业活动是中国湿地变化的主要驱动因素,湿地变化在中国分为三大不同特征区域,即西部三省/自治区(西藏、新疆和青海)、北部两省/自治区(黑龙江和内蒙古)和其他省市区.其中西部区域尤其是青藏高原,湿地变化的驱动因子以气候增温为主;新疆湿地由于气候增温和农业活动共同作用造成变化不大.北部省/自治区的湿地变化则主要由农业活动引起;而其他省市区的湿地变化几乎完全受控于人类的农业经济活动