• Journal of Ocean Engineering and Technology
• /
• v.34 no.6
• /
• pp.406-418
• /
• 2020

In this study, a nonlinear wave simulation code is developed using a higher-order spectral (HOS) method. The HOS method is very efficient because it can determine the solution of the boundary value problem using fast Fourier transform (FFT) without matrix operation. Based on the HOS order, the vertical velocity of the free surface boundary was estimated and applied to the nonlinear free surface boundary condition. Time integration was carried out using the fourth order Runge-Kutta method, which is known to be stable for nonlinear free-surface problems. Numerical stability against the aliasing effect was guaranteed by using the zero-padding method. In addition to simulating the initial wave field distribution, a nonlinear adjusted region for wave generation and a damping region for wave absorption were introduced for wave generation simulation. To validate the developed simulation code, the adjusted simulation was carried out and its results were compared to the eighth order Stokes theory. Long-time simulations were carried out on the irregular wave field distribution, and nonlinear wave propagation characteristics were observed from the results of the simulations. Nonlinear adjusted and damping regions were introduced to implement a numerical wave tank that successfully generated nonlinear regular waves. According to the variation in the mean wave steepness, irregular wave simulations were carried out in the numerical wave tank. The simulation results indicated an increase in the nonlinear interaction between the wave components, which was numerically verified as the mean wave steepness. The results of this study demonstrate that the HOS method is an accurate and efficient method for predicting the nonlinear interaction between waves, which increases with wave steepness.

• Journal of the Korean GEO-environmental Society
• /
• v.15 no.3
• /
• pp.65-74
• /
• 2014

Borehole investigation which is mainly used to figure out geotechnical characterizations at construction work has the benefit that it provides a clear and convincing geotechnical information. But it has limitations to get the overall information of the construction site because it is performed at point location. In contrast, geophysical measurements like seismic survey has the advantage that the geological stratum information of a large area can be characterized in a continuous cross-section but the result from geophysics survey has wide range of values and is not suitable to determine the geotechnical design values directly. Therefore it is essential to combine borehole data and geophysics data complementally. Accordingly, in this study, a three-dimensional spatial interpolation of the cross-sectional distribution of seismic refraction was performed using digitizing and geostatistical method (krigring). In the process, digital map were used to increase the trustworthiness of method. Using this map, errors of ground height which are broken out in measurement from boring investigation and geophysical measurements can be revised. After that, average seismic velocity are derived by comparing borehole data with geophysical speed distribution data of each soil layer. During this process, outlier analysis is adapted. On the basis of the average seismic velocity, integrated analysis techniques to determine the three-dimensional geological stratum information is established. Finally, this analysis system is applied to dam construction field.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.3
• /
• pp.38-46
• /
• 1990

A numerical solution method of the boundary integral equation for axisymmetric potential flows is presented. Those are represented by ring source and ring vorticity distribution. Strengths of ring source and ring vorticity are approximated by linear functions of a parameter $\zeta$ on a segment. The geometry of the body is represented by a cubic B-spline. Limiting integral expressions as the field point tends to the surface having ring source and ring vorticity distribution are derived upto the order of ${\zeta}ln{\zeta}$. In numerical calculations, the principal value integrals over the adjacent segments cancel each other exactly. Thus the singular part proportional to $\(\frac{1}{\zeta}\)$ can be subtracted off in the calculation of the induced velocity by singularities. And the terms proportional to $ln{\zeta}$ and ${\zeta}ln{\zeta}$ can be integrated analytically. Thus those are subtracted off in the numerical calculations and the numerical value obtained from the analytic integrations for $ln{\zeta}$ and ${\zeta}ln{\zeta}$ are added to the induced velocity. The four point Gaussian Quadrature formula was used to evaluate the higher order terms than ${\zeta}ln{\zeta}$ in the integration over the adjacent segments to the field points and the integral over the segments off the field points. The root mean square errors, $E_2$, are examined as a function of the number of nodes to determine convergence rates. The convergence rate of this method approaches 2.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.20 no.4
• /
• pp.245-256
• /
• 2016

In the companion paper (I - Database and Site Response Analyses), site-specific response analyses were performed at more than 300 domestic sites. In this study, a new site classification system and design response spectra are proposed using results of the site-specific response analyses. Depth to bedrock (H) and average shear wave velocity of soil above the bedrock ($V_{S,Soil}$) were adopted as parameters to classify the sites into sub-categories because these two factors mostly affect site amplification, especially for shallow bedrock region. The 20 m of depth to bedrock was selected as the initial parameter for site classification based on the trend of site coefficients obtained from the site-specific response analyses. The sites having less than 20 m of depth to bedrock (H1 sites) are sub-divided into two site classes using 260 m/s of $V_{S,Soil}$ while the sites having greater than 20 m of depth to bedrock (H2 sites) are sub-divided into two site classes at $V_{S,Soil}$ equal to 180 m/s. The integration interval of 0.4 ~ 1.5 sec period range was adopted to calculate the long-period site coefficients ($F_v$) for reflecting the amplification characteristics of Korean geological condition. In addition, the frequency distribution of depth to bedrock reported for Korean sites was also considered in calculating the site coefficients for H2 sites to incorporate sites having greater than 30 m of depth to bedrock. The relationships between the site coefficients and rock shaking intensity were proposed and then subsequently compared with the site coefficients of similar site classes suggested in other codes.

• 한국신재생에너지학회:학술대회논문집
• /
• 2009.06a
• /
• pp.353-356
• /
• 2009

A pivotal mechanical balance of plant for 75kW class molten carbonate fuel cells comprise of a catalytic burner and an ejector which has been designed and tested in KEPRI(Korea Electric Power Research Institute). The catalytic burner, which oxidizes residual fuel in the anode tail gas, was operated at several conditions. Some problems arose due to local overheating or auto-ignition, which could limit the catalyst life. The catalytic burner was designed by considering both gas mixing and gas velocity. Test results showed that the temperature distribution is very uniform. In addition, an ejector is a fluid machinery to be utilized for mixing fluids, maintaining vacuum, and transporting them. The ejector is placed at mixing point between the anode off gas and the cathode off gas or the fresh air Several ejectors were designed and tested to form a suction on the fuel tail gas and balance the differential pressures between anode and cathode over a range of operating conditions. The tests showed that the design of the nozzle and throat played an important role in balancing the anode tail and cathode inlet gas pressures. The 75kW MCFC system built in our ejector and catalytic burner was successfully operated from Novembe, 2008 to April, 2009. It recorded the voltage of 104V at the current of 754A and reached the maximum generating power of 78.5kW DC. The results for both stand-alone and integration into another balance of plant are discussed.

• Geophysics and Geophysical Exploration
• /
• v.4 no.3
• /
• pp.80-83
• /
• 2001

Four seismic reflection horizons in 3-D seismic data, coherence derived from the seismic data, and 38 well logs from the Boonsville Gas Filed in Texas were tried to be integrated and visualized in 3 dimensions. Time surface was constructed from pick times of the reflection horizons. Average velocities to each horizon at 38 well locations were calculated based on depth markers from the well logs and time picks from the 3-D seismic data. The time surface was transformed to depth surface through velocity interpolation. Coherence was calculated on the 3-D seismic data by semblance method. Spatial distribution of the coherence is captured easily in 3-D visualization. Comparing to a time-slice of seismic data, distinctive stratigraphic features could be correctly recognized on the 3-D visualization.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.4
• /
• pp.73-81
• /
• 1994

The source and source/dipole distribution methods using 3-dimensional panel method have been widely used for ship motion analysis. When these methods are used, large errors in the predicted hydrodynamic coefficients are introduced around the irregular frequencies caused by the resonance of imaginary internal flow. Therefore, the irregular frequencies need to be removed for an accurate prediction of ship motion. This paper adopts 3-dimensional translating and oscillating Green function derived by Wu. The adaptive integration method, stretching transform and stationary phase method are used for the calculation of the calculation of Green function and the integral equation is derived by distributing the Green function n ship surface and inner free-surface. The condition of zero normal velocity, that is, wall condition on inner free-surface has been successfully used for the removal of irregular frequencies in oscillating problems. The calculations are carried out for series 60($C_B=0.7$) vessel and the results are compared with those of other theoretical analyses and experiment.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.3 no.1
• /
• pp.16-24
• /
• 1998

In this study, Ocean General Circulation Model (OGCM) has been developed as a counterpart of Atmospheric General Circulation (AGCM) for the study of coupled ocean-atmosphere climate system. The oceanic responses to given atmospheric boundary conditions have been investigated using the OGCM. In an integration carried out over 100 simulated years with climatological monthly mean data (EXP 1), most parts of the model reached a quasi-equilibrium climate reproducing many of the observed large-scale oceanic features remarkably well. Some observed narrow currents, however, such as North Equatorial Counter Current, were inevitably distorted due to the model's relatively coarse resolution. The seasonal changes in sea ice cover over the southern oceans around Antarctica were also simulated. In an experiment (EXP 2) under boundary condition of 10-year monthly data (1982-1991) from NCEP/NCAR Reanalysis Project model properly reproduced major oceanic changes during the period, including El Ni$\tilde{n}$os of 1982-1983 and 1986-87. During the ENSO periods, the experiment showed eastward expansion of warm surface waters and a negative vertical velocity anomalies along' the equator in response to expansion of westerly current velocity anomalies as westerly wind anomalies propagated eastward. Simulated anomalous distribution and the time behavior in response to El Ni$\tilde{n}$o events is consistent with that of the observations. These experiments showed that the model has an ability to reproduce major mean and anomalous oceanic features and can be effectively used for the study of ocean-atmosphere coupling system.

• Journal of Korean Society of Transportation
• /
• v.19 no.4
• /
• pp.35-47
• /
• 2001

Traffic represents one of the largest sources of primary air pollutants in urban area. As a consequence. numerous abatement strategies are being pursued to decrease the ambient concentration of pollutants. A characteristic of most of the these strategies is a requirement for accurate data on both the quantity and spatial distribution of emissions to air in the form of an atmospheric emission inventory database. In the case of traffic pollution, such an inventory must be compiled using activity statistics and emission factors for vehicle types. The majority of inventories are compiled using passive data from either surveys or transportation models and by their very nature tend to be out-of-date by the time they are compiled. The study of current trends are towards integrating urban traffic control systems and assessments of the environmental effects of motor vehicles. In this study, a methodology of motor vehicle emission calculation by using real-time traffic data was studied. A methodology for estimating emissions of CO at a test area in Seoul. Traffic data, which are required on a street-by-street basis, is obtained from induction loops of traffic control system. It was calculated speed-related mass of CO emission from traffic tail pipe of data from traffic system, and parameters are considered, volume, composition, average velocity, link length. And, the result was compared with that of a method of emission calculation by VKT(Vehicle Kilometer Travelled) of vehicles of category.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.21 no.2
• /
• pp.131-136
• /
• 1985

The calculation of the resulting fluid motion is an important problem of ship hydrodynamics. For a partially immersed body the condition of constant pressure at the free surface can be linearized. The resulting linear boundary-value problem for the velocity potential is the Neumann-Kelvin problem. The two-dimensional Neumann-Kelvin problem is studied for the half-immersed circular cylinder by Ursell. Maruo introduced a slender body approach to simplify the Neumann-Kelvin problem in such a way that the integral equation which determines the singularity distribution over the hull surface can be solved by a marching procedure of step by step integration starting at bow. In the present pater for the two-dimensional Neumann-Kelvin problem, it has been suggested that any solution of the problem must have singularities in the corners between the body surface and free surface. There can be infinitely many solutions depending on the singularities in the coroners.