Unsteady rotor flow analysis using a diagonally implicit harmonic balance method and an overset mesh topology

The diagonally implicit harmonic balance method is developed in an overset mesh topology and applied to unsteady rotor flows analysis. Its efficiency is by reducing the complexity of a fully implicit harmonic balance method which becomes more flexible in handling the higher harmonics of the flow solutions. Applied to the overset mesh topology, the efficiency of the method becomes greater by reducing the number of solution interpolations required during the entire solution procedure as the method reduces the unsteady computation into periodic steady state. To verify the accuracy and efficiency of the method, both hovering and unsteady forward flight of Caradonna and Tung and AH-1G rotors are solved. Compared with wind-tunnel experiments, the numerical results demonstrate good agreements at computational cost an order of magnitude more efficient than the conventional time-accurate computation method. The proposed method has great potential in other engineering applications, including flapping wing vehicles, turbo-machinery, wind-turbines, etc.     

谐波平衡法在动导数快速预测中的应用研究

谐波平衡法以傅里叶级数展开为基础,将周期性非定常流场的非定常求解过程转化为几个定常流场的耦合求解过程,并通过重建得到整个流场的非定常过程. 建立了基于谐波平衡法的动导数快速预测方法,数值模拟了超声速带翼导弹俯仰的动态流场,并通过积分法获取了俯仰动导数,与实验结果吻合很好;且在同等计算精度下,谐波平衡法的计算效率是双时间步方法的13 倍. 应用谐波平衡法研究了较大范围内减缩频率对俯仰动导数的影响规律. 研究发现,对于本外形,当减缩频率降低到一定值后,俯仰动导数的值迅速变化,甚至发生变号;对此现象产生的原因进行了深入分析,并通过对导弹自激俯仰运动的数值模拟验证了该结果. 此外,针对大攻角条件下动态流场非线性强的特点,开展了谐波平衡法在大攻角下的适用性研究. 结果表明,谐波平衡法在大攻角下也能取得很好的计算结果.      

基于CFD的气动力建模及其在气动弹性中的应用

CFD技术为带有气动力非线性的气动弹性分析提供了一种研究途径,但是基于CFD的气动弹性直接数值模拟方法的计算量很大,不便开展定性分析和参数设计.基于CFD的非定常气动力模型的降阶技术为缓解计算效率与计算精度之间的矛盾以及系统的复杂性与易分析、易设计性之间的矛盾提供了行之有效的技术途径.综述了近年来发展的两类基于CFD技术的非定常气动力降阶技术和一种非线性气动力的谐波平衡方法,以及这些方法在非线性气动弹性研究中的运用.对比了各种方法的优越性并作了进一步的展望.      

Cascade aerodynamic gust response including steady loading effects

To predict the unsteady convected gust aerodynamic response of a cascade comprised of arbitrary thick and cambered aerofoils in an incompressible, inviscid, flow field, a complete first-order model is formulated. The flow is analysed by considering a periodic flow channel. The velocity potential is separated into steady and unsteady harmonic components, each described by a Laplace equation. The strong dependence of the unsteady aerodynamics on the steady effects of aerofoil and cascade geometry and incidence angle is manifested in the coupling of the unsteady and steady flow fields through the unsteady boundary conditions. Analytical solutions in individual grid elements of a body-fitted computational grid are then determined, with the complete solution obtained by assembly of these local solutions. The validity and capabilities of this model and solution technique are then demonstrated by analysing the steady and unsteady aerodynamics of both theoretical and experimental cascade configurations.     

高阶谐波平衡方法中非物理解来源分析及改进方法研究

对于周期性非定常问题,高阶谐波平衡（High-order Harmonic Balance, HOHB）方法将非定常方程的解用Fourier 级数展开至一定阶次,从而消除其中的时间导数项,大大降低了计算消耗. 本文以达芬振子方程为例,探讨了HOHB 方法中非物理解的来源,分析结果表明:非物理解出现的原因是在推导过程中非线性项的简化处理导致方程左右两边并不严格相等. 根据非线性项的特点,在其处理过程中扩充子时间层上的时域解,并将非线性项中出现的更高阶谐波截断,使方程左右两边严格相等. 通过对达芬振子方程进行数值模拟发现:改进方法在消除非物理解的同时,也显著减少了计算所需谐波数. 对比参考文献发现,同阶改进方法的精度和原始谐波平衡方法基本相当,证明了本方法的可行性. 最后将本方法应用于具有立方刚度非线性的气动弹性系统中,验证本方法的工程适用性. 但是,当方程中非线性项较多时,本方法所需要的计算消耗会有所增加.      

Equivalent damping of aeroelastic system of an airfoil with cubic stiffness

This paper presents a mathematical study on the subsonic aerodynamics acting on an airfoil with a cubic stiffness. One portion of aerodynamics is assumed as and replaced by an equivalent damping. Using the harmonic balance method, an equivalent system is deduced and studied by a numerical integration method. Numerical examples show the validity and feasibility of the proposed mathematical treatment of the aerodynamics. It reveals that the unsteady aerodynamics acting on the airfoil can not only affect the pitch stiffness but also result in additional damping.     

An analytical approximate technique for a class of strongly non-linear oscillators

An analytical approximate technique for large amplitude oscillations of a class of conservative single degree-of-freedom systems with odd non-linearity is proposed. The method incorporates salient features of both Newton's method and the harmonic balance method. Unlike the classical harmonic balance method, accurate analytical approximate solutions are possible because linearization of the governing differential equation by Newton's method is conducted prior to harmonic balancing. The approach yields simple linear algebraic equations instead of non-linear algebraic equations without analytical solution. With carefully constructed iterations, only a few iterations can provide very accurate analytical approximate solutions for the whole range of oscillation amplitude beyond the domain of possible solution by the conventional perturbation methods or harmonic balance method. Three examples including cubic-quintic Duffing oscillators are presented to illustrate the usefulness and effectiveness of the proposed technique.     

Simultaneous quantitative flow measurement using PIV on both sides of the air–water interface for breaking waves

An experimental method for simultaneously measuring the velocity fields on the air and water side of unsteady breaking waves is presented. The method includes a novel technique for seeding the air flow such that the air velocity can be resolved in the absence of wind. Low density particles that have large Stokes drag and ability to respond to high-frequency flow fluctuations are used to seed the air flow. Multi-camera, multi-laser particle image velocimetry setups are applied to small-scale unsteady breaking waves, yielding fully time-resolved velocity fields. The surface tension of the fluid is altered and controlled to form spilling breaking waves. Results for the velocity and vorticity fields of representative spilling breakers, which show shedding of an air-side vortex and well-documented generation of water-side vorticity, are presented and discussed.     

The Fourier pseudo-spectral method for the SRLW equation

In this paper we present a Fourier pseudo-spectral method with a restraint operator for the SRLW equation. We prove the stability of the schemes and give optimum error estimates.     

非定常流动模拟的时间离散方法

对于不同非定常流动问题, 采用合适的时间离散方法,可有效提高数值精度和计算效率. 本文在总结传统时间离散方法的基础上,对近些年发展的非线性频域法、谐波平衡法、经典时间谱方法、时间谱元法、时间有限差分法等进行了系统地总结.根据离散形式的不同,将上述方法分为时域推进法、频域谐波法、时域配点法和混合方法4大类.首先简要介绍了各类方法的数学思想以及研究进展,并重点比较了(准)周期性非定常流动计算中各方法的精度、效率以及适用范围.然后, 对各种时间离散格式的特点进行总结,并就不同的非定常流动问题如何选择合适的时间离散方法给予了建议.最后, 对这些新型时间离散格式在工程中的应用进行了简要介绍,并对其发展方向进行展望.      

Time discretization methods in the computation of unsteady flow

对于不同非定常流动问题, 采用合适的时间离散方法,可有效提高数值精度和计算效率. 本文在总结传统时间离散方法的基础上,对近些年发展的非线性频域法、谐波平衡法、经典时间谱方法、时间谱元法、时间有限差分法等进行了系统地总结.根据离散形式的不同,将上述方法分为时域推进法、频域谐波法、时域配点法和混合方法4大类.首先简要介绍了各类方法的数学思想以及研究进展,并重点比较了(准)周期性非定常流动计算中各方法的精度、效率以及适用范围.然后, 对各种时间离散格式的特点进行总结,并就不同的非定常流动问题如何选择合适的时间离散方法给予了建议.最后, 对这些新型时间离散格式在工程中的应用进行了简要介绍,并对其发展方向进行展望.     

Effect of Unsteady Wind Flow on Galloping of Tall Prismatic Structures

The galloping of tall prismatic cantilevertype structures due to unsteady wind is solvedanalytically. The unsteady wind was considered byadding a time varying wind speed component to the meanwind speed component. In reality, the time-varyingwind speed component is a random phenomenon that can bemodeled as a series of harmonic terms using thetransformation of the unsteady wind speed spectruminto the time domain. In doing this it is apparentthat the structure is subjected to multiharmonicexternal and parametric excitations due to theunsteady wind in addition to the nonlinearself-excited wind forces due to the steady wind speedcomponent. To have a clear insight into the unsteadywind effect, only one harmonic term is considered outof all the harmonic terms. The multiple-scale methodis used to study the effect of primary and secondaryresonances on the galloping response of the structure. Comparisons between the analytical results obtainedfrom the method of multiple scales and the numericalsolutions obtained from numerical integration indicate the accuracy of the analysis and thecomprehensive information obtained from the analyticalsolutions.     

Comparison of Symmetric and Asymmetric Schemes with Arithmetic and Harmonic Averaging for Fracture Flow on Cartesian Grids

Performance of four finite-difference schemes for fluid flow in rough-walled fractures on regular Cartesian grids is evaluated numerically. The four schemes are an asymmetric scheme with arithmetic averaging, an asymmetric scheme with harmonic averaging, a symmetric scheme with arithmetic averaging, and a symmetric scheme with harmonic averaging. The schemes are compared with respect to their simulated hydraulic aperture and the mass balance error. 1320 flow simulations with different grid sizes, mean fracture aperture and root mean square (RMS)/mean aperture ratio are completed. The asymmetric scheme with arithmetic averaging arises naturally, without any extra assumptions about the correct transmissivity averaging procedure, when one uses second-order finite differences to approximate the generalized Laplace operator expanded as a derivative of a product. Hydraulic apertures obtained with harmonic averaging are found to usually be smaller than those obtained with arithmetic averaging, especially when the ratio of aperture RMS to the mean aperture is larger. The traditionally used asymmetric schemes are found to be superior to symmetric schemes in terms of mass balance accuracy.

     

Unsteady stagnation-point flow of a Newtonian fluid in the presence of a magnetic field

The unsteady stagnation-point flow of a viscous fluid impinging on an infinite plate in the presence of a transverse magnetic field is examined and solutions are obtained. It is assumed that the infinite plate at y=0 is making harmonic oscillations in its own plane. A finite difference technique is employed and solutions for small and large frequencies of the oscillations are obtained for various values of the Hartmann's number.     

The Simplest Resonance Capture Problem,Using Harmonic Balance Based Averaging

We study the dynamics of capture into, or escape from, resonance in a strongly nonlinear oscillator with weak damping and forcing, using harmonic balance based averaging (HBBA). This system provides the simplest example of resonance capture that we know of. The HBBA technique, here adapted to tackle nonlinear resonances, provides a harmonic balance assisted approximation to the underlying, asymptotically correct, averaged dynamics. Allowing the harmonic balance approximation makes a variety of systems analytically tractable which might otherwise be intractable. The evolution equations for amplitude and phase of oscillations are derived first. Restricting attention near the primary resonance, the slow flow equations are approximately averaged. The resulting flow transparently shows the stable and unstable primary resonant solutions, as well as the trajectories that get captured into resonance and the ones that escape. Good agreement with numerics is obtained, showing the utility of HBBA near resonance manifolds.     

The application of reduced space harmonic balance method for the nonlinear vibration problem in rotor dynamics

The reduced space harmonic balance method is demonstrated to find the maximum vibration responses of rotor systems. Within the reduced space SQP method, transition from the high dimensional optimization space to the desired reduced space is accomplished by resorting to the null space decomposition technique, resulting in the elimination of the harmonic balance constraints. Numerical examples of rotor systems are presented to show the applicability of the proposed methodology.     

Residue harmonic balance approach to limit cycles of non-linear jerk equations

A residue harmonic balance is established for accurately determining limit cycles to parity- and time-reversal invariant general non-linear jerk equations with cubic non-linearities. The new technique incorporates the salient features of both methods of harmonic balance and parameter bookkeeping to minimize the total residue. The residue is separated into two parts in each step; one conforms to the present order of approximation and the remaining part for use in the next order. The corrections are governed by a set of linear ordinary differential equations that can be solved easily. Three specific cases of non-linear jerk equations are given to illustrate the validity and efficiency. The approximations to the angular frequency and the limit cycle are obtained and compared. The results show that the approximations obtained are in excellent agreement with the exact solutions for a wide range of initial velocities. The new technique is simple in principle and can be applied to other non-linear oscillating systems.     

基于时频变换的激波风洞天平信号分析与处理

在激波风洞中开展测力试验时,测力系统在风洞流场起动瞬间会受到冲击激励,从而对天平的输出信号产生惯性干扰.天平输出信号中叠加有动态气动力信号和惯性振动信号,有可能无法直接分辨出气动力信号的规律性,信号处理结果与真实气动力之间会产生较大的误差,导致处理结果不可靠.由于模型测力天平系统结构的复杂性,在极短的有效试验时间(毫秒...     

On the harmonic balance with linearization for asymmetric single degree of freedom non-linear oscillators

This paper deals with analytical approximation of non-linear oscillations of conservative asymmetric single degree of freedom systems, using the method of harmonic balance with linearization. This technique which consists of linearizing the governing equations prior to harmonic balance permits us to avoid solving complicated non-linear algebraic equations. But it could be applied only to symmetric oscillations for which it proves to be very simple and effective. This restriction is due to the fact that the method requires an appropriate initial approximate solution as input. Such a solution could not be readily identified for nonsymmetric oscillations, contrary the symmetric case where the fundamental harmonic works well. For these nonsymmetric oscillations, we propose in this paper to consider an initial approximation which consists of a small bias plus the fundamental harmonic. By expanding the corresponding harmonic balance equations respectively to first and second order in the bias, we are able to easily determine the bias and thus the required initial approximate solution that yields consistent solution at higher order. We use three examples to illustrate the proposed approach and reveal its simplicity and its very good convergence.     

Implicit solution of harmonic balance equation system using the LU-SGS method and one-step Jacobi/Gauss-Seidel iteration

ABSTRACT

This paper presents efficient alternative numerical methods for an implicit solution of the harmonic balance equation system for analysing temporal periodic unsteady flows. The proposed method employs approximate factorisation to decouple the common residual term and the time spectral source term of a harmonic balance equation system when it is discretised implicitly. With this approximate factorisation, the complexity of implicit solution of the discrete system is greatly reduced. The common residual term can be dealt with using a lower-upper symmetric-Gauss-Seidel (LU-SGS) method and the time spectral source term is integrated using a Jacobi iteration (JI) or one step Gauss-Seidel (GS) iteration, leading to the LU-SGS/JI method or LU-SGS/GS method. The NASA stage 35 compressor and the 1.5 stage Aachen turbine were used to demonstrate the effectiveness of the proposed methods in stabilising solution and its advantages in comparison with the existing lower-upper symmetric-Gauss-Seidel/block Jacobi (LU-SGS/BJ) method. The LU-SGS/GS method and the LU-SGS/JI method are more robust than the LU-SGS/BJ method in stabilising solution. The LU-SGS/GS method also has faster and tighter convergence and lower memory consumption in comparison with the LU-SGS/BJ method.