• Journal of KIISE:Computer Systems and Theory
• /
• v.32 no.9
• /
• pp.445-456
• /
• 2005

We propose a fast and accurate fluid solver of the Wavier-Stokes equations for the physics-based fluid simulations. Our method utilizes the solution of the Stokes equation as an initial guess for the velocity of the nonlinear term in the Wavier-Stokes equations. By guessing the initial velocity close to the exact solution of the given nonlinear differential equations, we can develop remarkably accurate and stable fluid solver. Our solver is based on the implicit scheme of finite difference methods, that makes it work well for large time steps. Since we employ the ADI method, our solver is also fast and has a uniform computation time. The experimental results show that our solver is excellent for fluids with high Reynolds numbers such as smoke and clouds.

• Journal of electromagnetic engineering and science
• /
• v.3 no.1
• /
• pp.39-44
• /
• 2003

An unconditionally stable compact 2D Alternating-Direction Implicit (ADI) FDTD method for calculating dispersion characteristics of waveguide structures is proposed. The numerical stability and numerical dispersion relation of the proposed method are also presented and discussed. Numerical wavelengths for the dominant and higher order modes in a hollow waveguide are obtained from numerical simulations and compared with those from the analytical dispersion relation. The numerical results show that the proposed scheme has the potential to successfully analyze a class of waveguides having locally fine geometry with reduced numerical costs.

• Journal of the Korean Mathematical Society
• /
• v.55 no.3
• /
• pp.719-734
• /
• 2018

An image restoration problem with Poisson noise arises in many applications of medical imaging, astronomy, and microscopy. To overcome ill-posedness, Total Variation (TV) model is commonly used owing to edge preserving property. Since staircase artifacts are observed in restored smooth regions, higher-order TV regularization is introduced. However, sharpness of edges in the image is also attenuated. To compromise benefits of TV and higher-order TV, the weighted sum of the non-convex TV and non-convex higher order TV is used as a regularizer in the proposed variational model. The proposed model is non-convex and non-smooth, and so it is very challenging to solve the model. We propose an iterative reweighted algorithm with the proximal linearized alternating direction method of multipliers to solve the proposed model and study convergence properties of the algorithm.

• The Transactions of the Korean Institute of Electrical Engineers P
• /
• v.59 no.3
• /
• pp.245-252
• /
• 2010

This paper proposes a constrained multi-area dispatch scheduling algorithm applicable to interconnected power system operations. The dispatch scheduling formulated as an MIP problem can be efficiently computed by GBD algorithm. GBD guarantees adequate computation speed and solution convergence by reducing the dimension of the dispatch scheduling problem. In addition, the regional decomposition technique based on ADM is introduced to obtain efficient inter-temporal OPF solution. It can find the most economic dispatch schedule incorporating power transactions without each regional utility's private information open.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.66 no.7
• /
• pp.1078-1082
• /
• 2017

This paper presents an L1 optimization based filtering technique which effectively eliminates the optical misalignment effects encountered in the squint guidance mode with strapdown seekers. We formulated a series of L1 optimization problems in order to separate the bias and the gradient components from the measured data, and solved them via the alternating direction method of multipliers (ADMM) and sparse matrix decomposition techniques. The proposed technique was able to rapidly detect arbitrary discontinuities and gradient changes from the measured signals, and was shown to effectively cancel the undesirable effects coming from the seeker misalignment angles. The technique was implemented on embedded flight computers and the real-time operational performance was verified via the hardware-in-the-loop simulation (HILS) tests in parallel with the automatic target recognition algorithms and the intra-red synthetic target images.

• Journal of Electrical Engineering and Technology
• /
• v.13 no.5
• /
• pp.1778-1790
• /
• 2018

Based on the requirements and characteristics of multi-zone autonomous decision-making in modern power system, fully distributed computing methods are needed to optimize the economic dispatch (ED) problem coordination of multi-regional power system on the basis of constructing decomposition and interaction mechanism. In this paper, four fully distributed methods based on alternating direction method of multipliers (ADMM) are used for solving the ED problem in distributed manner. By duplicating variables, the 2-block classical ADMM can be directly used to solve ED problem fully distributed. The second method is employing ADMM to solve the dual problem of ED in fully distributed manner. N-block methods based on ADMM including Alternating Direction Method with Gaussian back substitution (ADM_G) and Exchange ADMM (E_ADMM) are employed also. These two methods all can solve ED problem in distributed manner. However, the former one cannot be carried out in parallel. In this paper, four fully distributed methods solve the ED problem in distributed collaborative manner. And we also discussed the difference of four algorithms from the aspects of algorithm convergence, calculation speed and parameter change. Some simulation results are reported to test the performance of these distributed algorithms in serial and parallel.

• Journal of the Institute of Electronics Engineers of Korea TC
• /
• v.39 no.10
• /
• pp.38-45
• /
• 2002

This paper presents the new Compact 2D ADI-FDTD(Alternating-Direction Implicit Finite-Difference Time-Domain) method, where the time step is no longer restricted by the numerical stability condition. This method is an accelerating algorithm for the conventional Compact 2D FDTD method. To validate this algorithm, we have analyzed the dispersion characteristics of the hollow rectangular waveguide and the shielded microstrip line. The results of the proposed method are very well agreed with those of both the conventional analytic method and the Compact 2D FDTD method. The CPU time for analysis of this method is very much reduced compared with the conventional Compact 2D FDTD method. The proposed method is valuable as a fast algorithm in the research of dispersion characteristics of the waveguide structures.

• Proceedings of the Korean Institute of Electrical and Electronic Material Engineers Conference
• /
• 2003.07a
• /
• pp.111-114
• /
• 2003

In order to improve the surface uniformity and the conduction properties of the fabricated YBCO thick films, a system that applies alternating field vertically to the EPD field has been developed for the first time and applied to the electrophoretic deposition process. The applied alternating electric field caused a force to be exerted on each YBCO particle and resulted in a shaking of the particle in the direction of applied electric field, accomplishing a uniform particle orientation. The usual commercial electrical power was used for the vertically applied alternating voltage and the induced electric field was 25-120 V/cm at 60Hz. The thick film fabricated by the method developed in this paper showed better surface uniformity without crack and porosity and improved film characteristics such as critical temperature ($T_{c,zero}$ : 90 K) and critical current density ($2354\;A/cm^2$). Therefore, it is expected that the shaky-aligned electrophoretic deposition method can be used to fabricate superconductor films through a simpler process and at less expense.

• Proceedings of the IEEK Conference
• /
• 2004.06a
• /
• pp.131-134
• /
• 2004

This paper has been studied a ADI-FDTD(Alternating Direction Implicit Finite Difference Time Domain ) algorithm using an alternating Direction time-stepping scheme for GPR( Ground-Penetrating Radar ) system simulation. We did the numerical formulations for three-dimensional ADI-FDTD algorithm and PML(Perfect Matched Layer), and made an simple experiment on a arbitrary cube with programed algorithms. And then we compared its computed results with those of conventional FDTD.

• Journal of the Korean Mathematical Society
• /
• v.61 no.2
• /
• pp.395-408
• /
• 2024

There have been numerous studies on the characteristics of the solutions of ordinary differential equations for optimization methods, including gradient descent methods and alternating direction methods of multipliers. To investigate computer simulation of ODE solutions, we need to trace pseudo-orbits by real orbits and it is called shadowing property in dynamics. In this paper, we demonstrate that the flow induced by the alternating direction methods of multipliers (ADMM) for a C² strongly convex objective function has the eventual shadowing property. For the converse, we partially answer that convexity with the eventual shadowing property guarantees a unique minimizer. In contrast, we show that the flow generated by a second-order ODE, which is related to the accelerated version of ADMM, does not have the eventual shadowing property.