• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.9
• /
• pp.1183-1192
• /
• 2001

The numerical simulation has been conducted for the investigation of flow and heat transfer characteristics of an oblique impingement jet injected to a flat plate. The finite volume method was used to discretize the governing equations based on the non-orthogonal coordinate with non-staggered variable arrangement. The $textsc{k}$-$\varepsilon$-ν(sup)'2 turbulence model was employed to consider the consider the anisotropic flow characteristics generated by the impingement jet flow. The predicted results were compared with the experimental data and those of the standard $textsc{k}$-$\varepsilon$ turbulence model. The results of the $textsc{k}$-$\varepsilon$-ν(sup)'2 model showed better agreement with the experimental data than those of the standard $textsc{k}$-$\varepsilon$ model. In order to get the optimum condition, the flow and temperature fields were calculated with a variation of inclined angle($\alpha$=30$^{\circ}$~90$^{\circ}$) and the distance between the jet exit and impingement plate-to-diameter (L/D=4~10) at a fixed Reynolds number(Re=20,000). For a small L/D, the near-peak Nusselt numbers were not significantly effected by the inclined angle. The near-peak Nusselt numbers were not significantly affected by the L/D in the case of a large $\alpha$. The overall shape of the local Nusselt numbers was influenced by both the jet orifice-to-plate spacing and the jet angle.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2010.09a
• /
• pp.800-807
• /
• 2010

In this study, a series of numerical analyses has been performed in order to evaluate the performance of a full-scale closed-loop vertical ground heat exchanger constructed in Wonju. The circulation pipe HDPE, borehole and surrounding ground were modeled using FLUENT, a finite-volume method (FVM) program, for analyzing the heat transfer process of the system. Two user-defined functions (UDFs) accounting for the difference in the temperatures of the circulating inflow and outflow water and the change of the surrounding ground temperature with depth were adopted in the FLUENT model. The thermal properties of materials estimated in laboratory were used in the numerical analyses to compare the thermal efficiency of the cement grout with that of the bentonite grout used in the construction. The results of the simulation provide a verification of the in situ thermal response test data. The numerical model with the ground thermal conductivity of 4W/mK yielded the simulation result closer to the in-situ thermal response test than with the ground thermal conductivity of 3W/mK. From the results of the numerical analyses, the effective thermal conductivities of the cement and bentonite grouts were obtained to be 3.32W/mK and 2.99 W/mK, respectively.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.05a
• /
• pp.128-133
• /
• 2002

The heat transfer characteristics of a turbulent flow in a ribbed two-dimensional channel have been investigated numerically. The fully elliptic governing equations, coupled with a four-equation turbulence model, $\kappa-\omega-\bar{t^2}-\epsilon_t$, are solved by a finite volume method of SIMPLE type. Calculations have been carried out for three rib cross-sections : square, triangular, and semicircular, with various rib pitches and Reynolds numbers. The procedure appears to be satisfactory as the results for the square rib compare favorably with available experimental data and earlier calculation. The optimal rib pitch that yields the maximum heat transfer has been identified. It is also found that the square rib is most effective in enhancing the heat transfer. The semicircular rib, on the other hand, incurs the least amount of pressure drop but the improvement in heat transfer is substantially lower.

• Imaging Science in Dentistry
• /
• v.38 no.2
• /
• pp.81-87
• /
• 2008

Purpose: To investigate the relationship between 3D bone architectural parameters and direction-related elastic moduli of cancellous bone of mandibular condyle. Materials and Methods: Two micro-pigs (Micro-$pig^R$, PWG Genetics Korea) were used. Each pig was about 12 months old and weighing around 44 kg. 31 cylindrical bone specimen were obtained from cancellous bone of condyles for 3D analysis and measured by micro-computed tomography. Six parameters were trabecular thickness (Tb. Th), bone specific surface (BS/BV), percent bone volume (BV/TV), structure model index (SMI), degree of anisotropy (DA) and 3-dimensional fractal dimension (3DFD). Elastic moduli of three orthogonal directions (superior-inferior (SI), medial-lateral (ML), andterior-posterior (AP) direction) were calculated through finite element analysis. Results: Elastic modulus of superior-inferior direction was higher than those of other directions. Elastic moduli of 3 orthogonal directions showed different correlation with 3D architectural parameters. Elastic moduli of SI and ML directions showed significant strong to moderate correlation with BV/TV, SMI and 3DFD. Conclusion: Elastic modulus of cancellous bone of pig mandibular condyle was highest in the SI direction and it was supposed that the change into plate-like structure of trabeculae was mainly affected by increase of trabeculae of SI and ML directions.

• Advances in Energy Research
• /
• v.1 no.1
• /
• pp.53-78
• /
• 2013

Methane carbon dioxide reforming (MCDR) is a promising way of utilizing greenhouse gas for hydrogen-rich fuel production. Compared with other types of reactors, Compact Reformers (CRs) are efficient for fuel processing. In a CR, a thin solid plate is placed between two porous catalyst layers to enable efficient heat transfer between the two catalyst layers. In this study, the physical and chemical processes of MCDR in a CR are studied numerically with a 2D numerical model. The model considers the multi-component gas transport and heat transfer in the fuel channel and the porous catalyst layer, and the MCDR reaction kinetics in the catalyst layer. The finite volume method (FVM) is used for discretizing the governing equations. The SIMPLEC algorithm is used to couple the pressure and the velocity. Parametrical simulations are conducted to analyze in detail the effects of various operating/structural parameters on the fuel processing behavior.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.12
• /
• pp.2917-2925
• /
• 2000

Tensile strength and failure process of composite materials depend on the variation in fiber strength, matrix properties and fiber-matrix interfacial shear strength. A Monte-Carlo simulation considering variation in these factors has been widely used to analyze such a complicated phenomenon as a strength and simulated the failure process of unidirectional composites. In this study, a Monte Carlo simulation using 2-D and 3-D(square and hexagonal array) model was performed on unidirectional graphite/epoxy and glass/polyester composites. The results simulated by using 3-D hexagonal array model have a good agreement with the experimental data which were tensile strength and failure process of unidirectional composites.

• Steel and Composite Structures
• /
• v.42 no.2
• /
• pp.243-256
• /
• 2022

A coupled finite element method (FEM)-boundary element method (BEM) for analyzing the hydroelastic response of functionally graded carbon nanotube-reinforced composite (FG-CNTRC) floating plates under moving loads is firstly introduced in this article. For that aim, the plate displacement field is described utilizing a generalized shear deformation theory (GSDT)-based FEM, meanwhile the linear water-wave theory (LWWT)-relied BEM is employed for the fluid hydrodynamic modeling. Both computational domains of the plate and fluid are coincidentally discretized into 4-node Hermite elements. Accordingly, the C1-continuous plate element model can be simply captured owing to the inherent feature of third-order Hermite polynomials. In addition, this model is also completely free from shear correction factors, although the shear deformation effects are still taken into account. While the fluid BEM can easily handle the free surface with a lower computational effort due to its boundary integral performance. Material properties through the plate thickness follow four specific CNT distributions. Outcomes gained by the present FEM-BEM are compared with those of previously released papers including analytical solutions and experimental data to validate its reliability. In addition, the influences of CNT volume fraction, different CNT configurations, water depth, and load speed on the hydroelastic behavior of FG-CNTRC plates are also examined.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.1
• /
• pp.60-66
• /
• 2000

In order to improve cavitation characteristics for a high-speed propeller, leading edge shape of a 2-D hydrofoil is investigated numerically and experimentally. For flowfield analysis around the leading edge, the incompressible Reynolds Averaged Navier-Stokes(RANS) equation is solved using the standard $k-\varepsilon$ turbulence model and a finite volume method(FVM). The cavitation thickness, which is occurred on hydrofoil surface, is predicted using the panel code. It is shown that the calculation codes predict the experimental trend fairly well. From these results, new hydrofoils are designed

• Wind and Structures
• /
• v.17 no.6
• /
• pp.647-670
• /
• 2013

The paper presents a numerical approach to study of fluid flow characteristics and to predict performance of wind turbines. The numerical model is based on Finite-volume method (FVM) discretization of unsteady Reynolds-averaged Navier-Stokes (URANS) equations. The movement of turbine blades is modeled using moving mesh technique. The turbulence is modeled using commonly used turbulence models: Renormalization Group (RNG) k-${\varepsilon}$ turbulence model and the standard k-${\varepsilon}$ and k-${\omega}$ turbulence models. The model is validated with the experimental data over a large range of tip-speed to wind ratio (TSR) and blade pitch angles. In order to demonstrate the use of numerical method as a tool for designing wind turbines, two dimensional (2-D) and three-dimensional (3-D) simulations are carried out to study the flow through a small scale Darrieus type H-rotor Vertical Axis Wind Turbine (VAWT). The flows predictions are used to determine the performance of the turbine. The turbine consists of 3-symmetrical NACA0022 blades. A number of simulations are performed for a range of approaching angles and wind speeds. This numerical study highlights the concerns with the self-starting capabilities of the present VAWT turbine. However results also indicate that self-starting capabilities of the turbine can be increased when the mounted angle of attack of the blades is increased. The 2-D simulations using the presented model can successfully be used at preliminary stage of turbine design to compare performance of the turbine for different design and operating parameters, whereas 3-D studies are preferred for the final design.

• Journal of Korea Water Resources Association
• /
• v.30 no.6
• /
• pp.597-610
• /
• 1997

One-dimensional numerical models using finite difference methods for unsteady sediment transport on alluvial river channel are developed. The Preissmann implicit scheme and the Lax-Wendroff two-step explicit scheme with the Method of Characteristics for water motion and a forward time centered space explicit scheme for sediment motion are developed to simulate the sediment transport rate and the variation of channel bed level. The program correctness of each model is successfully verified using volume conservation tests. The sensitivity studies show that higher peak stage level, steeper channel slope and longer flooding duration produce more channel bed erosion. and median grain size, $D_{50}=0.4mm$ give maximum volume loss in this study. Finally, the numerical models are found to produce reasonable results from the various sensitivity tests which reveal that the numerical models have properly responded to the changes of each model parameter.