• Title/Summary/Keyword: Iteration Method

Search Result 1,145, Processing Time 0.036 seconds

Advanced two-level CMFD acceleration method for the 3D whole-core high-fidelity neutron adjoint transport calculation

  • Zhu, Kaijie;Hao, Chen;Xu, Yunlin
    • Nuclear Engineering and Technology
    • /
    • v.53 no.1
    • /
    • pp.30-43
    • /
    • 2021
  • In the 2D/1D method, a global adjoint CMFD based on the generalized equivalence theory is built to synthesize the 2D radial MOC adjoint and 1D axial NEM adjoint calculation and also to accelerate the iteration convergence of 3D whole-core adjoint transport calculation. Even more important, an advanced yet accurate two-level (TL) CMFD acceleration technique is proposed, in which an equivalent one-group adjoint CMFD is established to accelerate the multi-group adjoint CMFD and then to accelerate the 3D whole-core adjoint transport calculation efficiently. Based on these method, a new code is developed to perform 3D adjoint neutron flux calculation. Then a set of VERA and C5G7 benchmark problems are chosen to verify the capability of the 3D adjoint calculations and the effectiveness of TL CMFD acceleration. The numerical results demonstrate that acceptable accuracy of 2D/1D adjoint calculations and superior acceleration of TL CMFD are achievable.

Numerical Calculation of λ-Mode of the Diffusion Equation (수치해법을 이용한 중성자 확산방정식의λ-Mode 계산)

  • Noh, T.W.;Oh, S.K.;Kim, S.Y.;Kim, C.H.
    • Nuclear Engineering and Technology
    • /
    • v.19 no.4
    • /
    • pp.310-316
    • /
    • 1987
  • A successive iteration method to calculate the λ-modes of the diffusion equation was developed. The 2-group, 3-dimensional computer code MOGEN was developed to implement this method, The accuracy of the method was demonstrated using 2-dimensional bare homogeneous rectangular reactor. The numerical solution shows good agreement with the analytic solution in terms of eigenvalue and eigenfunction As for the standard CANDU-600 reactor, the 2-dimensional modes were generated and these represent the conventional mode characteristics well. Finally, application of theλ-mode in reactor engineering problems is described briefly.

  • PDF

Inverse Airfoil Design for Wind Turbine (역설계 기법을 이용한 풍력터빈 에어포일 형상 설계)

  • Ryu, Ki-Wahn;Park, Myoung-Ho
    • Journal of Wind Energy
    • /
    • v.4 no.2
    • /
    • pp.55-60
    • /
    • 2013
  • The mathematical implementation for inverse airfoil design of wind turbines is presented using vortex panel method based on assumptions of the two-dimensional incompressible potential flow. The vortex panel method employs linear distribution of the vortex strength to obtain the well converged solution. Stream function is adopted to get the basic formula for the inverse airfoil design, and a symmetric seed airfoil is given for initial data of the iteration approach. The final airfoil shape has been compared with the original airfoil shape for validation of the mathematical procedure.

Modified Single Loop Single Vector Method for Stability and Efficiency Improvement in Reliability-Based Design Optimization (신뢰성기반 최적설계에서 수치적 안정성과 효율성의 개선을 위해 수정된 Single Loop Single Vector 방법)

  • Kim, Bong-Jae;Lee, Jae-Ohk;Yang, Young-Soon
    • Journal of the Computational Structural Engineering Institute of Korea
    • /
    • v.18 no.1
    • /
    • pp.51-59
    • /
    • 2005
  • SLSV (single loop single vector) method is to solve the excessive computational cost problem in RBDO (reliability-based design optimization) by decoupling the nested iteration loops. However, the practical use of SLSV method to RBDO case is limited by the instability or inaccuracy of the method since it often diverges or converges to a wrong solution. Thus, in this paper, a new modified SLSV method is proposed. This method improves its convergence capability effectively by utilizing Inactive Design and Active MPP Design together with modified HMV (hybrid mean value) method. The usefulness of the proposed method is also verified through numerical examples.

A Development of Program on the Hydraulic Calculation in Sprinkler System Based on the Piping Network Analysis Method (배관망 해석 방법을 이용한 스프링클러 시스템의 수리계산 프로그램 개발)

  • 송철강;이명호;강계명
    • Fire Science and Engineering
    • /
    • v.16 no.1
    • /
    • pp.24-29
    • /
    • 2002
  • The purpose of this study is developing the computer program for hydraulic design sprinkler systems have performed the means for the general use of network analysis method. The computer program is based on the theoretical concepts of the related Hazen-Williams equations, a modified Bernoulli equations, and the Hardy Cross method of pipe network analysis. Looped piping calculations are solved by using either the Hardy Cross method or the other iteration methods. While the other methods are solved using simultaneous equations, the Hardy Cross method is concerned with one loop at a time using reiterative process. Due to its simplicity the Hardy Cross method will be the primary method described in this thesis. The purpose of this study is to develope hydraulic calculation program by using algorithm for network analysis method. The development of computer program for the hydraulic design of sprinkler systems will perform the means in the performance-based sprinkler system design.

CT Image Reconstruction of Wood Using Ultrasound Velocities II - Determination of the Initial Model Function of the SIRT Method -

  • Kim, Kwang-Mo;Lee, Jun-Jae
    • Journal of the Korean Wood Science and Technology
    • /
    • v.33 no.5 s.133
    • /
    • pp.29-37
    • /
    • 2005
  • A previous study verified that the SIRT (simultaneous iterative reconstruction technique) method is more efficient than the back-projection method as a CT algorithm for wood. However, it was expected that the determination of the initial model function of the SIRT method would influence the quality of CT image. Therefore, in this study, we intended to develop a technique that could be used to determine an adequate initial model function. For this purpose, we proposed several techniques, and for each technique we examined the effects of the initial model function on the average errors and the CT image at each iteration. Through this study, it was shown that the average error was decreased and the image quality was improved using the proposed techniques. This tendency was most pronounced when the back-projection method was used to determine the initial model function. From the results of this study, we drew the following conclusions: 1) The initial model function of the SIRT method should be determined with careful attention, and 2) the back-projection method efficiently determines the initial model function of the SIRT method.

Waveform Relaxation Method for Reactor Transient Analysis (원자로 천이해석을 위한 파형완화법)

  • Park, Keon-Woo;Co, Nam-Zin
    • Nuclear Engineering and Technology
    • /
    • v.27 no.6
    • /
    • pp.845-852
    • /
    • 1995
  • We investigate the concurrent solution of differential equations by the waveform relaxation (WR) method, an iterative method for analyzing linear and nonlinear dynamical systems in the time do-main. The method, at each iteration, decomposes the dynamical system into several subsystems, each of which is analyzed for the entire given time interval. The method, when efficiently implemented, results in algorithms with a highly parallelizable concurrent fraction. In this paper, the waveform relaxation method is introduced and applied to two types of reactor dynamics problems. It is concluded that the U method can be applied to reactor dynamics equations, but that its parallel performance on the KMRR dynamics is only modest.

  • PDF

A parameter calibration method for PFC simulation: Development and a case study of limestone

  • Xu, Z.H.;Wang, W.Y.;Lin, P.;Xiong, Y.;Liu, Z.Y.;He, S.J.
    • Geomechanics and Engineering
    • /
    • v.22 no.1
    • /
    • pp.97-108
    • /
    • 2020
  • The time-consuming and less objectivity are the main problems of conventional micromechanical parameters calibration method of Particle Flow Code simulations. Thus this study aims to address these two limitation of the conventional "trial-and-error" method. A new calibration method for the linear parallel bond model (CM-LPBM) is proposed. First, numerical simulations are conducted based on the results of the uniaxial compression tests on limestone. The macroscopic response of the numerical model agrees well with the results of the uniaxial compression tests. To reduce the number of the independent micromechanical parameters, numerical simulations are then carried out. Based on the results of the orthogonal experiments and the multi-factor variance analysis, main micromechanical parameters affecting the macro parameters of rocks are proposed. The macro-micro parameter functions are ultimately established using multiple linear regression, and the iteration correction formulas of the micromechanical parameters are obtained. To further verify the validity of the proposed method, a case study is carried out. The error between the macro mechanical response and the numerical results is less than 5%. Hence the calibration method, i.e., the CM-LPBM, is reliable for obtaining the micromechanical parameters quickly and accurately, providing reference for the calibration of micromechanical parameters.

Diffusion synthetic acceleration with the fine mesh rebalance of the subcell balance method with tetrahedral meshes for SN transport calculations

  • Muhammad, Habib;Hong, Ser Gi
    • Nuclear Engineering and Technology
    • /
    • v.52 no.3
    • /
    • pp.485-498
    • /
    • 2020
  • A diffusion synthetic acceleration (DSA) technique for the SN transport equation discretized with the linear discontinuous expansion method with subcell balance (LDEM-SCB) on unstructured tetrahedral meshes is presented. The LDEM-SCB scheme solves the transport equation with the discrete ordinates method by using the subcell balances and linear discontinuous expansion of the flux. Discretized DSA equations are derived by consistently discretizing the continuous diffusion equation with the LDEM-SCB method, however, the discretized diffusion equations are not fully consistent with the discretized transport equations. In addition, a fine mesh rebalance (FMR) method is devised to accelerate the discretized diffusion equation coupled with the preconditioned conjugate gradient (CG) method. The DSA method is applied to various test problems to show its effectiveness in speeding up the iterative convergence of the transport equation. The results show that the DSA method gives small spectral radii for the tetrahedral meshes having various minimum aspect ratios even in highly scattering dominant mediums for the homogeneous test problems. The numerical tests for the homogeneous and heterogeneous problems show that DSA with FMR (with preconditioned CG) gives significantly higher speedups and robustness than the one with the Gauss-Seidel-like iteration.

Resonance Elastic Scattering and Interference Effects Treatments in Subgroup Method

  • Li, Yunzhao;He, Qingming;Cao, Liangzhi;Wu, Hongchun;Zu, Tiejun
    • Nuclear Engineering and Technology
    • /
    • v.48 no.2
    • /
    • pp.339-350
    • /
    • 2016
  • Based on the resonance integral (RI) tables produced by the NJOY program, the conventional subgroup method usually ignores both the resonance elastic scattering and the resonance interference effects. In this paper, on one hand, to correct the resonance elastic scattering effect, RI tables are regenerated by using the Monte Carlo code, OpenMC, which employs the Doppler broadening rejection correction method for the resonance elastic scattering. On the other hand, a fast resonance interference factor method is proposed to efficiently handle the resonance interference effect. Encouraging conclusions have been indicated by the numerical results. (1) For a hot full power pressurized water reactor fuel pin-cell, an error of about +200 percent mille could be introduced by neglecting the resonance elastic scattering effect. By contrast, the approach employed in this paper can eliminate the error. (2) The fast resonance interference factor method possesses higher precision and higher efficiency than the conventional Bondarenko iteration method. Correspondingly, if the fast resonance interference factor method proposed in this paper is employed, the $k_{inf}$ can be improved by ~100 percent mille with a speedup of about 4.56.