• Advances in Energy Research
• /
• v.5 no.2
• /
• pp.91-105
• /
• 2017

Thermal hydraulic (TH) analysis of nuclear power reactors is utmost important. In this way, the numerical codes that preparing TH data in reactor core are essential. In this paper, a subchannel analysis of a Russian pressurized water reactor (WWER1000) core with enhanced numerical code is carried out. For this, in fluid domain, the mass, axial and lateral momentum and energy conservation equations for desired control volume are solved, numerically. In the solid domain, the cylindrical heat transfer equation for calculation of radial temperature profile in fuel, gap and clad with finite difference and finite element solvers are considered. The dependence of material properties to fuel burnup with Calza-Bini fuel-gap model is implemented. This model is coupled with Isotope Generation and Depletion Code (ORIGEN2.1). The possibility of central hole consideration in fuel pellet is another advantage of this work. In addition, subchannel to subchannel and subchannel to rod connection data in hexagonal fuel assembly geometry could be prepared, automatically. For a demonstration of code capability, the steady state TH analysis of a the WWER1000 core is compromised with Thermal-hydraulic analysis code (COBRA-EN). By thermal hydraulic parameters averaging Fuel Assembly-to-Fuel Assembly method, the one sixth (symmetry) of the Boushehr Nuclear Power Plant (BNPP) core with regular subchannels are modeled. Comparison between the results of the work and COBRA-EN demonstrates some advantages of the presented code. Using the code the thermal modeling of the fuel rods with considering the fission gas generation would be possible. In addition, this code is compatible with neutronic codes for coupling. This method is faster and more accurate for symmetrical simulation of the core with acceptable results.

• Journal of Korea Water Resources Association
• /
• v.48 no.8
• /
• pp.625-634
• /
• 2015

The diverse patterns of separation zone according to the marked bed discordance by dredging at confluence in addition to the confluence angle of tributary and discharge ratio between tributary and main channels have been analyzed by CCHE2D model simulation. The separation zone is defined by inside of zero velocity boundary at down-stream of confluence. The separation zone dose not formed at the $30^{\circ}$ of confluence angle of tributary. The size of separation zone increases as the discharge ratio and confluence angle increase in general. The separation zone decreases as the dredging depth increases which shows the relative momentum reduction compared by the flow volume increasing by dredging at confluence. The contraction factor with the variation of confluence angle and discharge ratio has been investigated and confirmed the corresponding conveyance decreasing results in backwater effect. The regression equation of shape factor with confluence angle and discharge and dredging depth ratios has been suggested.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.3 no.1
• /
• pp.1-8
• /
• 1999

In the present paper an attempt has been made to simulate the secondary injection-primary flow interaction in the conical rocket nozzle and to derive the performance of secondary injection thrust vector control(SITVC) system. Complex three-dimensional flowfield induced by the secondary injection is numerically analyzed by solving unsteady three-dimensional Euler equation with Beam and Warming's implicit approximate factorization method. Emphasized in the present study is the effect of secondary injection such as secondary mass flow rates and the momentum of secondary/primary nozzle flow mass rates upon the gross system performance parameters such as thrust ratio, specific impulse ratio and deflection angle. The results obtained in terms of system performance parameters show that lower secondary mass flow rate is advantageous for to reduce secondary specific impulse loss. It is further found that the nozzle with secondary jet injected downstream and interacting with fast primary flow is preferable for efficient and stable SITVC over the wide range of use with the penalty of side specific impulse loss.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.327-330
• /
• 2006

This paper represents the numerical analysis on effects of radius ratio in a concentric annulus with a rotating inner cylinder. The numerical model consisted of two cylinder which inner cylinder is rotating and outer cylinder is fix, and the axial direction is used the cyclic condition because of the length for axial direction is assumed infinite. The diameter of inner cylinder is assumed 86.8mm, the numerical parameters are angular velocity and radius ratio. Also, the whole walls of numerical model have no-slip and the working fluid is used water at $20^{\circ}C$. The numerical analysis is assumed the transient state to observe the flow variations by time and the 3-D cylindrical coordinate system. The calculation grid adopted a non-constant grid for dense arrangement near the wall side of cylinder, the standard $k-{\omega}$ high Reynolds number model to consider the effect of turbulence flow and wall, the fully implicit method for time term and the quick scheme for momentum equation. The numerical method is compared with the experimental results by Wereley and Lueptow, and the results are very good agreement. As the results, TVF isn't appeared when Re is small because of the initial flow instability is disappear by effect of the centrifugal force and viscosity. The vortex size is from 0.8 to 1.1 for TVF at various $\eta$, and the traveling distance for wavy vortex have the critical traveling distance for each case.

• Journal of Radiation Protection and Research
• /
• v.23 no.1
• /
• pp.1-6
• /
• 1998

A diffusion model of radioactive liquid effluents is developed and applied for YGN NPP's site, based on the Gaussian plume type model. Due to the complexity of oceanic diffusion characteristics of YGN site, a simple and reliable statistical model based on Reg. Guide 1.113 is developed. Also, a computer code package to calculate dilution factors as a function of plant operation conditions and pathway of radioactive materials. A liquid effluents diffusion model is developed by dividing the diffusion range into two categories, i. e, a near field mixing region and a far field mixing region. In the near field, the initial mixing is affected by a buoyance force, a high initial turbulence and momentum which is characterized by a plant operation condition and environmental conditions. The far field mixing is similar to gaseous effluents diffusion. So, beyond the near field region, wellknown Gaussian plume model was adopted. A different area averages of Gaussian plume equation was taken for each radioactive exposure pathway. As a result, we can get different dilution factors for different pathways. Results shows that present dilution factors used for YGN ODCM is too much overestimated compared with dilution factors calculated with the developed model.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.6 no.1
• /
• pp.87-94
• /
• 1986

Under natural condition, many rivers had shallow and gently curved shape in plane. A two dimensional mathematical model of the flow was a very attractive one. The flow characteristics in bended open channels were analyzed. The mathematical model based on the mass and the momentum equation of the two-dimensional unsteady flow was developed by introducing finite difference method and the double sweep algorithm. For the purpose of the verification of this model, the modeling results were applied to the L.F.M flume and the I.I.H.R flume. The results had a good agreement with the experimental data of the flumes. The results could be more close to the experimental data by controlling Chezy Coefficients in order to reduce the effect of friction around side wall, and be studied the importance of the convective term. The water surface profile, the direction and scale of depth average mean velocity and the path of the thread of maximum velocity in bended open channels could be computed.

• Fire Science and Engineering
• /
• v.15 no.1
• /
• pp.7-15
• /
• 2001

This paper describes the smoke filling process of a fire field model based on a self-developed SMEP (Smoke Movement Estimating Program) codo to the simulation of fire induced flows in the atrium space (SIVANS atrium at Japan) containing smoke radiation effect. The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k-$\varepsilon$ turbulence model with buoyancy term. Also it solves the radiation equation using the discrete ordinates method. The result of the calculated smoke temperature containing radiation effect has shown a better prediction than the result calculated by only convection effect in comparison with the experimental data. This seems to come from the radiation effect of $H_2$O and $CO_2$ gas under smoke productions. Thus, the consideration of the radiation effect under smoke in fire should be necessary in order to get more realistic result. Also the numerical results indicated that the smoke layer is developing at a rate of about 0.1 m/s. It would take about 450 seconds after starting the ultra fast fire of 560 kW that the smoke layer move down to 1.5m above the escape level.

• Journal of computational fluids engineering
• /
• v.18 no.3
• /
• pp.72-78
• /
• 2013

The turbulent flow behavior of air supply and exhaustion in the Shin-gum-ho subway station is analyzed for ordinary and emergency state. The depth of Shin-gum-ho station is 43.6m which consists of the island-type platform(8th floor in underground) and a two-story lobby (first & second floor in underground). An emergency stairway connects between the platform and the lobby. Ventilation operation mode for ordinary state is set up as a combination of air supply and exhaustion in the lobby and platform, while for emergency state it is set up as a full air supply in the lobby and a full exhaustion in the platform. The entire station is covered for simulation. The ventilation diffusers are modeled as 95 square shapes of $0.6m{\times}0.6m$ in the lobby and as 222 square shapes of $0.6m{\times}0.6m$ and 4 rectangular shapes of $1.2m{\times}0.8m$ in the platform. The total of 7.5million grids are generated and whole domain is divided to 22 blocks for MPI efficiency of calculation. Large eddy simulation(LES) is applied to solve the momentum equation and Smagorinsky model($C_s$=0.2) is used as SGS(subgrid scale) model. The time-averaged velocity fields are compared to experimental data and show a good agreement with it.

• Journal of the korean society of oceanography
• /
• v.37 no.3
• /
• pp.169-177
• /
• 2002

An eddy-resolving free-surface primitive-equation model with nonlinear terrain-following coordinates is established to study the exchange of water masses among the Asian marginal seas and their adjacent waters. A curvilinear coordinate system is used to generate the horizontal grid with a variable resolution for the regional oceans from $5^{\circ}$S to $45^{\circ}$N and $100^{\circ}$E to $155^{\circ}$E. The higher resolution region has about a 10 km by 10 km grid covering the complex geometry of the coastal marginal seas, while the lower resolution region has about a 30 km by 30 km grid covering the eastern Pacific. The model is initialized by the Levitus annual climitology and forced by the monthly mean air-sea fluxes of momentum, heat, and freshwater derived from the Comprehensive Ocean-Atmosphere Data Set. High-resolution and low-viscosity are identified as the key factors for a better representation of the exchange of waters through narrow straits and passages between the marginal seas and their adjacent waters. The dynamics of the loop currents and eddies in the South China Sea and Celebes Sea are examined in detail. It has found that the anticyclonic loop and detached eddies from the Kuroshio through the Luzon Strait play an important role in transporting warm and salty water into the South China Sea, while the cyclonic circulation of the Mindanao Current in the Celebes Sea plays a role in contributing cold water to the Indonesian throughflow. The deep undercurrent of the western Pacific is shown to provide fresher water to the South China Sea and Celebes Sea. These modeling results suggest that the exchange processes via the narrow straits and passages are of fundamental importance to the maintenance of water masses for the marginal sea region.

• Fire Science and Engineering
• /
• v.14 no.3
• /
• pp.6-12
• /
• 2000

This paper describes the smoke movement of a fire field model based on a self-developed SMEP(Smoke Movement Estimating Program) code to the simulation of fire induced flows in the two types of compartment space containing the radiation effect under smoke movement in room fire. The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k-$\varepsilon $ turbulence model with buoyancy term. Also it solves the radiation equation using the discrete ordinates method. The result of the calculated smoke temperature containing radiation effect has shown reasonable agreement compared with the experimental data. On the other hand, a difference of a lot was found between the temperature predicted by the SMEP with only convection effect and obtained by the experimental result. This seems to come from the radiation effect of $H_2$O and $CO_2$ gas under smoke productions. Thus, the consideration of the radiation effect under smoke in fire may be necessary in order to produce more realistic result.