• Nuclear Engineering and Technology
• /
• 제52권2호
• /
• pp.258-270
• /
• 2020

A new task of using the Jacobian-Free-Newton-Krylov (JFNK) method for the PWR core transient simulations involving neutronics, thermal hydraulics and mechanics is conducted. For the transient scenario of PWR, normally the Picard iteration of the coupled coarse-mesh nodal equations and parallel channel TH equations is performed to get the transient solution. In order to solve the coupled equations faster and more stable, the Newton Krylov (NK) method based on the explicit matrix was studied. However, the NK method is hard to be extended to the cases with more physics phenomenon coupled, thus the JFNK based iteration scheme is developed for the nodal method and parallel-channel TH method. The local gap conductance is sensitive to the gap width and will influence the temperature distribution in the fuel rod significantly. To further consider the local gap conductance during the transient scenario, a 1D mechanics model is coupled into the JFNK scheme to account for the fuel thermal expansion effect. To improve the efficiency, the physics-based precondition and scaling technique are developed for the JFNK iteration. Numerical tests show good convergence behavior of the iterations and demonstrate the influence of the fuel thermal expansion effect during the rod ejection problems.

• 해양환경안전학회지
• /
• 제27권1호
• /
• pp.193-200
• /
• 2021

본 연구에서는 국부 비선형 정적 해석을 효율적으로 수행하기 위하여 강성응축(Static condensation)을 활용한 해석기법을 제시하였다. 강성응축기법은 자유도 기반의 유한요소 모델 축소기법이며, 해석 모델을 관심 대상(Target) 부분과 응축되어 생략될(Omitted) 부분으로 구분한다. 본 연구에서는, 관심 대상 부분에는 비선형 영역, 생략될 부분에는 선형 영역으로 지정하였고, 선형 영역에 대응되는 강성 행렬 및 하중 벡터를 비선형 영역, 즉 관심 대상 부분으로 모두 응축하였다. 모델 응축 후에는 비선형 영역에 대한 강성 행렬 및 하중 벡터만으로 이루어진 축소 모델을 구성하였으며, 이 축소 모델만을 뉴턴-랩슨 반복(Newton-Raphson iteration)을 통해 갱신하여 효율적으로 비선형 해석을 수행하였다. 끝으로, 제안된 기법을 다양한 수치 예제에 적용하여 해석기법의 효율성과 신뢰성을 제시하였다.

• 호남수학학술지
• /
• 제35권3호
• /
• pp.417-433
• /
• 2013

We consider the iterative solution of a quadratic matrix equation with special coefficient matrices which arises in the quasibirth and death problem. In this paper, we show that the elementwise minimal positive solvent of the quadratic matrix equations can be obtained using Newton's method if there exists a positive solvent and the convergence rate of the Newton iteration is quadratic if the Fr$\acute{e}$chet derivative at the elementwise minimal positive solvent is nonsingular. Although the Fr$\acute{e}$chet derivative is singular, the convergence rate is at least linear. Numerical experiments of the convergence rate are given.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2000년도 춘계 학술대회논문집
• /
• pp.178-183
• /
• 2000

An efficient Newton-GMRES algorithm is presented for computing two-dimensional steady compressible viscous flows on unstructured hybrid meshes. The scheme is designed on cell-centered finite volume method which accepts general polygonal meshes. Steady-state solution is obtained with pseudo-transient continuation strategy. The preconditioned, restarted general minimum residual(GMRES) method is employed in matrix-free form to solve the linear system arising at each Newton iteration. The incomplete LU fartorization is employed for the preconditioning of linear system. The Spalart-Allmars one equation turbulence model is fully coupled with the flow equations to simulate turbulence effect. The accuracy, efficiency and robustness of the presently developed method are demonstrated on various test problems including laminar and turbulent flows over flat plate and airfoils.

• Journal of applied mathematics & informatics
• /
• 제16권1_2호
• /
• pp.195-205
• /
• 2004

To the unconstrained programme of non-convex function, this article give a modified BFGS algorithm. The idea of the algorithm is to modify the approximate Hessian matrix for obtaining the descent direction and guaranteeing the efficacious of the quasi-Newton iteration pattern. We prove the global convergence properties of the algorithm associating with the general form of line search, and prove the quadratic convergence rate of the algorithm under some conditions.

• 대한전기학회:학술대회논문집
• /
• 대한전기학회 2003년도 학술회의 논문집 정보 및 제어부문 B
• /
• pp.875-877
• /
• 2003

This paper use fault location algorithm for single-phase-to-ground faults on the teed circuit of a parallel transmission line that use only local end voltage and current information. When Newton-Raphson iteration method is used, the Initial value may cause error or cause not suitable result. Suggested new calculation model uses NVP methodology, which is one of the fault tolerance technology to solve this problem. EMTP simulation result has shown effectiveness of the algorithm under various conditions.

• 응용통계연구
• /
• 제6권2호
• /
• pp.319-328
• /
• 1993

뉴톤-랩슨 반복법이 최우추정량에 접근하는 비율이 초기값에 따라 가속화함을 보았다. 그러 므로 최우추정량을 구하기 어려운 경우에 통계적 목적 - Bahadur 효율, 콰지(Quasi) 우도비 검정 통계량의 점근분포, Bartlett 정정계수(correction factor)등 - 에 따라 뉴톤-랩슨 반복 의 횟수를 정하여 쓸 수 있다.

• 충청수학회지
• /
• 제25권4호
• /
• pp.669-676
• /
• 2012

In this note we extend one well-known iteration formula to derive new third-order methods for solving a nonlinear equation. The comparison with other methods is given.

• Journal of applied mathematics & informatics
• /
• 제4권2호
• /
• pp.341-378
• /
• 1997

We consider three classes of numerical methods for solv-ing the semi-explicit differential-algebraic equations of index 1 and higher. These methods use implicit multistep fixed stepsize methods and several iterative processes including simple iteration, full a2nd modified Newton iteration. For these methods we prove convergence theorems and derive error estimates. We consider different ways of choosing initial approximations for these iterative methods and in-vestigate their efficiency in theory and practice.

• Nuclear Engineering and Technology
• /
• 제54권8호
• /
• pp.3059-3072
• /
• 2022

A Jacobian-Free Newton Krylov Two-Nodal Coarse Mesh Finite Difference algorithm based on Nodal Expansion Method (NEM_TNCMFD_JFNK) is successfully developed and proposed to solve the three-dimensional (3D) and multi-group reactor physics models. In the NEM_TNCMFD_JFNK method, the efficient JFNK method with the Modified Incomplete LU (MILU) preconditioner is integrated and applied into the discrete systems of the NEM-based two-node CMFD method by constructing the residual functions of only the nodal average fluxes and the eigenvalue. All the nonlinear corrective nodal coupling coefficients are updated on the basis of two-nodal NEM formulation including the discontinuity factor in every few newton steps. All the expansion coefficients and interface currents of the two-node NEM need not be chosen as the solution variables to evaluate the residual functions of the NEM_TNCMFD_JFNK method, therefore, the NEM_TNCMFD_JFNK method can greatly reduce the number of solution variables and the computational cost compared with the JFNK based on the conventional NEM. Finally the NEM_TNCMFD_JFNK code is developed and then analyzed by simulating the representative PWR MOX/UO₂ core benchmark, the popular NEACRP 3D core benchmark and the complicated full-core pin-by-pin homogenous core model. Numerical solutions show that the proposed NEM_TNCMFD_JFNK method with the MILU preconditioner has the good numerical accuracy and can obtain higher computational efficiency than the NEM-based two-node CMFD algorithm with the power method in the outer iteration and the Krylov method using the MILU preconditioner in the inner iteration, which indicates the NEM_TNCMFD_JFNK method can serve as a potential and efficient numerical tool for reactor neutron diffusion analysis module in the JFNK-based multiphysics coupling application.