• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.4
• /
• pp.313-325
• /
• 2011

Finite element grid system using h-refinement on the Yellow Sea was constructed based on previous study (Suh, 1999b) from 14 K to 210 K and special attention was concentrated on refining the coastal zone. In grid generation, depth change between adjacent points and non-dimensional tidal wave length ratio were considered. As a result approximately a quarter of the total nodes are located nearby 5 m of shallow area. Accurate bathymetry data using 30's and ETOPO1 with open boundary conditions of 8 major tidal constituents extracted automatically from FES2004 have been applied. In tidal simulation a 3-dimensional nonlinear harmonic model was setup and tidal amplification due to changes in vertical turbulent and bottom friction were simulated. In this study not only 8 major tidal constituents but also nonlinear shallow tides $M_4,$, $MS_4$ and long period $M_f,$, $M_{sf}$ were reproduced. It is found that implication of spatial variation of friction coefficient plays a very important role in reproduction of astronomical and shallow tides which are computed by iterative computation of nonlinear terms. Also it should be considered differently with respect to tidal periods. To understand the distribution of tidal asymmetry, amplitude ratio of $M_4/M_2$ and phase differences $2g(M_2)-g(M_4)$ were calculated. Tidal distortion ratio marks up to 0.2 on the west coast showing shallow coastal characteristics and somewhat wide range of ebb-dominances in front of Mokpo area are reproduced.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.3
• /
• pp.163-171
• /
• 2024

Structural health monitoring for ships and offshore structures is important in various aspects. Ships and offshore structures are continuously exposed to various environmental conditions, such as waves, wind, and currents. In the event of an accident, immense economic losses, environmental pollution, and safety problems can occur, so it is necessary to detect structural damage or defects early. In this study, structural response data of multi-linked floating offshore structures under various wave load conditions was calculated by performing fluid-structure coupled analysis. Furthermore, the order reduction method with distortion base mode was applied to the structures for predicting the structural response by using the results of numerical analysis. The distortion base mode order reduction method can predict the structural response of a desired area with high accuracy, but prediction performance is affected by sensor arrangement. Optimization based on a genetic algorithm was performed to search for optimal sensor arrangement and improve the prediction performance of the distortion base mode-based reduced-order model. Consequently, a sensor arrangement that predicted the structural response with an error of about 84.0% less than the initial sensor arrangement was derived based on the root mean squared error, which is a prediction performance evaluation index. The computational cost was reduced by about 8 times compared to evaluating the prediction performance of reduced-order models for a total of 43,758 sensor arrangement combinations. and the expected performance was overturned to approximately 84.0% based on sensor placement, including the largest square root error.

• The Journal of the Acoustical Society of Korea
• /
• v.13 no.6
• /
• pp.75-82
• /
• 1994

The ray paths and travel times of sound wave in the ocean depend on the physical properties of the propagating media. Ocean Acoustic Tomography(OAT), which is inversely estimate the travel time variations between fixed sources and receivers the physical properties of the corresponding media can he understood. To apply ocean survey technology by using the OAT, the tomographic procedure requires forward problem that variation of the travel times be identified with the variability of the medium. Also, received signals must be satisfied the necessary conditions of ray path stability, identification and resolution in order for OAT to work. The canonical ocean has been determined based on the historical data and its travel time and ray path are used as reference values. The sound speed of canonical ocean in the East Sea is about 1523 m/s at the surface and 1458 m/s at the sound channel axis(400m). Sound speeds in the East Sea are perturbed by warm eddy whose horizontal extension is more than 100 km with deeper than 200 m in depth scale. In this study, an acoustic source and receiver are placed at the depth above the sound channel axis, 350 m, and are separated by 200 km range. Ray paths are identified by the ray theory methed in a range dependent medium whose sound speeds are functions of a range and depth. The eigenray information obtained from interpolation between the rays bracketing the receiver are used to simulate the received signal by convolution of source signal with the eigenray informations. The source signal is taken as a 400 Hz rectangular pulse signal, bandwidth is 16 Hz and pulse length is 64 ms. According to the analysis of the received signal and identified ray path by using numerical model of underwater sound propagation, simulated signals satisfy the necessary conditions of OAT, applied in the East Sea.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.30 no.6
• /
• pp.306-318
• /
• 2018

In order to evaluate the seismic capacity of massive vertical type breakwaters which have intensively been deployed along the coast of South Korea over the last two decades, we carry out the preliminary numerical simulation against the PoHang, GyeongJu, Hachinohe 1, Hachinohe 2, Ofunato, and artificial seismic waves based on the measured time series of ground acceleration. Numerical result shows that significant sliding can be resulted in once non-negligible portion of seismic energy is shifted toward the longer period during its propagation process toward the ground surface in a form of shear wave. It is well known that during these propagation process, shear waves due to the seismic activity would be amplified, and non-negligible portion of seismic energy be shifted toward the longer period. Among these, the shift of seismic energy toward the longer period is induced by the viscosity and internal friction intrinsic in the soil. On the other hand, the amplification of shear waves can be attributed to the fact that the shear modulus is getting smaller toward the ground surface following the descending effective stress toward the ground surface. And the weakened intensity of soil as the number of attacking shear waves are accumulated can also contribute these phenomenon (Das, 1993). In this rationale, we constitute the numerical model using the model by Hardin and Drnevich (1972) for the weakened shear modulus as shear waves go on, and shear wave equation, in the numerical integration of which $Newmark-{\beta}$ method and Modified Newton-Raphson method are evoked to take nonlinear stress-strain relationship into account. It is shown that the numerical model proposed in this study could duplicate the well known features of seismic shear waves such as that a great deal of probability mass is shifted toward the larger amplitude and longer period when shear waves propagate toward the ground surface.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.3
• /
• pp.101-114
• /
• 2019

In this study, we carry out the numerical simulation to trace the yearly shoreline change of Mang-Bang beach, which is suffering from erosion problem. We obtain the basic equation (One Line Model for shoreline) for the numerical simulation by assuming that the amount of shoreline retreat or advance is balanced by the net influx of longshore and cross-shore sediment into the unit discretized shoreline segment. In doing so, the energy flux model for the longshore sediment transport rate is also evoked. For the case of cross sediment transport, the modified Bailard's model (1981) by Cho and Kim (2019) is utilized. At each time step of the numerical simulation, we adjust a closure depth according to pertinent wave conditions based on the Hallermeier's analytical model (1978) having its roots on the Shield's parameter. Numerical results show that from 2017.4.26 to 2017.10.15 during which swells are prevailing, a shoreline advances due to the sustained supply of cross-shore sediment. It is also shown that a shoreline temporarily retreats due to the erosion by the yearly highest waves sequentially occurring from mid-October to the end of October, and is followed by gradual recovery of shoreline as high waves subdue and swells prevail. It is worth mentioning that great yearly circulation of shoreline completes when a shoreline retreats due to the erosion by the higher waves occurring from mid-March to the end of March. The great yearly circulation of shoreline mentioned above can also be found in the measured locations of shoreline on 2017.4.5, 2017.9.7, 2017.11.7, 2018.3.14. However, numerically simulated amount of shoreline retreat or advance is more significant than the physically measured one, and it should be noted that these discrepancies become more substantial for the case of RUN II where a closure depth is sustained to be as in the most morphology models like the Genesis (Hanson and Kraus, 1989).

• Geophysics and Geophysical Exploration
• /
• v.25 no.4
• /
• pp.227-241
• /
• 2022

Full waveform inversion (FWI) in the field of seismic data processing is an inversion technique that is used to estimate the velocity model of the subsurface for oil and gas exploration. Recently, deep learning (DL) technology has been increasingly used for seismic data processing, and its combination with FWI has attracted remarkable research efforts. For example, DL-based data processing techniques have been utilized for preprocessing input data for FWI, enabling the direct implementation of FWI through DL technology. DL-based FWI can be divided into the following methods: pure data-based, physics-based neural network, encoder-decoder, reparameterized FWI, and physics-informed neural network. In this review, we describe the theory and characteristics of the methods by systematizing them in the order of advancements. In the early days of DL-based FWI, the DL model predicted the velocity model by preparing a large training data set to adopt faithfully the basic principles of data science and apply a pure data-based prediction model. The current research trend is to supplement the shortcomings of the pure data-based approach using the loss function consisting of seismic data or physical information from the wave equation itself in deep neural networks. Based on these developments, DL-based FWI has evolved to not require a large amount of learning data, alleviating the cycle-skipping problem, which is an intrinsic limitation of FWI, and reducing computation times dramatically. The value of DL-based FWI is expected to increase continually in the processing of seismic data.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.219-235
• /
• 2014

With the development of the technology of underwater moving bodies, the need for developing the knowledge of surface effect interaction of free surface and underwater moving bodies is increased. Hence, the two-phase flow is a subject which is interesting for many researchers all around the world. In this paper, the non-linear free surface deformations which occur during the water-exit of a circular cylinder due to its buoyancy are solved using finite volume discretization based code, and using Volume of Fluid (VOF) scheme for solving two phase flow. Dynamic mesh model is used to simulate dynamic motion of the cylinder. In addition, the effect of cylinder mass in presence of an external force is studied. Moreover, the oblique exit and entry of a circular cylinder with two exit angles is simulated. At last, water-exit of a circular cylinder in six degrees of freedom is simulated in 3D using parallel processing. The simulation errors of present work (using VOF method) for maximum velocity and height of a circular cylinder are less than the corresponding errors of level set method reported by previous researchers. Oblique exit shows interesting results; formation of waves caused by exit of the cylinder, wave motion in horizontal direction and the air trapped between the waves are observable. In 3D simulation the visualization of water motion on the top surface of the cylinder and the free surface breaking on the front and back faces of the 3D cylinder at the exit phase are observed which cannot be seen in 2D simulation. Comparing the results, 3D simulation shows better agreement with experimental data, specially in the maximum height position of the cylinder.

• The Journal of the Acoustical Society of Korea
• /
• v.39 no.1
• /
• pp.1-7
• /
• 2020

The acoustic parabolic equation method in the ocean is an efficient technique to calculate the acoustic field in the range-dependent environment, emanating from a point source. However, we often need to use the directional source with a main beam in the practical problem. In this paper, we present two methods to implement the directional source in the acoustic parabolic equation code easily. One is simply to filter the Delta function idealized as an omni-directional point source. Another method is based on the rational filtering of the self-starter solution. It has a limitation not to separate the up-going and the down-going wave for the depth, but would be useful in implementing the mode propagation. Numerical examples for validation are given in the Pekeris environment and the deep sea environment.

• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.16 no.7
• /
• pp.1344-1349
• /
• 2012

The Korean DGPS station transmits the 200 bps GPS enhancement signal using the MSK modulation in frequency range of 283.5 kHz to 325 kHz. The land-based stations of 6 sites provide the service area of 80 km with the output power of 500 W. The ocean-based stations of 11 sites provide the output power of 300 W, which provide the DGPS service to 185 kM. Some places are serviced from two or three DGPS stations. The interferences among the DGPS stations using the high power can be occurred. Also, the performances of the user terminasl in dual service area can be degraded. In this paper, the protection ratios for the DGPS service are defined. Using the MF wave propagation model, the interferences among the DGPS stations and the adjacent wireless ground stations are analyzed. Also, the performances of DGPS user terminals are analyzed in the viewpoint of interference.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.4 no.1
• /
• pp.47-62
• /
• 2001

An experimental study was peformed to provide the numerical models to predict the physical fate of marine debris with the key information concerning their trifling and sinkage on the sea. For hygroscopic materials their floating times in a pseudo ocean environment were measured, and for non-hygroscopic objects their drifting limes in a certain wave condition were measured. For sinking materials or objects became to be sinkable free falling experiment was performed to evaluate their falling speed. Summarized information was condensed into a numerical model to predict the most possible region of a significant marine debris accumulation, and the seas connected to the Han river and the Nak-Dong river were considered for the simulation.