• Journal of Korean Society for Atmospheric Environment
• /
• v.13 no.5
• /
• pp.383-396
• /
• 1997

Numerical prediction of the pollutant dispersion over a two-dimensional hilly terrain is presented. The dispersion model used in the present work is based on the gradient diffusion theory and the finite-volume method on a non-orthogonal boundary-fitted grid system. The numerical model is validated by comparing the results with the available experimental data for the flat-floor dispersion within a turbulent boundary-layer. The numerical error analysis is performed based on the guideline of Kasibhatla et al.(1988) for the elevated-source dispersion in the flat-floor boundary layer having a power-law velocity and linear eddy-diffusivity profile. The influences of the two-dimensional hilly terrain on the dispersion from a continuously released source are numerically investigated by changing the emission locations and heights. It is found that the distributions of ground-level concentration are strongly influenced by the source location and the emission height. Hence, the terrain amplification factor is greatly enhanced when the pollutant source is located within a flow separation region. Dispersion from a source of short duration is also simulated and the duration time of the pollutant is compared at several downstream locations on a hilly terrain. The results of the numerical prediction are applied to the evaluation of environmental impacts due to the automobile exhausts at the seashore highway with a parallel mountain range.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.10
• /
• pp.1940-1954
• /
• 1992

Fine grid computations were attempted to analyze the turbulent flows in the near wall low Reynolds number region and the numerical analyses were incorporated by a finite-volume discretization with full find grid system and low Reynolds number k-.epsilon. model was employed in this region. For the improvement of low Reynolds number k-.epsilon. model, modification coefficient of eddy viscosity $f_{\mu}$ was derived as a function of turbulent Reynolds number $R_{+}$ and nondimensional length $y^{+}$ from the concept of two length scales of dissipation rate of turbulent kinetic energy. The modification coefficient $f_{\epsilon}$ in .epsilon. transport equation was also derived theoretically. In the turbulent kinetic energy equation, pressure diffusion term was added in order to consider low Reynolds number region effect. The main characteristics of this low Reynolds number k-.epsilon. model were founded as : (1) In high Reynolds number region, the present model has limiting behavior which approaches to the high Reynolds number model. (2) Present low Reynolds number k-.epsilon. model dose not need additional empirical constants for the transport equations of turbulent kinetic energy and dissipation of turbulent kinetic energy in order to consider wall effect. Present low Reynolds number turbulence model was tested in the pipe flow and obtained improved results in velocity profiles and Reynolds stress distributions compared with those from other k-.epsilon. models.s.s.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.33 no.4
• /
• pp.360-369
• /
• 1997

We calculated the residual current forced by buoyancy, wind stress, and tidal stress in the Chinhae Bay using a three-dimensional diagnostic model. The calculated current was also compared with the observation. The flow directs outward from the central area of the Bay in the upper layer, and also forms eddy-shape stucture in the upper and middle layers. The flow of bottom layer shows an opposite pattern compared to those of top and middle layers. The maximum speed was 6.05 em/see (September) and 3.49 cm/sec (November) in the upper layer, and 4.39 cm/sec on both month in the middle layer. The Kinetic energy of November (8.39xlO' W) was larger than that of September (1.24xlO 'W), mainly resulting from larger buoyancy effect in September. In spite of the roughness of the grid size(1 km) and wind date, the calculated flow shows eorrelation(r=0.71) with the observation. We expect that the correlation be increased by increased by adopting the fine grid and the variable coefficients of diffusion and viscosity.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.51-57
• /
• 2001

The eddy viscosity turbulence models were applied to predict the flows through a cascade, and the prediction performances of turbulence models were assessed by comparing with the experimental results for a controlled diffusion(CD) compressor blade. The original $\kappa-\omega$ turbulence model and $\kappa-\omega$ shear stress transport(SST) turbulence model were used as two-equation turbulence model which were enhanced for a low Reynolds number flow and the Baldwin-Lomax turbulence model was used as algebraic turbulence model. Farve averaged Wavier-Stokes equations in a two-dimensional, curvilinear coordinate system were solved by an implicit, cell-centered finite-volume computer code. The turbulence quantities are obtained by lagging when the men flow equations have been updated. The numerical analysis was made to the flows of CD compressor blade in a cascade at three different incidence angles (40. 43.4. 46 degrees). We found the reversion in the prediction performance of original $\kappa-\omega$ turbulence model and $\kappa-\omega$ SST turbulence model when the incidence angie increased. And the algebraic Baldwin-Lomax turbulence model showed inferiority to two-equation turbulence models.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.16 no.2
• /
• pp.92-102
• /
• 2004

In this paper, an analytical solution of calculating the excess temperature field due to a point heat source is presented in the presence of spring-neap modulation of convective alongshore flow. The basic form of the solution is identical to that given by Jung et al. (2003) but the convective term in the exponential kernel function is extended and a spring-neap variation in the horizontal eddy diffusivity is newly introduced. A set of calculations have been performed to examine the sensitivity of the heat build-up to the change in current fields and turbulent dispersion. Results indicate that the excess temperature field is confined within the tidal excursion distance, while the excess temperature field beyond the distance is mainly controlled by the horizontal diffusion. The heat build-up within the distance is considerably affected by the spring-neap variation in the horizontal eddy diffusivity; the relatively high excess temperature more than 1$^{\circ}C$ is extended further when the eddy diffusivity has spring-neap modulation.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.386-391
• /
• 2001

A finite element code based on P2P1 tetra element has been developed for the large eddy simulation (LES) of turbulent flows around a complex geometry. Fractional 4-step algorithm is employed to obtain time accurate solution since it is less expensive than the integrated formulation, in which the velocity and pressure fields are solved at the same time. Crank-Nicolson method is used for second order temporal discretization and Galerkin method is adopted for spatial discretization. For very high Reynolds number flows, which would require a formidable number of nodes to resolve the flow field, SUPG (Streamline Upwind Petrov-Galerkin) method is applied to the quadratic interpolation function for velocity variables, Noting that the calculation of intrinsic time scale is very complicated when using SUPG for quadratic tetra element of velocity variables, the present study uses a unique intrinsic time scale proposed by Codina et al. since it makes the present three-dimensional unstructured code much simpler in terms of implementing SUPG. In order to see the effect of numerical diffusion caused by using an upwind scheme (SUPG), those obtained from P2P1 Galerkin method and P2P1 Petrov-Galerkin approach are compared for the flow around a sphere at some Reynolds number. Smagorinsky model is adopted as subgrid scale models in the context of P2P1 finite element method. As a benchmark problem for code validation, turbulent flows around a sphere and a MIRA model have been studied at various Reynolds numbers.

• KEPCO Journal on Electric Power and Energy
• /
• v.5 no.3
• /
• pp.189-195
• /
• 2019

In this work, the flow and combustion characteristics using a 3-D numerical simulation was evaluated for a shale gas fueled combustor in a commercial class gas turbine. The Standard k-e turbulence model, 2 step methane oxidation mechanism, Finite rate/Eddy dissipation reaction model, DTRM radiation model were employed and validated well at the baseline condition (Natural Gas, Pilot Ratio 0.2). Based on the validated models, the combustion characteristics of shale gas was evaluated for three pilot ratios cases. It was found that NOx concentrations for all shale gas cases were less than the that for city gas, which imply that, at the selected PRs, the condition for combustion stability is satisfied. In addition, for higher PR, whereas the average temperatures at the exit are the same, the NOx increases. It means that diffusion combustion portion increases due to the higher PR.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.9
• /
• pp.1073-1084
• /
• 1999

The effects of secondary flow on the structure of a turbulent wake generated by a flat plate was investigated experimentally. The secondary flow was induced In a $90^{\circ}$ curved duct in which the flat plate wake generator was installed. The wake generator was installed in such a way that the wake velocity gradient exists in the span wise direction of the curved duct. Measurements were made in the plane containing the mean radius of curvature where pressure gradient and curvature effects were small compared with the secondary flow effect. All six components of the Reynolds stresses were measured in the curved duct. Turbulence intensities in the curved wake are higher than those in the straight wake due to an increase of the turbulent kinetic energy production by the secondary flow. In the inner wake region, shear stress and strain in the plane containing the velocity gradient of the wake show opposite signs with respect to each other, so that eddy viscosity Is negative in this region. This indicates that gradient-diffusion type turbulence models are not appropriate to simulate this type of flow.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.5
• /
• pp.962-978
• /
• 1989

A new model of the combustion rates of soot particle in turbulent flames has been suggested. This model applies the combustion rate of soot particles in laminar flames and uses local time-averaged quantities in order to consider the effect of the chemical reaction on the soot combustion in turbulent flames. The proposed rate equation has been tested for two propane-air turbulent round-jet diffusion flames and gives better predictions for the soot concentration field of two flames than the model previously used, especially in low temperature regions. A modified Monte carlo Method for analyzing radiative heat transfer of a flame also has been suggested and tested, which reveals good results.

• Asian Journal of Atmospheric Environment
• /
• v.10 no.1
• /
• pp.42-50
• /
• 2016

This study presents a performance evaluation of four different land surface models (LSM) available in Weather Forecast Research (WRF). The research site was located in Haean Basin in South Korea. The basin is very unique by its geomorphology and topography. For a better representation of the complex terrain in the mesoscale model were used a high resolution topography data with a spatial resolution of 30 meters. Additionally, land-use layer was corrected by ground mapping data-sets. The observation equipments used in the study were an ultrasonic anemometer with a gas analyzer, an automatic weather station and a tethered balloon sonde. The model simulation covers a four-day period during autumn. The result shows significant impact of LSM on meteorological simulation. The best agreement between observation and simulation was found in the case of WRF with Noah LSM (WRF-Noah). The WRF with Rapid Update Cycle LSM (WRF-RUC) has a very good agreement with temperature profiles due to successfully predicted fog which appeared during measurements and affected the radiation budget at the basin floor. The WRF with Pleim and Xiu LSM (WRF-PX) and WRF with Thermal Diffusion LSM (WRF-TD) performed insufficiently for simulation of heat fluxes. Both overestimated the sensible and underestimated the latent heat fluxes during the daytime.