• Nuclear Engineering and Technology
• /
• v.48 no.2
• /
• pp.404-418
• /
• 2016

This paper presents a documentation and development method to facilitate the certification of scientific computing software used in the safety analysis of nuclear facilities. To study the problems faced during quality assurance and certification activities, a case study was performed on legacy software used for thermal analysis of a fuelpin in a nuclear reactor. Although no errors were uncovered in the code, 27 issues of incompleteness and inconsistency were found with the documentation. This work proposes that software documentation follow a rational process, which includes a software requirements specification following a template that is reusable, maintainable, and understandable. To develop the design and implementation, this paper suggests literate programming as an alternative to traditional structured programming. Literate programming allows for documenting of numerical algorithms and code together in what is termed the literate programmer's manual. This manual is developed with explicit traceability to the software requirements specification. The traceability between the theory, numerical algorithms, and implementation facilitates achieving completeness and consistency, as well as simplifies the process of verification and the associated certification.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2008.05a
• /
• pp.387-390
• /
• 2008

Electromagnetic forming (EMF) technology, which is one of the high speed forming methods, has been used for the forming process in various industry fields. Numerical approach by finite element simulation of the EMF process is presented in this study. The implicit code is used to obtain the numerical model of the time-varying currents that are discharged through the coil in order to obtain the transient magnetic forces. In addition, the body forces generated in the workpiece are used as the loading condition to analyze deformation of thin sheet metal workpiece using explicit code. Numerical approach for a dimpled shape by EMF process is carried out and the simulated results of the dimpled shape by EMF are reviewed in view of the deformed shape and formability evaluation.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05a
• /
• pp.179-184
• /
• 1996

Power density distribution and criticality of a pressurized water reactor are calculated with a Monte Carlo calculation using the MCNP code. The MCNP model is based on one-eighth core symmetry. Individual fuel assemblies are modeled with fullscope three dimensional description except grid spacer. The fuel rod is divided into eight axial segments. Core internals above and below the active fuel region is represented as coolant. After 400 cycle calculations, the system converges to a k value of 1.09151$\pm$0.00066. Fission reaction rate in each rod is also calculated to use as the source term in pressure vessel fluence calculation.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05a
• /
• pp.40-45
• /
• 1996

A new concept of equivalent homogenization is proposed. The concept employs new set of homogenized parameters: homogenized cross sections (XS) and interface matrix (IM), which relates partial currents at the celt interfaces. The idea of interface matrix generalizes the idea of discontinuity factors (DFs), proposed and developed by K.Koebke and K.Smith. The method of K.Smith can be simulated within framework of new method, while the new method approximates heterogeneous cell better in case of the steep flux gradients at the cell interfaces. The attractive shapes of new concept are: improved accuracy, simplicity of incorporation in the existing codes, equal numerical expenses in comparison to the K.Smith's approach. The new concept is useful for: (a) explicit reflector/baffle simulation; (b) control blades simulation; (c) mixed UO2/MOX core simulation, The offered model has been incorporated in the finite difference code and in the nodal code PANBOX. The numerical results show good accuracy of core calculations and insensitivity of homogenized parameters with respect to in- core conditions.

• Journal of the Korean Geotechnical Society
• /
• v.29 no.3
• /
• pp.29-40
• /
• 2013

Seismic response of single degree of freedom system supported by shallow foundation was analyzed by using nonlinear explicit finite difference element code. Numerical analysis results were verified with dynamic centrifuge test results of the same soil profile and structural dimensions with the numerical analysis model at a centrifugal acceleration of 20 g. Differences between the analysis and the test results induced by the boundary conditions of control points can be reduced by adding additional local damping to the natural born cyclic hysteretic damping of the soil strata. The analysis results show good agreement with the test results in terms of both time histories and response spectra. Thus, it can be concluded that the nonlinear explicit finite difference element code will be a useful technique for estimating seismic residual displacement, earthpressure etc. which are difficult to measure during laboratory tests and real earthquake.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.11
• /
• pp.1749-1755
• /
• 2010

The performance of power window system was decided by driving characteristics of the window regulator part and reaction by the glass run. The performance of power window system usually has been predicted by experimental methods. In this paper, an analytical method using the explicit code was suggested to overcome the limit of the experimental methods. The friction coefficient of glass run was obtained by the friction test at various conditions and the Mooney-Rivlin model was used. Also, a mechanism of window regulator consisted of the fast belt system and the slip ring elements. And, we conducted the analysis considering characteristic of a motor and obtained the lifting speed of automotive glass with high reliability

• Journal of The Korean Astronomical Society
• /
• v.38 no.2
• /
• pp.59-66
• /
• 2005

A new ion transport code for planetary ionospheric studies has been developed with consideration of velocity differences among ion species involving ion-ion collision. Most of previous planetary ionosphere models assumed that ions diffuse through non-moving ion and neutral background in order to consolidate continuity and momentum equations for ions into a simple set of diffusion equations. The simplification may result in unreliable density profiles of ions at high altitudes where ion velocities are fast and their velocity differences are significant enough to cause inaccuracy when computing ion-ion collision. A new code solves explicitly one-dimensional continuity and momentum equations for ion densities and velocities by utilizing divided Jacobian matrices in matrix inversion necessary to the Newton iteration procedure. The code has been applied to Martian nightside ionosphere models, as an example computation. The computed density profiles of $O^+,\;OH^+$, and $HCO^+$ differ by more than a factor of 2 at altitudes higher than 200 km from a simple diffusion model, whereas the density profile of the dominant ion, $O_2^+$, changes little. Especially, the density profile of $HCO^+$ is reduced by a factor of about 10 and its peak altitude is lowered by about 40 km relative to a simple diffusion model in which $HCO^+$ ions are assumed to diffuse through non-moving ion background, $O_2^+$. The computed effects of the new code on the Martian nightside models are explained readily in terms of ion velocities that were solved together with ion densities, which were not available from diffusion models. The new code should thus be expected as a significantly improved tool for planetary ionosphere modelling.

• Nuclear Engineering and Technology
• /
• v.56 no.3
• /
• pp.959-965
• /
• 2024

The High Temperature Test Unit (HTTU) was an experimental set-up to conduct separate and integral effects tests of the Pebble Bed Modular Reactor (PBMR) core. The annular core consisted of a randomly packed bed of uniform spheres. Natural convection tests using both nitrogen and helium, and forced convection tests using nitrogen, were conducted. The maximum material temperature achieved during forced convection testing was 1200 ℃. This paper presents the numerical analysis of the flow and temperature distribution for a forced convection test using 3D CFD as well as a 1D systems-CFD computer code. Several modelling approaches are possible, ranging from a fully explicit to a semi-implicit method that relies on correlations of their associated phenomena. For the comparison between codes, the analysis was performed using a porous media approach, where the conduction and radiative heat transfer were lumped together as an effective thermal conductivity and the convective heat transfer was correlated between the solid and gas phases. The results from both codes were validated against the experimental measurements. Favourable results were obtained, in particular by the systems-CFD code with minimal computational and time requirements.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.114-125
• /
• 2008

Current state and perspective of DNS of turbulence and turbulent combustion are discussed with feature trend of the fastest supercomputer in the world. Based on the perspective of DNS of turbulent combustion, possibility of perfect simulations of IC engine is shown. In 2020, the perfect simulation will be realized with 30 billion grid points by 1EXAFlops supercomputer, which requires 4 months CPU time. The CPU time will be reduced to about 4 days if several developments were achieved in the current fundamental researches. To shorten CPU time required for DNS of turbulent combustion, two numerical methods are introduced to full-explicit full-compressible DNS code. One is compact finite difference filter to reduce spatial resolution requirements and numerical oscillations in small scales, and another is well-known point-implicit scheme to avoid quite small time integration of the order of nanosecond for fully explicit DNS. Availability and accuracy of these numerical methods have been confirmed carefully for auto-ignition, planar laminar flame and turbulent premixed flames. To realize DNS of IC engine with realistic kinetic mechanism, several DNS of elemental combustion process in IC engines has been conducted.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2005.04a
• /
• pp.172-175
• /
• 2005

The numerical simulation is performed for the acoustic emission and the wave propagation due to fiber breakage in single fiber composite plates by the finite element transient analysis. The acoustic emission and the following wave motions from a fiber breakage under a static loading is simulated to investigate the applicability of the explicit finite element method and the equivalent volume force model as a simulation tool of wave propagation and a modeling technique of an acoustic emission. For such a simple case of the damage event under static loading, various parameters affecting the wave motion are investigated for reliable simulations of the impact damage event. The high velocity and the small wave length of the acoustic emission require a refined analysis with dense distribution of the finite element and a small time step. In order to fulfill the requirement for capturing the exact wave propagation and to cover the 3-D simulation, we utilize the parallel FE transient analysis code and the parallel computing technology.