• 한국수자원학회논문집
• /
• 제30권2호
• /
• pp.143-154
• /
• 1997

이송항에는 5차 보간다항식을 사용하는 Holly-Pressmann 기법을, 확산항에는 Hobson 등이 제안한 양해법을 사용하는 연산자 분리기법을 사용하여 1차원 이송-확산방정식의 수치모형을 제안하였다. 제안된 모형을 검정하기 위하여 일정한 유속과 종확산계수를 갖는 순간적으로 부하된 오염원의 경우와 상류단에 연속적인 오염원을 갖는 경우에 대하여 본 모형의 해를 해석해와 기존의 모형으로부터 구한 해를 비교검토하였다. Courant 수와 Peclet 수를 가진 경우에 대한 수치해석을 통하여, 본 모형이 Courant 수가 1보다 큰 경우에 대해서도 안정된 해를 제공함을 알 수 있었으며, 해석해가 존재하는 경우에 본 모형을 적용하여 얻은 수치해와 비교한 바 전반적으로 잘 일치하였다. 본 모형의 확산항에 사용된 양해법에서는 일반적인 양해법의 단점인 계산시간간격의 제약이 상당히 완화되어 상대적으로 큰 계산시간간격에 대해서도 양호한 결과를 보였다.

• 대한수학회지
• /
• 제38권2호
• /
• pp.311-319
• /
• 2001

We show that it is possible to associate diffusion processes to second order perturbations of the Ornstein-Uhlenbeck operator L on the Wiener space of the form L = L + 1/2∑L$^2$(sub)ξ(sub)$\kappa$ where the ξ(sub)$\kappa$ are "tangent processes" (i.e., semimartingales with antisymmetric diffusion coefficients).

• Journal of applied mathematics & informatics
• /
• 제15권1_2호
• /
• pp.415-423
• /
• 2004

In allelic model X = ($x_1,\;x_2,...x_{d}$), $M_f(t)$= f(p(t)) - ${{\int}^{t}}_0$Lf(p(t))ds is a P-martingale for diffusion operator L under the certain conditions. In this note, we can show uniqueness of martingale problem associated with mean vector and obtain a complete description of ergodic property by using of the semigroup method.

• Journal of applied mathematics & informatics
• /
• 제13권1_2호
• /
• pp.393-404
• /
• 2003

In general, a differential operator can be used as a tool of treating the local properties of given function. However, when the given function is varied with high frequency and has irregular form with non-stationary evolution it may not act its role sufficiently as in case of nowhere differentiable curves. In this paper we introduce a multi-scale derivative as a form of weakened global derivative so that it may explain its semi global diffusion properties as well as local ones for the various irregular diffusion phenomena.

• Nuclear Engineering and Technology
• /
• 제54권4호
• /
• pp.1195-1205
• /
• 2022

Recent developments on spectral diffusion algorithms, i.e., algorithms which exploit the projection of the solution on the eigenfunctions of the Laplacian operator, demonstrated their effective applicability in fast transient conditions. Nevertheless, the numerical error introduced by these algorithms, together with the uncertainties associated with model parameters, may impact the reliability of the predictions on short-lived volatile fission product release from nuclear fuel. In this work, we provide an upper bound on the numerical error introduced by the presented spectral diffusion algorithm, in both constant and time-varying conditions, depending on the number of modes and on the time discretization. The definition of this upper bound allows introducing a methodology to a priori bound the numerical error on short-lived volatile fission product retention.

• 대한수학회지
• /
• 제34권3호
• /
• pp.731-770
• /
• 1997

We study Diriclet forms and related subjects for the Gibbs measures of classical unbounded sping systems interacting via potentials which are superstable and regular. For any Gibbs measure $\mu$, we construct a Dirichlet form and the associated diffusion process on $L^2(\Omega, d\mu), where \Omega = (R^d)^Z^\nu$. Under appropriate conditions on the potential we show that the Dirichlet operator associated to a Gibbs measure $\mu$ is essentially self-adjoint on the space of smooth bounded cylinder functions. Under the condition of uniform log-concavity, the Gibbs measure exists uniquely and there exists a mass gap in the lower end of the spectrum of the Dirichlet operator. We also show that under the condition of uniform log-concavity, the unique Gibbs measure satisfies the log-Sobolev inequality. We utilize the general scheme of the previous works on the theory in infinite dimensional spaces developed by e.g., Albeverio, Antonjuk, Hoegh-Krohn, Kondratiev, Rockner, and Kusuoka, etc, and also use the equilibrium condition and the regularity of Gibbs measures extensively.

• Nuclear Engineering and Technology
• /
• 제41권7호
• /
• pp.885-892
• /
• 2009

This manuscript will discuss a numerical method where the six equations of two-phase flow, the solid heat conduction equations, and the two equations that describe neutron diffusion and precursor concentration are solved together in a tightly coupled, nonlinear fashion for a simplified model of a nuclear reactor core. This approach has two important advantages. The first advantage is a higher level of accuracy. Because the equations are solved together in a single nonlinear system, the solution is more accurate than the traditional "operator split" approach where the two-phase flow equations are solved first, the heat conduction is solved second and the neutron diffusion is solved third, limiting the temporal accuracy to $1^{st}$ order because the nonlinear coupling between the physics is handled explicitly. The second advantage of the method described in this manuscript is that the time step control in the fully implicit system can be based on the timescale of the solution rather than a stability-based time step restriction like the material Courant limit required of operator-split methods. In this work, a pilot code was used which employs this tightly coupled, fully implicit method to simulate a reactor core. Results are presented from a simulated control rod movement which show $2^{nd}$ order accuracy in time. Also described in this paper is a simulated rod ejection demonstrating how the fastest timescale of the problem can change between the state variables of neutronics, conduction and two-phase flow during the course of a transient.

• 대한수학회지
• /
• 제33권4호
• /
• pp.823-855
• /
• 1996

We study Dirichlet forms and the associated diffusion processes for the Gibbs measures related to the quantum unbounded spin systems (lattice boson systems) interacting via superstable and regular potentials. This work is a continuation of the author's previous study on the classical systems [LPY] to the quantum cases. In [LPY], we constructed Dirichlet forms and the associated diffusion processes for the Gibbs measures of classical unbounded spin systems. Furthermore, we also showed the essential self-adjointness of the Dirichlet operator and the log-Sobolev inequality for any Gibbs measure under appropriate conditions on the potentials. In this atudy we try to extend the results of the classical systems to the quantum cases. Because of some technical difficulties, we are only able to construct a Dirichlet form and the associated diffusion process for any Gibbs measure of the quantum systems. We utilize the general scheme of the previous work on the theory in infinite dimensional spaces [AH-K1-2, AKR, AR1-2, Kus, MR, Ro, Sch] and the ideas we employed in our study of the calssical systems ]LPY].

• Journal of applied mathematics & informatics
• /
• 제17권1_2_3호
• /
• pp.283-296
• /
• 2005

It is well known that the ordering of the unknowns can have a significant effect on the convergence of a preconditioned iterative method and on its implementation on a parallel computer. To do so, we introduce a block red-black coloring to increase the degree of parallelism in the application of the block ILU preconditioner for solving sparse matrices, arising from convection-diffusion equations discretized using the finite difference scheme (five-point operator). We study the preconditioned PGMRES iterative method for solving these linear systems.

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

In the vortex particle method based on the vorticity-velocity formulation for solving the Wavier-Stokes equations, the unsteady, incompressible, viscous laminar flow over a NACA 0012 foil is simulated. By applying an operator-splitting method, the 'convection' and 'diffusion' equations are solved sequentially at each time step. The convection equation is solved using the vortex particle method, and the diffusion equation using the particle strength exchange(PSE) scheme which is modified to avoid a spurious vorticity flux. The scheme is improved for variety body shape using one image layer scheme. For a validation of the present method, we illustrate the early development of the viscous flow about an impulsively started NACA 0012 foil for Reynolds number 550.