• Korean Journal of Mathematics
• /
• v.13 no.1
• /
• pp.113-122
• /
• 2005

We consider the iterative schemes for the large sparse linear system to solve partial differential equations. Using spectral radius of iteration matrices, the optimal relaxation parameters and good parameters can be obtained. With those parameters we compare the effectiveness of the SOR and SSOR algorithms. Applying Crank-Nicolson approximation, we observe the error distribution according to domain decomposition. The number of processors due to domain decomposition affects time and error. Numerical experiments show that effectiveness of SOR and SSOR can be reversed as time size varies, which is not the usual case. Finally, these phenomena suggest conjectures about equilibrium time grid for SOR and SSOR.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.13 no.2 s.54
• /
• pp.231-242
• /
• 2009

The baseline distribution of a structure represents the statistical distribution of dynamic response feature from the healthy state of the structure. Generally, damage-sensitive dynamic response feature of a structure manifest themselves near the tail of a baseline statistical distribution. In this regard, some researchers have paid attention to extreme value distribution for modeling the tail of a baseline distribution. However, few researches have been conducted to theoretically understand the extreme value distribution from a perspective of statistical damage assessment. This study investigates the asymptotic convergence of domain of attraction in extreme value distribution through parameter estimation, which is needed for reliable statistical damage assessment. In particular, the asymptotic convergence of a domain of attraction is quantified with respect to the sample size out of which each extreme value is extracted. The effect of the sample size on false positive alarms in statistical damage assessment is quantitatively investigated as well. The validity of the proposed method is demonstrated through numerically simulated acceleration data on a two span continuous truss bridge.

• Proceedings of the Korean Institute of Electrical and Electronic Material Engineers Conference
• /
• 2001.07a
• /
• pp.1036-1038
• /
• 2001

The sapphire orientation dependence of the crystallization of ZnO thin films has been studied using real-time synchrotron x-ray scattering. The amorphous ZnO thin films with a 2400-${\AA}$-thick were grown on sapphire(110) and sapphire(001) substrates by radio frequency magnetron sputtering at room temperature. The amorphous ZnO films were crystallization into epitaxial ZnO(002) grains both on the sapphire(110) and on the sapphire(001) substrates. The epitaxial quality, such as mosaic distribution and crystal domain size, of the ZnO grains on the sapphire(110) is high, similar to that of the ZnO grains on the sapphire(001). With increasing the annealing temperature to 600$^{\circ}C$, the mosaic distribution and the crystal domain size of ZnO(002) grains in the film normal direction was improved and decreased, respectively.

• Journal of the Korean Statistical Society
• /
• v.26 no.3
• /
• pp.289-297
• /
• 1997

It is well known that the spacings, the differences of two successive order statistics, in a random sample of size n from a distribution function F are independent and exponentially distributed if F is itself the exponential distribution. In this paper we obtain an asymptotically similar result on a fixed number of upper spacings as n .to. .infty. for a general F under the assumption that F is in the domain of attraction of some extreme value distribution. For a heavy or short tailed F, appropriate log transformations of the sample should be proceded to get the result. As a by-product, we also get that each upper spacing diverges in probability to .infty. and converges in probability to 0 as n .to. .infty. for a heavy and short tailed F, respectively, which is fully expected.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.700-710
• /
• 2008

A Bubble size distribution model has been developed for the numerical simulation of cryogenic high-speed cavitating flow of the turbo-pumps in the liquid fuel rocket engine. The new model is based on the previous one proposed by the authors, in which the bubble number density was solved as a function of bubble size at each grid point of the calculation domain by means of Eulerian framework with respect to the bubble size coordinate. In the previous model, the growth/decay of bubbles due to pressure difference between bubble and liquid was solved exactly based on Rayleigh-Plesset equation. However, the unsteady heat transfer between liquid and bubble, which controls the evaporation/condensation rate, was approximated by a theoretical solution of unsteady heat conduction under a constant temperature difference. In the present study, the unsteady temperature field in the liquid around a bubble is also solved exactly in order to establish an accurate and efficient numerical simulation code for cavitating flows. The growth/decay of a single bubble and growth of bubbles with nucleation were successfully simulated by the proposed model.

• Wind and Structures
• /
• v.37 no.4
• /
• pp.255-274
• /
• 2023

The aim of this paper is to enhance the accuracy of predicting time-averaged external surface pressures on low-rise buildings by utilizing Computational Fluid Dynamics (CFD) simulations. To achieve this, benchmark studies of the Silsoe cube and the Texas Tech University (TTU) experimental building are employed for comparison with simulation results. The paper is structured into three main sections. In the initial part, an appropriate domain size is selected based on the precision of mean pressure coefficients on the windward face of the cube, utilizing Reynolds Averaged Navier-Stokes (RANS) turbulence models. Subsequently, recommendations regarding the optimal computational domain size for an isolated building are provided based on revised findings. Moving on to the second part, the Silsoe cube model is examined within a horizontally homogeneous computational domain using more accurate turbulence models, such as Large Eddy Simulation (LES) and hybrid RANS-LES models. For computational efficiency, transient simulation settings are employed, building upon previous studies by the authors at the Windstorm Impact, Science, and Engineering (WISE) Lab, Louisiana State University (LSU). An optimal meshing strategy is determined for LES based on a grid convergence study. Three hybrid RANS-LES cases are investigated to achieve desired enhancements in the distribution of mean pressure coefficients on the Silsoe cube. In the final part, a 1:10 scale model of the TTU building is studied, incorporating the insights gained from the second part. The generated flow characteristics, including vertical profiles of mean velocity, turbulence intensity, and velocity spectra (small and large eddies), exhibit good agreement with full-scale (TTU) measurements. The results indicate promising roof pressures achieved through the careful consideration of meshing strategy, time step, domain size, inflow turbulence, near-wall treatment, and turbulence models. Moreover, this paper demonstrates an improvement in mean roof pressures compared to other state-of-the-art studies, thus highlighting the significance of CFD simulations in building aerodynamics.

• Journal of the Korean GEO-environmental Society
• /
• v.15 no.4
• /
• pp.13-27
• /
• 2014

Catchments soil erosion, one of the most serious problems in the mountainous environment of the world, consists of a complex phenomenon involving the detachment of individual soil particles from the soil mass and their transport, storage and overland flow of rainfall, and infiltration. Sediment size distribution during erosion processes appear to depend on many factors such as rainfall characteristics, vegetation cover, hydraulic flow, soil properties and slope. This study involved laboratory flume experiments carried out under simulated rainfall in a 3.0 m long ${\times}$ 0.8 m wide ${\times}$ 0.7 m deep flume, set at $17^{\circ}$ slope. Five experimental cases, consisting of twelve experiments using three different sediments with two different rainfall conditions, are reported. The experiments consisted of detailed observations of particle size distribution of the out-flow sediment. Sediment water mixture out-flow hydrograph and sediment mass out-flow rate over time, moisture profiles at different points within the soil domain, and seepage outflow were also reported. Moisture profiles, seepage outflow, and movement of overland flow were clearly found to be controlled by water retention function and hydraulic function of the soil. The difference of grain size distribution of original soil bed and the out-flow sediment was found to be insignificant in the cases of uniform sediment used experiments. However, in the cases of non-uniform sediment used experiments the outflow sediment was found to be coarser than the original soil domain. The results indicated that the sediment transport mechanism is the combination of particle segregation, suspension/saltation and rolling along the travel distance.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.1 s.232
• /
• pp.22-29
• /
• 2005

A material cloud method, which is a new topology optimization method, is presented. In MCM, an optimal structure can be found out by manipulating sizes and positions of material clouds, which are lumps of material with specified properties. A numerical analysis for a specific distribution of material clouds is carried out using fixed background finite element mesh. Optimal material distribution can be element-wisely extracted from material clouds' distribution. In MCM, an expansion-reduction procedure of design domain for finding out better optimal solution can be naturally realized. Also the convergence of material distribution is faster and well-defined material distribution with fewer intermediate densities can be obtained. In addition, the control of minimum-member sizes in the material distribution can be realized to some extent. In this paper, basic concept of MCM is introduced, and formulation and optimization results of MCM are compared with those of the traditional density distribution method(DDM).

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.4 s.175
• /
• pp.1016-1023
• /
• 2000

The present work is concerned with three-dimensional finite element analysis of the hollow section extrusion process using a porthole die. The effects of related design parameters are discussed through the finite element simulation for extrusion of a triply-connected rectangular tubular section. For economic computation, mismatching refinement, an efficient domain decomposition method with different mesh density for each subdomain, is implemented. In order to obtain the uniform flow at the outlet, design parameters such as the hole size and the hole position are investigated and compared through the numerical analysis. Comparing the velocity distribution with that of the original design, it is concluded that the design modification enables more uniform flow characteristics. The analysis results are then successfully reflected on the industrial porthole die design.

• Journal of the Korean Ceramic Society
• /
• v.54 no.6
• /
• pp.506-510
• /
• 2017

One of the characteristics of polycrystalline graphene that determines its material properties is grain size. Mechanical properties such as Young's modulus, yield strain and tensile strength depend on the grain size and show a reverse Hall-Petch effect at small grain size limit for some properties under certain conditions. While there is agreement on the grain size effect for Young's modulus and yield strain, certain MD simulations have led to disagreement for tensile strength. Song et al. showed a decreasing behavior for tensile strength, that is, a pseudo Hall-Petch effect for the small grain size domain up to 5 nm. On the other hand, Sha et al. showed an increasing behavior, a reverse Hall-Petch effect, for grain size domain up to 10 nm. Mortazavi et al. also showed results similar to those of Sha et al. We suspect that the main difference of these two inconsistent results is due to the different modeling. The modeling of polycrystalline graphene with regular size and (hexagonal) shape shows the pseudo Hall-Petch effect, while the modeling with random size and shape shows the reverse Hall-Petch effect. Therefore, this study is conducted to confirm that different modeling is the main reason for the different behavior of tensile strength of the polycrystalline structures. We conducted MD simulations with models derived from the Voronoi tessellation for two types of grain size distributions. One type is grains of relatively similar sizes; the other is grains of random sizes. We found that the pseudo Hall-Petch effect and the reverse Hall-Petch effect of tensile strength were consistently shown for the two different models. We suspect that this result comes from the different crack paths, which are related to the grain patterns in the models.