• Nuclear Engineering and Technology
• /
• v.53 no.12
• /
• pp.3861-3878
• /
• 2021

The steady-state simplified spherical harmonics equations (SP_N equations) are a higher order approximation to the neutron transport equations than the neutron diffusion equation that also have reasonable computational demands. This work extends these results for the analysis of transients by comparing of two formulations of time-dependent SP_N equations considering different treatments for the time derivatives of the field moments. The first is the full system of equations and the second is a diffusive approximation of these equations that neglects the time derivatives of the odd moments. The spatial discretization of these methodologies is made by using a high order finite element method. For the time discretization, a semi-implicit Euler method is used. Numerical results show that the diffusive formulation for the time-dependent simplified spherical harmonics equations does not present a relevant loss of accuracy while being more computationally efficient than the full system.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.3
• /
• pp.1021-1029
• /
• 1991

본 연구에서는 2차원 정사각형 밀폐공간내에서 흡수 및 방사하는 회기체에 대 한 자연대류-복사 열전달을 P-1 및 P-3 근사법을 이용하고 수치해석을 통하여 유동 및 열전달 특성을 연구하였고 Plank 수, 광학두께 및 벽방사율의 영향을 조사하였다. 또한 P-3 근사해와 비교함으로써 P-1 근사해의 적용범위를 고찰하였다.

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

• Steel and Composite Structures
• /
• v.17 no.5
• /
• pp.561-572
• /
• 2014

In this paper, we present the scattering characteristics of infinite and finite array using method of moment (MoM) with Floquet harmonics and asymptotic waveform evaluation (AWE) technique. First, infinite cylindrical dipole array is analyzed using the MoM with entire domain basis function and cylindrical Floquet harmonics. To provide the validity of results, we fabricated the cylindrical dipole array and measured the transmission characteristics. The results show good agreements. Second, we analyzed the scattering characteristics of finite array. A large simulation time is needed to obtain the scattering characteristics of finite array over wide frequency range because Floquet harmonics can't be applied. So, we used the MoM with AWE technique using Taylor series and Pade approximation to overcome the shortcomings of conventional MoM. We calculated the radar cross section (RCS) as scattering characteristics using the proposed method in this paper and the conventional MoM for finite planar slot array, finite spherical slot array, and finite cylindrical dipole array, respectively. The compared results agree well and show that the proposed method in this paper is good for electromagnetic analysis of finite FSS.

• Journal of the korean Society of Automotive Engineers
• /
• v.14 no.2
• /
• pp.54-64
• /
• 1992

During combustion process in a diesel engine radiation heat transfer is the same order of magnitude as the convection heat transfer. An approximation of heat and momentum source distributions is applied at a level consistent with those used in modelling the soot distribution and the turbulence instead of modelling the fuel spray and the chemical kinetics. This paper illustrates a use of the third order spherical harmonics approximation to the radiative transfer equation and delta-Eddington approximation to the scattering phase function for droplets in the flow. Results are obtained numerically by a time marching finite difference scheme. This study aims to compare heat transfer with convection heat transfer and to investigate the importance of scattering by fuel droplets and of accounting for spatial variations in the extinction coefficient on the radiative heat flux distributions at the walls of a disc shaped diesel engine.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.11
• /
• pp.2999-3007
• /
• 1994

The unsteady numerical simulations have been presented for the laminar natural convection in a partially open compartment. Computations were performed within the domain of the compartment in order to show the thermal radiation and the partially opening effects on the flow fields and heat transfer characteristics. The results were shown for different Planck numbers(0.05~5) and opening ratios(0.25~0.75) being fixed with Ra=$10^5$ and Pr=0.71. Considering the flow which is buoyancy driven from the heated wall, and the buoyancy is not much affected by the further outside region from the opening, the numerical computations have been performed without an outer region by the particular boundary treatments on the flow velocity and temperature at the different partial openings. The confined numerical domain reduced the CPU time and the memory of computer. P-1 approximation of radiative transfer equation was employed with Marshak type boundary conditions along with the pseudo-black body approximation at the partial openings. The numerical results clearly show that the natural convective flow and heat transfer are much affected by increase of thermal radiation particularly from the initial state. When thermal radiation is not much affecting the flow ($PL{\le}1$), it was found that thermal radiation effects are almost negligible.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.6
• /
• pp.1652-1660
• /
• 1990

Combined radiative-convective heat transfer in a hot gas tube flow has been investigated numerically and experimentally. In the numerical analysis, a standard k-.epsilon. model is used for the evaluation of turbulent shear stresses and spherical harmonics method with the Weighted Sum of Gray Gases Model for the solution of radiative transfer equation. In the experimental study measured are the velocity and temperature of the hot gas flow generated by the propane gas combustion, and tude wall heat flux distribution. Numerical results are compared with experimental ones and it is confirmed that P-3 provides quite reliable results in the analysis of the combined radiation-convection system.

• Journal of Advanced Marine Engineering and Technology
• /
• v.19 no.2
• /
• pp.10-19
• /
• 1995

Heat trnasfer for absorbing and emitting media in laminar layer along the cylinders has been analyzed. Governing equation are transformed to local nonsimilarity equations by the dimensional analysis. The effects of the Stark number, Prandtl number, Optical radius and wall emissivity are mainly investigated. For the formal solution a numerical integration is performed and the results are compared with those obtained by P-1 and P-3 approximation. The results show that boundary layers consist of conduction-convection-radiation layer near the wall and convection-radiation layer far from the wall. As the Stark number of wall emissivity increases the local radiative heat flux is increased. The Pradtl number or curvature variations do not affect the radiative heat flux from the wall, but The Prandtl number or wall emissivity variations affect the conduction heat flux. Consequently the total heat flux from the wall are affected by the Prandtl number or wall emissivity variation.

• Journal of The Korean Astronomical Society
• /
• v.56 no.2
• /
• pp.187-194
• /
• 2023

The remote sensing technique of measuring the magnetic field was applied first to sunspots by Hale (1908). Later Babcock (1961) showed that the solar surface magnetic field on a global scale is a dipole in first-order approximation and that this dipole field reverses once every solar cycle. The Wilcox Solar Observatory (WSO) supplies the spherical harmonics coefficients of the solar corona magnetic field of each Carrington Rotation, calculated based on the remotely-sensed photospheric magnetic field of the solar surface. To infer the internal current system producing the global solar coronal magnetic field structure and evolution of the Sun, we calculate the multipole components of the solar magnetic field using the WSO data from 1976 to 2019. The prominent cycle components over the last 4 solar activity cycles are axis-symmetric fields of the dipole and octupole. This implies that the current inversion driving the solar magnetic field reversal originates from the equatorial region and spreads to the whole globe. Thus, a more accurate solar dynamo model must include an explanation of the origin and evolution of such solar internal current dynamics.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.7 no.4
• /
• pp.556-565
• /
• 1995

The physical model of interest is based upon the concentric cylinder, where the outside cylinder is filled with optically thick and high temperature phase change material(PCM). The fluid is flowing through the inside cylinder to transfer the appropriate energy. The fluid is flowing through the inside cylinder to transfer the appropriate energy. The governing equations for the phase change material including internal thermal radiation and for the turbulent transfer fluid have been employed and numerically solved. The optically thick phase change justifies the P-l spherical harmonics approximation, which is believed to be appropriate choice particularly for the much coupled problem like in this study. The solid/liquid interface, temperature distribution within the PCM and the heat flux from the PCM to the transfer fluid have been obtained and compared with those of laminar transfer fluid. The numerical results show that the turbulent transfer fluid accelerates the solid/liquid interface and results in the increase of heat transfer rate from the PCM. The internal thermal radiation within the PCM, however, does not always playa role to increase the heat transfer rate throughout the inside cylinder. It is believed that the combined heat flux has been picked up more in the inflowing area than in the pure conductive phase change material.