• Journal of Korean Tunnelling and Underground Space Association
• /
• v.8 no.2
• /
• pp.183-196
• /
• 2006

Rock joint roughness obtained from the camera-type 3D scanner was classified into waviness and unevenness. The classification criteria were established in the previous study; digital filtering was used to distinguish one from another. The classified and original profiles were used to produce metal moulds. For accurate machining of the moulds, the WEDM(Wire-cut Electric Discharge Machining) was adopted. Specimens were cast using high strength gypsum, and joint shear tests were performed by varying normal stress from low value to high one. Roughness mobilization characteristics depending on the asperity scale and the applied normal stress were investigated. A new equation was proposed to predict shear strength of rock joint, which can consider the characteristics of roughness mobilization and roughness parameters. The roughness quantification composed of waviness and unevenness was found to be a useful method to predict the joint shear strength.

• Journal of Korea Water Resources Association
• /
• v.30 no.2
• /
• pp.143-154
• /
• 1997

A numerical model for solving advection-diffusion equation is presented by splitoperator method combining the Holly-Preissmann scheme with a fifth-degree interpolating polynomial for advection operator and the explicit scheme porposed by Hobson et al. for diffusion operator. To examine the developed model, the obtained numerical solutions are compared with both the analytic solution and those from the existing models for the instantaneous source (Gaussian hill) and the continuous source (advanced front) at upstream boundary with constant velocity and diffusivity condition. For the various cases having different Courant and Peclet numbers, it is shown that the present study provides stable solutions even for Courant numbers exceeding one. The result obtained by the present study also agree well with existing analytical solutions for both cases. The proposed explicit scheme somewhat releases the conventional restriction of explicit schemes for determining the time step size and provides satisfactory results for relatively large time step size.

• Journal of the Korean Society of Systems Engineering
• /
• v.16 no.1
• /
• pp.1-8
• /
• 2020

In this paper we proposed to develop the reliability growth planning for the One-shot system using the PM2-Discrete model. The PM2-Discrete is the methodology specifically developed for discrete systems and is the first quantitative method available for formulating detailed plans in the discrete usage domain. First, the parameters RG, RI, T, MS and d of the PM2-Discrete model are set. Second, the case analysis was performed on One-shot system A. Third, the input parameter values were applied to drive the R(t) equation. Finally, using RGA 11 Software, the reliability Growth Planning Curve of One-shot system A was constructed. Also, the sensitivity analyses are performed for the changes of model parameters. The results of this study can be usefully used in establishing the reliability growth planning curve of the One-shot system.

• Journal of IKEEE
• /
• v.6 no.2 s.11
• /
• pp.136-145
• /
• 2002

The simulator for the analysis of the lattice temperature under the steady-state condition is developed. The heat flow equation using the Slotboom variables is discretized and the integration method of the thermal conductivity without using the numerical analysis method is presented. The simulations are executed on the $N^+P$ junction diode and BJT to verify the proposed method. The average relative error of the lattice temperature of $N^+P$ diode compared with DAVINCI is 2% when 1.4[V] forward bias is applied and the average relative error of the lattice temperature of BJT compared with MEDICI is 3% when 5.0[V] is applied to the collector contact and 0.5[V] is applied to the base contact. BANDIS using the proposed method of integration of thermal conductivity needs 3.45 times of matrix solution to solve one bias step and DAVINCI needs 5.1 times of matrix solution MEDICI needs 4.3 times of matrix solution.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.3
• /
• pp.241-248
• /
• 1987

The theoretical backgrounds of the several methods were surveyed and reviewed to fin out the adequate one to determine equivalent modal damping values in solving the dynamic problem of soil-structure interaction by mode superposition method. Furthermore the rigorous damping matrix of equation of motion was obtained through component mode synthesis technique and used in direct integration of the equation. The analytical results by direct integration method were compared with those of mode superposition approach using the various sets of equivalent modal damping values calculated by the methods to be reviewed. Two types of superstructures and four kinds of subsurface conditions were considered and combined to make soil-structure coupled models. It was realized that dissipating energy method gives the equivalent modal damping values which lead the most similar results to direct integration ones. In case of fixed base, the responses of all methods except stiffness weighted approach are almost equal to those of direct integration method.

• Nuclear Engineering and Technology
• /
• v.20 no.1
• /
• pp.9-18
• /
• 1988

Previously determining the fuel loading pattern is based on the trial and error method. For a candidate pattern, the core analysis is performed and the pattern is examined whether it satisfies the imposed constraints such as the power peaking or not. The pattern, then, is revised by the shuffling of assemblies and the revision is repeated until all of the conditions are met. This method unavoidably requires many iterative diffusion calculations, computing times and accumulated experiences. To overcome these disadvantages, a new method which is called backward diffusion calculation is introduced. If the most desirable power distribution is already known, the optimal loading pattern can be obtained by solving the backward diffusion equation with simple calculation. In this study, the basic equation for the backward diffusion calculation is derived and the optimal power and fuel distributions are searched in one-dimensional cylindrical geometry by using the proposed method. In addition, the basis to determine the optimal power and fuel distributions is suggested for the real core geometry.

• Nuclear Engineering and Technology
• /
• v.43 no.3
• /
• pp.271-278
• /
• 2011

When performing fire PSA in a nuclear power plant, an event mapping method, using an internal event PSA model, is widely used to reduce the resources used by fire PSA model development. Feasible initiating events and component failure events due to fire are identified to transform the fault tree (FT) for an internal event PSA into one for a fire PSA using the event mapping method. A surrogate event or damage term method is used to condition the FT of the internal PSA. The surrogate event or the damage term plays the role of flagging whether the system/component in a fire compartment is damaged or not, depending on the fire being initiated from a specified compartment. These methods usually require explicit states of all compartments to be modeled in a fire area. Fire event scenarios, when using explicit identification, such as surrogate or damage terms, have two problems: (1) there is no consideration of multiple fire propagation beyond a single propagation to an adjacent compartment, and (2) there is no consideration of simultaneous fire propagations in which an initiating fire event is propagated to multiple paths simultaneously. The present paper suggests a fire propagation equation to identify all possible fire event scenarios for an explicitly treated fire event scenario in the fire PSA. Also, a method for separating fire events was developed to make all fire events a set of mutually exclusive events, which can facilitate arithmetic summation in fire risk quantification. A simple example is given to confirm the applicability of the present method for a $2{\times}3$ rectangular fire area. Also, a feasible asymptotic approach is discussed to reduce the computational burden for fire risk quantification.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.55 no.6
• /
• pp.135-144
• /
• 2013

Ventilation efficiency has an important role in agricultural facilities such as greenhouse and livestock house to keep internally optimum environmental condition. Age-of-air concept allows to assess the ventilation efficiency of an agricultural facility according to estimating the ability of fresh air supply and contaminants emission using LMA and LMR. Most of these methods use a tracer gas method which has some limitations in experiment like dealing unstable and invisible gas. Therefore, the aim of this study was to develop a straightforward method to calculate age-of-air values with CFD simulation which has the advantage of saving computational time and resources and these method can solve the limitations in experiment using tracer gas method. The main idea of LMA computation is to solve the passive scalar transport equation with the assumption that the production of the time scalar throughout the room is uniform. In case of LMR calculation, the transport of the time scalar was reversed compulsively using UDF. The methodology to validate the results of this study was established by comparing with preceding research that had performed a computing LMA and LMR value by laboratory experiments and CFD simulations using tracer gas. As a result, the error was presented similarly level of results of preceding research. Some big errors could be caused by stagnated area and incongruity turbulence model. while the computational time was reduced to almost one fourth of that by preceding research.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.29 no.3
• /
• pp.138-144
• /
• 1987

This study was conducted to pursue the normalization of frequency distribution by making approach the coefficient of skewness to nearly zero tbrough SMEMAX transformation and to get probable minimum flows can be acquired by means of transforrnation equation which has been derivated by SMEMAX method to the annual minimum flow series of five watersheds along Geum river basin. The results obtained through SMEMAX method were compared with probable minimum flows according to return periods by Type III extremal distribution which has been determined as the best fitted one among probablility distributions for the analysis of minimum flow. All the results obtained through this study are summarized as follows. 1.SMEMAX transformation based on median value was proved to be the best method when the coefficient of skewness has less reliability because of the short duration for the observation and were not affected by accidental outliers. 2.SMEMAX transformation has found to be the best one for the coefficient of skewness to be made nearly zero in comparison with log and cubic root transformation. 3.Probable minimum flows according to the return periods were derivated by transformation equations obtained through theoretical analysis of SMEM AX transformation. 4.In general, probable minimum flows by SMEMAX method were appeared as higher values in the range of five and twenty years and as lower ones in the range of below than five and more than fifty years in return periods respectively, in comparison with the results of type III extremal distribution. 5.Relative errors in the probable minimum flows of SMEMAX method to the results of type III extremal distribution were shown to be within ten percent except those of one hundred years in return periods. 6.SMEMAX method was also confirmed to be useful for the analysis of minimum flow frequency as well as flood frequency analysis.

• The Journal of the Acoustical Society of Korea
• /
• v.31 no.1
• /
• pp.19-28
• /
• 2012

The time-domain solution from the Kirchhoff integral equation for an exterior problem is not unique at certain eigen-frequencies associated with the fictitious internal modes as happening in frequency-domain analysis. One of the solution methods is the CHIEF (Combined Helmholtz Integral Equation Formulation) approach, which is based on employing additional zero-pressure constraints at some interior points inside the body. Although this method has been widely used in frequency-domain boundary element method due to its simplicity, it was not used in time-domain analysis. In this work, the CHIEF approach is formulated appropriately for time-domain acoustic boundary element method by constraining the unknown surface pressure distribution at the current time, which was obtained by setting the pressure at the interior point to be zero considering the shortest retarded time between boundary nodes and interior point. Sound radiation of a pulsating sphere was used as a test example. By applying the CHIEF method, the low-order fictitious modes could be damped down satisfactorily, thus solving the non-uniqueness problem. However, it was observed that the instability due to high-order fictitious modes, which were beyond the effective frequency, was increased.