• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.5
• /
• pp.1158-1174
• /
• 1988

A numerical technique is developed for the solution of fully developed turbulent recirculating flow in the passage of variable area using the non-orthogonal coordinate transformation. In the numerical analysis, primitive pressure-velocity finite difference equations were solved by SIMPLER algorithm with 2-equation turbulence model and algebraic stress model (ASM). QUICK scheme on the differencing of convective terms which is free from the inaccuracies of numerical diffusion has been applied to the variable grids and the results compared with those from HYBRID scheme. In order to test the effect of streamline curvatures on turbulent diffusion Lee and Choi streamline curvature correction model which has been obtained by modifying the Leschziner and Rodi's model is testes. The ASM was also employed and the results are compared to those from another turbulence model. The results show that difference of convective differencing schemes and turbulence models give significant differences in the prediction of velocity fields in the expansion region and outlet region of the combustor, however show little differences in the parallel flow region.

• 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.23 no.8
• /
• pp.997-1009
• /
• 1999

Silica particle formation and growth process including chemical reaction, coagulation and sintering was studied in a counterflow diffusion flame burner. The counterflow geometry provides a one dimensional flow field, along the stagnation point streamline, which greatly simplifies interpretation of the particle growth characteristics. $SiCl_4$ has been used as the source of silicon in hydrogen/oxygen/argon flames. The temperature profiles obtained by calculation showed a good agreement with experiment data. Using one and two dimensional sectional method, aerosol dynamics equation in a flame was solved, and these two results were compared. The two dimensional section method can consider sintering effect and growth of primary particle during synthesis, thus it showed evolution of morphology of non-spherical particles (aggregates) using surface fractal dimension. The effects of flame temperature and chemical loading on particle dynamics were studied. Geometric mean diameter based on surface area and total number concentration followed the trend of experiment results, especially, the change of diameters showed the sintering effect in high temperature environment.

• Design & Manufacturing
• /
• v.14 no.2
• /
• pp.28-33
• /
• 2020

Recently, as the income level and quality of life have improved, the desire for a pleasant environment has increased, and a deodorization plan is required through thorough prevention and diffusion of odorous substances in waste treatment facilities recognized as hateful facilities, appropriate collection, and selection of the right prevention facilities. In this study, a waste disposal facility was modeled and computerized analysis for odor and ventilation analysis was conducted. Numerical analysis of the waste treatment facility was performed at the size of the actual plant. CATIA V5 R16 for numerical model generation and ANSYS FLUENT V.13 for general purpose flow analysis were used as analysis tools. The average air-age of the internal was 329 seconds, and the air-flow velocity was 0.384m/s. The odor diffusion analysis inside the underground pump room showed congestion-free air circulation through streamline distribution and air-age distribution. This satisfies the ASHRAE criteria. In addition, the results of diffusion analysis of odorous substances such as ammonia, hydrogen sulfide, methyl mercaptan and dimethyl sulfide were all expected to satisfy the regulatory standards. Particularly in the case of the waste loading area, the air-flow velocity was 0.297m/s, and the result of meeting the regulatory standards with 0.167ppm of ammonia, 0.00548ppm of hydrogen sulfide, 0.003ppm of methyl mercaptan, and 0.003ppm of dimethyl sulfide was found.

• Journal of applied mathematics & informatics
• /
• v.26 no.5_6
• /
• pp.1057-1069
• /
• 2008

We consider singularly perturbed Boundary Value Problems (BVPs) for third and fourth order Ordinary Differential Equations(ODEs) of convection-diffusion type with discontinuous source term and a small positive parameter multiplying the highest derivative. Because of the type of Boundary Conditions(BCs) imposed on these equations these problems can be transformed into weakly coupled systems. In this system, the first equation does not have the small parameter but the second contains it. In this paper a computational method named as 'An asymptotic finite element method' for solving these systems is presented. In this method we first find an zero order asymptotic approximation to the solution and then the system is decoupled by replacing the first component of the solution by this approximation in the second equation. Then the second equation is independently solved by a fitted mesh Finite Element Method (FEM). Numerical experiments support our theoritical results.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.6
• /
• pp.1640-1647
• /
• 1993

The numerical diffusion arising from streamline-to-grid skewness produces a deteriorating effect on the numerical accuracy. The QUICKER scheme to reduce the numerical diffusion requires more computational effort than the HYBRID scheme. This paper deals with the relative computational efficiencies of adopting QUICKER scheme with a coarser grid system and of adopting HYBRID scheme with a denser grid system. Laminar driven cavity flow with Re=400, 1000 is used as a test problem. It is found that QUICKER scheme with a coarser grid system is more efficient than the HYBRID scheme with a denser grid system.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.4
• /
• pp.865-875
• /
• 1988

A two-dimensional stepped wall jet was numerically investigated by applying three different models : One is the standard k-.epsilon. and the other is the modified k-.epsilon. model which takes account of the streamline curvature effect by modifying the Reynolds shear stress and a source term in the dissipation equation, and a third is curvature dependent third-order correlation model. In order to test the influences of the numerical result, both the upwind scheme and the skew-upwind scheme were sued for the computations. By comparing the numerical results with available experiments, it was found that the modified k-.epsilon. model gives best overall prediction accuracy only when the numerical diffusion is eliminated by using the skew-upwind scheme. The numerical scheme was found to have more pronounced effect on the accuracy of the turbulence computation than the turbulence models.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.12
• /
• pp.1663-1673
• /
• 2002

In the present paper, flow and heat transfer characteristics of confined impinging slot jets have been numerically investigated using a SIMPLE-based segregated SUPG finite element method. For laminar jets, it is shown that the skin friction coefficient obtained from the present SUPG formulation approaches the grid-independent Galerkin solution inducing negligible false diffusion in the flow field when a moderate number of grid points are used. For turbulent jets, the k-$\omega$turbulence model is adopted. The streamwise mean velocity and the heat transfer coefficient respectively agree very well with existing experimental data within limited ranges of parameters.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.82-88
• /
• 1998

As a par of a project related to the development of the design algorithm of a compact heat exchanger for the application of the electronic home appliances, the effect of the discreteness of the airflow boundary generated on the cooling fin surface on the heat transfer and pressure drop characteristics of the heat exchanger was studied numerically. In general, there are two critical design parameters seriously considered in the design of the heat exchanger; heat transfer rate(Q) and pressure drop coefficient(C/sub p/). Even though the higher heat transfer rate with lower pressure drop characteristics is required in a design of the heat exchanger, it is not an easy job to satisfy both conditions at the same time because these two parameters are phenomenally inversely proportional. To control the boundary layer thickness and its length along the streamline, the surface of the flat fin was modified to accelerate the heat transfer rate on the fin surface. To understand the effect of the discreted fin size(S/sub w/) and its location(S/sub h/) on the performance of the heat exchanger in the airflow field, the flat fin was modified as shown in Fig. 1. From this study, it was found that the smaller and more number of slits on the fin surface showed the higher energy diffusion rate. It means that the discreteness of the boundary layer is quite important on the heat transfer rate of the heat exchanger. On the other hand, if the fin surface configuration is very complex than needed, higher static pressure drop occurs than required in a system and it may be a reason of the induced aerodynamic noise in the heat exchanger.

• Journal of Mechanical Science and Technology
• /
• v.15 no.7
• /
• pp.895-905
• /
• 2001

This paper presents vector fields, three dimensional mean velocities, turbulent intensities, turbulent kinetic energy and Reynolds shear stresses measured in the X-Y plane of the gas swirl burner with a cone type baffle plate by using an X-type hot-wire probe. This experiment is carried out at the flow rates of 350 and 450ℓ/min which are equivalent to the combustion air flow rate necessary to release 15,000 kcal/hr in a gas furnace. The results show that the maximum axial mean velocity component exists around the narrow slits situated radially on the edge of a burner. Therefore, there is some entrainment of ambient air in the outer region of a burner. The maximum values of turbulent intensities occur around the narrow slits and in front of a burner up to X/R=1.5. Moreover, the turbulent intensity components show a relatively large value in the inner region due to the flow diffusion and mixing processes between the inclined baffle plate and the swirl vane. Consequently, the combustion reaction is expected to occur actively near these regions.