• Journal of the Optical Society of Korea
• /
• v.13 no.2
• /
• pp.286-293
• /
• 2009

In this paper, we present an optimal shape design method for a dielectric microlens which is used to focus an incoming infrared plane wave in wideband, by exploiting the finite difference time domain (FDTD) technique and the topology optimization technique. Topology optimization is a scheme to search an optimal shape by adjusting the material properties, which are design variables, within the design space. And by introducing the adjoint variable method, we can effectively calculate a derivative of the objective function with respect to the design variable. To verify the proposed method, a shape design problem of a dielectric microlens is tested when illuminated by a transverse electric (TE)-polarized infrared plane wave. In this problem, the design variable is the dielectric constant within the design space of a dielectric microlens. The design objective is to maximally focus the incoming magnetic field at a specific point in wideband.

• Proceedings of the KIEE Conference
• /
• 2002.11d
• /
• pp.124-127
• /
• 2002

This paper presents a 3D shape optimization algorithm for electromagnetic devices using the design sensitivity analysis with finite element method. The structural deformation analysis based on the deformation theory of the elastic body under stress is used for mesh renewing. The design sensitivity and adjoint variable formulae are derived for the 3D nonlinear finite element method with edge element. The proposed algorithm is applied to the shape optimization of 3D electromagnet to get a uniform flux density at the air gap.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.2 no.1
• /
• pp.27-39
• /
• 1998

We describe a fast numerical method for non-separable elliptic equations in self-adjoin form on irregular adaptive domains. One of the most successful results in numerical PDE is developing rapid elliptic solvers for separable EPDEs, for example, Fourier transformation methods for Poisson problem on a square, however, it is known that there is no rapid elliptic solvers capable of solving a general nonseparable problems. It is the purpose of this paper to present an iterative solver for linear EPDEs in self-adjoint form. The scheme discussed in this paper solves a given non-separable equation using a sequence of solutions of Poisson equations, therefore, the most important key for such a method is having a good Poison solver. High performance is achieved by using a fast high-order adaptive Poisson solver which requires only about 500 floating point operations per gridpoint in order to obtain machine precision for both the computed solution and its partial derivatives. A few numerical examples have been presented.

• Journal of Magnetics
• /
• v.18 no.3
• /
• pp.283-288
• /
• 2013

This paper presents structural topology optimization that is being applied for the design of electromagnetic vibration energy harvester. The design goal is to maximize the root-mean-square value of output voltage generated by external vibration leading structures. To calculate the output voltage, the magnetic field analysis is performed by using the finite element method, and the obtained magnetic flux linkage is interpolated by using Lagrange polynomials. To achieve the design goal, permanent magnet is designed by using topology optimization. The analytical design sensitivity is derived from the adjoint variable method, and the formulated optimization problem is solved through the method of moving asymptotes (MMA). As optimization results, the optimal location and shape of the permanent magnet are provided when the magnetization direction is fixed. In addition, the optimization results including the design of magnetization direction are provided.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.1
• /
• pp.1-8
• /
• 2014

Using the level set and the meshfree methods, we develop a topological shape optimization method applied to linear elasticity problems. Design gradients are computed using an efficient adjoint design sensitivity analysis(DSA) method. The boundaries are represented by an implicit moving boundary(IMB) embedded in the level set function obtainable from the "Hamilton-Jacobi type" equation with the "Up-wind scheme". Then, using the implicit function, explicit boundaries are generated to obtain the response and sensitivity of the structures. Global nodal shape function derived on a basis of the reproducing kernel(RK) method is employed to discretize the displacement field in the governing continuum equation. Thus, the material points can be located everywhere in the continuum domain, which enables to generate the explicit boundaries and leads to a precise design result. The developed method defines a Lagrangian functional for the constrained optimization. It minimizes the compliance, satisfying the constraint of allowable volume through the variations of boundary. During the optimization, the velocity to integrate the Hamilton-Jacobi equation is obtained from the optimality condition for the Lagrangian functional. Compared with the conventional shape optimization method, the developed one can easily represent the topological shape variations.

• Journal of Soil and Groundwater Environment
• /
• v.11 no.1
• /
• pp.37-44
• /
• 2006

This study summarizes advantages and disadvantages of numerical methods and compares ELLAM and LEZOOMPC to develop an efficient numerical modeling technique on contaminant transport. Eulerian-Lagrangian method and Eulerian method are commonly used numerical techniques. However Eulerian-Lagrangian method does not conserve mass globally and fails to treat boundary in a straightforward manner. Also, Eulerian method has restrictions on the size of Courant number and mesh Peclet number because of time truncation error. ELLAM (Eulerian Lagrangian Localized Adjoint Method) which has been popularly used for past 10 years in numerical modeling, is known for overcoming these numerical problems of Eulerian-Lagrangian method and Eulerian method. However, this study investigates advantages and disadvantages of ELLAM and suggests a change for the better. To figure out the disadvantages of ELLAM, the results of ELLAM, LEZOOMPC (Lagrangian-Eulerian ZOOMing Peak and valley Capturing), and visual MODFLOW are compared for four examples having different mesh Peclet numbers. The result of ELLAM generates numerical oscillation at infinite of mesh Peclet number, but that of LEZOOMPC yields accurate simulations. The simulation results suggest that the numerical error of ELLAM could be alleviated by adopting some schemes in LEZOOMPC. In other words, the numerical model which combines ELLAM with backward particle tracking, forward particle tracking, adaptively local zooming, and peak/valley capturing of LEZOOMPC can be developed for not only overcoming the numerical error of ELLAM, but also keeping the numerical advantage of ELLAM.

• Journal of Energy Engineering
• /
• v.16 no.4
• /
• pp.155-163
• /
• 2007

To insure nuclear reactor safety, the reactivity of control rods should be calculated by measuring the criticality of reactor core and it is regularly performed during the annual physics test period. Also, the core criticality should be monitored during the start-up operation to avoid reactivity induced accidents. Many research works on control rod reactivity measurement and subcriticality measurement have been accomplished throughout the world for decades and recently a new method named "Modified Neutron Source Multiplication Method (MNSM)" was proposed in Japan which is known to be improved overcoming limitations of traditional Neutron Source Multiplication Method (NSM). In this study, MNSM was tested in calculation of subcriticalities and in evaluation of application validity using the educational reactor in Kyung Hee University, AGN-201. For this study, a revised nuclear data library and a neutron transport code system TRANSX - PARTISN were established. Correction factors for various control rod positions were produced using the k-effective values and the corresponding flux distributions and adjoint flux distributions. Experimental values of the core criticality were obtained using the neutron count rates of the BF3 proportional counters. The results showed that the expected reactivity worth of control rods by MNSM agreed well with the theoretical values and the correction factors contributed much for this purpose.

• Proceedings of the KIEE Conference
• /
• 2009.07a
• /
• pp.658_659
• /
• 2009

This paper deals with the numerical implementation of the material sensitivity analysis, which is used to efficiently determine an optimal electric current distribution on a ship hull by corrosion of ship. A material sensitivity formula for the forward problem formulated in terms of the equivalent current method is analytically derived. Then, the components of the adjoint system including the electric pseudo-source are thoroughly investigated in order to obtain the first-order gradient information of an objective function with respect to the electric current. The purposed method has been successfully applied to a model ship and the predicted result on the underwater electric field due to corrosion of ship have been compared to that computed by FNREMUS software.

• Proceedings of the KIEE Conference
• /
• 2001.07b
• /
• pp.768-770
• /
• 2001

The topology optimizations of electromagnetic systems with the nonlinear and the eddy current effects are studied using the finite element method. The topology design sensitivity formulations of nonlinear magnetostatics and eddy current systems are derived using the adjoint variable method. A computer program is developed using object orient programming and applied to the topology optimization of a C-core actuator. A numerical study shows the effects of saturation and eddy current by comparing results of topology optimizations.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.100-110
• /
• 2007

A new design approach for a delicate treatment of complex geometries such as a wing/body configuration is arranged using overset mesh technique under large scale computing environment for turbulent viscous flow. Various pre- and post-processing techniques which are required of overset flow analysis and sensitivity analysis codes are discussed for design optimization problems based on gradient based optimization method (GBOM). The overset flow analysis code is validated by comparing with the experimental data of a wing/body configuration (DLR-F4) from the 1st Drag Prediction Workshop (DPW-I). In order to examine the applicability of the present design tools, careful design works for the drag minimization problem of a wing/body configuration are carried out by using the developed aerodynamic shape optimization tools for the viscous flow over multiple-body aircraft geometries.