• Journal of the Korean Mathematical Society
• /
• v.34 no.1
• /
• pp.87-118
• /
• 1997

This paper is concerned with a mixed Lagrangian formulation of the wiscous, stationary, incompressible Navier-Stokes equations $$ (1.1) -\nu\Delta u + (u \cdot \nabla)u + \nabla_p = f in \Omega $$ and $$ (1.2) \nubla \cdot u = 0 in \Omega $$ along with inhomogeneous Dirichlet boundary conditions on a portion of the boundary $$ (1.3) u = ^{0 on \Gamma_0 _{g on \Gamma_g, $$ where $\Omega$ is a bounded open domain in $R^d, d = 2 or 3$, or with a boundary $\Gamma = \partial\Omega$, which is composed of two disjoint parts $\Gamma_0$ and $\Gamma_g$.

• Coupled systems mechanics
• /
• v.7 no.5
• /
• pp.583-610
• /
• 2018

This paper discusses the issues associated with modeling frictional contact between solid bodies undergoing large deformations. The most common model for friction on contact interfaces in solid mechanics is the Coulomb friction model, in which two distinct responses are possible: stick and slip. Handling the transition between these two phases computationally has been a source of algorithmic instability, lack of convergence and non-unique solutions, particularly in the presence of large deformations. Most computational models for frictional contact have used penalty or updated Lagrangian approaches to enforce frictional contact conditions. These two approaches, however, present some computational challenges due to conditioning issues in penalty-type implementations and the iterative nature of the updated Lagrangian formulation, which, particularly in large simulations, may lead to relatively slow convergence. Alternatively, a plasticity-inspired implementation of frictional contact has been shown to handle the stick-slip conditions in a local, algorithmically efficient manner that substantially reduces computational cost and successfully avoids the issues of instability and lack of convergence often reported with other methods (Laursen and Simo 1993). The formulation of this approach, however, has been limited to the small deformations realm, a fact that severely limited its application to contact problems where large deformations are expected. In this paper, we present an algorithmically consistent formulation of this method that preserves its key advantages, while extending its application to the realm of large-deformation contact problems. We show that the method produces results similar to the augmented Lagrangian formulation at a reduced computational cost.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.2
• /
• pp.177-189
• /
• 2015

This paper presents an all-direct domain decomposition approach for large-scale structural analysis. The proposed approach achieves computational robustness and efficiency by enforcing the compatibility of the displacement field across the sub-domain boundaries via local Lagrange multipliers and augmented Lagrangian formulation (ALF). The proposed domain decomposition approach was compared to the existing FETI approach in terms of the computational time and memory usage. The parallel implementation of the proposed algorithm was described in detail. Finally, a preliminary validation was attempted for the proposed approach, and the numerical results of two- and three-dimensional problems were compared to those obtained through a dual-primal FETI approach. The results indicate an improvement in the performance as a result of the implementing the proposed approach.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.12 no.3
• /
• pp.272-281
• /
• 2009

In this paper, three-dimensional optimal evasive maneuver patterns for air-to-surface attack missiles against proportionally navigated anti-air defense missiles were investigated. An interception error of the defense missile is produced by an evasive maneuver of the attack missile. It is assumed that the defense missiles are continuously launched during the flight of attack missile. The performance index to be minimized is then defined as the negative square integral of the interception errors. The direct parameter optimization technique based on SQP and a co-evolution method based on the augmented Lagrangian formulation are adopted to get the attack missile's optimal evasive maneuver patterns. The overall shape of the resultant optimal evasive maneuver is represented as a deformed barrel-roll.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.6
• /
• pp.469-476
• /
• 2012

In this paper, a finite element domain decomposition method using local and mixed Lagrange multipliers for a large scal structural analysis is presented. The proposed algorithms use local and mixed Lagrange multipliers to improve computational efficiency. In the original FETI method, classical Lagrange multiplier technique was used. In the dual-primal FETI method, the interface nodes are used at the corner nodes of each sub-domain. On the other hand, the proposed FETI-local analysis adopts localized Lagrange multipliers and the proposed FETI-mixed analysis uses both global and local Lagrange multipliers. The numerical analysis results by the proposed algorithms are compared with those obtained by dual-primal FETI method.

• Journal of Magnetics
• /
• v.16 no.1
• /
• pp.77-82
• /
• 2011

This paper presents a new self-adjoint material sensitivity formulation for optimal designs and inverse problems in the high frequency domain. The proposed method is based on the continuum approach using the augmented Lagrangian method. Using the self-adjoint formulation, there is no need to solve the adjoint system additionally when the goal function is a function of the S-parameter. In addition, the algorithm is more general than most previous approaches because it is independent of specific analysis methods or gridding techniques, thereby enabling the use of commercial EM simulators and various custom solvers. For verification, the method was applied to the several numerical examples of dielectric material reconstruction problems in the high frequency domain, and the results were compared with those calculated using the conventional method.

• Computers and Concrete
• /
• v.1 no.2
• /
• pp.131-150
• /
• 2004

We present in this work a coupled phenomenological chemo-mechanical model that represents the degradation of concrete-like materials. The chemical behaviour is described by the nowadays well known simplified calcium leaching approach. And the mechanical damage behaviour is described by a continuum damage model which involves the gradient of the damage quantity. The coupled nonlinear problem at hand is addressed within the context of the finite element method. For the equation governing the calcium dissolution-diffusion part of the problem, special care is taken to treat the highly nonlinear calcium conductivity and solid calcium functions. The algorithmic design is based on a Newton-type iterative scheme where use is made of a recently proposed relaxed linearization procedure. And for the equation governing the damage part of the problem, an augmented Lagrangian formulation is used to take into account the damage irreversibility constraint. Finally, numerical simulations are compared with experimental results on cement paste.

• Proceedings of the KIEE Conference
• /
• 2007.11b
• /
• pp.226-228
• /
• 2007

This paper proposes Nonlinear Interior Point Method (NIPM) including discrete control variables optimal power flow formulations. The algorithm utilizes the robustness in terms of starting point and fast convergence for large scale power system of NIPM and an introduction of rounding penalty function which is augmented in the Lagrangian function to handle discrete control variables. The derived formulation shows a simplified approach to deal with discrete control problems which is implementable in real large scale systems.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.18 no.1
• /
• pp.17-26
• /
• 2014

This paper presents development of an improved domain decomposition method for large scale structural problem that aims to provide high computational efficiency. In the previous researches, we developed the domain decomposition algorithm based on augmented Lagrangian formulation and proved numerical efficiency under both serial and parallel computing environment. In this paper, new computational analysis by the proposed domain decomposition method is performed. For this purpose, reduction in computational time achieved by the proposed algorithm is compared with that obtained by the dual-primal FETI method under serial computing condition. It is found that the proposed methods significantly accelerate the computational speed for a linear structural problem.

• Journal of Mechanical Science and Technology
• /
• v.19 no.spc1
• /
• pp.388-394
• /
• 2005

Recently, the authors have developed a method for real-time dynamics of multibody systems, which combines a semi-recursive formulation to derive the equations of motion in dependent relative coordinates, along with an augmented Lagrangian technique to impose the loop closure conditions. The following numerical integration procedures, which can be grouped into the so-called structural integrators, were tested : trapezoidal rule, Newmark dissipative schemes, HHT rule, and the Generalized-${\alpha}$ family. It was shown that, for large multi body systems, Newmark dissipative was the best election since, provided that the adequate parameters were chosen, excellent behavior was achieved in terms of efficiency and robustness with acceptable levels of accuracy. In the present paper, the performance of the described method in combination with another group of integrators, the Implicit Runge-Kutta family (IRK), is analyzed. The purpose is to clarify which kind of IRK algorithms can be more suitable for real-time applications, and to see whether they can be competitive with the already tested structural family of integrators. The final objective of the work is to provide some practical criteria for those interested in achieving real-time performance for large and complex multibody systems.