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

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.2
• /
• pp.453-462
• /
• 1990

The interaction of forced convection-conduction with thermal radiation in laminar boundary layer over a plate fin is studied numerically. The analysis is based on complete solution whereby the heat conduction equation for the fin is solved simultaneously with the conservation equations for mass, momentum and energy in the fluid boundary layer adjacent to the fin. The fluid is a gray medium and diffusion(Rosseland) approximation is used to describe the radiative heat flux in the energy equation. The resulting boundary value problem are convection-conduction parameter N$_{c}$ and radiation-conduction parameter m, Prandtl number Pr. Numerical results are presented for gases with the Prandtl numbers of 0.7 & 5 with values of N$_{c}$ and M ranging from 0 to 10 respectively. The object of this study is to provide the first results on forced convection-radiation interaction in boundary layer flow over a semi-infinite flay plate which can be used for comparisons with future studies that will consider a more accurate expression for the radiative heat flux. The agreement of the results from the complete solution presented by E. M. Sparrow and those from this paper for the special case of M=0 is good. The overall rate of heat transfer from the fin considering radiative effect is higher than that from the fin neglecting radiative effect. The local heat transfer coefficient with radiative effect is higher than that without radiative effect. In the direction from tip to base, those coefficients decrease at first, attain minimum, and then increase. The larger values of N$_{c}$ M, Pr give rise to larger fin temperature variations and the fin temperature without radiative effect is always higher than that with radiative effect.

• Atmosphere
• /
• v.25 no.3
• /
• pp.501-510
• /
• 2015

In this study, the characteristics of flows around building groups are investigated using a computational fluid dynamics (CFD) model. For this, building groups with different volumetric ratios in a fixed area are considered. As the volumetric ratio of the building group increases, the region affected by the building group is widened. However, the wind-speed reduced area rather decreases with the volumetric ratio near the ground bottom (z ${\lesssim}$ 0.7H, here, H is the height of the building group) and, above 0.7H, it increases. As the volumetric ratio decreases (that is, space between buildings was widened), the size of recirculation region decreases but flow recovery is delayed, resulting in the wider wind-speed reduced area. The increase in the volumetric ratio results in larger drag force on the flow above the roof level, consequently reducing wind speed above the roof level. However, above z ${\gtrsim}$ 1.7H, wind speed increases with the volumetric ratio for satisfying mass conservation, resultantly increasing turbulent kinetic energy there. Inside the building groups, wind speed decreased with the volumetric ratio and averaged wind speed is parameterized in terms of the volumetric ratio and background flow speed. The parameterization method is applied to producing averaged wind speed for 80 urban areas in 7 cities in Korea, showing relatively good performance.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2003.06a
• /
• pp.1643-1646
• /
• 2003

Drawing is a mechanical operation that attenuates thick material to an appropriate thickness for the next processing or end usage. When the input material has the form of a bundle or bundles made of very thin and long shaped wire or fibers, this attenuation operation is called "bundle drawing" or "drafting" Drafting is being used widely in manufacturing staple yarns. which is indispensable for the textile industry. However, the bundle processed by this operation undertake more or less defects in the evenness of linear density. Such irregularities cause many problems not only for the product quality but also for the efficiency of the next successive processes. Since long there have been many researches tying to find out factors affecting the irregularity of linear desity, to obtain optimal drafting conditions, to develop efficient measuring and analysis methods of linear density of bundle, etc., but there exists yet no fundamental equation describing the dynamic behavior of the flowing bundle during processing. In this research a mathematical model for the dynamic behavior of the bundle fluid is to be set up on the basis of general physical lows representing physical variables, i.e. linear density and velocity as the dynamic state of bundle. The conservation of mass and momentum balance was applied to the fluid field of bundle. while the movement of′ individual material was taken into account. The constitutive model relating the surface force and the deformation of bundle was introduced by considering a representative prodedure that stands for the bundle movement. Then a fundamental equations system could be simplified considering a steady state of the process. On the basis of the simplified model, the simulation was performed and the results could be confirmed by the experiments under various conditions.

• Korean Journal of Poultry Science
• /
• v.28 no.2
• /
• pp.143-153
• /
• 2001

Using bioinformatic tools for searching the massive genome databases, it is possible to Identify new genes in few minutes for initial discoveries based on evolutionary conservation, domain homology, and tissue expression patterns, followed by further verification and characterization using the bench-top works. The development of high-density two-dimensional arrays has allowed the analysis of the expression of thousands of genes simultaneously in the humans, mice, rats, yeast, and bacteria to elucidate the genes and pathways involved in physiological processes. In addition, rapid and automated protein identification is being achieved by searching protein and nucleotide sequence databases directly with data generated from mass spectrometry. Recently, analysis at the bio-chemical level such as biochemical screening and metabolic profiling (Biochemical genomics) has been introduced as an additional approach for categorical assignment of gene function. To make advantage of recent achievements in computational approaches for facilitated gene discoveries in the avian model, chicken expression sequence tags (ESTs) have been reported and deposited in the international databases. By searching EST databases, a chicken heparanase gene was identified and functionally confirmed by subsequent experiments. Using combination of sub-tractive hybridization assay and Genbank database searches, a chicken heme -binding protein family (cSOUL/HBP) was isolated in the retina and pineal gland of domestic chicken and verified by Northern blot analysis. Microarrays have identified several host genes whose expression levels are elevated following infection of chicken embryo fibroblasts (CEF) with Marek's disease virus (MDV). The ongoing process of chicken genome projects and new discoveries and breakthroughs in genomics and proteomics will no doubt reveal new and exciting information and advances in the avian research.