Search

A term-by-term comparison has been conducted between the residuals of the continuous and discrete adjoint methods for quasi-1D flow equations. This comparison is done to identify the origin of numerical instabilities near boundaries of the discrete adjoint method. The results show that there is one-to-one correspondence between the terms of the two residuals. Furthermore, the adjoint boundary conditions are rigorously analyzed. It turns out that improvement can be achieved by replacing some of the terms of the discrete adjoint residual with the counterparts in the continuous adjoint method.

Unsteady turbulent flow simulations were performed based on the Reynolds-averaged Navier–Stokes (RANS) equations to investigate flow oscillation due to three-dimensional (3D) configuration of a Mach 1.5 supersonic open cavity flow with a length-to-depth ratio of 3. Two-dimensional (2D) and 3D unsteady simulation results were analyzed and compared with experimental data and Rossiter’s empirical prediction data. The three-dimensional cavity width-to-depth ratio (W/D) was 1, 3.8 and 7.6. Computational results indicated that pressure oscillation in the 2D flow was generated by a single-flow structure, whereas a multiple-flow structure generated multiple oscillation peaks in the 3D flow. The flow structure in the 3D cavity was investigated. For the 2D flow case, the cavity internal pressure wave was directly synchronized with the free shear layer. In the 3D flow case, an unstable spanwise flow due to the sidewall was observed. This spanwise fluctuation produced additional pressure oscillations coupled with the streamwise internal pressure wave. The numerical results indicate that the spanwise flow reduces the propagation speed of the internal pressure waves and the intensity of the corresponding pressure fluctuation.

In this paper, we present a computational aeroelastic analysis of flexible flapping wings in the vicinity of solid walls. The wall effects change the aerodynamic forces and moments of the wings, and thus the aeroelastic behavior. The numerical simulation is carried out using a fluid–structure interaction framework by coupling the computational fluid dynamics and computational structural dynamics. A preconditioned Navier–Stokes solver based on a finite volume method is used for the aerodynamic analysis. The structural analysis is performed using a nonlinear structural model based on a geometrically exact beam formulation. The method is validated using previous numerical and experimental results. The aeroelastic characteristics of the flexible wings with and without the walls are computed and compared.

The RANS equations are widely used to analyze complex flows over aircraft. The equations require a turbulence model for turbulent flow analyses. A suitable turbulence must be selected for accurate predictions of aircraft aerodynamic characteristics. In this study, numerical analyses of three-dimensional aircraft are performed to compare the results of various turbulence models for the prediction of aircraft aerodynamic characteristics. A 3-D RANS solver, MSAPv, is used for the aerodynamic analysis. The four turbulence models compared are the Sparlart–Allmaras (SA) model, Coakley’s $$q-\omega $$ model, Huang and Coakley’s $$k-\varepsilon $$ model, and Menter’s $$k-\omega $$ SST model. Four aircrafts are considered: an ARA-M100, DLR-F6 wing–body, DLR-F6 wing–body–nacelle–pylon from the second drag prediction workshop, and a high wing aircraft with nacelles. The CFD results are compared with experimental data and other published computational results. The details of separation patterns, shock positions, and $$C_{p}$$ distributions are discussed to find the characteristics of the turbulence models.

We propose a numerical scheme to simulate the time-domain echo signals at tracking radar for a realistic scenario where an EAD (expendable active decoy) and an airborne target are both in dynamic states. On various scenarios where the target takes different maneuvers, the trajectories of the EAD ejected from the target are accurately calculated by solving 6-DOF (Degreeof-Freedom) equations of the motion for the EAD. At each sampling time of the echo signal, the locations of the EAD and the target are assumed to be fixed. Thus, the echo power from the EAD can be simply calculated by using the Friis transmission formula. The returned power from the target can be computed based on the pre-calculated scattering matrix of the target. In this paper, an IPO (iterative physical optics) method is used to construct the scattering matrix database of the target. The sinc function-interpolation formulation (sampling theorem) is applied to compute the scattering at any incidence angle from the database. A simulator is developed based on the proposed scheme to estimate the echo signals, which can consider the movement of the airborne target and EAD, also the scattering of the target and the RF specifications of the EAD. For applications, we consider the detection probability of the target in the presence of the EAD based on Monte Carlo simulation.

A general purpose aircraft mission performance analysis program has been developed. The program can be used in design mode or in analysis mode. Fuel weight for a given mission profile can be estimated when the design mode is chosen, while mission time or mission range for a given fuel can be estimated when the analysis mode is chosen. The mission analysis program is written with JAVA and includes GUI(Graphic User Interface) for users’ conveniences. With a proper combination of databases for propulsion, aerodynamics and weight, the program can be configured to compute the performance of any type of aircraft. The program is validated by comparing its results with the results of a well known performance analysis program by ADD(Agency for Defense Development)

This paper investigates the convergence characteristics of the modified artificial compressibility method proposed by Turkel. In particular, a focus is mode on the convergence characteristics due to variation of the preconditioning factor (αu) and the artificial compressibility ( β) in conjunction with an upwind method. For the investigations, a code using the modified artificial compressibility is developed. The code solves the axisymmetric incompressible Reynolds averaged Navier-Stokes equations. The cell-centered finite volume method is used in conjunction with Roe’s approximate Riemann solver for the inviscid flux, and the central difference discretization is used for the viscous flux. Time marching is accomplished by the approximated factorizationalternate direction implicit method. In addition, Menter’s k-ω shear stress transport turbulence model is adopted for analysis of turbulent flows. Inviscid, laminar, and turbulent flows are solved to investigate the accuracy of solutions and convergence behavior in the modified artificial compressibility method. The possible reason for loss of robustness of the modified artificial compressibility method with αu >1.0 is given.

A program package for store separation analysis, SAFESEP, has been developed. The package includes CTS system software, 3-Dimensional visualization software, off-line store trajectory software and aerodynamic database software. In this paper, we describe the functions and the capabilities of the store separation analysis program package.