• Journal of electromagnetic engineering and science
• /
• v.10 no.2
• /
• pp.45-49
• /
• 2010

This paper presents an analytical formulation of the numerical reflection coefficient of the ID-FDTD scheme at the planar dielectric boundary for a TM wave incidence. The reflection coefficient is formulated in an approximate manner, and the accuracy of this method is numerically verified. The effective dielectric constant for a grid on the interface is obtained, and then reduced to that of the Yee scheme for a small cell size.

• Geomechanics and Engineering
• /
• v.16 no.1
• /
• pp.97-104
• /
• 2018

The finite difference software Flac3D is used to study the influence of tunnel burial depth, tunnel diameter and lateral pressure coefficient of original rock stress on the stress and deformation of tunnel surrounding rock under sandstone condition. The results show that the maximum shear stress, the radius of the plastic zone and the maximum displacement in the surrounding rock increase with the increase of the diameter of the tunnel. When the lateral pressure coefficient is 1, it is most favorable for surrounding rock and lining structure, with the increase or decrease of lateral pressure coefficient, the maximum principal stress, surrounding displacement and plastic zone range of surrounding rock and lining show a sharp increase trend, the plastic zone on the lining increases with the increase of buried depth.

• Journal of Korea Foundry Society
• /
• v.10 no.6
• /
• pp.495-502
• /
• 1990

A basic heat flow model has been developed to estimate the heat transfer coefficient at the casting/mold interface during squeeze casting. Based on the measured temperature profiles in squeeze casting of Al-4.5%Si alloy, heat transfer coefficients which vary with time were calculated by numerical method. The influences of the load and the amount of fraction solid on the heat transfer coefficient have also been studied. Using the calculated heat transfer coefficient two dimensional solidification analysis in the squeeze casting process was carried out by the finite difference method, and the results were in good agreement with the experiments. It may be concluded that heat flow analysis in the squeeze casting process with accurate heat transfer coefficient at the casting /mold interface is important for a proper design of cooling in die and finally for improving productivity and die life as well.

• Journal of Power System Engineering
• /
• v.11 no.1
• /
• pp.20-27
• /
• 2007

In this paper, finite element modeling was performed for vibration analysis of a rotor system installed in sunroof motor, and analysis process was developed for natural frequency and unbalance response analysis. At the same time, to reduce analysis modeling error caused by the difference between analysis and measured values, finite element model updating was conducted using an optimization algorithm, i.e. hybrid genetic algorithm and simulated annealing (HGASA) method. For this end experimental modal test was carried out and by using the measured frequency response function (FRF), model updating was performed considering both cases where core coil was removed and included. And acceptable result was obtained. Also, dynamic property coefficient of bush bearing which influences vibration response of the rotor system was estimated.

• Journal of electromagnetic engineering and science
• /
• v.18 no.3
• /
• pp.212-214
• /
• 2018

The finite-difference time-domain (FDTD) model was developed to analyze electromagnetic (EM) wave propagation in an inhomogeneous ionosphere. The EM analysis of ionosphere is complicated, owing to various propagation environments that are significantly influenced by plasma frequency, cyclotron frequency, and collision frequency. Based on the simple auxiliary differential equation (ADE) technique, we present an accurate FDTD algorithm suitable for the EM analysis of complex phenomena in the ionosphere under arbitrary-direction geomagnetic field. Numerical examples are used to validate our FDTD model in terms of the reflection coefficient of a single magnetized plasma slab. Based on the FDTD formulation developed here, we investigate EM wave propagation characteristics in the ionosphere using realistic ionospheric data for South Korea.

• Geomechanics and Engineering
• /
• v.17 no.4
• /
• pp.317-322
• /
• 2019

This paper presents the results of parametric analyses to compute the axial capacity of a suction pile using 2D and 3D finite element approaches. The study is intended to simplify the process of analyzing suction piles from 3D to 2D model. The research focuses on obtaining the coefficient to be applied into the 2D model in order to obtain results that are as close as possible to the 3D model. Two 2D models were used in the analysis, namely the plane strain and axisymmetric models. The analyses were performed using two actual offshore soil data of the North and West Java Indonesia. The study reveals that the simplification of model through 2D Finite Element is achievable by applying the appropriate coefficient to the stiffness parameters. The results show that the simplified model of the 2D FEA provides more conservative results (with the difference between 2% to 7%) than the 3D FEA.

• Ocean Systems Engineering
• /
• v.10 no.2
• /
• pp.201-226
• /
• 2020

A systematic numerical comparative study of the performance of semicircular and rectangular submerged breakwaters interacting with solitary waves is the basis of this paper. To accomplish this task, Nwogu's extended Boussinesq model equations are employed to simulate the interaction of the wave with breakwaters. The finite difference technique has been used to discretize the spatial terms while a fourth-order predictor-corrector method is employed for time discretization in our numerical model. The proposed computational scheme uses a staggered-grid system where the first-order spatial derivatives have been discretized with fourth-order accuracy. For validation purposes, five test cases are considered and numerical results have been successfully compared with the existing analytical and experimental results. The performances of the rectangular and semicircular breakwaters have been examined in terms of the wave reflection, transmission, and dissipation coefficients (RTD coefficients) denoted by K_R, K_T, K_D. The latter coefficient K_D emerges due to the non-energy conserving K_R and K_T. Our computational results and graphical illustrations show that the rectangular breakwater has higher reflection coefficients than semicircular breakwater for a fixed crest height, but as the wave height increases, the two reflection coefficients approach each other. un the other hand, the rectangular breakwater has larger dissipation coefficients compared to that of the semicircular breakwater and the difference between them increases as the height of the crest increases. However, the transmission coefficient for the semicircular breakwater is greater than that of the rectangular breakwater and the difference in their transmission coefficients increases with the crest height. Quantitatively, for rectangular breakwaters the reflection coefficients K_R are 5-15% higher while the diffusion coefficients K_D are 3-23% higher than that for the semicircular breakwaters, respectively. The transmission coefficients K_T for rectangular breakwater shows the better performance up to 2.47% than that for the semicircular breakwaters. Based on our computational results, one may conclude that the rectangular breakwater has a better overall performance than the semicircular breakwater. Although the model equations are non-dissipative, the non-energy conserving transmission and reflection coefficients due to wave-breakwater interactions lead to dissipation type contribution.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.7
• /
• pp.849-861
• /
• 1997

In this study, the effect of hole geometry of the cooling system on the flow and temperature field was numerically calculated. The finite volume method was employed to discretize the governing equation based on the non-orthogonal coordinate with non-staggered variable arrangement. The standard k-.epsilon. turbulence model was used and also the predicted results were compared with the experimental data to validate numerical modeling. The predicted results showed good agreement in all cases. To analyze the effect of the discharge coefficient for slots of different length to width, the inlet chamfering and radiusing holes were considered. The discharge coefficient was increased with increment of the chamfering ratio, radiusing ratio and slot length to width and also the effect of radiusing showed better result than chamfering in all cases. In order to analyze the difference between the predicted results with plenum region and without plenum region, the velocity profiles of jet exit region for a various flow conditions were calculated. The normal velocity components of jet exit showed big difference for the low slot length to width and high blowing rate cases. To analyze the flow phenomena injected from a row of inclined holes in a real turbine blade, three dimensional flow and temperature distribution of the region including plenum, hole and cross stream with flow conditions were numerically calculated. The results have shown three-dimensional flow characteristics, such as the development of counter rotating vortices, jetting effect and low momentum region within the hole in addition to counter rotating vortex structure in the cross stream.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.26 no.1
• /
• pp.58-66
• /
• 2022

In this study, to examine the effect of the transverse thermal expansion behavior of FRP reinforcing bars and concrete on the concrete cover thickness, based on 20℃, when the temperature changes from -70℃ to 80℃, the behavior of concrete was studied theoretically and numerically. Theoretical elastic analysis and nonlinear finite element analysis were performed on FRP reinforced concrete with different diameters and cover thicknesses of FRP reinforcement. As a result, at a negative temperature difference, concrete was compressed, and the theoretical strain result and the finite element result were similar, but at a positive temperature difference, tensile stress and further cracks occurred in the concrete, which was 1.2 to 1.4 times larger than the theoretical result. The ratio of the diameter of the FRP reinforcing bar to the thickness of the concrete cover (c/db) is closely related to the occurrence of cracks. Since the transverse thermal expansion coefficient of FRP reinforcing bars is three times greater than that of concrete, it is necessary to consider this in design.

• Geomechanics and Engineering
• /
• v.18 no.3
• /
• pp.225-234
• /
• 2019

A numerical stepwise approach for cavity expansion problem in strain-softening rock or soil mass is investigated, which is compatible with Mohr-Coulomb and generalized Hoek-Brown failure criteria. Based on finite difference method, plastic region is divided into a finite number of concentric rings whose thicknesses are determined internally to satisfy the equilibrium and compatibility equations, the material parameters of the rock or soil mass are assumed to be the same in each ring. For the strain-softening behavior, the strength parameters are assumed to be a linear function of deviatoric plastic strain (${\gamma}p^*$) for each ring. Increments of stress and strain for each ring are calculated with the finite difference method. Assumptions of large-strain for soil mass and small-strain for rock mass are adopted, respectively. A new numerical stepwise approach for limited pressure and plastic radius are obtained. Comparisons are conducted to validate the correctness of the proposed approach with Vesic's solution (1972). The results show that the perfectly elasto-plastic model may underestimate the displacement and stresses in cavity expansion than strain-softening coefficient considered. The results of limit expansion pressure based on the generalised H-B failure criterion are less than those obtained based on the M-C failure criterion.