• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.1
• /
• pp.1-18
• /
• 2019

In contrast to the well-developed convex splitting schemes for gradient flows of two-component system, there were few efforts on applying the convex splitting idea to gradient flows of multi-component system, such as the vector-valued Cahn-Hilliard (vCH) equation. In the case of the vCH equation, one need to consider not only the convex splitting idea but also a specific method to manage the partition of unity constraint to design an unconditionally energy stable scheme. In this paper, we propose a constrained Convex Splitting (cCS) scheme for the vCH equation, which is based on a convex splitting of the energy functional for the vCH equation under the constraint. We show analytically that the cCS scheme is mass conserving and unconditionally uniquely solvable. And it satisfies the constraint at the next time level for any time step thus is unconditionally energy stable. Numerical experiments are presented demonstrating the accuracy, energy stability, and efficiency of the proposed cCS scheme.

• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.8 no.8
• /
• pp.1818-1827
• /
• 2004

This paper carried out numerical calculation about the directivity synthesis problem in two dimension of adaptive ultrasonic transducers. Adaptive system for directivity synthesis is constructed by multi-electrode array on the surface of only one piezoelectric ceramic plate combined with optimal algorithm. In order to realize the desired directivity that is established arbitrarily, the optimal vibrational displacement is calculated by optimal algorithm (DFP method). Secondly, the optimal voltage of electrodes which is correspond to the calculated vibrational displacement is calculated used by finite element method and DFP method. And, the vibration analysis in two dimension of piezoelectric transducer carried out by means of finite element method.

• Journal of the Korea Society of Computer and Information
• /
• v.8 no.3
• /
• pp.66-72
• /
• 2003

An ordinary differential equation in analytical numerics is utilized to some applications, for example, physics, mechanical engineering, electrical engineering, thermodynamics and etc. But this equation has problems a lots to process in the real time processing by software method. This paper is proposed a systolic Arrays architecture for solving the Runge-Kutta method. it is one of method for solving an ordinary differential equation. the proposed its architecture is very high speed and regular. this hardware proposed in this paper may be part of the mathematical problem solver's tool kit in the future and may be available to many applications in the engineering.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.21 no.1
• /
• pp.1-16
• /
• 2017

The Allen-Cahn equation is solved numerically by operator splitting Fourier spectral methods. The basic idea of the operator splitting method is to decompose the original problem into sub-equations and compose the approximate solution of the original equation using the solutions of the subproblems. The purpose of this paper is to characterize higher order operator splitting schemes and propose several higher order methods. Unlike the first and the second order methods, each of the heat and the free-energy evolution operators has at least one backward evaluation in higher order methods. We investigate the effect of negative time steps on a general form of third order schemes and suggest three third order methods for better stability and accuracy. Two fourth order methods are also presented. The traveling wave solution and a spinodal decomposition problem are used to demonstrate numerical properties and the order of convergence of the proposed methods.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.58 no.2
• /
• pp.225-231
• /
• 2009

In this paper, it is shown that Carson's equation can still be applied for the calculation of the series reactance of transmission lines with no ground return current as well as the one with ground return. It is proved in the following method. First two voltage drop equations for three-phase three wire transmission line are derived, one without considering ground return and the other using Carson's equation. The impedance matrix of the two equations are different from each other. But if we put the condition of zero ground current, $I_a+I_b+I_c=0$, those two equations becomes the identical equations. Therefore even a transmission line is not grounded, its line parameters can still be obtained using the Carson's equation. It has been confused whether or not Carson's equation can be used for an ungrounded system. It is because where ever Carson's equation is shown in the book, it also says that the system has ground return current paths as a premise. It is also verified with EMTP studies on the test circuit.

• Proceedings of the KIEE Conference
• /
• 2006.07c
• /
• pp.1494-1495
• /
• 2006

A complete finite element analysis method for discharge onset process, which is governed and coupled by charge transport equation and electric field equation, was presented. The charge transport equation of first order was transformed into a second-order one by utilizing the artificial diffusion scheme. The two second-order equations were analyzed by the finite element formulation which is well-developed for second-order ones. The Fowler-Nordheim injection boundary condition was adopted for charge transport equation. After verifying the numerical results by comparing to the analytic solutions using parallel plane electrodes with one carrier system, we extended the result to blade-plane electrodes in 2D xy geometry with three carriers system. Radius of the sharp tip was taken to be 50 ${\mu}m$. When this sharp geometry was solved by utilizing the space discretizing methods, the very sharp tip was found to cause a singularity in electric field and space charge distribution around the tip. To avoid these numerical difficulties in the FEM, finer meshes, a higher order shape function, and artificial diffusion scheme were employed.

• Proceedings of the KSR Conference
• /
• 2008.11b
• /
• pp.169-178
• /
• 2008

In this paper, it is shown that Carson's equation can still be applied for the calculation of the series reactance of transmission lines with no ground return current as well as the one with ground return. It is proved in the following method. First two voltage drop equations for three-phase three wire transmission line are derived, one without considering ground return and the other using Carson's equation. The impedance matrix of the two equations are different from each other. But if we put the condition of zero ground current, $I_a+I_b+I_c=0$, those two equations becomes the identical equations. Therefore even a transmission line is not grounded, its line parameters can still be obtained using the Carson's equation. It has been confused whether or not Carson's equation can be used for an ungrounded system. It is because where ever Carson's equation is shown in the book, it also says that the system has ground return current paths as a premise. It is also verified with EMTP studies on the test circuit.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2001.09a
• /
• pp.117-124
• /
• 2001

A simplified method fur the eigenpair sensitivities of damped system with multiple eigenvalues is presented. This approach employs a reduced equation to determine the sensitivities of eigenpairs of the damped vibratory systems with multiple natural frequencies. In the proposed method, adjacent eigenvectors and orthonormal conditions are used to compute an algebraic equation whose order is (n+m)x(n+m), where n is the number of coordinates and m the number of multiplicity of multiple natural frequencies. The proposed method is an improved Lee and Jung's method which was developed previously. Two equations are used to find eigenvalue derivatives and eigenvector derivatives in Lee and Jung's method. A significant advantage of this approach over Lee and Jung's method is that one algebraic equation newly developed is enough to compute such eigenvalue derivatives and eigenvector derivatives. This method can be consistently applied to both structural systems with structural design parameters and mechanical systems with lumped design parameters. To demonstrate the theory of the proposed method and its possibilities in the case of multiple eigenvalues, the finite element model of the cantilever beam and 5-DOF mechanical system in the case of a non-proportionally damped system are considered as numerical examples. The design parameter of the cantilever beam is its height. and that of the 5-DOF mechanical system is a spring.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.3
• /
• pp.237-244
• /
• 2016

This paper presents a nonlinear Moving Least Squares(MLS) difference method for material nonlinearity problem. The MLS difference method, which employs strong formulation involving the fast derivative approximation, discretizes governing partial differential equation based on a node model. However, the conventional MLS difference method cannot explicitly handle constitutive equation since it solves solid mechanics problems by using the Navier's equation that unifies unknowns into one variable, displacement. In this study, a double derivative approximation is devised to treat the constitutive equation of inelastic material in the framework of strong formulation; in fact, it manipulates the first order derivative approximation two times. The equilibrium equation described by the divergence of stress tensor is directly discretized and is linearized by the Newton method; as a result, an iterative procedure is developed to find convergent solution. Stresses and internal variables are calculated and updated by the return mapping algorithm. Effectiveness and stability of the iterative procedure is improved by using algorithmic tangent modulus. The consistency of the double derivative approximation was shown by the reproducing property test. Also, accuracy and stability of the procedure were verified by analyzing inelastic beam under incremental tensile loading.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.106-112
• /
• 1999

In this article, combination of the FAS-FMG multi-grid method and the Krylov subspace method was presented in solving two dimensional driven-cavity flows. Three algorithms of the Krylov subspace method, CG, CGSTAB(Bi-CG Stabilized) and GMRES method were tested with MILU preconditioner. As a smoother of the pressure correction equation, the MILU-CG is recommended rather than MILU-GMRES(k) or MILU-CGSTAB, since the MILU-GMRES(k) preconditioner has too much computation on the coarse grid compared to the MILU-CG one. As for the momentum equation, relatively cheap smoother like SIP solver may be sufficient.