• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.9
• /
• pp.34-40
• /
• 2005

In present study, a two dimensional unsteady Navier-Stokes solver has been developed using the Diagonalized ADI (DADI) method and implicit dual time stepping method. The jacobian matrices in steady state Navier-Stokes equations are introduced from inviscid flux terms. The implicit treatment of artificial dissipation terms results in a block penta-diagonal matrix system and it becomes a scalar penta-diagonal matrix by diagonalization. In steady state equations about fictitious time, a new residual including a real time derivative term is introduced. From a converged solution about fictitious time, a real time unsteady solution can be obtained, which is called 'implicit dual time stepping method'. For code validation, an oscillating flat plate, a regular Karman vortices past a circular cylinder and shock buffeting around a bicircular airfoil problems are numerically solved. And they are compared with a theoretical solution, experiments and other researcher's computations.

• Structural Engineering and Mechanics
• /
• v.22 no.5
• /
• pp.541-562
• /
• 2006

This paper compares the performance of axially loaded concrete filled steel tube (CFST) columns cast using a conventionally vibrated normal concrete (NC) and a novel self-consolidating concrete (SCC) made with a new viscosity modifying admixture (VMA). A total of sixteen columns with a standard compressive strength of about 50 MPa for both SCC and NC were tested by applying concentric axial load through the concrete core. Columns were fabricated without and with longitudinal and hoop reinforcement (Series I and Series II, respectively) in addition to the tube confinement. The slenderness of the columns expressed as height to diameter ratio (H/D) ranged between 4.8 and 9.5 for Series CI and between 3.1 and 6.5 for Series CII. The strength and ductility of SCC columns were found comparable to those of their NC counterparts as the maximum strength enhancement in NC columns ranged between 1.1% and 7.5% only. No significant difference in strain development was found due to the presence of SCC or NC or due to the presence of longitudinal and hoop reinforcement. Biaxial stress development in the steel tube as per von Mises yield criterion showed similar characteristics for both SCC and NC columns. The confined strength ($f^{\prime}_{cc}$) of SCC was found to be lower than that of NC and $f^{\prime}_{cc}$ also decreased with the increase of slenderness of the columns. Analytical models for the prediction of confined concrete strength and axial strength of CFST columns were developed and their performance was validated through test results. The proposed models were found to predict the axial strength of CFST columns better than existing models and Code based design procedures.

• 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 the Korea institute for structural maintenance and inspection
• /
• v.16 no.4
• /
• pp.107-115
• /
• 2012

The present study evaluates the validity of different equations specified in code provisions and proposed by the existing researchers to predict the concrete tensile capacities (direct tensile strength, splitting tensile strength and modulus of rupture) using a comprehensible database including 361 lightweight concrete (LWC), 1,335 normal-weight concrete (NWC) and 221 heavy-weight concrete (HWC) specimens. Most of the equations express the concrete tensile strengths as a function of its compressive strength based on the limited NWC concrete test data. However, the present database shows that the concrete tensile capacities are significantly affected by its unit weight as well. As a result, the inconsistency between experiments and predictions by the different models increases when the concrete unit weight is below 2,100 kg/$m^3$ and concrete compressive strength is above 50 MPa. On the other hand, new models proposed by the present study considering the concrete unit weight predict the tensile strengths of concrete with more accuracy.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.8 no.3
• /
• pp.437-450
• /
• 1996

A set of stability equations is formulated for natural convection flows adjacent to a vertical isothermal surface melting in cold pure water. It takes account of the nonparallelism of the base flows. The melting rate is regarded as a blowing velocity at the ice surface. The numerical solutions of the linear stability equations which constitute a two-point boundary value problem are accurately obtained for various values of the density extremum parameter $R=(T_m-T_{\infty})/(T_0-T_{\infty})$ in the range $0.3{\leq}R{\leq}0.6$, by using a computer code COLNEW. The blowing effects on the base flow becomes more significant as ambient temperature ($T_{\infty}$) increases to $T_{\infty}=10^{\circ}C$. The maximum decrease of heat transfer rate is about 6.4 percent. The stability results show that the melting at surface causes the critical Grashof number $G^*$ and the maximum frequency of disturbances to decrease. In comparision with the results for the conventional parallel flow model, the nonparallel flow model has a higher critical Grashof number but has lower amplification rates of disturbances than does the parallel flow model. The spatial amplification contours exhibit that the selective frequency $B_0$ of the nonparallel flow model is higher than that of the parallel flow model and that the effects of melting are rather small. The present study also indicates that the selective frequency $B_0$ can be easily predicted by the value of the frequency parameter $B^*$ at $G^*$, which comes from the neutral stability results of the nonparallel flow model.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.671-686
• /
• 2004

The shock-induced combustion ramjet (shcramjet) is a hypersonic airbreathing propulsion concept which over-comes the drawbacks of the long, massive combustors present in the scramjet by using a standing oblique detonation wave (a coupled shock-combustion front) as a means of nearly instantaneous heat addition. A novel shcramjet combustor design that makes use of wedge-shaped flameholders to avoid detonation wave-wall interactions is proposed and analyzed with computational fluid dynamics (CFD) simulations in this study. The laminar, two-dimensional Navier-Stokes equations coupled with a non-equilibrium hydrogen-air combustion model based on chemical kinetics are used to represent the physical system. The equations are solved with the WARP (window-allocatable resolver for propulsion) CFD code (see: Parent, B. and Sislian, J. P., “The Use of Domain Decomposition in Accelerating the Convergence of Quasihyperbolic Systems”, J. of Comp. Physics, Vol. 179, No. 1,2002, pages 140-169). The solver was validated with experimental results found in the literature. A series of steady-state numerical simulations was conducted using WARP and it was deter-mined by means of thrust potential calculations that this combustor design is a viable one for shcramjet propulsion: assuming a shcramjet flight Mach number of twelve at an altitude of 36,000 m, the geometrical dimensions used for the combustor give rise to an operational range for combustor inlet Mach numbers between six and eight. Different shcramjet flight Mach numbers would require different combustor dimensions and hence a variable geometry system in or-der to be viable.

• Journal of Ship and Ocean Technology
• /
• v.10 no.4
• /
• pp.12-23
• /
• 2006

In this paper, a numerical study is carried out to investigate the turbulent flow around a twin-skeg container ship model with rudders including propeller effects. A commercial CFD code, FLUENT is used with body forces distributed on the propeller disk to simulate the ship stem and wake flows with the propeller in operation. A multi-block, matching, structured grid system has been generated for the container ship hull with twin-skegs in consideration of rudders and body-force propeller disks. The RANS equations for incompressible fluid flows are solved numerically by using a finite volume method. For the turbulence closure, a Reynolds stress model is used in conjunction with a wall function. Computations are carried out for the bare hull as well as the hull with appendages of a twin-skeg container ship model. For the bare hull, the computational results are compared with experimental data and show generally a good agreement. For the hull with appendages, the changes of the stem flow by the rudders and the propellers have been analyzed based on the computed result since there is no experimental data available for comparison. It is found the flow incoming to the rudders has an angle of attack due to the influence of the skegs and thereby the hull surface pressure and the limiting streamlines are changed slightly by the rudders. The axial velocity of the propeller disk is found to be accelerated overall by about 35% due to the propeller operation with the rudders. The area and the magnitude of low pressure on the hull surface enlarge with the flow acceleration caused by the propeller. The propellers are found to have an effect on up to the position where the skeg begins. The propeller slipstream is disturbed strongly by the rudders and the flow is accelerated further and the transverse velocity vectors are weakened due to the flow rectifying effect of the rudder.

• Nuclear Engineering and Technology
• /
• v.22 no.1
• /
• pp.1-11
• /
• 1990

Heavy water is used as a moderator and a coolant in the pressurized heavy water reactor Because of the high cost of heavy water, downgraded heavy water generated in the reactor system is recycled to the reactor after being concentrated up to 99.8% or more in heavy water upgraders. This study investigates the process of upgraders and then suggests a theoretical model. The relations between process variables are derived from tower packing characteristics, vapour-liquid equilibria, and mass-heat balance equations at a steady state operation of the upgrader h computer program UPGR is developed, using the algorithm that solves the nonlinear equations step by step. It shows that the results of computer simulation are in good agreement with the operating data of the Wolsung upgrader. Thus, this computer code offers the optimum operating guide and is now applied to manage the performance of upgraders for the effective operation of the heavy water upgraders.

• Journal of the Institute of Convergence Signal Processing
• /
• v.24 no.2
• /
• pp.97-102
• /
• 2023

Numerical simulations have been performed to investigate the hydrodynamic characteristics of the rudder since they play an important role in naval architecture fields. Although some values such as hydrodynamics forces can be measured easily in the towing tanks, it is difficult to obtain the detailed information of the flow fields such as pressure distribution, velocity distribution, vortex generation from experiments. In the present study, the effects of hydrodynamic coefficients and Reynolds number acting on the rudder were studied by using Computational Fluid Dynamics(CFD). Ansys fluent, one of commercial CFD solvers, solves the Navier-Stokes equations and the k-epsilon turbulence model is selected for the viscous model to solve RANS equations. At first, drag coefficients and lift coefficient for different angle of attack are obtained by using a CFD commercial code for KCS rudder. Secondly, the 2-D lift coefficients and drag coefficients are compared with 3-D coefficients at the same conditions. Thirdly, the effects of Reynolds number on the hydrodynamic forces are investigated.

• Structural Engineering and Mechanics
• /
• v.89 no.2
• /
• pp.181-197
• /
• 2024

A numerical method is presented in this paper, for buckling analysis of thin arbitrary stiffened composite cylindrical shells under axial compression. The stiffeners can be placed inside and outside of the shell. The shell and stiffeners are operated as discrete elements, and their interactions are taking place through the compatibility conditions along their intersecting lines. The governing equations of motion are obtained based on Koiter's theory and solved by utilizing the principle of the minimum potential energy. Then, the buckling load coefficient and the critical buckling load are computed by solving characteristic equations. In this formulation, the elastic and geometric stiffness matrices of a single curved strip of the shell and stiffeners can be located anywhere within the shell element and in any direction are provided. Moreover, five stiffened composite shell specimens are made and tested under axial compression loading. The reliability of the presented method is validated by comparing its numerical results with those of commercial software, experiments, and other published numerical results. In addition, by using the ANSYS code, a 3-D finite element model that takes the exact geometric arrangement and the properties of the stiffeners and the shell into consideration is built. Finally, the effects of Poisson's ratio, shell length-to-radius ratio, shell thickness, cross-sectional area, angle, eccentricity, torsional stiffness, numbers and geometric configuration of stiffeners on the buckling of stiffened composite shells with various end conditions are computed. The results gained can be used as a meaningful benchmark for researchers to validate their analytical and numerical methods.