• Title/Summary/Keyword: Predictor-Corrector scheme

Search Result 34, Processing Time 0.036 seconds

A Grid Generation Technique for the External Flow Fields Utilizing the Predictor-Corrector Scheme (Predictor-Corrector를 활용한 외부 유동장 격자 생성 기법)

  • Kim B. S.
    • Journal of computational fluids engineering
    • /
    • v.2 no.1
    • /
    • pp.84-92
    • /
    • 1997
  • In this paper a new structured grid generation technique is introduced. This new technique utilizes predictor-corrector approach, and is a marching scheme in the global sense as the hyperbolic scheme is. In the predictor step, one layer of grid cells is obtained by using Modified Advancing Front Method which generates a collection of quadrilateral cells simultaneously. In the corrector step, the layer of grid cells that is calculated in the predictor step is adjusted by solving Laplace equations to prevent grid lines from skewing and overlapping in highly curved configurations. It is shown that the resultant algorithm, named a MAP scheme, which combines the Modified Advancing Front Method as a Predictor with an elliptic scheme as a corrector can be used to generate globally smooth and locally near-orthogonal grids for external flow fields even for highly curved configurations. Examples of grid generations for external flow fields about several configurations by use of the present approach are given, and its applicability and flexibility have been demonstrated and discussed.

  • PDF

A PREDICTOR-CORRECTOR SCHEME FOR THE NUMERICAL SOLUTION OF THE BOUSSINESQ EQUATION

  • Ismail, M.S.;Bratsos, A.G.
    • Journal of applied mathematics & informatics
    • /
    • v.13 no.1_2
    • /
    • pp.11-27
    • /
    • 2003
  • A fourth order in time and second order in space scheme using a finite-difference method is developed for the non-linear Boussinesq equation. For the solution of the resulting non-linear system a predictor-corrector pair is proposed. The method is analyzed for local truncation error and stability. The results of a number of numerical experiments for both the single and the double-soliton waves are given.

Numerical Simulation of Mold Filling Processes of Castings by using of Predictor-two step Corrector-VOF

  • Xun, Sun;Junqing, Wang;Hwang, Ho-Young;Choi, Jeong-Kil
    • Journal of Korea Foundry Society
    • /
    • v.22 no.6
    • /
    • pp.299-303
    • /
    • 2002
  • 수정된 SIMPLE법과 VOF의 결합으로 predictor-two step corrector-VOF라고 불리는 새로운 알고리즘이 주조 시 용탕 충전과정을 해석하기 위해 개발되었다. 운동량보존으로부터 유도된 새 2단계 속도 경계조건 처리법은 용탕의 자유표면을 추적하는 데 사용되었다. 본 연구에서는 2개의 예제 계산을 통해 계산정확도와 속도에 대한 Courant 수의 영향을 살펴보았다. 그 결과 적당한 Courant 수의 증가는 계산 정확도의 감소 없이 용탕 계산 속도를 향상시킬 수 있는 것으로 나타났다. 또한 만족할 만한 계산 정확도와 효율이 이 알고리즘의 실제 제품 해석을 통해 얻어졌다.

A spent nuclear fuel source term calculation code BESNA with a new modified predictor-corrector scheme

  • Duy Long Ta ;Ser Gi Hong ;Dae Sik Yook
    • Nuclear Engineering and Technology
    • /
    • v.54 no.12
    • /
    • pp.4722-4730
    • /
    • 2022
  • This paper introduces a new point depletion-based source term calculation code named BESNA (Bateman Equation Solver for Nuclear Applications), which is aimed to estimate nuclide inventories and source terms from spent nuclear fuels. The BESNA code employs a new modified CE/CM (Constant Extrapolation - Constant Midpoint) predictor-corrector scheme in depletion calculations for improving computational efficiency. In this modified CE/CM scheme, the decay components leading to the large norm of the depletion matrix are excluded in the corrector, and hence the corrector calculation involves only the reaction components, which can be efficiently solved with the Talyor Expansion Method (TEM). The numerical test shows that the new scheme substantially reduces computing time without loss of accuracy in comparison with the conventional scheme using CRAM (Chebyshev Rational Approximation Method), especially when the substep calculations are applied. The depletion calculation and source term estimation capability of BESNA are verified and validated through several problems, where results from BESNA are compared with those calculated by other codes as well as measured data. The analysis results show the computational efficiency of the new modified scheme and the reliability of BESNA in both isotopic predictions and source term estimations.

Development of Finite Element Method for the Extended Boussinesq Equations (확장형 Boussinesq 방정식의 유한요소모형 개발)

  • Woo, Seung-Buhm;Choi, Young-Kwang;Yoon, Byung-Il
    • The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
    • /
    • v.12 no.3
    • /
    • pp.133-141
    • /
    • 2007
  • A finite element model is developed for the extended Boussinesq equations that is capable of simulating the dynamics of long and short waves. Galerkin weighted residual method and the introduction of auxiliary variables for 3rd spatial derivative terms in the governing equations are used for the model development. The Adams-Bashforth-Moulton Predictor Corrector scheme is used as a time integration scheme for the extended Boussinesq finite element model so that the truncation error would not produce any non-physical dispersion or dissipation. This developed model is applied to the problems of solitary wave propagation. Predicted results is compared to available analytical solutions and laboratory measurements. A good agreement is observed.

The Efficient Algorithm for Simulating the Multiphase Flow

  • Yoon Seong Y;Yabe T.
    • Journal of computational fluids engineering
    • /
    • v.9 no.1
    • /
    • pp.18-24
    • /
    • 2004
  • The unified simulation for the multiphase flow by predictor-corrector scheme based on CIP method is introduced. In this algorithm, the interface between different phases is identified by a density function and tracked by solving an advection equation. Solid body motion is modeled by the translation and angular motion. The mathematical formulation and numerical results are also described. To verify the efficiency, accuracy and capability of proposed algorithm, two dimensional incompressible cavity flow, the motion of a floating ball into water and a single rising bubble by buoyancy force are numerically simulated by the present scheme. As results, it is confirmed that the present scheme gives an efficient, stable and reasonable solution in the multiphase flow problem.

Analysis of the Flow in Square Cavity Using CIP Method in Non-staggered Grid Arrangement (비엇갈림 격자계에서 CIP법을 이용한 캐비티내의 유동해석)

  • Lee, Jung-Hee;Kang, Joon;Im, Do-Kiun;Kim, Charn-Jung
    • Proceedings of the KSME Conference
    • /
    • 2003.04a
    • /
    • pp.1821-1826
    • /
    • 2003
  • In this study, we discuss CIP method, which can treat compressible/incompressible problem and multi-phase problem. We can apply this method to the general equations by using CCUP. In this paper, non-staggered grid arrangement and predictor-corrector method are used to enhance later development and the solution accuracy and convergence performance. To validate the numerical algorithm proposed in this paper, the two-dimensional unsteady flow in the driven cavity is simulated. The numerical results of this subject using the present algorithm are compared with other numerical results. As a result, it is prived that the present scheme gives stable and improved solutions, and the results even coarse grid are in good agreement with other result using a fine grid arrangement.

  • PDF

Application of TVD-McCormack Scheme to Analysis of Dam-Break Problems (댐붕괴 문제의 해석에 관한 TVD-McCormack기법의 적용)

  • Lee, Jong-Kyu;Kim, Tae-Kwan
    • Journal of Korea Water Resources Association
    • /
    • v.36 no.3 s.134
    • /
    • pp.365-374
    • /
    • 2003
  • This is a study on application of a TVD-Mccormack scheme for the computation of one-dimensional dam-break flows. The TVD scheme not only has the ability to damp out oscillations, but also does not contain terms with adjustable parameters. Moreover, the TVD-McCormack scheme does not cause any additional difficulty when dealing with the source term of the equation and retains second-order accuracy in both space and time. In this study, by appropriately designing the limiter functions, the TVD property can be achieved, and numerical oscillations near a jump discontinuities can be eliminated or reduced. Also, this numerical scheme has less computational errors when the direction of the predictor-corrector step is in the same direction as the shock wane propagation.

Nonlinear Analysis of Underwater Towed Cable Using Robust Nodal Position Finite Element Method (강건 절점위치 유한요소법을 이용한 수중 예인 케이블의 비선형 거동해석)

  • Lee, Euntaek;Go, Gwangsoo;Ahn, Hyung Taek;Kim, Seongil;Chun, Seung Yong;Kim, Jung Suk;Lee, Byeong Hee
    • Journal of the Society of Naval Architects of Korea
    • /
    • v.53 no.5
    • /
    • pp.388-399
    • /
    • 2016
  • A motion analysis of an underwater towed cable is a complex task due to its nonlinear nature of the problem. The major source of the nonlinearity of the underwater cable analysis is that the motion of the cable involves large rigid-body motion. This large rigid-body motion makes difficult to use standard displacement-based finite element method. In this paper, the authors apply recently developed nodal position-based finite element method which can deal with the geometric nonlinearity due to the large rigid-body motion. In order to enhance the stability of the large-scale nonlinear cable motion simulation, an efficient time-integration scheme is proposed, namely predictor/multi-corrector Newmark scheme. Three different predictors are introduced, and the best predictor in terms of stability and robustness for impulsive cable motion analysis is proposed. As a result, the nonlinear motion of underwater cable is predicted in a very efficient manner compared to the classical finite element of finite difference methods. The efficacy of the method is demonstrated with several test cases, involving static and dynamic motion of a single cable element, and also under water towed cable composed of multiple cable elements.

Strongly coupled partitioned six degree-of-freedom rigid body motion solver with Aitken's dynamic under-relaxation

  • Chow, Jeng Hei;Ng, E.Y.K.
    • International Journal of Naval Architecture and Ocean Engineering
    • /
    • v.8 no.4
    • /
    • pp.320-329
    • /
    • 2016
  • An implicit method of solving the six degree-of-freedom rigid body motion equations based on the second order Adams-Bashforth-Moulten method was utilised as an improvement over the leapfrog scheme by making modifications to the rigid body motion solver libraries directly. The implementation will depend on predictor-corrector steps still residing within the hybrid Pressure Implicit with Splitting of Operators - Semi-Implicit Method for Pressure Linked Equations (PIMPLE) outer corrector loops to ensure strong coupling between fluid and motion. Aitken's under-relaxation is also introduced in this study to optimise the convergence rate and stability of the coupled solver. The resulting coupled solver ran on a free floating object tutorial test case when converged matches the original solver. It further allows a varying 70%-80% reduction in simulation times compared using a fixed under-relaxation to achieve the required stability.