• Proceeding of EDISON Challenge
• /
• 2016.11a
• /
• pp.49-54
• /
• 2016

본 연구에서는 유/무인 우주발사체의 재진입모듈에 적용되는 Heat-Shield(열차단막)의 형상에 따른 유동특성을 수치해석으로 분석했다. 재진입모듈이 지구의 대기권을 다시 진입하는 환경(고도 70km)을 기준으로 해석을 진행했다. 열차단막의 형상은 평판, 곡률이 다른 타원으로 나누었으며, 각 형상별 유동특성을 확인했다. 결론적으로 재진입 모듈의 형상에 따라 압력과 속도분포를 기준으로 열적인 분포를 예상했으며, 계산된 항력계수를 비교했다. 단순한 유동으로 열적 분포를 예상한 것에 한계가 있지만 대기권 재진입 모듈의 2차원 설계에 도움이 될 형상기준을 제시했다. 수치해석은 모두 Edison_CFD에서 제공하는 툴을 이용해 수행했으며, 전처리에는 e-mega (structured)을 후처리에는 e-dava를 이용했다. 해석 solver는 '정렬격자기반 2차원 압축성 유동 범용해석 소프트웨어'를 이용했다.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.7
• /
• pp.51-60
• /
• 2002

The Monte-Carlo simulation of statistical analysis is used to investigate the final conditions of states as well as the footprint boundaries resulting from the atmospheric re-entry dispersions. The re-entry dispersions in this paper are specified by a $7\times7$ covariance matrix of latitude, longitude, altitude, bank angle, flight path angle, heading error, and range at entry velocity. The error sources that affect these at re-entry for a deboost are the uncertainties associated with atmospheric density and temperature, initial errors, wind, and estimation error of aerodynamic coefficients. Using $3{\sigma}_n$ deviations of these errors and a nominal flight trajectory, the covariance matrix of state variables can be determined by performing a trajectory error analysis. Major considerations in the application of the Monte-Carlo method are the simulation of perturbed trajectories, bank reversal, and determination of the impact points for each of these trajectories. This paper analyzes the results of uncertainties from the viewpoint of aero-coefficients and bank reversal.

• Journal of Power System Engineering
• /
• v.7 no.4
• /
• pp.67-73
• /
• 2003

The entry guidance design involves trajectory optimization and generation of a drag acceleration profile as the satisfaction of trajectory conditions during the entry flight. The reference trajectory is parameterized and optimized as piecewise linear functions of the velocity. A regularization technique is employed to achieve desired properties of the optimal drag profile. The regularized problem has smoothness properties and the minimization of performance index then prevents the drag acceleration from varying too fast, thus eliminating discontinuities. This paper shows the trajectory control using the simple control law as well as the information of reference drag acceleration.

• Journal of the Korean Society for Precision Engineering
• /
• v.17 no.6
• /
• pp.111-118
• /
• 2000

The design of reference trajectory with respect to drag acceleration is necessary to decelerate from hypersonic speed safely after atmosphere re-entry of space vehicle. The re-entry guidance design involves trajectory optimization, generation of a reference drag acceleration profile with the satisfaction of 6 trajectory constraints during the re-entry flight. This reference drag acceleration profile can be considered as the reference trajectory. The cost function is composed of the accumulated total heating on vehicle due to the reduction of weight. And a regularization is needed to prevent optimal drag profile from varying too fast and achieve realized trajectory. This paper shows the relations between velocity, drag acceleration and altitude in drag acceleration profile, and how to determine the reference trajectory.

• Journal of Astronomy and Space Sciences
• /
• v.22 no.2
• /
• pp.147-174
• /
• 2005

To understand the physical processes that control the high-latitude lower thermospheric dynamics, we quantify the forces that are mainly responsible for maintaining the high-latitude lower thermospheric wind system with the aid of the National Center for Atmospheric Research Thermosphere-Ionosphere Electrodynamics General Circulation Model (NCAR-TIEGCM). Momentum forcing is statistically analyzed in magnetic coordinates, and its behavior with respect to the magnitude and orientation of the interplanetary magnetic field (IMF) is further examined. By subtracting the values with zero IMF from those with non-zero IMF, we obtained the difference winds and forces in the high-latitude 1ower thermosphere(<180 km). They show a simple structure over the polar cap and auroral regions for positive($B_y$ > 0.8|$\overline{B}_z$ |) or negative($B_y$ < -0.8|$\overline{B}_z$|) IMF-$\overline{B}_y$ conditions, with maximum values appearing around -80$^{\circ}$ magnetic latitude. Difference winds and difference forces for negative and positive $\overline{B}_y$ have an opposite sign and similar strength each other. For positive($B_z$ > 0.3125|$\overline{B}_y$|) or negative($B_z$ < -0.3125|$\overline{B}_y$|) IMF-$\overline{B}_z$ conditions the difference winds and difference forces are noted to subauroral latitudes. Difference winds and difference forces for negative $\overline{B}_z$ have an opposite sign to positive $\overline{B}_z$ condition. Those for negative $\overline{B}_z$ are stronger than those for positive indicating that negative $\overline{B}_z$ has a stronger effect on the winds and momentum forces than does positive $\overline{B}_z$ At higher altitudes(>125 km) the primary forces that determine the variations of tile neutral winds are the pressure gradient, Coriolis and rotational Pedersen ion drag forces; however, at various locations and times significant contributions can be made by the horizontal advection force. On the other hand, at lower altitudes(108-125 km) the pressure gradient, Coriolis and non-rotational Hall ion drag forces determine the variations of the neutral winds. At lower altitudes(<108 km) it tends to generate a geostrophic motion with the balance between the pressure gradient and Coriolis forces. The northward component of IMF By-dependent average momentum forces act more significantly on the neutral motion except for the ion drag. At lower altitudes(108-425 km) for negative IMF-$\overline{B}_y$ condition the ion drag force tends to generate a warm clockwise circulation with downward vertical motion associated with the adiabatic compress heating in the polar cap region. For positive IMF-$\overline{B}_y$ condition it tends to generate a cold anticlockwise circulation with upward vertical motion associated with the adiabatic expansion cooling in the polar cap region. For negative IMF-$\overline{B}_z$ the ion drag force tends to generate a cold anticlockwise circulation with upward vertical motion in the dawn sector. For positive IMF-$\overline{B}_z$ it tends to generate a warm clockwise circulation with downward vertical motion in the dawn sector.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.2
• /
• pp.141-149
• /
• 2010

This paper contains impulsive maneuver which considers fuel consumption balance of chief satellite and deputy satellite in satellite formation flying. Thrust input is obtained by Lagrange' Multiplier method which is constructed by cost function with weight parameter of each satellite. Energy matching constraint is applied for boundedness of relative orbit, and theoretical solutions are verified by simulation results. Simulations are divided into two scenarios, with or without air-drag effect. This paper's results are expected to be used in real satellite formation flying, when fuel-balancing impulsive maneuver for relative orbit boundedness is needed.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.6
• /
• pp.530-536
• /
• 2010

In-flight icing is a critical technical issue for aircraft safety and, in particular, ice accretions on aircraft surfaces can drastically impair aerodynamic performances and control authority. In order to investigate icing effects on the aerodynamic characteristics of KC-100 aircraft, a state-of-the-art CFD code, FENSAP-ICE, was used. A main wing section and full configuration of KC-100 aircraft were considered for the icing analysis. Also, shapes of iced area were calculated for the design of anti-/de-icing devices. The iced areas around leading edge of main wing and horizontal tail wing were observed maximum 7.07% and 11.2% of the chord length of wing section, respectively. In case of wind shield, 16.7% of its area turned out to be covered by ice. The lift of KC-100 aircraft were decreased to 64.3%, while the drag was increased to 55.2%.

• Journal of Astronomy and Space Sciences
• /
• v.25 no.4
• /
• pp.415-424
• /
• 2008

Kwak et al. (2008) found that the mean neutral wind pattern in the high-latitude lower thermosphere is dominated by rotational flow than by divergent flow. As an extension of the our previous work (Kwak et al. 2008), we performed a term analysis of vorticity equation that describes the driving forces for the rotational component of the horizontal wind in order to determine key processes that causes strong rotational flow in the high-latitude lower thermospheric winds. For this study the National Center for Atmospheric Research Thermosphere-Ionosphere Electrodynamics General Circulation Model (NCAR-TIEGCM) is used. The primary forces that determine variations of the vorticity are the ion drag term and the horizontal advection term. Significant contributions, however, can be made by the stretching term. The effects of IMF on the vorticity forces are seen down to around 105-110km.

• Atmosphere
• /
• v.16 no.4
• /
• pp.303-318
• /
• 2006

A parameterization of gravity wave drag induced by cumulus convection (GWDC) proposed by Chun and Baik is implemented in the KMA operational global NWP model (GDAPS), and effects of the GWDC on the forecast for July 2005 by GDAPS are investigated. The forecast result is compared with NCEP final analyses data (FNL) and model's own analysis data. Cloud-top gravity wave stresses are concentrated in the tropical region, and the resultant forcing by the GWDC is strong in the tropical upper troposphere and lower stratosphere. Nevertheless, the effect of the GWDC is strong in the mid- to high latitudes of Southern Hemisphere and high latitudes of Northern Hemisphere. By examining the effect of the GWDC on the amplitude of the geopotential height perturbation with zonal wavenumbers 1-3, it is found that impact of the GWDC is extended to the high latitudes through the change of planetary wave activity, which is maximum in the winter hemisphere. The GWDC reduces the amplitude of zonal wavenumber 1 but increases wavenumber 2 in the winter hemisphere. This change alleviates model biases in the zonal wind not only in the lower stratosphere where the GWDC is imposed, but also in the whole troposphere, especially in the mid- to high latitudes of Southern Hemisphere. By examining root mean square error, it is found that the GWDC parameterization improves GDAPS forecast skill in the Southern Hemisphere before 7 days and partially in the Northern Hemisphere after about 5 days.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.4
• /
• pp.281-289
• /
• 2016

Accurate prediction of the trajectory and time of a time-varying mass parachute system remains essential in the mission requiring a precision airdrop to the ground. In this study, we investigate the altitude-varying behavior of a cross-type parachute system designed to deliver a time-varying mass object like flare. The dynamics of the descending parachute system was analyzed based on the Runge-Kutta method of the ordinary differential system. The drag coefficients of the cross-type parachute and flare were calculated by a CFD code based on the incompressible Navier-Stokes equation. Finally, by using a simplified gust wind model in troposphere, the combined effects of gust wind and time-varying mass were examined in detail.