• Journal of Radiation Protection and Research
• /
• v.17 no.1
• /
• pp.43-56
• /
• 1992

As one of the methods to ameliorate the ray effects which are the nature of anomalous computational effects due to the discretization of the angular variable in discrete ordinates approximations, a computational program, named TWODET (TWO dimensional Discrete Element Transport), has developed in 2 dimensional cartesian coordinates system using the discrete elements method, in which the discrete angle quadratures are steered by the spatially dependent angular fluxes. The results of the TWODET calculation with K-2, L-3 discrete angular quadratures, in the problem of a centrally located, isotropically emitting flat source in an absorbing square, are shown to be more accurate than that of the DOT 4.3 calculation with S-10 full symmetry angular quadratures, in remedy of the ray effect at the edge flux distributions of the square. But the computing time of the TWODET is about 4 times more than that of the DOT 4.3. In the problem of vacuum boundaries just outside of the source region in an absorbing square, the results of the TWODET calculation are shown severely anomalous ray effects, due to the sudden discontinuity between the source and the vacuum, like as the results of the DOT 4.3 calculation. In the probelm of an external source in an absorbing square in which a highly absorbing medium is added, the results of the TWODET calculation with K-3, L-4 show a good ones like as, somewhat more than, that of the DOT 4.3 calculation with S-10.

• Journal of the Korean Institute of Illuminating and Electrical Installation Engineers
• /
• v.17 no.3
• /
• pp.110-118
• /
• 2003

Spacer dampers maintain the constant gaps between each conductor in a bundle conductor-transmission line, and are installed at proper intervals to keep a line from all sorts of damages derived from the vibration energy caused by mechanical or electrical external factors. It is most important to embody a technology which considers difficulties of maintenance and repair, and has optimum elements in order to prevent accidents such as destruction by fire or the snapping of a wire by the effect of vibration phenomenon coming from transmission lint In the resent thesis, therefore, the analysis of vibratory characteristics of spacer damper is set up by analytical methods such as the analysis of conductor motion's governing equation, the equation of spacer damper's motion, spacer damper-fastened wire's motion in a span and the numerical analysis of finite difference method. Furthermore, the installation distance between spacer dampers was scrutinized by simulations of various vibration phenomena which change at any time as actual conditions do, and hereafter difference method. Furthermore, the installation distance between spacer dampers was scrutinized by simulations of various vibration phenomena which change at any time as actual conditions do, and hereafter we will be able to analyze all kinds of vibration phenomena coming from a boltless spacer damper with 6 bundle conductor for 765kV transmission line based on new analytical methods.

• The Journal of Engineering Geology
• /
• v.18 no.4
• /
• pp.381-391
• /
• 2008

The eigenvalue technique is introduced to overcome the problem of truncation errors caused by temporal discretization of numerical analysis. The eigenvalue technique is different from simulation in that only the space is discretized. The spatially discretized equation is diagonized and the linear dynamic system is then decoupled. The time integration can be done independently and continuously for any nodal point at any time. The results of eigenvalue technique are compared with the exact solution and FEM numerical solution. The eigenvalue technique is more efficient than the FEM to the computation time and the computer storage in the same conditions. This technique is applied to the solute transport analysis in nonuniform flow fields around underground storage caverns. This method can be very useful for time consuming simulations. So, a sensitivity analysis is carried out by using this method to analyze the safety of caverns from nearly located contaminant sources. According to the simulations, the reaching time from source to the nearest cavern may takes 50 years with longitudinal dispersivity of 50 m and transversal dispersivity of 5 m, respectively.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.20 no.5
• /
• pp.612-620
• /
• 2019

The increasing use of curved beams in buildings, vehicles, ships, and aircraft has results in considerable effort being directed toward developing an accurate method for analyzing the dynamic behavior of such structures. The stability behavior of elastic curved beams has been the subject of a large number of investigations. Solutions of the relevant differential equations have traditionally been obtained by the standard finite difference. These techniques require a great deal of computer time as the number of discrete nodes becomes relatively large under conditions of complex geometry and loading. One of the efficient procedures for the solution of partial differential equations is the method of differential quadrature. The differential quadrature method(DQM) has been applied to a large number of cases to overcome the difficulties of the complex algorithms of programming for the computer, as well as excessive use of storage due to conditions of complex geometry and loading. In this study, the in-plane extensional vibration for asymmetric curved beams with linearly varying cross-section is analyzed using the DQM. Fundamental frequency parameters are calculated for the member with various parameter ratios, boundary conditions, and opening angles. The results are compared with the result by other methods for cases in which they are available. According to the analysis of the solutions, the DQM, used only a limited number of grid points, gives results which agree very well with the exact ones.

• The Journal of Engineering Geology
• /
• v.24 no.2
• /
• pp.261-271
• /
• 2014

An almost permanent anchor (friction type) is resistant to ground deformation due to the friction between the soil and grout at a fixed length from the anchor body. The purpose of this study is to calculate the force of bearing resistance for a bearing anchor in enlarged boreholes. We conducted analytical and numerical analyses, along with laboratory testing, to find the quantities of bearing resistance prior to grouting in EBA (Enlarged Bearing Anchor) construction. The force of bearing resistance from the analytical method was defined as a function of general borehole diameter, expanded borehole diameter, and soil unconfined compressive strength. We also employed the Flac 3D finite difference numerical modeling code to analyze the bearing resistance of the soil conditions. We then created a laboratory experimental model to measure bearing resistance and carried out a pull-out test. The results of these three analyses are presented here, and a regression analysis was performed between bearing resistance and uniaxial compression strength. The laboratory results yield the strongest bearing resistance, with reinforcement 28.5 times greater than the uniaxial compression strength; the analytical and numerical analyses yielded values of 13.3 and 9.9, respectively. This results means that bearing resistance of laboratory test appears to be affected by skin friction resistance. To improve the reliability of these results, a comparison field study is needed to verify which results (analytical, numerical, or laboratory) best represent field observations.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.1
• /
• pp.127-140
• /
• 1988

The management of sewage and rainfall runoff becomes an emerging problem with the growth of urban communities. From the uncontrollable excess intensity or amount of rainfall, the conditions of sewer surcharge or manhole overflow could be generated in the combined sewer network where municipal or industrial wastewaters and rainfall runoff flow. The predictive model far the prevention of property and human life losses from this inundation was studied in this research. In the development of a mathematical flow model for the combined sewer surcharge and overflow, the Preissmann Slot concept and the four-point implicit method of finite difference were utilized. For the usage in personal computer, the overlapping segment method that required less memory storage was adopted. Through the simulation of hypothetical sewer network, the conservation of discharge volume was checked, and the usefulness of the Preissmann Slot was assured from the temporal distribution of discharge and depth along the sewer network. Also the possible field application for the correction of sewer diameters and slopes in the design of sewer network which has no surcharge/overflow condition was suggested.

• Nuclear Engineering and Technology
• /
• v.22 no.2
• /
• pp.116-127
• /
• 1990

The fission gas release model used In the SPEAR-BETA fuel performance code was modified by use of effective thermal conductivity for cracked fuel and by laking Into account axial fission-gas mixing between the fuel-clad gap and the plenum. With use of this modified model the fission gas release was analyzed under various power ramping conditions of P$_{max}$ and $\Delta$.fP. Effective fuel thermal conductivity that accounts for the effect of fuel tracking was used in calculation of the fuel temperature distribution and the Internal gas pressure under power ramping conditions. Mixing and dilution effects due to axial gas flow were also considered in computing the width and the thermal conductivity of the gap. The effect of axial gas flow w3s solved by the Crank-Nicholson method. The finite difference method was used to save running time in the calculation. The present modified fission-gas release model was validated by comparing its predicted results with experimental data from various lamping tests In the literature and calculated results with use of the models used In the SPEAR-BETA and FEMAXI-IV codes. Results obtained with use of the present modified model showed better agreement with experimental data reported in the literature than those results with use of the latter codes. The fuel centerline temperature calculated with introduction of effective thermal conductivity for centerline temperature calculated with Introduction of effective thermal conductivity for cracked fuel was 200 higher fission gas release predicted with use of the modified model was nearly 6% larger on the average than that calculated by use of the unmodified model used in the SPEAR-BETA code.e SPEAR-BETA code.e.

• The KIPS Transactions:PartA
• /
• v.15A no.3
• /
• pp.141-149
• /
• 2008

In this research we present a novel method which combines and visualizes the design model and the FDM-based simulation result of solidification. Moreover we employ VR displays and visualize stereoscopic images to provide an effective analysis environment. First we reconstruct the solidification simulation result to a rectangular mesh model using a conventional simulation software. Then each point color of the reconstructed model represents a temperature value of its position. Next we map the two models by finding the nearest point of the reconstructed model for each point of the design model and then assign the point color of the design model as that of the reconstructed model. Before this mapping we apply mesh subdivision because the design model is composed of minimum number of points and that makes the point distribution of the design model not uniform compared with the reconstructed model. In this process the original shape is preserved in the manner that points are added to the mesh edge which length is longer than a predefined threshold value. The implemented system visualizes the solidification simulation data on the design model, which allows the user to understand the object geometry precisely. The immersive and realistic working environment constructed with use of VR display can support the user to discover the defect occurrence faster and more effectively.

• Tunnel and Underground Space
• /
• v.14 no.5
• /
• pp.327-338
• /
• 2004

The purpose of this research is to establish a concept of the factor of safety of tunnels which is a quantitative estimate of the stability of tunnels. Based on this concept, a numerical technique which calculates the factor of safety of tunnels was developed. To obtain the safety factor of a tunnel, the strength reduction technique in which a series of analyses are repeated with reduced ground strength until the tunnel collapses were employed. With this technique, the failure plane, as well as the factor of safety, can be obtained without prescribing the trial failure plane. Analyses were conducted with FLA $C^{2D}$(ver3.3), a geotechnical analysis program which is based on the finite difference method. From the result, the location of plastic zones, the maximum convergence and the maximum stress generated in the support system were also analyzed. The result shows that factors of safety are higher for the 1st and 2nd rock classes, and lower for the lower rock classes. Furthermore, factor of safety is higher when $K_{0}$ =0.5 compared to at in case of $K_{0}$ =2.0. Through this research, it is found that the factor of safety defined in this research can be used as a good quantitative index representing the stability of tunnels. Also, close examination of the results can help adjustment of the quantity and location of additional supports.s.

• Journal of the Korean Geotechnical Society
• /
• v.32 no.4
• /
• pp.5-14
• /
• 2016

Numerical simulation of dynamic soil-pile-structure interaction embedded in a dry sand was carried out. 3D model of the dynamic centrifuge model tests was formulated in a time domain to consider nonlinear behavior of soil using the finite difference method program, FLAC3D. As a modeling methodology, Mohr-Coulomb criteria was adopted as soil constitutive model. Soil nonlinearity was considered by adopting the hysteretic damping model, and an interface model which can simulate separation and slip between soil and pile was adopted. Simplified continuum modeling (Kim et al., 2012) was used as boundary condition to reduce analysis time. Calibration process for numerical modeling results and test results was performed through the parametric study. Verification process was then performed by comparing numerical modeling results with another test results. Based on the calibration and validation procedure, it is identified that proposed modeling method can properly simulate dynamic behavior of soil-pile system in dry condition.