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

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

• Communications of the Korean Mathematical Society
• /
• v.22 no.1
• /
• pp.145-155
• /
• 2007

In this paper, we show the convergence of approximate solutions to the convective porous media equation using methodology developed in [8]. First, we obtain the approximate transport equation for the given convective porous media equation. Then using the averaging lemma, we obtain the convergence.

• 한국방재학회:학술대회논문집
• /
• 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.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.

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

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

• Journal of Energy Engineering
• /
• v.16 no.1 s.49
• /
• pp.46-52
• /
• 2007

The piecewise constant angular approximation is developed to replace the conventional angular quadrature sets in the solution of the second-order, multi-dimensional $S_{N}$ neutron transport equations. The newly generated quadrature sets by this method substantially mitigate ray effects and can be used in the same manner as the conventional quadrature sets are used. The discrete-ordinates and the piecewise-constant approximations are applied to both the first-order Boltzmann and the second-order form of neutron transport equations in treating angular variables. The result is that the mitigation of ray effects is only achieved by the piecewise-constant method, in which new angular quadratures are generated by integrating angle variables over the specified region. In other sense, the newly generated angular quadratures turn out to decrease the contribution of mixed-derivative terms in the even-parity equation that is one of the second-order neutron transport equation. This result can be interpreted as the entire elimination or substantial mitigation of ray effect are possible in the simplified even-parity equation which has no mixed-derivative terms.

• Journal of Korea Water Resources Association
• /
• v.45 no.5
• /
• pp.473-480
• /
• 2012

A two-dimensional depth-integrated hydrodynamic and a depth-averaged passive scalar transport models were developed by using a Compact Finite Volume Method (CFVM) which can assure a higher order accuracy. A typical wave current interaction experimental data set was compared with the computed results by the proposed CFVM model, and resonable agreements were observed from the comparisons. One and two dimensional scalar advection tests were conducted, and very close agreements were observed with very little numerical diffusion. Finally, a turbulent mixing simulation was done in an open channel flow, and a reasonable similarity with LES data was observed.