• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.9
• /
• pp.1058-1065
• /
• 2004

In the present study, characteristics of steady state laminar flows of a straight duct connected to a 180$^{\circ}$ curved duct were examined in the entrance region through experimental and numerical analyses. For the analysis, the governing equations of laminar flows in the Cartesian coordinate system were applied. Flow characteristics such as velocity profiles, and secondary flows were investigated numerically and experimentally in a square cross-sectional straight duct by the PIV system and a CFD code(STAR CD). For the PIV measurement, working fluid produced from mosquito coils smoke. The experimental data were obtained at 9 points dividing the test sections by 400 mm. Experimental and numerical results can be summarized as follows. Critical Reynolds number, Recr which indicates transition from laminar steady flow to transition steady flow was 2,150. As Reynolds number, Re, was increased, dimensionless velocity profiles at the outer wall were increased due to the effect of the centrifugal force and the secondary flows. The intensity of a secondary flow became stronger at the inner wall rather than the outer wall regardless of Reynolds number.

• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.247-254
• /
• 2000

The base pressure of vacuum vessel of the KSTAR (Korea Superconducting Tokamak Advanced Research) Tokamak is to be a ultra high vacuum, $10^{-6}{\sim}10^{-7}Pa$, to produce clean plasma with low impurity containments. For this purpose, the KSTAR vacuum vessel and plasma facing components need to be baked up to at least $250^{\circ}C,\;350^{\circ}C$ respectively, within 24 hours by hot nitrogen gas from a separate baking/cooling line system to remove impurities from the plasma-material interaction surfaces before plasma operation. Here by applying the implicit numerical method to the heat balance equations of the system, overall temperature distributions of the KSTAR vacuum vessel and plasma facing components are obtained during the whole baking process. The model for 2-dimensional baking analysis are segmented into 9 imaginary sectors corresponding to each plasma facing component and has up-down symmetry. Under the resulting combined loads including dead weight, baking gas pressure, vacuum pressure and thermal loads, thermal stresses in the vacuum vessel during bakeout are calculated by using the ANSYS code. It is found that the vacuum vessel and its supports are structurally rigid based on the thermal stress analyses.

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.46-53
• /
• 1997

The flow field around a three dimensional minivan-like body has been simulated. This study solves 3-D unsteady incompressible Navier-Stokes equations on a non-orthogonal curvilinear coordinate system using second-order accurate schemes for the time derivatives, and third/second-order scheme for the spatial derivatives. The Marker-and-Cell concept is applied to efficiently solve continuity equation. A H-H type of multi-block grid system is generated around a three dimensional minivan-like body. Turbulent flows have been modeled by the Baldwin-Lomax turbulent model. To validate present procedure, the flows around the Ahmed body with 12.5° of slant angle are simulated. A good agreement with other numerical results is achived. After code validation, the flows around a mimivan-like body are simulated. The simulation shows three dimensional vortex-pair just behind body. The flow separation is also observed on the rear of the body. It has concluded that the results of present study properly agreed with physical flow phenomena.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.155-158
• /
• 2002

In general, TVC(Thermal Vapor Compressor) is used to boost/compress a low pressure vapor to a higher pressure for further utilization. The one-dimensional method is simple and reasonably accurate, but cannot realize the detail as like the back flow and recirculation in the mixing chamber, viscous shear effect, and etc. In this study, the axisymmetric How simulations have been performed to reveal the detailed flow characteristics for the various ejector shapes. The Navier-Stokes and energy equations are solved together with the continuity equation In the compressible flow fields. The standard $k-{\epsilon}$ model is selected for the turbulence modeling. The commercial computational fluid dynamic code FLUENT software is used for the simulation. The results contain the entrainment ratio under the various motive, suction and discharge pressure conditions. The numerical results are compared with the experimental data, and the comparison shows the good agreement. The three different flow regimes (double chocking, single chocking and back flow) have been clearly distinguished according to each boundary pressure values. Also the effects of the various shape variables (nozzle position, nozzle outlet diameter, mixing tube diameter, mixing tube converging angle, and etc.) are quantitatively discussed.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.115-118
• /
• 2002

The present study deals with CFD analysis of a plastic heat exchanger with corrugated wall. This exchanger has sinusoidal corrugations, and the flow through the exchanger is three dimensional. In addition, CFX-5.4, a commercial code utilizing unstructured mesh, was used as a computational method for solving RANS(Reynolds-Averaged Navier-Stokes) equations, and the applied turbulence model is $k-{\varepsilon}$ model. The factors to affect the efficiency of a plastic heat exchanger are heat conductivity, flow characteristics and so on. For those two factors, heat conductivity is fixed by the wall material. Therefore, the How along the corrugation affects the efficiency more, provided the same material. In conclusion, the heat transfer enhancement of a plastic heat exchanger with corrugated wall can be recognized from the flow characteristics such as velocity streamline, local heat transfer coefficient, velocity contour, and pressure contour. To confirm the results, both of the measured and the computational data for pressure loss were compared with each other, and they were identical.

• Proceedings of the IEEK Conference
• /
• 2001.10a
• /
• pp.693-698
• /
• 2001

This paper establishes a coding method system including the liberation and size distribution of recycling materials in the shredder operation. Every particle in the shredded product becomes a code number using the liberation model and size distribution equation transforming of weight percentage into particles number percentage. One set of database can be obtained after all particles have been coded. This database is suitable for the size reduction operation in the process simulation of waste recycling. Coupling with the developed air classification, sizing and separating operations, the whole process simulation will be completely established for diversified application. A typical simulation for the rolling cutting shredder product of waste TV had been demonstrated under this coding system. The breakage size distribution of Gaudin and Schumann equation were selected for the shredding operation simulation. The Gaudin's liberation model was suitable fur the liberation simulation. Both of these equations were transformed weight percentage into particles distribution for the necessary of particle coding method. A better recycling operation for this shredded solid waste can be concluded from the comparison of simulation results with their sorted grade, recovery or economic of materials in different processes.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2006.06a
• /
• pp.265-270
• /
• 2006

High-frequency electromagnetic (EM) wave propagation phenomena associated with borehole ground-penetrating radar (GPR) surveys are complex. To improve the understanding of governing physical processes, we present a finite-difference time-domain solution of Maxwell's equations in cylindrical coordinates. This approach allows us to model the full EM wavefield associated with borehole GPR surveys. The algorithm can be easily implemented perfectly matched layers for absorbing boundaries, frequency-dependent media, and finite-length transmitter antenna.

• Journal of The Korean Astronomical Society
• /
• v.40 no.4
• /
• pp.165-169
• /
• 2007

The numerical simulations associated with the interaction of High Velocity Clouds (HVC) with the Magnetized Galactic Interstellar Medium (ISM) are a powerful tool to describe the evolution of the interaction of these objects in our Galaxy. In this work we present a new project referred to as Theoretical Virtual i Observatories. It is oriented toward to perform numerical simulations in real time through a Web page. This is a powerful astrophysical computational tool that consists of an intuitive graphical user interface (GUI) and a database produced by numerical calculations. In this Website the user can make use of the existing numerical simulations from the database or run a new simulation introducing initial conditions such as temperatures, densities, velocities, and magnetic field intensities for both the ISM and HVC. The prototype is programmed using Linux, Apache, MySQL, and PHP (LAMP), based on the open source philosophy. All simulations were performed with the MHD code ZEUS-3D, which solves the ideal MHD equations by finite differences on a fixed Eulerian mesh. Finally, we present typical results that can be obtained with this tool.

• Structural Engineering and Mechanics
• /
• v.65 no.4
• /
• pp.369-380
• /
• 2018

In this paper, a multi-layered finite element model for laminated glass plates is introduced. A layer-wise theory is applied to the analysis of laminated glass due to the combination of stiff and soft layers; the independent layers are connected via Lagrange multipliers. The von $K{\acute{a}}rm{\acute{a}}n$ large deflection plate theory and the constant Poisson ratio for constitutive equations are assumed to capture the possible effects of geometric nonlinearity and the time/temperature-dependent response of the plastic foil. The linear viscoelastic behavior of a polymer foil is included by the generalized Maxwell model. The proposed layer-wise model was implemented into the MATLAB code and verified against detailed three-dimensional models in ADINA solver using different hexahedral finite elements. The effects of temperature, load duration, and creep/relaxation are demonstrated by examples.

• Journal of computational fluids engineering
• /
• v.4 no.2
• /
• pp.17-30
• /
• 1999

The unsteady transitional boundary layer due to rotor-stator interaction was studied at two operation points, the design and one off design points. The off design point leads to lower blade loading and lower Reynolds number. A Navier-Stokes code developed in the previous study was parallelized to expedite computations. A low Reynolds number turbulence model was used to close the momentum equations. All computations show good agreement with experimental data. The wake induced transitional strip on the suction side of the stator is clearly captured at design point operation. There is no noticeable change in shape and phase angle of the wake induced strip even in the laminar sublayer. The wake induced transitional strip at off design point shows more complex structure. The wake induced transitional strip is observed only in the turbulent layer, and becomes obscure in the laminar sublayer and buffer layer. This behavior is probably consequent upon that the transition is governed by both wake induced strip and natural transition mechanism by Tollmien-Schlichting wave.