• Journal of Korean Port Research
• /
• v.14 no.4
• /
• pp.451-461
• /
• 2000

The design and maintenance of navigation channel and water facilities of an harbor which is located at the mouth of river or at the estuary area are difficult due to the complexity of estuarial water and sediment circulation. Effects of deepening navigable waterways, of changing coastline configurations, or of discharging dredged material to the open sea are necessary to be investigated and predicted in terms of water quality and possible physical changes to the coastal environment. A borad analysis of the transport mechanism in the estuary area was made in terms of sediment property, falling velocity, concentration and flow characteristics. In order to simulate the transport processes, a two-dimensional finite element model is developed, which includes erosion, transport and deposition mechanism of suspended sediments. Galerkin’s weighted residual method is used to solve the transient convection-diffusion equation. The fluid domain is subdivided into a series of triangular elements in which a quadratic approximation is made for suspended sediment concentration. Model could deal with a continuous aggregation by stipulating the settling velocity of the flocs in each element. The model provides suspended sediment concentration, bed shear stress, erosion versus deposition rate and bed profile at the given time step.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.2
• /
• pp.245-252
• /
• 2002

Numerical analysis was conducted to characterize particle deposition on a horizontal rotating disk with thermophorectic effect under laminar flow field. The particle transport mechanisms considered were convection, Brownian diffusion, gravitational settling and thermophoresis. The averaged particle deposition velocities and their radial distributions for the upper surface of the disk were calculated from the particle concentration equation in a Eulerian frame of reference for rotating speeds of 0∼1000rpm and temperature differences of 0∼5K. It was observed from the numerical results that the rotation effect of disk increased the averaged deposition velocities, and enhanced the uniformity of local deposition velocities on the upper surface compared with those of the disk at rest. It was also shown that the heating of the disk with ΔT=5K decreased deposition velocity over a fairly broad range of particle sizes. Finally, an approximate deposition velocity model for the rotating disk was suggested. The comparison of the present numerical results with the results of the approximate model and the available experimental results showed relatively good agreement between them.

• Nuclear Engineering and Technology
• /
• v.56 no.4
• /
• pp.1213-1224
• /
• 2024

The liquid lead-bismuth cooled fast reactor has been in a single-phase, low-pressure, and high-temperature state for a long time during operation. Considering the requirement of calculation efficiency for long-term transient accident calculation, based on a homogeneous hydrodynamic model, one-dimensional heat conduction model, coolant flow and heat transfer model, neutron kinetics model, coolant and material properties model, this study used the fully implicit difference scheme algorithm of the convection-diffusion term to solve the basic conservation equation, to develop the transient analysis program NUSOL-LMR 2.0 for the lead-bismuth fast reactor system. The steady-state and typical design basis accidents (including reactivity introduction, loss of flow caused by main pump idling, excessive cooling, and plant power outage accidents) for the ABR have been analyzed. The results are compared with the international system analysis software ATHENA. The results indicate that the developed program can stably, accurately, and efficiently predict the transient accident response and safety characteristics of the lead-bismuth fast reactor system.

• Journal of Korea Soil Environment Society
• /
• v.2 no.2
• /
• pp.81-94
• /
• 1997

In many solute transport studies, either flux or resident concentration has been used. Choice of the concentration mode was dependent on the monitoring device in solute displacement experiments. It has been accepted that no priority exists in the selection of concentration mode in the study of solute transport. It would be questionable, however, to accept the equivalency in the solute transport parameters between flux and resident concentrations in structured soils exhibiting preferential movement of solute. In this study, we investigate how they differ in the monitored breakthrough curves (BTCs) and transport parameters for a given boundary and flow condition by performing solute displacement experiments on a number of undisturbed soil columns. Both flux and resident concentrations have been simultaneously obtained by monitoring the effluent and resistance of the horizontally-positioned TDR probes. Two different solute transport models namely, convection-dispersion equation (CDE) and convective lognormal transfer function (CLT) models, were fitted to the observed breakthrough data in order to quantify the difference between two concentration modes. The study reveals that soil columns having relatively high flux densities exhibited great differences in the degree of peak concentration and travel time of peak between flux and resident concentrations. The peak concentration in flux mode was several times higher than that in resident one. Accordingly, the estimated parameters of flux mode differed greatly from those of resident mode and the difference was more pronounced in CDE than CLT model. Especially in CDE model, the parameters of flux mode were much higher than those of resident mode. This was mainly due to the bypassing of solute through soil macropores and failure of the equilibrium CDE model to adequate description of solute transport in studied soils. In the domain of the relationship between the ratio of hydrodynamic dispersion to molecular diffusion and the peclet number, both concentrations fall on a zone of predominant mechanical dispersion. However, it appears that more molecular diffusion contributes to the solute spreading in the matrix region than the macropore region due to the nonliearity present in the pore water velocity and dispersion coefficient relationship.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.2
• /
• pp.181-188
• /
• 2012

The present study on nitrogen-diluted non-premixed counterflow flames with finite burner diameters experimentally investigates the important role of the outer edge flame in flame extinction. Flame stability diagrams mapping the flame extinction response of nitrogen-diluted non-premixed counterflow flames to varying global strain rates in terms of the burner diameter, burner gap, and velocity ratio are explored. There exists a critical nitrogen mole fraction beyond which the flame cannot be sustained, and also the curves of the critical nitrogen mole fraction versus the global strain rate have C-shapes in terms of burner diameter, burner gap, and velocity ratio. In flames with sufficiently high strain rates, the curves of the critical nitrogen mole fractions versus global strain rate collapse into one curve, and the flames can have the 1-D flame response of typical diffusion flames. Three flame extinction modes are identified: flame extinctions through the shrinkage of the outer edge flame with and without an oscillation of the outer edge flame prior to the extinction and flame extinction through a flame hole at the flame center. The measured flame surface temperature and a numerical evaluation of the fractional contribution of each term in the energy equation show that the radial conductive heat loss at the flame edge destabilizes the outer edge flame, and the conductive and convection heat addition to the outer edge from the trailing diffusion flame stabilizes the outer edge flame. The radial conductive heat loss at the flame edge is the dominant extinction mechanism acting through the shrinkage of the outer edge flame.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.5
• /
• pp.1319-1332
• /
• 1995

Numerical analysis was performed to characterize the particle deposition behavior on a horizontal free-standing wafer with thermophoretic effect under the turbulent flow field. A low Reynolds number k-.epsilon. turbulence model was used to analyze the turbulent flow field around the wafer, and the temperature field for the calculation of the thermophoretic effect was predicted from the energy equation introducing the eddy diffusivity concept. The deposition mechanisms considered were convection, diffusion, sedimentation, turbulence and thermophoresis. For both the upper and lower surfaces of the wafer, the averaged particle deposition velocities and their radial distributions were calculated and compared with the laminar flow results and available experimental data. It was shown by the calculated averaged particle deposition velocities on the upper surface of the wafer that the deposition-free zone, where the deposition velocite is lower than 10$^{-5}$ cm/s, exists between 0.096 .mu.m and 1.6 .mu.m through the influence of thermophoresis with positive temperature difference of 10 K between the wafer and the ambient air. As for the calsulated local deposition velocities, for small particle sizes d$_{p}$<0.05 .mu.m, the deposition velocity is higher at the center of the wafer than at the wafer edge, whereas for particle size of d$_{p}$ = 2.0 .mu.m the deposition takes place mainly on the inside area of the wafer. Finally, an approximate model for calculating the deposition velocities was recommended and the calculated deposition velocity results were compared with the present numerical solutions, those of Schmidt et al.'s model and the experimental data of Opiolka et al.. It is shown by the comparison that the results of the recommended model agree better with the numerical solutions and Opiolka et al.'s data than those of Schmidt's simple model.