• Magazine of the Korean Society of Agricultural Engineers
• /
• v.34 no.1
• /
• pp.87-99
• /
• 1992

A mathematical model, UNSATR which predicts the seepage flow through the body of dike especially under the tidal fluctuation has been developed. This model has been revised from UNSAT2 model which was developed on the basis of the saturated-unsaturated theory by Neuman. UNSATR has been verified and applied to the hydraulic model in order to estimated the seepage quantity, the formation of free water surface etc. The results lead to the following conclusions : 1. Seepage rates between the mathematical model and hydraulic model experiment are very similar to each other both in constant and transient water level conditions. 2. The lapsed time to be steady state of the free water surface becomes late as the tidal levels are relatively low mainly due to the seepage flow from the unsaturated zone of the body of dike. 3. Under the transient state of water levels, owing to the flow from the unsaturated domain, streamlines crossing to the free water surface are found and time lag during a falling tide may allow the free water surface inside the body of dike to stand at a high level than the outside water level. 4. The utility and validity of UNSATR model are convinced when the analyses on seepage problems through the porous embankment of the soil structures on the conditions of the steady and unsteady states are carried out.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2744-2749
• /
• 2008

Numerical calculations are carried out for the natural convection induced by temperature difference between a cold outer square cylinder and a hot inner circular cylinder for Rayleigh number of $Ra=10^7$. This study investigates the effect of the inner cylinder location on the heat transfer and fluid flow. The location of inner circular cylinder ($\delta$) is changed vertically along the center-line of square enclosure. The natural convection bifurcates from unsteady to steady state according to $\delta$. Two critical positions of ${\delta}_{C,L}$ and ${\delta}_{C,U}$ as a lower bound and an upper bound are ${\delta}_{C,L}=0.05$ and ${\delta}_{C,U}=0.18$, respectively. Within the defined bounds, the thermal and flow fields are steady state. When the inner cylinder locates at ${\delta}{\geq}{\delta}_{C,U}$, the space between the upper surface of inner cylinder and the top surface of the enclosure forms a relatively shallow layer where the natural convection characterized as the pure Rayleigh-Benard convection forms alternately the upwelling and downwelling plums, as a result that a series of cells known as Benard cells is derived.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.62-68
• /
• 2004

A comprehensive numerical analysis has been carried out for both non-reacting and reacting flows in a scramjet engine combustor with and without a cavity. The theoretical formulation treats the complete conservation equations of chemically reacting flows with finite-rate chemistry of hydrogen-air. Turbulence closure is achieved by means of a k-$\omega$ two-equation model. The governing equations are discretized using a MUSCL-type TVD scheme, and temporally integrated by a second-order accurate implicit scheme. Transverse injection of hydrogen is considered over a broad range of injection pressure. The corresponding equivalence ratio of the overall fuel/air mixture ranges from 0.167 to 0.50. The work features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous studies. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the .underlying physical mechanisms. Much of the flow unsteadiness is related not only to the cavity, but also to the intrinsic unsteadiness in the flow-field. The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The roles of the cavity, injection pressure, and heat release in determining the flow dynamics are examined systematically.

• Journal of the Korean GEO-environmental Society
• /
• v.2 no.1
• /
• pp.81-89
• /
• 2001

Soil-bentonite mixture is often used for barrier wall to prevent seawater intrusion. In this study, the effect of seawater on the permeability of soil-bentonite mixture is examined, and the effect of permeability change on the seawater intrusion is investigated. Seawater intrusion in coastal areas was modeled using a finite element method. Seawater intrusion in the seawater-contaminated zone was determined by considering the hydraulic conductivity changes using the residual flow procedure (RFP) in the simulation model. Steady state and unsteady state conditions with variations in ground water levels in an inland area were investigated. The interface between fresh water and seawater, found by the proposed method, was located lower at the seawater side and the level at the fresh water side is higher than those by conventional methods.

• Wind and Structures
• /
• v.29 no.4
• /
• pp.225-234
• /
• 2019

Conductor galloping is a common disaster for the transmission lines. Among the existing analytical methods, the wind tunnel test is highlighted as the most effective approach to obtain the aerodynamic coefficients. In this paper, the aerodynamic coefficients of 12 conductor models covered with the crescent-shaped ice, which were fabricated considering the surface roughness of the iced conductor, were obtained based on the wind tunnel test. The influence of the Reynolds number and the shape parameter ${\beta}$, defined as the ratio of ice thickness to the diameter, were investigated. In addition, the effect of surface roughness of the iced conductor was discussed. Subsequently, unsteady areas of conductor galloping were calculated according to the Den Hartog criterion and the Nigol criterion. The results indicate that the aerodynamic coefficients of iced conductors change sharply at the attack angles of $20^{\circ}$ and $170^{\circ}$ with the increase of ${\beta}$. The surface roughness of iced conductors changed the range of attack angle, which was influenced by the increase of the Reynolds number. The experimental results can provide insights for preventing and controlling galloping.

• Nuclear Engineering and Technology
• /
• v.51 no.7
• /
• pp.1842-1852
• /
• 2019

The AP1000 reactor coolant pump is a vertical shielded-mixed flow pump, is the most important coolant power supply and energy exchange equipment in nuclear reactor primary circuit system, whose steadystate and transient performance affect the safety of the whole nuclear island. Moreover, safety demonstration of reactor coolant pump is the most important step to judge whether it can be practiced, among which software simulation is the first step of theoretical verification. This paper mainly introduces the fluid-solid coupling simulation method applied to reactor coolant pump, studying the feasibility of simulation results based on workbench fluid-solid coupling technology. The study found that: for the unsteady calculations of the pure liquid media, the average head of the reactor coolant pump with bidirectional fluid-solid coupling decreases to a certain extent. And the coupling result is closer to the real experimental value. The large stress and deformation of rotor under different flow conditions are mainly distributed on impeller and idler, and the stress concentration mainly occurs at the junction of front cover plate and blade outlet. Among the factors that affect the dynamic stress change of rotor, the pressure load takes a dominant position.

• Ocean Systems Engineering
• /
• v.4 no.3
• /
• pp.185-200
• /
• 2014

The objective of the present simulations is to evaluate the applicability of the standard $k-{\varepsilon}$ turbulence model in engineering practice in the subcritical to supercritical flow regimes. Two-dimensional numerical simulations of flow around a circular cylinder at $Re=1{\times}10^5$, $5{\times}10^5$ and $1{\times}10^6$, had been performed using Unsteady Reynolds-Averaged Navier Stokes (URANS) equations with the standard $k-{\varepsilon}$ turbulence model. Solution verification had been studied by evaluating grid and time step size convergence. For each Reynolds number, several meshes with different grid and time step size resolutions were chosen to calculate the hydrodynamic quantities such as the time-averaged drag coefficient, root-mean square value of lift coefficient, Strouhal number, the coefficient of pressure on the downstream point of the cylinder, the separation angle. By comparing the values of these quantities of adjacent grid or time step size resolutions, convergence study has been performed. Solution validation is obtained by comparing the converged results with published numerical and experimental data. The deviations of the values of present simulated quantities from those corresponding experimental data become smaller as Reynolds numbers increases from $1{\times}10^5$ to $1{\times}10^6$. This may show that the standard $k-{\varepsilon}$ model with enhanced wall treatment appears to be applicable for higher Reynolds number turbulence flow.

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

When the heat release and acoustic pressure fluctuations are generated in the combustor by irregular combustion, these fluctuations affect the mass flow rate of the propellants injected through the injectors. Also, the variations of the mass flow rate by these fluctuations again bring about irregular combustion and furthermore that is related with combustion instability. Therefore, it is very important to identify the mass variation for the pressure fluctuation on the injector and to investigate its transfer function. So, we first have studied quantifying the variation of mass flow rate generated in simplex swirl injector by injection pressure fluctuation. To acquire the transient mass flow rate in orifice with time, we have tried to measure of the flow axial velocity and liquid film thickness in orifice. The axial velocity is acquired through theoretical approach after measuring the pressure in orifice and the flow area in the orifice is measured by electric conductance method. As results, mass flow rate calculated by axial velocity and liquid film thickness measuring in orifice accorded with mass flow rate acquired by direct measuring method in the small error range within 1 percents in steady state and within 6 percents as average mass flow rate in pulsated state. Hence this method can be used to measure the mass flow rate not only in steady state but also in unsteady state because the mass flow rate in the orifice can acquire with time and this method shows very high accuracy based on the experimental results.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.4 no.1
• /
• pp.46-56
• /
• 1994

The characteristics of flows, temperatures, concentrations of the boron are numerically studied when uniform axial magnetic fields are applied in the Czechralski crucible. The to governing factors to the flow regimes are buoyancy, thermocapillarity, centrifugal forces, magnetic forces, diffusion coefficient and segregation coefficient of the boron. Since the concentration of the boron is so low that buoyancy effects are negligible, it cannot affect the flow and temperature fields. From the fact that the flow fields are rotationally symmetric, two velocity components in the meridional plane and the circumferential velocity are calculated together with the temperature in the steady state. Based on the known velocity and temperature distributions the unsteady concentration distributions of the boron are calculated. As the strength of the magnetic is increased, the flow velocities are decreased. Circumferential velocities are large near the crucible side-wall and in the region below the rotating crystal. Steep temperatures gradient near the edge of the rotating crystal causes the Marangoni convection. It has been found out that the convection characteristics affects the unsteady transport phenomena of the boron.

• The Journal of Engineering Geology
• /
• v.19 no.2
• /
• pp.205-215
• /
• 2009

The conventional numerical models to analyze flow in subsurface porous media under the transient state usually generate numerical oscillation and unstability due to local flux domain for critical cases such as infiltration into initially dry soil during rainfall period. In this case, it is required refined mesh and small time step, but it decrease efficiency of computation. In this study, numerical unstability in discontinuity domain is removed by applying particle tracking algorithm to simulate unsteady subsurface flow with inflow boundary condition. Finally the hybrid LE FEM improving numerical stability is proposed. The hypothetical domains with unsteady uniform and nonuniform flow field were used to demonstrated algorithm verification. In comparison with analytic solution, we obtained reasonable results and conducted simulation of hypothetical 3-D recharge/pumping area. The proposed algorithm can simulate saturated/unsaturated porous media with more practical problems and will greatly contribute to accuracy and stability of numerical computation.