• Nuclear Engineering and Technology
• /
• v.45 no.2
• /
• pp.191-202
• /
• 2013

A numerical analysis of thermal stratification in the upper plenum of the MONJU fast breeder reactor was performed. Calculations were performed for a 1/6 simplified model of the MONJU reactor using the commercial code, CFX-13. To better resolve the geometrically complex upper core structure of the MONJU reactor, the porous media approach was adopted for the simulation. First, a steady state solution was obtained and the transient solutions were then obtained for the turbine trip test conducted in December 1995. The time dependent inlet conditions for the mass flow rate and temperature were provided by JAEA. Good agreement with the experimental data was observed for steady state solution. The numerical solution of the transient analysis shows the formation of thermal stratification within the upper plenum of the reactor vessel during the turbine trip test. The temporal variations of temperature were predicted accurately by the present method in the initial rapid coastdown period (~300 seconds). However, transient numerical solutions show a faster thermal mixing than that observed in the experiment after the initial coastdown period. A nearly homogenization of the temperature field in the upper plenum is predicted after about 900 seconds, which is a much shorter-term thermal stratification than the experimental data indicates. This discrepancy may be due to the shortcoming of the turbulence models available in the CFX-13 code for a natural convection flow with thermal stratification.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2004.04a
• /
• pp.57-62
• /
• 2004

Streamline simulation researches have been extensively accomplished due to the swiftness of computation and the reduction of numerical dispersion. In this study, we developed a streamline simulation model using a semianalytical solution of ID transport equation. To validate accuracy of the developed model, we compared simulation results of contaminant transport, which were acquired by streamline simulation models using an analytical solution, a numerical solution, and a semianalytical solution. The developed model using the semianalytical solution matched well with the model using an analytical solution. However, streamline simulation model using a numerical solution showed numerical dispersion. For an advection-dominant flow, there was little difference in the simulation results between the developed model and tile analytical model, but the differences between the analytical model and the numerical model were cleary shown. From the comparison of computing time we know that the streamline simulation using the semianalytical solution is 2-60 times as fast as the streamline simulation using the numerical solution.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.2 no.1
• /
• pp.27-39
• /
• 1998

We describe a fast numerical method for non-separable elliptic equations in self-adjoin form on irregular adaptive domains. One of the most successful results in numerical PDE is developing rapid elliptic solvers for separable EPDEs, for example, Fourier transformation methods for Poisson problem on a square, however, it is known that there is no rapid elliptic solvers capable of solving a general nonseparable problems. It is the purpose of this paper to present an iterative solver for linear EPDEs in self-adjoint form. The scheme discussed in this paper solves a given non-separable equation using a sequence of solutions of Poisson equations, therefore, the most important key for such a method is having a good Poison solver. High performance is achieved by using a fast high-order adaptive Poisson solver which requires only about 500 floating point operations per gridpoint in order to obtain machine precision for both the computed solution and its partial derivatives. A few numerical examples have been presented.

• Journal of computational fluids engineering
• /
• v.17 no.4
• /
• pp.41-48
• /
• 2012

A numerical analysis of thermal stratification in the upper plenum of the MONJU fast breeder reactor was performed. Calculations were performed for a 1/6 simplified model of the MONJU reactor using the commercial code, CFX-13. To better resolve the geometrically complex upper core structure of the MONJU reactor, the porous media approach was adopted for the simulation. First, a steady state solution was obtained and the transient solutions were then obtained for the turbine trip test conducted in December 1995. The time dependent inlet conditions for the mass flow rate and temperature were provided by JAEA. Good agreement with the experimental data was observed for steady state solution. The numerical solution of the transient analysis shows the formation of thermal stratification within the upper plenum of the reactor vessel during the turbine trip test. The temporal variations of temperature were predicted accurately by the present method in the initial rapid coastdown period (~300 seconds). However, transient numerical solutions show a faster thermal mixing than that observed in the experiment after the initial coastdown period. A nearly homogenization of the temperature field in the upper plenum is predicted after about 900 seconds, which is a much shorter-term thermal stratification than the experimental data indicates. This discrepancy is due to the shortcoming of the turbulence models available in the CFX-13 code for a natural convection flow with thermal stratification.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.7 no.2
• /
• pp.375-389
• /
• 2015

An improved boundary element method is presented for numerical analysis of hydrodynamic behavior of marine structures. A new algorithm for numerical solution of the finite depth free-surface Green function in three dimensions is developed based on multiple series representations. The whole range of the key parameter R/h is divided into four regions, within which different representation is used to achieve fast convergence. The well-known epsilon algorithm is also adopted to accelerate the convergence. The critical convergence criteria for each representation are investigated and provided. The proposed method is validated by several well-documented benchmark problems.

• Journal of Institute of Control, Robotics and Systems
• /
• v.18 no.9
• /
• pp.823-829
• /
• 2012

An excitation table is commonly used for vibration and ride tests for parts or assemblies of automobiles, aircrafts, or other heavy systems. The authors have analyzed several kinematic properties of an excitation table that is under development for heavy transport vehicles. It consists of one table and 7 linear hydraulic actuators. The authors have performed mobility analysis, inverse kinematics, forward kinematics, and singularity analysis. Especially, we have proposed a fast forward kinematic solution considering the limited motion of the excitation table. On the assumption that the motion variables such as rotation angles and displacements are small, the forward kinematic problem is converted to the observer problem of a linear system. This provides a fast solution. Also we have verified that there are no singularity points in the working range by numerical analysis.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.19-22
• /
• 1998

Automatic mesh generation procedures applied to industrial now problems lead to complex mesh topologies where usually no special considerations to mesh resolution are taken. In the present study a fast and flexible solution algorithm in combination with generalized higher order discretization schemes is presented and its application to intake port calculation is demonstrated.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.14 no.4
• /
• pp.201-210
• /
• 2010

We propose a new robust and accurate method for the numerical solution of medical image segmentation. The modified Allen-Cahn equation is used to model the boundaries of the image regions. Its numerical algorithm is based on operator splitting techniques. In the first step of the splitting scheme, we implicitly solve the heat equation with the variable diffusive coefficient and a source term. Then, in the second step, using a closed-form solution for the nonlinear equation, we get an analytic solution. We overcome the time step constraint associated with most numerical implementations of geometric active contours. We demonstrate performance of the proposed image segmentation algorithm on several artificial as well as real image examples.

• Smart Structures and Systems
• /
• v.18 no.1
• /
• pp.183-196
• /
• 2016

In this paper a nonlinear, bistable, single degree of freedom system is considered. It consists of a Duffing oscillator externally excited by a non-resonant, harmonic force. A customized perturbation scheme is proposed to achieve an approximate expression for periodic solutions. It is based on the evaluation of the quasi-steady (slow) solution, and then on a variable change followed by two perturbation steps which aim to capture the fast, decaying contribution of the response. The reconstructed solution, given by the sum of the slow and fast contributions, is in a good agreement with the one obtained by numerical integration.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.12
• /
• pp.2040-2047
• /
• 2001

In this paper, a singularity problem of the state transition matrix is investigated in the spatial propagation when the spatial matrix differential equation is constructed via time finite element analysis. A parametric study shows that the degree of singularity of the state transition matrix depends on the degree of flexibility of the structures. As an alternative to avoid the numerical problems due to the singularity, an analytic solution fur spatial propagation of the flexible structures is proposed. In the proposed method, the spatial properties of the structure are analytically expressed by a combination of transcendental functions. The analytic solution serves fast and accurate results by eliminating the possibility of the error accumulation caused by the boundary condition. Several numerical examples are shown to validate the effectiveness of the proposed methods.