• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.15 no.10 s.89
• /
• pp.977-982
• /
• 2004

In this paper, a dispersive anisotropic perfectly matched layer(APML) is proposed using piecewise linear recursive convolution(PLRC) for finite difference time domain(FDTD) methods. This proposed APML can be utilized for the analysis of a nonlinear dispersive medium as absorbing boundary condition(ABC). The formulation is simple modification to the original AMPL and can be easily implemented. Also it has some advantages of the PLRC approach-fast speed, low memory cost, and easy formulation of multiple pole susceptibility. We applied this APML to 2-D propagation problems in dispersive media such as Debye and Lorentz media The results showed good absorption at boundaries.

• Smart Structures and Systems
• /
• v.28 no.6
• /
• pp.745-760
• /
• 2021

Dynamic buckling of structure is one of the failure modes that needs to be considered since it may result in catastrophic failure of the structure in a short period of time. For a thin fiber-reinforced polymer (FRP) plate under compression, buckling is an inherent hazard which will be intensified by the existence of defects like holes, cracks, and delamination. On the other hand, the growth of the delamination is another prime concern for thin FRP plates. In the current paper, reinforcing the plates against buckling is realized by using SMA wires in the form of stitches. A numerical framework is proposed to simulate the dynamic instability emphasizing the effect of the SMA stitches in suppressing delamination growth. The suggested algorithm is more accurate than the other methods when considering the transformation point of the SMA wires and the modeling of the cohesive zone using simple and yet reliable technique. The computational design of the method by producing the line by line orders leads to a simple algorithm for simulating the super-elastic behavior. The Lagoudas constitutive model of the SMA material is implemented in the form of user material subroutines (VUMAT). The normal bilinear spring model is used to reproduce the cohesive zone behavior. The nonlinear finite element formulation is programmed into FORTRAN using the Newmark-beta numerical time-integration approach. The obtained results are compared with the results obtained by the finite element method using ABAQUS/Explicit solver. The obtained results by the proposed algorithm and those by ABAQUS are in good agreement.

• Journal of Mechanical Science and Technology
• /
• v.20 no.1
• /
• pp.110-124
• /
• 2006

The multi scale wavelet-Galerkin method implemented in an adaptive manner has an advantage of obtaining accurate solutions with a substantially reduced number of interpolation points. The method is becoming popular, but its numerical efficiency still needs improvement. The objectives of this investigation are to present a new numerical scheme to improve the performance of the multi scale adaptive wavelet-Galerkin method and to give detailed implementation procedure. Specifically, the subdomain technique suitable for multiscale methods is developed and implemented. When the standard wavelet-Galerkin method is implemented without domain subdivision, the interaction between very long scale wavelets and very short scale wavelets leads to a poorly-sparse system matrix, which considerably worsens numerical efficiency for large-sized problems. The performance of the developed strategy is checked in terms of numerical costs such as the CPU time and memory size. Since the detailed implementation procedure including preprocessing and stiffness matrix construction is given, researchers having experiences in standard finite element implementation may be able to extend the multi scale method further or utilize some features of the multiscale method in their own applications.

• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.52-61
• /
• 2014

The present paper deals with the efficient computation of higher-order CFD methods for compressible flow using graphics processing units (GPU). The higher-order CFD methods, such as discontinuous Galerkin (DG) methods and correction procedure via reconstruction (CPR) methods, can realize arbitrary higher-order accuracy with compact stencil on unstructured mesh. However, they require much more computational costs compared to the widely used finite volume methods (FVM). Graphics processing unit, consisting of hundreds or thousands small cores, is apt to massive parallel computations of compressible flow based on the higher-order CFD methods and can reduce computational time greatly. Higher-order multi-dimensional limiting process (MLP) is applied for the robust control of numerical oscillations around shock discontinuity and implemented efficiently on GPU. The program is written and optimized in CUDA library offered from NVIDIA. The whole algorithms are implemented to guarantee accurate and efficient computations for parallel programming on shared-memory model of GPU. The extensive numerical experiments validates that the GPU successfully accelerates computing compressible flow using higher-order method.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.2
• /
• pp.177-189
• /
• 2015

This paper presents an all-direct domain decomposition approach for large-scale structural analysis. The proposed approach achieves computational robustness and efficiency by enforcing the compatibility of the displacement field across the sub-domain boundaries via local Lagrange multipliers and augmented Lagrangian formulation (ALF). The proposed domain decomposition approach was compared to the existing FETI approach in terms of the computational time and memory usage. The parallel implementation of the proposed algorithm was described in detail. Finally, a preliminary validation was attempted for the proposed approach, and the numerical results of two- and three-dimensional problems were compared to those obtained through a dual-primal FETI approach. The results indicate an improvement in the performance as a result of the implementing the proposed approach.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• v.10 no.9
• /
• pp.4307-4325
• /
• 2016

With the aim of solving the problems of long processing times, high memory consumption and low event throughput in the current processing approaches in out-of-order RFID event streams, an efficient complex event processing method based on INFA-HTS (Improved Nondeterministic Finite Automaton-Hash Table Structure) is presented in this paper. The contribution of this paper lies in the fact that we use INFA and HTS to successfully realize the detection of complex events for out-of-order RFID event streams. Specifically, in our scheme, to detect the disorder of out-of-order event streams, we expand the traditional NFA model into a new INFA model to capture the related RFID primitive events from the out-of-order event stream. To high-efficiently manage the large intermediate capturing results, we use the HTS to store and process them. As a result, these problems in the existing methods can be effectively solved by our scheme. The simulation results of our experiments show that our proposed method in this paper outperforms some of the current general processing approaches used to process out-of-order RFID event streams.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.11a
• /
• pp.840-846
• /
• 2007

This paper investigates the dynamic behavior of a HDD spindle system due to the imperfect roundness of a rotating shaft. The shaft of a spindle motor rotates with eccentricity by the unbalanced mass of the rotating part. The eccentricity generates the run-out of a spindle motor which results in the eccentric motion of a rotating part. Roundness of a shaft affects this motion which limits the memory capacity of a HDD. This research proposes a modified Reynolds equation for the coupled journal and thrust FDBs to include the variable film thickness due to the roundness. Finite element method is used to solve the Reynolds equation for the pressure distribution. Reaction forces and friction torque are obtained by integrating the pressure and shear stress, respectively. The dynamic behavior is determined by solving the equations of a motion of a HDD spindle system in six degrees of freedom with the Runge-Kutta method to characterize the motion of a rotating part. This research shows that the roundness of a rotating shaft causes the excitation frequency with integer multiple of a rotating frequency.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2002.05a
• /
• pp.348-352
• /
• 2002

a basis such as IT(Information Technology), NT(Nano Technology) and BT(Bio Technology). Recently, NT is applied to various fields that are composed of science, industry, media and semiconductor-micro technology. It has need of IT that is ultra-precision positioning technology with strokes of many hundreds mm and maintenance of nm precision in fields of ultra micro process, ultra precision measurement, photo communication part and photo magnetic memory. This thesis represents optimal design on ultra-precision positioning with single plane X-Y stage and development of artificial control system for adequacy of industrial demand. Also, dynamic simulation on global stage is performed by using ADAMS (Automated Dynamic Analysis of Mechanical System) for the purpose of grasping dynamic characteristic on user designed X-Y global stage. The error between displacements from micro stage and from FEM(Finite Element Method) is 3.53% by verifications of stability on micro stage and control performance. As maximum Von-mises stress on hinge of micro stage is 5.981kg/mm$^2$ that is 1.5% of yield stress, stability on hinge is secured. Preparing previous results, optimal design of micro stage can be possible, and reliance of results with FEM can be secured.

• Journal of computational fluids engineering
• /
• v.17 no.2
• /
• pp.35-41
• /
• 2012

This paper presents incompressible Navier-Stokes solution algorithm for 2D Free-surface flow problems on the Cartesian mesh, which was implemented to run on Graphics Processing Units(GPU). The INS solver utilizes the variable arrangement on the Cartesian mesh, Finite Volume discretization along Constrained Interpolation Profile-Conservative Semi-Lagrangian(CIP-CSL). Solution procedure of incompressible Navier-Stokes equations for free-surface flow takes considerable amount of computation time and memory space even in modern multi-core computing architecture based on Central Processing Units(CPUs). By the recent development of computer architecture technology, Graphics Processing Unit(GPU)'s scientific computing performance outperforms that of CPU's. This paper focus on the utilization of GPU's high performance computing capability, and presents an efficient solution algorithm for free surface flow simulation. The performance of the GPU implementations with double precision accuracy is compared to that of the CPU code using an representative free-surface flow problem, namely. dam-break problem.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.9
• /
• pp.2173-2180
• /
• 1995

This paper presents free vibration analysis method using the characteristics of the rotationally periodic structures and includes the analysis results of a tire as an example. The normal modes of the rotationally periodic structures are the kind of standing waves, so all sectors have the same deflection shapes, and only different phases. This property makes it possible to derive the analysis method called sector method. The sector method can give the accurate natural frequencies and the corresponding mode shapes of the rotationally periodic structure with information of only one sector. When the free vibration analysis is performed to find the dynamic characteristics of the rotationally periodic structure by using the sector method, the computer memory spaces and the CPU times can be saved. We obtained much economic benefits by using the sector method in the analysis of dynamic characteristics of a tire made of non-linear materials.