• 방사성폐기물학회지
• /
• 제16권2호
• /
• pp.211-221
• /
• 2018

중성자 확산방정식에 대해 개발된 노달 확산이론을 단순 우성 중성자 수송방정식에 적용할 수 있는 노달 수송이론을 제시한다. 노달이론으로 다항식전개 노달법과 해석함수전개 노달법을 채택하였고 단순 우성 수송방정식은 수송방정식에 대한 합리적 근사이며 기존의 노달해법이 방향 차분된 단순 우성 수송방정식에 정확히 적용될 수 있음을 수치적으로 확인하였다. 본 연구에서는 방법론 개발이 목적이므로 노드 당 최소한의 미지수를 정의하여 사용했지만 미지수를 추가함으로써 정확도를 증가시킬 수 있음은 기존의 노달 확산이론의 경우와 같다. 즉 중성자 수송방정식에 대해 노달이론을 적용하여 소격격자에 대해 계산 정확성이 확보되고 이는 결국 계산 효율성 증대로 나타난다.

• Nuclear Engineering and Technology
• /
• 제32권6호
• /
• pp.605-617
• /
• 2000

Nodal transport methods are studied for the solution of two dimensional discrete-ordinates, simplified even-parity transport equation(SEP) which is known to be an approximation to the true transport equation. The polynomial expansion nodal method(PEN) and the analytic function expansion nodal method(AFEN）which have been developed for the diffusion theory are used for the solution of the discrete-ordinates form of SEP equation. Our study shows that while the PEN method in diffusion theory can directly be converted without complication, the AFEN method requires a theoretical modification due to the nonhomogeneous property of the transport equation. The numerical results show that the proposed two methods work well with the SEP transport equation with higher accuracies compared with the conventional finite difference method.

• Nuclear Engineering and Technology
• /
• 제37권6호
• /
• pp.525-536
• /
• 2005

The interfacial area transport equation dynamically models the changes in interfacial structures along the flow field by mechanistically modeling the creation and destruction of dispersed phase. Hence, when employed in the numerical thermal-hydraulic system analysis codes, it eliminates artificial bifurcations stemming from the use of the static flow regime transition criteria. Accounting for the substantial differences in the transport mechanism for various sizes of bubbles, the transport equation is formulated for two characteristic groups of bubbles. The group 1 equation describes the transport of small-dispersed bubbles, whereas the group 2 equation describes the transport of large cap, slug or chum-turbulent bubbles. To evaluate the feasibility and reliability of interfacial area transport equation available at present, it is benchmarked by an extensive database established in various two-phase flow configurations spanning from bubbly to chum-turbulent flow regimes. The geometrical effect in interfacial area transport is examined by the data acquired in vertical fir-water two-phase flow through round pipes of various sizes and a confined flow duct, and by those acquired In vertical co-current downward air-water two-phase flow through round pipes of two different sizes.

• 한국지하수토양환경학회지:지하수토양환경
• /
• 제8권3호
• /
• pp.8-12
• /
• 2003

The centrifuge test result on capped sediment was compared to the advection- dispersion equation proposed for one layered to predict contaminant transport parameters. The fitted contaminant transport parameters for the centrifuge test results were one to three orders of magnitude greater than the estimated parameters from the advection-dispersion equation. This indicates that the centrifuge model over estimated the contaminant transport phenomena. Thus, the centrifuge provides a non-conservative approach to modeling contaminant transport. It should be also noted that the advection-dispersion equation used in this study is a one layered model. Two layered modeling approaches are more appropriate for modeling this data since there are two layers with different partitioning coefficients. Further research is required to model the centrifuge test using two-layered advection-dispersion models.

• Nuclear Engineering and Technology
• /
• 제52권3호
• /
• pp.485-498
• /
• 2020

A diffusion synthetic acceleration (DSA) technique for the S_N transport equation discretized with the linear discontinuous expansion method with subcell balance (LDEM-SCB) on unstructured tetrahedral meshes is presented. The LDEM-SCB scheme solves the transport equation with the discrete ordinates method by using the subcell balances and linear discontinuous expansion of the flux. Discretized DSA equations are derived by consistently discretizing the continuous diffusion equation with the LDEM-SCB method, however, the discretized diffusion equations are not fully consistent with the discretized transport equations. In addition, a fine mesh rebalance (FMR) method is devised to accelerate the discretized diffusion equation coupled with the preconditioned conjugate gradient (CG) method. The DSA method is applied to various test problems to show its effectiveness in speeding up the iterative convergence of the transport equation. The results show that the DSA method gives small spectral radii for the tetrahedral meshes having various minimum aspect ratios even in highly scattering dominant mediums for the homogeneous test problems. The numerical tests for the homogeneous and heterogeneous problems show that DSA with FMR (with preconditioned CG) gives significantly higher speedups and robustness than the one with the Gauss-Seidel-like iteration.

• Nuclear Engineering and Technology
• /
• 제54권9호
• /
• pp.3181-3204
• /
• 2022

The variational nodal method for solving the neutron transport equation has evolved over 40 years. Based on a functional form of the Boltzmann neutron transport equation, the method now comprises a complete set of variants that can be employed for different problems. This paper presents an extensive review of the development of the variational nodal method. The emphasis is on summarizing the whole theoretical system rather than validating the methodologies. The paper covers the variational nodal formulation of the Boltzmann neutron transport equation, the Ritz procedure for various application purposes, the derivation of boundary conditions, the extension for adjoint and perturbation calculations, and treatments for anisotropic scattering sources. Acceleration approaches for constructing response matrices and solving the resulting system of algebraic equations are also presented.

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

• 대한기계학회논문집B
• /
• 제22권1호
• /
• pp.1-11
• /
• 1998

For the analysis of phonon heat transfer within short time and spatial scales, conventional macroscopic heat conduction equations with jump boundary conditions are tried and the results are compared to those of equation of phonon radiative transport(EPRT), which is one of microscopic transport equation. In transient state the macroscopic temperatures show far different behavior from EPRT. In steady state the hyperbolic temperatures with temperature jump at the wall from time relaxation model agrees well with EPRT temperatures. Since EPRT is also an approximate form of microscopic transport equation and there are no experimental results to verify the proposed model in this study, we can not conclude whether the approaching method from this study is valid or not. To the authors' knowledge, there are no experimental results available which can be used to test the validity of these models. Such an experiment, while difficult to conduct, would be invaluable.

• 한국전기전자재료학회:학술대회논문집
• /
• 한국전기전자재료학회 2001년도 춘계학술대회 논문집 유기절연재료 전자세라믹 방전플라즈마 연구회
• /
• pp.164-167
• /
• 2001

The electron transport coefficients, the electron drift velocity W, the longitudinal diffusion coefficient $ND_L$ and $D_L/{\mu}$, in pure $CO_2$ were calculated over the wide E/N range from 0.01 to 500 Td at 1 Torr by two-term approximation of the Boltzmann equation for determination of electron collision cross sections set and for quantitative characteristic analysis of $CO_2$ molecular gas. And for propriety of two-term approximation of Boltzmann equation analysis, the calculated results compared with the electron transport coefficients measured by Nakamura.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1995년도 가을 학술발표회 논문집
• /
• pp.76-83
• /
• 1995

Estimation of blasting behavior of explosives is prerequisite in the numerical analysis of blasting works. In this study, blasting pressure is estimated by the finite difference method using the Flux-Corrected Transport Algorithm. To formulate the behavior of blasting gas, the mass conservation equation, the moment conservation equation, the energy conservation equation and the ideal gas state equation are used. The simplified species conservation equation is included to simulate the behavior of reacting explosives. To verify the calculation, the Sod's shock tube problem, the strong shock problem and the reacting problem we used. Numerical results show that the shock wave can be captured by means of the FCT algorithm in the reacting and nonreacting states.