• Asian Journal of Atmospheric Environment
• /
• v.11 no.4
• /
• pp.270-282
• /
• 2017

To evaluate the deposition amount on a ground surface, mesoscale numerical models coupled with atmospheric chemistry are widely used for larger horizontal domains ranging from a few to several hundreds of kilometers; however, these models are rarely applied to high-resolution simulations. In this study, the performance of a dry and wet deposition model is investigated to estimate the amount of deposition via computational fluid dynamics (CFD) models with high grid resolution. Reynolds-averaged Navier-Stokes (RANS) simulations are implemented for a cone and a two-dimensional ridge to estimate the dry deposition rate, and a constant deposition velocity is used to obtain the dry deposition flux. The results show that the dry deposition rate of RANS generally corresponds to that observed in wind-tunnel experiments. For the wet deposition model, the transport equation of a new scalar concentration scavenged by rain droplets is developed and used instead of the scalar concentration scavenged by raindrops falling to the ground surface just below the scavenging point, which is normally used in mesoscale numerical models. A sensitivity analysis of the proposed wet deposition procedure is implemented. The result indicates the applicability of RANS for high-resolution grids considering the effect of terrains on the wet deposition.

• Nuclear Engineering and Technology
• /
• v.49 no.7
• /
• pp.1379-1387
• /
• 2017

This paper presents a comparison of homogenization techniques implemented in the APOLLO3 platform for transport core calculations: standard scalar flux weighting and new flux-moment homogenization, in different combinations with (or without) leakage models. Besides the historical B1-homogeneous model, a new B-heterogeneous one has indeed been implemented recently in the two/three-dimensional-transport solver using the method of characteristics. First analyses have been performed on a very simple Sodium Fast Reactor core with a regular hexagonal lattice. They show that using the heterogeneous leakage model in association with flux-moment homogenization strongly improves the prediction of $k_{eff}$ and void reactivity effects. These good results are confirmed when the application is done to the fissile assemblies of the more complex CFV (Low Void Effect) core of the ASTRID (Advanced Sodium Technological Reactor for Industrial Demonstration) project of sodium-cooled fast breeder reactor (Generation IV).

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

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.7 no.2
• /
• pp.193-203
• /
• 1983

The mixed convective heat transfer problems are characterized by the relatively significant contribution of buoyancy force to the transport processes of momentum and heat. Past analytical studies on this kind of problems have been carried out by employing either the conventional R-.epsilon. turbulence model which includes constant turbulent Prandtl number .sigma.$_{+}$ 1 or an extended R-.epsilon. turbulence model which takes account of the buoyancy effect in appropriate length scale equations. But in the latter case, the temperature variance .the+a.$^{2}$ over bar is approximated by a model under local equilibrium condition and the time scale ratio between velocity and temperature is assumed to be constant. These approximation is known to break down when the buoyancy effect is dominant. The present study is aimed at development of new computational turbulence closure level which can be applied to this rather complex turbulent process. The temperature variance is obtained directly by solving its dynamic transport equation and the time scale ratio which is variable in space is computed by a solution of a dynamic equation for the rate of scalar dissipation .epsilon.$_{\thetod}$ It was found that the computational results are in good agreement with available experimental data of wide range of unstable conditions.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1263-1267
• /
• 2007

A more complete two-dimensional vertical numerical model has been developed to describe the saltwater intrusion in an estuary. The model is based on the previous studies in order to obtain a better accuracy. The non-linear terms of the governing equations are analyzed and the $\sigma$-coordinate system is employed in the vertical direction with full transformation which is recently issued in several studies because numerical errors can be generated during the coordinate transformation of the diffusion term. The advection terms of the governing equations are discretized by an upwind scheme in second-order of accuracy. By employing 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. In previous researches, some terms induced from the transformation have been intentionally excluded since they are asked the complicate discretization of the numerical model. However, the lack of these terms introduces significant errors during the numerical simulation of scalar transport problems, such as saltwater intrusion and sediment transport in an estuary. The numerical accuracy attributable to the full transformation is verified by comparing results with a previous model in a simply sloped topography. The numerical model is applied to the Han River estuary. Very reasonable agreements for salinity intrusion are observed.

• Proceedings of the Korean Society of Coastal and Ocean Engineers Conference
• /
• 2003.08a
• /
• pp.208-215
• /
• 2003

Convective-diffusion equations appear in various disciplines such as hydrology, chemical engineering and oceanography dealing with the transport problem of scalar quantities. Since it is nonlinear, numerical methods are generally used to obtain its solution. Very limited number of analytical solutions are available usually in cases when the convective velocity is constant or has a simple functional form (for some collection of the solutions, see Noye, 1987). There is however a continuing need to develop analytical solutions because of its practical importance. Analytical solutions of the convection-diffusion equation are valuable not only for the better understanding on the transport process but the verification of numerical schemes. (omitted)

• Nuclear Engineering and Technology
• /
• v.54 no.9
• /
• pp.3478-3487
• /
• 2022

In this study, we incorporate an anisotropic scattering scheme involving spherical harmonics into the method of characteristics (MOC). The neutron transport solution in a light water reactor can be significantly improved because of the impact of an anisotropic scattering source with the MOC flat source approximation. Several problems are selected to verify the proposed scheme and investigate its effects and accuracy. The MOC anisotropic scattering source is based on the expansion of spherical harmonics with Legendre polynomial functions. The angular flux, scattering source, and cross section are expanded in terms of the surface spherical harmonics. Later, the polynomial is expanded to achieve the odd and even parity of the source components. Ultimately, the MOC angular and scalar fluxes are calculated from a combination of two sources. This paper presents various numerical examples that represent the hot and cold conditions of a reactor core with boron concentration, burnable absorbers, and control rod materials, with and without a reflector or baffle. Moreover, a small critical core problem is considered which involves significant neutron leakage at room temperature. We demonstrate that an anisotropic scattering source significantly improves solution accuracy for the small core high-leakage problem, as well as for practical large core analyses.

• 한국연소학회:학술대회논문집
• /
• 2014.11a
• /
• pp.183-184
• /
• 2014

The multi-environment probability density function model has been applied to simulate the turbulent stratified premixed flames. The direct quadrature method of moments (DQMOM) has been adopted to solve the transport PDF equation due to its computational efficiency and robustness. Computations are made for the non-swirling turbulent stratified premixed flames including SWB1, SWB5 and SWB9. The numerical results obtained in this study are precisely compared with experimental data in terms of axial velocity, unconditional means and conditional means for scalar field including temperature and species mass fraction.

• Journal of computational fluids engineering
• /
• v.10 no.3 s.30
• /
• pp.48-54
• /
• 2005

Three-dimensional incompressible flow past an open cavity in a channel is investigated using Detached Eddy Simulation(DES). The length to depth ratio of the cavity is 2 and the Reynolds number defined with the cavity depth is 3,360. The DES methods are based on the Menter's SST model. In the present work, two types of inflow conditions are used: one is RANS profile, the other is LES inflow from another Large Eddy Simulation(LES) of fully developed channel flow. The results are compared with experimental data and LES results in terms of the mean statistics, temporal physics and scalar transport phenomenon of the flow.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.239-243
• /
• 2005

Fully three-dimensional Large Eddy Simulation calculations of the flow past 2D cavity are conducted to study the purging of neutrally buoyant or dense miscible contaminants introduced instantaneously inside the cavity. The length to depth ratio(L/D) is 2 and Reynolds number based on the depth is 3,360. Fully developed turbulent inflow are fed at the inlet from precursor simulation of channel flow. Mean flow pattern and unsteady features are investigated based on the experimental data of Pereira and Sousa. From the study of mass exchange processes, it is found that the mechanism of removal of the contaminant is very different between the non-buoyant and buoyant cases. In the buoyant case, internal wave motion which interacts with a strong cavity vortex is dominant in the ejection mechanism of the contaminants.