• 제목/요약/키워드: Newton linearization

검색결과 16건 처리시간 0.024초

Newton-Krylov Method for Compressible Euler Equations on Unstructured Grids

  • Kim Sungho;Kwon Jang Hyuk
    • 한국전산유체공학회:학술대회논문집
    • /
    • 한국전산유체공학회 1998년도 추계 학술대회논문집
    • /
    • pp.153-159
    • /
    • 1998
  • The Newton-Krylov method on the unstructured grid flow solver using the cell-centered spatial discretization oi compressible Euler equations is presented. This flow solver uses the reconstructed primitive variables to get the higher order solutions. To get the quadratic convergence of Newton method with this solver, the careful linearization of face flux is performed with the reconstructed flow variables. The GMRES method is used to solve large sparse matrix and to improve the performance ILU preconditioner is adopted and vectorized with level scheduling algorithm. To get the quadratic convergence with the higher order schemes and to reduce the memory storage. the matrix-free implementation and Barth's matrix-vector method are implemented and compared with the traditional matrix-vector method. The convergence and computing times are compared with each other.

  • PDF

자기동조 가중최소자승법을 이용한 AOA 측위 알고리즘 개발 (Development of an AOA Location Method Using Self-tuning Weighted Least Square)

  • 이성호;김동혁;노기홍;박경순;성태경
    • 제어로봇시스템학회논문지
    • /
    • 제13권7호
    • /
    • pp.683-687
    • /
    • 2007
  • In last decades, several linearization methods for the AOA measurements have been proposed, for example, Gauss-Newton method and Closed-Form solution. Gauss-Newton method can achieve high accuracy, but the convergence of the iterative process is not always ensured if the initial guess is not accurate enough. Closed-Form solution provides a non-iterative solution and it is less computational. It does not suffer from convergence problem, but estimation error is somewhat larger. This paper proposes a Self-Tuning Weighted Least Square AOA algorithm that is a modified version of the conventional Closed-Form solution. In order to estimate the error covariance matrix as a weight, a two-step estimation technique is used. Simulation results show that the proposed method has smaller positioning error compared to the existing methods.

ON THE LINEARIZATION OF DEFECT-CORRECTION METHOD FOR THE STEADY NAVIER-STOKES EQUATIONS

  • Shang, Yueqiang;Kim, Do Wan;Jo, Tae-Chang
    • 대한수학회지
    • /
    • 제50권5호
    • /
    • pp.1129-1163
    • /
    • 2013
  • Based on finite element discretization, two linearization approaches to the defect-correction method for the steady incompressible Navier-Stokes equations are discussed and investigated. By applying $m$ times of Newton and Picard iterations to solve an artificial viscosity stabilized nonlinear Navier-Stokes problem, respectively, and then correcting the solution by solving a linear problem, two linearized defect-correction algorithms are proposed and analyzed. Error estimates with respect to the mesh size $h$, the kinematic viscosity ${\nu}$, the stability factor ${\alpha}$ and the number of nonlinear iterations $m$ for the discrete solution are derived for the linearized one-step defect-correction algorithms. Efficient stopping criteria for the nonlinear iterations are derived. The influence of the linearizations on the accuracy of the approximate solutions are also investigated. Finally, numerical experiments on a problem with known analytical solution, the lid-driven cavity flow, and the flow over a backward-facing step are performed to verify the theoretical results and demonstrate the effectiveness of the proposed defect-correction algorithms.

Error Analysis of time-based and angle-based location methods

  • Kim, Dong-Hyouk;Song, Seung-Hun;Sung, Tae-Kyung
    • 한국항해항만학회:학술대회논문집
    • /
    • 한국항해항만학회 2006년도 International Symposium on GPS/GNSS Vol.1
    • /
    • pp.479-483
    • /
    • 2006
  • Indoor positioning is recently highlighted and various kinds of indoor positioning systems are under developments. Since positioning systems have their own characteristics, proper positioning scheme should be chosen according to the required specifications. Positioning methods are often classified into time-based and angle-based one, and this paper presents the error analysis of these location methods. Because measurement equations of these methods are nonlinear, linearization is usually needed to get the position estimate. In this paper, Gauss-Newton method is used in the linearization. To analyze the position error, we investigate the error ellipse parameters that include eccentricity, rotation angle, and the size of ellipse. Simulation results show that the major axes of error ellipses of TOA and AOA method lie in different quadrants at most region of workspace, especially where the geometry is poor. When the TOA/AOA hybrid scheme is employed, it is found that the error ellipse is reduced to the intersection of ellipses of TOA and AOA method.

  • PDF

시간기반과 각도기반의 측위방식 성능비교 및 오차 특성 분석 (Error Analysis of Time-Based and Angle-Based Location Methods)

  • 김동혁;송승헌;성태경
    • 제어로봇시스템학회논문지
    • /
    • 제12권10호
    • /
    • pp.962-967
    • /
    • 2006
  • Indoor positioning is highlighted recently, and various kinds of indoor positioning systems are under developments. Since positioning systems have their own characteristics, proper positioning scheme should be chosen according to the required specifications. Positioning methods are classified into time-based and angle-based one. This paper presents the error analysis of time-based and angle-based location methods. Because measurements of these methods are nonlinear, linearizations are needed in both cases to estimate the user position. In the linearization, Gauss-Newton method is used in both cases. To analyze the position error, we investigate the error ellipse parameters that include eccentricity, rotation angle, and the size of ellipse. Simulation results show that the major axes of TOA and AOA method lie in different quadrants at most region of workspace, especially where the geometry is poor. When the TOA/AOA hybrid is employed, it is found that the error ellipse is reduced to the intersection of ellipses of TOA and AOA.

Development of an AOA Location Method Using Covariance Estimation

  • Lee, Sung-Ho;Roh, Gi-Hong;Sung, Tae-Kyung
    • 한국항해항만학회:학술대회논문집
    • /
    • 한국항해항만학회 2006년도 International Symposium on GPS/GNSS Vol.1
    • /
    • pp.485-489
    • /
    • 2006
  • In last decades, several linearization methods for the AOA measurements have been proposed, for example, Gauss-Newton method and closed-form solution. Gauss-Newton method can achieve high accuracy, but the convergence of the iterative process is not always ensured if the initial guess is not accurate enough. Closed-form solution provides a non-iterative solution and it is less computational. It does not suffer from convergence problem, but estimation error is somewhat larger. This paper proposes a self-tuning weighted least square AOA algorithm that is a modified version of the conventional closed-form solution. In order to estimate the error covariance matrix as a weight, two-step estimation technique is used. Simulation results show that the proposed method has smaller positioning error compared to the existing methods.

  • PDF

고차 판 유한요소의 기하학적 비선형 해석 (Geometrically Nonlinear Analysis of Higher Order Plate Bending Finite Element)

  • 신영식
    • 대한토목학회논문집
    • /
    • 제8권3호
    • /
    • pp.1-10
    • /
    • 1988
  • 본 연구에서는 고차 판 유한요소의 판의 기하학적 비선형 해석에의 적용성을 고찰한다. 고차판요소는 3 차원 연속체로부터 Total Lagrangian 형태로 나타낸 운동방정식을 이산화하고 고차 판이론을 도입하여 유도한다. 유한변형을 고려한 기하학적 비션형 방정식은 Newton-Raphson반복법으로 내력벡터를 선형화하여 강도매트릭스를 반복계산하여 푼다. 요소매트릭스는 shear locking 현상을 피하기 위하여 Gauss 적분법을 이용한 선택적 감차적분으로 계산한다. 여러가지 예제해석을 통하여 고차 판요소의 효율성과 정확도를 고찰하였다.

  • PDF

Analytical approximate solution for Initial post-buckling behavior of pipes in oil and gas wells

  • Yu, Yongping;Sun, Youhong;Han, Yucen
    • Coupled systems mechanics
    • /
    • 제1권2호
    • /
    • pp.155-163
    • /
    • 2012
  • This paper presents analytical approximate solutions for the initial post-buckling deformation of the pipes in oil and gas wells. The governing differential equation with sinusoidal nonlinearity can be reduced to form a third-order-polynomial nonlinear equation, by coupling of the well-known Maclaurin series expansion and orthogonal Chebyshev polynomials. Analytical approximations to the resulting boundary condition problem are established by combining the Newton's method with the method of harmonic balance. The linearization is performed prior to proceeding with harmonic balancing thus resulting in a set of linear algebraic equations instead of one of non-linear algebraic equations, unlike the classical method of harmonic balance. We are hence able to establish analytical approximate solutions. The approximate formulae for load along axis, and periodic solution are established for derivative of the helix angle at the end of the pipe. Illustrative examples are selected and compared to "reference" solution obtained by the shooting method to substantiate the accuracy and correctness of the approximate analytical approach.

Two-scale approaches for fracture in fluid-saturated porous media

  • de Borst, Rene;Rethore, Julien;Abellan, Marie-Angele
    • Interaction and multiscale mechanics
    • /
    • 제1권1호
    • /
    • pp.83-101
    • /
    • 2008
  • A derivation is given of two-scale models that are able to describe deformation and flow in a fluid-saturated and progressively fracturing porous medium. From the micromechanics of the flow in the cavity, identities are derived that couple the local momentum and the mass balances to the governing equations for a fluid-saturated porous medium, which are assumed to hold on the macroscopic scale. By exploiting the partition-of-unity property of the finite element shape functions, the position and direction of the fractures are independent from the underlying discretization. The finite element equations are derived for this two-scale approach and integrated over time. The resulting discrete equations are nonlinear due to the cohesive crack model and the nonlinearity of the coupling terms. A consistent linearization is given for use within a Newton-Raphson iterative procedure. Finally, examples are given to show the versatility and the efficiency of the approach.