Search

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 program that predicts trajectories of projectiles influenced by the interference flow field of helicopters is developed. The interference flow field are computed using a compressible inviscid solver in conjunction with an actuator disc model. The trajectories are predicted using 6-DOF (Degree of Freedom) equations as well as an alternative form of modified point mass equations of motion. The method for the interference flow field prediction method are validated with ROBIN(ROtor Body INteraction) model. A Sierra international bullet and a 105mm projectile are used to validate the trajectory method. Trajectories of a Sierra International bullet and a HYDRA 70 rocket firing from a helicopter are predicted.

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)

In this paper, numerical simulations are performed on the lock-in phenomena of vortex induced vibration(VIV) of a two dimensional cylinder. A deforming grid as well as a rigidly moving grid are used to simulate the movement of the cylinder. The grid deformation is accomplished by the linear spring analogy. Converged solutions, which are obtained by controling the grid size and the non-dimensional time step, are used for comparison and validation of the analysis results. Moreover, the efficiency and the accuracy of the coupling methods for fluid-structure interaction are examined.

This paper presents the computations of the longitudinal and directional stability derivatives for the SDM(Standard Dynamic Model). The static and dynamic derivatives are evaluated at once using forced harmonic oscillations in the pitch and yaw directions. For the numerical simulations, a 3-D Euler solver that uses a dual time stepping method for unsteady time accurate simulations is applied. This work investigates the variation of the derivatives in terms of the Mach number and the several motion parameters. Good agreement of the pitch and yaw stability derivatives with previously published numerical results and experimental results are observed.

A computer program that can estimate static, dynamic stability and control derivatives using a subsonic-supersonic panel method is developed. The panel method uses subsonic-supersonic source and elementary horse shoe vortex distributions, and their strengths are determined by solving the boundary condition approximated with a thin body assumption. In addition, quasi-steady analysis on the body fixed coordinate system allows the estimation of damping coefficients of aircraft 3 axes. The code is validated by comparing the neutral point, roll and pitch damping of delta wings with published analysis results. Finally, the static, dynamic stability and control derivatives of F-18 are compared with experimental data as well as other numerical results to show the accuracy and the usefulness of the code. 초 록 아음속-초음속 패널법(panel method)을 이용하여 항공기의 정적 안정성 미계수와 동적 안정성 미계수 및 조종미계수를 예측할 수 있는 프로그램을 개발하였다. 사용된 코드는 아음속-초음속 소스(source)와 말굽 와류(elementary horse shoe vortex)의 분포를 사용하고, 그 분포의 크기는 얇은 물체 근사(thin body approximation)을 적용하여 간략히 한 경계 조건을 이용하여 계산하였다. 항공기에 부착된 물체 좌표계에서 준정상(quasi-steady) 해석 을 통해서 항공기 3축의 댐핑 계수를 예측하였다. 개발된 코드는 삼각날개(delta-wing)의 중립점(neutral point), 롤, 피치 댐핑 계수의 이론치와 비교하여 검증하였다. 마지막으로 F-18의 정적, 동적 안정성 미계수와 조종 미계수를 풍동 시험치와 계산치에 비교하여 개발 한 코드의 정확성과 유용성을 확인하였다.

This paper presents computations of the dynamic derivatives of three dimensional flight vehicle configurations using CFD. The pitch dynamic derivatives are computed from the pitch sinusoidal motion, while the roll damping is computed based on steady state calculation using a non-inertial frame method. The Basic Finner and the SDM(Standard Dynamic Model) are chosen for the benchmark tests against other numerical and experimental results. For the flow calculations, a 3-D Euler solver that can be run both on the non-inertial frame and on the inertial frame is developed. A dual-time stepping method is applied for the unsteady time accurate simulations. A good agreement of pitch-roll dynamic derivatives with previously published numerical results and the experimental results is observed.

In this study, Constraint Force Equation Methodology (CFEM) is used to develop a multi-body dynamic analysis program with constraints. Seven constraint models are implemented to analyze constraint motions of multiple bodies. The augmented equations with the constraints are solved with the 4th order Runge-Kutta method for higher degree of accuracy. The analysis code is verified by comparing the analysis results of the motion of bodies with various constraints to published results.

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.

In this study, fairing configurations for an aircraft are designed and the aerodynamic analyses of the fairings are performed to find the best choice for the aircraft. Fairings considered are wing-tip fairing and wing-body fairing. Wing alone analyses are done for the wing-tip faring selection, while wing-body-tail analyses are done for the wing-body fairing selection. A 3-D RANS solver with Menter`s SST turbulence model are used for the aerodynamic analyses. The effects on the drag of the aircraft are examined by comparing the analysis results with and without the farings.

Off-line 6-DOF simulation program for store separation analysis has been developed. The developed program enables to predict a trajectory of a store from the database which was constructed by wind tunnel testing or CFD analysis. The flow angle method was applied to the program for predicting aerodynamic coefficients from the database and the ejector forces and constraints were enabled to incorporate the equations of motion for computing the trajectory. Using the program, the trajectories were calculated and the results are compared with the CTS results.

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.

Euler codes or Navier-Stokes codes for compressible flows suffer severe degradation in convergence as Mach number approaches zero. The convergence problem arose from the wide disparity in characteristic speeds can be solved using preconditioning methods without large modifications. In this paper, a preconditioned RANS(Reynolds Averaged Navier-Stokes) solver is developed for analysis of low Mach number flows. In order to validate the method, computational examples are chosen and the results are compared with the experimental data and the existing computed results showing a good accuracy and convergence characteristics for steady inviscid, laminar and turbulent flows at low Mach number.

Aerodynamic data measured in a wind tunnel has inevitable errors due to the presence of the wind tunnel walls. These unwanted interference effects must be corrected for the wall interference free aerodynamic data. Streamline curvature effects are caused by straightening of streamlines due to wind tunnel walls. Classical Glauert`s correction method that is a standard method for fixed wing aircraft is not suitable for rotary wing aircraft. In this paper, Heyson`s correction method of which wake model is compatible with rotors is used to correct the rotor shaft angle as well as the dynamic pressure. The results of Heyson`s method are compared with Glauert`s correction method.

6-DOF simulation program is developed in order to increase the efficiency of the store separation analysis. This S/W is much faster than a method based on CFD(Computational Fluid Dynamics) technology, and allows the simulation of stores with fixed shape as well as with extensible wings, because it uses aerodynamic databases which are prepared beforehand. In this paper, aerodynamic databases of stores are obtained with MSAP(Multi-body Separation Analysis Program), and unsteady damping coefficients are modeled with Missile Datcom. These databases and the 6-DOF simulation program are used to predict the trajectory of an external store, while its wings are being deployed. The analysis results indicate that the safe separations of the store can be achieved not only with the wing fixed but with the wings being deployed.

In this paper, we summarize the results of free drop wind tunnel test, separation analysis and flight test in order to verify the safety during the separations of an external fuel tank and the LAU-131 rocket launcher from XKO-l. The wind tunnel test was conducted to show the safety in free drop of the stores and to gather the trajectory data for fine tune of MSAP(Multi-body Separation Analysis Program). The enhanced MSAP was then used to predict the trajectories of the stores with and without the ejector forces. A correlation of MSAP results for free drop case was also made to show the safety of jettison with the free drop type bomb rack. Moreover, the flight test was conducted. and its results were compared to analysis results. Finally, the safe jettison boundary was determined from the flight test.

A study on the test, design and fabrication of wind tunnel model for measurement of air load distribution on wing surfaces is presented. 447 pressure taps are installed normal to the wing surfaces, and measured by PSI-8400 system using total 8 ESPs modules installed in the model. The test was performed at 50 m/sec constant speed in the low speed wind tunnel of Agency for Defense Development. Tests were carried out to determine effects of angle of attack, angle of sideslip and flap and stores for the load distribution of wing. The test results in this paper can be applied to the design optimization of structure and validation of computational fluid dynamics.