• Economic and Environmental Geology
• /
• v.54 no.4
• /
• pp.465-473
• /
• 2021

This numerical study aims at analyzing the effect of flow characteristics, caused by geometrical features such as intersecting angles, on solute mixing at fracture junctions. It showed that not only Pe, the ratio of advection to diffusion, but also the intersecting angles played an important role in solute mixing at the junction. For the intersection angles less than 90°, the fluid flowed to the outlet in the same direction as the injected flow direction, which increased the contact at the junction with the streamlines coming from the different inlets. On the other hand, for the intersecting angles greater than 90°, the fluid flowed out to the outlet opposite to the flow direction in the inlet, leading to minimizing the contact at the junction. Therefore, in the former case, solute mixing occurred even at high Pe, and in the latter case, solutes transport along the streamlines even at low Pe. For Pe < 1, the complete mixing model was known to occur, but for the intersecting angle greater than 150°, no complete solute mixing occurred. Overall, the transition from the complete mixing model to the streamline-routing model occurred for Pe = 0.1 - 100, but it highly depended on the intersecting angles. Specifically, the transition occurred at Pe = 0.1 - 10 for intersecting angles ≧ 150° and at Pe = 10 - 100 for intersecting angles ≦ 30°. For Pe > 100, the streamline-routing model was dominant regardless of intersecting angles. For Pe > 1,000, the complete streamline-routing model appeared only for the intersecting angles greater than 150°. For the intersecting angles less than 150°, the streamline-routing model dominated over the complete solute mixing, but solute mixing still occurred at the fracture junction.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.1
• /
• pp.12-24
• /
• 1998

The paper presents the mechanism of secondary flows and the associated total pressure losses occurring in turbine cascades with turning angle of about 127 and 77 degree. Velocity and pressure measurements are taken in seven traverse planes through the cascade passage using a prism type five hole probe. Oil-film flow visualization is also conducted on blade and endwall surfaces. The characteristics of the limiting streamlines show that the three dimensional separation is an important flow feature of endwall and blade surfaces. The larger turning results in much stronger contribution of the secondary flows to the loss developing mechanism. A large part of the endwall loss region at downstream pressure side is found to be very thin when compared to that of the cascade inlet and suction side endwall. Evolution of overall loss starts quite early within the cascade and the rate of the loss growth is much larger in the blade of large turning angle than in the blade of small turning angle.

• Journal of Advanced Marine Engineering and Technology
• /
• v.33 no.1
• /
• pp.62-70
• /
• 2009

In this paper, the fluid flow behavior past a pair of rectangular cylinders placed in a two dimensional horizontal channel has been investigated using Lattice-Boltzmann Method(LBM). The LBM has built up on the D2Q9 model and the single relaxation time method called the Lattice-BGK(Bhatnagar-Gross-Krook)model. Streamlines, velocity, vorticity and pressure contours are provided to analyze the important characteristics of the flow field for a wide range of non dimensional parameters that present in our simulation. Special attention is paid to the effect of spacing(d) between two cylinders and the blockage ratio A(=h/H), where H is the channel height and h is the rectangular cylinder height. for different Reynolds numbers. The first cylinder is called upstream cylinder and the second one as downstream cylinder. The downstream fluid flow fields have been more influenced by its blockage ratios(A) and Reynolds numbers(Re) whereas the upstream flow patterns(in front of downstream cylinder) by the gap length(d) between two cylinders. Moreover, it is observed that after a certain gap, both upstream and downstream flow patterns are almost similar size and shape. The simulation result has been compared with analytical solution and it is found to be in excellent agreement.

• Aerospace Engineering and Technology
• /
• v.7 no.2
• /
• pp.131-143
• /
• 2008

Slag mass deposition was required to predict performance accurately of KSLV-I kick motor(KM) system. The validation of the numerical analysis was performed with mass flow rate measured at 4th ground test of the KM. The study described here included internal flow field of KM at various time steps during burning. Slag mass accumulation was computed through the aluminum oxide particle paths to deviate from the gas flow streamlines in flight. These numerical analysis was performed with Fluent 6.3 program The effects for the acceleration, origins and diameters of the aluminum oxide particles was analyzed, finally the total slag mass accumulation was acquired. We confirmed that the slag mass deposition was agreement well with predicted slag mass based on kick motor the grounded test.

• Korea-Australia Rheology Journal
• /
• v.19 no.4
• /
• pp.211-219
• /
• 2007

Here we demonstrate complex transient behavior of viscoelastic liquid described numerically with the Leonov model in straight and contraction channel flow domains. Finite element and implicit Euler time integration methods are employed for spatial discretization and time marching. In order to stabilize the computational procedure, the tensor-logarithmic formulation of the constitutive equation with SUPG and DEVSS algorithms is implemented. For completeness of numerical formulation, the so called traction boundaries are assigned for flow inlet and outlet boundaries. At the inlet, finite traction force in the flow direction with stress free condition is allocated whereas the traction free boundary is assigned at the outlet. The numerical result has illustrated severe forward-backward fluctuations of overall flow rate in inertial straight channel flow ultimately followed by steady state of forward flow. When the flow reversal occurs, the flow patterns exhibit quite complicated time variation of streamlines. In the inertialess flow, it takes much more time to reach the steady state in the contraction flow than in the straight pipe flow. Even in the inertialess case during startup contraction flow, quite distinctly altering flow patterns with the lapse of time have been observed such as appearing and vanishing of lip vortices, coexistence of multiple vortices at the contraction comer and their merging into one.

• Water for future
• /
• v.21 no.3
• /
• pp.291-291
• /
• 1988

The two dimensional time dependent flow past a circular cylinder is analyzed numerically. In the analysis, equations of conservation of mass and momentum are transformed to equations of stream function-vorticity and vorticity transport, and nondimensionalized by nondimensional parameters representing flow characteristics, The resulting stream function-vorticity equstion and vorticity transport equation are solved by successive over relaxation scheme and alternating direct implicit scheme. Numerical experments are performed for the flow in the range of Reynolds number 125 to 275. The time dependent streamlines, vorticities, pressure on cylinder surface, separation angle, and drag and lift coefficients are calculated, and the method for estimation of pressure on cylinder surface and the outer boundary limit are developed.

• Journal of computational fluids engineering
• /
• v.21 no.4
• /
• pp.96-101
• /
• 2016

A two-dimensional Stokes flow past a circular disk in a circular tube is analyzed. The circular disk is located coaxially with the circular tube and the Hagen-Poiseuille flow exists at upstream and downstream far from the circular disk. The Stokes approximation is used and the flow is investigated analytically by using the method of eigenfunction expansion and the method of least square. From the analysis, the stream function and the pressure of the flow field are obtained, and the streamlines and pressure distribution are shown. Also, the pressure and shear stress distributions on the circular disk and circular tube wall are calculated, and shown for some typical radii of the circular disk. The additional pressure drop induced by the disk and the drag force exerted on the disk are compared as functions of the radius of the circular disk, and it is shown that the shear force on the wall of the tube increases due to the disk.

• Journal of Mechanical Science and Technology
• /
• v.16 no.10
• /
• pp.1346-1354
• /
• 2002

The purpose of this work is to study the effects of specularly reflecting wall under the combined radiative and laminar free convective heat transfer in an infinite square duct. An absorbing and emitting gray medium is enclosed by the opaque and diffusely emitting walls. The walls may reflect diffusely or specularly. Boussinesq approximation is used for the buoyancy term. The radiative heat transfer is evaluated using the direct discrete ordinates method. The parameters under considerations are Rayleigh number, conduction to radiation parameter, optical thickness, wall emissivity and reflection mode. The differences caused by the reflection mode on the stream line, and temperature distribution and wall heat fluxes are studied. Some differences are observed for the categories mentioned above if the order of the conduction to radiation parameter is less than order of 10$\^$-3/ fer the range of Rayleigh number studied. The differences at the side wall heat flux distributions are observed as long as the medium is optically thin. As the top wall emissivity decreases, the differences between these two modes are increased. As the optical thickness decreases at the fixed wall emissivity, the differences also increase. The difference of the streamlines or the temperature contours is not as distinct as the side wall heat flux distributions. The specular reflection may alter the fluid motion.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.12
• /
• pp.3135-3147
• /
• 1993

A numerical method is presented which can solve the steady flow and heat transfer from a rotating and heated circular cylinder in a uniform flow for a range of Reynolds number form 5 to 100. The steady response of the flow and heat transfer is simulated for various spin parameter. The effects on the flow field and heat transfer characteristics known as lift, drag and heat transfer coefficient are analyzed and the streamlines, velocity vectors, vorticity, temperature distributions around it were scrutinized numerically. As spin parameter increases the region of separation vortex becomes smaller than upper one and the lower region will vanish. The lift force, a large part is due to the pressure force, increases as the Reynolds number and it increases linearly as spin parameter increases. The pressure coefficient changes rapidly with spin parameter on the lower surface of the cylinder and the vorticity is sensitive to the spin parameter near separation region. As spin parameter increases the maximum heat coefficient and the thin thermal layer on front region are moved to direction of rotation. However, with balance between the local increase and decrease, the overal heat transfer coefficient is almost unaffected by rotation.

• Journal of the Society of Naval Architects of Korea
• /
• v.43 no.4 s.148
• /
• pp.460-466
• /
• 2006

The three-dimensional RANS equations were applied to analyze the flow field of a submarine. To validate the code, the DARPA SUBOFF bare hull and an eliipsoid at angles of attack of $10^{\circ}\;and\;30^{\circ}$ were simulated and good agreement with experiments was obtained. After the code validation, the flows over the full configuration of DARPA SUBOFF model having a fairwater and four stern appendages were simulated at four angles of attack $(0^{\circ},\;10^{\circ},\;20^{\circ},\;30^{\circ})$ and three yaw angles $(10^{\circ},\;20^{\circ},\;30^{\circ})$ Specifically, the pressure contours and streamlines of fairwater and stern appendage were compared as the angle of attack and yaw angle changed. The variations of hydrodynamic forces were also calculated.