基于CFD/CSD分区耦合的气动弹性数值模拟

提出了一种改进的常体积转换局部插值方法和判断两套网格之间宿主-受体关系的高效算法。结构计算采用带面内旋转自由度的四节点四边形平板壳单元模型,气动计算基于有限体积离散的任意拉格朗日-欧拉(ALE)格式N-S方程,并采用Delaunay图映射法来处理网格移动。根据分区耦合原理,建立一套基于计算流体力学(CFD)/计算结构力学(CSD)的气动弹性计算方法,并应用于AGARD445.6机翼颤振速度系数和颤振频率系数、切尖三角翼的翼尖颤振幅值和颤振频率的分析中。与已有文献相比,数值计算结果误差分别为7.85%和9.46%、7.28%和1.03%,相对较小,验证了该气动弹性分析方法和程序的可靠性和可行性。     

基于聚合多重网格方法的跨音速颤振的数值模拟

采用流固耦合方法对跨音速颤振进行了数值模拟。流体方面在非结构网格上用有限体积方法求解了Euler方程;结构方面则求解了后掠机翼典型剖面的结构模态方程。时间推进采用双时间步长：对每一真实时间步,都通过基于聚合多重网格方法的伪时间步推进,对流体和结构方程交替迭代．得到一个稳态的流固耦合的解。文章最后给出了NACA64A010翼型剖面的跨音速颤振边界．与相关文献的计算结果符合良好。     

基于ALE有限元法的流固耦合强耦合数值模拟

针对不同流固耦合问题,提出一种基于任意拉格朗日-欧拉(ALE)有限元技术的分区强耦合算法.运用半隐式特征线分裂算法求解ALE描述下的不可压缩黏性流体Navier-Stokes方程.分别考虑一般平面运动刚体和几何非线性固体,采用复合隐式时间积分法推进结构运动方程,故可选用较大时间步长;进一步应用单元型光滑有限元法求解几何非线性固体大变形,获得更精确结构解且不影响计算效率.运用子块移动技术结合正交-半扭转弹簧近似法高效更新流体动网格;同时将一质量源项引入压力泊松方程满足几何守恒律,无需复杂构造网格速度差分格式.采用简单高效的固定点法配合Aitken动态松弛技术实现各场耦合,可灵活选择先进单场求解技术,具备较好程序模块性.运用本文算法分别模拟了H型桥梁截面颤振问题和均匀管道流内节气阀涡激振动问题.研究表明,数值结果与已有文献数据吻合,计算精度和求解效率均令人满意.     

探测器对超音速刚性盘-缝-带型降落伞系统的影响

研究了流体初始马赫数为 2.0 时, 探测器的存在与否对刚性盘-缝-带型降落伞系统气动减速性能以及流场流体结构特性的影响. 对于非定常可压缩流体的数值模拟, 流场采用了三层块结构自适应网格加密技术, 配合混合形式的TCD (tuned center difference)和WENO (weighted essentially non-oscillatory)计算格式以及基于拉伸涡亚格子模型的大涡模拟方法来处理超音速流中的激波以及大尺度湍流旋涡结构等. 结果表明: 无探测器时, 降落伞系统的流场结构稳定, 扰动较小; 有探测器存在时, 探测器后端的湍流尾迹和伞衣内部逆向运动溢出的流体与伞衣前端的弓形激波周期性的相互作用, 使得激波位置发生前移、激波倾角变小, 伞衣内部流场难以达到平衡稳定状态. 这加剧了降落伞系统的气动阻力振荡脉动变化, 降低了降落伞系统气动阻力系数, 同时也使得降落伞系统流场尾迹结构更加复杂.      

非线性流体-刚体结构相互作用问题的一种数值模拟方法

给出了一种模拟非线性流体-刚体结构相互作用问题的数值方法．文中假定结构承受大的刚体运动,流体流动受非线性有粘或无粘的场方程支配并满足自由表面和两相耦合界面上的非线性边界条件,利用任意拉氏-欧氏（ALE）网格系统构造了数值模型．采用所探讨的多块数值格式,允许可动重造网格间有独立的相对运动,从而克服了流体网格与固体大运动匹配的困难．通过数值离散化,导出了描述非线性流固耦合动力学的数值方程并应用耦合迭代过程对其作了求解．通过算例,说明了所提出数值模型的应用．     

大跨度桥梁多模态耦合颤振的TMD控制

推导了装有TMD的结构在气动自激力作用下的动力微分方程,基于模态空间中多模态耦合颤振分析手段,运用考虑安装TMD的多模态自动分析法对结构-TMD系统进行了颤振分析和TMD控制分析,使多模态自动分析法能适用于TMD颤振控制分析,避免了双参数搜索和迭代计算,提高了计算效率。对某在建三塔悬索桥进行了原结构颤振分析和TMD-结构颤振频域分析,探讨TMD控制参数对颤振临界风速的影响。     

计算气动弹性力学中的界面映射方法研究

非线性气动弹性体振动研究中,涉及到非线性的结构动力学和非线性的流体动力学耦合问题,在耦合边界上要满足两个系统的连续性相容条件,必须在边界处进行数据的交换。本文针对非线性气动弹性问题的计算流体动力学(CFD)和计算结构动力学(CSD)的耦合计算方法,在常体积转换法(CVT)的基础上,发展了一种耦合界面的数据映射矩阵(IMM)。该方法仅需要局部的网格信息,将耦合边界上载荷信息和位移信息的转换放在同一个映射矩阵中来处理,并且该矩阵可以通用求解CFD/CSD的耦合问题,克服了占用大量CPU时间和内存的弊端。最后将该界面映射方法应用于柔性大展弦比机翼的气动弹性计算和AGARD445.6机翼的颤振预测中,结果表明该方法能够高效、高精度地处理不同网格体系间的数据交换,并具有处理复杂非规则几何体信息转换的能力。     

基于ALE有限元法的流固耦合强耦合数值模拟

针对不同流固耦合问题,提出一种基于任意拉格朗日--欧拉(ALE)有限元技术的分区强耦合算法. 运用半隐式特征线分裂算法求解ALE描述下的不可压缩黏性流体Navier-Stokes方程. 分别考虑一般平面运动刚体和几何非线性固体,采用复合隐式时间积分法推进结构运动方程,故可选用较大时间步长;进一步应用单元型光滑有限元法求解几何非线性固体大变形,获得更精确结构解且不影响计算效率. 运用子块移动技术结合正 交--半扭转弹簧近似法高效更新流体动网格;同时将一质量源项引入压力泊松方程满足几何守恒律,无需复杂构造网格速度差分格式. 采用简单高效的固定点法配合Aitken动态松弛技术实现各场耦合,可灵活选择先进单场求解技术,具备较好程序模块性. 运用本文算法分别模拟了H型桥梁截面颤振问题和均匀管道流内节气阀涡激振动问题. 研究表明,数值结果与已有文献数据吻合,计算精度和求解效率均令人满意.      

流固耦合分析的一种改进CBS有限元算法

针对流固耦合问题, 发展了一种基于任意拉格朗日-欧拉(ALE)描述有限元法的弱耦合分区算法. 运用半隐式特征线分裂算法求解Navier-Stokes方程, 在压力Poisson 方程中引入质量源项以满足几何守恒律; 运用子块移动技术更新动态网格, 并配以光滑处理防止网格质量下降; 采用Newmark-β 法求解结构运动方程. 为保持流体-结构界面处速度和动量守恒, 利用修正结合界面边界条件方法求解界面处速度通量和动量通量. 运用本方法分别模拟了不同雷诺数下单圆柱横向和两向流致振动、串列双圆柱两向流致振动. 计算表明, 本文方法计算效率高, 计算结果与已有实验和数值计算数据吻合.     

Developments of numerical methods for linear and nonlinear fluid-solid interaction dynamics with applications

本文综述了线性与非线性流固耦合问题数值方法的进展及工程应用. 讨论了四种数值分析方法: (1) 混合有限元-子结构-子区域数值模型, 以求解有限域线性流固耦合问题, 如流体晃动, 声腔-结构耦合, 流体中的压力波, 化工容器的地震响应,坝水耦合等; (2) 混合有限元-边界元数值模型, 以求解涉及无限域的线性流固耦合问题, 如大型浮体承受飞机降落冲击, 船舰的炮击回应等; (3) 混合有限元-有限差分(体积) 数值模型, 以求解不涉及破浪和两相分离的非线性流固耦合问题; (4) 混合有限元-光滑粒子数值模型, 以求解涉及破浪和两相分离的非线性流固耦合问题. 文中推荐分区迭代求解过程, 以便应用现有的固体及流体求解器, 于毎一时间步长分别求解固体及流体的方程, 通过耦合迭代收敛, 向前推进以达问题求解. 文中选用的工程应用例子包含气-液-壳三相耦合, 液化天然气船水晃动, 人体步行冲击引起的声腔-建筑结构耦合, 大型浮体承受飞机降落冲击的瞬态动力回应, 涉及破浪和两相分离的气-翼耦合及结构于水上降落的冲击. 数值分析结果与可用的实验或计算结果作了比较, 以说明所述方法的精度及工程应用价值. 文中列出了基于流固耦合的波能采积装置模型, 以应用线性系统的共振及非线性系统的周期解原理, 有效地采积波能. 本文列出了231 篇参考文献, 以便读者进一步研讨所感兴趣方法.     

DEVELOPMENT OF A THREE-DIMENSIONAL VISCOUS AEROELASTIC SOLVER FOR NONLINEAR PANEL FLUTTER

A new three-dimensional (3-D) viscous aeroelastic solver for nonlinear panel flutter is developed in this paper. A well-validated full Navier–Stokes code is coupled with a finite-difference procedure for the von Karman plate equations. A subiteration strategy is employed to eliminate lagging errors between the fluid and structural solvers. This approach eliminates the need for the development of a specialized, tightly coupled algorithm for the fluid/structure interaction problem. The new computational scheme is applied to the solution of inviscid two-dimensional panel flutter problems for subsonic and supersonic Mach numbers. Supersonic results are shown to be consistent with the work of previous researchers. Multiple solutions at subsonic Mach numbers are discussed. Viscous effects are shown to raise the flutter dynamic pressure for the supersonic case. For the subsonic viscous case, a different type of flutter behavior occurs for the downward deflected solution with oscillations occurring about a mean deflected position of the panel. This flutter phenomenon results from a true fluid/structure interaction between the flexible panel and the viscous flow above the surface. Initial computations have also been performed for inviscid, 3-D panel flutter for both supersonic and subsonic Mach numbers.     

Computational Investigation of Wing Flutter

Numerical aeroelastic simulation of a high-aspect-ratio transport type wing model in transonic region is presented. The aeroelastic responses of the wing are extracted by integrating compressible thin-layer Navier-Stokes equations coupled with the equations of motion of the wing structure, in a time dependent manner. The Yee-Harten implicit TVD scheme and the Wilson's θ method are employed to integrate these equations, respectively. Flutter boundaries were found for Mach number range, 0.7 to 0.85 and the results were compared with experimental flutter boundaries. Futhermore, Limit Cycle Oscillations were found and the characteristics of the flutter and limit cycle oscillations are investigated and discussed.     

Validation of CFD‐CSD coupling interface methodology using commercial codes

Coupling interface between computational fluid dynamics (CFD) and computational structural dynamics (CSD) is required to provide exchange of information for the simulation of fluid–structure interaction (FSI) phenomena. Accuracy and consistency of information exchanged through coupling interface between the independent CFD and CSD solvers plays a central role in the simulation and prediction of FSI phenomenon, like flutter. In this paper validation of an implemented coupling interface methodology is presented for subsonic, transonic and near supersonic mach regime. The test case chosen for this purpose is the flutter of AGARD445.6 standard I‐wing weakened model configuration for subsonic to near transonic flow regime. Gambit^® and Fluent^® are used for CFD grid generation and solution of fluid dynamic equations, respectively. CSD modeling and simulation are provided by numerical time integration of modal dynamic equations derived through the finite element modeling in ANSYS^® environment. Copyright © 2009 John Wiley & Sons, Ltd.     

A hybrid immersed boundary and material point method for simulating 3D fluid–structure interaction problems

A numerical method is developed for solving the 3D, unsteady, incompressible Navier–Stokes equations in curvilinear coordinates containing immersed boundaries (IBs) of arbitrary geometrical complexity moving and deforming under forces acting on the body. Since simulations of flow in complex geometries with deformable surfaces require special treatment, the present approach combines a hybrid immersed boundary method (HIBM) for handling complex moving boundaries and a material point method (MPM) for resolving structural stresses and movement. This combined HIBM & MPM approach is presented as an effective approach for solving fluid–structure interaction (FSI) problems. In the HIBM, a curvilinear grid is defined and the variable values at grid points adjacent to a boundary are forced or interpolated to satisfy the boundary conditions. The MPM is used for solving the equations of solid structure and communicates with the fluid through appropriate interface‐boundary conditions. The governing flow equations are discretized on a non‐staggered grid layout using second‐order accurate finite‐difference formulas. The discrete equations are integrated in time via a second‐order accurate dual time stepping, artificial compressibility scheme. Unstructured, triangular meshes are employed to discretize the complex surface of the IBs. The nodes of the surface mesh constitute a set of Lagrangian control points used for tracking the motion of the flexible body. The equations of the solid body are integrated in time via the MPM. At every instant in time, the influence of the body on the flow is accounted for by applying boundary conditions at stationary curvilinear grid nodes located in the exterior but in the immediate vicinity of the body by reconstructing the solution along the local normal to the body surface. The influence of the fluid on the body is defined through pressure and shear stresses acting on the surface of the body. The HIBM & MPM approach is validated for FSI problems by solving for a falling rigid and flexible sphere in a fluid‐filled channel. The behavior of a capsule in a shear flow was also examined. Agreement with the published results is excellent. Copyright © 2007 John Wiley & Sons, Ltd.     

A strongly coupled,embedded-boundary method for fluid–structure interactions of elastically mounted rigid bodies

In the present paper, an embedded-boundary formulation that is applicable to fluid–structure interaction problems is presented. The Navier–Stokes equations for incompressible flow are solved on a Cartesian grid which is not aligned with the boundaries of a body that undergoes large-angle/large-displacement rigid body motions through the fixed grid. A strong-coupling scheme is adopted, where the fluid and the structure are treated as elements of a single dynamical system, and all of the governing equations are integrated simultaneously and interactively in the time domain. A demonstration of the accuracy and efficiency of the method is given for a variety of fluid–structure interaction problems.     

Die-swell,splashing drop and a numerical technique for solving the Oldroyd B model for axisymmetric free surface flows

This work deals with the development of a numerical method for simulating viscoelastic axisymmetric free surface flow of an Oldroyd B fluid. A novel formulation is developed for the computation of the non-Newtonian extra-stress components on rigid boundaries and on the symmetry axis. The full free surface stress conditions are employed. The resulting governing equations are solved by finite differences on a Marker-and-cell (MAC) type grid. Validation is provided by simulating a pipe flow problem. The classical die-swell problem is solved and swelling ratios are provided. The height of the splash caused by a falling liquid drop for various Reynolds and Weissenberg numbers is then studied, and the height of the splash is shown to diminish with increasing viscoelasticity.     

A flexible rod in annular leakage flow: Influence of turbulence and equilibrium offset,and analysis of instability mechanisms

The linear stability of a flexible, cylindrical rod subjected to annular leakage flow is studied. The mathematical models developed by Li, Kaneko, and Hayama in 2002 and Fujita and Shintani in 2001 are bridged and extended, to account for a flexible rod with equilibrium offset (eccentricity) in laminar or turbulent leakage flow. Stability characteristics are analyzed numerically for a variety of configurations. It is found that simply supported rods may become unstable at a certain critical flow speed by either divergence or flutter, depending on dimensions and fluid/solid properties. It is furthermore found that the critical flow speed is quite insensitive to use of a laminar friction model at high Reynolds numbers in cases of divergence, but sensitive to it in cases of flutter. These findings are verified analytically though analysis of an energy equation. This equation shows that (i) divergence instability is independent of fluid friction; (ii) flutter instability is caused solely by fluid friction. It also suggests a possible explanation to the question of why a ‘wrong’ fluid friction assumption gives a too large critical flow speed in cases of flutter instability at a high Reynolds number.     

Cyber-physical flexible wing for aeroelastic investigations of stall and classical flutter

The fluid–structure interactions of a finite aspect ratio, cantilevered, flexible wing were investigated using a cyber-physical system to virtually augment the torsional dynamics of the wing. Cyber-physical systems (CPS), which have recently been pursued by a small number of research groups, have proven to be a very useful mechanism to interrogate the fluid–structure interaction parameter space. The premise of a CPS is to use dynamic feedback control to make a system behave according to desired equations of motion. Systems are composed of embedded hardware and software coupled with real-time computing to give the user the flexibility to quickly explore a range of structural parameters. With the advancement of modern control theory, robotics and embedded systems, CPSs integrate both simulation and physical properties into a smart structure which can be used to push the boundaries of research investigations.The CPS in this work allows for the investigation of dynamic aeroelastic instabilities of a three-dimensional, flexible, rectangular planform wing. Two dynamic instability regimes are observed: first, stall flutter, in which the torsional or pitch mode is excited through the dynamic stall process, and second, coupled (or classical) flutter, in which the pitch mode couples with the bending mode. By varying the torsional stiffness and therefore the frequency of torsional versus bending oscillations of the wing, both of these regimes can be attained at the same aerodynamic conditions using the CPS.     

Flutter and resonances of a flag near a free surface

We investigate the effects of a nearby free surface on the stability of a flexible plate in axial flow. Confinement by rigid boundaries is known to affect flag flutter thresholds and fluttering dynamics significantly, and this work considers the effects of a more general confinement involving a deformable free surface. To this end, a local linear stability is proposed for a flag in axial uniform flow and parallel to a free surface, using one-dimensional beam and potential flow models to revisit this classical fluid–structure interaction problem. The physical behaviour of the confining free surface is characterized by the Froude number, corresponding to the ratio of the incoming flow velocity to that of the gravity waves. After presenting the simplified limit of infinite span (i.e. two-dimensional problem), the results are generalized to include finite-span and lateral confinement effects. In both cases, three unstable regimes are identified for varying Froude number. Rigidly-confined flutter is observed for low Froude number, i.e. when the free surface behaves as a rigid wall, and is equivalent to the classical problem of the confined flag. When the flow and wave velocities are comparable, a new instability is observed before the onset of flutter (i.e. at lower reduced flow speed) and results from the resonance of a structural bending wave and one of the fundamental modes of surface gravity waves. Finally, for large Froude number (low effect of gravity), flutter is observed with significant but passive deformation of the free surface in response of the flag’s displacement.     

Ideal gas flows with separation zones and time-dependent contact discontinuities of complicated shape

Some examples of flows with separation zones andmovable contact discontinuities obtained as a result of the numerical integration of the time-dependent equations for an ideal gas are presented. The examples concern a steady annular separation zone on the blunt nose of a body in a supersonic flow, periodic shedding of unsteady discontinuities from a cylinder in a steady uniform subsonic flow with a supercritical Mach number, and the complicated deformation of a contact (tangential) discontinuity, namely, the boundaries of a two-dimensional jet, either subsonic or supersonic, flowing into a cocurrent subsonic low-velocity flow. A multiple increase in the difference grid capacity in the numerical integration of the Euler equations indicates the absence of a noticeable scheme viscosity effect in the examples calculated. The inviscid nature of the separation flows obtained is also confirmed by their well-known counterparts constructed in the ideal incompressible fluid approximation. The time-average velocity fields of the two-dimensional jet and the intensity of its sound field are in reasonable agreement with the available data.