• Journal of Environmental Science International
• /
• v.7 no.6
• /
• pp.745-751
• /
• 1998

The well-mixed room model has been traditionally used to predict the concentrations of contaminants in indoor environments. However, this is inappropriate because the flow fields in many indoor environments distribute contaminants non-uniformly, due to imperfect air mixing. Thus, some means used to describe an imperfectly mixed room are needed. The simplest model that accounts for imperfect air mixing is a two-zone model. Therefore, this study on development of computer program far the two-zone model is carried out to propose techniques of estimating the concentration of contaminants in the room. To do this, an important consideration is to divide a room into two-zone, i.e. the lower and upper zone assuming that the air and contaminants are well mixed within each zone. And between the zones the air recirculation is characterized through the air exchange parameter. By this basic assumption, the equations for the conservation of mass are derived for each zone. These equations are solved by using the computational technique. The language used to develope the program is a VISUAL BASIC. The value of air exchange coefficient($f_12$) is the most difficult to forecast when the concentrations of contaminants in an imperfectly mixed room are estimated by the two-zone model. But, as the value of $f_12$ increases, the air exchange between each zone increases. When the value of $f_12$ is approximately 15, the concentrations in both zone approach each other, and the entire room may be approximately treated as a single well-mixed room. Therefore, this study is available for designing of the ventilation to improve the air quality of indoor environments. Also, the two-zone model produces the theoretical base which may be extended to the theory for the multi-zone model, that will be contributed to estimate the air pollution in large enclosures, such as shopping malls, atria buildings, atria terminals, and covered sports stadia.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.39-46
• /
• 2011

A numerical study of a turbulent natural convection in an enclosure with the lattice Boltzmann method (LBM) is presented. The primary emphasis of the present study is placed on investigation of accuracy and numerical stability of the LBM for the turbulent natural convection flow. A HYBRID method in which the thermal equation is solved by the conventional Reynolds averaged Navier-Stokes equation method while the conservation of mass and momentum equations are resolved by the LBM is employed in the present study. The elliptic-relaxation model is employed for the turbulence model and the turbulent heat fluxes are treated by the algebraic flux model. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with the deferred correction way to ensure accuracy and stability of solutions. The present LBM is applied to the prediction of a turbulent natural convection in a rectangular cavity and the computed results are compared with the experimental data commonly used for the validation of turbulence models and those by the conventional finite-volume method. It is shown that the LBM with the present HYBRID thermal model predicts the mean velocity components and turbulent quantities which are as good as those by the conventional finite-volume method. It is also found that the accuracy and stability of the solution is significantly affected by the treatment of the convection term, especially near the wall.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.5
• /
• pp.353-363
• /
• 2017

The dynamics of a water-steam system in a once-through boiler was simulated based on the physics-based modeling approach, representing the system in response to large load change or scale disturbance simulations. The modeling considered the mass, energy conservation, and momentum equation in the water pipe and the focus was limited to the sub-critical pressure region. An evaporator tube modeling was validated against the reference data. A simplified boiler system consisting of economizer, evaporator, and superheater was constructed to match a 500 MW power boiler. The dynamic response of the system following a disturbance was discussed along with the quantitative response characteristics. The dynamic response of the boiler system was further evaluated by checking the case of an off-design point operation of the feedwater-to-fuel supply ratio. The results re-emphasized the significance of controlling the feedwater-to-fuel supply ratio and additional design requirements of the water-steam separator and spray attemperator.

• Journal of Bio-Environment Control
• /
• v.9 no.4
• /
• pp.212-222
• /
• 2000

This study was performed to investigate the performance of heat recovery device attached to exhaust gas flue connected to combustion chamber of greenhouse heating system. The experimental heat recovery system is mainly consisted of LPG combustion chamber and two heat recovery units; unit-A is attached directly to the exhaust gas flue, and unit-B is connected with unit-A. Heat recovery performance was evaluated by estimating total energy amounts by using enthalpy difference between two measurement points together with mass flow rate of gas and/or air passing through each heat recovery unit depending on 5 different flow rates controlled by voltage meter. The results of this experimental study, such as heat exchange behavior of supply air tubes and exhaust air passages crossing the tubes, pressure drop between inlet and outlet, heat recovery performance of exchange unit, etc., will be used as fundamental data for designing optimum heat recovery device to be used for fuel saving purpose by reducing heat loss amounts mostly wasted outside of greenhouse through flue.

• 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 the Korean Society for Aeronautical & Space Sciences
• /
• v.39 no.12
• /
• pp.1166-1173
• /
• 2011

Electro-Mechanical Ignition Safety Device(EMISD) for solid rocket motor is designed and manufactured. The EMISD utilizes a true rotary solenoid for arming mechanism and an electric squib(initiator) for generating ignition energy. In order to prove the ignition capability of the EMISD, 10-cc Closed Bomb Test(CBT) is performed, which measures the pressure built by high temperature and high pressure gas generated by operating EMISD. The pressure built in the free volume of 10-cc closed bomb and the opening time of the ignition gas outlet are calculated using one dimensional gas dynamic model which is composed of the ideal gas equation and mass-energy conservation equation. Comparing the test result with model prediction, it is realized that the pressure built in the free volume of closed bomb due to the firing of EMISD, has the efficiency ratio of about 34%.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2010.03a
• /
• pp.1409-1416
• /
• 2010

It is essential to predict a blasting-induced excavation damage zone (EDZ) beyond the proposed excavation line of a tunnel because the unwanted damage area requires extra support system for tunnel safety. Complicated blasting process which may hinder a proper characterization of the damage zone can be effectively represented by two loading mechanisms. The one is a dynamic impulsive load generating stress waves outwards immediately after detonation. The other is a gas pressure that remains for a relatively long time. Since the gas pressure reopens up the arrested cracks and continues to extend some cracks, it contributes to the final formation of EDZ induced by blasting. This paper presents the simple method to evaluate EDZ induced by gas pressure during blasting in rock. The EDZ is characterized by analyzing crack propagation from the blasthole. To do this, a model of the blasthole with a number of radial cracks of equal length in an infinite elastic plane is considered. In this model, the crack propagation is simulated by using three conditions, the crack propagation criterion, the mass conservation of the gas, and the adiabatic condition. As a result, the stress intensity factor of the crack generally decreases as crack propagates from the blasthole so that the length of the crack is determined. In addition, the effect of rock properties, initial number of cracks, and the adiabatic exponent are investigated.

• Journal of the Korean Geotechnical Society
• /
• v.29 no.6
• /
• pp.77-86
• /
• 2013

Embankments on soft ground experience significant deformation during time-dependent consolidation settlement, as well as an initial undrained settlement. Since infinitesimal strain theory assumes no configuration change and minute strain during deformation, finite strain analysis is required for better prediction of geotechnical problems involving large strain and geometric change induced by imposed loadings. Updated Lagrangian formulation is developed for time-dependent consolidation combining both force equilibrium and mass conservation of fluid, and mechanical constitutive equation is written in Janumann stress rate. Numerical convergence during Newton's iteration in large deformation analysis is improved by Nagtegaal's approach of considering the effect of rotation in mechanical constitutive relationship. Numerical simulations are conducted to discuss numerical reliability and applicability of developed numerical code: deformation of cantilever beam, two-dimensional consolidation. The numerical results show that developed formulation can efficiently describe large deformation problems. Proposed formulation is expected to facilitate the upgrading of a numerical code based on infinitesimal strain theory to that based on finite strain analysis.

• Journal of the Korea Concrete Institute
• /
• v.20 no.3
• /
• pp.317-323
• /
• 2008

In the current study, A model for the flow analysis of fresh and highly flowable concrete has been developed using a particle method, the moving particle semi-implicit (MPS) method. The phenomena on the flow of concrete has been considered as a visco-plastic flow problem, and the basic governing equation of concrete particle dynamics has been based on the Navier-Stokes equation in Lagrangian form and the conservation of mass. In order to formulate a visco-plastic flow constitutive law of fresh concrete, concrete is modeled as a highly viscous material in the state of non-flow and as a visco-plastic material in the state of flow after reaching the yield stress of fresh concrete. A flow test of fresh concrete in the L-box was simulated and the predicted flow was well matched with the experimental result. The developed method was well showed the flow motion of concrete particles because it was formulated to be based on the motion of visco-plastic fluid dynamics.

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

The conventional segregated finite element formulation produces a small and simple matrix at each step than in an integrated formulation. And the memory and cost requirements of computations are significantly reduced because the pressure equation for the mass conservation of the Navier-Stokes equations is constructed only once if the mesh is fixed. However, segregated finite element formulation solves Poisson equation of elliptic type so that it always needs a pressure boundary condition along a boundary even when physical information on pressure is not provided. On the other hand, the conventional integrated finite element formulation in which the governing equations are simultaneously treated has an advantage over a segregated formulation in the sense that it can give a more robust convergence behavior because all variables are implicitly combined. Further it needs a very small number of iterations to achieve convergence. However, the saddle-paint-type matrix (SPTM) in the integrated formulation is assembled and preconditioned every time step, so that it needs a large memory and computing time. Therefore, we newly proposed the P2PI semi-segregation formulation. In order to utilize the fact that the pressure equation is assembled and preconditioned only once in the segregated finite element formulation, a fixed symmetric SPTM has been obtained for the continuity constraint of the present semi-segregation finite element formulation. The momentum equation in the semi-segregation finite element formulation will be separated from the continuity equation so that the saddle-point-type matrix is assembled and preconditioned only once during the whole computation as long as the mesh does not change. For a comparison of the CPU time, accuracy and condition number between the two methods, they have been applied to the well-known benchmark problem. It is shown that the newly proposed semi-segregation finite element formulation performs better than the conventional integrated finite element formulation in terms of the computation time.