• Title/Summary/Keyword: preconditioning

Search Result 261, Processing Time 0.022 seconds

PRECONDITIONING FOR THE p-VERSION BOUNDARY ELEMENT METHOD IN THREE DIMENSIONS WITH TRIANGULAR ELEMENTS

  • Cao, Wei-Ming;Guo, Benqi
    • Journal of the Korean Mathematical Society
    • /
    • v.41 no.2
    • /
    • pp.345-368
    • /
    • 2004
  • A preconditioning algorithm is developed in this paper for the iterative solution of the linear system of equations resulting from the p-version boundary element approximation of the three dimensional integral equation with hypersingular operators. The preconditioner is derived by first making the nodal and side basis functions locally orthogonal to the element internal bases, and then by decoupling the nodal and side bases from the internal bases. Its implementation consists of solving a global problem on the wire-basket and a series of local problems defined on a single element. Moreover, the condition number of the preconditioned system is shown to be of order $O((1+ln/p)^{7})$. This technique can be applied to discretization with triangular elements and with general basis functions.

ANALYSIS OF FLOW FIELD AROUND NON-LIFTING FORWARD FLIGHT ROTOR USING LOW MACH NUMBER PRECONDITIONING (저마하수 예조건화 기법을 이용한 무양력 전진 비행 로터 주위 유동장 해석)

  • Kim, Jee-Woong;Park, Soo-Hyung
    • 한국전산유체공학회:학술대회논문집
    • /
    • 2008.03a
    • /
    • pp.251-255
    • /
    • 2008
  • Flow field around helicopter involves incompressible flow near the blade root and compressible flow at the blade tip. A problem occurs for low Mach number flow due to the stiffness of the governing equations. Time-derivative preconditioning techniques have been incorporated to reduce the stiffness that occurs at low speed region. The preconditioned form of the compressible Navier-Stokes and Euler equations is used. Computations are performed for the Caradonna-Tung's hovering and non-lifting forward flight case. Computational results are in good agreement with the experimental data.

  • PDF

ANALYSIS OF FLOW FIELD AROUND NON-LIFTING FORWARD FLIGHT ROTOR USING LOW MACH NUMBER PRECONDITIONING (저마하수 예조건화 기법을 이용한 무양력 전진 비행 로터 주위 유동장 해석)

  • Kim, Jee-Woong;Park, Soo-Hyung
    • 한국전산유체공학회:학술대회논문집
    • /
    • 2008.10a
    • /
    • pp.251-255
    • /
    • 2008
  • Flow field around helicopter involves incompressible flow near the blade root and compressible flow at the blade tip. A problem occurs for low Mach number flow due to the stiffness of the governing equations. Time-derivative preconditioning techniques have been incorporated to reduce the stiffness that occurs at low speed region. The preconditioned form of the compressible Navier-Stokes and Euler equations is used. Computations are performed for the Caradonna-Tung's hovering and non-lifting forward flight case. Computational results are in good agreement with the experimental data.

  • PDF

An Approach to Alleviate Cancellation Problem of Temperature Preconditioned Navier-Stokes Equations (온도예조건화 나비어스톡스 방정식의 계산오차 문제 완화 방법 연구)

  • Lee, Sang-Hyeon
    • Journal of the Korean Society of Propulsion Engineers
    • /
    • v.14 no.1
    • /
    • pp.11-19
    • /
    • 2010
  • An approach to reduce cancellation errors of the temperature preconditioned Navier-Stokes equations is proposed. This approach is also applied to the conventional preconditioning methods. Adiabatic laminar viscous flows around a circular cylinder are calculated at different Mach numbers. It is shown that a redefinition of total enthalpy for reducing magnitude of the enthalpy remarkably alleviates cancellation problems of the temperature preconditioning.

HYBRID REORDERING STRATEGIES FOR ILU PRECONDITIONING OF INDEFINITE SPARSE MATRICES

  • Lee Eun-Joo;Zgang Jun
    • Journal of applied mathematics & informatics
    • /
    • v.22 no.1_2
    • /
    • pp.307-316
    • /
    • 2006
  • Incomplete LU factorization preconditioning techniques often have difficulty on indefinite sparse matrices. We present hybrid reordering strategies to deal with such matrices, which include new diagonal reorderings that are in conjunction with a symmetric nondecreasing degree algorithm. We first use the diagonal reorderings to efficiently search for entries of single element rows and columns and/or the maximum absolute value to be placed on the diagonal for computing a nonsymmetric permutation. To augment the effectiveness of the diagonal reorderings, a nondecreasing degree algorithm is applied to reduce the amount of fill-in during the ILU factorization. With the reordered matrices, we achieve a noticeable improvement in enhancing the stability of incomplete LU factorizations. Consequently, we reduce the convergence cost of the preconditioned Krylov subspace methods on solving the reordered indefinite matrices.

Study on Preconditioning of the clavier-Stokes Equations Using 3-Dimensional Unstructured Meshes (3차원 비정렬격자계를 이용한 Navier-Stokes해의 Preconditioning에 관한 연구)

  • Nam, Young-Sok;Choi, Hyoung-Gwon;Yoo, Jung-Yul
    • Transactions of the Korean Society of Mechanical Engineers B
    • /
    • v.25 no.11
    • /
    • pp.1581-1593
    • /
    • 2001
  • An efficient variable-reordering method for finite element meshes is used and the effect of variable-reordering is investigated. For the element renumbering of unstructured meshes, Cuthill-McKee ordering is adopted. The newsy reordered global matrix has a much narrower bandwidth than the original one, making the ILU preconditioner perform bolter. The effect of variable reordering on the convergence behaviour of saddle point type matrix it studied, which results from P2/P1 element discretization of the Navier-Stokes equations. We also propose and test 'level(0) preconditioner'and 'level(2) ILU preconditioner', which are another versions of the existing 'level(1) ILU preconditioner', for the global matrix generated by P2/P1 finite element method of incompressible Navier-Stokes equations. We show that 'level(2) ILU preconditioner'performs much better than the others only with a little extra computations.

EIGENVALUES FOR THE SEMI-CIRCULANT PRECONDITIONING OF ELLIPTIC OPERATORS WITH THE VARIABLE COEFFICIENTS

  • Kim, Hoi-Sub;Kim, Sang-Dong;Lee, Yong-Hun
    • Journal of the Korean Mathematical Society
    • /
    • v.44 no.3
    • /
    • pp.627-645
    • /
    • 2007
  • We investigate the eigenvalues of the semi-circulant preconditioned matrix for the finite difference scheme corresponding to the second-order elliptic operator with the variable coefficients given by $L_vu\;:=-{\Delta}u+a(x,\;y)u_x+b(x,\;y)u_y+d(x,\;y)u$, where a and b are continuously differentiable functions and d is a positive bounded function. The semi-circulant preconditioning operator $L_cu$ is constructed by using the leading term of $L_vu$ plus the constant reaction term such that $L_cu\;:=-{\Delta}u+d_cu$. Using the field of values arguments, we show that the eigenvalues of the preconditioned matrix are clustered at some number. Some numerical evidences are also provided.

Study on Robustness of Incomplete Cholesky Factorization using Preconditioning for Conjugate Gradient Method (불완전분해법을 전처리로 하는 공액구배법의 안정화에 대한 연구)

  • Ko, Jin-Hwan;Lee, Byung-Chai
    • Transactions of the Korean Society of Mechanical Engineers A
    • /
    • v.27 no.2
    • /
    • pp.276-284
    • /
    • 2003
  • The preconditioned conjugate gradient method is an efficient iterative solution scheme for large size finite element problems. As preconditioning method, we choose an incomplete Cholesky factorization which has efficiency and easiness in implementation in this paper. The incomplete Cholesky factorization mettled sometimes leads to breakdown of the computational procedure that means pivots in the matrix become minus during factorization. So, it is inevitable that a reduction process fur stabilizing and this process will guarantee robustness of the algorithm at the cost of a little computation. Recently incomplete factorization that enhances robustness through increasing diagonal dominancy instead of reduction process has been developed. This method has better efficiency for the problem that has rotational degree of freedom but is sensitive to parameters and the breakdown can be occurred occasionally. Therefore, this paper presents new method that guarantees robustness for this method. Numerical experiment shows that the present method guarantees robustness without further efficiency loss.

Numerical Characteristics of Upwind Schemes for Preconditioned Compressible Navier-Stokes Equations (예조건화된 압축성유동 수치기법에서의 풍상차분법의 수치특성 검토)

  • Gill J. H.;Lee D. H.;Choi Y. H.;Kwon J. H.;Lee S. S.
    • 한국전산유체공학회:학술대회논문집
    • /
    • 2002.10a
    • /
    • pp.95-102
    • /
    • 2002
  • Studies of the numerical characteristics of implicit upwind schemes, such as upwind ADI, Line Gauss-Seidel(LGS) and Point Gauss-Seidel(LU) algorithms, for preconditioned Navier-Stokes equations ate performed. All the algorithms are expressed in approximate factorization form and Von Neumann stability analysis and convergence studies are made. Preconditioning is applied for efficient convergence at low Mach numbers and low Reynolds numbers. For high aspect ratio computations, the ADI and LGS algorithms show efficient and uniform convergence up to moderate aspect ratio if we adopt viscous preconditioning based on min- CFL/max- VNN time-step definition. The LU algorithm, on the other hand, shows serious deterioration in convergence rate as the grid aspect ratio increases. Computations for practical applications also verify these results.

  • PDF

Acupoint Electrical Stimulation Preconditioning to Reduce Hemorrhoid Pain: A Case Report (경혈 전기 자극 전처치요법을 통해 치료한 치핵 통증: 증례보고)

  • Lee, Seungmin Kathy
    • Korean Journal of Acupuncture
    • /
    • v.37 no.3
    • /
    • pp.198-202
    • /
    • 2020
  • A 30-year-old woman with Grade III hemorrhoid complained of excruciating pain that continued for several hours, especially with defecation. She was not able to frequent the clinic due to COVID-19 shutdowns, therefore additional treatment using acupoint electrical stimulation (AES) was self-administered. She administered AES bilaterally on BL57 and LR10 for fifteen minutes before each defecation as a preconditioning treatment. She assessed her pain using a Numerical Rating Scale (NRS) during defecation, and 3 hours later. The patient initially complained of pain rating 9 on the NRS. After the first session of AES, the pain dropped to 5. On one-month follow-up, the pain was at 3 and the patient was able to terminate all treatment. Self-administered AES preconditioning at BL57 and LR10 can be used to reduce extreme cases of hemorrhoid pain.