• Journal of the Korean Society of Propulsion Engineers
• /
• v.17 no.1
• /
• pp.18-25
• /
• 2013

Assessment and validation of RANS turbulence models are conducted for the optimal analysis of supersonic converging-diverging nozzle through the comparison between computational results and experimental data. One/two equation turbulence closures such as Spalart-Allmaras, RNG k-${\varepsilon}$, and k-${\omega}$ SST are employed to simulate the two-dimensional nozzle flow. Computational results with the turbulence models mentioned fairly well predict shock structure of the nozzle-inside and pressure distribution along the wall. Especially, SST model among the employed ones shows the best agreement to experimental results.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.675-678
• /
• 2008

The effect of wave and current interactions on jet-like current flowing against waves was investigated based on numerical simulations. The numerical simulations are conducted by a combined model system of REF/DIF(a wave model) plus SHORECIRC(a current model) and a Boussinesq equation model, FUNWAVE. In the simulations, regular and irregular waves refracted due to the jet-like opposing current were focused along the core region of current, and the jet-like current was earlier spreaded when the waves had larger wave heights. The numerical results show that the rapid change of wave height distribution in transverse direction near current inlet plays a significant role to spread the jet-like current. In other words, the gradients of radiation stress forcing in transverse direction have a more significant effect on the jet-like current than its relatively small gradients forcing in flowing direction, which tend to accelerate the current, do. In conclusion, it is indispensible to take into account the interaction effect of wave transformation and current characteristics when waves meet the opposing jet-like current such as river mouth.

• Journal of Korean Port Research
• /
• v.15 no.1
• /
• pp.47-56
• /
• 2001

A numerical model for practical use based on the 1-line theory is presented to simulate shoreline changes due to construction of offshore structures. The shoreline change model calculates the longshore sediment transport rate using breaking waves. Before the shoreline change model execution, a wave model, adopting the modified Boussinesq equation including the breaking parameters and bottom friction term, was used to provide the longshore distribution of the breaking waves. The contents of present model are outlined first. Then to examine the characteristics of this model, the effects of the parameters contained in this model are clarified through the calculations of shoreline changes for simple cases. Finally, as the guides for practical application of this model, several comments are made on the parameters used in the model, such as transport parameter, average beach slope, breaking height variation alongshore, depth of closure, etc. with the presentation of typical examples of 3-dimensional movable bed experimental results for application of this model. Here, beach change behind the offshore structures is represented by the movement of the shoreline position. Analysis gives that the transport parameters should be taken as site specific parameters in terms of time scale for the shoreline change and adjusted to achieve the best agreement between the calculated and the observed near the structures.

• International Journal of Aerospace System Engineering
• /
• v.2 no.2
• /
• pp.67-72
• /
• 2015

The main concern of this study is to integrate the EARSM into an industrial RANS solver in conjunction with the $k-{\omega}$ model, as proposed by Hellsten (EARSMKO2005). In order to improve the robustness, particular limiters are introduced to turbulent conservative variables, and a suitable full-approximation storage (FAS) multi-grid (MG) strategy is designed to incorporate turbulence model equations. The present limiters and MG strategy improve both robustness and efficiency significantly but without degenerating accuracy. Two discretization approachs for velocity gradient on cell interfaces are implemented and compared with each other. Numerical results of a three-dimensional supersonic square duct flow show that the proper discretization of velocity gradient improves the accuracy essentially. To assess the capability of the resulting EARSM $k-{\omega}$ model to predict complex engineering flow, the case of Common Research Model (CRM, Wing-Body) is performed. All the numerical results demonstrate that the resulting model performs well and is comparable to the standard two-equation models such as SST $k-{\omega}$ model in terms of computational effort, thus it is suitable for industrial applications.

• Journal of Navigation and Port Research
• /
• v.43 no.5
• /
• pp.320-327
• /
• 2019

This study examined the port expansion plan for the fishery port at the east coast of Korea, in accordance with permission conditions for coastal ports such as a limit on the cargo volume and passenger demand for the coastal tourism belt formation. The site was chosen as a municipal coastal port attracting the new ocean industry and building waterfront zone as a hub of new marine tourism. Two different numerical models (Swan and Bouss 2D) were used. Before applying to the target sea area, some numerical tests were conducted for the variation according to Bouss-2D's strong/weak and nonlinear technique compared to the irregular diffraction of semi-infinite breakwater with a theoretical solution. As a result, there was a difference in strong nonlinearity with breaking waves and it was necessary to experiment with a strong nonlinear analysis technique for the actual site. Two numerical models were applied to the fishery port site and the tranquility of some alternatives were analyzed. The numerical results show the most suitable plan was ALT-1, with satisfied harbor tranquility and reasonable economic sense. The extension of the east breakwater and enlarged turning basin of the F-Land plan have brought generally more stable harbor tranquility than the ALT-1. The result can be used as a reference for the port expansion plan in the future.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.4B
• /
• pp.243-251
• /
• 2012

To investigate the characteristics of rip current occurring at Haeundae beach, observations obtained from a buoy and a CCTV were analyzed and numerical experiments were conducted. During observed rip-current events, the CCTV images showed that a couple of wave-trains, which are close to regular waves with slightly different directions, propagated to the beach, and wavelet analyses of data from the buoy showed very narrow-banded spectra with a peak frequency. From the evidences, it was inferred that a known mechanism of generating rip current due to the nodal line area of honeycomb-patterned wave crest was one of the significant factors of rip current occurrences of Haeundae beach. The mechanism has been explained by the following: When two wave-trains with slightly different directions propagate to a beach, wave crests of the incident wave-trains form honeycomb pattern due to nonlinear interaction. The nodal lines of honeycomb pattern are developed in the cross-shore direction. And longshore currents flow toward the nodal line area which has very low wave energy. Consequently their mass flux is expelled through the area toward the sea direction. To confirm the generation, numerical experiments were performed using a nonlinear Boussinesq equation model. In the cases with two incident wave-trains with slightly different directions and with a monochromatic wave propagating over a submerged shoal, it was seen that the honeycomb pattern of wave crests was well developed, and thus rip currents were evolved along the nodal lines.

• Journal of Korea Water Resources Association
• /
• v.45 no.6
• /
• pp.545-555
• /
• 2012

A numerical modeling has become increasingly popular and more important to the study of water waves with a rapid advancement of computer technology. However, different types of problems are induced during simulating wave motion. One of the key problems is re-reflection to a computation domain at the incident boundary. The internal wave generating-absorbing boundary conditions have been commonly used in numerical wave models to prevent re-reflection. For the Navier-Stokes equations model, the internal wave maker using a mass source function of the continuity equation has been used to generate various types of waves. Nonetheless, almost every numerical experiment is performed in two dimensions and only a few tests have been expanded to three dimensions. More recently, a momentum source function of the Boussinesq equations is applied to generate essentially directional waves in the three dimensional Navier-Stokes equations model. In this study, the internal wave maker using a momentum source function is employed to generate targeted linear waves in the three-dimensional LES model.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2003.09a
• /
• pp.574-578
• /
• 2003

There are so many methods to estimate the groundwater recharge. These methods can be categorized into four groups. First groupis related to the water balance analysis, second group is concerned with baseflow/springflow recession, and third group is interested in some types of tracers; environmental tracers and/or temperature profile. The limitation of these types of methods is that the estimated results of recharge are presented in the form of an average over some time period. Forth group has a little different approach. They use the time series data of hydraulic head and specific yield evaluated from field test, and the results of estimation are described in the sequential form. But their approach has a serious problem. The estimated results in forth typeof methods are generally underestimated because they cannot consider the discharge phase of water table fluctuation coupled with the recharge phase. Ketchum el. at. (2000) proposed calibrated method, considering recharge- and discharge-coupled water table fluctuation. But the dischargeis considered just as the areal average with discharge rate. On the other hand, there are many methods to estimate the source wavelet with observed data set in geophysics/signal processing and geophysical methods are rarely applied to the estimation of groundwater recharge. The purpose this study is the evaluation of the applicability of one of the geophysical method in the estimation of sequential recharge rate. The applied geophysical method is called minimum entropy deconvolution (MED). For this purpose, numerical modeling with linearized Boussinesq equation was applied. Using the synthesized hydraulic head through the numerical modeling, the relative sequenceof recharge is calculated inversely. Estimated results are very concordant with the applied recharge sequence. Cross-correlations between applied recharge sequence and the estimated results are above 0.985 in all study cases. Through the numerical test, the availability of MED in the estimation of the recharge sequence to groundwater was investigated

• Journal of the Korean GEO-environmental Society
• /
• v.5 no.2
• /
• pp.51-62
• /
• 2004

Granular compaction piles have the load bearing capacity of the soft ground increase and have the settlement of foundation built on the reinforced soil reduce. The granular compaction group piles also have the consolidation of the soft ground accelerate and have the liquefaction caused by earthquake prevent using the granular materials such as sand, gravel, stone etc. However, this method is one of unuseful methods in Korea. The Granular compaction piles are constructed by grouping it with a raft system. The confining pressure at the center of bulging failure depth is a major variable in relation to estimate for the ultimate bearing capacity of the granular compaction piles. Therefore, a share of loading is determined considering the effect of load concentration ratio between the granular compaction piles and surrounding soils, and varies the magnitude of the confining pressure. In this study, method for the determination of the ultimate bearing capacity is proposed to apply a change of the horizontal pressure considering bulging failure depth, surcharge and loaded area. Also, the ultimate bearing capacity of the granular compaction piles is evaluated on the basis of previous study on the estimation of the ultimate bearing capacity and compared with the results obtained from laboratory scale model tests. And using the result from laboratory model tests, it is studied increase effect of the bearing capacity on the granular compaction piles and variance of coefficient of consolidation for the ground.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.4
• /
• pp.23-32
• /
• 1988

The rates of infiltration contributed to the flow fo water in an unconfined aquifer under the partially penetrated stream at an ungaged station and the corresponding base flow in channel are coupled by using the hydraulic and/or hydrologic characteristics obtained from the geomorphologic and soil maps. For the determination of groundwater flow, the linearized model which is originally Boussinesq's nonlinear equation is applied in this study. Also, a stream flow routing model for base flow in channel is based on a simplification of the Saint-venant. The distributed runoff model with piecewise spatial uniformity is presented for obtaining its solution based on a finite difference technique of the kinematic wave equations. The method developed in this study was tested to the Bocheong watershed(area : $475.5km^2$) of the natural stream basin which is one of tributaries in Geum River basin in Korea. As a result, it is suggested that the rationality of hydro-graph separation according to a wide variability in hydrogeologic properties be worked out as developing the physically based subsurface model. The results of the present model are shown to be possible to simulate a base flow due to an arbitrary rate of infiltration for ungaged basins.