• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.98-104
• /
• 2011

Numerical simulation in exploration geophysics provides important insights into subsurface wave propagation phenomena. Although elastic wave simulations take longer to compute than acoustic simulations, an elastic simulator can construct more realistic wavefields including shear components. Therefore, it is suitable for exploration of the responses of elastic bodies. To overcome the long duration of the calculations, we use a Graphic Processing Unit (GPU) to accelerate the elastic wave simulation. Because a GPU has many processors and a wide memory bandwidth, we can use it in a parallelised computing architecture. The GPU board used in this study is an NVIDIA Tesla C1060, which has 240 processors and a 102 GB/s memory bandwidth. Despite the availability of a parallel computing architecture (CUDA), developed by NVIDIA, we must optimise the usage of the different types of memory on the GPU device, and the sequence of calculations, to obtain a significant speedup of the computation. In this study, we simulate two- (2D) and threedimensional (3D) elastic wave propagation using the Finite-Difference Time-Domain (FDTD) method on GPUs. In the wave propagation simulation, we adopt the staggered-grid method, which is one of the conventional FD schemes, since this method can achieve sufficient accuracy for use in numerical modelling in geophysics. Our simulator optimises the usage of memory on the GPU device to reduce data access times, and uses faster memory as much as possible. This is a key factor in GPU computing. By using one GPU device and optimising its memory usage, we improved the computation time by more than 14 times in the 2D simulation, and over six times in the 3D simulation, compared with one CPU. Furthermore, by using three GPUs, we succeeded in accelerating the 3D simulation 10 times.

• 한국해양학회지
• /
• v.30 no.6
• /
• pp.565-583
• /
• 1995

We have used a nonlinear quasi-geostrophic model to study effects of lateral friction and bottom topography on the motion of warm eddies. The two empirical orthogonal functions of the stream function, accounting for the vertical structure, represent the barotropic and first baroclinic dynamic modes. This model is integrated 360 days on a 1000 km ${\times}$ 1000 km domain with a resolution of 10 km ${\times}$ 10 km including both the thermocline and idealized topography of the East Sea. Prescribed inflow through the Korea Strait is compensated by outflow through the Tsugaru Strait. The balance between the nonlinear advection term and the planetary ${\beta}$-effect tends to make northward movement of warm eddy over a flat bottom. The motion of a warm eddy over a sloping topography can be dominated by the nonlinear advection, while nonlinearity plays a secondary role over a flat topography. For eddies dispersing over topography, the nonlinear tendency is a function of time. For a strong warm eddy, northward propagation can occur. For intermediate strength of eddies one might expect a balance between the nonlinear term and the topographic ${\beta}$-effect. As nonlinearity decreases with eddy dispersion, southward motion along the slope may occur by such as a topographic Rossby wave. Our numerical simulations have confirmed the importance of lateral friction on eddy motions, in such a way that the northward penetration of the warm eddy increases drastically by the decrease of the lateral friction. The northward motion of warm eddy can be prevented by reducing the Reynolds number sufficiently. We have also demonstrated the crucial role of topographic effects in the eddy motion process.

• Journal of Navigation and Port Research
• /
• v.29 no.10 s.106
• /
• pp.853-858
• /
• 2005

Recently, various novel numerical models based on Navier-Stokes equation have been developed for calculating wave motions in the sea with coastal or ocean structures. Among those models, Volume Of Fluid (VOF) method might be the most popular one, and it has been used for numerical simulations of wave motions including complicated phenomena of wave breakings. VOF method, however, needs enormous computation time and large computational storage memories in general, thus it is practically difficult to use this method for calculations in the case of random waves because long and stable computation (e.g for more than 100 significant wave periods) is required to obtain statistically meaningful results. On the other hand if the wave motion is potential motion, Boundary Element Method (BEM), which is a much faster and more accurate method than VOF method, can be effectively used. The aim of this study is to develop a new efficient model applicable to calculations of wave motion and/or wave-structure interactions under random waves. To achieve this, a strictly combined BEM-VOF model has been developed by making the best use of both methods' merits; VOF method is used in a restricted fluid domain around a structure where complicated phenomena of wave breakings may exist, and BEM is used in the other domains far from the disturbance where the wave motion may be assumed to be potential. The verification of the model was performed with numerical results for Stokes' 5th order wave propagation and a random wave propagation.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.8 no.3
• /
• pp.23-27
• /
• 2008

The new conjunctive surface-subsurface flow model at a large scale was developed by using a 1-D Diffusion Wave (DW) model for surface flow interacting with the 3-D Volume Averaged Soil-moisture Transport (VAST) model for subsurface flow for the comprehensive terrestrial water and energy predictions in Land Surface Models (LSMs). A selection of numerical implementation schemes is employed for each flow component. The 3-D VAST model is implemented using a time splitting scheme applying an explicit method for lateral flow after a fully implicit method for vertical flow. The 1-D DW model is then solved by MacCormack finite difference scheme. This new conjunctive flow model is substituted for the existing 1-D hydrologic scheme in Common Land Model (CLM), one of the state-of-the-art LSMs. The new conjunctive flow model coupled to CLM is tested for a study domain around the Ohio Valley. The simulation results show that the interaction between surface flow and subsurface flow associated with the flow routing scheme matches the runoff prediction with the observations more closely in the new coupled CLM simulations. This improved terrestrial hydrologic module will be coupled to the Climate extension of the next-generation Weather Research and Forecasting (CWRF) model for advanced regional, continental, and global hydroclimatological studies and the prevention of disasters caused by climate changes.

• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• 2005.10a
• /
• pp.153-159
• /
• 2005

Recently, various novel numerical models based on Navier-Stokes equation rave been developed for calculating wave motions in the sea with coastal or ocean structures. Among those models, Volume Of Fluid (VOF) method might be the most popular one, and it has been used for numerical simulations of wave motions including complicated phenomena of wave breakings. VOF method, however, needs enormous computation time and large computational storage memories in general, thus it is practically difficult to use VOF method for calculations in the case of random waves because long and stable computation ( e.g. for more than 100 significant wave periods) is required to obtain statistically meaningful results. On the other hand of the wave motion is potential motion, Boundary Element Method (BEM), which is a much faster and more accurate method than VOF method, am be effectively used. The aim of this study is to develop a new efficient model applicable to calculations of wave motion and/or wave-structure interactions under random waves. To achieve this, a strictly combined BEM-VOF model has been developed by making the best use of both methods' merits; VOF method is used in a restricted fluid domain around a structure where complicated phenomena of wave breakings may exist, and BEM is used in the other domains far from the disturbance where the wave motion may be assumed to be potential. The verification of the model was performed with numerical results for Stokes'5th order wave propagation and a random wave propagation.

• Journal of the Korean Society of Radiology
• /
• v.15 no.6
• /
• pp.803-811
• /
• 2021

In this study, the exposure amount of IASCC test worker was evaluated by applying the process simulation technology. Using DELMIA Version 5, a commercial process simulation code, IASCC test facility, hot cells, and workers were prepared, and IASCC test activities were implemented, and the cumulative exposure of workers passing through the dose-distributed space could be evaluated through user coding. In order to simulate behavior of workers, human manikins with a degree of freedom of 200 or more imitating the human musculoskeletal system were applied. In order to calculate the worker's exposure, the coordinates, start time, and retention period for each posture were extracted by accessing the sub-information of the human manikin task, and the cumulative exposure was calculated by multiplying the spatial dose value by the posture retention time. The spatial dose for the exposure evaluation was calculated using MCNP6 Version 1.0, and the calculated spatial dose was embedded into the process simulation domain. As a result of comparing and analyzing the results of exposure evaluation by process simulation and typical exposure evaluation, the annual exposure to daily test work in the regular entrance was predicted at similar levels, 0.388 mSv/year and 1.334 mSv/year, respectively. Exposure assessment was also performed on special tasks performed in areas with high spatial doses, and tasks with high exposure could be easily identified, and work improvement plans could be derived intuitively through human manikin posture and spatial dose visualization of the tasks.