• 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.

• Nuclear Engineering and Technology
• /
• v.26 no.4
• /
• pp.536-547
• /
• 1994

The state feedback optimal control techniques are used in designing the reactor control system. The mathematical plant model with the temperature feedback effects is established from the one delayed neutron group point kinetics equation and the singly lumped thermal-hydraulic balance equations, and is expressed in terms of state variables. The LQR (Linear Quadratic Regulator) control system is designed, being followed by the LQG (Linear Quadratic Gaussian) design to determine the optimal conditions of rod movement for the desired reactor power responses. And two different servo control schemes, the ordinary feedback system and the order increased regulating system, are proposed for the purpose of input tacking. The general control characteristics such as stability margins and output responses are discussed. Comparing each other, it is found that the order increased regulating system has far better control characteristics than the ordinary feedback system.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.6
• /
• pp.1385-1394
• /
• 1994

A three-dimensional numerical model of surface buoyant jets with variable density was established. The model uses fully nonlinear, time-dependent, three-dimensional, ${\sigma}$-transformed equations of motion and equation of heat transport. A semi-implicit numerical scheme in time has been adopted for computational efficiency. The model was applied for thermal jets discharging into a stagnant water and the simulated results were compared with a hydraulic experimental data set showing good agreement. Comparative studies of exchange coefficients and stability functions indicated that spatial variation of exchange coefficients should be considered and the existing stability functions should be modified to simulate surface buoyant jets accurately.

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

In this study, a numerical model was established and applied to simulate the steady-state groundwater and heat flow in an isotropic, heterogeneous, three dimensional aquifer system with uniform thermal properties and no change of state. This model was developed as an aid in screening large groundwater-flow systems as prospects for underground waste storage. Driving forces on the system are external hydrologic conditions of recharge from precipitation and fixed hydraulic head boundaries. Heat flux includes geothermal heat-flow, conduction to the land surface, advection from recharge, and advection to or from fixed-head boundaries. The model uses an iterative procedure that alternately solves the groundwater-flow and heat-flow equations, updating advective flux after solution of the groundwater-flow equation, and updating hydraulic conductivity after solution of the heat-flow equation. Dierect solution is used for each equation. Travel time is determined by particle tracking through the modeled space. Velocities within blocks are linear interpolations of velocities at block faces. Applying this model to the groundwater-flow system located in Jigyung-ri. Songla-myun, Youngil-gun. Kyungsangbuk-do, the groundwater-flow system including distribution of head, temperature and travel time and flow line, is analyzed.