• Nuclear Engineering and Technology
• /
• v.53 no.2
• /
• pp.455-465
• /
• 2021

The Molten Salt Reactor (MSR), one of the six advanced reactor types of the 4th generation nuclear energy systems, has many impressive features including economic advantages, inherent safety and nuclear non-proliferation. This paper introduces a system analysis code named TREND, which is developed and used for the steady and transient simulation of MSRs. The TREND code calculates the distributions of pressure, velocity and temperature of single-phase flows by solving the conservation equations of mass, momentum and energy, along with a fluid state equation. Heat structures coupled with the fluid dynamics model is sufficient to meet the demands of modeling MSR system-level thermal-hydraulics. The core power is based on the point reactor neutron kinetics model calculated by the typical Runge-Kutta method. An incremental PID controller is inserted to adjust the operation behaviors. The verification and validation of the TREND code have been carried out in two aspects: detailed code-to-code comparison with established thermal-hydraulic system codes such as RELAP5, and validation with the experimental data from MSRE and the CIET facility (the University of California, Berkeley's Compact Integral Effects Test facility).The results indicate that TREND can be used in analyzing the transient behaviors of MSRs and will be improved by validating with more experimental results with the support of SINAP.

• Advances in aircraft and spacecraft science
• /
• v.9 no.2
• /
• pp.95-102
• /
• 2022

The paper investigates the process of pulsation of a spherical cavity (bubble) in a liquid under the influence of a source of ultrasonic vibrations. The process of pulsation of a cavitation pocket in liquid is investigated. The Kirkwood-Bethe model was used to describe the motion. A numerical solution algorithm based on the Runge-Kutta-Felberg method of 4-5th order with adaptive selection of the integration step has been developed and implemented. It was revealed that if the initial bubble radius exceeds a certain value, then the bubble will perform several pulsations until the moment of collapse. The same applies to the case of exceeding the amplitude of ultrasonic vibrations of a certain value. The proposed algorithm makes it possible to fully describe the process of cavitation pulsations, to carry out comprehensive parametric studies and to evaluate the influence of various process parameters on the intensity of cavitation.

• Journal of Soil and Groundwater Environment
• /
• v.13 no.6
• /
• pp.50-61
• /
• 2008

In this study, new integrated particle tracking algorithm was developed to reduce the inherent problem of Eulerian- Lagrangian method, or adverse effect of particle tracking error on mass balance error. The new integrated particle tracking algorithm includes numerical method, semi-analytical method, and analytical method which consider both temporal and spatial changes of velocity field during time step. Detail of mathematical derivations is well illustrated and four examples are made to verify through the comparison of the new integrated particle tracking with analytical solution or Runge-Kutta method. Additionally, It was shown that the there is better superiority of the new integrated particle tracking algorithm over other existing particle tracking method such as Lu's method.

• Journal of Ocean Engineering and Technology
• /
• v.34 no.6
• /
• pp.406-418
• /
• 2020

In this study, a nonlinear wave simulation code is developed using a higher-order spectral (HOS) method. The HOS method is very efficient because it can determine the solution of the boundary value problem using fast Fourier transform (FFT) without matrix operation. Based on the HOS order, the vertical velocity of the free surface boundary was estimated and applied to the nonlinear free surface boundary condition. Time integration was carried out using the fourth order Runge-Kutta method, which is known to be stable for nonlinear free-surface problems. Numerical stability against the aliasing effect was guaranteed by using the zero-padding method. In addition to simulating the initial wave field distribution, a nonlinear adjusted region for wave generation and a damping region for wave absorption were introduced for wave generation simulation. To validate the developed simulation code, the adjusted simulation was carried out and its results were compared to the eighth order Stokes theory. Long-time simulations were carried out on the irregular wave field distribution, and nonlinear wave propagation characteristics were observed from the results of the simulations. Nonlinear adjusted and damping regions were introduced to implement a numerical wave tank that successfully generated nonlinear regular waves. According to the variation in the mean wave steepness, irregular wave simulations were carried out in the numerical wave tank. The simulation results indicated an increase in the nonlinear interaction between the wave components, which was numerically verified as the mean wave steepness. The results of this study demonstrate that the HOS method is an accurate and efficient method for predicting the nonlinear interaction between waves, which increases with wave steepness.

• Journal of the Korean Institute of Intelligent Systems
• /
• v.13 no.3
• /
• pp.259-259
• /
• 2003

Genetic algorithms using a Process of genetic evolution of an organism are appropriate for hard problems that have not been solved by any deterministic method. Up to now, the controller design method has been made with the frequency dependent specification but the design method with the time specification has gotten little progress. In this paper, we study the controller design to satisfy the performance of a plant using the generalized Manabe standard form. When dealing with a controller design in the case of two parameter configurations, there are some situations that neither a known pseudo inverse technique nor the inverse method can be applicable. In this case, we propose two methods of designing a controller by the gradient algorithm and the new pseudo inverse method so that the desired closed polynomials are either equalized to or approximated to the designed polynomial. Design methods of the proposed controller are implemented in Java.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.6 no.2
• /
• pp.12-20
• /
• 1998

Exter ballistics of a typical high-speed projectile is studied through a flow-visualization experiment and an unstructured grid Navier-Srokes computation. Experiment produced a schlieren photograph that adequately shows the characteristic features of this complex flow, namely two kinds of oblique cone shocks and turbulent wake developing into the downstream. A hybrid scheme of finite volume-element method is used to simulate the compressible Reynolds-Averaged Navier-Stok- es solution on unstructured grids. Osher's approximate Riemann solver is used to discretize the cinvection term. Higher-order spatial accuracy is obtained by MUSCL extension and van Albada ty- pe flux limiter is used to stabilize the numerical oscillation near the solution discontinuity. Accurate Gakerkin method is used to discretize the viscous term. Explict fourth-order Runge-Kutta method is used for the time-stepping, which simplifies the application of MUSCL extension. A two-layer k-$\varepsilon$ turbulence model is used to simulate the turbulent wakes accurately. Axisymmetric folw and two-dimensional flow with an angle of attack have been computed. Grid-dependency is also checked by carrying out the computation with doubled meshes. 2-D calculation shows that effect of angle of attack on the flow field is negligible. Axi-symmetric results of the computation agrees well with the flow visualization. Primary oblique shock is represented within 2-3 meshes in numerical results, and the varicose mode of the vortex shedding is clearly captured in the turbulent wake region.

• Progress in Superconductivity
• /
• v.14 no.1
• /
• pp.1-10
• /
• 2012

We fabricated a closed coils by using resistive-joint method and the joint resistance of the coils were estimated by field decay technique in liquid nitrogen. We used the Runge-kutta method for the numerical analysis to calculate the decay properties. The closed coil was wound by $(Bi,Pb)_2Sr_2Ca_2Cu_3O_x$/Ag tape. Both ends the tape were overlapped and soldered to each other. The current was induced in a closed coils by external magnetic flux density. Its decay characteristic was observed by means of measuring the magnetic flux density generated by induced current at the center of the closed coil with hall sensor. The joint resistance was calculated as the ratio of the inductance of the loop to the time constants. The joint resistances were evaluated as a function of critical current of loop, contact length, sweep time, and external magnetic flux density in a contact length of 7 cm. It was observed that joint resistance was dependent on contact length of a closed coil, but independent of critical current, sweep time, and external magnetic flux density. The joint resistance was measured to be higher for a standard four-probe method, compared with that for the field decay technique. This implies that noise of measurement in a standard four-probe method is larger than that of field decay technique. It was estimated that joint resistance was $8.0{\times}10^{-9}{\Omega}$ to $11.4{\times}10^{-9}{\Omega}$ for coils of contact length for 7 cm. It was found that 40Pb/60Sn solder are unsuitable for persistent mode.

• Structural Engineering and Mechanics
• /
• v.89 no.1
• /
• pp.33-46
• /
• 2024

The present research delves into the analysis of nonlinear free and forced vibrations of porous functionally graded (FG) shallow shells reinforced with oblique stiffeners, which are embedded in a nonlinear elastic foundation (NEF) subjected to external excitation. Two distinct types of PFG shallow shells, characterized by even and uneven porosity distribution along the thickness direction, are considered in the research. In order to model the stiffeners, Lekhnitskii's smeared stiffeners technique is implemented. With the stress function and first-order shear deformation theory (FSDT), the nonlinear model of the oblique stiffened shallow shells is established. The strain-displacement relationships for the system are derived via the FSDT and utilization of the von-Kármán's geometric assumptions. To discretize the nonlinear governing equations, the Galerkin method is employed. The model such developed allows analysis of the effects of the stiffeners with various angles as desired, in addition to the quantitative investigation on the influence of the surrounding nonlinear elastic foundations. To numerically solve the problem of vibrations, the 4th-order P-T method is used, as this method, known for its enhanced accuracy and reliability, proves to be an effective choice. The validation of the present research findings includes a comprehensive comparison with outcomes documented in existing literature. Additionally, a comparative analysis of the numerical results against those obtained using the 4th Runge-Kutta method is performed. The impact of stiffeners with varying angles and material parameters on the vibration characteristics of the present system is also explored. The researchers and engineers working in this field may use the results of this study as benchmarks in their design and research for the considered shell systems.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.69-69
• /
• 2011

본 논문에서는 비정상 상태의 비압축성 유동장을 해석하기 위하여 물체맞춤격자방법이 아닌 가상경계법을 사용하였다. 가상경계법은 구조격자를 사용하여 구조물 경계면에서 Momentum Forceing을 사용하여 가상의 경계를 만들어 유동장을 해석하는 방법이다. Navier-Stoke 방정식의 수치 이산화 방법으로 Kim et al(1985)이 사용한 Fractional Step Method(FSM)을 사용하였다. 시간에 대하여 semi-implicit FSM를 사용하였고, 확산항에 대해서는 2차 정확도의 Crank-Nicolson Method를 대류항은 3차 정확도의 Runge-Kutta Method를 사용 하였다. 본 연구에서는 가상경계법을 이용한 유동장 해석이 교량 단면에 대하여 수치해석이 가능한지 검토하였다. 가상경계법은 현재 많은 연구가 유선형의 구조물에 대하여 수행되어 오고 있다. 교량 단면과 같은 각 진 구조물에 대한 검토는 아직 미비한 실정이다. 가상경계법에서 다루고 있는 구조물 경계면에서의 Momentum Forcing 방법이 유선형의 구조물에 맞추어 연구가 진행되었기 때문이다. 먼저 본 연구의 프로그램을 검증하기 위하여 원형 실린더에 대하여 가상경계법을 적용한 결과 Re 수 200에서 Strouhal Number, 양력계수, 항력계수를 이전 연구 결과와 비교하였다. Williamson(1988)과 Zhang(1995)의 연구결과와 유사한 결과를 얻을 수 있다. 그리고 교량의 단면과 같은 각진 구조물(Bluff Body)에 대하여 가상경계법 적용하였다. 본 논문의 연구에서 평가 대상으로 하고 있는 2차원 교량 단면에 대하여 유동장 해석을 하였다. 본 논문에서 정량적인 유체력과 유동장에 대한 비교 및 검토가 이루어지지 못했지만 압력장과 유선의 형태가 이론적인 값을 벗어나지 않고 있는 것으로 확인 되었다. Re 수 2700에서 전산 해석을 수행하였으며, 교량 단면 주위의 압력계수와 박리현상 그리고 후류에서의 Vortex shedding 현상이 모두 적절한 분포가 나타나는 것을 확인할 수 있었다. 따라서 가상경계법을 이용하여 각진 구조물에 대한 주위 유동장해석에 대한 가능성을 확인하였으며, 풍동실험과의 결과비교를 통하여 가상경계법을 이용하여 교량 단면 주위의 유동장 해석 결과를 정량적으로 비교할 것이다.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.59-72
• /
• 2022

In this work the numerical results of the flow around a 5:1 rectangular cylinder at Reynolds numbers 3 000 and 40 000, zero angle of attack and smooth incoming flow condition are presented. Implicit Large Eddy Simulations (ILES) have been performed with a high-order accurate spatial scheme and an implicit high-order accurate time integration method. The spatial approximation is based on a discontinuous Galerkin (dG) method, while the time integration exploits a linearly-implicit Rosenbrock-type Runge-Kutta scheme. The aim of this work is to show the feasibility of high-fidelity flow simulations with a moderate number of DOFs and large time step sizes. Moreover, the effect of different parameters, i.e., dimension of the computational domain, mesh type, grid resolution, boundary conditions, time step size and polynomial approximation, on the results accuracy is investigated. Our best dG result at Re=3 000 perfectly agrees with a reference DNS obtained using Nek5000 and about 40 times more degrees of freedom. The Re=40 000 computations, which are strongly under-resolved, show a reasonable correspondence with the experimental data of Mannini et al. (2017) and the LES of Zhang and Xu (2020).