• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.10
• /
• pp.1368-1379
• /
• 2000

A numerical study on swirling pulverized coal combustion in an axisymmetric enclosure is carried out by applying the 2-phase weighted sum of gray gases model (WSGGM) approach with the discrete ordinate method (DOM) to model the radiative heat transfer equation. In the radiative transfer equation, the same polynomial equation and coefficients for weighting factors as those for gas are adopted for the coal/char particles as a function of partial pressure and particle temperature. The Eulerian balance equations for mass, momentum, energy, and species mass fractions are adopted with the standard and RNG k-${\varepsilon}$ turbulence model, whereas the Lagrangian approach is used for the particulate phase. The eddy-dissipation model is employed for the reaction rate for gaseous mixture, and the single-step and two-step first-order reaction model for the devolatilization process for coal. Special attention is given to establish the thermal boundary conditions on radiative transfer equation By comparing the numerical results with experimental ones, the radiation model used here is confirmed and found to provide an alternative for simulating the radiative transfer.

• International Journal of Fluid Machinery and Systems
• /
• v.10 no.3
• /
• pp.254-263
• /
• 2017

The cavitation erosion remains an industrial issue for many applications. This paper deals with the cavitation intensity, which can be described as the fluid mechanical loading leading to cavitation damage. The estimation of this quantity is a challenging problem both in terms of modeling the cavitating flow and predicting the erosion due to cavitation. For this purpose, a numerical methodology was proposed to estimate cavitation intensity from 3D unsteady cavitating flow simulations. CFD calculations were carried out using Code_Saturne, which enables U-RANS equations resolution for a homogeneous fluid mixture using the Merkle's model, coupled to a $k-{\varepsilon}$ turbulence model with the Reboud's correction. A post-process cavitation intensity prediction model was developed based on pressure and void fraction derivatives. This model is applied on a flow around a hydrofoil using different physical (inlet velocities) and numerical (meshes and time steps) parameters. The article presents the cavitation intensity model as well as the comparison of this model with experimental results. The numerical predictions of cavitation damage are in good agreement with experimental results obtained by pitting test.

• Journal of Ship and Ocean Technology
• /
• v.7 no.4
• /
• pp.16-31
• /
• 2003

The finite volume based multi-block RANS code, WAVIS developed at KRISO, is used to simulate the turbulent flows around a submarine with the realizable $\textsc{k}-\varepsilon$ turbulence model. RANS methods are verified and validated at the level of validation uncertainty 1.54% of the stagnation pressure coefficient for the solution of the turbulent flows around SUBOFF submarine model without appendages. Another SUBOFF configuration, axisymmetric body with four identical stem appendages, is also computed and validated with the experimental data of the nominal wake and hydrodynamic coefficients. The hydrodynamic forces and moments for SUBOFF model and a practical submarine are predicted at several drift and pitch angles. The computed results are in extremely good agreement with experimental data. Furthermore, it is noteworthy that all the computations at the present study were carried out in a PC and the CPU time required for 2.8 million grids was about 20 hours to get fully converged solution. The current study shows that CFD can be a very useful and cost effective tool for the prediction of the hydrodynamic performance of a submarine in the basic design stage.

• Macromolecular Research
• /
• v.16 no.2
• /
• pp.139-148
• /
• 2008

Biodegradable and elastic poly(L-lactide-co-$\varepsilon$-caprolactone) (PLCL) was electrospun to prepare nanofibers, and N-isopropylacrylamide (NIPAAm) was then grafted onto their surfaces under aqueous conditions using $^{60}Co-{\gamma}$ irradiation. The graft yield increased with increasing irradiation dose from 5 to 10 kGy and the nanofibers showed a greater graft yield compared with the firms. SEM confirmed that the PLCL nanofibers maintained an interconnected pore structure after grafting with NIPAAm. However, overdoses of irradiation led to the excessive formation of homopolymer gels on the surface of thc PLCL nanofibers. The equilibrium swelling and deswelling ratio of the PNIPAAm-g-PLCL nanofibers (prepared with 10 kGy) was the highest among the samples, which was consistent with the graft yield results. The phase-separation characteristics of PNIPAAm in aqueous conditions conferred a unique temperature-responsive swelling behavior of PNIPAAm-g-PLCL nanofibers, showing the ability to absorb a large amount of water at < $32^{\circ}C$, and abrupt collapse when the temperature was increased to $40^{\circ}C$. In accordance with the temperature-dependent changes in swelling behavior, the release rate of indomethacin and FITC-BSA loaded in PNIPAAm-g-PLCL nanofibers by a diffusion-mediated process was regulated by the change in temperature. Both model drugs demonstrated greater release rate at $40^{\circ}C$ relative to that at $25^{\circ}C$. This approach of the temperature-controlled release of drugs from PNIPAAm-g-PLCL nanofibers using gamma-ray irradiation may be used to design drugs and protein delivery carriers in various biomedical applications.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.1
• /
• pp.1-10
• /
• 1997

Viscous flow around actual ship is calculated by an use of RANS equations. The propriety of this computing method, usefulness to hull form design and the scale effect which is the effect of viscous flow depending on the scale of ship model are investigated. Reynolds stress is modelled by using k-${\varepsilon}$ turbulence model and the law of wall is applied near the body. Body fitted coordinates are introduced for the treatment of the arbitrary 3-dimensional shape of the ship hull form. The transformed equations in the computational domain are numerically solved by an employment of FVM. In the calculation of pressure, SIMPLE method is adopted and the solution of the discretized equation is obtained by the line-by-line method with the use of TDMA The calculations of two ships, 4410 TEU container carrier and 50,000 DWT class bulk carrier, are performed at model and actual ship scale. The results are compared and discussed with the model test results which are viscous resistance, nominal wake distribution at propeller plane and limiting streamline on the hull surface. They describe the effect of stem form and the scale effect very well. In particular, the calculated nominal wake distribution and limiting streamline are agreed qualitatively with the experiments and the viscous resistance values are estimated within ${\pm}5%$ difference from the resistance tests.

• International Journal of Air-Conditioning and Refrigeration
• /
• v.10 no.3
• /
• pp.117-128
• /
• 2002

This paper describes the smoke filling process of a fire field model based on a self-deve-loped 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-$\varepsilon$ turbulence model with buoyancy production term. Also it solves the radiation equation using the discrete ordinates method. Compressibility is assumed and the perfect gas law is used. Comparison of the calculated upper-layer average tempera-ture and smoke layer clear height with the zone models has shown reasonable agreement. The zone models used are the CFAST and the NBTC one-room. For atrium fires with ceiling glass the ceiling heat flux by solar heat causes a high smoke temperature near the ceiling. However, it has no effect on the smoke movement such as the smoke layer clear heights that are important in fire safety. In conclusion, the smoke layer clear heights that are important in evacuation activity except the early of a fire were not as sensitive as the smoke layer tem-perature to the nature of ceiling heat flux condition. Thus, a fire sensor in atrium with ceiling glass has to consider these phenomena.

• Journal of Biomedical Engineering Research
• /
• v.13 no.2
• /
• pp.87-96
• /
• 1992

In thIns paper, a numerical simulation of steady laminar and turbulent flow in a two dimensional model for the total artificial heart is'presented. A trlleaflet polyurethane valve was simulated at the outflow orifice while the Inflow orifice had a trileaflet or a flap valve. The finite analytic numerical method was employed to obtain solutions to the governing equations in the Cartesian coordinates. The closure for turbulence model was achieved by employing the k-$\varepsilon$-E model. The SIMPLER algo rithm was used to solve the problem in primitive variables. The numerical solutions of the slulated model show that regions of relative stasis and trapped vortices were smaller within the ventricular chamber with the flap valve at the Inflow orifice than that with the trileaflet valve. The predicted Reynolds stresses distal to the inflow valve within the ventricular chamber were also found to be smaller wlth the flap valve than with the trlleaflet valve. These resu1ts also suggest a correlation be- tween high turbulent stresses and the presence of thrombus In the vicinity of the valves in the total artificial hearts. The computed velocity vectors and trubulent stresses were comparable with previ ously reported in vitro measurements in artificial heart chambers. Analysis of the numerical solo talons suggests that geometries similar to the flap valve(or a tilting disc valve) results in a better flow dynamics within the total artificial heart chamber compared to a trileaflet valve.

• Geomechanics and Engineering
• /
• v.10 no.4
• /
• pp.387-403
• /
• 2016

In this paper, an integrated model for the wave (current)-induced seabed response is presented. The present model consists of two parts: hydrodynamic model for wave-current interactions and poro-elastic seabed model for pore accumulations. In the wave-current model, based on the fifth-order wave theory, ocean waves were generated by adding a source function into the mass conservation equation. Then, currents were simulated through imposing a steady inlet velocity on one domain and pressure outlet on the other side. In addition, both of the Reynolds-Averaged Navier-Stokers (RANS) Equations and $k-{\varepsilon}$ turbulence model would be applied in the fluid field. Once the wave pressures on the seabed calculated through the wave-current interaction model, it would be applied to be boundary conditions on the seabed model. In the seabed model, the poro-elastic theory would be imposed to simulate the seabed soil response. After comparing with the experimental data, the effect of currents on the seabed response would be examined by emphasize on the residual mechanisms of the pore pressure inside the soil. The build-up of the pore water pressure and the resulted liquefaction phenomenon will be fully investigated. A parametric study will also be conducted to examine the effects of waves and currents as well as soil properties on the pore pressure accumulation.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.2
• /
• pp.150-162
• /
• 2010

Assessment of hydrodynamic performance of a ship hull has been focused on a model ship rather than a full-scale ship. In order to design the propeller of a ship, model-scale wake is often extended to full-scale based upon an empirical method or designer's experience, since wake measurement data for a full-scale ship is very rare. Recently modern CFD tools made some success in reproducing wake field of a model ship, which implicates that there are some possibilities of the accurate prediction of full-scale wakes. In this paper firstly the evaluation of model-scale wake obtained by Fluent package was performed. It was found that CFD calculation with the Reynolds-stress model (RSM) provided much better agreement with wake measurement in the towing tank than with the realizable k-$\varepsilon$ model (RKE). In the next full-scale wake was calculated using the same package to find out the difference between model and full-scale wakes. Three hull forms of KLNG, KCS, KVLCC2 having measurement data open for the public, were chosen for the comparison of resistance, form factor, and propeller plane wake between model ships and full-scale ships.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.21 no.1
• /
• pp.72-83
• /
• 2017

A thermo-fluid dynamic analysis was performed through the numerical simulation of a missile canister. Calculation was made in a fixed analytical volume and fully evaporated water was used as a coolant. To analyze the interaction among the hot gas, coolant, and mixture flow, Realizable $k-{\varepsilon}$ turbulence and VOF(Volume Of Fluid) model were chosen and parametric study was performed with the change of coolant flow rate. It could be found that the pressure on the canister top nonlinearly increased with the increase of coolant flow rate. Temperature and coolant distribution were closely related to the flow behavior in canister. Temperature on the canister bottom indicated a decrease being proportional to coolant flow rate in early times but after a specific time, the temperature increased with the tendency being reversed. In addition, the early part of temperature showed a fluctuating phenomenon because of the overall circulatory flow of mixture gas.