Based on modified Leishman-Beddoes （L-B） state space model at low Mach number （lower than 0.3）, the airfoil aeroelastic system is presented in this paper. The main modifications for L-B model include a new dynamic...Based on modified Leishman-Beddoes （L-B） state space model at low Mach number （lower than 0.3）, the airfoil aeroelastic system is presented in this paper. The main modifications for L-B model include a new dynamic stall criterion and revisions of normal force and pitching moment coefficient. The bifurcation diagrams, the limit cycle oscillation （LCO） phase plane plots and the time domain response figures are applied to investigating the stall flutter bifurcation behavior of airfoil aeroelastic systems with symmetry or asymmetry. It is shown that the symmetric periodical oscillation happens after subcritical bifurcation caused by dynamic stall, and the asymmetric periodical oscillation, which is caused by the interaction of dynamic stall and static divergence, only happens in the airfoil aeroelastic system with asymmetry. Validations of the modified L-B model and the airfoil aeroelastic system are presented with the experimental airload data of NACA0012 and OA207 and experimental stall flutter data of NACA0012 respectively. Results demonstrate that the airfoil aeroelastic system presented in this paper is effective and accurate, which can be applied to the investigation of airfoil stall flutter at low Mach number.展开更多
A novel extension to SMAC scheme is proposed for variable density flows under low Mach number approximation. The algorithm is based on a predictor—corrector time integration scheme that employs a projection method fo...A novel extension to SMAC scheme is proposed for variable density flows under low Mach number approximation. The algorithm is based on a predictor—corrector time integration scheme that employs a projection method for the momentum equation. A constant-coefficient Poisson equation is solved for the pressure following both the predictor and corrector steps to satisfy the continuity equation at each time step. The proposed algorithm has second order centrally differenced convective fluxes with upwinding based on Cell Peclet number while diffusive flux are viscous fourth order accurate. Spatial discretization is performed on a collocated grid system that offers computational simplicity and straight forward extension to curvilinear coordinate systems. The algorithm is kinetic energy preserving. Further in this paper robustness and accuracy are demonstrated by performing test on channel flow with non-Boussinesq condition on different temperature ratios.展开更多
基金Aeronautical Science Foundation of China (08A52003)Science and Technology Foundation of Rotorcraft Aeromechanics Laboratory (9140C4001010901)
文摘Based on modified Leishman-Beddoes （L-B） state space model at low Mach number （lower than 0.3）, the airfoil aeroelastic system is presented in this paper. The main modifications for L-B model include a new dynamic stall criterion and revisions of normal force and pitching moment coefficient. The bifurcation diagrams, the limit cycle oscillation （LCO） phase plane plots and the time domain response figures are applied to investigating the stall flutter bifurcation behavior of airfoil aeroelastic systems with symmetry or asymmetry. It is shown that the symmetric periodical oscillation happens after subcritical bifurcation caused by dynamic stall, and the asymmetric periodical oscillation, which is caused by the interaction of dynamic stall and static divergence, only happens in the airfoil aeroelastic system with asymmetry. Validations of the modified L-B model and the airfoil aeroelastic system are presented with the experimental airload data of NACA0012 and OA207 and experimental stall flutter data of NACA0012 respectively. Results demonstrate that the airfoil aeroelastic system presented in this paper is effective and accurate, which can be applied to the investigation of airfoil stall flutter at low Mach number.
文摘A novel extension to SMAC scheme is proposed for variable density flows under low Mach number approximation. The algorithm is based on a predictor—corrector time integration scheme that employs a projection method for the momentum equation. A constant-coefficient Poisson equation is solved for the pressure following both the predictor and corrector steps to satisfy the continuity equation at each time step. The proposed algorithm has second order centrally differenced convective fluxes with upwinding based on Cell Peclet number while diffusive flux are viscous fourth order accurate. Spatial discretization is performed on a collocated grid system that offers computational simplicity and straight forward extension to curvilinear coordinate systems. The algorithm is kinetic energy preserving. Further in this paper robustness and accuracy are demonstrated by performing test on channel flow with non-Boussinesq condition on different temperature ratios.
