Experimental benchmark of a free plunging wing with imposed flap oscillations

The validation of fluid–structure interaction solvers is difficult since there is a lack of experimental data. Therefore, in this work an aeroelastic experiment is presented. The focus is on the temporal coupling between fluid and structure dynamics. Issues in the spatial coupling are eliminated by using a rigid wing. The wing, with a harmonically actuated 0.2c trailing edge flap, has a degree of freedom in the plunge (vertical) direction. The wing has a chord of 0.5 m and is suspended with springs. The wing motion is constrained by a vertical rail system.For simplicity attached flow is desired and therefore the set angle of attack is α=0°. The Reynolds number is approximately Re=700 000 and the flap deflects over a range of about ±2°. The damped natural frequency of the structure expressed as a reduced frequency is about k=0.194 and measurements are performed for reduced flap frequencies ranging from k=0.1 to k=0.3. Displacements and time dependent aerodynamic forces are measured and for k=0.198 2-D PIV measurements are performed. The planar PIV measurements are used to intrinsically determine the unsteady loads using Noca׳s method.As expected the aeroelastic problem shows similarities with a viscously damped mass–damper–spring, meaning the maximum excursion of the wing is found near the system eigenfrequency. The lift is dominated by the flap motion and the effective angle of attack due to the motion introduces phase shifts of the lift signal with respect to the flap phase angle.The experiment has been set up and executed with the necessary precision, but small ambiguities are found in the lift and drag disqualifying the data for validation. Nevertheless the data set provides a clear insight into typical loads and motions and can be used for comparative studies. It can also be used to (re)design future experiments to improve the quality of the data to the desired level of accuracy for validation.     

A novel hybrid solid-like fluid-like (SLFL) method for the simulation of dry granular flows

This paper proposes a novel hybrid method to simulate the dry granular flow of materials over a wide range of inertial numbers that simultaneously covers the quasi-static and dense granular flow regimes. To overcome the lack of incremental objectivity whenever large deformations occur in solid-like regimes and to remove computational singularities in fluid-like regimes close to rest, the elastic–perfectly plastic theory based on the Drucker–Prager yield criterion is combined with the theory of dense granular flows. By implementing some new modifications at the boundaries and removing all ghost particles, smoothed particle hydrodynamics (SPH) is used as the framework for the method. A number of benchmark problems have been solved to show the capabilities of the new modified SPH method. Precise prediction of both location and pressure makes the modifications comparable with the previous works on SPH. Finally, the method is used to solve the classic 2D dry granular cliff collapse problem and to model dry granular material flow inside a rotary drum. The outcomes of the numerical simulation show good agreement with tabletop experiments and published results.     

Stress induced phase transitions in silicon

Silicon has a tremendous importance as an electronic, structural and optical material. Modeling the interaction of a silicon surface with a pointed asperity at room temperature is a major step towards the understanding of various phenomena related to brittle as well as ductile regime machining of this semiconductor. If subjected to pressure or contact loading, silicon undergoes a series of stress-driven phase transitions accompanied by large volume changes. In order to understand the material's response for complex non-hydrostatic loading situations, dedicated constitutive models are required. While a significant body of literature exists for the dislocation dominated high-temperature deformation regime, the constitutive laws used for the technologically relevant rapid low-temperature loading have severe limitations, as they do not account for the relevant phase transitions. We developed a novel finite deformation constitutive model set within the framework of thermodynamics with internal variables that captures the stress induced semiconductor-to-metal ($\mathrm{cd-Si}\to \mathrm{\beta }\mathrm{-Si}$), metal-to-amorphous ($\mathrm{\beta }\mathrm{-Si}\to \mathrm{a-Si}$) as well as amorphous-to-amorphous ($\mathrm{a-Si}\to \mathrm{hda-Si}$, $\mathrm{hda-Si}\to \mathrm{a-Si}$) transitions. The model parameters were identified in part directly from diamond anvil cell data and in part from instrumented indentation by the solution of an inverse problem. The constitutive model was verified by successfully predicting the transformation stress under uniaxial compression and load–displacement curves for different indenters for single loading–unloading cycles as well as repeated indentation. To the authors' knowledge this is the first constitutive model that is able to adequately describe cyclic indentation in silicon.     

Retrieval of aerosol size distribution using improved quantum-behaved particle swarm optimization on spectral extinction measurements

An improved quantum-behaved particle swarm optimization (IQPSO) algorithm is employed to determine aerosol size distribution (ASD). The direct problem is solved using the anomalous diffraction approximation and Lambert–Beer's Law. Compared with the standard particle swarm optimization algorithm, the stochastic particle size optimization algorithm and the original QPSO, our IQPSO has faster convergence speed and higher accuracy within a smaller number of generations. Optimization parameters for the IQPSO were also evaluated; we recommend using four measurement wavelengths and 50 particles. Size distributions of various aerosol types were estimated using the IQPSO under dependent and independent models. Finally, experimental ASDs at different locations in Harbin were recovered using the IQPSO. All our results confirm that the IQPSO algorithm is an effective and reliable technique for estimating ASD.     

Estimating the tensile strength of super hard brittle materials using truncated spheroidal specimens

New approaches need to be introduced to measure the tensile capacity of super hard materials since the standard methods are not effective. To pursue this objective, a series of laboratory tests were constructed to replicate the fracture mechanism of diamond-based materials. Experiments indicate that under a certain compressive test condition, stresses normal to the axisymmetric line in truncated spheroidal specimens (bullet-shaped specimens) are in tension contributing to the tensile fracture of the material. From experimental and numerical studies, it is concluded that semi-prolate spheroidal specimens can be used to determine precisely the tensile strength of brittle stiff diamond-like composites.     

On the spray drying of uniform functional microparticles

Spray drying is a typical method to produce particles in dry powder forms at industrial scale. Most spray-dried products often show a wide range of particle properties even within the same batch. At Monash University, we utilise a microfluidic spray drying approach to generate uniform microparticles with tightly controlled characteristics and sizes in a scalable, almost waste-free process. The technique is useful to correlate the effects of formulation and spray drying conditions on the properties of spray-dried particles, and can be used to test new formulations for targeted applications such as encapsulation and release of active ingredients. The synthesis route can be applied to other self-assembling systems, including mesoporous, crystalline, and hierarchically structured microparticles. As spray drying is commonly used in commercial scales, the understanding of how functional particles are formed in relation to formulations and process conditions could assist in developing a cost effective, energy and material-efficient route to produce powders with better properties and ease of handling for more advanced applications such as selective adsorption and bio-separation.     

Relating geologic units and mobility system kinematics contributing to Curiosity wheel damage at Gale Crater,Mars

Curiosity landed on plains to the north of Mount Sharp in August 2012. By June 2016 the rover had traversed 12.9 km to the southwest, encountering extensive strata that were deposited in a fluvial-deltaic-lacustrine system. Initial drives across sharp sandstone outcrops initiated an unacceptably high rate of punctures and cracks in the thin aluminum wheel skin structures. Initial damage was found to be related to the drive control mode of the six wheel drive actuators and the kinematics of the rocker-bogie suspension. Wheels leading a suspension pivot were forced onto sharp, immobile surfaces by the other wheels as they maintained their commanded angular velocities. Wheel damage mechanisms such as geometry-induced stress concentration cracking and low-cycle fatigue were then exacerbated. A geomorphic map was generated to assist in planning traverses that would minimize further wheel damage. A steady increase in punctures and cracks between landing and June 2016 was due in part because of drives across the sharp sandstone outcrops that could not be avoided. Wheel lifetime estimates show that with careful path planning the wheels will be operational for an additional ten kilometers or more, allowing the rover to reach key strata exposed on the slopes of Mount Sharp.     

Effect of adhesive thickness,adhesive type and scarf angle on the mechanical properties of scarf adhesive joints

The effects of adhesive thickness, adhesive type and scarf angle, which are determined as the main control parameters by the dimensional analysis, on the mechanical properties of a scarf adhesive joint (SJ) subjected to uniaxial tensile loading are examined using a mixed-mode cohesive zone model (CZM) with a bilinear shape to govern the interface separation. Particularly, the adhesive-dependence of the vital cohesive parameters of CZM, which mainly include initial stiffness, total fracture energy and separation strength, is introduced emphatically. The numerical results demonstrate that the ultimate tensile loading increases as the adhesive thickness decreases. Cross the ultimate tension, the joint loses the load-bearing capacity when adopting the brittle adhesive but sustains partial load-bearing capacity while selecting the ductile adhesive. In addition, for the joint with the ductile adhesive, the maximum applied displacement until the complete failure of it is directly proportional to the adhesive thickness, which is different from the case using the brittle adhesive. Taking the combination of the ultimate loading and applied displacement into account, failure energy is employed to evaluate the joint performances. The results show that the failure energy of the joint with the brittle adhesive increases as the adhesive thickness decreases. Conversely, the situation of the joint using the ductile adhesive is vice versa. Moreover, the effect of the adhesive thickness becomes more noticeable with decreasing the scarf angle owing to the variation of the proportion of each component of the mixed-mode. Furthermore, all the characteristic parameters (the ultimate tensile loading, the maximum applied displacement and the failure energy) that adopted to describe the performances of SJ increase as the scarf angle decreases. Finally, the numerical method employed in this study is validated by comparing with existing experimental results.     

Modeling dynamic brittle behavior of materials with circular flaws or pores

Compressive failure of brittle materials is driven primarily by crack growth from pre-existing flaws in the material. These flaws, such as grain boundaries, pores, preexisting cracks, inclusions and missing grains, are randomly spaced and have a range of possible shapes and sizes. The current work proposes a micromechanics-based model for compressive dynamic failure of brittle materials with circular pore flaws, which incorporates both the number density and the size distribution of flaws. Results show that the distribution of flaw sizes is very important, particularly at moderate strain rate, since analyses based solely on the mean flaw size overpredict strength. Therefore, in order to increase dynamic strength at low to moderate strain rates, it is most effective to control the presence of large flaws. At very high strain rates, however, crack growth is activated even in small flaws and therefore controlling the total number density rather than the size of the flaws is effective for increasing dynamic strength. Finally, the model shows that neglecting very small flaws in the pore population may not have significant effects on the results in many cases, suggesting that the model is a useful tool for identifying a minimum resolution required for experimental characterization of microstructure.     

Offshore floating wind turbine and its dynamic problems

Green energy sources and ocean wind power are plentiful in deep sea. More and more offshore wind power plants are constructed in the deep water over hundred meters below the surface. While offshore floating wind turbine system is working, wind turbine, floating foundation, and mooring system affect each other with wind, waves, and currents acting on them. Various offshore floating wind turbine systems and the encountered environmental loads are briefly reviewed and discussed. It is difficult and crucial to comprehensively analyze the aerodynamic-hydrodynamic-service system-structure under the coupling effect of offshore floating wind turbine system. The environmental flow field, structure scale, and rational applications of theories and approaches should be well considered in advance.     

Aeroelastic analysis of wind turbines using a tightly coupled CFD–CSD method

This paper presents aeroelastic analyses of wind turbines, using the compressible flow Helicopter Multi-Block (HMB2) solver of Liverpool University, coupled with a Computational Structural Dynamics method. For this study, the MEXICO and NREL Phase VI wind turbines were employed. A static aeroelastic method was first employed for the analysis of the MEXICO blade and the effect of the torsional stiffness was studied at 10, 15 and 24 m/s axial wind speeds. The torsional deformations showed strong dependency on this parameter and the blade region from mid-span to the tip was the most susceptible to aeroelastic effects. The work progressed by studying both the static and dynamic response on the NREL wind turbine, where the nacelle and the tower were considered. Mean deflections between the static and dynamic methods showed consistency and, due to the structural properties of this blade, flapping modes were dominant. The dynamic aeroelastic method enabled an assessment of the effect of flapping on the blade loads, in conjunction with the effect of tower. Aeroelastic effects were found to be secondary for the MEXICO blade, but had a stronger effect on the larger NREL Phase VI blade.     

An experimental investigation on the characteristics of fluid–structure interactions of a wind turbine model sited in microburst-like winds

An experimental investigation is performed to assess the characteristics of the fluid–structure interactions and microburst-induced wind loads acting on a wind turbine model sited in microburst-liked winds. The experiment study was conducted with a scaled wind turbine model placed in microburst-like winds generated by using an impinging-jet-typed microburst simulator. In addition to quantifying complex flow features of microburst-like winds, the resultant wind loads acting on the turbine model were measured by using a high-sensitive force–moment sensor as the turbine model was mounted at different radial locations and with different orientation angles with respect to the oncoming microburst-like winds. The measurement results reveal clearly that, the microburst-induced wind loads acting on the turbine model were distinctly different from those in a conventional atmospheric boundary layer (ABL) wind. With the scales of the wind turbine model and the microburst-like wind used in the present study, the dynamic wind loadings acting on the turbine model were found to be significantly higher (i.e., up to 4 times higher for the mean loads, and up to 10 times higher for the fluctuation amplitudes) than those with the same turbine model sited in ABL winds. Both the mean values and fluctuation amplitudes of the microburst-induced wind loads were found to vary significantly with the changes of the mounted site of the turbine model, the operating status (i.e., with the turbine blades stationary or freely rotating), and the orientation angle of the turbine model with respect to the oncoming microburst-like wind. The dynamic wind load measurements were correlated to the flow characteristics of the microburst-like winds to elucidate underlying physics. The findings of the present study are helpful to gain further insight into the potential damage caused by the violent microbursts to wind turbines to ensure safer and more efficient operation of the wind turbines in thunderstorm-prone areas.     

Aeroelastic analysis of large horizontal wind turbine baldes

A nonlinear aeroelastic analysis method for large horizontal wind turbines is described. A vortex wake method and a nonlinear finite element method (FEM) are coupled in the approach. The vortex wake method is used to predict wind turbine aerodynamic loads of a wind turbine, and a three-dimensional (3D) shell model is built for the rotor. Average aerodynamic forces along the azimuth are applied to the structural model, and the nonlinear static aeroelastic behaviors are computed. The wind rotor modes are obtained at the static aeroelastic status by linearizing the coupled equations. The static aeroelastic performance and dynamic aeroelastic responses are calculated for the NH1500 wind turbine. The results show that structural geometrical nonlinearities significantly reduce displacements and vibration amplitudes of the wind turbine blades. Therefore, structural geometrical nonlinearities cannot be neglected both in the static aeroelastic analysis and dynamic aeroelastic analysis.     

复杂工况下的大型风力机气动性能和尾迹研究Investigation of aerodynamic performance and wake of the large scale wind turbine in complicated conditions

在复杂工况下,大型风力机非定常特性会更严重,导致风力机气动性能变化和尾迹预测更加复杂。本文主要针对稳态偏航、动态偏航、风剪切和随机风速场等复杂工况,基于自由涡尾迹方法,嵌入复杂工况的模块,加入了动态失速模型和三维旋转效应模型修正,实现了复杂工况数值模拟计算,比较了不同复杂工况的气动载荷和尾迹形状。最后,得出了风力机在复杂工况下的气动性能、载荷和尾迹叶尖涡线特性,并计算出风力机在复杂工况下的气动载荷超调量和迟滞时间。对推进自由涡尾迹方法应用于风力机工程的大批工况载荷计算,提高大型风力机的载荷计算精度和设计水平等具有重要意义。     

考虑流固耦合作用的充气膜结构风压分布研究

充气膜结构是典型的风敏感型柔性结构,风荷载经常起关键的控制作用。本文利用ANSYS14.5程序中的workbench平台,考虑流固耦合作用,研究矩形平面气承式充气膜结构的风压系数分布。其中,选用基于雷诺时均模拟法的RNGk-ε湍流模型进行风场模拟,采用弱耦合分析方法模拟流固耦合风荷载效应。分析的参数选择风向角、结构内压、矢跨比和平面长宽比。针对矢跨比分别为1/4,1/3和1/2,长宽比分别为5/3,2/1和3/1的柔性充气膜结构模型,计算不同内压及不同风向角作用下的结构响应。结果表明,考虑流固耦合作用时,充气膜结构的风压体型系数比不考虑流固耦合作用的刚性模型明显偏大,其影响因子在1.25~1.5之间;充气膜结构的风压系数分布受风向角、内压、长宽比及矢跨比的影响较大。     

浮式垂直轴风机的动力学建模、仿真与实验研究

考虑气动力和水动力的耦合研究浮式垂直轴风机系统的运动响应,将固定式垂直轴风机的气动载荷计算方法进一步推广到海上浮式垂直轴风机的气动载荷计算.考虑阻尼力、波浪力、风载荷、系泊力等,建立了浮式垂直轴风机系统的纵荡-垂荡-纵摇运动方程.考虑动态失速和浮式基础运动,基于双致动盘多流管理论,推导了风机叶片气动载荷计算公式,编制了数值计算程序.以Sandia 17 m风机为例,验证了气动载荷计算程序的正确性.最后进行了模型实验,其中模型的风机为Φ型达里厄垂直轴风机,支撑基础为桁架式Spar型浮式基础,将模型实验结果与数值计算结果进行了对比,验证了耦合计算程序.结果表明,数值计算得到的风机系统的垂荡、纵摇运动的RAO(幅值响应算子)曲线与模型实验结果吻合较好,验证了耦合程序的正确性.然而,由于数值计算与模型实验在运动自由度、阻尼、风载荷等方面存在差别,数值计算结果与模型实验结果仍有一定的差异.     

漂浮运动对风力机气动特性的影响分析

海上浮式风机为蕴藏丰富的深海风能开发提供了有效的解决方案,浮式基础存在的大幅度纵荡、纵摇和艏摇运动可以改变风轮与流场间的相互作用,从而影响风力机的气动特性。基于叶素-动量理论及其修正方法,以NREL 5MW风力机为研究对象,考虑浮式基础运动对叶片不同径向位置处相对入流风速的影响,提出了风载荷的计算模型,通过编程计算获得了叶轮转矩和风力机功率,并比较了不同运动形式对风力机功率波动的影响。结果表明,纵摇对其功率特性影响最大,这为海上浮式风机的优化设计提供理论依据与数据基础。     

漂浮运动对风力机气动特性的影响分析Effect of floating foundation motion on the aerodynamic characteristics of the floating offshore wind turbine

海上浮式风机为蕴藏丰富的深海风能开发提供了有效的解决方案,浮式基础存在的大幅度纵荡、纵摇和艏摇运动可以改变风轮与流场间的相互作用,从而影响风力机的气动特性。基于叶素-动量理论及其修正方法,以NREL 5MW风力机为研究对象,考虑浮式基础运动对叶片不同径向位置处相对入流风速的影响,提出了风载荷的计算模型,通过编程计算获得了叶轮转矩和风力机功率,并比较了不同运动形式对风力机功率波动的影响。结果表明,纵摇对其功率特性影响最大,这为海上浮式风机的优化设计提供理论依据与数据基础。     

Recent progresses in studies on wave-current loads on foundation structure with piles and slab

The foundation structure with piles and slab is widely used in offshore wind farm construction in shallow water. Experimental studies on the hydrodynamic loads acting on the piles and slab under irregular waves and currents are summarized with discussion on the effects of pile grouping on the wave forces and wave impact loads on the slab locating near the free surface. By applying the theoretical solution of the wave diffracted by the slab and using the Morison equation to evaluate the wave force on the piles, the effects of the slab on the wave forces acting on the piles are analyzed. Based on the Reynolds-averaged Navier-Stokes (RANS) equations and the volume of fluid (VOF) method, a numerical wave basin is developed to simulate the wave-structure interaction. The computed maximum wave force on the foundation structure with piles and slab agrees well with the measured data. The violent deformation, breaking, and run-up of the wave around the structure are presented and discussed. Further work on the turbulent flow structures and large deformation of the free surface due to interaction of the waves and foundation structures of offshore wind farms needs more efficient approaches for evaluating hydrodynamic loads under the effects of nonlinear waves and currents.     

随机风速下高速列车的运行安全可靠性

基于可靠性理论提出了随机风速作用下高速列车风致安全分析的新方法, 这种方法可以对随机风速作用下高速列车的运行安全可靠性进行有效评估.首先基于Cooper理论和谐波叠加法计算随车移动点的脉动风速, 建立随机风速作用下高速列车非定常气动载荷的计算方法, 并通过数值仿真得到气动载荷系数的标准差随侧偏角的变化规律. 然后建立高速列车车辆系统动力学模型, 并对计算模型的正确性进行验证.最后以随机风速、侧力系数、升力系数、侧滚力矩系数、摇头力矩系数和点头力矩系数为基本随机变量, 研究随机风速作用下高速列车的运行安全可靠性和可靠性灵敏度, 给出随机风速作用下高速列车的概率特征风速曲线.研究结果表明:随着车速和风速的增大, 系统的失效概率增大;通过可靠性灵敏度分析发现侧力系数和侧滚力矩系数对高速列车的运行安全影响最大, 应该特别注意这两个参数的变化对高速列车运行安全性的影响;传统确定性方法得到的高速列车的安全域曲线偏于保守, 基于可靠性的方法可得到更合理的安全域曲线.