• 유체기계공업학회:학술대회논문집
• /
• 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.

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

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

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

• Journal of the Korean Society of Food Science and Nutrition
• /
• v.21 no.2
• /
• pp.213-218
• /
• 1992

The reaching time to the freezing point was to be fast in the order of 2% agar gel, 5% agar gel, 20% gelatin gel, pork, respectively. The freezing time and the passing time through the zone of the maximum ice crystal formation had linear relationship with the coolant temperature. The average diameter d$_{p}$ of ice crystal in a soybean protein gel and the moving of freezing front were represented an inverse proportion, and the moving velocity of freezing front was shown as 3.4$\times$10$^{-6}$ $\textrm{cm}^2$/sec from predicted theoretical formula. This value was very close to experimental results. The storage temperature did not give any influences for the growth of ice crystal in inside soybean protein gels during freezing conservation. The relationship between freezing condition and structure of freezing front was as follows : (moving velocity of freezing front) : (mass transfer rate of water at freezing point)$\times$(surface area of freezing front).

• Fire Science and Engineering
• /
• v.13 no.4
• /
• pp.20-29
• /
• 1999

This paper describes the smoke filling process 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 atrium space containing a ceiling heat flux. The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k- epsilon turbulence model with buoyancy term. Compressibility is assumed and the perfect gas law is used. The results of the calculated upper-layer average temperature and smoke layer interface height has shown reasonable agreement compared with the zone models. The zone models used are the CFAST developed at the Building and Fire Research Laboratory NIST U.S.A. and the NBTC one-room of FIRECALC developed at CSIRO, Australia. The smoke layer interface heights that are important in fire safety were not as sensitive as the smoke layer temperature to the nature of ceiling heat flux condition.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.15 no.4
• /
• pp.18-25
• /
• 2011

Solid propulsion systems have simple structures compared to other propulsion systems and are suitable for long-term storage. However the systems generally have limits on control of thrust levels. In this paper we suggest control algorithms for combustion chamber pressure of variable thrust solid propulsion systems using special nozzles such as pintle valve. For the pressure control within the chamber, we use a simple pressure change model by considering only mass conservation within the combustion chamber, design a classical algorithm and also a nonlinear controller using the feedback linearization technique. Also we derive the equation of the thrust for an under-expanded one-dimensional nozzle and then design a proportional-intergral controller after linearizing the thrust model for an operating point. Finally, we demonstrate the performance of the controller through a numerical simulation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.10
• /
• pp.1089-1095
• /
• 2011

Numerical simulation is performed for nucleate boiling on a micro-finned surface, which has been widely used to enhance heat transfer, by solving the equations governing the conservation of mass, momentum, and energy in the liquid and vapor phases. The bubble motion is determined by a sharp-interface level-set method, which is modified to include the effect of phase change and to treat the no-slip and contact-angle conditions, as well as the evaporative heat flux from the liquid microlayer on immersed solid surfaces such as micro fins and cavities. The numerical results for bubble formation, growth, and departure on a microstructured surface including fins and cavities show that the bubble behavior during nucleate boiling is significantly influenced by the fin-cavity arrangement and the fin-fin spacing.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.20 no.6
• /
• pp.746-751
• /
• 2014

In this paper, ship performance prediction solver was developed using open source computational fluid dynamics (CFD) libraries. The mass- and momentum-conservation equations and turbulent model with a wall function for the turbulent closer were considered. The volume fraction transport equation with a high-resolution interface capturing scheme were selected for free-surface simulation. The predicted wave pattern around KRISO container ship (KCS) using developed program showed good agreement against existing experimental data. For the revolution of a propeller in the propulsive test, general grid interface (GGI) library was used. The predicted propulsive performance showed 7 % difference against experimental data. Two-phase mixture model was developed to simulate a cavitation and applied to a propeller. The sheet cavitation on the propeller was predicted well. From results, the potential of the numerical tank developed by open source libraries was verified by applying it to KCS.

• Journal of Plant Biotechnology
• /
• v.40 no.3
• /
• pp.147-155
• /
• 2013

In this experiment, we report for the first time mass propagation by in vitro mircrotuberization of Codonopsis lanceolata. We first examined the effect of cytokinins on multiple shoot induction. 2.0 $mg{\cdot}L^{-1}$ of kinetin not only gave the highest rate of shoot induction (19.1%) but also the elongation of shoot (17.1 mm). Secondly, we investigated the effect of sugars on in vitro microtuberization from nodal segments. The diameter of tuberous roots was enlarged in the half-strength MS medium supplemented with 145.9 mM sucrose. Histological analysis revealed that the number of parenchymatous cell containing starch grains increased in the tuberous roots. In addition, unlike in non-tuberous root, vascular bundles were scattered inner cortex layer. Thirdly, in order to preserve and stimulate the germination, microtubers were stored at $4^{\circ}C$ refrigerator during 9 months and then transplanted to the artificial soils (vermiculrite : peatmoss = 1:1 v/v), resulting that the rates of survival and germination were 75% and 70%, respectively. These results indicated that mass propagation of C. lanceolata was achieved by in vitro microtuber formation, suggesting that this protocol might be applied for not only the propagation of elite clones but also conservation of C. lanceolata germplasm.