• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.2
• /
• pp.146-158
• /
• 1991

The flow around a two-dimensional foil section Is measured by a LDV(Laser Doppler Velocimetry) system which is capable of measuring the datailed flow field without interfering the original flow field. A 2-color 3-beam LDV system, which is capable of mea,;tiring 2 velocity components simultaneously and uses 2W Ar-Ion laser source, is used to measure the flow field around an NACA0012 foil section. The measured flow velocities are analysed iii order to study the boundary layer characteristics, flow separation and the detail structure of the flow near the trailing edge of the foil. The boundary layer characteristics are compared with the results by the head's momentum integral method. For the case of small angle of attack at relatively higher Reynolds number, both results show good agreements. The measured data of the velocity field around an NACA0012 foil section would be valuable data to validate the CFD(Computational Fluid Dynamic) calculation results. The developed experimental technique to evaluate the characteristics of two-dimensional foil sections is essential tool to develope new blade sections which have good lift characteristics and better cavitation performances.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.1
• /
• pp.94-101
• /
• 1994

This paper describes a short-term prediction model for the movement of an oil slick in shallow waters. Under the assumption that the initial movement of the oil slick is governed by spreading and advection, the model has been developed and applied to Kyungki-Bay near Incheon Harbor. The initial spreading is estimated by using an empirical formula. The depth-averaged momentum equations are solved numerically for the volume transport velocities, in which the $M_2$ tide is the main driving source. A staggered grid system is adopted fur spatial discretization and the half-time method is implemented for time marching. The numerical result is visualized with the help of animation and thus the contaminated area is displayed on a monitor in time sequence. The input data are the time, the location and the volume of spill accident as well as environmental data such as md and $M_2$ tide.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.24 no.2
• /
• pp.187-207
• /
• 2019

In order to develop a high performance ocean model, we used Julia, a Just-In-Time compile language, and to obtain the solution of the momentum equation, we made the code to solve the Poisson equation by the Successive Over-Relaxation method. And then we made two models to test Julia calculation codes. First, a simple channel form is modeled to test constant source/sink conditions. Second, the simplified Yellow Sea was modeled to test tidal forcing, Coriolis forces, and the effect of vertical eddy diffusivity coefficients. The model has been tested with a total of eight cases in the two scenarios. As a result of the test, the depth-averaged current speed of the three cases in Scenario 1 converged perfectly to the theoretical value, and that showed well a vertical flow velocity gradient due to the bottom friction. Also, the result of Scenario 2 represented well the amphidromic points of Yellow Sea and the tidal characteristics of mid-western and southwestern coast of Korea. Therefore, it is considered that the ocean model using Julia language has developed successfully, this suggests that the ocean model has come to the stage of successful transition from a classical compile language to a Just-In-Time compile language.

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

In recent years, countries like Europe and Japan have been involved in many researches on the Low-Crested Structure (LCS) which is the method to protect beach erosion and it is regarded as an alternative to the submerged breakwaters, and compiled its results and released the design manual. In the past, studies on LCS have focused on two-dimensional wave transmission and calculating required weight of armor units, and these were mainly examined and discussed based on experiments. In this study, three-dimensional numerical analysis is performed on permeable LCS. The open-source CFD code olaFlow based on the Navier-Stokes momentum equations is applied to the numerical analysis, which is a strongly nonlinear analysis method that enables breaking and turbulence analysis. As a result, the distribution characteristics of the LCS such as water level, water flow, and turbulent kinetic energy were examined and discussed, then they were carefully compared and examined in the case of submerged breakwaters. The study results indicate that there is a difference between the flow patterns of longshore current near the shoreline, the spatial distribution of longshore and on-offshore directions of mean turbulent kinetic energy in case of submerged breakwaters and LCS. It is predicted that the difference in these results leads to the difference in sand movement.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.1
• /
• pp.11-22
• /
• 1987

A plane buoyant jet discharged vertically upward into a crossflow is analyzed by numerical solution of the governing equations of continuity, momentum and constituent transport. The turbulent transport is modelled by the Prandtl's mixing length theory. In the numerical solution procedure, the governing equations are transformed by stream function and vorticity transport, non-dimensionalyzed by discharge velocity, slot width, and parameters representing flow characteristics, and solved by Gauss-Seidel iteration method with successive underrelaxation. The numerical experiments were performed for the region of established flow of buoyant jet in the range of discharge densimetric Froude number of 4 to 32 and in the range of velocity ratio of 8 to 15, which is the ratio of discharge velocity to crossflow velocity. Variations of velocities and temperatures, flow patterns and vorticity patterns of receiving water due to buoyant jet were investigated. Also investigated are the effects of velocity ratio and discharge densimetric Froude number on the trajectories of buoyant jet. Computed are velocities, temperatures and local densimetric Froude numbers along the trajectory of the buoyant jet. Spreading rate and dispersion ratio were analyzed in terms of discharge densimetric Froude number, local densimetric Froude number and distance from the source along the jet trajectory. It was noted that the similarity law holds in both the profiles of velocity and temperatures across the jet trajectory and the integral type analysis of Gaussian distribution is applicable.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2014.02a
• /
• pp.221-221
• /
• 2014

The present work proposes an improved numerical simulator for design and modification of large area capacitively coupled plasma (CCP) processing chamber. CCP, as notoriously well-known, demands the tremendously huge computational cost for carrying out transient analyses in realistic multi-dimensional models, because electron dissociations take place in a much smaller time scale (${\Delta}t{\approx}10-8{\sim}10-10$) than time scale of those happened between neutrals (${\Delta}t{\approx}10-1{\sim}10-3$), due to the rf drive frequencies of external electric field. And also, for spatial discretization of electron flux (Je), exponential scheme such as Scharfetter-Gummel method needs to be used in order to alleviate the numerical stiffness and resolve exponential change of spatial distribution of electron temperature (Te) and electron number density (Ne) in the vicinity of electrodes. Due to such computational intractability, it is prohibited to simulate CCP deposition in a three-dimension within acceptable calculation runtimes (<24 h). Under the situation where process conditions require thickness non-uniformity below 5%, however, detailed flow features of reactive gases induced from three-dimensional geometric effects such as gas distribution through the perforated plates (showerhead) should be considered. Without considering plasma chemistry, we therefore simulated flow, temperature and species fields in three-dimensional geometry first, and then, based on that data, boundary conditions of two-dimensional plasma discharge model are set. In the particular case of SiH4-NH3-N2-He CCP discharge to produce deposition of SiNxHy thin film, a cylindrical showerhead electrode reactor was studied by numerical modeling of mass, momentum and energy transports for charged particles in an axi-symmetric geometry. By solving transport equations of electron and radicals simultaneously, we observed that the way how source gases are consumed in the non-isothermal flow field and such consequences on active species production were outlined as playing the leading parts in the processes. As an example of application of the model for the prediction of the deposited thickness uniformity in a 300 mm wafer plasma processing chamber, the results were compared with the experimentally measured deposition profiles along the radius of the wafer varying inter-electrode gap. The simulation results were in good agreement with experimental data.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.5B
• /
• pp.473-482
• /
• 2009

This study numerically discusses wave forces acting on a vertical wall such as breakwaters or revetments, subjected to incident undular or turbulent bores. Due to the complex hydrodynamics of bore, its wave forces have been predicted, mainly through laboratory experiments. Numerical simulations in this paper were carried out by CADMAS-SURF(CDIT, 2001), which is based on Navier-Stokes momentum equations and VOF method (Hirt and Nichols, 1981) for tracking free water surface. Its original source code was also partly revised to generate bore in the numerical water channel. Numerical raw data computed by CADMAS-SURF included great strong spike phenomena that show the abrupt jumps of wave loads. To resolve this undesired noise of raw data, the band-pass filter with the frequency of 5Hz was utilized. The filtered results showed reasonable agreements with the experimental results performed by Matsutomi (1991) and Ramsden (1996). It was confirmed that CADMASSURF can be applied to the design of coastal structures against tsunami bores. In addition, the transformation process and propagation speed of bores in the same 2-d water channel were discussed by the variations of water level for time and space. The numerical results indicated that the propagation speed of bore was changed due to the nonlinear interactions between negative and reflected waves.