• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2401-2405
• /
• 2008

In the present study, a three-dimensional least square/level set based two-phase flow code was developed for the simulation of three-dimensional sloshing problems using finite element discretization. The present method can be utilized for the analysis of a free surface flow problem in a complex geometry due to the feature of FEM. Since the finite element method is employed for the spatial discretization of governing equations, an unstructured mesh can be naturally adopted for the level set simulation of a free surface flow without an additional load for the code development except that solution methods of the hyperbolic type redistancing and advection equations of the level set function should be devised in order to give a bounded solution on the unstructured mesh. From the numerical experiments of the present study, it is shown that the proposed method is both robust and accurate for the simulation of three-dimensional sloshing problems.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.9 s.240
• /
• pp.1049-1056
• /
• 2005

A conservative pressure-based finite-volume numerical method has been developed for computing flow and heat transfer by using an unstructured grid system. The method admits arbitrary convex polyhedra. Care is taken in the discretization and solution procedures to avoid formulations that are cell-shape-specific. A collocated variable arrangement formulation is developed, i.e. all dependent variables such as pressure and velocity are stored at cell centers. Gradients required for the evaluation of diffusion fluxes and for second-order-accurate convective operators are found by a novel second-order accurate spatial discretization. Momentum interpolation is used to prevent pressure checkerboarding and the SIMPLE algorithm is used for pressure-velocity coupling. The resulting set of coupled nonlinear algebraic equations is solved by employing a segregated approach, leading to a decoupled set of linear algebraic equations fer each dependent variable, with a sparse diagonally dominant coefficient matrix. These equations are solved by an iterative preconditioned conjugate gradient solver which retains the sparsity of the coefficient matrix, thus achieving a very efficient use of computer resources.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.8 no.5
• /
• pp.416-424
• /
• 1984

A numerical and experimental investigation was carried out on the single-phase natural circulation system. This study is concerned with the multiloop system which is relevant to the primary system of the pressurized water reactor. For numerical analysis, five time-dependent governing equations were derived using the one-dimensional lumped parameter model. These equations were discretized by the space-time integration technique, and a simplified computer program, SIMFARS, was developed to solve those discretized equations. Experiments were performed for two purposes-one is to validate the developed code, and the other is to understand the qualitative behavior of the natural circulation loop. Comparison of the computational results with experiments, and several experimental and numerical results are presented in this article.

• Nuclear Engineering and Technology
• /
• v.19 no.3
• /
• pp.198-204
• /
• 1987

A plant simulation code, MCSIM (Micro-Computer SIMulator), has been developed to simulate plant transient accidents for pressurized water reactors. Reactor coolant system is modeled using decoupled energy and momentum equations, drift flux two-phase flow model and integral momentum equation. A two-fluid pressurizer model is used to simulate the pressurizer dynamics. Pot Boiler model is used for steam generator, steady-state decoupled energy and momentum equations for secondary side system, and point kinetics equations for nuclear power calculation. For test of the present version of MCSIM, complete loss of flow and RCCA withdrawal accidents are calculated with MCSIM. The results are compared with those in FSAR of KNU 5 & 6.

• Computers and Concrete
• /
• v.17 no.6
• /
• pp.739-760
• /
• 2016

This paper reports on punching shear behavior of reinforced concrete panels, investigated experimentally and through finite element simulation. The aim of the study was to examine the punching shear of high strength concrete panels incorporating different types of aggregate and silica fume, in order to assess the validity of the existing code models with respect to the role of compressive and tensile strength of high strength concrete. The variables in concrete mix design include three types of coarse aggregates and three water-cementitious ratios, and ten-percent replacement of silica fume. The experimental results were compared with the results produced by empirical prediction equations of a number of widely used codes of practice. The prediction of the punching shear capacity of high strength concrete using the equations listed in this study, pointed to a potential unsafe design in some of them. This may be a reflection of the overestimation of the contribution of compressive strength and the negligence of the role of flexural reinforcement. The overall findings clearly indicated that the extrapolation of the relationships that were developed for normal strength concrete are not valid for high strength concrete within the scope of this study and that finite element simulation can provide a better alternative to empirical code Equations.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.55-58
• /
• 2002

Parallel finite element code has been recently developed for the analysis of the incompressible Wavier-Stokes equations using domain decomposition method. Metis and MPI libraries are used for the domain partitioning of an unstructured mesh and the data communication between sub-domains, respectively. For unsteady computation of the incompressible Navier-Stokes equations, 4-step splitting method is combined with P1P1 finite element formulation. Smagorinsky and dynamic model are implemented for the simulation of turbulent flows. For the validation performance-estimation of the developed parallel code, three-dimensional Laplace equation has been solved. It has been found that the speed-up of 40 has been obtained from the present parallel code fir the bench mark problem. Lastly, the turbulent flows around the MIRA model and Tiburon model have been solved using 32 processors on IBM SMP cluster and unstructured mesh. The computed drag coefficient agrees better with the existing experiment as the mesh resolution of the region increases, where the variation of pressure is severe.

• Journal of computational fluids engineering
• /
• v.11 no.1 s.32
• /
• pp.57-66
• /
• 2006

An optimal shape design approach is presented for a subsonic S-shaped intake using aerodynamic sensitivity analysis. Two-equation turbulence model is employed to capture strong counter vortices in the S-shaped duct more precisely. Sensitivity analysis is performed for the three-dimensional Navier-Stokes equations coupled with two-equation turbulence models using a discrete adjoint method For code validation, the result of the flow solver is compared with experiment data and other computational results of bench marking test. To study the influence oj turbulence models and grid refinement on the duct flow analysis, the results from several turbulence models are compared with one another and the minimum number of grid points, which can yield an accurate solution is investigated The adjoint variable code is validated by comparing the complex step derivative results. To realize a sufficient and flexible design space, NURBS equations are introduced as a geometric representation and a new grid modification technique, Least Square NURBS Grid Approximation is applied With the verified flow solver, the sensitivity analysis code and the geometric modification technique, the optimization of S-shaped intake is carried out and the enhancement of overall intake performance is achieved The designed S-shaped duct is tested in several off-design conditions to confirm the robustness of the current design approach. As a result, the capability and the efficiency of the present design tools are successfully demonstrated in three-dimensional highly turbulent internal flow design and off-design conditions.

• Architectural research
• /
• v.19 no.4
• /
• pp.117-124
• /
• 2017

In this research, we carried out finite element analysis depends on the variations such as the strength of the main bar, concrete, shear-span ratio(a/d) and existence of shear reinforcing bar. Throughout the results of FEM analysis, we were able to figure out how each variation can effect on shear performance. As the strength of concrete increased, the maximum shear force enhancement effect of each specimen was evaluated. As a result, the shear strengthening effect was 51~97% for shear reinforced specimens, and 26~44% for non-shear reinforced specimens. As the yield strength of reinforcing bars increases, the shear reinforcement effect of the specimen the specimens without shear reinforcement were 3%~6% higher than those with shear reinforcement. Theoretical and analytical values were compared using the design equations obtained from the CEB code. Theoretical and analytical values were compared using the design equations obtained from the CEB code. As a result, the error rate was the highest at 3.64 in the S1.0-C0 series and the lowest at 1.46 in the S1.7-C1 series. Therefore, the design equation of the CEB code is estimated to underestimate the actual shear strength of deep beams that are not subjected to shear reinforcement.

• Journal of the Korean Recycled Construction Resources Institute
• /
• v.7 no.2
• /
• pp.145-150
• /
• 2019

This paper proposes a model for the calculation of the ultimate torsional strength in normal-strength and high-strength concrete beams which include the concrete contribution strength and use a reasonable thickness of shear flow. The adequacy of the proposed model is evaluated by comparing the calculated torsional strength with the experimentally observed results from 104 test specimens reported in the literature. The results are also compared with the calculations of the KCI and the ACI building code equations, and those of other model which include the concrete contribution strength. The comparisons show that the ultimate torsional strengths calculated by the proposed equation and Rahal's equation are closer to the experimentally observed results than those calculated by the code equations.

• Nuclear Engineering and Technology
• /
• v.49 no.4
• /
• pp.838-849
• /
• 2017

This paper proposes penalty factor equations that take into consideration the weld strength over-match given in the classified form similar to the revised equations presented in the Code Case N-779 via cyclic elastic-plastic finite element analysis. It was found that the $K_e$ analysis data reflecting elastic follow-up can be consolidated by normalizing the primary-plus-secondary stress intensity ranges excluding the nonlinear thermal stress intensity component, $S_n$ to over-match degree of yield strength, $M_F$. For the effect of over-match on $K_n{\times}K_{\nu}$, dispersion of the $K_n{\times}K_{\nu}$ analysis data can be sharply reduced by dividing total stress intensity range, excluding local thermal stresses, $S_{p-lt}$ by $M_F$. Finally, the proposed equations were applied to the weld between the safe end and the piping of a pressurizer surge nozzle in pressurized water reactors in order to calculate a cumulative usage factor. The cumulative usage factor was then compared with those derived by the previous $K_e$ factor equations. The result shows that application of the proposed equations can significantly reduce conservatism of fatigue assessment using the previous $K_e$ factor equations.