• 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.

• Bulletin of the Korean Mathematical Society
• /
• v.46 no.4
• /
• pp.627-644
• /
• 2009

The gauge-Uzawa method which has been constructed in [11] is a projection type method to solve the evolution Navier-Stokes equations. The method overcomes many shortcomings of projection methods and displays superior numerical performance [11, 12, 15, 16]. However, we have obtained only suboptimal accuracy via the energy estimate in [11]. In this paper, we study semi-discrete gauge-Uzawa method to prove optimal accuracy via energy estimate. The main key in this proof is to construct the intermediate equation which is formed to gauge-Uzawa algorithm. We will estimate velocity errors via comparing with the intermediate equation and then evaluate pressure errors via subtracting gauge-Uzawa algorithm from Navier-Stokes equations.

• Bulletin of the Korean Mathematical Society
• /
• v.52 no.5
• /
• pp.1669-1681
• /
• 2015

In this paper, we provide semi-convergence results of the parameterized inexact Uzawa method with singular preconditioners for solving singular saddle point problems. We also provide numerical experiments to examine the effectiveness of the parameterized inexact Uzawa method with singular preconditioners.

• Korean Journal of Mathematics
• /
• v.22 no.4
• /
• pp.645-658
• /
• 2014

Many projection methods have been progressively constructed to find more accurate and efficient solution of the Navier-Stokes equations. In this paper, we consider most recently constructed projection methods: the pressure correction method, the gauge method, the consistent splitting method, the Gauge-Uzawa method, and the stabilized Gauge-Uzawa method. Each method has different background and theoretical proof. We prove equivalentness of the pressure correction method and the stabilized Gauge-Uzawa method. Also we will obtain that the Gauge-Uzawa method is equivalent to the gauge method and the consistent splitting method. We gather theoretical results of them and conclude that the results are also valid on other equivalent methods.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.15 no.3
• /
• pp.233-251
• /
• 2011

The Gauge-Uzawa method [GUM] in [9] which is a projection type algorithm to solve evolution Navier-Stokes equations has many advantages and superior performance. But this method has been studied for backward Euler time discrete scheme which is the first order technique, because the classical second order GUM requests rather strong stability condition. Recently, the second order time discrete GUM was modified to be unconditionally stable and estimated errors in [12]. In this paper, we contemplate several GUMs which can be derived by the same manner within [12], and we dig out properties of them for both stability and accuracy. In addition, we evaluate an stability condition for the classical GUM to construct an adaptive GUM for time to make free from strong stability condition of the classical GUM.

• Korean Journal of Mathematics
• /
• v.21 no.2
• /
• pp.125-150
• /
• 2013

The stabilized Gauge-Uzawa method (SGUM), which is a second order projection type algorithm to solve the time-dependent Navier-Stokes equations, has been newly constructed in 2013 Pyo's paper. The accuracy of SGUM has been proved only for time discrete scheme in the same paper, but it is crucial to study for fully discrete scheme, because the numerical errors depend on discretizations for both space and time, and because discrete spaces between velocity and pressure can not be chosen arbitrary. In this paper, we find out properties of the fully discrete SGUM and estimate its errors and stability to solve the evolution Navier-Stokes equations. The main difficulty in this estimation arises from losing some cancellation laws due to failing divergence free condition of the discrete velocity function. This result will be extended to Boussinesq equations in the continuous research (part II) and is essential in the study of part II.

• Korean Journal of Mathematics
• /
• v.27 no.1
• /
• pp.193-219
• /
• 2019

The Navier-Stokes equations with variable density are challenging problems in numerical analysis community. We recently built the 2nd order stabilized Gauge-Uzawa method [SGUM] to solve the Navier-Stokes equations with constant density and have estimated theoretically optimal accuracy. Also we proved that SGUM is unconditionally stable. In this paper, we apply SGUM to the Navier-Stokes equations with nonconstant variable density and find out the stability condition of the algorithms. Because the condition is rather strong to apply to real problems, we consider Allen-Cahn scheme to construct unconditionally stable scheme.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.24 no.2
• /
• pp.121-141
• /
• 2020

Nanofluids is the fluids mixed with nanoscale particles and the mixed nano size materials affect heat transport. Researchers in this field has been focused on modeling and numerical computation by engineers In this paper, we analyze stability constraint of the dominant equations and check validate of the condition for most kinds of materials. So we mathematically analyze stability of the system. Also we apply Gauge-Uzawa algorithm to solve the system and prove stability of the 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.

• Bulletin of the Korean Mathematical Society
• /
• v.61 no.1
• /
• pp.93-115
• /
• 2024

In this paper, the Gauge-Uzawa methods for the Darcy-Brinkman equations driven by temperature and salt concentration (DBTC) are proposed. The first order backward difference formula is adopted to approximate the time derivative term, and the linear term is treated implicitly, the nonlinear terms are treated semi-implicit. In each time step, the coupling elliptic problems of velocity, temperature and salt concentration are solved, and then the pressure is solved. The unconditional stability and error estimations of the first order semi-discrete scheme are derived, at the same time, the unconditional stability of the first order fully discrete scheme is obtained. Some numerical experiments verify the theoretical prediction and show the effectiveness of the proposed methods.