• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.10
• /
• pp.954-961
• /
• 2008

An experimental study has been conducted to investigate the flow field around the wing having a circular damage hole. The damage was represented by a circular hole passing through the model with 10% airfoil chord diameter and normal to the chord. The hole was centered at quarter or half chord. The PIV flow fields and static pressure measurements on the wing upper and lower surface were carried out at Rec=2.85×105 based on the chord length. The PIV results showed the two types of flow structures around a damage hole were formed. The first one was a weak jet that formed an attached wake behind the damage hole. The second one resulted from increased incidence; this was a strong jet where the flow through the hole penetrates into the free-stream resulting in extensive separation of oncoming boundary layer flow and development of a separated wake with reverse flow. The surface pressure data showed a big pressure alteration near the circular damage hole. The severity of pressure alteration was increased as a damage hole located nearer to the leading edge.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.16 no.3
• /
• pp.178-189
• /
• 2004

For designing a reliable harbor, a methodology for estimating design waves of 97.5% operable harbor condition is suggested using long-term wave data. For a practical application of the methodology, a marine police harbor was selected as a site. Wave data used were collected from February 1993 to December 2003 at Jodo wave gage station in front of Pusan harbor. Joint distributions of significant wave height and significant wave period for specified wave directions were obtained and used to feed as input waves for parabolic mild-slope wave model. Results showed that input waves with significant wave height of 1.75 m, significant wave period off sec and wave direction E yield design waves height of 1.06 m at the site of interests, which is a 97.5% operable harbor condition. Wind waves generated inside harbor showed to be no effect on the design wave condition. Swells propagated from deep water into harbor are shown to be dominant effects on the design waves of operable harbor condition.

• Journal of Korean Society for Atmospheric Environment
• /
• v.32 no.5
• /
• pp.513-525
• /
• 2016

To propose an effective ensemble methods in predicting $PM_{10}$ concentration, six experiments were designed by different ensemble average methods (e.g., non-weighted, single weighted, and cluster weighted methods). The single weighted method was calculated the weighted value using both multiple regression analysis and singular value decomposition and the cluster weighted method was estimated the weighted value based on temperature, relative humidity, and wind component using multiple regression analysis. The effects of ensemble average methods were significantly better in weighted average than non-weight. The results of ensemble experiments using weighted average methods were distinguished according to methods calculating the weighted value. The single weighted average method using multiple regression analysis showed the highest accuracy for hourly $PM_{10}$ concentration, and the cluster weighted average method based on relative humidity showed the highest accuracy for daily mean $PM_{10}$ concentration. However, the result of ensemble spread analysis showed better reliability in the single weighted average method than the cluster weighted average method based on relative humidity. Thus, the single weighted average method was the most effective method in this study case.

• Journal of Astronomy and Space Sciences
• /
• v.21 no.4
• /
• pp.303-314
• /
• 2004

In this paper we report some properties of inner magnetospheric structure inferred from the T01_s code, one of the latest magnetospheric models by Tsyganenko. We have constructed three average storms representing moderate, strong, and severe intensity storms using 95 actual storms. The three storms are then modelled by the T01_s code to examine differences in magnetic structure among them. We find that the magnetic structure of intense storms is strikingly different from the normal structure. First, when the storm intensity is large, the field lines anchored at dayside longitudinal sectors become warped tailward to align to the solar wind direction. This is particularly so for the field lines anchored at the longitudinal sectors from postnoon through dusk. Also while for the moderate storm the equatorial magnetic field near geosynchronous altitude is found to be weakest near midnight sector, this depression region expands into even late afternoon sector during the severe storm. Accordingly the field line curvature radius at the equator in the premidnight geosynchronous region becomes unusually small, reaching down to a value less than 500 km. We attribute this strong depression and the dawn-dusk asymmetry to the combined effect from the enhanced tail current and the westward expansion/rotation of the partial ring current.

• Journal of Astronomy and Space Sciences
• /
• v.27 no.3
• /
• pp.205-212
• /
• 2010

To better understand the physical processes that maintain the high-latitude lower thermospheric dynamics, we have identified relative contributions of the momentum forcing and the heating to the high-latitude lower thermospheric winds depending on the interplanetary magnetic field (IMF) and altitude. For this study, we performed a term analysis of the potential vorticity equation for the high-latitude neutral wind field in the lower thermosphere during the southern summertime for different IMF conditions, with the aid of the National Center for Atmospheric Research Thermosphere-Ionosphere Electrodynamics General Circulation Model (NCAR-TIEGCM). Difference potential vorticity forcing and heating terms, obtained by subtracting values with zero IMF from those with non-zero IMF, are influenced by the IMF conditions. The difference forcing is more significant for strong IMF $B_y$ condition than for strong IMF $B_z$ condition. For negative or positive $B_y$ conditions, the difference forcings in the polar cap are larger by a factor of about 2 than those in the auroral region. The difference heating is the most significant for negative IMF $B_z$ condition, and the difference heatings in the auroral region are larger by a factor of about 1.5 than those in the polar cap region. The magnitudes of the difference forcing and heating decrease rapidly with descending altitudes. It is confirmed that the contribution of the forcing to the high-latitude lower thermospheric dynamics is stronger than the contribution of the heating to it. Especially, it is obvious that the contribution of the forcing to the dynamics is much larger in the polar cap region than in the auroral region and at higher altitude than at lower altitude. It is evident that when $B_z$ is negative condition the contribution of the forcing is the lowest and the contribution of the heating is the highest among the different IMF conditions.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.14 no.1
• /
• pp.81-88
• /
• 1996

Remote sensing technique is a probable means to estimate distribution of actual evapotranspiration over wide area in connection with regional characteristics of vegetation and landuse. Factors controlling evapotranspiration from ground are air temperature, humidity, wind, radiation, soil moisture and so on. Not only the vegetation influences directly the evapotranspiration, but also these factors strongly influnce the vegetation at the area. Therefore we can expect high correlation between the evapotranspiration and the vegetation. To grasp the state of vegetation at any point, NDVI calculated from NOAA/AVHRR data is utilized. It can be considered that evapotranspiration at a forest region is linearly proportional to the NDVI. Here, a model which adopts a direct method to estimate actual evapotranspiration is developed by using the relationship between NDVI and evapotranspiration. This method makes possible to estimate evapotranspiration of Korean Peninsula including North Korea where enough meteorological and hydrological data are unavailable.

• Atmosphere
• /
• v.23 no.1
• /
• pp.1-11
• /
• 2013

The angular momentum transport of an idealized axisymmetric vortex in the developing stage was investigated using the Weather Research and Forecast (WRF) model. The balanced axisymmetric vortex was constructed based on an empirical function for tangential wind, and the temperature, geopotential, and surface pressure were obtained from the balanced equation. The numerical simulation was carried out for 6 days on the f-plane with the Sea Surface Temperature (SST) set as constant. The weak vortex at initial time was intensified with time, and reached the strength of tropical cyclone in a couple of days. The Absolute Angular Momentum (AAM) was transported along with the secondary circulation of the vortex. Total AAM integrated over a cylinder of radius of 2000 km decreased with simulation time, but total kinetic energy increased rapidly. From the budget analysis, it was found that the surface friction is mainly responsible for the decrease of total AAM. Also, contribution of the surface friction to the AAM loss was about 90% while that of horizontal advection was as small as 8%. The trajectory of neutral numerical tracers following the secondary circulation was presented for the Lagrangian viewpoint of the transports of absolute angular momentum. From the analysis using the trajectory of tracers it was found that the air parcel was under the influence of the surface friction continuously until it leaves the boundary layer near the core. Then the air parcel with reduced amount of angular momentum compared to its original amount was transported from boundary layer to upper level of the vortex and contributed to form the anti-cyclone. These results suggest that the tropical cyclone loses angular momentum as it develops, which is due to the dissipation of angular momentum by the surface friction.

• Ocean Systems Engineering
• /
• v.3 no.4
• /
• pp.257-273
• /
• 2013

Top tensioned riser (TTR) is often used in a floating oil/gas production system deployed in deep water for oil/gas transport. This study focuses on the extension of the existing numerical code, known as CABLE3D, to allow for static and dynamic simulation of a TTR connected to a floating structure through a tensioner system or buoyancy can, and restrained by riser guides at different elevations. A tensioner system usually consists of three to six cylindrical tensioners. Although the stiffness of individual tensioner is assumed to be linear, the resultant stiffness of a tensioner system may be nonlinear. The vertical friction between a TTR and the hull at its riser guide is neglected assuming rollers are installed there. Near the water surface, a TTR is forced to move horizontally due to the motion of the upper deck of a floating structure as well as related riser guides. The extended CABLE3D is then integrated into a numerical code, known as COUPLE, for the simulation of the dynamic interaction among the hull of a floating structure, such as spar or TLP, its mooring system and riser system under the impact of wind, current and waves. To demonstrate the application of the extended CABLE3D and its integration with COUPLE, the numerical simulation is made for a truss spar under the impact of Hurricane "Ike". The mooring system of the spar consists of nine mooring lines and the riser system consists of six TTRs and two steel catenary risers (SCRs).

• Wind and Structures
• /
• v.30 no.5
• /
• pp.533-546
• /
• 2020

Large Eddy Simulation (LES) is used to study the effects of steady slot suction on the aerodynamic forces of and flow around a wall-mounted finite-length square cylinder. The aspect ratio H/d of the tested cylinder is 5, where H and d are the cylinder height and width, respectively. The Reynolds number based on free-stream oncoming flow velocity U_∞ and d is 2.78×10⁴. The suction slot locates near the leading edge of the free end, with a width of 0.025d and a length of 0.9d. The suction coefficient Q (= U_s/U_∞) is varied as Q = 0, 1 and 3, where U_s is the velocity at the entrance of the suction slot. It is found that the free-end steady slot suction can effectively suppress the aerodynamic forces of the model. The maximum reduction of aerodynamic forces occurs at Q = 1, with the time-mean drag, fluctuating drag, and fluctuating lift reduced by 3.75%, 19.08%, 40.91%, respectively. For Q = 3, all aerodynamic forces are still smaller than those for Q = 0 (uncontrolled case), but obviously higher than those for Q = 1. The involved control mechanism is successfully revealed, based on the comparison of the flow around cylinder free end and the near wake for the three tested Q values.

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

A model of advection-dominated accretion flows has been highlighted in the last decade. Most of calculations are based on self-similar solutions of equations governing the accreting flows. We revisit self-similar solutions of the simplest form of advection-dominated accretion flows. We explore the parameter space thoroughly and seek another category of self-similar solutions. In this study we allow the parameter f less than zero, which denotes the fraction of energy transported through advection. We have found followings: 1. For f > 0, in real ADAF solutions the ratio of specific heats ${\gamma}$ satisfies 1 < ${\gamma}$ < 5/3 for O ${\leq}$ f ${\leq}$ 1. On the other hands, in wind solutions a rotating disk does not exist. 2. For f < 0, in real ADAF solutions with ${\epsilon}$ greater than zero ${\gamma}$ requires rather exotic range, that is, ${\gamma}$ < 1 or ${\gamma}$ > 5/3. When -5/2 < ${\epsilon}$' < 0, however, allowable ${\gamma}$ can be found in ${\gamma}$ < 5/3, in which case 4{\Omega}_0$,_ is imaginary. 3. For a negative $u_0$,+ with f > 0, solutions are only allowed with exotic ${\gamma}$, that is, 1 < ${\gamma}$ or ${\gamma}$ > (5f/2-5/3)/(5f/2-1)when O < f < 2/5, (5f/2-5/3)/(5f/2-1) < ${\gamma}$ < 1 when f > 2/5. Since ${\epsilon}$' is negative, 4{\Omega}_0$,+ is again an imaginary quantity. For a negative $u_0$,+ with f < 0, ${\gamma}$ is allowed in 1 < 7 < (5|f|/2 + 5/3)/(5|f|/2 + 1). We briefly discuss physical implications of what we have found.