Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

飞行器跨声速气动弹性数值分析

将流体和结构运动方程分别构造为含子迭代的计算格式，发展了一种紧耦合气动弹性分析方法．其中流体计算的空间离散采用改进的HLLEW（Harten-Lax-van Leer-Einfeldt-Wada）格式．TFI（transfinite interpolation）方法用于生成随结构变形的自适应多块动网格．利用所发展的方法，对一翼-身-尾气动外形，数值预测了马赫数在0．3-1．3范围内的气动颤振边界．并详细研究了时间步长、子迭代步数、初始流场、耦合方法、疏密网格对颤振计算结果的影响．

Numerical Analyses of Transonic Fluid/Structure interaction in Aircraft Engineering

Through the coupling between aerodynamic and structural governing equations, a fully implicit multiblock aeroelastic solver was developed for transonic fluid/stricture interaction. The Navier-Stokes fluid equations are solved based on LU-SGS (lower-upper symmetric Gauss-Seidel) Time-marching subiteration scheme and HLLEW (Harten-Lax-van Leer-Einfeldt-Wada) spacing discretization scheme and the same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Transfinite interpolation (TFI) is used for the grid deformation of blocks neighboring the flexible surfaces. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between fluid and structure. The developed code was fort validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. In the subsonic and transonic range, the calculated flutter speeds and frequencies agree well with experimental data, however, in the supersonic range, the present calculation overpredicts the experimental flutter points similar to other computations. Then the flutter character of a complete aircraft configuration is analyzed through the calculation of the change of structural stiffness. Finally, the phenomenon of aileron buzz is simulated for the weakened model of a supersonic transport wing/body model at Mach numbers of 0.98 and l.05. The calculated unsteady flow shows, on the upper surface, the shock wave becomes stronger as the aileron deflects downward, and the flow behaves just contrary on the lower surface of the wing. Corresponding to general theoretical analysis, the flow instability referred to as aileron buzz is induced by a stronger shock alternately moving on the upper and lower surfaces of wing. For the rigid structural model, the flow is stable at all calculated Mach numbers as observed in experiment

Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

基于并行CFD方法的颤振特性计算

将流体和结构运动方程分别构造为含子迭代的计算格式,发展一种紧耦合气动弹性分析方法。其中流体计算的空间离散采用改进的HLLEW(Harten-Lax-Van Leer-Einfeldt-Wada)格式,同时采用MPI通信实现并行计算。弹性体网格变形法用于生成随结构变形的四面体非结构网格。利用所发展的方法,对一标模445.6翼型进行了跨声速气动颤振特性研究,分析其颤振特性并对比实验结果。同时还运用所发展的方法对某机带控制舵垂尾的颤振特性进行了研究,确立了所发展方法对实际工程颤振问题的适应性。

System size dependence of associated yields in hadron-triggered jets STAR Collaboration

We present results on the system size dependence of high transverse momentum di-hadron correlations at root s(NN) = 200 GeV as measured by STAR at RHIC. Measurements in d + Au, Cu + Cu and Au + Au collisions reveal similar jet-like near-side correlation yields (correlations at small angular separation Delta phi similar to 0, Delta eta similar to 0) for all systems and centralities. Previous measurements have shown Chat the away-side (Delta phi similar to pi) yield is suppressed in heavy-ion collisions. We present measurements of the away-side Suppression as a function of transverse momentum and centrality in Cu + Cu and Au + Au collisions. The suppression is found to be similar in Cu + Cu and An + An collisions at a similar number of participants. The results are compared to theoretical calculations based on the patron quenching model and the modified fragmentation model. The observed differences between data and theory indicate that the correlated yields presented here will further constrain dynamic energy loss models and provide information about the dynamic density profile in heavy-ion collisions. (C) 2009 Elsevier B.V. All rights reserved.

Measurement of the Bottom Quark Contribution to Nonphotonic Electron Production in p plus p Collisions at root s=200 GeV

The contribution of B meson decays to nonphotonic electrons, which are mainly produced by the semileptonic decays of heavy-flavor mesons, in p + p collisions at root s = 200 GeV has been measured using azimuthal correlations between nonphotonic electrons and hadrons. The extracted B decay contribution is approximately 50% at a transverse momentum of p(T) >= 5 GeV/c. These measurements constrain the nuclear modification factor for electrons from B and D meson decays. The result indicates that B meson production in heavy ion collisions is also suppressed at high p(T).

Higher Moments of Net Proton Multiplicity Distributions at RHIC

200 GeV corresponding to baryon chemical potentials (mu(B)) between 200 and 20 MeV. Our measurements of the products kappa sigma(2) and S sigma, which can be related to theoretical calculations sensitive to baryon number susceptibilities and long-range correlations, are constant as functions of collision centrality. We compare these products with results from lattice QCD and various models without a critical point and study the root s(NN) dependence of kappa sigma(2). From the measurements at the three beam energies, we find no evidence for a critical point in the QCD phase diagram for mu(B) below 200 MeV.

Three-Particle Coincidence of the Long Range Pseudorapidity Correlation in High Energy Nucleus-Nucleus Collisions

We report the first three-particle coincidence measurement in pseudorapidity (Delta eta) between a high transverse momentum (p(perpendicular to)) trigger particle and two lower p(perpendicular to) associated particles within azimuth |Delta phi| < 0.7 in root s(NN) = 200 GeV d + Au and Au + Au collisions. Charge ordering properties are exploited to separate the jetlike component and the ridge (long range Delta eta correlation). The results indicate that the correlation of ridge particles are uniform not only with respect to the trigger particle but also between themselves event by event in our measured Delta eta. In addition, the production of the ridge appears to be uncorrelated to the presence of the narrow jetlike component.

gamma cross section in p plus p collisions at root s=200 GeV

We report on a measurement of the gamma(1S + 2S + 3S) -> e(+)e(-) cross section at midrapidity in p + p collisions at root s = 200 GeV. We find the cross section to be 114 +/- 38(stat + fit)(-24)(+23)(syst) pb. Perturbative QCD calculations at next-to-leading order in the color evaporation model are in agreement with our measurement, while calculations in the color singlet model underestimate it by 2 sigma. Our result is consistent with the trend seen in world data as a function of the center-of-mass energy of the collision and extends the availability of gamma data to RHIC energies. The dielectron continuum in the invariant-mass range near the gamma is also studied to obtain a combined yield of e(+)e(-) pairs from the sum of the Drell-Yan process and b-(b) over bar production.

Power law behavior of the isotope yield distributions in the multifragmentation regime of heavy ion reactions

Isotope yield distributions in the multifragmentation regime were studied with high-quality isotope identification, focusing on the intermediate mass fragments (IMFs) produced in semiviolent collisions. The yields were analyzed within the framework of a modified Fisher model. Using the ratio of the mass-dependent symmetry energy coefficient relative to the temperature, a(sym)/T, extracted in previous work and that of the pairing term, a(p)/T, extracted from this work, and assuming that both reflect secondary decay processes, the experimentally observed isotope yields were corrected for these effects. For a given I = N - Z value, the corrected yields of isotopes relative to the yield of C-12 show a power law distribution Y (N, Z)/Y(C-12) similar to A(-tau) in the mass range 1 <= A <= 30, and the distributions are almost identical for the different reactions studied. The observed power law distributions change systematically when I of the isotopes changes and the extracted tau value decreases from 3.9 to 1.0 as I increases from -1 to 3. These observations are well reproduced by a simple deexcitation model, with which the power law distribution of the primary isotopes is determined to be tau(prim) = 2.4 +/- 0.2, suggesting that the disassembling system at the time of the fragment formation is indeed at, or very near, the critical point.

Parton energy loss in heavy-ion collisions via direct-photon and charged-particle azimuthal correlations

Charged-particle spectra associated with direct photon (gamma(dir)) and pi(0) are measured in p + p and Au + Au collisions at center-of-mass energy root(S)(NN) = 200 GeV with the STAR detector at the Relativistic Heavy Ion Collider. A shower-shape analysis is used to partially discriminate between gamma(dir) and pi(0). Assuming no associated charged particles in the gamma(dir) direction ( near side) and small contribution from fragmentation photons (gamma(frag)), the associated charged-particle yields opposite to gamma(dir) (away side) are extracted. In central Au + Au collisions, the charged-particle yields at midrapidity (vertical bar eta vertical bar < 1) and high transverse momentum (3 < (assoc)(PT) < 16 GeV/c) associated with gamma(dir) and pi(0) (vertical bar eta vertical bar < 0.9, 8 < (trig)(PT) < 16 GeV/c) are suppressed by a factor of 3-5 compared with p + p collisions. The observed suppression of the associated charged particles is similar for gamma(dir) and pi(0) and independent of the gamma(dir) energy within uncertainties. These measurements indicate that, in the kinematic range covered and within our current experimental uncertainties, the parton energy loss shows no sensitivity to the parton initial energy, path length, or color charge.

Azimuthal di-hadron correlations in d plus Au and Au plus Au collisions at root s(NN)=200 GeV measured at the STAR detector

Yields, correlation shapes, and mean transverse momenta p(T) of charged particles associated with intermediate-to high-p(T) trigger particles (2.5 < p(T) < 10 GeV/c) in d + Au and Au + Au collisions at root s(NN) = 200 GeV are presented. For associated particles at higher p(T) greater than or similar to 2.5 GeV/c, narrow correlation peaks are seen in d + Au and Au + Au, indicating that the main production mechanism is jet fragmentation. At lower associated particle pT < 2 GeV/c, a large enhancement of the near- (Delta phi similar to 0) and away-side (Delta phi similar to pi) associated yields is found, together with a strong broadening of the away-side azimuthal distributions in Au + Au collisions compared to d + Au measurements, suggesting that other particle production mechanisms play a role. This is further supported by the observed significant softening of the away-side associated particle yield distribution at Delta phi similar to pi in central Au + Au collisions.

Balance functions from Au+Au, d+Au, and p+p collisions at root s(NN)=200 GeV

Balance functions have been measured for charged-particle pairs, identified charged-pion pairs, and identified charged-kaon pairs in Au + Au, d + Au, and p + p collisions at root s(NN) = 200 GeV at the Relativistic Heavy Ion Collider using the STAR detector. These balance functions are presented in terms of relative pseudorapidity, Delta eta, relative rapidity, Delta y, relative azimuthal angle, Delta phi, and invariant relative momentum, q(inv). For charged-particle pairs, the width of the balance function in terms of Delta eta scales smoothly with the number of participating nucleons, while HIJING and UrQMD model calculations show no dependence on centrality or system size. For charged-particle and charged-pion pairs, the balance functions widths in terms of Delta eta and Delta y are narrower in central Au + Au collisions than in peripheral collisions. The width for central collisions is consistent with thermal blast-wave models where the balancing charges are highly correlated in coordinate space at breakup. This strong correlation might be explained by either delayed hadronization or limited diffusion during the reaction. Furthermore, the narrowing trend is consistent with the lower kinetic temperatures inherent to more central collisions. In contrast, the width of the balance function for charged-kaon pairs in terms of Delta y shows little centrality dependence, which may signal a different production mechanism for kaons. The widths of the balance functions for charged pions and kaons in terms of q(inv) narrow in central collisions compared to peripheral collisions, which may be driven by the change in the kinetic temperature.