• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.2
• /
• pp.495-506
• /
• 2013

This study developed the Time-split Mixing Model (TMM) which can represent the pollutant mixing process on a three-dimensional open channel through constructing the conceptual model based on Taylor's assumption (1954) that the shear flow dispersion is the result of combination of shear advection and diffusion by turbulence. The developed model splits the 2-D mixing process into longitudinal mixing and transverse mixing, and it represents the 2-D advection-dispersion by the repetitive calculation of concentration separation by the vertical non-uniformity of flow velocity and then vertical mixing by turbulent diffusion sequentially. The simulation results indicated that the proposed model explains the effect of concentration overlapping by boundary walls, and the simulated concentration was in good agreement with the analytical solution of the 2-D advection-dispersion equation in Taylor period (Chatwin, 1970). The proposed model could explain the correlation between hydraulic factors and the dispersion coefficient to provide the physical insight about the dispersion behavior. The longitudinal dispersion coefficient calculated by the TMM varied with the mixing time unlike the constant value suggested by Elder (1959), whereas the transverse dispersion coefficient was similar with the coefficient evaluated by experiments of Sayre and Chang (1968), Fischer et al. (1979).

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.240-240
• /
• 2015

As Simpson et al. (1990) emphasized the importance of the straining process in the stratification and mixing in the estuarine circulation process, various researches have investigated on the relative contribution of each process to the overall potential energy anomaly dynamics. However, many numerical works have done only for two dimensional modeling along channel or the short distance cross sectional three dimensional simulations as Burchard et al. (2008) and the estuarine channel was not simulated so far. But, in the study on the physics of shallow coastal seas, spatial dimension in the three dimensional way affects significantly on results of a particular numerical model. Therefore, the comparison of two and three dimensional models is important to understand the real physics of mixing and stratification in an estuary. Also, as Geyer and MacCready (2013) pointed out that the lateral process seems to be important in determining the periodic stratifications, to study such process the three dimensional modeling must be required. The present study uses a numerical model to show the signification roles of each term of the time-dependent dynamic equation for the potential energy anomaly (PEA) in controlling along and lateral channel flows and different stratification structures. Moreover, we present the relationships between the ${\Phi}$-advection, the depth mean straining, vertical mixing and vertical advection can explain well how water level, salinity distribution and across velocity 2D model are slightly different from 3D.

• Proceedings of the KSRS Conference
• /
• v.1
• /
• pp.21-24
• /
• 2006

Processes controlling the interannual variation of mixed layer temperature (MLT) averaged over the NINO3 domain ($150-90^{\circ}W$, $5^{\circ}N-5^{\circ}S$) are studied using an ocean data assimilation product that covers the period of 1993 to 2003. Advective tendencies are estimated here as the temperature fluxes through the domain's boundaries, with the boundary temperature referenced to the domain-averaged temperature to remove the dependence on temperature scale. The overall balance is such that surface heat flux opposes the MLT change but horizontal advection and subsurface processes assist the change. The zonal advective tendency is caused primarily by large-scale advection of warm-pool water through the western boundary of the domain. The meridional advective tendency is contributed mostly by Ekman current advecting large-scale temperature anomalies though the southern boundary of the domain. Unlike many previous studies, we explicitly evaluate the subsurface processes that consist of vertical mixing and entrainment. In particular, a rigorous method to estimate entrainment allows an exact budget closure. The vertical mixing across the mixed layer (ML) base has a contribution in phase with the MLT change. The entrainment tendency due to temporal change in ML depth is negligible comparing to other subsurface processes. The entrainment tendency by vertical advection across the ML base is dominated by large-scale changes in wind-driven upwelling and temperature of upwelling water. Tropical instability waves (TIWs) result in smaller-scale vertical advection that warms the domain during La Ni? cooling events. When the advective tendencies are evaluated by spatially averaging the conventional local advective tendencies of temperature, the apparent effects of currents with spatial scales smaller than the domain (such as TIWs) become very important as they redistribute heat within the NINO3 domain. However, such internal redistribution of heat does not represent external processes that control the domain-averaged MLT.

• Proceedings of KOSOMES biannual meeting
• /
• 2001.10a
• /
• pp.137-144
• /
• 2001

A Plankton ecosystem model was developed to investigate effects of hydrodynamic processes including advection and diffusion on size-structured phytoplankton dynamics in the mesohaline zone of the York River estuarine system, Virginia, USA. The model included 12 state variables representing the distribution of carbon and nutrients in the surface mixed layer. Groupings of autotrophs and heterotrophs were based on cell site and ecological hierarchy Forcing functions included incident radiation, temperature, wind stress, mean How and tide which includes advective transport and turbulent mixing. The ecosystem model was developed in FORTRAN using differential equations that were solved using the 4th order Runge-Kutta technique. The model showed that microphytoplankton blooms during winter-spring resulted from a combination of vertical advection and diffusion of phytoplankton cells rather than in-situ production in the lower York River estuary.

• Journal of Wetlands Research
• /
• v.15 no.1
• /
• pp.59-69
• /
• 2013

The numerical simulations were performed to analyze the advection-diffusion processes of dredging-induced turbidity plume using RMA2 and RMA4 models in Bunam reservoir, Seosan, Chungnam. Field survey was also performed to measure the turbidity using the multi water quality monitoring system (YSI6600EDS). In the field survey, the vertical and horizontal distributions of the turbidity were measured during the dredging operation in Bunam reservoir. RMA2 model was used to simulate the velocity distributions in both the whole domain and the 2nd part of Bunam reservoir. RMA4 model was also used to simulate the concentration distribution in only the 2nd part of Bunam reservoir, where the dredging work were conducted. The comparison of the simulation results with the field data for the advection-diffusion of the turbidity plume using the concentration ratio concepts shows that the numerical model can be applied to analyze the environmental impact of dredging works.

• 한국방재학회:학술대회논문집
• /
• 2008.02a
• /
• pp.233-236
• /
• 2008

To analyze the saline intrusion in the place, such as an estuary, the three-dimensional numerical model is developed. In this study, the advection terms of the governing equations are discretized by upwind scheme. By using an explicit scheme for the longitudinal direction and an implicit scheme for the vertical direction, the numerical model is free from the restriction of temporal step size caused by a relatively small grid ratio. The equation of state is used to consider the density, and the scalar transport equation for salinity is employed the third order TVD to scheme to prevent unphysical oscillation near discontinuity. In order to verify saline intrusion, the numerical model is conducted to compare the previous model in the lock exchange. The present model generally show a good agreement with the previous one.

• Water for future
• /
• v.29 no.4
• /
• pp.187-198
• /
• 1996

A numerical model for unsteady dispersion of horizontal line source in turbulent shear flow is developed. A fractional step finite difference method is used which splits the unsteady two-dimensional advective diffusion equation into the longitudinal advection and the vertical diffusion equations, and solves them alternately for half time intervals by the Holly-Preissmann scheme and the Crank-Nicholson scheme, respectively. The developed numerical model is verified using a semi-analytic solution for steady dispersion in turbulent shear flow. Dispersion of an instantaneous plane source in turbulent shear flow is analyzed using the model. The degree of mixing at the same dimensionless time is almost the same regardless of the friction factor, and the travel distance required to reach a certain degree of mixing is inversely proportional to the square root of the friction factor.

• Journal of Environmental Science International
• /
• v.10 no.5
• /
• pp.329-336
• /
• 2001

Characgcteristics of nocturnal boundary layer(NBL) were analyzed by the upper-air observations data using with the airsonde and pilot balloons from 1994 to 1999 in Kyungpook province. The automate weather boundary layer can become stably stratified when the surface is cooler than the air. Stable nocturnal boundary layer height were estimated from the top of surface stable layer where the vertical gradient of temperature and mixing ratio tend to zero or negative. The depth of the stable nocturnal boundary layer depended largely on the thermal effect rather than the wind effect at nighttime. The NBL was more developed on the land than on the coastal region. The stability index (bulk Richardson number) showed that the NBL was stable when the wind was weak and the vertical gradient of the temperature was strong. The heat budget in the NBL was studied by considering the effect of the radiative and the cooled by both the longwave radiative flux and the divergence of the heat flux, while NBL under the cloudy sky the longwave radiative flux played a role of the warming. It was noted that the heat was not conserved in both cases. To complete the heat budget in the NBL the warming/cooling by advection and subsidence must be considered.

• Journal of Environmental Science International
• /
• v.11 no.5
• /
• pp.419-436
• /
• 2002

The current paper reports on the enhancement of O$_3$, CO, NO$_2$, and aerosols during the Asian dust event that occurred over Korea on 1 May 1999. To confirm the origin and net flux of the O$_3$, CO, NO$_2$, and aerosols, the meteorological parameters of the weather conditions were investigated using Mesoscale Meteorological Model 5(MM5) and the TOMS total ozone and aerosol index, the back trajectory was identified using the Hybrid Single-Particle Lagrangian Integrated Trajectory Model(HYSPLIT), and the ozone and ozone precursor concentrations were determined using the Urban Ashed Model(UAM). In the presence of sufficiently large concentrations of NO$\sub$x/, the oxidation of CO led to O$_3$ formation with OH, HO$_2$, NO, and NO$_2$ acting as catalysts. The sudden enhancement of O$_3$, CO, NO$_2$ and aerosols was also found to be associated with a deepening cut-off low connected with a surface cyclone and surface anticyclone located to the south of Korea during the Asian dust event. The wave pattern of the upper trough/cut-off low and total ozone level remained stationary when they came into contact with a surface cyclone during the Asian dust event. A typical example of a stratosphere-troposphere exchange(STE) of ozone was demonstrated by tropopause folding due to the jet stream. As such, the secondary maxima of ozone above 80 ppbv that occurred at night in Busan, Korea on 1 May 2001 were considered to result from vertical mixing and advection from a free troposphere-boundary layer exchange in connection with an STE in the upper troposphere. Whereas the sudden enhancement of ozone above 100 ppbv during the day was explained by the catalytic reaction of ozone precursors and transport of ozone from a slow-moving anticyclone area that included a high level of ozone and its precursors coming from China to the south of Korea. The aerosols identified in the free troposphere over Busan, Korea on 1 May 1999 originated from the Taklamakan and Gobi deserts across the Yellow River. In particular, the 1000m profile indicated that the source of the air parcels was from an anticyclone located to the south of Korea. The net flux due to the first invasion of ozone between 0000 LST and 0600 LST on 1 May 1999 agreed with the observed ground-based background concentration of ozone. From 0600 LST to 1200 LST, the net flux of the second invasion of ozone was twice as much as the day before. In this case, a change in the horizontal wind direction may have been responsible for the ozone increase.

• 한국방재학회:학술대회논문집
• /
• 2008.02a
• /
• pp.799-802
• /
• 2008

General behaviors based on hydraulic characteristics of natural streams and channels have been recently analyzed and developed via various numerical models. However in the states of natural hydraulics, an experimental research must be performed simultaneously with the mathematical analysis due to effects of hydraulic properties such as meander, sediment, and so on. In this study based on 2-D advection-dispersion equation, flow and tracer experiments were performed in the S-curved meandering laboratory channel with a rectangular cross-section. The channel was equipped with instrument carriages which was equipped with an auto-traversing system to be used with velocity measuring sensors throughout the depth and breadth of the flow field. To measure concentration distribution of the salt solution was adjusted to that of the flume water by adding methanol and a red dye (KMnO4) was added to aid the visualization of the tracer cloud, the tracer was instantaneously injected into the flow as a full-depth vertical line source by the instantaneous injector and the initial concentration of the tracer was 100,000 mg/l. The secondary current as well as the primary flow pattern was analyzed to investigate the flow distribution in the meandering channels. The velocity distribution of the primary flow for all cases skewed toward the inner bank at the first bend, and was almost symmetric at the crossovers, and then shifted toward the inner bank again at the next alternating bend. Thus, one can clearly notice that the maximum velocity occurs taking the shortest course along the channel, irrespective of the flow conditions. The result of the tracer tests shows that pollutant clouds are spreading following the maximum velocity lines in each cases with various mixing patterns like superposition, separation, and stagnation of pollutant clouds. Flow characteristics in each cases performed in this study can be compared with tracer dispersion characteristics with using evaluation of longitudinal and transverse dispersion coefficients(LDC, TDC). As expected, LDC and TDC in meandering parts have been evaluated with increasing distribution and straight parts have effected to evaluate minimum of LDC and TDC due to symmetric flow patterns and attenuations of secondary flow.