• Journal of the Korean Society of Marine Environment & Safety
• /
• v.9 no.2
• /
• pp.53-57
• /
• 2003

Low salinity water appears by outflow of fresh water from the Yangtze River in southwestern sea of the Korean peninsula. The water volume discharged form the Yangtze River is not constant with year, according to the time series data recorded in the past, the maximum value of the discharged volume is two times as high asr minimum value. Therefore, the distribution of low salinity water over study area with different discharged fresh water volume is studied using numerical model, Princeton Ocean Model (POM). POM is three dimensional circulation model selecting a $\sigma$­coordinate. According to the result of numerical simulation by the Model, current velocity on the continental slope is faster than those in other regions, current which flows toward the central part of the Yellow Sea through western part of Jeju exists, and also, southward flow along the coastal region exists. the greater discharged volume from the Yangrze River is. the lower salinity water appears closer to Jeju.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.8
• /
• pp.20-28
• /
• 2004

For the prediction on the onset of oblique shock wave-induced vortex breakdown, computational studies on the Oblique Shock wave/Vortex Interaction (OSVI) are conducted and compared with both experimental results and analytic mode1. A Shock-stable numerical scheme, the Roe scheme with Mach number-based function (RoeM), and a two-equation eddy viscosity-transport approach arc used for three-dimensional turbulent flow computations. The computational configuration is identical to available experiment, and we attempt to ascertain the effect of parameters such as a vortex strength, streamwise velocity deficit, and shock strength at a freestream Mach number of 2.49. Numerical simulations using the k-w SST turbulence model and suitably modeled vortex profiles are able to accurately reproduce many fine features through a direct comparison with experimental observations. The present computational approach to determine the criterion on the onset of oblique shock wave-induced vortex breakdown is found to be in good agreement with both the experimental result and the analytic prediction.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.15 no.4
• /
• pp.273-280
• /
• 2012

The numerical simulation using the Lattice Boltzmann Method in the field of computational fluid dynamics becomes wider in the engineering applications because of its simplicity of update rules compared to the conventional Navier-Stokes solvers. Here, a two-dimensional D2Q9 LB model is numerically tested with a few new computational treatment on the free surface. The single relaxation time is applied under the gravitational field where applied only in the higher density fluid because of its big density difference. At the free surface, the reconstruction techniques in combination with boundary conditions is adopted in order to get some distribution function coming into the fluid site from the air one, and surface tension, early stable test for the gravitional field is considered in it. With the implementation of the gravitational profile, conserving the overall mass and grid dependency are observed during the calculations and freesurface advance track is well captured with an experiment.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.4
• /
• pp.243-258
• /
• 2005

A fully-nonlinear two-dimensional dynamically interfaced nested tidal model has been developed with improved efficiency, usability, relocatability and also adopting different time steps for economizing computational times, which is based on two-dimensional tidal model of Choi(1990) with resolution of $1/15^{\circ}$ inlatitude by $1/12^{\circ}$ in longitude covering the whole Yellow and East China Seas continental shelf. To illustrate the validity and applicability of the present model, numerical model experiment was carried out for tidal distribution in the Yellow Sea and the East China Sea emphasizing the Saemangeum area with enhanced grid system. The results from comparison with observed value showed that the present model gives reasonable agreement. The present modeling system is a relocatable tide forecasting system designed to predict tides at any location within the Yellow Sea and the East China Sea continental shelf sea.

• Journal of Korea Water Resources Association
• /
• v.54 no.6
• /
• pp.407-418
• /
• 2021

In this study, the two-dimensional flow model, Nays2DH, and driftwood dynamics model were combined to analyze the flow and driftwood behavior depending on the characteristics of the inflow of driftwood and the length of the driftwood stem. In particular, the Dashpot-spring model was added to the driftwood dynamics model to simulate the collision motion of the driftwood, and the wood jam characteristics by the collision of the driftwood were compared. As a result of the simulation, the pass rate of the obstacle section, the travel distance of wood jam, and the mean position of the wood pieces were respondent sensitively by the length of the driftwood stem, but the cohort size of the driftwood supply was insignificant excepting for the pass rate. Through this study, we could understand the interaction between hydraulic structures and driftwood, and through this, it is believed that it will be helpful in establishing a durable maintenance plan for hydraulic structures by predicting the transport and jam formation phenomena of driftwood in advance.

• Geophysics and Geophysical Exploration
• /
• v.9 no.3
• /
• pp.241-249
• /
• 2006

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms commonly used, however, may not give sufficiently precise computational results of traveltime data particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. Considering the complex geology of Korea, we assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution. The performance of the algorithm developed in this study is demonstrated by the comparison of the analytic and numerical solutions for the homogeneous anisotropic earth as well as through the numerical experiment for the two layer model whose anisotropic properties are greatly different each other. We expect that the developed modeling algorithm can be used in the development of processing and inversion schemes of seismic data acquired in strongly anisotropic environment, such as migration, velocity analysis, cross-well tomography and so on.

• 한국연소학회:학술대회논문집
• /
• 2002.06a
• /
• pp.52-60
• /
• 2002

A two-dimensional direct numerical simulation is performed to investigate the flame structure of $CH_4/N_2-Air$ counterflow nonpremixed flame interacting with a single vortex. The detailed transport properties and a modified 16-step augmented reduced mechanism based on Miller and Bowman's detailed reaction mechanism are adopted in this calculation. To quantify the strain on flame induced by a vortex, a scalar dissipation rate (SDR) is introduced. Results show that the fuel and air-side vortex cause an unsteady extinction. In this case, the flame interacting with a vortex is extinguished at much larger SDR than steady flame. It is also found that air-side vortex extinguishes a flame more rapidly than fuel-side vortex. The unsteady effect induced by flame-vortex interaction does not lead to a transient OH overshoot of the maximum steady concentration observed in experiment, while $HO_2$ radical increases more than the maximum steady concentration with increasing SDR. In addition, it is seen that NO and $NO_2$ are not sensitive to the unsteady variation of SDR.

• Journal of the Korean Society of Combustion
• /
• v.9 no.2
• /
• pp.10-17
• /
• 2004

A two-dimensional direct numerical simulation is performed to investigate the flame structure of $CH_4/N_2-Air$ counterflow nonpremixed flame interacting with a single vortex. The detailed transport properties and a modified 16-step augmented reduced mechanism based on Miller and Bowman#s detailed reaction mechanism are adopted in this calculation. To quantify the strain on flame induced by a vortex, a scalar dissipation rate (SDR) is introduced. The results show that fuel-side and air-side vortex cause an unsteady extinction. In this case, the flame interacting with a vortex is extinguished at much larger SDR than steady flame. It is also found that air-side vortex extinguishes a flame more rapidly than fuel-side vortex. The unsteady effect induced by flame-vortex interaction does not lead to a transient OH overshoot of the maximum steady concentration observed in experiment, while $HO_2$ radical increases more than the maximum steady concentration with increasing SDR. In addition, it is seen that NO and $NO_2$ are not sensitive to the unsteady variation of SDR.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2003.04a
• /
• pp.184-187
• /
• 2003

Complete prediction of second order permeability tensor for three dimensional circular braided preform is critical to understand the resin transfer molding process of composites. The permeability can be predicted by considering resin flow through the multi-axial fiber structure. In this study, permeability tensor for a 3-D circular braided preform is calculated by solving a boundary problem of a periodic unit cell. Flow field through the unit cell is obtained by using a 3-D finite volume method (FVM) and Darcy's law is utilized to obtain permeability tensor. Flow analysis for two cases that a fiber tow is regarded as impermeable solid and permeable porous medium is carried out respectively. It is found that the flow within the intra-tow region of the braided preform is negligible if inter-tow porosity is relatively high but the flow through the tow must be considered when the porosity is low. To avoid checkerboard pressure field and improve the efficiency of numerical computation, a new interpolation function for velocity variation is proposed on the basis of analytic solutions. Permeability of the braided preform is measured through a radial flow experiment and compared with the permeability predicted numerically.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.195-201
• /
• 2008

In the jet engines on the aircrafts cruising at high altitude over 20 km and subsonic speed, the Reynolds number in terms of the compressor blades becomes very low. In such an operating condition with low Reynolds number, it is widely reported that total pressure loss of the air flow through the compressor cascades increases dramatically due to separation of the boundary layer and the secondary-flow. But the detail of flow mechanisms causes the total pressure loss has not been fully understood yet. In the present study, two series of numerical investigations were conducted to study the effects of Reynolds number on the aerodynamic characteristics of compressor cascades. At first, the incompressible flow fields in the two-dimensional compressor cascade composed of C4 airfoils were numerically simulated with various values of Reynolds number. Compared with the corresponding experimental data, the numerically estimated trend of total pressure loss as a function of Reynolds number showed good agreement with that of experiment. From the visualized numerical results, the thickness of boundary layer and wake were found to increase with the decrease of Reynolds number. Especially at very low Reynolds number, the separation of boundary layer and vortex shedding were observed. The other series, as the preparatory investigation, the flow fields in the transonic compressor, NASA Rotor 37, were simulated under the several conditions, which corresponded to the operation at sea level static and at 10 km of altitude with low density and temperature. It was found that, in the case of operation at high altitude, the separation region on the blade surface became lager, and that the radial and reverse flow around the trailing edge become stronger than those under sea level static condition.