• Nuclear Engineering and Technology
• /
• v.51 no.1
• /
• pp.163-169
• /
• 2019

The robust controller synthesis and analysis of the water level process in the U-tube system generator (UTSG) is addressed in this paper. The parameter uncertainties of the steam generator (SG) are modeled as multiplicative perturbations which are normalized by designing suitable weighting functions. The relative errors of the nominal SG model with respect to the other operating power level models are employed to specify the weighting functions for normalizing the plant uncertainties. Then, a robust controller is designed based on ${\mu}$-synthesis and D-K iteration, and its stability robustness is verified over the whole range of power operations. A gain-scheduled controller with $H_{\infty}$-synthesis is also designed to compare its robustness with the proposed controller. The stability analysis is accomplished and compared with the previous QFT design. The ${\mu}$-analysis of the system shows that the proposed controller has a favorable stability robustness for the whole range of operating power conditions. The proposed controller response is simulated against the power level deviation in start-up and shutdown stages and compared with the other concerning controllers.

• Nuclear Engineering and Technology
• /
• v.54 no.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.

• Journal of Biomedical Engineering Research
• /
• v.38 no.5
• /
• pp.219-226
• /
• 2017

Penalized-likelihood (PL) reconstruction methods for transmission tomography are known to provide improved image quality for reduced dose level by efficiently smoothing out noise while preserving edges. Unfortunately, however, most of the edge-preserving penalty functions used in conventional PL methods contain at least one free parameter which controls the shape of a non-quadratic penalty function to adjust the sensitivity of edge preservation. In this work, to avoid difficulties in finding a proper value of the free parameter involved in a non-quadratic penalty function, we propose a new adaptive method of space-variant smoothing with a simple quadratic penalty function. In this method, the smoothing parameter is adaptively selected for each pixel location at each iteration by using the image roughness measured by a pixel-wise standard deviation image calculated from the previous iteration. The experimental results demonstrate that our new method not only preserves edges, but also suppresses noise well in monotonic regions without requiring additional processes to select free parameters that may otherwise be included in a non-quadratic penalty function.

• Proceedings of the KSEEG Conference
• /
• 2002.04a
• /
• pp.167-169
• /
• 2002

The extended Born, or localized nonlinear approximation of integral equation (IE) solution has been applied to inverting single-hole electromagnetic (EM) data using a cylindrically symmetric model. The extended Born approximation is less accurate than a full solution but much superior to the simple Born approximation. When applied to the cylindrically symmetric model with a vertical magnetic dipole source, however, the accuracy of the extended Born approximation is greatly improved because the electric field is scalar and continuous everywhere. One of the most important steps in the inversion is the selection of a proper regularization parameter for stability. Occam's inversion (Constable et al., 1987) is an excellent method for obtaining a stable inverse solution. It is extremely slow when combined with a differential equation method because many forward simulations are needed but suitable for the extended Born solution because the Green's functions, the most time consuming part in IE methods, are repeatedly re-usable throughout the inversion. In addition, the If formulation also readily contains a sensitivity matrix, which can be revised at each iteration at little expense. The inversion algorithm developed in this study is quite stable and fast even if the optimum regularization parameter Is sought at each iteration step. Tn this paper we show inversion results using synthetic data obtained from a finite-element method and field data as well.

• The Bulletin of The Korean Astronomical Society
• /
• v.44 no.2
• /
• pp.71.3-71.3
• /
• 2019

Previously we developed a method of coronal force-free field construction using vector potentials. In this method, the boundary normal component of the vector potential should be adjusted at every iteration step to implement the boundary normal current density, which is provided by observations. The method was a variational method in the sense that the excessive kinetic energy is removed from the system at every iteration step. The boundary condition imposing the normal current density, however, is not compatible with the variational procedure seeking for the minimum energy state, which is employed by most force-free field solvers currently being used. To resolve this problem, we have developed a totally new method of force-free field construction. Our new method uses a unique magnetic field description using two scalar functions. Our procedure is non-variational and can impose the boundary normal current density exactly. We have tested the new force-free solver for standard Low & Lou fields and Titov-Demoulin flux ropes. Our code excels others in both examples, especially in Titov-Demoulin flux ropes, for which most codes available now yield poor results. Application to a real active region will also be presented.

• Smart Structures and Systems
• /
• v.20 no.1
• /
• pp.99-114
• /
• 2017

In cognitive science, it is illustrated how the collective opinions of a group of individuals answers to questions involving quantity estimation. One example of this approach is introduced in this article as Star Graph (SG) algorithm. This graph describes the details of communication among individuals to share their information and make a new decision. A new labyrinthine network of neighbors is defined in the decision-making process of the algorithm. In order to prevent getting trapped in local optima, the neighboring networks are regenerated in each iteration of the algorithm. In this algorithm, the normal distribution is utilized for a group of agents with the best results (guidance group) to replace the existing infeasible solutions. Here, some new functions are introduced to provide a high convergence for the method. These functions not only increase the local and global search capabilities but also require less computational effort. Various benchmark functions and engineering problems are examined and the results are compared with those of some other algorithms to show the capability and performance of the presented method.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.34 no.3
• /
• pp.243-251
• /
• 2016

LiDAR data processing steps include noise removal, filtering, classification, segmentation, shape recognition, modeling, and quality assessment. This paper focuses on modeling and accuracy evaluation of 3D objects with curved surfaces. The appropriate modeling functions were determined by analyzing surface patch shape. Existing methods for modeling curved surface features require linearization, initial approximation, and iteration of the non-linear functions. However, proposed method could directly estimate the unknown parameters of the modeling functions. The results demonstrate feasibility of the proposed method. The proposed method was applied to the simulated and real building data of hemi-spherical and semi-cylindrical surfaces. The parameters and accuracy of the modeling functions were estimated. It is expected that the proposed method would contribute to automatic modeling of various objects.

• Computers and Concrete
• /
• v.21 no.4
• /
• pp.441-449
• /
• 2018

FE models for complex or large-scaled structures that need detailed modeling of structural components are usually constructed using commercial analysis softwares. Updating of such FE model by conventional sensitivity-based methods is difficult since repeated computation for perturbed parameters and manual calculations are needed to obtain sensitivity matrix in each iteration. In this study, an FE model updating procedure avoiding such difficulties by using response surface (RS) method and a Pareto-based multiobjective optimization (MOO) was formulated and applied to FE models constructed with a commercial analysis package. The test building is a low-rise reinforced concrete building that has been seismically retrofitted. Dynamic properties of the building were extracted from vibration tests performed before and after the seismic retrofits, respectively. The elastic modulus of concrete and masonry, and spring constants for the expansion joint were updated. Two RS functions representing the errors in the natural frequencies and mode shape, respectively, were obtained and used as the objective functions for MOO. Among the Pareto solutions, the best compromise solution was determined using the TOPSIS (Technique for Order of Preference by Similarity to Ideal Solution) procedure. A similar task was performed for retrofitted building by taking the updating parameters as the stiffness of modified or added members. Obtained parameters of the existing building were reasonably comparable with the current code provisions. However, the stiffness of added concrete shear walls and steel section jacketed members were considerably lower than expectation. Such low values are seemingly because the bond between new and existing concrete was not as good as the monolithically casted members, even though they were connected by the anchoring bars.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.21 no.3
• /
• pp.167-180
• /
• 2017

According to the fact that the low convergence level of the complete elliptic integral of the first kind for the modulus which having values approach to one. In this paper we propose novelty of the complete elliptic integral which having new infinite series that consists of new modulus introduced as own modulus function. We apply scheme of iteration by substituting the common modulus with own modulus function into the new infinite series. We obtained so many new exact formulas of the complete elliptic integral derived from this method correspond to the number of iterations. On the other hand, it has been also obtained a lot of new transformation functions with the corresponding own modulus functions. The calculation results show that the enhancement of the number of significant figures of the new infinite series of the complete elliptic integral of the first kind corresponds to the level of quadratic convergence.

• Korean Journal of Optics and Photonics
• /
• v.6 no.3
• /
• pp.245-256
• /
• 1995

Rugate filters of inhomogeneous refractive index were designed using the Fourier transform and the effect of reflectance, stop bandwidth, optical thickness, and Q function on the rugate filter was investigated. An iterative correction process using a merit function was employed to fit an initial design to the target spectrum. Three Q functions derived by Sossi, Bovard, and Fabricius, respectively, were compared in terms of the number of iteration, merit function, and optimum optical thickness. The result shows that after a number of iterations the Q functions by Bovard and Fabricius produce high rejection rugate filters closer to the target spectrum than the Sossi's Q function and the optimal optical thickness is determined by the stop-band width of the rugate filter. ilter.