No abstract

SUMMARYThe overset mesh method chimera is popular within the rotorcraft research community, because the use of multiple, non‐matching grids make the CFD simulations of bodies in relative motion much simpler. Consequently, the relative motion between the helicopter blades and fuselage can be accurately accounted for. In this paper, the method for treating overset grids within CFD codes is presented. It is compatible with multi‐block, structured‐grid solvers. The proposed method is based on hierarchy of overset, non‐matching grids, whose cells are automatically identified as computational or non‐computational and localised with respect to all grids they overlap with. The efficiency of the method relies on the hierarchical, multi‐step approach, for the overset mesh localisation and the use of a tree search. Because of the high efficiency of the algorithm, the search for overlapping cells can be carried out on‐the‐fly, during time‐marching of the unsteady, implicit CFD solver. In addition, the algorithm is suitable for parallel execution. The method has been demonstrated for several flows, ranging from simple aerofoils to rotor‐body interaction. The paper presents and demonstrates the method and shows that it has a low CPU overhead. It also highlights the limitations of the method and suggests remedies for improvement. Copyright © 2013 John Wiley & Sons, Ltd.

The computation of the aerodynamic forces arising from forced periodic motions is required for the generation of dynamic terms in models for flight simulation. The periodicity can be used to avoid using fully unsteady calculations by using the harmonic balance method. The current paper develops an implicit solver for the harmonic balance equations. The method is tested on two transonic test cases and evaluation is made against the unsteady simulation results. The first case is for the pitching NACA 0012 aerofoil. The second is for forced pitching of the F-5 wing with a wing tip launcher and missile. A reduction in computational time by one order of magnitude compared with the unsteady solver is obtained. Nomenclature A = matrix in frequency domain equation c = chord D = matrix in harmonic balance equation E = transformation matrix between frequency and time domains e = energy F, G, H = convective fluxes I = residual of semidiscrete system I = identity matrix k = reduced frequency n H = number of harmonics p = pressure R = residual vector T = period t = time u, v, w = Cartesian velocity components W = conserved variables = angle of attack t = pseudo time step t = real time step = density ! = frequency Subscripts m = mean value hb = harmonic balance 1 = freestream value 0 = amplitude Superscripts n = time level = Fourier coefficient

No abstract

No abstract