• 한국전산유체공학회지
• /
• 제9권2호
• /
• pp.43-51
• /
• 2004

Interaction between fluid flow and thermal radiation has received considerable attention due to its numerous applications in engineering field. In this case the thermofluid properties of radiating fluid vary with the variation of temperature field caused by absorption and emission of radiant heat. To analyze the radiation heat transfer in radiating fluid, the simultaneous solution of the radiative transfer equation (RTE) and the fluid dynamics equations is required. This means that the numerical procedure used for the RTE must be computationally efficient to permit its inclusion in the other submodels, and must be compatible with the other transport equations. The finite volume method (FVM) and the discrete ordinates method (DOM) are usually employed to simulate radiation problems in generalized coordinates. These two representative methods are examined and compared, especially in view of the numerical integration of the radiation intensity over solid angle. The FVM shows better accuracy than the DOM owing to less constraints of the selection of control angle.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2004년도 춘계 학술대회논문집
• /
• pp.62-70
• /
• 2004

Interaction between fluid flow and thermal radiation has received considerable attention due to its numerous applications in engineering field. In this case the thermofluid properties of radiating fluid vary with the variation of temperature field caused by absorption and emission of radiant heat. To analyze the radiation heat transfer in radiating fluid, the simultaneous solution of the radiative transfer equation (RTE) and the fluid dynamics equations is required. This means that the numerical procedure used for the RTE must be computationally efficient to permit its inclusion in the other submodels, and must be compatible with the other transport equations. The finite volume method (FVM) and the discrete ordinates method (DOM) are usually employed to simulate radiation problems in generalized coordinates. These two representative methods are examined and compared, especially in view of the numerical integration of the radiation intensity over solid angle. The FVM shows better accuracy than the DOM owing to less constraints of the selection of control angle.

• Nuclear Engineering and Technology
• /
• 제28권5호
• /
• pp.488-499
• /
• 1996

We have extended the source projection analytic nodal discrete ordinates method (SPANDOM) for more flexible applicability in analysis of hexagonal assembly cores. The method (SPANDOM-FH) does not invoke transverse integration but instead solves the discrete ordinates equation analytically after the source term is projected and represented in hybrid form of high-order polynomials and exponential functions. SPANDOM-FH which treats a hexagonal node as one node is applied to two fast reactor benchmark problems and compared with TWOHEX. The results of comparison indicate that the present method SPANDOM-FH predicts accurately $k_eff$ and flux distributions in hexagonal assembly cores. In addition, SPANDOM-FH gives the continuous two dimensional intranodal scalar flux distributions in a hexagonal node. The reentering models between TWOHEX and SPANDOM were also compared and it was confirmed that SPANDOM's model is more realistic. Through the results of benchmark problems, we conclude that SPANDOM-FH has the sufficient accuracy for the nuclear design of fast breeder reactor (FBR) cores with hexagonal assemblies.

• 에너지공학
• /
• 제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)을 사용하는 경우 광첨두 현상이 완전 제거 또는 극단적으로 감소되었던 이전의 결과를 이론적으로 설명한다.

• 한국연소학회:학술대회논문집
• /
• 한국연소학회 2004년도 제29회 KOSCI SYMPOSIUM 논문집
• /
• pp.49-54
• /
• 2004

Detailed flame structures of the opposed flames formed for different oxidant compositions are studied numerically. The detailed chemical reactions are modeled by using the CHEMKIN code. Only the $CO_2$ and $H_2O$ are assumed to participate by absorbing the radiative energy while all other gases are assumed to be transparent. The discrete ordinates method and the narrow band based WSGGM with a gray gas regrouping technique are applied for modeling the radiative transfer through non-homogeneous and non-isothermal combustion gas mixtures generated by the opposed flow flames. The results show that the different radiation model can cause different results for flame structures and the WSGGM with gray gas regrouping is successful in modeling the opposed flames with non-gray gas mixture. The numerical results show that the increases in $CO_2$ and $H_2O$ compositions cause to reduce the flame temperature and the NO formation.

• 설비공학논문집
• /
• 제7권2호
• /
• pp.266-275
• /
• 1995

A numerical investigation on the conduction-natural convection-surface radiation conjugate heat transfer in the enclosure having substrate and chips has been performed. A 2-dimensional simulation model is developed by considering heat transfer by conduction, convection and radiation. The solutions to the equation of radiative transfer are obtained by the discrete ordinates method using S-4 quadrature. The effects of Rayleigh number and the substrate-fluid thermal conductivity ratio on the cooling of chip are analyzed. The result shows that radiation is the dominant heat transfer mode in the enclosure.

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 2001년도 추계학술발표회요약집
• /
• pp.55-55
• /
• 2001

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 2001년도 추계학술발표회요약집
• /
• pp.58.2-58
• /
• 2001

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 2002년도 추계학술발표회요약집
• /
• pp.51.2-51
• /
• 2002

• Nuclear Engineering and Technology
• /
• 제49권6호
• /
• pp.1114-1124
• /
• 2017

Acceleration/preconditioning strategies available in the SCEPTRE radiation transport code are described. A flexible transport synthetic acceleration (TSA) algorithm that uses a low-order discrete-ordinates ($S_N$) or spherical-harmonics ($P_N$) solve to accelerate convergence of a high-order $S_N$ source-iteration (SI) solve is described. Convergence of the low-order solves can be further accelerated by applying off-the-shelf incomplete-factorization or algebraic-multigrid methods. Also available is an algorithm that uses a generalized minimum residual (GMRES) iterative method rather than SI for convergence, using a parallel sweep-based solver to build up a Krylov subspace. TSA has been applied as a preconditioner to accelerate the convergence of the GMRES iterations. The methods are applied to several problems involving electron transport and problems with artificial cross sections with large scattering ratios. These methods were compared and evaluated by considering material discontinuities and scattering anisotropy. Observed accelerations obtained are highly problem dependent, but speedup factors around 10 have been observed in typical applications.