• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 1996년도 춘계학술발표회논문집(1)
• /
• pp.131-136
• /
• 1996

An analytic nodal expansion method has been derived for the multigroup neutron diffusion equation in 2-D cylindrical(R-Z) coordinate. In this method we used the second order Legendre polynomials for source, and transverse leakage, and then the diffusion eqaution was solved analytically. This formalism has been applied to 2-D LWR model. $textsc{k}$$_{eff}$, power distribution, and computing time have been compared with those of ADEP code(finite difference method). The benchmark showed that the analytic nodal expansion method in R-Z coordinate has good accuracy and quite faster than the finite difference method. This is another merit of using R-Z coordinate in that the transverse integration over surfaces is better than the linear integration over length. This makes the discontinuity factor useless.s.

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

• 한국경영과학회지
• /
• 제38권1호
• /
• pp.201-216
• /
• 2013

Diffusion is a basic mathematical tool for modern financial engineering. The theory of the estimation methods for diffusion models is an important topic of the financial engineering. Many researches have been tried to apply the likelihood estimation method for estimating diffusion models. However, the likelihood estimation method for diffusion is complicated and needs much amount of computing. In this paper we develop the estimation methods which are simple enough to be compared to the Euler approximation method, and efficient enough statistically to be compared to the likelihood estimation method. We devise pseudo-likelihood and propose the maximum pseudo-likelihood estimation methods. The pseudo-likelihoods are obtained by approximating the transition density with normal distributions. The means and the variances of the distributions are obtained from the delta expansion suggested by Lee, Song and Lee (2012). We compare the newly suggested estimators with other existing estimators by simulation study. From the simulation study we find the maximum pseudo-likelihood estimator has very similar properties with the maximum likelihood estimator. Also the maximum pseudo-likelihood estimator is easy to apply to general diffusion models, and can be obtained by simple numerical steps.

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 1995년도 추계학술발표회논문집(1)
• /
• pp.29-34
• /
• 1995

The good features of the analytic function expansion nodal (AFEN) method are utilized to develop a practical scheme jot the multigroup diffusion problems, in combination with the polynomial expansion nodal (PEN) method. The thermal group fluxes exhibiting strong gradients are solved by the AFEN method[1-6], while the fast group fluxes that are smoother than the thermal group fuzes are solved by the PEN method[7-9]. The scheme is applied to a MOX-fuel loaded core with good results.

• 지구물리와물리탐사
• /
• 제17권3호
• /
• pp.129-136
• /
• 2014

전자 탐사는 신호원의 파형에 따라 주파수 영역과 시간 영역법으로 나누어진다. 주파수 영역과 시간 영역은 수학적으로 Fourier 변환 관계에 있으므로, 주파수 영역 자료를 Fourier 변환하여 시간 영역 자료를 얻어낼 수 있다. 즉, 시간 영역 전자 탐사의 모델링 자료는 주파수 영역에서 수행한 모델링 자료의 적절한 변환을 통해 얻어질 수 있다. 따라서 주파수-시간 영역 변환은 전자 탐사에서 매우 중요한 부분이다. 분산 전개법(DEM)은 신속하고 효과적인 주파수-시간 영역 변환 기법 중의 하나이다. 분산 전개법에서는 전자기장은 분산 함수와 분산 시간의 급수로 전개하며, 분산 시간은 주어진 주파수 자료에 의해 결정된다. 특히 적정 분산 시간의 설정은 분산 전개법의 정확성을 결정하는 주요 요소이다. 이 연구에서는 급수 전개에 의해 얻어진 주파수 영역 자료의 오차를 최소화하는 방법을 사용하여 적정 분산 시간의 설정 방법을 개발하였다. 반무한 공간 및 2층 구조 모델에 대하여 이 방법을 적용한 결과, 분산 전개법은 상당히 넓은 시간 대역에서 정확한 결과를 나타냄을 확인하였다.

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

• Journal of applied mathematics & informatics
• /
• 제35권3_4호
• /
• pp.277-302
• /
• 2017

In this paper, we have proposed a numerical method for Singularly Perturbed Boundary Value Problems (SPBVPs) of reaction-diffusion type of third order Ordinary Differential Equations (ODEs). The SPBVP is reduced into a weakly coupled system of one first order and one second order ODEs, one without the parameter and the other with the parameter ${\varepsilon}$ multiplying the highest derivative subject to suitable initial and boundary conditions, respectively. The numerical method combines boundary value technique, asymptotic expansion approximation, shooting method and finite difference scheme. The weakly coupled system is decoupled by replacing one of the unknowns by its zero-order asymptotic expansion. Finally the present numerical method is applied to the decoupled system. In order to get a numerical solution for the derivative of the solution, the domain is divided into three regions namely two inner regions and one outer region. The Shooting method is applied to two inner regions whereas for the outer region, standard finite difference (FD) scheme is applied. Necessary error estimates are derived for the method. Computational efficiency and accuracy are verified through numerical examples. The method is easy to implement and suitable for parallel computing. The main advantage of this method is that due to decoupling the system, the computation time is very much reduced.

• Korean Chemical Engineering Research
• /
• 제57권5호
• /
• pp.652-665
• /
• 2019

Analytical solutions of the reactant concentration inside porous spherical catalytic particles were obtained from unsteady reaction-diffusion equation by applying eigenfunction expansion method. Various surface concentrations as exponentially decaying or oscillating function were considered as boundary conditions to solve the unsteady partial differential equation as a function of radial distance and time. Dirac delta function was also used for the instantaneous injection of the reactant as the surface boundary condition to calculate average reactant concentration inside the particles as a function of time by Laplace transform. Besides spherical morphology, other geometries of particles, such as cylinder or slab, were considered to obtain the solution of the reaction-diffusion equation, and the results were compared with the solution in spherical coordinate. The concentration inside the particles based on calculation was compared with the bulk concentration of the reactant molecules measured by photocatalytic decomposition as a function of time.

• 한국세라믹학회지
• /
• 제24권5호
• /
• pp.484-490
• /
• 1987

Diffusion coupling experiment was done to study expansion of body and soild reaction in BaCO3-TiO2 system. Specimen of BaCO3 and TiO2 was formed with Pt-mark's method. Each specimen was fired at interval of 25℃ from 900℃ to 1000℃ for 2hrs. After that, specimen was fixed with resin and polished. Product layers of specimen were observed with SEM and EDS. The result were following; 1. Diffusion component is Ba2+, which diffuse toward TiO2. 2. Large crack between layer of BaCO3 and Ba2TiO4 was generated because of difference of thermal expansion coefficient. 3. Ba2TiO4 is formed to TiO2 body by the reaction of BaTiO3 and BaO and its structure is very porous. 4. BaTiO3 changes immediately to Ba2TiO4 by the reaction of BaO. But BaTiO3 which formed by the reaction of TiO2 and Ba2TiO4 exsists as layer because the diffusion distance of Ba2+ is far.

• Fisheries and Aquatic Sciences
• /
• 제7권2호
• /
• pp.90-95
• /
• 2004

An analytical model for predicting the convection-diffusion of solute dumped in a homogeneous open sea of constant water depth has been developed in a time-integral form. The model incorporates spatially uniform, uni-directional, mean and oscillatory currents for horizontal convection, the settling velocity for the vertical convection, and the anisotropic turbulent diffusion. Two transformations were introduced to reduce the convection-diffusion equation to the Fickian type diffusion equation, and then the Galerkin method was then applied via the expansion of eigenfunctions over the water column derived from the Sturm-Liouville problem. A series of calculations has been performed to demonstrate the applicability of the model.