• Progress in Superconductivity and Cryogenics
• /
• v.17 no.2
• /
• pp.45-49
• /
• 2015

It has been well known that a toroid is the inevitable shape for a high temperature superconducting (HTS) coil as a component of a large scale superconducting magnetic energy storage system (SMES) because it is the best option to minimize a magnetic field intensity applied perpendicularly to the HTS wires. Even though a perfect toroid coil does not have a perpendicular magnetic field, for a practical toroid coil composed of many HTS pancake coils, some type of perpendicular magnetic field cannot be avoided, which is a major cause of degradation of the HTS wires. In order to suggest an optimum design solution for an HTS SMES system, we need an accurate, fast, and effective calculation for the magnetic field, mechanical stresses, and stored energy. As a calculation method for these criteria, a numerical calculation such as an finite element method (FEM) has usually been adopted. However, a 3-dimensional FEM can involve complicated calculation and can be relatively time consuming, which leads to very inefficient iterations for an optimal design process. In this paper, we suggested an intuitive and effective way to determine the maximum magnetic field intensity in the HTS coil by using an analytic and statistical calculation method. We were able to achieve a remarkable reduction of the calculation time by using this method. The calculation results using this method for sample model coils were compared with those obtained by conventional numerical method to verify the accuracy and availability of this proposed method. After the successful substitution of this calculation method for the proposed design program, a similar method of determining the maximum mechanical stress in the HTS coil will also be studied as a future work.

• Composites Research
• /
• v.26 no.6
• /
• pp.363-372
• /
• 2013

Numerical simulation for Resin Transfer Molding manufacturing process is attempted by using the eXtended Finite Element Method (XFEM) combined with the level set method. XFEM allows to obtaining a good numerical precision of the pressure near the resin flow front, where its gradient is discontinuous. The enriched shape functions of XFEM are derived by using the level set values so as to correctly describe the interpolation with the resin flow front. In addition, the level set method is used to transport the resin flow front at each time step during the mold filling. The level set values are calculated by an implicit characteristic Galerkin FEM. The multi-frontal solver of IPSAP is adopted to solve the system. This work is validated by comparing the obtained results with analytic solutions. Moreover, a localization method of XFEM and level set method is proposed to increase the computing efficiency. The computation domain is reduced to the small region near the resin flow front. Therefore, the total computing time is strongly reduced by it. The efficiency test is made with a simple channel flow model. Several application examples are analyzed to demonstrate ability of this method.

• Structural Engineering and Mechanics
• /
• v.90 no.6
• /
• pp.577-590
• /
• 2024

Considering the Column-Supported Group Silos (CSGSs) often arranged by rows in practical applications, earthquake responses will be affected by group effect. Since group effect presenting uncertainties, establishing the analytic model and evaluating characteristics of CSGSs seems necessary. This study aimed at providing a simplified method to evaluate seismic performances of the CSGSs. Firstly, the CSGSs with different storage granule heights are used as numerical examples to derive the base shear formula for three-particle dynamic analytical model. Then, the base shear distribution coefficient is defined as the group effect index. The simplified calculation method of the group silos based on the distribution coefficients is proposed. Finally, based on the empty, half, and full granular storage conditions, the empirical design parameters for the group silos system are given by combining finite element simulation with shaking table test. The group effect of storage granule heights of group silos on its frequency and base shear are studied by comparative analysis between group silos and independent single silo. The results show that the frequency of CSGSs decreases with the increasing weight of the stored granule. The connection between the column top and silo bottom plate is vulnerable, and structural measures should be strengthened to improve its damage resistance. In case of different storage granule heights, distribution coefficients are effective to reconstruction the group effect. The complex calculations of seismic response for CSGSs can be avoided by adopting the empirical distribution coefficients obtained in this study. The proposed method provides a theoretical reference for evaluation on the seismic performances of the CSGSs.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.7 no.4
• /
• pp.192-198
• /
• 2004

A numerical method is developed to solve a two-dimensional diffraction problem for a body located in a sediment pocket where a heavier muddy water is trapped. In the present study, the wave exciting forces acting on a submerged body in the water-sediment interface by an incident wave is investigate. It is assumed that the heavier mud is trapped locally in a sediment pocket. A mathematical formulation is made in the scope of the potential theory. The fluid is assumed to be inviscid, incompressible and its motion irrotational. The boundary conditions on the unknown free surface and interface are linearized. As a method of solution, the localized finite-element method is adopted. In the method, the computation domain is reduced by utilizing the complete set of analytic solutions known in the infinite subdomain to be truncated by introduction of an appropriate juncture conditions. The main advantage of this method is that any complex geometry of the boundaries can be easily accommodated. Computations are carried out for mono-chromatic plane progressive surface waves normally incident on the domain. Numerical results are compared with those obtained by Lassiter based on Schwingers variational method. Good Agreements are obtained in general. Another numerical computations are made for the cases with and without a body in the sediment pocket.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.43-53
• /
• 1999

Numerical modeling and inversion for electromagnetic exploration methods are essential to understand behaviour of electromagnetic fields in complex subsurface. In this study, a finite element method was adopted as a numerical scheme for the 2.5-dimensional forward problem. And a finite element equation considering linear conductivity variation was proposed, when 2.5-dimensional differential equation to couple eletric and magnetic field was implemented. Model parameters were investigated for near-field with large source effects and far-field with responses dominantly by homogeneous half-space. Numerical responses by this study were compared with analytic solutions in homogeneous half-space. Blocky inversion model was modified to be applied to the forward calculation in this study and it was also adopted in the inversion algorithm. Resolution for isolated bodies were investigated to confirm possibility and limitation of inversion for electromagnetic exploration data.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.3 no.1
• /
• pp.41-54
• /
• 1999

For free vibration of non-symmetric thin-walled circular arches including restrained warping effect, the elastic strain and kinetic energy is derived by introducing displacement fields of circular arches in which all displacement parameters are defined at the centroid axis. The cubic Hermitian polynomials are utilized as shape functions for development of the curved thin-walled beam element having eight degrees of freedom. Analytical solution for in-plane free vibration behaviors of simply supported thin-walled curved beams with monosymmetric cross-sections is newly derived. Also, a finite element formulation using two noded curved beams element is presented by evaluating elastic stiffness and mass matrices. In order to illustrate the accuracy and practical usefulness of this study, analytical and numerical solutions for free vibration of circular arches are presented and compared with solutions analyzed by the straight beam element and the ABAQUS's shell element.

• Geophysics and Geophysical Exploration
• /
• v.19 no.1
• /
• pp.45-50
• /
• 2016

Generally, a point source has been routinely used in the electrical resistivity measurements because of easy installation. If steel-cased wells are used as long electrodes, we can expect the better depth of investigation. However, the resistivity data with long electrodes can not be processed with a conventional inversion algorithm because a long electrode produces the different primary potential distribution compared with the point source. In this study, we proposed a new technique to process the electrical resistivity data with long electrodes by replacing the long electrode with a sequence of point electrodes. Comparing the potentials obtained from the technique with the analytic/numerical solution, we ensure that the proposed technique can be used for the numerical resistivity modeling based on the finite difference or finite element method.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.105-115
• /
• 2000

As large underground projects such as radioactive waste disposal, hot water and heat storage, and geothermal energy become influential, the study, which consider all aspects of thermics, hydraulics and mechanics would be needed. Thermo Hydro-Mechanical coupling analysis is one of the most complex numerical technique because it should be implemented with the combined three governing equations to analyze the behavior of rock mass. In this study, finite element code, which is based on Biot's consolidation theory, was developed to analyze the thermo-hydro-mechanical coupling in continuum rock mass. To verify the implemented program, one-dimensional consolidation model under the isothermal and non-isothermal conditions was analyzed and was compared with the analytic solution. The parametric study on two-dimensional consolidation was also performed and the effects of several factors such as poisson's ratio and hydraulic anisotropy on rock mass behavior were investigated. In the future, this program would be revised to be used for analysis of general discontinuous media with incorporating discrete joint model.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.355-365
• /
• 2000

As large underground projects such as radioactive waste disposal, hot water and heat storage, and geothermal energy become influential, the study, which consider all aspects of thermics, hydraulics and mechanics would be needed. Thermo-Hydro-Mechanical coupling analysis is one of the most complex numerical technique because it should be implemented with the combined three governing equations to analyze the behavior of rock mass. In this study, finite element code, which is based on Biot's consolidation theory, was developed to analyze the thermo-hydro-mechanical coupling in continuum rock mass. To verify the implemented program, one-dimensional consolidation model under the isothermal and non-isothermal conditions was analyzed and was compared with the analytic solution. The parametric study on two-dimensional consolidation was also performed and the effects of several factors such as poisson's ratio and hydraulic anisotropy on rock mass behavior were investigated. In the future, this program would be revised to be used for analysis of general discontinuous media with incorporating discrete joint model.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.6
• /
• pp.714-725
• /
• 2003

The theory of the current-type electromagnetic flowmeter for a high temporal resolution was developed for two-phase flow measurements. To predict the output of the current-type flowmeter, the three-dimensional virtual potential distribution C and the newly introduced flow pattern coefficient f were derived and computed. The output of flowmeter depends on the liquid conductivity (sensitive to temperature) and flow configurations of the two-phase flow with the sinusoidal excitation over 100 Hz. The flow pattern coefficient was specially devised to separate the dependency on the flow configuration of the two-phase flow from that on the liquid conductivity which can be expressed with the calibration of single-phase flow. Using the finite difference method, the three-dimensional virtual potential distributions were computed for the electrode of finite size. By taking derivative of the virtual potential, the weight functions were evaluated and compared with existing analytic series solution for the point-electrode. There was a reasonable correspondence between the present and existing results. In addition, the flow pattern coefficients were evaluated for annular flows with various film thicknesses, and compared with the experimental results by the impedance spectroscopy. The numerical results agreed well with the experimental data.