• Ocean and Polar Research
• /
• v.25 no.spc3
• /
• pp.361-371
• /
• 2003

The effect of bottom friction on the steady wind-driven circulation in the Yellow Sea and the East China Sea (YSECS) has been studied using a two-dimensional numerical model with and without tidal forcing. Upwind flow experiment in YSECS has also been carried out with a schematic time variation in the wind field. The surface water setup and circulation pattern due to steady wind forcing are found to be very sensitive to the bottom friction. When the effects of tidal currents are neglected, the overall current velocities are overestimated and eddies of various sizes appear, upwind flow is formed within the deep trough of the Yellow Sea, forming a part of the topographic gyre on the side of Korea. When tidal forcing is taken into account, the wind-induced surface elevations are smoothed out due to the strong tide-induced bottom friction, which is aligned almost normal to the wind stresses; weak upwind flow is farmed in the deep trough of the Yellow Sea, west and south of Jeju. Calculation with wind forcing only through a parameterized linear bottom friction produces almost same results from the calculation with $M_2$ tidal forcing and wind forcing using a quadratic bottom friction, supporting Hunter (1975)'s linearization of bottom friction which includes the effect of tidal current, can be applied to the simulation of wind-driven circulation in YSECS. The results show that steady wind forcing is not a dominant factor to the winter-time upwind flow in YSECS. Upwind flow experiment which considers the relaxation of pressure gradient (Huesh et al. 1986) shows that 1) a downwind flow is dominant over the whole YSECS when the northerly wind reaches a maximum speed; 2) a trend of upwind flow near the trough is found during relaxation when the wind abates; 3) a northward flow dominates over the YSECS after the wind stops. The results also show that the upwind flow in the trough of Yellow Sea is forced by a wind-induced longitudinal surface elevation gradient.

• Journal of Korea Water Resources Association
• /
• v.52 no.2
• /
• pp.149-161
• /
• 2019

In this study, the effect of the fluid flow characteristics on the length distribution of the fracture elements composing the fracture network is analyzed numerically using the 3D fracture crack network model. The truncated power-law distribution is applied to generate the length distribution of the fracture elements and the simulations of fluid flow are carried out with the exponent ${\beta}_l$ from 1.0 to 6.0. As a result of simulations, when the exponent ${\beta}_l$ increases, the length distribution of the fracture elements gradually decreases, and the connectivity between the fracture elements affecting the permeability of the fracture network becomes weak. When we analyzed the distributions of flow rate calculated at each fracture element with the exponent ${\beta}_l$, the mean flow rate at ${\beta}_l=1.0$ was estimated to be about 447 times larger than that at ${\beta}_l=6.0$ and for the flow calculated at the outflow boundary of the fracture network, the case of ${\beta}_l=1.0$ was estimated to be 6,440 times larger than that of ${\beta}_l=6.0$.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.17 no.4
• /
• pp.384-398
• /
• 2015

In this paper, the high-resolution Weather Research and Forecasting/Noah-MultiParameterization (WRF/Noah-MP) modeling system is configured for the Cheongmicheon Farmland site in Korea (CFK), and its performance in land and atmospheric simulation is evaluated using the observed data at CFK during the 2014 special observation period (21 August-10 September). In order to explore the usefulness of turning on Noah-MP dynamic vegetation in midterm simulations of surface and atmospheric variables, two numerical experiments are conducted without dynamic vegetation and with dynamic vegetation (referred to as CTL and DVG experiments, respectively). The main results are as following. 1) CTL showed a tendency of overestimating daytime net shortwave radiation, thereby surface heat fluxes and Bowen ratio. The CTL experiment showed reasonable magnitudes and timing of air temperature at 2 m and 10 m; especially the small error in simulating minimum air temperature showed high potential for predicting frost and leaf wetness duration. The CTL experiment overestimated 10-m wind and precipitation, but the beginning and ending time of precipitation were well captured. 2) When the dynamic vegetation was turned on, the WRF/Noah-MP system showed more realistic values of leaf area index (LAI), net shortwave radiation, surface heat fluxes, Bowen ratio, air temperature, wind and precipitation. The DVG experiment, where LAI is a prognostic variable, produced larger LAI than CTL, and the larger LAI showed better agreement with the observed. The simulated Bowen ratio got closer to the observed ratio, indicating reasonable surface energy partition. The DVG experiment showed patterns similar to CTL, with differences for maximum air temperature. Both experiments showed faster rising of 10-m air temperature during the morning growth hours, presumably due to the rapid growth of daytime mixed layers in the Yonsei University (YSU) boundary layer scheme. The DVG experiment decreased errors in simulating 10-m wind and precipitation. 3) As horizontal resolution increases, the models did not show practical improvement in simulation performance for surface fluxes, air temperature, wind and precipitation, and required three-dimensional observation for more agricultural land spots as well as consistency in model topography and land cover data.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.28 no.3
• /
• pp.279-293
• /
• 1995

The most obvious and easily recognizable sources of potential water pollution are point sources such as domestic and industrial wastes. But recently, the potential effects of nonpoint sources on water quality have been increased apparently. In order to evaluate the characteristics and the effects of nonpoint sources on water quality, this study was performed in Suyeong Bay from May, 1992 to July, 1992. The depth-averaged 2-dimensional numerical model, which consists of the hydrodynamic model and the diffusion model was applied to simulate the water quality in Suyeong Bay. When flowrate was $65.736m^3/s,$ the concentration of pollutants (COD, TSS and VSS) at Oncheon stream (Sebeong bridge) during second flush were very high as much as 121.4mg/l of COD, 1148.0mg/l of TSS and 262.0mg/1 of VSS. When flowrate was 4.686m^3/s, the concentration of pollutants $(TIN,\;NH_4\;^+-\;N,\;NO_2\;^--N\;and\;PO_4\;^{3-}-P)$ during the first flush were very high as much as 20.306mg/1 of TIN, 14.154mg/1 of $NH_4\;^+-N$, 9.571mg/l of $NO_2\;^--N$ and l.785mg/l of $PO_2\;^{3-}-P$ As results of the hydrodynamic model simulation, the computed maximum velocity of tidal currents in Suyeong Bay was 0.3m/s and their direction was clockwise flow for ebb tide and counter clockwise flow for Hood tide. Four different methods were applied for the diffusion simulation in Suyeong Bay. There were the effects for the water quality due to point loads, annual nonpoint loads and nonpoint loads during the wet weather and the investigation period, respectively. The efforts of annual nonpoint loads and nonpoint loads during the wet weather seem to be slightly deteriorated in comparison with the effects of point loads. However, the bay was significantly polluted by the nonpoint loads during the investigation period. In this case, COD and SS concentrations ranged 2.0-30.0mg/l, 7.0- 200.0mg/l in ebb tide, respectively. From these results, it can be emphasized that the large amount of pollutants caused by nonpoint sources during the wet weather were discharged into the bay, and affected significantly to both the water quality and the marine ecosystem. Therefore, it is necessary to consider the loadings of nonpoint pollutants to plan wastewater treatment plant.

• International Journal of Highway Engineering
• /
• v.11 no.1
• /
• pp.165-175
• /
• 2009

With the current move towards adopting mechanistic-empirical concepts in the design of pavement structures, state-of-the-art mechanistic analysis methodologies are needed to determine accurate pavement responses, such as stress, strain, and deformation. Previous laboratory studies of pavement foundation geomaterials, i.e., unbound granular materials used in base/subbase layers and fine-grained soils of a prepared subgrade, have shown that the resilient responses followed by nonlinear, stress-dependent behavior under repeated wheel loading. This nonlinear behavior is commonly characterized by stress-dependent resilient modulus material models that need to be incorporated into finite element (FE) based mechanistic pavement analysis methods to predict more realistically predict pavement responses for a mechanistic pavement analysis. Developed user material subroutine using aforementioned resilient model with nonlinear solution technique and convergence scheme with proven performance were successfully employed in general-purpose FE program, ABAQUS. This numerical analysis was investigated in predicted critical responses and domain selection with specific mesh generation was implemented to evaluate better prediction of pavement responses. Results obtained from both axisymmetric and three-dimensional (3D) nonlinear FE analyses were compared and remarkable findings were described for nonlinear FE analysis. The UMAT subroutine performance was also validated with the instrumented full scale pavement test section study results from the Federal Aviation Administration's National Airport Pavement Test Facility (FAA's NAPTF).

• Journal of the Korean Institute of Gas
• /
• v.7 no.4 s.21
• /
• pp.44-52
• /
• 2003

In this work dynamic heat transfer in a CPFS (cable penetration fire stop) system built in the firewall of nuclear power plants is three-dimensionally investigated to develop a test-simulator that can be used to verify effectiveness of the sealant. Dynamic heat transfer in the fire stop system is formulated in a parabolic PDE (partial differential equation) subjected to a set of initial and boundary conditions. First, the PDE model is divided into two parts; one corresponding to heat transfer in the axial direction and the other corresponding to heat transfer on the vertical planes. The first PDE is converted to a series of ODEs (ordinary differential equations) at finite discrete axial points for applying the numerical method of SOR (successive over-relaxation) to the problem. The ODEs are solved by using an ODE solver In such manner, the axial heat flux can be calculated at least at the finite discrete points. After that, all the planes are separated into finite elements, where the time and spatial functions are assumed to be of orthogonal collocation state at each element. The initial condition of each finite element can be obtained from the above solution. The heat fluxes on the vertical planes are calculated by the Galerkin FEM (finite element method). The CPFS system was modeled, simulated, and analyzed here. The simulation results were illustrated in three-dimensional graphics. Through simulation, it was shown clearly that the temperature distribution was influenced very much by the number, position, and temperature of the cable stream, and that dynamic heat transfer through the cable stream was one of the most dominant factors, and that the feature of heat conduction could be understood as an unsteady-state process.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.36 no.3
• /
• pp.185-192
• /
• 2023

In this study, fluctuating wind velocity for time history analysis is simulated by a single variate, single-dimensional random process using the KBC2022 spectrum about across-wind direction. This study analyzed and obtained the inelastic dynamic response for structures modeled as a single-degree-of-freedom system. It is assumed that the wind response is excellent in the primary mode, the change in vibration owing to plasticization is minor, along-wind vibration and across-wind vibration are independent, and the effect of torsional vibration is small. The numerical results, obtained by the Newmark-𝛽 method, shows the time-history responses and trends of maximum displacements. As a result of analyzing the inelastic dynamic response of the structure with the second stiffness ratio(𝛼) and yield displacement ratio (𝛽) as variables, it is identified that as the yield displacement ratio (𝛽) increases when the second stiffness ratio is constant, the maximum displacement ratio decreases, then reaches a minimum value, and then increases. When the stiffness ratio is greater than 0.5, there is a yield point ratio at which the maximum displacement ratio is less than 1, indicating that the maximum deformation is reduced compared to the elastically designed building even if the inelastic behavior is permitted in the inelastic wind design.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.25 no.2
• /
• pp.63-75
• /
• 2013

The Inner Port Phase 2 area of the Pyeongtaek-Dangjin Port is enclosed by a total of three permeable sea-walls, and the disposal site to the east of the Inner Port Phase 2 is also enclosed by two permeable sea-walls. The maximum tidal range measured in the Inner Port Phase 2 and in the disposal site in May 2010 is 4.70 and 2.32 m, respectively. It reaches up to 54 and 27%, respectively of 8.74 m measured simultaneously in the exterior. Regression formulas between the difference of hydraulic head and the rate of interior water volume change, are induced. A three-dimensional numerical hydrodynamic model for the Asan Bay is constructed incorporating a module to compute water discharge through the permeable sea-walls at each computation time step by employing the formulas. Hydrodynamics for the period from 13th to 27th May, 2010 is simulated by driving forces of real-time reconstructed tide with major five constituents($M_2$, $S_2$, $K_1$, $O_1$ and $N_2$) and freshwater discharges from Asan, Sapkyo, Namyang and Seokmoon Sea dikes. The skill scores of modeled mean high waters, mean sea levels and mean low waters are excellent to be 96 to 100% in the interior of permeable sea-walls. Compared with the results of simulation to obstruct the flow through the permeable sea-walls, the maximum current speed increases by 0.05 to 0.10 m/s along the main channel and by 0.1 to 0.2 m/s locally in the exterior of the Outer Sea-wall of Inner Port. The maximum bottom shear stress is also intensified by 0.1 to 0.4 $N/m^2$ in the main channel and by more than 0.4 $N/m^2$ locally around the arched Outer Sea-wall. The module developed to compute the flow through impermeable seawalls can be practically applied to simulate and predict the advection and dispersion of materials, the erosion or deposion of sediments, and the local scouring around coastal structures where large-scale permeable sea-walls are maintained.

• The Journal of Engineering Geology
• /
• v.27 no.1
• /
• pp.81-90
• /
• 2017

Preceding displacements in tunnel are difficult to predict since the measurements of displacements after excavation can not be performed immediately. In the present study, The longitudinal displacements which can be measured immediately after excavation are used to predict the crown settlements occurring before excavation only if fault is located at the tunnel crown. Three-dimensional finite element analysis was conducted using 28 numerical models with various fault attitudes to analyze the correlation between the longitudinal displacements on tunnel face and preceding crown settlements. The results, $L_{face}/C$ ratio show 2~12% in the drives with dip models and 2~13% in the drives against dip models individually. In addition, each model has a certain $L_{face}/C$ ratio. The result of the regression analysis show that the coefficient of determination is over 0.8 in most models. Therefore, crown settlements occurring before excavation can be predicted by analyzing the longitudinal displacements occurring on tunnel faces.

• Korean Journal of Remote Sensing
• /
• v.34 no.5
• /
• pp.749-762
• /
• 2018

In this study, the changes of the sunshine environment due to the construction of buildings were analyzed by season. Using a geographic information system (GIS) data, the topography and the buildings were constructed around Pukyong National University (PKNU) in Busan. The numerical model was performed for a week in each season, before and after the construction of buildings. Even before the construction of the high-rise building complex, the area of sunshine block is wide due to the apartment complex located in the southeast of the PKNU campus and the mountain in the south. After the construction, the sunshine-blocked area became wider after the sunrise and before the sunset. The area of sunshine block after 1 hour at sunrise increased by 1.60%, 1.50%, 1.58% and 1.36% in the vernal equinox, summer solstice, autumnal equinox and winter solstice, respectively. The building complex in the east (south) of the PKNU campus made shadow in more than 1,000 m (750 m) toward the west (east) just before the sunrise (sunset). Especially, the sunshine duration in PKNU campus decreased by 46.61%, 22.75%, 58.56%, and 11.31% in the vernal equinox, summer solstice, autumnal equinox and winter solstice, respectively. The analysis of the sunshine duration for a dormitory building showed that the construction of the building complex in the south of the PKNU campus reduced the sunshine area of the western (southern) wall of the dormitory by 30.91% (49.45%) for a winter week.