• Journal of Mechanical Science and Technology
• /
• v.19 no.3
• /
• pp.913-923
• /
• 2005

The dynamic mixed model (DMM) combined with a box filter of Zang et. al. (1993) has been generalized for passive scalar transport and applied to large eddy simulation of turbulent channel flows with Prandtl number up to 10. Results from a priori test showed that DMM is capable of predicting both subgrid-scale (SGS) scalar flux and dissipation rather accurately for the Prandtl numbers considered. This would suggest that the favorable feature of DMM, originally developed for the velocity field, works equally well for scalar transport problem. The validity of the DMM has also been tested a posteriori. The results of the large eddy simulation showed that DMM is superior to the dynamic Smagorinsky model in the prediction of scalar field and the model performance of DMM depends to a lesser degree on the ratio of test to grid filter widths, unlike in the a priori test.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.1B
• /
• pp.73-81
• /
• 2009

This paper presents a numerical investigation of the impact of the secondary currents on solute transport in open-channel flows. The RANS model with Reynolds stress model is used for flow modeling, and the GGDH(generalized gradient diffusion hypothesis) model is used to close the scalar transport equation. Using the developed model, the impact of secondary currents on solute transport in open channel flows over smooth-rough strip is investigated. Through numerical experiments, the secondary currents are found to affect the solute spreading, leading a movement of the position of the peak concentration and a skewed distribution of solute concentration. Due to the lateral flow of secondary currents near the free surface, the concentration at the rough strip is found to be larger than that at the smooth strip bed. The solute at the rough strip is more rapidly transported than smooth bed. A magnitude analysis of the solute transport rate in scalar transport equation is also carried out to investigate the effect of secondary currents and scalar flux on the concentration distribution.

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

In this study, turbulent flows in a planar combustor which has a square rib-type flame holder are numerically investigated by Large Eddy Simulation(LES). Firstly, the flow fields with or without jet injection downstream of the flame-holder are examined using uniform inlet velocity. Comparison of the present LES results with experimental one shows a good agreement. Secondly, to investigate mixing of oxidizer(air) and fuel injected behind the flame holder, the scalar-transport equation is introduced and solved. From the instantaneous flow and scalar fields, complex and intense mixing phenomena between fuel and jet are observed. It is shown that the ratio of jet to blocked air velocity is an important factor to determine the flow structure. Especially, when the ratio is large enough, the fuel jet penetrates the main vortices shed from the flame holder, resulting in significant changes in the flow and scalar fields.

• Nuclear Engineering and Technology
• /
• v.24 no.3
• /
• pp.319-328
• /
• 1992

The Galerkin formulation of the finite element method is applied to the integral law of the first-order form of the one-group neutron transport equation in one-dimensional spherical geometry. Piecewise linear or quadratic Lagrange polynomials are utilized in the integral law for the angular flux to establish a set of linear algebraic equations. Numerical analyses are performed for the scalar flux distribution in a heterogeneous sphere as well as for the criticality problem in a uniform sphere. For the criticality problems in the uniform sphere, the results of the finite element method, with the use of continuous finite elements in space and angle, are compared with the exact solutions. In the heterogeneous problem, the scalar flux distribution obtained by using discontinuous angular and spatical finite elements is in good agreement with that from the ANISN code calculation.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.55 no.6
• /
• pp.135-144
• /
• 2013

Ventilation efficiency has an important role in agricultural facilities such as greenhouse and livestock house to keep internally optimum environmental condition. Age-of-air concept allows to assess the ventilation efficiency of an agricultural facility according to estimating the ability of fresh air supply and contaminants emission using LMA and LMR. Most of these methods use a tracer gas method which has some limitations in experiment like dealing unstable and invisible gas. Therefore, the aim of this study was to develop a straightforward method to calculate age-of-air values with CFD simulation which has the advantage of saving computational time and resources and these method can solve the limitations in experiment using tracer gas method. The main idea of LMA computation is to solve the passive scalar transport equation with the assumption that the production of the time scalar throughout the room is uniform. In case of LMR calculation, the transport of the time scalar was reversed compulsively using UDF. The methodology to validate the results of this study was established by comparing with preceding research that had performed a computing LMA and LMR value by laboratory experiments and CFD simulations using tracer gas. As a result, the error was presented similarly level of results of preceding research. Some big errors could be caused by stagnated area and incongruity turbulence model. while the computational time was reduced to almost one fourth of that by preceding research.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.11
• /
• pp.1407-1416
• /
• 2004

A two-dimensional direct numerical simulation was performed to investigate the flame structure of CH$_4$/$N_2$-Air counterflow nonpremixed flame interacting with a single vortex. The detailed transport properties and a modified 16-step augmented reduced mechanism based on Miller and Bowman's detailed chemistry were adopted in this computation. The results showed that an initially flat stagnation plane, on which an axial velocity was zero, was deformed into a complex-shaped plane, and an initial stagnation point was moved far away from a vortex head when the counterflow field was perturbed by the vortex. It was noted that the movement of stagnation point could alter the species transport mechanism to the flame surface. It was also identified that the altered species transport mechanism affected the distributions of the mixture fraction and the scalar dissipation rate.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1972-1979
• /
• 2003

Transport of a scalar quantity, such as chemical concentration or temperature, is important in many engineering applications and environmental flows. Here we report on results obtained from the large eddy simulations of flow and concentration fields inside the tank performed using a spectral multi-domain technique. The computations were driven by specifying the impeller-induced flow at the blade tip radius (Yoon et al. $^{(1)}$). This study focused on the concentration development at different molecular diffusivities in a stirred tank operated under turbulent conditions. The main objective of the work presented here is to study the large-scale mixing structure at different molecular diffusivities in a stirred tank by using the large eddy simulation. The time sequence of concentration and flow fields shows the flow dependency of the concentration development. The presence of spatial inhomogenieties is detailed by observing the time variation of local concentration at different positions.

• 한국방재학회:학술대회논문집
• /
• 2008.02a
• /
• pp.233-236
• /
• 2008

To analyze the saline intrusion in the place, such as an estuary, the three-dimensional numerical model is developed. In this study, the advection terms of the governing equations are discretized by upwind scheme. By using an explicit scheme for the longitudinal direction and an implicit scheme for the vertical direction, the numerical model is free from the restriction of temporal step size caused by a relatively small grid ratio. The equation of state is used to consider the density, and the scalar transport equation for salinity is employed the third order TVD to scheme to prevent unphysical oscillation near discontinuity. In order to verify saline intrusion, the numerical model is conducted to compare the previous model in the lock exchange. The present model generally show a good agreement with the previous one.

• Nuclear Engineering and Technology
• /
• v.43 no.3
• /
• pp.243-256
• /
• 2011

This paper considers the concept of analog computing based on a cellular neural network (CNN) paradigm to simulate nuclear reactor dynamics using a time-dependent second order form of the neutron transport equation. Instead of solving nuclear reactor dynamic equations numerically, which is time-consuming and suffers from such weaknesses as vulnerability to transient phenomena, accumulation of round-off errors and floating-point overflows, use is made of a new method based on a cellular neural network. The state-of-the-art shows the CNN as being an alternative solution to the conventional numerical computation method. Indeed CNN is an analog computing paradigm that performs ultra-fast calculations and provides accurate results. In this study use is made of the CNN model to simulate the space-time response of scalar flux distribution in steady state and transient conditions. The CNN model also is used to simulate step perturbation in the core. The accuracy and capability of the CNN model are examined in 2D Cartesian geometry for two fixed source problems, a mini-BWR assembly, and a TWIGL Seed/Blanket problem. We also use the CNN model concurrently for a typical small PWR assembly to simulate the effect of temperature feedback, poisons, and control rods on the scalar flux distribution.

• Nuclear Engineering and Technology
• /
• v.32 no.6
• /
• pp.605-617
• /
• 2000

Nodal transport methods are studied for the solution of two dimensional discrete-ordinates, simplified even-parity transport equation(SEP) which is known to be an approximation to the true transport equation. The polynomial expansion nodal method(PEN) and the analytic function expansion nodal method(AFEN）which have been developed for the diffusion theory are used for the solution of the discrete-ordinates form of SEP equation. Our study shows that while the PEN method in diffusion theory can directly be converted without complication, the AFEN method requires a theoretical modification due to the nonhomogeneous property of the transport equation. The numerical results show that the proposed two methods work well with the SEP transport equation with higher accuracies compared with the conventional finite difference method.