• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.29 no.10
• /
• pp.538-550
• /
• 2017

In a PWR (Pressurized Water Reactor), the appropriate heat removal from the surface of fuel rod bundle is important for ensuring thermal margins and safety. Although many CFD (Computational Fluid Dynamics) software have been used to predict complex flows inside fuel assemblies with mixing vanes, there is no domestic regulatory guideline for the comprehensive evaluation of CFD software. Therefore, from the nuclear regulatory perspective, it is necessary to perform the systematic assessment and prepare the domestic regulatory guideline for checking whether valid CFD software is used for nuclear safety problems. In this study, to provide systematic evaluation and guidance on the applicability of CFD software to the domestic nuclear safety area, the results of the sensitivity analysis for the effect of the discretization scheme accuracy for the convection terms and turbulence models, which are main factors that contribute to the uncertainty in the calculation of the nuclear safety problems, on the prediction performance for the turbulent flow distribution inside the fuel assembly with split-type mixing vanes were explained.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.7 s.250
• /
• pp.833-840
• /
• 2006

The paper presents a numerical solution to the problem of a hot rigid indenter sliding over a thermoelastic Winkler foundation with a thermal contact resistance at constant speed. It is shown analytically that no steady-state solution can exist for sufficiently high temperature or sufficiently small normal load or speed, regardless of the thermal contact resistance. However the steady state solution may exist in the same situation if the thermal contact resistance is considered. This means that the effect of the large values of temperature difference and small value of force or velocity which occur at no steady state can be lessened due to the thermal contact resistance. When there is no steady state, the predicted transient behavior involves regions of transient stationary contact interspersed with regions of separation regardless of the thermal contact resistance. Initially, the system typically exhibits a small number of relatively large contact and separation regions, but after the initial transient, the trailing edge of the contact area is only established and the leading edge loses contact, reducing the total extent of contact considerably. As time progresses, larger and larger numbers of small contact areas are established, unlit eventually the accuracy of the algorithm is limited by the discretization used.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.11 no.4
• /
• pp.91-102
• /
• 1991

Traditionally, a dynamic programming(DP) technique has been used to the multireservoir control system. The algorithm has inherent problem to increase computational requirements exponentially due to discretization of variables and expanding the dimension of the system. To solve this problem, this paper describes transforming the optimal control system into a quadratic programming(QP), using Laguerre polynomials(LP) and its properties. The objective function of the proposed QP is independent of time variable. The solution of the QP is obtained by nonlinear programming(NLP) using augmented Lagrangian multiplier method. The numerical experiment shows that the water level of reservoirs is higher than Lee's and the evaluated benefit value is about the same as other researcher's.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.05a
• /
• pp.330-333
• /
• 2005

The paper presents a numerical solution to the problem of a hot rigid indenter siding over a thermoelastic Winkler foundation with a thermal contact resistance at constant speed. It is shown analytically that no steady-state solution can exist for sufficiently high temperature or sufficiently small normal load or speed regardless of the thermal contact resistance. However, the steady state solution may exist in the same situation if the thermal contact resistance is considered. This means that the effect of the large values of temperature difference and small value of force or velocity which occur at no steady state can be lessened due to the thermal contact resistance. When there is no steady-state the predicted transient behavior involves regions of transient stationary contact interspersed with regions of separation regardless of the thermal contact resistance. Initially, the system typically exhibits a small number of relatively large contact and separation regions, but after the initial transient the trailing edge of the contact area is only established and the leading edge loses contact, reducing the total extent of contact considerably. As time progresses, larger and larger number of small contact areas are established, until eventually the accuracy of the algorithm is limited by the discretization used.

• Structural Engineering and Mechanics
• /
• v.40 no.3
• /
• pp.393-410
• /
• 2011

The eight-node 3D solid element is one of the most extensively used elements in computational mechanics. This is due to its simple shape and easy of discretization. However, due to the parasitic shear locking, it should not be used to simulate the behaviour of structural members in bending dominant conditions. Previous researches have indicated that the introduction of incompatible mode into the displacement field of the solid element could significantly reduce the shear locking phenomenon. In this study, an incompatible mode eight-node solid element, which considers both geometric and material nonlinearities, is developed for modelling of structural members at elevated temperatures. An algorithm is developed to extend the state determination procedure at ambient temperature to elevated temperatures overcoming initially converged stress locking when the external load is kept constant. Numerical studies show that this incompatible element is superior in terms of convergence, mesh insensitivity and reducing shear locking. It is also showed that the solid element model developed in this paper can be used to model structural behaviour at both ambient and elevated temperatures.

• Structural Engineering and Mechanics
• /
• v.61 no.2
• /
• pp.193-207
• /
• 2017

In this study, the vibration of an electrostatically actuated micro cantilever beam is analyzed in which a viscoelastic layer covers a portion of the micro beam length. This proposed model is considered as the main element of mass and pollutant micro sensors. The nonlinear motion equation is extracted by means of Hamilton principle, considering nonlinear shortening effect for Euler-Bernoulli beam. The non-linear effects of electrostatic excitation, geometry and inertia have been taken into account. The viscoelastic model is assumed as Kelvin-Voigt model. The motion equation is discretized by Galerkin approach. The linear free vibration mode shapes of non-uniform micro beam i.e. the linear mode shape of the system by considering the geometric and inertia effects of viscoelastic layer, have been employed as comparison function in the process of the motion equation discretization. The discretized equation of motion is solved by the use of multiple scale method of perturbation theory and the results are compared with the results of numerical Runge-Kutta approach. The frequency response variations for different lengths and thicknesses of the viscoelastic layer have been founded. The results indicate that if a constant volume of viscoelastic layer is to be deposited on the micro beam for mass or gas sensor applications, then a modified configuration may be found by using the analysis of this paper.

• Water for future
• /
• v.29 no.5
• /
• pp.173-184
• /
• 1996

In this study, a two-dimensional shallow-water equation was used to develop the mathematical model for computing water levels and flow distribution. In the discretization equations, based on the finite volume method (FVM), the third order Runge-Kutta method and the third order upwind scheme were introduced to handle the unsteady and vconvective terms in the governing equations. To determine the accuracy of the developed model, it was applied to the rectangular horizontal channel in a frictionless flow. The water depth and velocity obtained by the numerical model were found to agree closely with the exact solution. The model was also applied to the rectangular channel with both the symmetric and the non symmetric constriction. The velocity distribution of the flow and the propagation of the flood wave were simulated and the results well described the flow characteristics.

• Korea-Australia Rheology Journal
• /
• v.21 no.1
• /
• pp.47-58
• /
• 2009

We have investigated the transient behavior of 1D fully developed Poiseuille viscoelastic flow under finite pressure gradient described by the Oldroyd-B and Leonov constitutive equations. For analysis we employ a simple $2^{nd}$ order discretization scheme such as central difference for space and the Crank-Nicolson for time approximation. For the analysis of the Oldroyd-B model, we also apply the analytical solution, which is obtained again in this work in terms of elementary solution procedure simpler than the previous one (Waters and King, 1970). Both models demonstrate qualitatively similar solutions, but their eventual steady flowrate exhibits noticeable difference due to the absence or presence of shear thinning behavior. In the inertialess flow, the flowrate instantaneously attains a large value corresponding to the Newtonian creeping flow and then decreases to its steady value when the applied pressure gradient is low. However with finite liquid density the flow field shows severe fluctuation even accompanying reversals of flow directions. As the assigned pressure gradient increases, the flowrate achieves its steady value significantly higher than its value during oscillations after quite long period of time. We have also illustrated comparison between 1D and 2D results and possible mechanism of complex 2D flow rearrangement employing a previous solution of [mite element computation. In addition, we discuss some mathematical points regarding missing boundary conditions in 2D modeling due to the change of the type of differential equations when varying from inertialess to inertial flow.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.1 no.1
• /
• pp.128-144
• /
• 1998

The BEM for linear elastic stress analysis is applied to the highly rotating axisymmetric body problem which also involves the thermoelastic effects due to steady-state thermal conduction. The axisymmetric BEM formulation is briefly summarized and an alternative approach for transforming the volume integrals associated with such body force kernels into equivalent boundary integrals is described in a way of using the concept of inner product and vector identity. A discretization scheme for higher order BE is outlined for numerical treatment of the resulting boundary integral equations, and it is consequently illustrated by determining the stress distributions of the turbine rotor disk of the small turbojet engine(ADD 500) for which a FEM stress solution has been furnished by author.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.3
• /
• pp.607-621
• /
• 1988

The paper describes the incorporation of an algebraic stress model(ASM) of turbulence in to a semi-elliptic solution procedure for the prediction of turbulent flow in passage around a 180.deg. square sectioned bend. The numerical results are obtained from a finite-volume discretization with applications of QUICK scheme and full find grid system without PSL approximation. Results show that the better agreements in velocity profiles with experimental data than those from k, $\varepsilon$ equation model with wall function and PSL are obtained. Predictions of Reynolds stresses also show good agreements with the experimental data.