• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.25 no.5
• /
• pp.291-300
• /
• 2013

Typhoon simulation method is widely used to estimate sea surface wind speeds during the typhoon periods. Holland (1980) model has been regarded to provide relatively better results for observed wind data. JTWC or RSMC best track data are available for typhoon modeling, but these data show slightly different because the data generation process are different. In this paper, a Newton-Raphson method is used to solve the two nonlinear equations based on the Holland model that is formed by the two typhoon parameters, i.e. the longest radius of 25 m/s and 15 m/s wind speeds, respectively. The solution is the radius of maximum wind speed which is of importance for typhoon modeling. This method is based on the typhoon wind profile of JMA and it shows that Holland model appears to fit better the characteristics of typhoons on the temporal and spatial changes than that of the other models. In case of using RSMC best track data, the method suggested in this study shows better and more reasonable results for the estimation of radius of maximum wind speed because the consistency of the input data is assured.

• Journal of Korean Ophthalmic Optics Society
• /
• v.7 no.2
• /
• pp.1-11
• /
• 2002

A mathematical model is proposed to analyze the force; acting on the hard contact lens fitted on the cornea. The model incorporates the nonlinear equations and their numerical solution program, based on the formulations of surface tension force arising from the capillary action in the tear-film layer between the lens and cornea. The model simulates how the adhesion between lens and cornea varies according to the base curves and diameters of the lenses. When the spherical lens is fitted on the spherical cornea it is to rotate downward due to the weight of lens itself until it reaches an equilibrium position along the cornea where the counter(upward) moment caused by net force between the upper and lower portion of the periphery of lens. It is found that both the adhesion and displacement of lens along the cornea, where the gravity of lens balances the capillary-induced upward force, increases rapidly as the base curve of lens increases, i.e., as the lens gets flatter, while the increase in the diameter of lenses has resulted in the less increase in the rotation and adhesion. With the base curve and diameters of lenses being remained constant the increase in surface tension of tear film yields the increase in the adhesion between the cornea and lens while the initial rotation of lens is inversely proportional to the surface tension of the tear film.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.1
• /
• pp.81-92
• /
• 2013

The ultimate flexural strength of HSB I-girders, considering the effect of local bucking, was investigated through a series of nonlinear finite element analysis. The girders were selected such that the inelastic local flange buckling or the plastic yielding of compression flanges governs the flexural strength. Both homogeneous sections fabricated from HSB600 or HSB800 steel and hybrid sections with HSB800 flanges and SM570-TMC web were considered. In the FE analysis, the flanges and web were modeled using thin shell elements and initial imperfections and residual stresses were imposed on the FE model. An elasto-plastic strain hardening material was used for steels. After establishing the validity of present FE analysis by comparing FE results with test results published in the literature, the effects of initial imperfection and residual stress on the inelastic flange local buckling behavior were assessed. The ultimate flexural strengths of 60 I-girders with various compression flange slenderness were obtained by FE analysis and compared with those calculated from the KHBDC, AASHTO LRFD and Eurocode 3 provisions. Based on the comparison, the applicability of design equations in these specifications for the flexural strength of I-girder considering flange local buckling was evaluated.

• The Journal of the Acoustical Society of Korea
• /
• v.14 no.1E
• /
• pp.56-64
• /
• 1995

The least mean fourth (LMF) adaptive algorithm is a stochastic gradient method that minimizes the error in the mean fourth sense. Despite its potential advantages, the algorithm is much less popular than the conventional least mean square (LMS) algorithm in practice. This seems partly because the analysis of the LMF algorithm is much more difficult than that of the LMS algorithm, and thus not much still has been known about the algorithm. In this paper, we explore the statistical convergence behavior of the LMF algorithm when the input to the adaptive filter is zero-mean, wide-sense stationary, and Gaussian. Under a system idenrification mode, a set of nonlinear evolution equations that characterizes the mean and mean-squared behavior of the algorithm is derived. A condition for the conbergence is then found, and it turns out that the conbergence of the LMF algorithm strongly depends on the choice of initial conditions. Performances of the LMF algorithm are compared with those of the LMS algorithm. It is observed that the mean convergence of the LMF algorithm is much faster than that of the LMS algorithm when the two algorithms are designed to achieve the same steady-state mean-squared estimation error.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.4
• /
• pp.401-412
• /
• 2008

Simulations of three-dimensional numerical wave tank are performed to investigate wave force acting on a large cylindrical structure and consequent wave deformation, which are induced by bore after breaking waves. The numerical model is based on the three-dimensional Navier-Stokes equations with a finite-difference method combined with a volume of fluid(VOF) method, which is capable of tracking the complex free surface, including wave breaking. In order to promote wave breaking of the incident wave, the approach slope was built seaward of the structure with a constant slope and a large cylindrical structure was installed on a flat bed. The incident waves were broken on the approach slope or flat bed by its wave height. In the present study, all waves acting on the large cylindrical structure were limited to breaking bore after wave breaking. The effects of the position of the structure and the incident wave height on the wave force and wave transformations were mainly investigated with the concern of wave breaking. Further, the relations between the variation of wave energy by wave propagation after wave breaking and wave force acting on the structure were discussed to give the understanding of the full-linear wave-structure interactions in three-dimensional wave fields.

• Journal of Broadcast Engineering
• /
• v.1 no.2
• /
• pp.176-187
• /
• 1996

This paper deals with the accurate estimation of 3- D pose (position and orientation) of a moving object with reference to the world frame (or robot base frame), based on a sequence of stereo images taken by cameras mounted on the end - effector of a robot manipulator. This work is an extension of the previous work[1]. Emphasis is given to the 3-D pose estimation relative to the world (or robot base) frame under the presence of not only the measurement noise in 2 - D images[ 1] but also the camera position errors due to the random noise involved in joint angles of a robot manipulator. To this end, a new set of discrete linear Kalman filter equations is derived, based on the following: 1) the orientation error of the object frame due to measurement noise in 2 - D images is modeled with reference to the camera frame by analyzing the noise propagation through 3- D reconstruction; 2) an extended Jacobian matrix is formulated by combining the result of 1) and the orientation error of the end-effector frame due to joint angle errors through robot differential kinematics; and 3) the rotational motion of an object, which is nonlinear in nature, is linearized based on quaternions. Motion parameters are computed from the estimated quaternions based on the iterated least-squares method. Simulation results show the significant reduction of estimation errors and also demonstrate an accurate convergence of the actual motion parameters to the true values.

• Ocean and Polar Research
• /
• v.30 no.4
• /
• pp.453-466
• /
• 2008

In order to examine the generation mechanism of long ocean waves along the west coast of Korea and to understand the amplification process of the long ocean waves, sea level, atmospheric pressure and wind data observed every minute from 2007 March 29 to 2007 April 1 were analyzed and onedimensional numerical ocean model experiments were performed. An atmospheric pressure jump propagated southeastward from Backryungdo to Yeonggwang along the west coast of Korea with speed of $13{\sim}27\;m/s$ between 2007 March 30 23:00 and 2007 April 1 1:30. Average magnitude of pressure jump was 4.2 hPa. As a moving atmospheric jump propagated from north to south along the coast, long ocean waves were generated and the sea level abnormally rose or fell at Anheung, Kunsan, Wido and Yeonggwang. Average amplitude of sea level rise (or fall) was about 113.6 cm. In a one-dimensional numerical ocean model, nonlinear shallow water equations were numerically integrated and a moving atmospheric pressure jump with traveling speed of 24 m/s was used as an external force. While the atmospheric pressure jump travels over 60 m depth ocean, a long ocean wave is generated. Because the propagation speed of the atmospheric jump is almost equal to that of the long ocean wave, Proudman resonance occurs and the long ocean wave amplifies. As the atmospheric pressure jump moves into the coastal area shallower than 60 m, the speed of the long ocean wave decreases and Proudman resonance effect decreases. However, the amplitude of the long ocean wave increases and wave length becomes shorter because of shoaling effect. When the long ocean wave hits the land boundary, amplitude of the long ocean wave drastically amplifies due to reflection. Data analysis and numerical experiments suggest that the southeastward propagation of an atmospheric pressure jump over the shallow ocean, which is a necessary condition for Proudaman resonance, generated the long ocean waves along the west coast of Korea on 2007 March 31 and the ocean waves amplified due to shoaling effect in the coastal area and reflection at the shore.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.60 no.6
• /
• pp.43-54
• /
• 2018

The accurate estimation of reference crop evapotranspiration ($ET_o$) is essential in irrigation water management to assess the time-dependent status of crop water use and irrigation scheduling. The importance of $ET_o$ has resulted in many direct and indirect methods to approximate its value and include pan evaporation, meteorological-based estimations, lysimetry, soil moisture depletion, and soil water balance equations. Artificial neural networks (ANNs) have been intensively implemented for process-based hydrologic modeling due to their superior performance using nonlinear modeling, pattern recognition, and classification. This study adapted two well-known ANN algorithms, Backpropagation neural network (BPNN) and Generalized regression neural network (GRNN), to evaluate their capability to accurately predict $ET_o$ using daily meteorological data. All data were obtained from two automated weather stations (Chupungryeong and Jangsu) located in the Yeongdong-gun (2002-2017) and Jangsu-gun (1988-2017), respectively. Daily $ET_o$ was calculated using the Penman-Monteith equation as the benchmark method. These calculated values of $ET_o$ and corresponding meteorological data were separated into training, validation and test datasets. The performance of each ANN algorithm was evaluated against $ET_o$ calculated from the benchmark method and multiple linear regression (MLR) model. The overall results showed that the BPNN algorithm performed best followed by the MLR and GRNN in a statistical sense and this could contribute to provide valuable information to farmers, water managers and policy makers for effective agricultural water governance.

• Journal of Astronomy and Space Sciences
• /
• v.22 no.3
• /
• pp.263-272
• /
• 2005

For the Communication, Ocean and Meteorological Satellite (COMS) which will be launched in 2008, an algorithm for finding the precise location of the sun-glint point on the ocean surface is studied. The precise locations of the sun-glint are estimated by considering azimuth and elevation angles of Sun-satellite-Earth geometric position and the law of reflection. The obtained nonlinear equations are solved by using the Newton-Raphson method. As a result, when COMS is located at $116.2^{\circ}E$ or $128.2^{\circ}E$ longitude, the sun-glint covers region of ${\pm}10^{\circ}(N-S)$ latitude and $80-150^{\circ}(E-W)$ longitude. The diurnal path of the sun-glint in the southern hemisphere is curved towards the North Pole, and the path in the northern hemisphere is forwards the south pole. The algorithm presented in this paper can be applied to predict the precise location of sun-glint region in any other geostationary satellites.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.2
• /
• pp.608-616
• /
• 2018

If a general nonlinear system is linearized by the successive multiplication of the 1st and 2nd order systems, then there are four types of poles in this linearized system: the pole of the 1st order system and the equal poles, two distinct real poles, and complex conjugate pair of poles of the 2nd order system. Linear Quadratic (LQ) control is a method of designing a control law that minimizes the quadratic performance index. It has the advantage of ensuring the stability of the system and the pole placement of the root of the system by weighted matrix adjustment. LQ control by the weighted matrix can move the position of the pole of the system arbitrarily, but it is difficult to set the weighting matrix by the trial and error method. This problem can be solved using the characteristic equations of the Hamiltonian system, and if the control weighting matrix is a symmetric matrix of constants, it is possible to move several poles of the system to the desired closed loop poles by applying the control law repeatedly. The paper presents a method of calculating the state weighting matrix and the control law for moving the equal poles with Jordan blocks to two real poles using the characteristic equation of the Hamiltonian system. We express this characteristic equation with a state weighting matrix by means of a trigonometric function, and we derive the relation function (${\rho},\;{\theta}$) between the equal poles and the state weighting matrix under the condition that the two real poles are the roots of the characteristic equation. Then, we obtain the moving-range of the two real poles under the condition that the state weighting matrix becomes a positive semi-finite matrix. We calculate the state weighting matrix and the control law by substituting the two real roots selected in the moving-range into the relational function. As an example, we apply the proposed method to a simple example 3rd order system.