• Journal of Korea Water Resources Association
• /
• v.32 no.2
• /
• pp.185-195
• /
• 1999

The finite difference method and the method of characteristics are frequently used for the numerical analysis of kinematic wave model. Truncation errors cause the peak discharge dissipated in the solution from the finite difference method. The peak discharge is conserved in the solution from the finite difference method. The peak discharge is conserved in the solution from the method of characteristics, however, the shock may deteriorates the numerical solution. In this paper, distinctive features of each scheme are investigated for the numerical analysis of kinematic wave model, and applicability of shock fitting algorithm such as Propagating Shock Fitting and Approximated Shock Fitting methods are studied. Propagating Shock Fitting method appears to treat shock properly, however, it failed to fit the shock appropriately when applied to a sudden inflow change in a long river. Approximate Shock Sitting method, which uses finer elements, is found to be more proper shock-fitting than the Propagating Shock Fitting method. Comparisons are made between two solution from the kinematic wave theory with shock fitting and full dynamic wave theory, and the results are discussed.

• Geophysics and Geophysical Exploration
• /
• v.7 no.2
• /
• pp.148-154
• /
• 2004

This article reviews the development of three-dimensional (3-D) magnetotelluric (MT) modeling. The 3-D modeling of electromagnetic fields is essential in understanding the physics of MT soundings, and in implementing an inversion method to reconstruct a 3-D resistivity image. Although various numerical schemes have been developed over the last two decades, practical methods have been quite limited. However, the recent rapid improvement in computer speed and memory, as well as the advance in iterative solution algorithms for a large system of equations, makes it possible to model the MT responses of complex 3-D structures, which have been very difficult to simulate before. The use of staggered grids in finite difference method has become popular, conserving a magnetic flux and an electric current and allowing for realistic discontinuous fields. The convergence of numerical solutions has been greatly accelerated by adopting Krylov subspace methods, proper preconditioning techniques, and static divergence corrections. The vector finite-element method using edge elements is also free from the discontinuity problem, and seems a natural choice for modeling complex structures including irregular topography because its flexibility allows one to capture full geometric complexity.

• The Journal of Korean Academy of Prosthodontics
• /
• v.50 no.1
• /
• pp.36-43
• /
• 2012

Purpose: To analyze the stress distribution of the implant and its supporting structures through 3D finite elements analysis for implants with different hexagon heights and to make the assessment of the mechanical stability and the effect of the elements. Materials and methods: Infinite elements modeling with CAD data was designed. The modeling was done as follows; an external connection type ${\phi}4.0mm{\times}11.5mm$ Osstem$^{(R)}$ USII (Osstem Co., Pusan, Korea) implant system was used, the implant was planted in the mandibular first molar region with appropriate prosthetic restoration, the hexagon (implant fixture's external connection) height of 0.0, 0.7, 1.2, and 1.5 mm were applied. ABAQUS 6.4 (ABAQUS, Inc., Providence, USA) was used to calculate the stress value. The force distribution via color distribution on each experimental group's implant fixture and titanium screw was studied based on the equivalent stress (von Mises stress). The maximum stress level of each element (crown, implant screw, implant fixture, cortical bone and cancellous bone) was compared. Results: The hexagonal height of the implant with external connection had an influence on the stress distribution of the fixture, screw and upper prosthesis and the surrounding supporting bone. As the hexagon height increased, the stress was well distributed and there was a decrease in the maximum stress value. If the height of the hexagon reached over 1.2mm, there was no significant influence on the stress distribution. Conclusion: For implants with external connections, a hexagon is vital for stress distribution. As the height of the hexagon increased, the more effective stress distribution was observed.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2009.10a
• /
• pp.797-803
• /
• 2009

The accurate prediction of dynamic characteristics of high speed rotors, such as gas turbines, is important to avoid the possibility of operating the machinery near the critical speeds or unstable speed regions. However, the dynamic analysis methods and softwares for gas turbines have been developed in the process of producing many gas turbines by manufacturers and most of them have seldom been disclosed to the public. Recently, commercial FEM softwares, such as SAMCEF, ANSYS and NASTRAN, started supporting some rotordynamics analysis modules based on 3-D finite elements. In this paper, the dynamic analysis method using commercial S/W, especially ANSYS, is attempted for the small-size gas turbine engine rotor, and the analysis capability and limitations of its rotordyamics module are evaluated for further improvement of the module. As the preliminary procedure, the rotordyamic analysis capability of ANSYS was tested and evaluated with the reference models of the well-known dynamics. The limitations in application of the rotordynamics module were then identified. Under the current capability and limitations of ANSYS, it is shown that Lee diagram, a new frequency-speed diagram enhanced with the concept of $H{\infty}$ in rotating machinery, can be indirectly obtained from FRFs computed from harmonic response analysis of ANSYS. Finally, it is demonstrated based on the modeling and analysis method developed in the process of the S/W verification that the conventional Campbell diagram, Lee diagram, mode shapes and critical speeds of the small-size gas turbine engine rotor can be computed using the ANSYS rotordynamics module.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.20 no.1
• /
• pp.36-44
• /
• 2010

The accurate prediction of dynamic characteristics of high speed rotors, such as gas turbines, is important to avoid the possibility of operating the machinery near the critical speeds or unstable speed regions. However, the dynamic analysis methods and softwares for gas turbines have been developed in the process of producing many gas turbines by manufacturers and most of them have seldom been disclosed to the public. Recently, commercial FEM softwares, such as SAMCEF, ANSYS and NASTRAN, started supporting some rotordynamics analysis modules based on 3-D finite elements. In this paper, the dynamic analysis method using commercial S/W, especially ANSYS, is attempted for the small-size gas turbine engine rotor, and the analysis capability and limitations of its rotordyamics module are evaluated for further improvement of the module. As the preliminary procedure, the rotordyamic analysis capability of ANSYS was tested and evaluated with the reference models of the well-known dynamics. The limitations in application of the rotordynamics module were then identified. Under the current capability and limitations of ANSYS, it is shown that Lee diagram, a new frequency-speed diagram enhanced with the concept of $H{\infty}$ in rotating machinery, can be indirectly obtained from FRFs computed from harmonic response analysis of ANSYS. Finally, it is demonstrated based on the modeling and analysis method developed in the process of the S/W verification that the conventional Campbell diagram, Lee diagram, mode shapes and critical speeds of the small-size gas turbine engine rotor can be computed using the ANSYS rotordynamics module.

• Journal of the Institute of Electronics Engineers of Korea SD
• /
• v.39 no.1
• /
• pp.1-11
• /
• 2002

In this paper, new simple method for the calculation of depletion region under complex geometry and general purpose numerical simulator that could handle this were developed and applied in the analysis of SCM with nanoscale tip, which is a promising tool for high resolution dopant profiling. Our simple depletion region seeking algorithm alternatively switches material of elements to align ionized element boundary with contour of zero potential. To prove the validity of our method we examined whether our results satisfy the definition of depletion region and compared those with known values of un junction and MOS structure. By modeling of capacitance based on the shape of depletion region and potential distribution, we could calculate the CV curve and dC/dV curve between silicon substrate and nanoscale SCM tip.

• Tribology and Lubricants
• /
• v.30 no.5
• /
• pp.303-310
• /
• 2014

In order to reduce friction and improve reliability, researchers have applied various surface texturing methods to highly sliding machine elements such as mechanical seals and piston rings. Despite extensive theoretical research on surface texturing, previous numerical results are only applicable to isothermal and iso-viscous conditions. Because the lubricant flow pattern of textured bearing surfaces is much more complicated than that for non-textured bearings, the Navier?Stokes equation is more suitable than the Reynolds equation for the former. This study carries out a thermohydrodynamic (THD) lubrication analysis to investigate the lubrication characteristics of a single micro-dimpled parallel thrust bearing cell. The analysis involves using the continuity, Navier?Stokes, energy, temperature?viscosity relation, and heat conduction equations with the commercial computational fluid dynamics (CFD) code FLUENT. This study discretizes these equations using the finite volume method and solves them using the SIMPLE algorithm. The results include finding the streamlines, pressure and temperature distributions, and variations in the friction force and leakage for various dimple radii and depths. Increasing the dimple radius and decreasing the depth causes a recirculation flow to form because of a strong vortex, and the oil temperature greatly increases compared with the non-textured case. The present numerical scheme and results are applicable to THD analysis of various surface-textured sliding bearings and can lead to further study.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.42 no.9
• /
• pp.745-751
• /
• 2014

For the multi-nozzle propulsion, the rupture pressure of nozzle closure has an effect on the initial strain rate of ignition. Moreover, the deviation of rupture pressure for each nozzle closure leads to side forces which can disturb the attitude control of rocket. When designed, it should be considered whether nozzle closures are ruptured equally and exactly in the intented pressure. In this paper, the rupture behavior is analyzed by analytical and experimental methods for plate and "+" notched nozzle closures. The rupture pressure and deviation for operating temperature, whether notched or not and notched directions are analyzed. This paper provides a comparison between rupture pressure prediction of finite elements method which tool is Abaqus/Explicit and results of the rupture test. Jonson-Cook shear failure model which corresponds to the damage initiation criterion were used in this simulation.

• Steel and Composite Structures
• /
• v.19 no.5
• /
• pp.1115-1144
• /
• 2015

Current climate conditions along with advances in technology make further design and verification methods for structural strength and reliability of wind turbine towers imperative. Along with the growing interest for "green" energy, the wind energy sector has been developed tremendously the past decades. To this end, the improvement of wind turbine towers in terms of structural detailing and performance result in more efficient, durable and robust structures that facilitate their wider application, thus leading to energy harvesting increase. The wind tower industry is set to expand to greater heights than before and tapered steel towers with a circular cross-section are widely used as more capable of carrying heavier loads. The present study focuses on the improvement of the structural response of steel wind turbine towers, by means of internal stiffening. A thorough investigation of the contribution of stiffening rings to the overall structural behavior of the tower is being carried out. These stiffening rings are placed along the tower height to reduce local buckling phenomena, thus increasing the buckling strength of steel wind energy towers and leading the structure to a behavior closer to the one provided by the beam theory. Additionally to ring stiffeners, vertical stiffening schemes are studied to eliminate the presence of short wavelength buckles due to bending. For the purposes of this research, finite element analysis is applied in order to describe and predict in an accurate way the structural response of a model tower stiffened by internal stiffeners. Moreover, a parametric study is being performed in order to investigate the effect of the stiffeners' number to the functionality of the aforementioned stiffening systems and the improved structural behavior of the overall wind converter.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.2
• /
• pp.139-149
• /
• 1994

In this paper, a finite element technique is applied to the prediction of the wave resonance phenomena in harbors. The mild-slope equation is used with a partial reflection boundary condition introduced to model the energy dissipating effects on the solid boundary. For an efficient modeling of the radiation condition at infinity, a new infinite element is developed. The shape function of the infinite element is derived from the asymptotic behavior of the first kind of the Hankel's function in the analytical boundary series solutions. For the computational efficiency, the system matrices of the element are constructed by performing the relevant integrations in the infinite direction analytically. Comparisons with the results from experiments and other solution methods show that the present model gives fairly good results. Numerical experiments are also carried out to determine the proper distance to the infinite elements from the mouth of the halter, which directly affect the accuracy and efficiency of the solution.