• Nuclear Engineering and Technology
• /
• v.41 no.7
• /
• pp.875-884
• /
• 2009

We present three nuclear/hydrogen-related R&D activities being performed at KAIST: air-ingressed LOCA analysis code development, gas turbine analysis tool development, and hydrogen-production system analysis model development. The ICE numerical technique widely used for the safety analysis of water-reactors is successfully implemented into GAMMA, with which we solve the basic equations for continuity, momentum conservation, energy conservation of the gas mixture, and mass conservation of 6 species (He, N2, O2, CO, CO2, and H2O). GAMMA has been extensively validated using data from 14 test facilities. We developed a tool to predict the characteristics of HTGR helium turbines based on the throughflow calculation with a Newton-Raphson method that overcomes the weakness of the conventional method based on the successive iteration scheme. It is found that the current method reaches stable and quick convergence even under the off-normal condition with the same degree of accuracy. The dynamic equations for the distillation column of HI process are described with 4 material components involved in the HI process: H2O, HI, I2, H2. For the HI process we improved the Neumann model based on the NRTL (Non-Random Two-Liquid) model. The improved Neumann model predicted a total pressure with 8.6% maximum relative deviation from the data and 2.5% mean relative deviation, and liquid-liquid-separation with 9.52% maximum relative deviation from the data.

• Journal of Korean Society of Steel Construction
• /
• v.21 no.3
• /
• pp.221-232
• /
• 2009

The load distribution factor (LDF) values of simplified composite H beam panel bridges (SCHPBs) that were subjected to one lane and two lane loads were investigated using three dimensional finite element analyses with the computer program ABAQUS (2007). This study considered some design parameters such as the slab thickness, the steel plate thickness, the span length, and the continuity of the SCHPBs in the development of new LDFs. The distribution values that were obtained from these analyses were compared with those from the AASHTO Standard, LRFD, and the equations presented by Tarhini and Frederick, Huo et al., Back and Shin, and Cai. The AASHTO Standard distribution factors for SCHPBs were found to be very conservative. Sometimes, the distribution values from the finite element analyses for interior girders were similar to the results of the AASHTO LRFD, whereas the values for exterior girders were conservative in most cases. The new distribution values that were presented in this study produced LDFs that are more conservative than those from the finite element method. For the simple application of the design to SCHPBs, bridge engineers can use 0.42 for the interior girder and 0.32 for the exterior girder. The proposed values improve the current design procedure for the LDF problem and increase SCHPB design efficiency.

• Structural Engineering and Mechanics
• /
• v.7 no.2
• /
• pp.181-202
• /
• 1999

Dynamic response of axisymmetric arbitrary laminated composite cylindrical shell of finite length, using three-dimensional elasticity equations are studied. The shell is simply supported at both ends. The highly coupled partial differential equations are reduced to ordinary differential equations (ODE) with variable coefficients by means of trigonometric function expansion in axial direction. For cylindrical shell under dynamic load, the resulting differential equations are solved by Galerkin finite element method, In this solution, the continuity conditions between any two layer is satisfied. It is found that the difference between elasticity solution (ES) and higher order shear deformation theory (HSD) become higher for a symmetric laminations than their unsymmetric counterpart. That is due to the effect of bending-streching coupling. It is also found that due to the discontinuity of inplane stresses at the interface of the laminate, the slope of transverse normal and shear stresses aren't continuous across the interface. For free vibration analysis, through dividing each layer into thin laminas, the variable coefficients in ODE become constants and the resulting equations can be solved exactly. It is shown that the natural frequency of symmetric angle-ply are generally higher than their antisymmetric counterpart. Also the results are in good agreement with similar results found in literatures.

• Journal of the Society of Naval Architects of Korea
• /
• v.40 no.5
• /
• pp.36-42
• /
• 2003

A CFD simulation technique has been developed to handle the unsteady body motion with large amplitude by use of overlapping multi-block grid system. The three-dimensional, viscous and incompressible flow around body is investigated by solving the Navier-Stokes equations, and the motion of body is represented by moving effect of the grid system. Composite grid system is employed in order to deal with both the body motion with large amplitude and the condition of numerical wave maker in convenience at the same time. The governing equations, Navier-Stokes (N-S) and continuity equations, are discretized by a finite volume method, in the framework of an O-H type boundary-fitted grid system (inner grid system including test model) and a rectangular grid system (outer grid system including simulation equipments for generation of wave environments). If this study, several flow configurations, such as an oscillating cylinder with large KC number, are studied in order to predict and evaluate the hydrodynamic forces. Furthermore, the motion simulation of a Series 60 model advancing in a uniform flow under the condition of enforced roll motion of angle 20$^{\circ}$ is performed in the developed numerical wave tank.

• International Journal of Automotive Technology
• /
• v.7 no.7
• /
• pp.867-872
• /
• 2006

Fuel cells are devices that convert chemical energy directly into electrical energy. Owing to the high efficiency of the fuel cells, a large number of research work have been done during these years. Among many kinds of the fuel cells, a polymer electrolyte membrane fuel cell is such kind of thing which works under low temperature. Because of the specialty, it stimulated intense global R&D competition. Most of the major world automakers are racing to develop polymer electrolyte membrane fuel cell passenger vehicles. Unfortunately, there are still many problems to be solved in order to make them into the commercial use, such as the thermal and water management in working process of PEMFCs. To solve the difficulites facing the researcher, the analysis of the inner mechanism of PEMFC should be implemented as much as possible and mathematical modeling is an important tool for the research of the fuel cell especially with the combination of experiment. By regarding some of the assumptions and simplifications, using the finite element technique, a two-dimensional electrochemical mode is presented in this paper for the further comparison with experimental data. Based on the principals of the problem, the equations of electronic charge conservation equation, gas-phase continuity equation, and mass balance equation are used in calculating. Finally, modeling results indicate some of the phenomenon in a unit cell, and the relationships between potential and current density.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.189-194
• /
• 2001

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.

• Structural Engineering and Mechanics
• /
• v.38 no.4
• /
• pp.385-403
• /
• 2011

Continuous deep girders which transmit the gravity load from the upper wall to the lower columns have frequently long end shear spans between the boundary of the upper wall and the face of the lower column. This paper presents the results of tests and analyses performed on three 1:2.5 scale specimens with long end shear spans, (the ratios of shear-span/total depth: 1.8 < a/h < 2.5): one designed by the conventional approach using the beam theory and two by the strut-and-tie approach. The conclusions are as follows: (1) the yielding strength of the continuous RC deep girders is controlled by the tensile yielding of the bottom longitudinal reinforcements, being much larger than the nominal strength predicted by using the section analysis of the girder section only or using the strut-and-tie model based on elastic-analysis stress distribution. (2) The ultimate strengths are 22% to 26% larger than the yielding strength. This additional strength derives from the strain hardening of yielded reinforcements and the shear resistance due to continuity with the adjacent span. (3) The pattern of shear force flow and failure mode in shear zone varies depending on the amount of vertical shear reinforcement. And (4) it is necessary to take into account the existence of the upper wall in the analysis and design of the deep continuous transfer girders that support the upper wall with a long end shear span.

• Proceedings of the KOSOMBE Conference
• /
• v.1997 no.05
• /
• pp.414-417
• /
• 1997

This paper presents vessel contour for extracting features and segmentating narrow blood vessels down to a diameter of two pixels in digital subtraction angiographic image. We present a new tracking algorithm for contour, mainly blood vessels in DSA image, and extracting properties such as their intensities, diameters, and center lines by exploiting spatial continuity. The proposed algorithm comes to detect blood vessel's boundary using difference edge detector one of homogeneity operator and find a next centerline position by direction vector of edge information. This algorithm enhanced variation of vessel's diameter compared to Sun's tracking algorithm and lessoned to compute as direction vector decide adaptively entire vessel's direction field. The processed images are intended to support radiologists in diagnosis, radiation therapy planning, and surgical planning. The algorithm should be useful for automating angiographic analyses of blood vessels.

• 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.

• Computers and Concrete
• /
• v.4 no.2
• /
• pp.135-156
• /
• 2007

Three 3D nonlinear finite-element models are developed to study the behavior of concrete beams and plates with and without external reinforcement by fibre-reinforced plastic (FRP). All three models are formulated based upon the 3D theory of elasticity. The stress model is modified from the element developed by Ramtekkar, et al. (2002) to incorporate material nonlinearity in the formulation. Both transverse stress and displacement components are used as nodal degrees-of-freedom to ensure the continuity of both stress and displacement components between the elements. The displacement model uses only displacement components as nodal degrees-of-freedom. The transition model has both stress and displacement components as nodal degrees-of-freedom on one surface, and only displacement components as nodal degrees-of-freedom on the opposite surface. The transition model serves as a connector between the stress and the displacement models. The developed models are validated by comparing the results of the analyses with an existing experimental result. Parametric studies of the effects of the externally reinforced FRP on the load capacity of reinforced concrete (RC) beams and concrete plates are performed to demonstrate the practicality and the efficiency of the proposed models.