• Journal of Korea Water Resources Association
• /
• v.44 no.6
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.126-134
• /
• 2008

We developed a three dimensional Navier-Stokes code based on the level set method to simulate two phase flows with high density ratio. The Navier-Stokes equations with consideration of the surface tension effects are solved by using SIMPLE algorithm on a non-staggered grid. The present code is validated by simulating two test problems. First one is to simulate a rising bubble inside a cube. The thickness of the interface of the bubble is shown to affect the pressure distribution around the interface. As the thickness decreases, the pressure field around the interface becomes more oscillatory. As the bubble rises, a ring vortex is shown to form around the interface and the bubble eventually develops into an ellipsoidal shape. Merge of two bubbles inside a container is secondly tested to show the robustness of the present code for two phase flow simulation. Numerical results show stable and reliable behavior during the process of merging of two bubbles. The velocity and pressure fields around the interface of bubbles are shown oscillation free during the merging of two bubbles.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.22 no.1
• /
• pp.1-7
• /
• 2013

Pneumatic cylinders are widely used to actuators in automatic equipments because they are relatively inexpensive, simple to install and maintain, offer robust design and operation, are available in a wide range of standard sizes and design alternatives. This paper presents a comparative study among the dynamic characteristics of meter-out and meter-in speed control of pneumatic cushion cylinders with a relief valve type cushion mechanism. Because of the nonlinear differential equations and a requirement for simultaneous iterative solution in a mathematical model of a double acting pneumatic cushion cylinder, a computer simulation is carried out to investigate pressure, temperature, mass flow rate in cushion chamber and displacement and velocity time histories of piston under various operating conditions. It is found that the piston velocity and pressure response in meter-in speed control are more oscillatory than with meter-out those when pneumatic cushion cylinders are driven at a high-speed. In meter-out speed control, the effective area of the flow control valve is larger than that of meter-in, and the supply pressure has to be much higher than the pressure required to move the load because it has also to overcome the back pressure in cushion chamber.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.267-277
• /
• 1993

A finite element method(FEM) is presented and applied to the two-dimensional bottom turbulent boundary layer. The time-dependent incompressible motion of a viscous fluid is formulated by using the well-known Navier-Stokes equations and vorticity equation in terms of the velocity and pressure fields. The general numerical formulation is based on Galerkin method and solved by introducing the mixing length theory of Prandtl for eddy kinematic viscosity of a turbulent flow field. Numerical computations of the transport of sediment on an arbitrary sea-bed due to wave motion in the turbulent boundary layer are carried out. The results obtained by the FEM made clear the difference in characteristic features between the boundary layer due to oscillatory flow and the boundary layer due to wave motion.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.7 no.3
• /
• pp.528-537
• /
• 1995

Natural convective flow and heat transfer characteristics in a partially opened enclosure fitted with a horizontal divider are investigated numerically. The enclosure is composed of a lower hot and a upper cold horizontal walls and adiabatic vertical walls. A divider is attached perpendicularly to the vertical insulated wall. The governing equations are solved by using the finite element method with Galerkin method. The computations have been carried out by varying the length of divider, the opening size, and the Rayleigh number based on the temperature difference between two horizontal walls and the enclosure height for air(Pr=0.71). As result, when the opening size is fixed, the intensity of the secondary flow is weaken as the length of divider increases. The maximum heat transfer rate over the upper cold wall occurs at a position bounded on the opening. However, when the length of divider is increased considerably, its maximum occurs at the right wall. The stability and frequency of oscillation are affected by the Rayleigh number and length of divider. The Nusselt number is increased with the increase of the opening size and the increase of Rayleigh number.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.3
• /
• pp.547-554
• /
• 1992

Response characteristics and false counts of laser and white light optical particle counters (OPC) have been studied as a function of particle size using monodisperse polystylen laterx (PSL) particles. Theoretical light scattering calculations for He-Ne laser based counter have been compared with the experimental results and thus good agreements have been found. The light scattering intensity in monochromatic light shows an oscillatory character for the transparent and spherical particles of PSL due to Mie resonance. Because of this effect, the response of the LAS-X OPC showed almost same responses in the diameter ranges of 0.4mu.m to 0.6mu.m and 0.7mu.m to 1.0mu.m for PSL particles. A laser optical particle counter with high flow rate applied for clean room has been studied to identify the noise sources. Three different manufacturer's clean room optical particle counters alos have been tested to measure the background noise level.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• v.y2005m4
• /
• pp.364-367
• /
• 2005

A series of computational simulations have been carried out for non-reacting and reacting flows in channel-type and divergent scramjet combustor configurations with and without a cavity. Transverse injection of hydrogen is considered with the injection pressure varying from 0.5 to 1.5 MPa. The detailed resolution of the oscillatory flow and flame dynamics are captured at sufficient scale to identify the underlying physical mechanisms, and the role of combustor configuration, cavity and amount of heat addition could be understood.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.11 no.8
• /
• pp.314-321
• /
• 2001

Static and oscillatory loss of stability of composite pipes conveying fluid is Investigated. The theory of than walled beams is applied and transverse shear. rotary inertia, primary and secondary warping effects are incorporated. The governing equations and the associated boundary conditions are derived through Hamilton's variational principle. The governing equations and the associated boundary conditions are transformed to an eigenvlaue problem which provides the Information about the dynamic characteristics of the system. Numerical analysis is performed by using extended Gelerkin method. Variation of critical velocity of fluid with fiber angles and mass patios of fluid to pipe Including fluid is investigated.

• Journal of Ocean Engineering and Technology
• /
• v.15 no.3
• /
• pp.9-19
• /
• 2001

In this paper we present the aspect of inertial oscillation typically observed in the spin-up of fluids at low Rossby numbers in a circular cylinder. Numerical computations for the quasi three-dimensional equation as well as one-dimensional equation are performed to estimate the predictability of the one-dimensional equation with Ekman pumping/suction models. It is assumed that the discrepancy between the two results may be attributed to the inertial oscillation The detailed analysis to the numerical results reveals that the axial plane is dominated by a comparatively strong oscillatory flows caused by the inertial oscillation. In view of the fact that the time-averaged flow field however agrees to the Taylor-Proudman theorem, it is recommended that further analysis is needed to obtain an improved one-dimensional model like the Reynolds-averaged Navier-Stokes equation for turbulent flows.

• Journal of the Korea Institute of Building Construction
• /
• v.11 no.1
• /
• pp.83-90
• /
• 2011

The purpose of this research is to investigate the effect of AFm formation on the stiffening process of cement paste. High and low alkali sulfate clinkers were used for the experiments. The flow and stiffening behavior of cement paste was investigated using modified ASTM C403 penetration resistance test and oscillatory shear rheology. X-ray powder diffraction (XRD) was used for phase identification associated with stiffening of the paste. It was found from the results that low alkali clinker mixture produced very strong premature stiffening whereas high alkali clinker mixture did not cause premature stiffening. This is because of the large amount of alkali sulfate present in the clinker. Addition of calcium and sodium chloride to the high alkali clinker mixture caused faster stiffening and set.