• Particle and aerosol research
• /
• v.9 no.3
• /
• pp.163-171
• /
• 2013

Soot particles emitted from combustion processes are often coated by non-absorbing organic materials, which enhance the global warming effect of soot particles. It is of importance to study the condensation characteristics of soot particles experimentally and theoretically to reduce the uncertainty of the climate impact of soot particles. In this study, the condensational growth of soot particles in a tubular coater was modeled by a one-dimensional (1D) plug flow model and a two-dimensional (2D) laminar flow model. The effects of 2D heat and mass transports on the predicted particle growth were investigated. The temperature and coating material vapor concentration distributions in radial direction, which the 1D model could not accounted for, affected substantially the particle growth in the coater. Under the simulated conditions, the differences between the temperatures and vapor concentrations near the wall and at the tube center were large. The neglect of these variations by the 1D model resulted in a large error in modeling the mass transfer and aerosol dynamics occurring in the coater. The 1D model predicted the average temperature and vapor concentration quite accurately but overestimated the average diameter of the growing particles considerably. At the outermost grid, at which condensation begins earliest due to the lowest temperature and saturation vapor concentration, condensing vapor was exhausted rapidly because of the competition between condensations on the wall and on the particle surface, decreasing the growth rate. At the center of the tube, on the other hand, the growth rate was low due to high temperature and saturation vapor concentration. The effects of Brownian diffusion and thermophoresis were not high enough to transport the coating material vapor quickly from the tube center to the wall. The 1D model based on perfect radial mixing could not take into account this phenomenon, resulting in a much higher growth rate than what the 2D model predicted. The result of this study indicates that contrary to a previous report for a thermodenuder, 2D heat and mass transports must be taken into account to model accurately the condensational particle growth in a coater.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.687-690
• /
• 2002

Multimode boundary-layer transition on a NACA0012 airfoil is experimentally investigated under periodically passing wakes and the moderate level of free-stream turbulence. The periodic wakes are generated by rotating circular cylinders clockwise or counterclockwise around the airfoil. The free-stream turbulence is produced by a grid upstream of the rotating cylinder, and its intensity(Tu) at the leading edge of the airfoil is $0.5\;or\;3.5\;{\%}$. The Reynolds number ($Re_c$) based on chord length (C) of the alrfoil is $2.0{\times}10^5$, and Strouhal number ($St_c$) of the passing wake is about 0.7. Time- and phase-averaged streamwise mean velocities and turbulence fluctuations are measured with a single hot-wire probe, and especially, the corresponding wall skin friction is evaluated using a computational Preston tube method. The wake-passing orientation changes pressure distribution on the airfoil in a different manner irrespective of the free-stream turbulence. Regardless of free-stream turbulence level, turbulent patches for the receding wakes propagate more rapidly than those for the approaching wake because adverse pressure gradient becomes larger. The patch under the high free-stream turbulence ($Tu=3.5{\%}$) grows more greatly in laminar-like regions compared with that under the low background turbulence ($Tu=0.5{\%}$) in laminar regions. The former, however, does not greatly change the original turbulence level in the very near-wall region while the latter does it. At further downstream, the former interacts vigorously with high environmental turbulence inside the pre-existing transitional boundary layer and gradually lose his identification, whereas the latter keep growing in the laminar boundary layer. The calmed region is more clearly observed under the lower free-stream turbulence level and for the receding wakes. The calmed region delays the breakdown further downstream and stabilizes more the boundary layer.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.9 s.252
• /
• pp.921-928
• /
• 2006

Wake-induced boundary-layer transition on a NACA0012 airfoil with zero angle of attack is experimentally investigated in periodically passing wakes under the moderate level of free-stream turbulence. The periodic wakes are generated by rotating circular cylinders clockwise or counterclockwise around the airfoil. The free-stream turbulence is produced by a grid upstream of the rotating cylinder, and its intensities $(Tu_{\infty})$ at the leading edge of the airfoil are 0.5 and 3.5%, respectively. The Reynolds number (Rec) based on chord length (C) of the airfoil is $2.0{\times}10^5$, and Strouhal number (Stc) of the passing wake is about 1.4. Time- and phase-averaged streamwise mean velocities and turbulence fluctuations are measured with a single hot-wire probe, and especially, the corresponding wall skin friction is evaluated using a computational Preston tube method. The patch under the high free-stream turbulence $(Tu_{\infty}=3.5%)$ grows more greatly in laminar-like regions compared with that under the low turbulence $(Tu_{\infty}=0.5%)$ in laminar regions. The former, however, does not greatly change the turbulence level in very near-wall region while the latter does it. At further downstream, the former interacts vigorously with high environmental turbulence inside the pre-existing transitional boundary layer and gradually loses its identification, whereas the latter keeps growing in the laminar boundary layer. The calmed region is more clearly observed under the lower free-stream turbulence level and with the receding wakes.

• Ecology and Resilient Infrastructure
• /
• v.4 no.1
• /
• pp.54-62
• /
• 2017

This study investigated compound open channel flow using OpenFOAM Large Eddy Simulation (LES). Large eddy simulations were carried out by solving the filtered continuity and momentum equations numerically. One equation LES and non-uniform grid were applied to capture the anisotropic turbulence and secondary flow near the wall. The results of large eddy simulations of turbulent flow in a compound open channel with deep and shallow flood plain depths are presented. These LESs are validated with experimental data, resulting in a good agreement between measured and calculated data. The role of anisotropic turbulence in generating secondary currents is illustrated.

• International Journal of Aeronautical and Space Sciences
• /
• v.12 no.2
• /
• pp.149-155
• /
• 2011

The viscous flow in an axisymmetric nozzle was analyzed while accounting for the mesh sizes in both in the free stream and the boundary layer. The Navier-Stokes equations were resolved using the finite volume method in order to determine the supersonic flow parameters at the exit of the converging-diverging nozzle. The numerical technique in the aforementioned method uses the flux vector splitting of Van Leer. An adequate time stepping parameter, along with the Courant, Friedrich, Lewis coefficient and mesh size level, was selected to ensure numerical convergence. The boundary layer thickness significantly affected the viscous flow parameters at the exit of the nozzle. The best solution was obtained using a very fine grid, especially near the wall at which a strong variation of velocity, temperature and shear stress was observed. This study confirmed that the boundary layer thickness can be obtained only if the size of the mesh is lower than a certain value. The nozzles are used at the exit of the shock tube in order to obtain supersonic flows for various tests. They also used in propulsion to obtain the thrust necessary to the displacement of the vehicles.

• Journal of computational fluids engineering
• /
• v.14 no.1
• /
• pp.9-17
• /
• 2009

In this study we investigated the characteristics of fluid flow and heat transfer within a coolant passage in the cylinder head assembly of an internal combustion engine by using a commercial CFD code, CFX The complex coolant passage of the cylinder head assembly was modelled by suitable choice of a grid system and careful attention was paid in the construction of meshes near the walls where significant cooling occurs. To treat the simultaneous heating and cooling of the combustion walls we invented a methodology allowing a heat source within the solid wall and the convective cooling at the interface between the solid and the fluid. We managed to reproduce the experimental results by adjusting parameters appropriately. We have found that high temperature was concentrated at the surface of the cylinder jacket. It turned out that the effect of oil cooling from the piston head was unexpectedly significant. On the other hand the effect of cooling from the ambient air is almost negligible. The CFD method proposed in this study is believed to be useful in the early stage of the design of the engine-cooling system.

• Nuclear Engineering and Technology
• /
• v.36 no.3
• /
• pp.248-262
• /
• 2004

Turbulence models are separately assessed for a three dimensional thermal-hydraulic analysis of the integral reactor SMART. Seven models (mixing length, k-l, standard $k-{\epsilon},\;k-{\epsilon}-f{\mu},\;k-{\epsilon}-v2$, RRSM, and ERRSM) are investigated for flat plate channel flow, rotating channel flow, and square sectioned U-bend duct flow. The results of these models are compared to the DNS data and experiment data. The results are assessed in terms of many aspects such as economical efficiency, accuracy, theorization, and applicability. The standard $k-{\epsilon}$ model (high Reynolds model), the $k-{\epsilon}-v2$ model, and the ERRSM (low Reynolds models) are selected from the assessment results. The standard $k-{\epsilon}$ model using small grid numbers predicts the channel flow with higher accuracy in comparison with the other eddy viscosity models in the logarithmic layer. The elliptic-relaxation type models, $k-{\epsilon}-v2$, and ERRSM have the advantage of application to complex geometries and show good prediction for near wall flows.

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

• Wind and Structures
• /
• v.37 no.4
• /
• pp.255-274
• /
• 2023

The aim of this paper is to enhance the accuracy of predicting time-averaged external surface pressures on low-rise buildings by utilizing Computational Fluid Dynamics (CFD) simulations. To achieve this, benchmark studies of the Silsoe cube and the Texas Tech University (TTU) experimental building are employed for comparison with simulation results. The paper is structured into three main sections. In the initial part, an appropriate domain size is selected based on the precision of mean pressure coefficients on the windward face of the cube, utilizing Reynolds Averaged Navier-Stokes (RANS) turbulence models. Subsequently, recommendations regarding the optimal computational domain size for an isolated building are provided based on revised findings. Moving on to the second part, the Silsoe cube model is examined within a horizontally homogeneous computational domain using more accurate turbulence models, such as Large Eddy Simulation (LES) and hybrid RANS-LES models. For computational efficiency, transient simulation settings are employed, building upon previous studies by the authors at the Windstorm Impact, Science, and Engineering (WISE) Lab, Louisiana State University (LSU). An optimal meshing strategy is determined for LES based on a grid convergence study. Three hybrid RANS-LES cases are investigated to achieve desired enhancements in the distribution of mean pressure coefficients on the Silsoe cube. In the final part, a 1:10 scale model of the TTU building is studied, incorporating the insights gained from the second part. The generated flow characteristics, including vertical profiles of mean velocity, turbulence intensity, and velocity spectra (small and large eddies), exhibit good agreement with full-scale (TTU) measurements. The results indicate promising roof pressures achieved through the careful consideration of meshing strategy, time step, domain size, inflow turbulence, near-wall treatment, and turbulence models. Moreover, this paper demonstrates an improvement in mean roof pressures compared to other state-of-the-art studies, thus highlighting the significance of CFD simulations in building aerodynamics.

• Journal of the Korean earth science society
• /
• v.28 no.2
• /
• pp.202-213
• /
• 2007

The Odaesan earthquake $(M_L=4.8)$ occurred near Mt. Odae, Jinbu-Myon, Pyongchang-Gun, Kangwon Province on January 20, 2007. It has a shallow focal depth about 10 km. Its felt area covers most of the southern peninsula except some southern and western inland area. The maximum MM intensity was VI in the areas including Jinbu, Doam, Kangreung, Jumunjin, and Pyongchang. In these areas, there was a very strong shaking that caused several cracks on the walls of buildings and houses, slates falling off the roof, tiles being off the wall, things falling off the desk, and rock falling from the mountains. In order to get fault plane solutions, grid searches were performed by fitting distributions of P-wave first-motion polarities and SH/P amplitude ratios for each event. The results showed that the main shock represented right-lateral strike-slip sense and two aftershocks, reverse sense. It seems that the seismogenic fault may be the NNE-SSW trending Weoljeongsa fault near the epicenter based on the distribution of epicenters (foreshock, main shock, and aftershocks), damage area, and fault plane solution. The distribution of the epicenters indicates that the length of the subsurface rupture is estimated to be about 2 km.