• Asian Journal of Atmospheric Environment
• /
• 제11권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
• /
• 제49권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).

• 에너지공학
• /
• 제16권1호
• /
• pp.46-52
• /
• 2007

특정한 방향에 대해 방향중성자속(angular neutron flux)을 정의하는 방향차분 방정식(discrete-ordinates or $S_{N}$ equation)과 달리 방향변수를 구분된 방향영역에 대하여 적분한 값을 사용하고, 해당 방향영역 내에서 방향중성자속이 일정하다고 가정하는 영역상수법(piecewise-constant method)을 개발하였다. 기존 방향차분법과 본 연구에서 개발된 영역상수법을 1계 수송방정식(1'st-order Boltzmann transport equation)과 2계 우성 방정식(even-parity equation)에 적용하여 방향차분 방정식인 $S_{N}$ 방정식과 유사 방향차분방정식($S_{N}$-like equation)인 $PC_{N}$ 방정식을 유도하였다. 우성 방정식에 영역상수법을 적용한 경우 기존 방향차분법의 단점인 광첨두 현상(ray effect)이 현저히 감소함을 확인하였는데 이는 우성 방정식의 혼합 미분항의 기여도가 작아지기 때문인 것으로 판단된다. 이러한 이론은 우성 방정식에서 혼합 미분항이 제거된 단순우성 방정식(simplified even-parity equation)을 사용하는 경우 광첨두 현상이 완전 제거 또는 극단적으로 감소되었던 이전의 결과를 이론적으로 설명한다.

• 대한기계학회논문집
• /
• 제7권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.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2007년도 학술발표회 논문집
• /
• 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.

• 한국해안해양공학회:학술대회논문집
• /
• 한국해안해양공학회 2003년도 한국해안해양공학발표논문집
• /
• 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
• /
• 제54권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년도 제49회 KOSCO SYMPOSIUM 초록집
• /
• 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.

• 한국전산유체공학회지
• /
• 제10권3호
• /
• 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년도 추계 학술대회논문집
• /
• 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.