• Journal of the Optical Society of Korea
• /
• v.15 no.2
• /
• pp.196-201
• /
• 2011

Near-field properties of light emanated from a subwavelength double slit of finite length in a thin metal film, which is essential for understanding fundamental physical mechanisms for near-field optical beam manipulations and various potential nanophotonic device applications, is investigated by using a three-dimensional finite-difference time-domain method. Near-field intensity distribution along the propagation direction of light after passing through the slit has been obtained from the phase relation of transverse electric and magnetic fields and the wave impedance. It is found that the near field of emerged light from the both slits is evanescent, that is consistent with conventional surface plasmon localization near the metal surface. Due to the finite of the slit, the amplitude of this evanescent field does not monotonically approach to than of the infinite slit as the slit length increases, i.e. the near-field of the longer slit along the center line can be weaker than that of the shorter one.

• Korean Journal of Computational Design and Engineering
• /
• v.12 no.4
• /
• pp.255-262
• /
• 2007

In the finite element analysis of forming process, objects are described with a finite number of elements and nodes and the approximated solutions can be obtained by the variational principle. One of the shortcomings of a finite element analysis is that the structure of mesh has become inefficient and unusable because discretization error increases as deformation proceeds due to severe distortion of elements. If the state of current mesh satisfies a certain remeshing criterion, analysis is stopped instantly and resumed with a reconstructed mesh. In the study, a new remeshing algorithm using tetrahedral elements has been developed, which is adapted to the desired mesh density. In order to reduce the discretization error, desired mesh sizes in each lesion of the workpiece are calculated using the Zinkiewicz and Zhu's a-posteriori error estimation scheme. The pre-constructed mesh is constructed based on the modified point insertion technique which is adapted to the density function. The object domain is divided into uniformly-sized sub-domains and the numbers of nodes in each sub-domain are redistributed, respectively. After finishing the redistribution process of nodes, a tetrahedral mesh is reconstructed with the redistributed nodes, which is adapted to the density map and resulting in good mesh quality. A goodness and adaptability of the constructed mesh is verified with a testing measure. The proposed remeshing technique is applied to the finite element analyses of forging processes.

• Structural Engineering and Mechanics
• /
• v.14 no.1
• /
• pp.99-118
• /
• 2002

The scaled-boundary finite element method is a novel semi-analytical technique, combining the advantages of the finite element and the boundary element methods with unique properties of its own. The method works by weakening the governing differential equations in one coordinate direction through the introduction of shape functions, then solving the weakened equations analytically in the other (radial) coordinate direction. These coordinate directions are defined by the geometry of the domain and a scaling centre. This paper presents a general development of the scaled boundary finite-element method for two-dimensional problems where two boundaries of the solution domain are similar. Unlike three-dimensional and axisymmetric problems of the same type, the use of logarithmic solutions of the weakened differential equations is found to be necessary. The accuracy and efficiency of the procedure is demonstrated through two examples. The first of these examples uses the standard finite element method to provide a comparable solution, while the second combines both solution techniques in a single analysis. One significant application of the new technique is the generation of transition super-elements requiring few degrees of freedom that can connect two regions of vastly different levels of discretisation.

• Geophysics and Geophysical Exploration
• /
• v.7 no.4
• /
• pp.234-244
• /
• 2004

In order to calculate 3-dimensional wavefield using finite-element method in frequency domain, we must factor so huge sparse impedance matrix. Because of difficulties of handling of this huge impedance matrix, 3-dimensional wave equation modeling is conducted mainly in time domain. In this study, we simulate the 3-D wavefield using finite-element method in Laplace domain by combining high-performance parallel finite-element solver and SWEET (Suppressed Wave Equation Estimation of Traveltime) algorithm which can calculate the traveltime and the amplitude. To verify this combination, we applied it to the SEG/EAGE 3D salt model in serial and parallel computing environments.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.8
• /
• pp.1253-1260
• /
• 2001

Transient thermal analysis of a three-dimensional axisymmetric automotive disk brake is presented in this paper. Temperature fields are obtained using a hybrid FFT-FEM scheme that combines Fourier transform techniques and finite element method. The use of a fast Fourier transform algorithm can avoid singularity problems and lead to inexpensive computing time. The transformed problem is solved with finite element scheme for each frequency domain. Inverse transforms are then performed for time domain solution. Numerical examples are presented for validation tests. Comparisons with analytical results show very good agreement. Also, a 3-D simulation, based upon an automotive brake disk model is performed.

• Bulletin of the Korean Mathematical Society
• /
• v.53 no.3
• /
• pp.853-874
• /
• 2016

Based on a particular overlapping domain decomposition technique, a parallel finite element discretization algorithm for the generalized Stokes equations is proposed and investigated. In this algorithm, each processor computes a local approximate solution in its own subdomain by solving a global problem on a mesh that is fine around its own subdomain and coarse elsewhere, and hence avoids communication with other processors in the process of computations. This algorithm has low communication complexity. It only requires the application of an existing sequential solver on the global meshes associated with each subdomain, and hence can reuse existing sequential software. Numerical results are given to demonstrate the effectiveness of the parallel algorithm.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1999.05a
• /
• pp.5-8
• /
• 1999

The present work is concerned with three-dimensional finite element analysis of the hollow section extrusion process using a porthole die. For economic computation, mismatching refinement, an efficient domain decomposition method with different mesh density for each subdomain, is implemented. The proposed method improves the computational efficiency significantly, especially fur the three-dimensional analysis of extrusion problems. As a numerical example, extrusion of the underframe part of a railroad vehicle are analyzed. For three-dimensional mesh generation of a complicated shape with hexahedral elements, a modified grid-based approach with the surface element layer is utilized. The analysis results are then successfully reflected on the industrial porthole die design.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.84-85
• /
• 2003

A parallelized FEM code based on domain decomposition method has been recently developed for a large scale computational fluid dynamics. A 4-step splitting finite element algorithm is adopted for unsteady computation of the incompressible Navier-Stokes equation, and Smagorinsky LES(Large Eddy Simulation) model is chosen for turbulent flow computation. Both METIS and MPI library are used for domain partitioning and data communication between processors respectively. Tiburon of Hyundai-motor is chosen as the computational model at $Re=7.5{\times}10^{5}$, which is based on the car height. It is confirmed that the drag under road condition is smaller than that of wind tunnel condition.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1621-1626
• /
• 2003

This study deals with the heat transfer analysis on phase change optical disc with land/groove recording by means of numerical method. Finite difference time domain(FDTD) method was used to obtain the amount of absorption of light propagating inside disc and finite difference element(FEM) method was used to calculate the temperature distribution. The calculated results present the detailed information of recording characteristics on the phase change optical disc. The temperature profiles are quite different between the land track and the groove track. The recorded mark shape on land track is smaller and more elliptic than that on groove track. It is shown that the thermal problem to the neighboring track takes place due to secondary peaks. It is found that the different write strategy should be applied to land and groove recording, respectively.

• Proceedings of the KIEE Conference
• /
• 1999.11b
• /
• pp.450-454
• /
• 1999

In order to evaluate the integrity of nuclear power plants, J-integral calculation is crucial. For this purpose, finite element method is popularly used to obtain J-integral. However, high cost time consuming preprocess should be performed to design the finite element model of a cracked structure. Also, the J-integral should be verified by alternative method since it may differ depending on the calculation method. The objective of this paper is to develop a three-dimensional elastic-plastic J-integral analysis system which is named as EPAS. The EPAS program consists of an automatic mesh generator for a through-wall crack and a surface crack, a solver based on ABAQUS program, and a J-integral calculation program which provides DI(Domain Integral) and EDI(Equivalent Domain Integral) based J-integral calculation. Using the EPAS program, an optimized finite element model for a cracked structure can be generated and corresponding J-integral can be obtained subsequently.