跨声速机翼阵风减缓研究

阵风响应分析是大型民用飞机设计必不可少的工作. 利用操纵面的主动偏转实现机翼阵风减缓是未来民用飞行器的一个关键技术. 基于CFD/CSD耦合的气动弹性仿真方法,将阵风视为输入,翼根弯矩作为输出,通过系统辨识方法建立跨音速阵风响应的状态空间分析模型. 而后将副翼作动位移视为系统输入,建立副翼作动对应的机翼响应分析模型. 耦合上述2个模型,通过最优控制方法设计副翼偏转的控制律,实现跨音速机翼的阵风响应减缓. 通过设计状态观测器得到最优控制反馈所需的状态量. 通过数值算例验证了所设计的阵风减缓控制律的有效性,能将翼根弯矩减少60%～80%.      

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

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

Modernization and optimization of a legacy open-source CFD code for high-performance computing architectures

Legacy codes remain a crucial element of today's simulation-based engineering ecosystem due to the extensive validation process and investment in such software. The rapid evolution of high-performance computing architectures necessitates the modernization of these codes. One approach to modernization is a complete overhaul of the code. However, this could require extensive investments, such as rewriting in modern languages, new data constructs, etc., which will necessitate systematic verification and validation to re-establish the credibility of the computational models. The current study advocates using a more incremental approach and is a culmination of several modernization efforts of the legacy code MFIX, which is an open-source computational fluid dynamics code that has evolved over several decades, widely used in multiphase flows and still being developed by the National Energy Technology Laboratory. Two different modernization approaches,‘bottom-up’ and ‘top-down’, are illustrated. Preliminary results show up to 8.5x improvement at the selected kernel level with the first approach, and up to 50% improvement in total simulated time with the latter were achieved for the demonstration cases and target HPC systems employed.     

Portable implementation model for CFD simulations. Application to hybrid CPU/GPU supercomputers

Nowadays, high performance computing (HPC) systems experience a disruptive moment with a variety of novel architectures and frameworks, without any clarity of which one is going to prevail. In this context, the portability of codes across different architectures is of major importance. This paper presents a portable implementation model based on an algebraic operational approach for direct numerical simulation (DNS) and large eddy simulation (LES) of incompressible turbulent flows using unstructured hybrid meshes. The strategy proposed consists in representing the whole time-integration algorithm using only three basic algebraic operations: sparse matrix–vector product, a linear combination of vectors and dot product. The main idea is based on decomposing the nonlinear operators into a concatenation of two SpMV operations. This provides high modularity and portability. An exhaustive analysis of the proposed implementation for hybrid CPU/GPU supercomputers has been conducted with tests using up to 128 GPUs. The main objective consists in understanding the challenges of implementing CFD codes on new architectures.     

混合驱动水下滑翔机水动力参数辨识

与传统水下滑翔机相比,混合驱动水下滑翔机可以分别利用头部的浮力驱动单元和机身尾部的螺旋桨推进单元进行驱动从而实现不同形式的运动,具备低功耗、长航程、良好机动性等特点,具有广泛的应用前景. 准确的动力学模型以及精确的水动力参数是实现混合驱动水下滑翔机控制系统设计以及精确导航的基础. 在混合驱动水下滑翔机动力学模型已知的前提下获得准确的水动力参数是本文主要研究的问题. 本文以天津大学研制的混合驱动水下滑翔机“海燕” 作为研究对象,提出一种在有限航行参数条件下,基于大数据统计分析的计算流体力学(computational fluid dynamics, CFD) 和参数辨识相结合来获取混合驱动水下滑翔机的水动力参数的方法. 即首先建立滑翔机的动力学模型,推导出稳态数据与所求水动力参数的关系;然后采用CFD 的方法得到其升力系数,根据大量稳态纯滑翔实验数据,结合大数据统计分析,辨识出其剩余水动力参数;最后,根据混合驱动模式下的实验数据辨识出与螺旋桨相关的参数,从而得到其整套的水动力参数. 该方法不仅结合了CFD 方法具有获取复杂外形结构航行水动力的特点,而且可以有效利用大量现场实验数据,因而能够更加准确地辨识其实际运动. 通过运动仿真与试验对比,验证了该辨识方法的正确性和有效性,对滑翔机的研究具有指导意义.      

基于医学影像的血流动力学分析

随着计算机成像技术的发展,计算机断层扫描(CT)和核磁共振(MRI)已经广泛应用于临床,特别是血管类疾病的诊断,比如动脉狭窄、血管瘤、血管畸形等。除了可以提供高分辨率的静态图像,先进的MRI技术还可以通过时间序列直接反映血流动力学的变化。而基于计算机成像和三维重建技术,血流动力学的参数又可以通过计算流体力学的方法进行详细的分析。如何将血流动力学参数和临床诊断相结合是近年来在转化医学领域研究的热点。文中将结合文献调研和作者自己的研究工作对基于医学图像的血流动力学分析进行综述,并探讨未来的研究方向。     

轴流泵参数化设计及应力与模态特征分析

以船用轴流式喷水推进泵为对象,探索了轴流泵参数化设计、水动力性能、静强度和结构声学特征分析的数值途径。轴流泵叶轮采用升力法设计,导叶采用流线法设计,叶片三维造型在NUMECA参数化设计平台中完成。轴流泵水动力性能校核由粘性CFD计算完成,CFD计算同时提取得到叶片分布式水动力载荷。叶片静强度校核由ANSYS有限元计算叶片应力和应变特征完成,应力分析时同时考虑水动力载荷、重力载荷和离心力载荷。叶片结构声学特征分析由NASTRAN有限元计算叶片模态振型和振型频率完成。计算结果表明,轴流泵扬程和功率满足设计指标,效率达87.13%;叶轮叶片形变相对于叶顶间隙来说为极小量,可忽略不计,叶片存在局部应力集中现象,最大应力小于许用应力,满足静强度要求;叶片前四阶振型特征与分析经验一致,且振型频率远离轴频、叶频及其谐频特征频率,能够避免共振产生。     

Parallel computation of turbulent flows using moment base lattice Boltzmann method

This paper describes parallel computing approach for simulating turbulent flows using a moment base lattice Boltzmann method. The distribution functions of the lattice Boltzmann method are expressed by corresponding moments. Choosing proper relaxation times for higher order moments, a minimum numerical dissipation is implicitly added to stabilise the method at high Reynolds numbers. Validation of the method is made by computing free decaying periodic turbulent flows and fully developed turbulent channel flows on a GPU platform. Though the present method requires additional work to calculate the higher order moments, it is shown that additional computational cost is negligible in the GPU computing. The numerical results stably obtained for the turbulent flows are in good agreement with those of a pseudo-spectral method and corresponding DNS database.     

Numerical study on the suppression of the vortex-induced vibration of an elastically mounted cylinder by a traveling wave wall

In the present paper, the commercial CFD code “Fluent” was employed to perform 2-D simulations of an entire process that included the flow around a fixed circular cylinder, the oscillating cylinder (vortex-induced vibration, VIV) and the oscillating cylinder subjected to shape control by a traveling wave wall (TWW) method. The study mainly focused on using the TWW control method to suppress the VIV of an elastically supported circular cylinder with two degrees of freedom at a low Reynolds number of 200. The cross flow (CF) and the inline flow (IL) displacements, the centroid motion trajectories and the lift and drag forces of the cylinder that changed with the frequency ratios were analyzed in detail. The results indicate that a series of small-scale vortices will be formed in the troughs of the traveling wave located on the rear part of the circular cylinder; these vortices can effectively control the flow separation from the cylinder surface, eliminate the oscillating wake and suppress the VIV of the cylinder. A TWW starting at the initial time or at some time halfway through the time interval can significantly suppress the CF and IL vibrations of the cylinder and can remarkably decrease the fluctuations of the lift coefficients and the average values of the drag coefficients; however, it will simultaneously dramatically increase the fluctuations of the drag coefficients.