• Journal of computational fluids engineering
• /
• v.9 no.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.03a
• /
• 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
• /
• v.28 no.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.

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

• 한국연소학회:학술대회논문집
• /
• 2004.11a
• /
• 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.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.7 no.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.

• Proceedings of the Korean Nuclear Society Conference
• /
• 2001.10a
• /
• pp.55-55
• /
• 2001

• Proceedings of the Korean Nuclear Society Conference
• /
• 2001.10a
• /
• pp.58.2-58
• /
• 2001

• Proceedings of the Korean Nuclear Society Conference
• /
• 2002.10a
• /
• pp.51.2-51
• /
• 2002

• Nuclear Engineering and Technology
• /
• v.49 no.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.