• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.11 no.1
• /
• pp.73-80
• /
• 1999

This paper describes a numerical study of three-dimensional buoyant turbulent flow in a stairwell model with three convective differencing schemes, which include the upwind differencing scheme, the hybrid scheme and QUICK scheme. The Reynolds-averaged Navier-Stokes and energy equations are solved with a two-equation turbulence model. The Boussinesq approximation is used to model buoyancy terms in the governing equations. Three-dimensional predictions of the velocity and temperature fields are presented and are compared with experimental data. Three-dimensional simulations with each scheme have predicted the overall features of the flow fairly satisfactorily. A better agreement with experimental is achieved with QUICK scheme.

• Nuclear Engineering and Technology
• /
• v.51 no.5
• /
• pp.1272-1278
• /
• 2019

A lumped kinetic modeling platform is developed to investigate the coupled nuclear/thermo-fluid features of the closed natural circulation loop in a low power lead cooled fast reactor. This coolant material serves a reliable choice with noticeable thermo-physical safety characteristics in terms of natural convection. Boussienesq approximation is resorted to appropriately reduce the governing partial differential equations (PDEs) for the fluid flow into a set of ordinary differential equations (ODEs). As a main contributing step, the coolant circulation speed is accordingly correlated to the loop operational power and temperature levels. Further temporal analysis and control synthesis activities may thus be carried out within a more consistent state space framework. Nyquist stability criterion is thereafter employed to carry out a sensitivity analysis for the system stability at various power and heat sink temperature levels and results confirm a widely stable natural circulation loop.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.33 no.6
• /
• pp.257-264
• /
• 2021

This study established a numerical model capable of calculating the wave overtopping rate of coastal structures by nonlinear irregular waves using the FUNWAVE-TVD model, a fully nonlinear Boussinesq equation model. Here, a numerical model was established by coding the mean value approach equations of EurOtop (2018) and empirical formula by Goda (2009), and adding them as subroutines of the FUNWAVE-TVD model. The verification of the model was performed by numerically calculating the wave overtopping rate of nonlinear irregular waves on vertical wall structures and comparing them with the experimental results presented in EurOtop (2018). As a result of the verification, the numerical calculation result according to the EurOtop equation of this model was very well matched with the experimental result in all relative freeboard (R_c/H_mo) range under non-impulsive wave conditions, and the numerical calculation result of empirical formula was evaluated slightly smaller than the experimental result in R_c/H_mo < 0.8 and slightly larger than the experimental result in R_c/H_mo > 0.8. The results of this model were well represented in both the exponential curve and the power curve under impulsive wave conditions. Therefore, it was confirmed that this numerical model can simulate the wave overtopping rate caused by nonlinear irregular waves in an vertical wall structure.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.3
• /
• pp.342-351
• /
• 1999

A jet injected normally into a cross flow has been found to have the cross section of a horseshoe shape. It occurs by a twin vortex motion in the region downstream of the jet injection. Such a flow is inherently and highly three-dimensional and numerical calculations should play an important role. The three-dimensional momentum equations with buoyancy effect and energy equation are solved to obtain the velocity distributions, center-line trajectories, cross sectional shape and entrainment. The density difference is sufficiently small, so that the Boussinesq approximation is considered to be valid. The SIMPLE algorithm is applied in a staggered grid system of a calculational domain for the numerical method.

• Proceedings of the KSME Conference
• /
• 2001.06d
• /
• pp.522-527
• /
• 2001

Comprehensive numerical computations are made of side-heated squire cavity which is exposed to zero mean g-jitter. Numerical solutions are acquires to the governing two-dimensional Navier-Stokes equations for a Boussinesq fluid. When the system is exposed to pure sinusoidal g-jitter inclined to the vertical axis, in spite of zero mean gravity there exist non zero net flow fields [8]. The resonance phenomenon are observed in moderate Rayleigh number. And, by comprehensive numerical work, unlike[5], it is found that they are related with the overshoot phenomenon of the sudden gravity up problem.

• Journal of Korean Society for Atmospheric Environment
• /
• v.16 no.5
• /
• pp.511-520
• /
• 2000

The urban atmosphere is characterized by th difference in surface and atmospheric environment between urban and more natural area. To investigate th climatic effect of land use type and anthropogenic heat of urban on wind field, numerical simulations were carried out under typical summer synoptic condition. The wind model PNU_MCM(Pusan National University Mesoscale Circulation Model) is based on the three-dimensional Boussinesq equations, taking into account the hydrostatic assumption . Since lane-use differs over every subdivision on Pusan the surface energy budget model includes sub0grid parameterization scheme which can calculate the total heat flux over a grid surface composed of different surfaces. The simulated surface wind agrees well with the observed value, and average over 6 days which represent typical summer lan-sea breeze days, August 1998, i.e. negligible gradient winds and almost clear skies. Urbanization makes sea-breeze enhance at day and reduce land-breeze at night. The results show that contribution of land-use type is much larger than that of anthropogenic heat in Pusan.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.5
• /
• pp.1308-1318
• /
• 1995

The influence of varying rotation speed of both crystal and crucible was numerically investigated for the Czochralski silicon-crystal growth. Based on a simplified model assuming flatness of free surfrae, the Navier-Stokes Boussinesq equations were employed to identify the flow pattern, temperature distribution as well as the shape of the melt/crystal interface. The present results showed that the interface shape was relatively convex with respect to the melt at lower pulling rate and tended to be concave as the pulling rate increased. In particular, the experimentally observed gull-winged shape of the interface was qualitatively in agreement with the predicted shape. The rotation of crystal alone little affected the growth system. When the rotation speed of the crucible was increased, there occurred inversion of the interface shape from convex to concave pattern. At rapid rotation of the crucible, an interesting channel formation was predictied primarily due to the assumption of laminar flow.

• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.556-561
• /
• 2000

Comprehensive numerical computations are made of side-heated squire cavity which is exposed to time dependent external mechanical forcing. Numerical solutions are acquires to the governing two-dimensional Navier-Stokes equations for a Boussinesq fluid. Time dependent heat transfer characteristics of interior fluid are analyzed to illustrate resonance phenomenon. When system is exposed to pure sinusoidal mechanical forcing, the numerical results disclose that the basic mechanism of resonance of mechanical forcing is same as that of thermal forcing of Ref. [3, 9]. In comparatively small amplitude of mechanical forcing, thermal characteristics of the system are similar to basic system(${\varepsilon}=0$).

• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.271-278
• /
• 2020

The investigation into a full-scale 27 m high, by 6 m wide, thermosyphon loop. The simulation model is based on a one-dimensional axially-symmetrical control volume approach, where the loop is divided into a series of discreet control volumes. The three conservation equations, namely, mass, momentum and energy, were applied to these control volumes and solved with an explicit numerical method. The flow is assumed to be quasi-static, implying that the mass-flow rate changes over time. However, at any instant in time the mass-flow rate is constant around the loop. The boussinesq approximation was invoked, and a reasonable correlation between the experimental and theoretical results was obtained. Experimental results are presented and the flow regimes of the working fluid inside the loop identified. The results indicate that a series of such thermosyphon loops can be used as a cavity cooling system and that the one-dimensional theoretical model can predict the internal temperature and mass-flow rate of the thermosyphon loop.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.45-50
• /
• 2000

A numerical investigation is made of three-dimensional buoyant convection of a Boussinesq-fluid in a vertical cylinder. The top and bottom endwalls are thermally insulated. Flow is driven by the substantial azimuthal variations in thermal boundary conditions. Comprehensive numerical solutions to the Navier-Stokes equations are obtained. The representative Rayleigh number is large, thus, the overall flow pattern is of boundary layer-type. Three-dimensional (low characteristics are described. Specially, the global flow and the heat transfer features are delineated when the severity of azimuthal variation of sidewall temperature n, is intensified. Temperature and velocity fields on the meridional planes and the planes of constant height are presented. The global flow weakens as n becomes large. The pattern of the local Nusselt number on the surface of cylinder is similar regardless of n. The convective gain in heat transfer activities is reduced as n increases.