• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.3
• /
• pp.253-266
• /
• 2019

In this article, a systematic solution based on the sequence of expansion method is planned to solve the time-fractional diffusion equation, time-fractional telegraphic equation and time-fractional wave equation in three dimensions using a current and valid approximate method, namely the ADM, VIM, and the NIM subject to the estimate initial condition. By using these three methods it is likely to find the exact solutions or a nearby approximate solution of fractional partial differential equations. The exactness, efficiency, and convergence of the method are demonstrated through the three numerical examples.

• Nuclear Engineering and Technology
• /
• v.56 no.4
• /
• pp.1296-1309
• /
• 2024

The aim of this work is to explore the effect of the double subdiffusion on the stability in BWRs. A BWR novel reduced order model with double subdiffusion effects: reduced order fractional model (DS-F-ROM) to describe the neutron and heat transfer processes was proposed for this study. The double subdiffusion was developed with a fractional-order two-equation model, and with different fractional-orders and relaxation times. The stability analysis was carried out using the root-locus method and change from the s to the W domain and were confirmed using the time-domain evolution of neutron flux for a unit step change in reactivity. The results obtained using the reduced fractional-order model are presented for different anomalous diffusion coefficient values. Results are compared with normal diffusion and P1 equations, which are obtained straightforwardly with DS-ROM when relaxation time tends to zero, and when the anomalous diffusion coefficient tends to one, respectively.

• Journal of the Korean Mathematical Society
• /
• v.53 no.4
• /
• pp.929-967
• /
• 2016

Let p(t, x) be the fundamental solution to the problem $${\partial}^{\alpha}_tu=-(-{\Delta})^{\beta}u,\;{\alpha}{\in}(0,2),\;{\beta}{\in}(0,{\infty})$$. If ${\alpha},{\beta}{\in}(0,1)$, then the kernel p(t, x) becomes the transition density of a Levy process delayed by an inverse subordinator. In this paper we provide the asymptotic behaviors and sharp upper bounds of p(t, x) and its space and time fractional derivatives $$D^n_x(-{\Delta}_x)^{\gamma}D^{\sigma}_tI^{\delta}_tp(t,x),\;{\forall}n{\in}{\mathbb{Z}}_+,\;{\gamma}{\in}[0,{\beta}],\;{\sigma},{\delta}{\in}[0,{\infty})$$, where $D^n_x$ x is a partial derivative of order n with respect to x, $(-{\Delta}_x)^{\gamma}$ is a fractional Laplace operator and $D^{\sigma}_t$ and $I^{\delta}_t$ are Riemann-Liouville fractional derivative and integral respectively.

• Nonlinear Functional Analysis and Applications
• /
• v.27 no.3
• /
• pp.679-689
• /
• 2022

In this paper, we establish sufficient conditions for the existence and uniqueness of solution for a class of initial boundary value problems with Dirichlet condition in regard to a category of fractional-order partial differential equations. The results are established by a method based on the theorem of Lax Milligram.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.16 no.3
• /
• pp.193-204
• /
• 2012

This paper is comparing numerical schemes for a differential equation with convection and fourth-order diffusion. Our model equation is $h_t+(h^2-h^3)_x=-(h^3h_{xxx})_x$, which arises in the context of thin film flow driven the competing effects of an induced surface tension gradient and gravity. These films arise in thin coating flows and are of great technical and scientific interest. Here we focus on the several numerical methods to apply the model equation and the comparison and analysis of the numerical results. The convection terms are treated with well known WENO methods and the diffusion term is treated implicitly. The diffusion and convection schemes are combined using a fractional step-splitting method.

• Nuclear Engineering and Technology
• /
• v.47 no.5
• /
• pp.608-616
• /
• 2015

The migration of silver (Ag) in silicon carbide (SiC) and $^{110m}Ag$ through SiC of irradiated tristructural isotropic (TRISO) fuel has been studied for the past three to four decades. However, there is no satisfactory explanation for the transport mechanism of Ag in SiC. In this work, the diffusion coefficients of Ag measured and/or estimated in previous studies were reviewed, and then pre-exponential factors and activation energies from the previous experiments were evaluated using Arrhenius equation. The activation energy is $247.4kJ{\cdot}mol^{-1}$ from Ag paste experiments between two SiC layers produced using fluidized-bed chemical vapor deposition (FBCVD), $125.3kJ{\cdot}mol^{-1}$ from integral release experiments (annealing of irradiated TRISO fuel), $121.8kJ{\cdot}mol^{-1}$ from fractional Ag release during irradiation of TRISO fuel in high flux reactor (HFR), and $274.8kJ{\cdot}mol^{-1}$ from Ag ion implantation experiments, respectively. The activation energy from ion implantation experiments is greater than that from Ag paste, fractional release and integral release, and the activation energy from Ag paste experiments is approximately two times greater than that from integral release experiments and fractional Ag release during the irradiation of TRISO fuel in HFR. The pre-exponential factors are also very different depending on the experimental methods and estimation. From a comparison of the pre-exponential factors and activation energies, it can be analogized that the diffusion mechanism of Ag using ion implantation experiment is different from other experiments, such as a Ag paste experiment, integral release experiments, and heating experiments after irradiating TRISO fuel in HFR. However, the results of this work do not support the long held assumption that Ag release from FBCVD-SiC, used for the coating layer in TRISO fuel, is dominated by grain boundary diffusion. In order to understand in detail the transport mechanism of Ag through the coating layer, FBCVD-SiC in TRISO fuel, a microstructural change caused by neutron irradiation during operation has to be fully considered.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1601-1606
• /
• 2004

In this paper, the unsteady behavior of the viscous flow field past an impulsively started elliptic cylinder is studied numerically. In order to analyze flow field, we introduce vortex particle method. The vorticity transport equation is solved by fractional step algorithm which splits into convection term and diffusion term. The convection term is calculated with Biot-Savart law, the no-through boundary condition is employed on solid boundaries. The diffusion term is modified based on the scheme of particle strength exchange. The particle redistributed scheme for general geometry is adapted. The flows around an elliptic cylinder are investigated for various attack angles at Reynolds number 200. The comparison between numerical results of present study and experimental data shows good agreements.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.7
• /
• pp.809-817
• /
• 2004

In this paper, the unsteady behavior of the viscous flow field past an impulsively started elliptic cylinder is studied numerically. In order to analyze flow field, we introduce vortex particle method. The vorticity transport equation is solved by fractional step algorithm which splits into convection term and diffusion term. The convection term is calculated with Biot-Savart law, the no-through boundary condition is employed on solid boundaries. The diffusion term is modified based on the scheme of particle strength exchange. The particle redistributed scheme for general geometry is adapted. The flows around an elliptic cylinder are investigated for various attack angles at Reynolds number 200. The comparison between numerical results of present study and experimental data shows good agreements.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.1-1
• /
• 2009

In this presentation, I talk about various fluid simulation methods that have been developed for computer graphics special effects since 1996. They are all based on CFD but sacrifice physical reality for visual plausability and time. But as the speed of computer increases rapidly and the capability of GPU (graphics processing unit) improves, methods for more physical realism have been tried. In this talk, I will focus on four aspects of fluid simulation methods for computer graphics: (1) particle level-set methods, (2) particle-based simulation, (3) methods for exact satisfaction of incompressibility constraint, and (4) GPU-based simulation. (1) Particle level-set methods evolve the surface of fluid by means of the zero-level set and a band of massless marker particles on both sides of it. The evolution of the zero-level set captures the surface in an approximate manner and the evolution of marker particles captures the fine details of the surface, and the zero-level set is modified based on the particle positions in each step of evolution. (2) Recently the particle-based Lagrangian approach to fluid simulation gains some popularity, because it automatically respects mass conservation and the difficulty of tracking the surface geometry has been somewhat addressed. (3) Until recently fluid simulation algorithm was dominated by approximate fractional step methods. They split the Navier-Stoke equation into two, so that the first one solves the equation without considering the incompressibility constraint and the second finds the pressure which satisfies the constraint. In this approach, the first step introduces error inevitably, producing numerical diffusion in solution. But recently exact fractional step methods without error have been developed by fluid mechanics scholars), and another method was introduced which satisfies the incompressibility constraint by formulating fluid in terms of vorticity field rather than velocity field (by computer graphics scholars). (4) Finally, I want to mention GPU implementation of fluid simulation, which takes advantage of the fact that discrete fluid equations can be solved in parallel.

• Journal of Korea Water Resources Association
• /
• v.32 no.6
• /
• pp.607-616
• /
• 1999

A fractional step finite difference model for the longitudinal dispersion of nonconservative contaminants is developed. It is based on splitting the longitudinal dispersion equation into a set of three equations each to be solved over a one-third time step. The fourth-order Holly-Preissmann scheme, an analytic solution, and the Crank-Nicholson scheme are used to solve the equations for the pure advection, the first-order decay, and the diffusion, respectively. To test the model, it is applied to simulate the longitudinal dispersion of continuous source released into a nonuniform flow field as well as the dispersion of an instantaneous source in a uniform flow field. The results are compared with the exact solution and those computed by an existing model. Compared to the existing model which uses Euler method for the first-order decay equation, the present model yield more accurate results as the decay coefficient increases.