• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.10
• /
• pp.961-969
• /
• 2013

In this study, the numerical analysis of a turbulent flow around both a staggered and an inline tube bundle was conducted using ANSYS CFX V.13, a commercial CFD software. The flow was assumed to be steady, incompressible, and isothermal. According to the CFD Best Practice Guideline, the sensitivity study for grid size, accuracy of the discretization scheme for convection term, and turbulence model was conducted, and its result was compared with the experimental data to estimate the applicability of the CFD Best Practice Guideline. It was concluded that the CFD Best Practice Guideline did not always guarantee an improvement in the prediction performance of the commercial CFD software in the field of tube bundle flow.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.3
• /
• pp.211-219
• /
• 2008

A three-dimensional inviscid flow solver has been developed based on unstructured meshes for the investigation of the effect of the external stores on the tail surfaces of a generic fighter aircraft. The numerical method is based on a vertex-centered finite-volume discretization and an implicit point Gauss-Seidel time integration. The calculations were made for a steady flow and the computed results were compared with experimental data to validate the flow solver. An unsteady time-accurate computation of the generic fighter aircraft with external stores at transonic flight conditions showed that the external stores cause unsteady loading on the horizontal tail surface due to the mutual interference between their wake and the horizontal tail surface. It was shown that downward deflection of the trailing edge flap significantly reduces the undesirable interference effect.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.5 no.6
• /
• pp.21-27
• /
• 2001

In this paper, a nonlinear finite element procedure is presented for the dynamic analysis of reinforced concrete shells. A computer program, named RCAHEST(reinforced concrete analysis in higher evaluation system technology), for the analysis of reinforced concrete structures was used. A 4-node flat shell element will drilling rotational stiffness is used for spatial discretization. The layered approach is used to discretize behavior of concrete and reinforcement through the thickness. Material nonlinearity is taken into account by comprising tensile, compressive and shear models of cracked concrete and a model of reinforcing steel. The smeared crack approach is incorporated. Solution of the equations of motion is obtained by numerical integration using Hilber-Hughes-Taylor(HHT) algorithm. The proposed numerical method for the nonlinear dynamic analysis of reinforced concrete shells is verified by comparison with reliable analytical results.

• Journal of Energy Engineering
• /
• v.5 no.1
• /
• pp.50-55
• /
• 1996

A numerical simulation on a primary calcinator of porcelain was performed with using Fluent to calculate the heat efficiency by studying velocity vector and temperature profile according to variables such as the location of outlet and porcelain. Control-Volume based Finite Difference Method and Up-wind scheme are used for discretization of differential equation. SIMPLEC Algorithm and standard k-$\varepsilon$ turbulent model are selected to resolve the pressure-velocity coupling and the turbulent. The result of simulation showed that the whole velocity vector field in a calcinator was varied greatly according to the location of outlet. But the whole temperature profile at each zone was still high regardless of the location of outlet because of the radiation. But the temperature of a case with a outlet at sidepart of preheating or cooling zone was little high compared to the case with a outlet on the top of preheating zone. The velocity vector field and temperature profile in a calcinator were almost not affected by the location of porcelain, but the temperature inside a porcelain was much affected according to the place where it was located. The heat efficiency in a calcinator was 44.6% and the gas temperature in the outlet was about 1000 K.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.5
• /
• pp.1-9
• /
• 2003

Detailed flow of a shrouded tail rotor in hover is studied by using a compressible inviscid flow solver on unstructured meshes. The numerical method is based on a cell-centered finite-volume discretization and an implicit Gauss-Seidel time integration. Numerical simulation is made for a single blade attached to the center body and guide by the duct by imposing a periodic boundary condition between adjacent rotor blades. The results show that the performance of an isolated rotor without shroud compares well with experiment. In case of a shrouded rotor, correction of the collective pitch angle is made such that the overall performance matches with experiment to account for the uncertainties of the experimental model configuration. Details of the flow field compare well with the experiment confirming the validity of the present method.

• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.6
• /
• pp.461-468
• /
• 2012

In this study, a Cartesian-grid method based on finite volume approach is applied to simulate the ship motions in large amplitude waves. Fractional step method is applied for pressure-velocity coupling and TVD limiter is used to interpolate the cell face value for the discretization of convective term. Water, air, and solid phases are identified by using the concept of volume-fraction function for each phase. In order to capture the interface between air and water, the tangent of hyperbola for interface capturing (THINC) scheme is used with weighed line interface calculation (WLIC) method which considers multidimensional information. The volume fraction of solid body embedded in the Cartesian grid system is calculated using a level-set based algorithm, and the body boundary condition is imposed by a volume weighted formula. Numerical simulations for the two-dimensional barge type model and Wigley hull in linear waves have been carried out to validate the newly developed code. To demonstrate the applicability for highly nonlinear wave-body interactions such as green water on the deck, numerical analysis on the large-amplitude motion of S175 containership is conducted and all computational results are compared with experimental data.

• Coupled systems mechanics
• /
• v.9 no.2
• /
• pp.163-182
• /
• 2020

Coupled thermo-mechanical analysis of reinforced concrete slab at elevated temperatures from a fire accounting for nonlinear thermal parameters is carried out. The main focus of the paper is put on a one-way continuous reinforced concrete slab exposed to fire from the single (bottom) side as the most typical working condition under fire loading. Although contemporary techniques alongside the fire protection measures are in constant development, in most cases it is not possible to avoid the material deterioration particularly nearby the exposed surface from a fire. Thereby the structural fire resistance of reinforced concrete slabs is mostly influenced by a relative distance between reinforcement and the exposed surface. A parametric study with variable concrete cover ranging from 15 mm to 35 mm is performed. As the first part of a one-way coupled thermo-mechanical analysis, transient nonlinear heat transfer analysis is performed by applying the net heat flux on the exposed surface. The solution of proposed heat analysis is obtained at certain time steps of interest by α-method using the explicit Euler time-integration scheme. Spatial discretization is done by the finite element method using a 1D 2-noded truss element with the temperature nodal values as unknowns. The obtained results in terms of temperature field inside the element are compared with available numerical and experimental results. A high level of agreement can be observed, implying the proposed model capable of describing the temperature field during a fire. Accompanying thermal analysis, mechanical analysis is performed in two ways. Firstly, using the guidelines given in Eurocode 2 - Part 1-2 resulting in the fire resistance rating for the aforementioned concrete cover values. The second way is a fully numerical coupled analysis carried out in general-purpose finite element software DIANA FEA. Both approaches indicate structural fire behavior similar to those observed in large-scale fire tests.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.9
• /
• pp.1562-1570
• /
• 2003

The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. The radial return mapping is one of the most robust integration scheme currently used. Nonlinear kinematic hardening model of Armstrong-Fredrick type has recovery term and the direction of kinematic hardening increment is not parallel to that of plastic strain increment. In this case, The conventional radial return mapping method cannot be applied directly. In this investigation, we expanded the radial return mapping method to consider the nonlinear kinematic hardening model and implemented this integration scheme into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using Newton method and bisection method. Using dynamic yield condition derived from linearization of flow rule, the integration scheme for elastoplastic and viscoplastic constitutive model was unified. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.1
• /
• pp.94-101
• /
• 1994

This paper describes a short-term prediction model for the movement of an oil slick in shallow waters. Under the assumption that the initial movement of the oil slick is governed by spreading and advection, the model has been developed and applied to Kyungki-Bay near Incheon Harbor. The initial spreading is estimated by using an empirical formula. The depth-averaged momentum equations are solved numerically for the volume transport velocities, in which the $M_2$ tide is the main driving source. A staggered grid system is adopted fur spatial discretization and the half-time method is implemented for time marching. The numerical result is visualized with the help of animation and thus the contaminated area is displayed on a monitor in time sequence. The input data are the time, the location and the volume of spill accident as well as environmental data such as md and $M_2$ tide.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2014.02a
• /
• pp.221-221
• /
• 2014

The present work proposes an improved numerical simulator for design and modification of large area capacitively coupled plasma (CCP) processing chamber. CCP, as notoriously well-known, demands the tremendously huge computational cost for carrying out transient analyses in realistic multi-dimensional models, because electron dissociations take place in a much smaller time scale (${\Delta}t{\approx}10-8{\sim}10-10$) than time scale of those happened between neutrals (${\Delta}t{\approx}10-1{\sim}10-3$), due to the rf drive frequencies of external electric field. And also, for spatial discretization of electron flux (Je), exponential scheme such as Scharfetter-Gummel method needs to be used in order to alleviate the numerical stiffness and resolve exponential change of spatial distribution of electron temperature (Te) and electron number density (Ne) in the vicinity of electrodes. Due to such computational intractability, it is prohibited to simulate CCP deposition in a three-dimension within acceptable calculation runtimes (<24 h). Under the situation where process conditions require thickness non-uniformity below 5%, however, detailed flow features of reactive gases induced from three-dimensional geometric effects such as gas distribution through the perforated plates (showerhead) should be considered. Without considering plasma chemistry, we therefore simulated flow, temperature and species fields in three-dimensional geometry first, and then, based on that data, boundary conditions of two-dimensional plasma discharge model are set. In the particular case of SiH4-NH3-N2-He CCP discharge to produce deposition of SiNxHy thin film, a cylindrical showerhead electrode reactor was studied by numerical modeling of mass, momentum and energy transports for charged particles in an axi-symmetric geometry. By solving transport equations of electron and radicals simultaneously, we observed that the way how source gases are consumed in the non-isothermal flow field and such consequences on active species production were outlined as playing the leading parts in the processes. As an example of application of the model for the prediction of the deposited thickness uniformity in a 300 mm wafer plasma processing chamber, the results were compared with the experimentally measured deposition profiles along the radius of the wafer varying inter-electrode gap. The simulation results were in good agreement with experimental data.