• Nuclear Engineering and Technology
• /
• v.55 no.7
• /
• pp.2712-2722
• /
• 2023

For an accurate and efficient time-dependent Monte Carlo (TDMC) neutron transport analysis, several advanced methods are newly developed and implemented in the Seoul National University Monte Carlo code, McCARD. For an efficient control of the neutron population, a dynamic weight window method is devised to adjust the weight bounds of the implicit capture in the time bin-by-bin TDMC simulations. A moving geometry module is developed to model a continuous insertion or withdrawal of a control rod. Especially, the history-based batch method for the TDMC calculations is developed to predict the unbiased variance of a bin-wise mean estimate. The developed methods are verified for three-dimensional problems in the C5G7-TD benchmark, showing good agreements with results from a deterministic neutron transport analysis code, nTRACER, within the statistical uncertainty bounds. In addition, the TDMC analysis capability implemented in McCARD is demonstrated to search the optimum detector positions for the pulsed-neutron-source experiments in the Kyoto University Critical Assembly and AGN201K.

• Korean Journal of Hydrosciences
• /
• v.2
• /
• pp.85-98
• /
• 1991

This paper is focused on the development of the RIVNET1 model, which is a dynamic wave model, for flood analysis in dendritic river networks with arbitrary cross-sections. This model adopted the $-point implicit RDM and utilized a relaxation algorithim in order to solve the governing equations. The double-sweep method was used to reduce the C.P.U. time to solve the matrix system of the model. This model is applied the analyze flood waves of the Ohid river in the U.S.A. and the Keum river in Korea. The results of analysis obtained from this model are compared with those of the DWOPER and observed data.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.3 no.3
• /
• pp.39-45
• /
• 2004

In rigid-plastic finite element method, there is a heavy computation time and convergence problem. In this study, static-explicit rigid-plastic finite element method will be introduced. This method is the way that restrict the convergence interval. In result, convergence problem and computation time due to large non-linearity in the existing numerical analysis method were no longer a critical problem. 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.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.4
• /
• pp.510-519
• /
• 2015

Concerning the issue of high-dimensions, hybrid uncertainties of randomness and intervals including implicit and highly nonlinear limit state function, reliability analysis based on the hybrid uncertainty reliability mode combining with back propagation neural network (HU-BP neural network) is proposed in this paper. Random variables and interval variables are as input layer of the neural network, after the training and approximation of the neural network, the response variables are obtained through the output layer. Reliability index is calculated by solving the optimization model of the most probable point (MPP) searching in the limit state band. Two numerical cases are used to demonstrate the method proposed in this paper, and finally the method is employed to solving an engineering problem of the aerospace friction plate. For this high nonlinear, small failure probability problem with interval variables, this method could achieve a good analysis result.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.8
• /
• pp.2524-2539
• /
• 1996

A design process and an axisymmetric boundary integral equation method for precise structural analysis of the aircraft gas turbine rotor disk were developed. This axisymmetric boundary integral equation method for stress and steady-state thermal analysis was improved in solution accuracy by appling an implicit technique for Cauchy principal value evaluation, a subelement technique for weak singular integral evaluation and a double exponentical integral technoque for internal point solution near boundary surfaces. Stresses, temperatures, low cycle fatigue lifes and critical speeds for the turbine rotor disk of the thrust 1421 N class turbojet engine were analysed in a pratical calculation model problem.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.1
• /
• pp.79-87
• /
• 2013

An isogeometric topological shape optimization method is developed using the level sets and Heaviside enrichments. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set functions, which facilitates to handle complicated topological shape changes. The Heaviside enrichment improves the isogeometric analysis by adding some enrichment functions to model the internal boundaries. The proposed topological shape optimization method has several benefits: exact geometric models can be obtained using the isogeometric approach and the limitation of tensor-product patches can be overcome using the Heaviside enrichments to represent the internal voids. Even in a single patch, discontinuous displacement fields as well as smooth stress field can be obtained. Since the level sets offer the implicit moving boundary inside the domain, it is easy to represent the topological shape variations in the isogeometric analysis using Heaviside enrichments.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.474-479
• /
• 2004

According to our previous study, it is confirmed that the Petrov-Galerkin natural element method (PGNEM) completely resolves the numerical integration inaccuracy in the conventional Bubnov-Galerkin natural element method (BG-NEM). This paper is an extension of PG-NEM to two-dimensional nonlinear dynamic problem. For the analysis, a constant average acceleration method and a linearized total Lagrangian formulation is introduced with the PG-NEM. At every time step, the grid points are updated and the shape functions are reproduced from the relocated nodal distribution. This process enables the PG-NEM to provide more accurate and robust approximations. The representative numerical experiments performed by the test Fortran program, and the numerical results confirmed that the PG-NEM effectively and accurately approximates the nonlinear dynamic problem.

• Wind and Structures
• /
• v.27 no.3
• /
• pp.149-161
• /
• 2018

A reliability analysis method is proposed in this paper based on the maximum entropy (MaxEnt) principle in which constraints are specified in terms of the fractional moments instead of integer moments. Then a multiplicative dimensional reduction method (M-DRM) is introduced to compute the fractional moments. The method is applicable for both explicit and implicit limit state functions of complex structures. After two examples illustrate the accuracy and efficiency of this method in comparison to the Monte Carlo simulation (MCS), the method is used to analyze the flutter reliability of suspension bridge. The results show that the empirical formula method in which the limit state function is explicitly represented as a function of variables is only a too conservative estimate for flutter reliability analysis but is not accurate adequately. So it is not suitable for reliability analysis of bridge flutter. The actual flutter reliability analysis should be conducted based on a finite element method in which limit state function is implicitly represented as a function of variables. The proposed M-DRM provide an alternate and efficient way to analyze a much more complicated flutter reliability of long span suspension bridge.

• Steel and Composite Structures
• /
• v.43 no.1
• /
• pp.1-17
• /
• 2022

The free and forced nonlinear dynamic behaviors of Porous Functionally Graded Material (PFGM) plates are examined by means of a High-Order Implicit Algorithm (HOIA). The formulation is developed using the Third-order Shear Deformation Theory (TSDT). Unlike previous works, the formulation is written without resorting to any homogenization technique neither rule of mixture nor considering FGM as a laminated composite, and the distribution of the porosity is assumed to be gradually variable through the thickness of the PFGM plates. Using the Hamilton principle, we establish the governing equations of motion. The Finite Element Method (FEM) is used to compute approximations of the resulting equations; FEM is adopted using a four-node quadrilateral finite element with seven Degrees Of Freedom (DOF) per node. Nonlinear equations are solved by a HOIA. The accuracy and the performance of the proposed approach are verified by presenting comparisons with literature results for vibration natural frequencies and dynamic response of PFGM plates under external loading. The influences of porosity volume fraction, porosity distribution, slenderness ratio and other parameters on the vibrations of PFGM plate are explored. The results demonstrate the significant impact of different physical and geometrical parameters on the vibration behavior of the PFGM plate.

• Nuclear Engineering and Technology
• /
• v.40 no.6
• /
• pp.477-488
• /
• 2008

The objective of the paper is to analyze the thermally induced density wave oscillations in water cooled boiling water reactors. A transient thermal hydraulic model is developed with a characteristics-based implicit finite-difference scheme to solve the nonlinear mass, momentum and energy conservation equations in a time-domain. A two-phase flow was simulated with a one-dimensional homogeneous equilibrium model. The model treats the boundary conditions naturally and takes into account the compressibility effect of the two-phase flow. The axial variation of the heat flux profile can also be handled with the model. Unlike the method of characteristics analysis, the present numerical model is computationally inexpensive in terms of time and works in a Eulerian coordinate system without the loss of accuracy. The model was validated against available benchmarks. The model was extended for the purpose of studying the flow-induced density wave oscillations in forced circulation and natural circulation boiling water reactors. Various parametric studies were undertaken to evaluate the model's performance under different operating conditions. Marginal stability boundaries were drawn for type-I and type-II instabilities in a dimensionless parameter space. The significance of adiabatic riser sections in different boiling reactors was analyzed in detail. The effect of the axial heat flux profile was also investigated for different boiling reactors.