Finite element solutions of axisymmetric Euler equations for an incompressible and inviscid fluid

In this paper we present a finite element method for the numerical solution of axisymmetric flows. The governing equations of the flow are the axisymmetric Euler equations. We use a streamfunction angular velocity and vorticity formulation of these equations, and we consider the non-stationary and the stationary problems. For industrial applications we have developed a general model which computes the flow past an annular aerofoil and a duct propeller. It is able to take into account jumps of angular velocity and vorticiy in order to model the flow in the presence of a propeller. Moreover, we compute the complete flow around the after-body of a ship and the interaction between a ducted propeller and the stern. In the stationary case we have developed a simple and efficient version of the characteristics/finite element method. Numerical tests have shown that this last method leads to a very fast solver for the Euler equations. The numerical results are in good agreement with experimental data.     

Three Dimensional Numerical Flow Simulation Around Parallel Rectangular Cylinders

Noises of flow around parallel rectangular cylinders are likely to be caused by wind flow. According to the wind tunnel experiment, it is known that there are three kinds of special noises. Two kinds of noises occur when the wind comes from a perpendicular direction to the cylinders. These noises are caused by the vortex excited oscillation. Another noise occurs when the wind comes from an inclined direction to the cylinders. This noise is very high frequency, and its mechanism is not clear. Therefore in this study, the three dimensional numerical flow simulation is performed to clarify the high frequency noise of flow around parallel rectangular cylinders.     

Wake impacts on aerodynamic and aeroelastic behaviors of a horizontal axis wind turbine blade for sheared and turbulent flow conditions

The magnitude and temporal variations of wind speed considerably influence aerodynamic and structural responses of MW-sized horizontal axis wind turbines. Thus, this paper investigates the variations in airloads and blade behavior of a wind turbine blade resulting from operations in sheared and turbulent flow conditions. First, in order to validate the present methods, comparisons of aerodynamic results were made among the blade element momentum method, free-wake method, and numerical results from the previous studies. Then, the validated methods were applied to a national renewable energy laboratory 5 MW reference wind turbine model for fluid–structure interaction analyses. From the numerical simulations, it can be clearly seen that unfavorable airloads and blade deformations occur due to the sheared and turbulent flow conditions. In addition, it is clear that wake impacts are not as substantial at those of high wind speeds; however, the effects obviously affect the aerodynamic and structural behaviors of the blade at lower wind speeds. Therefore, it is concluded that the numerical results markedly indicate the demand for accurate assessment of wake dynamics for accurate estimations of the aerodynamic and structural responses for sheared and turbulent flow environments.     

某调距导管桨系泊工况水力性能和流激直发声预报

采用经验证的计算流体力学方法,对某艏辅推调距导管桨设计螺距和系泊工况螺距的水动力性能进行了有效预报,并对系泊工况装船桨流激噪声进行了分析。系泊工况下,由于导管桨的抽吸作用在导管外壁近壁面区域存在与导管内部流动方向相反的逆向流动,且导管桨尾流场速度梯度分布不均匀、流动紊乱,此时桨叶与导管的推力之比约为1.2∶1。系泊工况船+桨的瞬态流场脉动信息表明,导管桨各部件噪声源强度均表现出从1倍到4倍叶频依次下降的规律,最强幅值集中在桨叶导边和导管内壁;在远场声源级频谱曲线中轴向测点线谱较高峰值位置体现出导管桨进流流场的流动特性。对比分析该艏辅推整体和各部件宽带声源级指向性,可知旋转部件(桨叶、桨榖)对总噪声级的贡献较大,静止部件是径向测点噪声的主要贡献源。     

SWT-120风洞稳定段的性能测量

在普通超音速风洞中,由于受到噪声干扰,很难进行有效的边界层转捩特性试验研究,针对超声速流动特点发展较低噪声风洞十分必要,而稳定段设计的好坏直接影响到下游试验段噪声水平。本文介绍一座低噪声风洞稳定段的结构和性能测量结果。先对稳定段的结构设计做了简单介绍,然后对实验结果进行分析,实验结果表明在大角度扩散段内装置孔锥,稳定段安装消音夹层和阻尼网组等部件后,气流的速度脉动和压力脉动明显降低,其中压力脉动降低一个量级,速度脉动为1%。进一步优化设计和改进工艺,速度脉动还可进一步降低。测量结果表明SWT-120稳定段的设计是成功的,对我们以后发展更高性能的静风洞有借鉴和参考价值。     

Wave envelope and infinite elements for acoustical radiation

Finite element models are presented for the calculation of near and far field acoustical radiation. These models are applied to the specific problem of fan noise radiation from axisymmetric turbofan inlets. In all cases conventional acoustic finite elements are used within an inner region close to the inlet. The far field is represented by infinite elements or wave envelope elements. Theory and results are presented for the case with zero mean flow. Comparisons of computed data with analytic solutions and measured values establish the utility of both the infinite element and wave envelope element schemes in determining the near field values of acoustical pressure. The wave envelope scheme is shown to be effective also in the far field. Both schemes use meshes an order of magnitude more sparse that would be required in conventional numerical discretizations, and may consequently be applied at modest computational cost.     

Third-Order Upwind Finite Element Simulation of Flow Around Two Square Cylinders

The aerodynamic behavior of the flow around two square cylinders is presented on the basis of the numerical simulation of the incompressible Navier-Stokes equations using a third-order upwind finite element scheme. It is well known that flow patterns around the two square cylinders are more complicated than flow patterns around one square cylinder because of interference between the Karman vortices behind the two square cylinders. In this paper, two kinds of cylinder arrangements are chosen as computational models. One type is that of two square cylinders arranged vertically to the direction of a uniform flow, and the other is arranged horizontally to the direction of a uniform flow.     

A finite element method for the three-dimensional non-steady Navier-Stokes equations

The algorithm for solving the three-dimensional non-steady Navier-Stokes equations by the explicit forward Euler method is shown and the Galerkin finite element formulation is presented. As a numerical example, an entrace flow in a square duct is illustrated.     

Production of oscillatory flow in wind tunnels

A method for producing oscillatory flow in open-circuit wind tunnels driven by centrifugal fans is described. Performance characteristics of a new device installed on two such tunnels of greatly differing size are presented. It is shown that sinusoidal variations of the working section flow, having peak-to-peak amplitudes up to at least 30% of the mean flow speed and frequencies up to, typically, that corresponding to the acoustic quarter-wave-length frequency determined by the tunnel size, can be obtained with negligible harmonic distortion or acoustic noise difficulties. A brief review of the various methods that have been used previously is included and the advantages and disadvantages of these different techniques are highlighted. The present technique seems to represent a significant improvement over many of them.     

Residual‐based stabilization of the finite element approximation to the acoustic perturbation equations for low Mach number aeroacoustics

The acoustic perturbation equations (APE) are suitable to predict aerodynamic noise in the presence of a non‐uniform mean flow. As for any hybrid computational aeroacoustics approach, a first computational fluid dynamics simulation is carried out from which the mean flow characteristics and acoustic sources are obtained. In a second step, the APE are solved to get the acoustic pressure and particle velocity fields. However, resorting to the finite element method (FEM) for that purpose is not straightforward. Whereas mixed finite elements satisfying an appropriate inf–sup compatibility condition can be built in the case of no mean flow, that is, for the standard wave equation in mixed form, these are difficult to implement and their good performance is yet to be checked for more complex wave operators. As a consequence, strong simplifying assumptions are usually considered when solving the APE with FEM. It is possible to avoid them by resorting to stabilized formulations. In this work, a residual‐based stabilized FEM is presented for the APE at low Mach numbers, which allows one to deal with the APE convective and reaction terms in its full extent. The key of the approach resides in the design of the matrix of stabilization parameters. The performance of the formulation and the contributions of the different terms in the equations are tested for an acoustic pulse propagating in sheared‐solenoidal mean flow, and for the aeolian tone generated by flow past a two‐dimensional cylinder. Copyright © 2016 John Wiley & Sons, Ltd.     

Prolonged simulation of near-free surface underwater explosion based on Eulerian finite element method

In the area of naval architecture and ocean engineering, the research about the underwater explosion problem is of great significance. To achieve prolonged simulation of near-free surface underwater explosion, the underwater explosion transient numerical model is established in this paper based on compressible Eulerian finite element method(EFEM). Compared with Geers–Hunter formula, EFEM is availably validated by simulating the free-field underwater explosion case. Then, the bubble pulsation and flow field dynamic characteristics of the cases with different underwater explosive depth are compared in this work. Lastly, the height of the water hump and the pressure of flow flied are analyzed quantitatively through the simulation results.     

基于非结构网格的TTGC有限元格式的实现及在超声速流动中的应用

基于非结构网格系统,实现了时空三阶精度的TTGC有限元格式,并在三阶TTGC格式上发展了基于人工粘性的激波捕捉技术。在非结构网格下,采用这种方法对若干典型的超声速流动问题(SOD激波管、马赫数为3的前台阶流动以及马赫数为8的高超声速圆柱流动)进行了验证计算。结果表明,TTGC格式分辨率高,在粗糙网格下能够准确的模拟超声速流场中的激波、接触间断等复杂流动现象,并且能有效的控制间断附近的数值色散现象。与传统的有限体积方法相比,本文实现的TTGC有限元格式在模拟超声速流动问题方面具有格式精度高、数值耗散小等优点。     

The non‐reflective interface: an innovative forcing technique for computational acoustic hybrid methods

The present article concerns a commonly used methodology for the numerical simulation of acoustic emission and propagation phenomena. We consider the so‐called multi‐stage hybrid acoustic approach, in which a given noise problem is simulated via a sequence of weakly coupled computations of noise generation and acoustic propagation stages, wherein the simulation of the propagation stage is based on advanced Computational AeroAcoustics (CAA) techniques. The paper introduces an original forcing technique, namely, the Non‐Reflective Interface (NRI), to enable the transfer of an acoustic signal from an a priori noise generation stage into a CAA‐based acoustic propagation phase. Unlike most existing forcing techniques, the NRI is non‐reflective (or anechoic) in nature and, therefore, can properly handle the backscattering effects arising during the noise propagation stage. This attribute makes the NRI‐based weak‐coupling procedure and the associated CAA‐based hybrid approach compatible with a larger variety of realistic noise problems (such as those involving installed configurations in wind tunnel experiments, for instance). The NRI technique is first validated via several test cases of increasing complexity and is then applied to two aerodynamic noise problems. Copyright © 2015 John Wiley & Sons, Ltd. Copyright © 2015 John Wiley & Sons, Ltd.     

水下声波传播过程的时域间断Galerkin有限元方法

本文构建了声压波动方程的改进时域间断Galerkin有限元方法.传统时域连续有限元方法在计算高梯度、强间断特征水中声波传播问题时往往会出现虚假数值振荡现象,这些数值振荡会影响正常波动的计算精度.为了解决这一问题,本文通过引入人工阻尼的方式构建了改进的时域间断Galerkin有限元方法,并针对具有高梯度、强间断特征的多障...     

两种不同弦长螺旋桨的风洞试验及分析

本文介绍了两种不同弦长螺旋桨的风洞试验结果。试验在西北工业大学NF-3风洞的三元试验段进行,试验风速分别为V=20m/s、30m/s和40m/s,每个风速下,螺旋桨旋转速度为:900r/min、1200r/min、1500r/min、1800r/min、2100r/min、2400r/min、2700r/min、3000r/min、3300r/min、3600和3900r/min。风洞试验结果表明:当两种螺旋桨的翼型相同、桨叶角沿径向分布相同,但弦长沿径向的分布不同时,它们的拉力、扭矩、功率、效率以及前进比会有显著区别。其中,1#螺旋的弦长大于2#螺旋桨,在相同的试验风速和螺旋桨转速下,1#螺旋桨的拉力、扭矩和功率高于2#螺旋桨,但效率低于2#螺旋桨。     

高速列车头型长细比对气动噪声的影响

高速列车的头尾车外形对气动噪声具有重要的影响.工程实践中随着车速的增加,车辆头部越来越细长,日本高速磁悬浮列车实践中甚至出现了具有极端长细比的头部形状.本文以讨论头型长细比对列车气动噪声的影响规律为出发点,应用非线性声学求解器(NLAS)和FW-H声学比拟法的混合算法,在3种运行速度下对基于CRH380A高速列车头型概化的4种不同头型长细比的模型车的气动噪声进行了数值模拟.给出了不同头型长细比列车的流场特征、气动阻力和气动噪声.结果表明,列车的气动总阻力随头型长细比的增大而减小,且头型长细比对列车总气动阻力的影响随运行速度的增加而增强.而头型长细比对气动噪声的影响呈现出较为复杂的影响,并不存在单调的影响关系;综合考虑气动阻力和气动噪声,长细比最大的头型综合性能较优,但差异并不显著,因此在不考虑微气压波等因素的条件下,简单增加车头长细比并不一定能带来明显的气动噪声性能提升.     

Numerical simulation of flow‐induced noise using LES/SAS and Lighthill's acoustic analogy

We present a finite element (FE) formulation of Lighthill's acoustic analogy for the hybrid computation of noise generated by turbulent flows. In the present approach, the flow field is computed using large eddy simulation and scale adaptive simulation turbulence models. The acoustic propagation is obtained by solving the variational formulation of Lighthill's acoustic analogy with the FE method. In order to preserve the acoustic energy, we compute the inhomogeneous part of Lighthill's wave equation by applying the FE formulation on the fine flow grid. The resulting acoustic nodal loads are then conservatively interpolated to the coarser acoustic grid. Subsequently, the radiated acoustic field can be solved in both time and frequency domains. In the latter case, an enhanced perfectly matched layer technique is employed, allowing one to truncate the computational domain in the acoustic near field, without compromising the numerical solution. Our hybrid approach is validated by comparing the numerical results of the acoustic field induced by a corotating vortex pair with the corresponding analytical solution. To demonstrate the applicability of our scheme, we present full 3D numerical results for the computed acoustic field generated by the turbulent flow around square cylinder geometries. The sound pressure levels obtained compare well with measured values. Copyright © 2009 John Wiley & Sons, Ltd.     

Finite element analysis of air flow around permeable sand fences

A numerical study of the turbulent air flow in a trench trap and the turbulent flow around a permeable sand fence is reported in this paper. The two-dimensional modified k–ε turbulence model proposed by Kato and Launder is used to predict the turbulent characteristics of the air flow. The discretization method for the governing equations is the three-step Taylor/Galerkin finite element method proposed by the authors. For the flow in a trench trap the numerical results are compared with experimental data obtained under realistic conditions using a large wind tunnel. For the air flow around a permeable sand fence a pressure loss model is used to represent the effect of the porosity of the fence on the flow field. © 1997 John Wiley & Sons, Ltd.     

Finite Difference Method for the Acoustic Radiation of an Elastic Plate Excited by a Turbulent Boundary Layer: A Spectral Domain Solution

A finite difference method is developed to study, on a two-dimensional model, the acoustic pressure radiated when a thin elastic plate, clamped at its boundaries, is excited by a turbulent boundary layer. Consider a homogeneous thin elastic plate clamped at its boundaries and extended to infinity by a plane, perfectly rigid, baffle. This plate closes a rectangular cavity. Both the cavity and the outside domain contain a perfect fluid. The fluid in the cavity is at rest. The fluid in the outside domain moves in the direction parallel to the system plate/baffle with a constant speed. A turbulent boundary layer develops at the interface baffle/plate. The wall pressure fluctuations in this boundary layer generates a vibration of the plate and an acoustic radiation in the two fluid domains. Modeling the wall pressure fluctuations spectrum in a turbulent boundary layer developed over a vibrating surface is a very complex and unresolved task. Ducan and Sirkis [1] proposed a model for the two-way interactions between a membrane and a turbulent flow of fluid. The excitation of the membrane is modeled by a potential flow randomly perturbed. This potential flow is modified by the displacement of the membrane. Howe [2] proposed a model for the turbulent wall pressure fluctuations power spectrum over an elastomeric material. The model presented in this article is based on a hypothesis of one-way interaction between the flow and the structure: the flow generates wall pressure fluctuations which are at the origin of the vibration of the plate, but the vibration of the plate does not modify the characteristics of the flow. A finite difference scheme that incorporates the vibration of the plate and the acoustic pressure inside the fluid cavity has been developed and coupled with a boundary element method that ensures the outside domain coupling. In this paper, we focus on the resolution of the coupled vibration/interior acoustic problem. We compare the results obtained with three numerical methods: (a) a finite difference representation for both the plate displacement and the acoustic pressure inside the cavity; (b) a coupled method involving a finite difference representation for the displacement of the plate and a boundary element method for the interior acoustic pressure; (c) a boundary element method for both the vibration of the plate and the interior acoustic pressure. A comparison of the numerical results obtained with two models of turbulent wall pressure fluctuations spectrums - the Corcos model [3] and the Chase model [4] - is proposed. A difference of 20 dB is found in the vibro-acoustic response of the structure. In [3], this difference is explained by calculating a wavenumber transfer function of the plate. In [6], coupled beam-cavity modes for similar geometry are calculated by the finite difference method. This revised version was published online in July 2006 with corrections to the Cover Date.