• Nuclear Engineering and Technology
• /
• v.54 no.2
• /
• pp.532-545
• /
• 2022

In-depth convergence analyses for neutronics/thermal-hydraulics (T/H) coupled calculations are performed to investigate the performance of nonlinear methods based on the Fixed-Point Iteration (FPI). A simplified neutronics-T/H coupled system consisting of a single fuel pin is derived to provide a testbed. The xenon equilibrium model is considered to investigate its impact during the nonlinear iteration. A problem set is organized to have a thousand different fuel temperature coefficients (FTC) and moderator temperature coefficients (MTC). The problem set is solved by the Jacobi and Gauss-Seidel (G-S) type FPI. The relaxation scheme and the Anderson acceleration are applied to improve the convergence rate of FPI. The performances of solution schemes are evaluated by comparing the number of iterations and the error reduction behavior. From those numerical investigations, it is demonstrated that the number of FPIs is increased as the feedback is stronger regardless of its sign. In addition, the Jacobi type FPIs generally shows a slower convergence rate than the G-S type FPI. It also turns out that the xenon equilibrium model can cause numerical instability for certain conditions. Lastly, it is figured out that the Anderson acceleration can effectively improve the convergence behaviors of FPI, compared to the conventional relaxation scheme.

• Advances in Energy Research
• /
• v.5 no.3
• /
• pp.255-275
• /
• 2017

The advantages of using natural circulation (NC) as a cooling system, has prompted the worldwide development to investigate this phenomenon more than before. The interesting application of the NC in low power experimental facilities and research reactors, highlights the obligation of study in these laminar flows. The inherent oscillations of NC between hot source and cold sink in low Grashof numbers necessitates stability analysis of cooling flow with experimental or numerical schemes. For this type of analysis, numerical methods could be implemented to desired mass, momentum and energy equations as an efficient instrument for predicting the behavior of the flow field. In this work, using the explicit, implicit and Crank-Nicolson methods, the fluid flow parameters in a natural circulation experimental test loop are obtained and the sensitivity of solving approaches are discussed. In this way, at first, the steady state and transient results from explicit are obtained and compared with experimental data. The implicit and crank-Nicolson scheme is investigated in next steps and in subsequent this research is focused on the numerical aspects of instability prediction for these schemes. In the following, the assessment of the flow behavior with coarse and fine mesh sizes and time-steps has been reported and the numerical schemes convergence are compared. For more detail research, the natural circulation of fluid was modeled by ANSYS-CFX software and results for the experimental loop are shown. Finally, the stability map for rectangular closed loop was obtained with employing the Nyquist criterion.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.6 no.3
• /
• pp.211-221
• /
• 1982

Numerical analysis has been made on the heat transfer of two dimensional turbulent channel with a slat type blockage. Especially the effects of the height of slat and Reynolds number on the heat transfer characteristics of channel wall have been investigated. The methods of accelerating the convergence of the numerical solution of governing differential equation have been also examined. Line-by-line iterative method shows higher convergence rate than point-by-point iterative method for solution of both momentum equation and energy equation. The results show that the ratio of heat transfer coefficient of the wall near the blockage to that of the fully developed flow increase with increasing the ratio of blockage to channel height and decreasing the Reynolds number. These trends of variation of heat transfer coefficient with respect to the height of slat and Reynolds number agree with those of Sparrow's experiment on the pipe flow with slat type blockage.

• Proceedings of the Korean Institute of Intelligent Systems Conference
• /
• 1995.10b
• /
• pp.71-77
• /
• 1995

In this paper are proposed robust iterative learning control(ILC) algorithms for both linear continuous time-invariant system and linear discrete-time system. In contrast to conventional methods, the proposed learning algorithms are constructed based on both time domain performance and iteration-domain performance. The convergence of the proposed learning algorithms is proved. Also, it is shown that the proposed method has robustness in the presence of external disturbances and the convergence accuracy can be improved. A numerical external disturbances and the convergence accuracy can be improved. A numerical example is provided to show the effectiveness of the proposed algorithm.

• Proceedings of the KSME Conference
• /
• 2001.06c
• /
• pp.604-609
• /
• 2001

In the sheet metal forming process, the drawbead is used to control the flow of material during the forming process. The drawbead provides proper restraining force to the material and prevents defects such as wrinkling or breakage. For these reasons, many studies for designing the effective drawbead have been conducted. For the analysis, the numerical method called the static-explicit finite element method was used. The finite element analysis code for this method has been developed and applied to the drawbead process problems. In result, convergence problem and computation time due to large non-linearity in the existing numerical analysis methods were no longer a critical problem. Futhermore, this approach could treat the contact friction problem easily by applying very small time intervals. It is expected that various results from the numerical analysis will give very useful information for the design of tools in sheet metal forming process.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.1
• /
• pp.22-30
• /
• 2008

This paper deals with the nonlinear optimal control approach to helicopter maneuver problems using the indirect method. We apply a penalty function to the deviation from a prescribed trajectory to convert the system optimality to an unconstrained optimal control problem. The resultant two-point boundary value problem has been solved by using the multiple-shooting method. This paper focuses on the effect of the number of shooting nodes and initialization methods on the numerical solution in order to define the minimum number of shooting nodes required for numerical convergence and to provide a method increasing convergence radius of the indirect method. The results of this study can provide an approach to improve numerical stability and convergence of the indirect method when we solve the optimal control problems of an inherently unstable helicopter system.

• Journal of applied mathematics & informatics
• /
• v.32 no.5_6
• /
• pp.609-619
• /
• 2014

Nonsmooth equations with finitely many maximum functions is often used in the study of complementarity problems, variational inequalities and many problems in engineering and mechanics. In this paper, we consider the global convergence methods for nonsmooth equations with finitely many maximum functions. The steepest decent method and the smoothing gradient method are used to solve the nonsmooth equations with finitely many maximum functions. In addition, the convergence analysis and the applications are also given. The numerical results for the smoothing gradient method indicate that the method works quite well in practice.

• Journal of Korean Society of Steel Construction
• /
• v.10 no.3 s.36
• /
• pp.451-461
• /
• 1998

The solution of anisotropic plate via the classical methods is limited to relatively load and boundary conditions. If these conditions are more complex, the analysis becomes increasingly tedious and even impossible. For many plate problems of considerable practical interest, analytic solutions to the governing differential equations cannot be found. Among the numerical techniques presently available, the finite difference method and the finite element method are powerful numerical methods. The objective of this paper is to compare with each numerical methods for the buckling load and modes of anisotropic composite laminated plates considering shear deformation. In applying numerical methods to solve differential equations of anisotropic plates, this study uses the finite difference method and the finite element method. In determining the eigenvalue by Finite Difference Method, this paper represent good convergence compared with Finite Element Method. Several numerical examples and buckling modes show the effectiveness of various numerical methods and they will give a guides in deciding minimum buckling load and various mode shapes.

• Steel and Composite Structures
• /
• v.29 no.6
• /
• pp.759-770
• /
• 2018

In the present investigation, by using the two numerical methods, free vibration analysis of laminated annular and annular sector plates have been studied. In order to obtain the main equations two different shell theories such as Love's shell theory and first-order shear deformation theory (FSDT) have been used for modeling. After obtaining the fundamental equations in briefly, the methods of harmonic differential quadrature (HDQ) and discrete singular convolution (DSC) are used to solve the equation of motion. Accuracy, convergence and reliability of the present HDQ and DSC methods were tested by comparing the existing results obtained by different methods in the literature. The effects of some geometric and material properties of the plates are investigated via these two methods. The advantages and accuracy of the HDQ and DSC methods have also been examined with different grid numbers and shell theory. Some results for laminated annular plates and laminated circular plates were also been supplied.

• Journal of Ocean Engineering and Technology
• /
• v.37 no.5
• /
• pp.215-224
• /
• 2023

While melting glaciers due to global warming have facilitated the development of polar routes, Arctic vessels require reliable anti-icing methods to prevent hull icing. Currently, the existing anti-icing method, i.e., the heating coil method, has disadvantages, such as disconnection and power inefficiency. Therefore, a carbon nanotube-based surface heating element method was developed to address these limitations. In this study, the numerical analysis of the surface heating element method was performed using ANSYS. The numerical analysis included conjugate heat transfer and computational fluid dynamics to consider the conduction solids and the effects of wind speed and temperature in cold environments. The numerical analysis method of the surface heating element method was validated by comparing the experimental results of the heating coil method with the numerical analysis results (under the -30 ℃ conditions). The surface heating element method demonstrated significantly higher efficiency, ranging from 56.65-80.17%, depending on the conditions compared to the heating coil method. Moreover, even under extreme environmental conditions (-45 ℃), the surface heating element method satisfied anti-icing requirements. The surface heating element method is more efficient and economical than the heating coil method. However, proper heat flux calculation for environmental conditions is required to prevent excessive design.