• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.4
• /
• pp.191-197
• /
• 2021

To calculate the induced charge of atoms in molecular dynamics, linear equations for the induced charges need to be solved. As induced charges are determined at each time step, the process involves considerable computational costs. Hence, an efficient method for calculating the induced charge distribution is required when analyzing large systems. This paper introduces the Uzawa method for solving saddle point problems, which occur in linear systems, for the solution of the Lagrange equation with constraints. We apply the accelerated Uzawa algorithm, which reduces computational costs noticeably using the Schur complement and preconditioned conjugate gradient methods, in order to overcome the drawback of the Uzawa parameter, which affects the convergence speed, and increase the efficiency of the matrix operation. Numerical models of molecular dynamics in which two gold nanoparticles are placed under external electric fields reveal that the proposed method provides improved results in terms of both convergence and efficiency. The computational cost was reduced by approximately 1/10 compared to that for the Gaussian elimination method, and fast convergence of the conjugate gradient, as compared to the basic Uzawa method, was verified.

• Journal of the Korean Mathematical Society
• /
• v.44 no.2
• /
• pp.455-466
• /
• 2007

The zero-distribution of the Fourier integral $${\int}^{\infty}_{-{\infty}}\;Q(u)e^{p(u)+^{izu}du$$, where P is a polynomial with leading term $-u^{2m}(m\;{\geq}\;1)$ and Q an arbitrary polynomial, is described. To this end, an asymptotic formula for the integral is established by applying the saddle point method.

• Journal of the Korean Mathematical Society
• /
• v.49 no.4
• /
• pp.843-854
• /
• 2012

In this paper, the iterative Uzawa method with a modified Lagrangian functional is investigated to seek a saddle point for the semicoercive variational Signorini inequality.

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

To understand the concept of dynamic motion in two degree nonlinear coupled system, free vibration not including damping and excitation is investigated with the concept of nonlinear normal mode. Stability analysis of a coupled system is conducted, and the theoretical analysis performed for the bifurcation phenomenon in the system. Bifurcation point is estimated using harmonic balance method. When the bifurcation occurs, the saddle point is always found on Poincare's map. Nonlinear phenomenon result in amplitude modulation near the saddle point and the internal resonance in the system making continuous interchange of energy. If the bifurcation in the normal mode is local, the motion remains stable for a long time even when the total energy is increased in the system. On the other hand, if the bifurcation is global, the motion in the normal mode disappears into the chaos range as the range becomes gradually large.

• Proceedings of the Korean Institute of Information and Commucation Sciences Conference
• /
• 2008.05a
• /
• pp.677-680
• /
• 2008

The paper deals with the analysis of radiation characteristics of antenna in the multi-layered media structures. The dyadic Green's function for three layer medium is complex because the Green's functions belonging to the kernel of the integral equation are expressed as Sommerfeld integrals, in which surface wave effects are automatically included. When certain condition are met, the integral can be evaluated approximated by the method of Saddle-point integration. In this study, we propose a method to calculate a radiation pattern for several antennas by using the method of Saddle-point integration. Numerical results show how the radiation characteristics are affected by parameter of dielectric media.

• The Transactions of the Korean Institute of Electrical Engineers A
• /
• v.48 no.8
• /
• pp.958-964
• /
• 1999

The distance in load parameter space to the closest saddle node bifurcation (CSNB) point provides the worst case power margin to voltage instability and the left eigenvector at CSNB identifies the most effective direction to steer the system to maximize voltage stability under contingency. This paper presents an improved direct method for computing CSNB: the order of nonlinear systems equations is reduced to about twice of the size of load flow equations in contrast to about three-times in Dobson's direct method; the initial guess for the direct method is computed efficiently and robustly by combined use of continuation power flow, a pair of multiple load flow solution with Lagrange interpolation. It is also shown that voltage stability may be enhanced significantly with shift of generations in the direction of the left eigenvector at CSNB.

• Journal of the Korean Mathematical Society
• /
• v.52 no.6
• /
• pp.1195-1207
• /
• 2015

Lagrange multiplier method for solving the contact problem in elasticity is considered. Based on lower semicontinuity of sensitivity functional we prove the convergence of modified dual scheme to corresponding saddle point.

• Proceedings of the ESK Conference
• /
• 1993.04a
• /
• pp.145-145
• /
• 1993

In bicycle design, saddle is the major part which effects the subjective comfort of rider. This study examines the correlation between the subjective comfort and elastic body pressure. The elastic body pressure measuring instrument for free-form surface such as saddle is developed by force sensor register, A/D converter and computer. The subjective comfort is measured quantit- atively by 11-point scale method and the elastic body pressure distribution is obtained through 3 different saddles at 4 postures. The pressure distribution is presented by computerized equi- pressure contour. While mean pressure, standard deviation of pressure, maximum pressure are inversely proportional to subjective comfort, the modified saddle-bearing weight which is the surface integral of pressure is directly proportional. Consequently, standard deviation of pressure is most important characteristic which affects variation of subjective comfort.

• Proceedings of the KIEE Conference
• /
• 1997.11a
• /
• pp.231-234
• /
• 1997

This paper presents improved direct method for calculating the closest saddle node bifurcation (CSNB) point, which is also applicable to the selection of appropriate load shedding, reactive power compensation point detection. The proposed method reduced dimension of nonlinear equation compared with that of Dobson's direct method. The improved direct method, utilizing Newton Iterative method converges very quickly. But it diverges if the initial guess is not very close to CSNB. So the direct method is performed with the initial values obtained by carrying out the iterative method twice, which is considered most efficient at this time. Since sparsity techniques can be employed, this method is a good choice to a large scale system on-line application. Proposed method has been tested for 5-bus, New England 30-bus system.

• Proceedings of the Korea Concrete Institute Conference
• /
• 1998.10a
• /
• pp.405-411
• /
• 1998

An iterative numerical computational algorithm is presented to design a plate of shell element subjected to membrane and flexural forces. Based on equilibrium consideration, equations for capacities of top and bottom reinforcements in two orthogonal directions have been derived. The amount of reinforcement is determined locally, i. e., for each sampling point, from the equilibrium between applied and internal forces. One case of design is performed for a hyperbolic paraboloid saddle shell (originally used by Lin and Scordelis) to check the design strength against a consistent design load, therefore, to verify the adequacy of design practice for reinforced concrete shells. Based on nonlinear analyses performed, the analytically calculated ultimate load exceeded the design ultimate load from 14-43% for an analysis with relatively low to high tension stiffening, ${\gamma}$ =5~20 cases. For these cases, the design method gives a lower bound on the ultimate load with respect to Lower bound theorem. This shows the adequacy of the current practice at least for this saddle shell case studied. To generalize the conclusion many more designs-analyses are performed with different shell configurations.