• The Korean Journal of Applied Statistics
• /
• v.23 no.1
• /
• pp.189-196
• /
• 2010

It is never an easy task to physically randomize the sequence of cards. For instance, US 1970 draft lottery resulted in a social turmoil since the outcome sequence of 366 birthday numbers showed a significant relationship with the input order (Wikipedia, "Draft Lottery 1969", Retrieved 2009/05/01). We are motivated by Laplace's 1825 book titled Philosophical Essay on Probabilities that says "Suppose that the numbers 1, 2, ..., 100 are placed, according to their natural ordering, in an urn, and suppose further that, after having shaken the urn, to shuffle the numbers, one draws one number. It is clear that if the shuffling has been properly done, each number will have the same chance of being drawn. But if we fear that there are small differences between them depending on the order in which the numbers were put into the urn, we can decrease these differences considerably by placing these numbers in a second urn in the order in which they are drawn from the first urn, and then shaking the second urn to shuffle the numbers. These differences, already imperceptible in the second urn, would be diminished more and more by using a third urn, a fourth urn, &c." (translated by Andrew 1. Dale, 1995, Springer. pp. 35-36). Laplace foresaw what would happen to us in 150 years later, and, even more, suggested the possible tool to handle the problem. But he did omit the detailed arguments for the solution. Thus we would like to write the supplement in modern terms for Laplace in this research note. We formulate the problem with a lottery box model, to which Markov chain theory can be applied. By applying Markov chains repeatedly, one expects the uniform distribution on k states as stationary distribution. Additionally, we show that the probability of even-number of successes in binomial distribution with trials and the success probability $\theta$ approaches to 0.5, as n increases to infinity. Our theory is illustrated to the cases of truncated geometric distribution and the US 1970 draft lottery.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.1
• /
• pp.1-7
• /
• 2011

In this paper, based on an isogeometric approach, we have developed a shape design optimization method for plane elasticity problems subjected to design-dependent loads. The conventional shape optimization using the finite element method has some difficulties in the parameterization of geometry. In an isogeometric analysis, however, the geometric properties are already embedded in the B-spline basis functions and control points so that it has potential capability to overcome the aforementioned difficulties. The solution space for the response analysis can be represented in terms of the same NURBS basis functions to represent the geometry, which yields a precise analysis model that exactly represents the normal and curvature depending on the applied loads. A continuum-based isogeometric adjoint sensitivity is extensively derived for the plane elasticity problems under the design-dependent loads. Through some numerical examples, the developed isogeometric sensitivity analysis method is verified to show excellent agreement with finite difference sensitivity.

• The Journal of Engineering Geology
• /
• v.3 no.2
• /
• pp.149-165
• /
• 1993

Nowadays, the field hydraulic test is still an only method to evaluate groundwater characteristics in subsurface. The results of hydraulic test are very important for the concept model of fracture hydrogeology as well as the geometric pattern of fractures. The hydraulic tests performed in Korea are generally analysed under such assumption as steady radial flow in homogeneous aquifer or along simple geometry of fractures. Also the transmissivity measured in a fixed interval length is equivalent to a sum of individual fracture transmissivities in test legth. The boundary effects of weH hydraulics and the geometry of flow paths are hardly obtained from the test results analysed by a steady flow method. To circumvent this problem, the flow dimensional analysis was attempted from the results of constant pressure injection test carried out in a fractured granite area. A comparison of the hydraulic conductivity values from the transient and steady analysis shows that the latter is about a factor of 2~3 higher than the former. However, it was possible to analyse a flow dimension of each test interval from flow rate variation with time. The upper part of the bedrock(<10m deep) indicates an open boundary and the flow dimension shows nearly steady states, while the lower part of the bedrock(>25m deep) is characterized as sublinear flow dimension with a dosed boundary. In one of the test sections(15m deep), the flow dimension was changed from linear flow to spherical flow. From the experience of this study, one of the immediate problems to be solved is to enhance the field testing equipments, i.e., an accurate flowmeter with autorecording and a pressure detecting device to be able to install in the test section.

• 제어로봇시스템학회:학술대회논문집
• /
• 1991.10b
• /
• pp.1649-1654
• /
• 1991

In 1995 the VSOP satellite, which is called MUSES-B in Japan, will be launched under the VLBI Space Observatory Programme(VSOP) promoted by ISAS(Institute of Space and Astronautical Science) of Japan. We are now developing the GPS Receiver(GPSR) and On-board Orbit Determination System. This paper describes the GPS(Global Positioning System), VSOP, GPSR(GPS Receiver system) configuration and the results of the GPS system analysis. The GPSR consists of three GPS antennas and 5 channel receiver package. In the receiver package, there are two 16 bits microprocessing units. The power consumption is 25 Watts in average and the weight is 8.5 kg. Three GPS antennas on board enable GPSR to receive GPS signals from any NAVSTARs(GPS satellites) which are visible. NAVSATR's visibility is described as follows. The VSOP satellite flies from 1, 000 km to 20, 000 km in height on the elliptical orbit around the earth. On the other hand, the orbit of NAVSTARs are nearly circular and about 20, 000 km in height. GPSR can't receive the GPS signals near the apogee, because NAVSTARs transmit the GPS signals through the NAVSTAR's narrow beam antennas directed toward the earth. However near the perigee, GPSR can receive from 12 to 15 GPS signals. More than 4 GPS signals can be received for 40 minutes, which are related to GDOP(Geometric Dillusion Of Precision of selected NAVSTARs). Because there are a lot of visible NAVSTARs, GDOP is small near the perigee. This is a favorqble condition for GPSR. Orbit determination system onboard VSOP satellite consists of a Kalman filter and a precise orbit propagator. Near the perigee, the Kalman filter can eliminate the orbit propagation error using the observed data by GPSR. Except a perigee, precise onboard orbit propagator propagates the orbit, taking into account accelerations such as gravities of the earth, the sun, the moon, and other acceleration caused by the solar pressure. But there remain some amount of calculation and integration errors. When VSOP satellite returns to the perigee, the Kalman filter eliminates the error of the orbit determined by the propagator. After the error is eliminated, VSOP satellite flies out towards an apogee again. The analysis of the orbit determination is performed by the covariance analysis method. Number of the states of the onboard filter is 8. As for a true model, we assume that it is based on the actual error dynamics that include the Selective Availability of GPS called 'SA', having 17 states. Analytical results for position and velocity are tabulated and illustrated, in the sequel. These show that the position and the velocity error are about 40 m and 0.008 m/sec at the perigee, and are about 110 m and 0.012 m/sec at the apogee, respectively.

• Korean Journal of Remote Sensing
• /
• v.30 no.6
• /
• pp.797-807
• /
• 2014

The Sea Surface Temperature (SST) is one of the most important oceanic environmental factors in determining the change of marine environments and ecological activities. Satellite thermal infrared images can be effective for understanding the global trend of sea surface temperature due to large scale. However, their low spatial resolution caused some limitations in some areas where complicated and refined coastal shapes due to many islands are present as in the Korean Peninsula. The coastal ocean is also very important because human activities interact with the environmental change of coastal area and most aqua farming is distributed in the coastal ocean. Thus, low-cost airborne thermal infrared remote sensing with high resolution capability is considered for verifying its possibility to extract SST and to monitor the changes of coastal environment. In this study, an airborne thermal infrared system was implemented using a low-cost and ground-based thermal infrared camera (FLIR), and more than 8 airborne acquisitions were carried out in the western coast of the Korean Peninsula during the periods between May 23, 2012 and December 7, 2013. The acquired thermal infrared images were radiometrically calibrated using an atmospheric radiative transfer model with a support from a temperature-humidity sensor, and geometrically calibrated using GPS and IMU sensors. In particular, the airborne sea surface temperature acquired in June 25, 2013 was compared and verified with satellite SST as well as ship-borne thermal infrared and in-situ SST data. As a result, the airborne thermal infrared sensor extracted SST with an accuracy of $1^{\circ}C$.

• Tunnel and Underground Space
• /
• v.15 no.1 s.54
• /
• pp.28-38
• /
• 2005

According to the recent development of measurement system utilizing one or a set of boreholes, visualization of the explored underground became to be a major issue. It induced even the introduction of monitoring apparatuses on the borehole wall with multi-function tool, but the usage of these was often limited by where is unfavorable rock condition and a few of engineers can approach. And so, a portable type of borehole camera with only the essential function has been investigated and a few of commercial models about this is recently being applied into the field condition. This paper was based on the monitoring results obtained using a commercial model by Dr. Nakagawa. Discontinuities in rock mass were the topic for the visualization, and it was studied how can visualize their three dimensional distribution and what a numerical formulation is needed and how to understand the visualization result. The numerical formulation was based on the geometric correlation between the dip direction / dip of discontinuous plane and the trend / plunge of borehole, a set of the equation of a plane was induced. As field application of this into two places, it is found that the above visualization methodology will be especially an useful geotechlical tool for analyzing the local distribution of discontinuities.

• Journal of Korean Society of Transportation
• /
• v.23 no.7 s.85
• /
• pp.77-86
• /
• 2005

Level of service(LOS) is a quantify measure describing operational conditions within a traffic stream, generally, in terms of such service measures as speed, travel time, freedom to measures, traffic interruptions, comfort and convenience. The LOS is leveled by highway facilities according to measure of effectiveness(MOE) and then used to evaluate performance capacity. The current evaluation of a urban road is performed by only a aspect of traffic operation without any concepts of safety. Therefore, this paper presents a method for evaluation of risk order for urban road with new MOE, user cost analysis, considering both smooth traffic operation(congestion) and traffic safety(accident). The user coat is included traffic accident cast by traffic safety and traffic congestion cost by traffic operation. First of all, a number of traffic accident and accident rate by highway geometric is inferred from urban road traffic accident prediction model (Poul Greibe(2001)) Secondly, a user cost is inferred as traffic accident cast and traffic congestion cost is putting together. Thirdly, a method for evaluation of a urban road is inferred by user cost analysis. Fourthly a accident rate by segment predict with traffic accidents and data related to the accidents in $1996{\sim}1998$ on 11 urban road segments, Gwang-Ju, predicted accident rate. Traffic accident cost predict using predicted accident rate, and, traffic congestion cost predict using predicted average traffic speed(KHCM). Fifthly, a risk order are presented by predicted user cost at each segment in urban roads. Finally, it si compared and evaluated that LOS of 11 urban road segments, Gwang-Ju, by only a aspect of traffic operation without any concepts of safety and risk order by a method for evaluation of urban road in this paper.

• Journal of Korea Water Resources Association
• /
• v.38 no.10 s.159
• /
• pp.801-810
• /
• 2005

The lag time is one of the most important factors for estimating a flood runoff from streams. It is well known to be under the influence of the morphometric properties of basins which could be expressed by catchment shape descriptors. In this paper, the notion of the geometric characteristics of an equivalent ellipse proposed by Moussa(2003) is applied for calculating the lag time of geomorphological instantaneous unit hydrograph(GIUH) at the basin outlet. The lag time is obtained from the observed data of rainfall and runoff by using the method of moments suggested by Nash(1957), and the procedure based on geomorphology is used for GIUH. The relationships between the basin morphometric properties and the hydrological response are discussed as applied to 3 catchments In Korea. Additionally, the shapes of equivalent ellipse are examined how then are transformed from upstream area to downstream one. As a result, the relationship between the hydrological response and descriptors is shown to be comparatively good, and the shape of ellipse is presented to approach a circle along the river downwards. These results may be expanded to the estimation of hydrological response of ungauged catchment.

• Proceedings of KOSOMES biannual meeting
• /
• 2006.11a
• /
• pp.181-186
• /
• 2006

Geostationary Orbit satellite, unlike other sun-synchronous polar-orbit satellites, will be able to take a picture of a large region several times a day (almost with everyone hour interval). For geostationary satellite, the target region is fixed though the location of sun is changed always. However, Sun-synchronous polar-orbit satellites able to take a picture of target region same time a everyday. Thus Ocean signal is almost same. Accordingly, the ocean signal of a given target point is largely dependent on time. In other words, the ocean signal detected by geostationary satellite sensor must translate to the signal of target when both sun and satellite are located in nadir, using another correction model. This correction is performed with a standardization of signal throughout relative geometric relationship among satellite-sun-target points. This relative ratio called bidirectional factor. To find relationship between time and $[L_w]_N$/Bidirectional Factor differences, we are calculate solar position, geometry parameters. And reflectance, total radiance at the top of atmosphere(). And water leaving radiance, normalized water leaving radiance. And calculate bidirectional factor, that is the ratio of $[L_w]_N$ between target region and aiming the point. Then, we can make the bidirectional factor lookup table for one year imaging. So, we suggested for necessary to simulation experiment bidirectional factor in more various condition(wavelength and ocean/air condition).

• Ocean Systems Engineering
• /
• v.1 no.4
• /
• pp.285-296
• /
• 2011

Spar platforms have several advantages for deploying wind turbines in offshore for depth beyond 120 m. The merit of spar platform is large range of topside payloads, favourable motions compared to other floating structures and minimum hull/deck interface. The main objective of this paper is to present the response analysis of the slack moored spar platform supporting 5MW wind turbine with bottom keel plates in regular and random waves, studied experimentally and numerically. A 1:100 scale model of the spar with sparD, sparCD and sparSD configuration was studied in the wave basin ($30{\times}30{\times}3m$) in Ocean engineering department in IIT Madras. In present study the effect of wind loading, blade dynamics and control, and tower elasticity are not considered. This paper presents the details of the studies carried out on a 16 m diameter and 100 m long spar buoy supporting a 90 m tall 5 MW wind turbine with 3600 kN weight of Nacelle and Rotor and 3500 kN weight of tower. The weight of the ballast and the draft of the spar are adjusted in such a way to keep the centre of gravity below the centre of buoyancy. The mooring lines are divided into four groups, each of which has four lines. The studies were carried out in regular and random waves. The operational significant wave height of 2.5 m and 10 s wave period and survival significant wave height of 6 m and 18 s wave period in 300 m water depth are considered. The wind speed corresponding to the operational wave height is about 22 knots and this wind speed is considered to be operating wind speed for turbines. The heave and surge accelerations at the top of spar platform were measured and are used for calculating the response. The geometric modeling of spar was carried out using Multisurf and this was directly exported to WAMIT for subsequent hydrodynamic and mooring system analysis. The numerical results were compared with experimental results and the comparison was found to be good. Parametric study was carried out to find out the effect of shape, size and spacing of keel plate and from the results obtained from present work ,it is recommended to use circular keel plate instead of square plate.