• Journal of applied mathematics & informatics
• /
• v.6 no.2
• /
• pp.455-466
• /
• 1999

In this note we present a numerical scheme for finding an approxximate solution of an equation which can be viewed as a model for spatial diffusion of age-depednent biological populations. Discretization of the model yields a linear system with a block tridi-agonal matrix. Our main concern will be discussion of stability for this scheme by examining the eigenvalues of the block tridiagonal matrix. Numerical results are presented.

• Computers and Concrete
• /
• v.10 no.5
• /
• pp.523-539
• /
• 2012

Silica fume has long been used as a mineral admixture to improve the durability and produce high strength and high performance concrete. And in marine and coastal environments, penetration of chloride ions is one of the main mechanisms causing concrete reinforcement corrosion. In this paper, we proposed a numerical procedure to predict the chloride diffusion in a hydrating silica fume blended concrete. This numerical procedure includes two parts: a hydration model and a chloride diffusion model. The hydration model starts with mix proportions of silica fume blended concrete and considers Portland cement hydration and silica fume reaction respectively. By using the hydration model, the evolution of properties of silica fume blended concrete is predicted as a function of curing age and these properties are adopted as input parameters for the chloride penetration model. Furthermore, based on the modeling of physicochemical processes of diffusion of chloride ion into concrete, the chloride distribution in silica fume blended concrete is evaluated. The prediction results agree well with experiment results of chloride ion concentrations in the hydrating concrete incorporating silica fume.

• Journal of Environmental Health Sciences
• /
• v.31 no.6
• /
• pp.489-497
• /
• 2005

A numerical method for the mathematical water modeling in 2-dimensional flow has been developed. The model based on a split operator technique, in which, the advection term is calculated using the upwind scheme. The diffusion term is one- dimensionalized and calculated using Crank-Nicholson's implicit finite difference scheme to reduce the numerical errors from large time steps and variable spacings. It also provides a relatively simple and economic method for more accurate simulation of pollutant dispersion. Water depths and flow velocities in the Boreyong reservoir during the normal water periods were predicted by numerical experiments with a 2-dimensional flow model so as to provide current field data for the study of advection and diffusion of pollutants. Developed 2-dimensional water quality model is applied to Boreyong reservoir to simulate a spatial and periodical changes of water quality.

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.43 no.7
• /
• pp.1127-1134
• /
• 1994

The analytical model for the induced Hall field in the magnetotransistor considering the diffusion of carriers has been proposed and verified by experiment and simulation. Previous models for the induced Hall field in the magnetotransistor do not consider the influence of the diffusion carrier transport. However, the carrier diffusion in the base region of magnetotransistors cannot be neglected and should be considered to evaluated the Hall field in the magnetotransistors accurately. We have measured the Hall voltage in the base region of the fabricated magnetotransistors. The measured values have been compared with the numerical results evaluated from our diffusion model as well as the calculated results from the conventional model. The evaluated Hall voltage from the diffusion model agrees well with the measured values while the sign of the Hall voltage calculated by the conventional model is opposite to that of the measured values in the saturation region. This discrepancy is due to the fact that the diffusion model considers the carrier diffusion while the conventional one does not. The Hall field model including the influence of carrier diffusion may be an important tool to optimize the device structure and to understand the operating principle of the magnetotransistor.

• Nuclear Engineering and Technology
• /
• v.54 no.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.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.25 no.12
• /
• pp.660-667
• /
• 2013

This study presents a numerical investigation of the heat and mass transfer kinetics of a fin-tube-type adsorption bed using a two-dimensional numerical model with silica-gel and water as the adsorbent and refrigerant pair. The performance is strongly affected by the heat and mass transfer in the adsorption bed, but the details of the mass transfer kinetics remain unclear. The validity of inter-particle models used to simulate mass-transfer kinetics were examined, such as a constant pressure model and non-constant pressure model, and the valid ranges of the diffusion ratio for each model are proposed. The COP and SCP have been numerically calculated as the performance indexes according to the diffusion ratio. The constant pressure model, which is commonly used in previous research, was found to be valid only in a limited range of diffusion ratio.

• The Pure and Applied Mathematics
• /
• v.22 no.2
• /
• pp.159-168
• /
• 2015

Abstract. We propose a fast and robust finite difference method for Merton's jump diffusion model, which is a partial integro-differential equation. To speed up a computational time, we compute a matrix so that we can calculate the non-local integral term fast by a simple matrix-vector operation. Also, we use non-uniform grids to increase efficiency. We present numerical experiments such as evaluation of the option prices and Greeks to demonstrate a performance of the proposed numerical method. The computational results are in good agreements with the exact solutions of the jump-diffusion model.

• Journal of computational fluids engineering
• /
• v.6 no.4
• /
• pp.15-25
• /
• 2001

Numerical prediction of the diffusion controlled transition in a turbine gas pass is important because it can change the local heat transfer rate over a turbine blade as much as three times. In this study, the gas flow over turbine blade is simplified to the flat plate boundary layer, and an adaptive grid scheme redistributing grid points within the computation domain is proposed with a great emphasis on the construction of the grid control function. The function is sensitized to the second invariant of the mean strain tensor, its spatial gradient, and the interaction of pressure gradient and flow deformation. The transition process is assumed to be described with a κ-ε turbulence model. An elliptic solver is employed to integrate governing equations. Numerical results show that the proposed adaptive grid scheme is very effective in obtaining grid independent numerical solution with a very low grid number. It is expected that present scheme is helpful in predicting actual flow within a turbine to improve computation efficiency.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.18 no.1
• /
• pp.27-41
• /
• 2014

In this paper, we present an efficient numerical method for multiphase image segmentation using a multiphase-field model. The method combines the vector-valued Allen-Cahn phase-field equation with initial data fitting terms containing prescribed interface width and fidelity constants. An efficient numerical solution is achieved using the recently developed hybrid operator splitting method for the vector-valued Allen-Cahn phase-field equation. We split the modified vector-valued Allen-Cahn equation into a nonlinear equation and a linear diffusion equation with a source term. The linear diffusion equation is discretized using an implicit scheme and the resulting implicit discrete system of equations is solved by a multigrid method. The nonlinear equation is solved semi-analytically using a closed-form solution. And by treating the source term of the linear diffusion equation explicitly, we solve the modified vector-valued Allen-Cahn equation in a decoupled way. By decoupling the governing equation, we can speed up the segmentation process with multiple phases. We perform some characteristic numerical experiments for multiphase image segmentation.

• Communications of the Korean Mathematical Society
• /
• v.26 no.2
• /
• pp.321-338
• /
• 2011

In this paper, a class of delayed epidemic model with diffusion is investigated. By analyzing the associated characteristic transcendental equation, its linear stability is investigated and Hopf bifurcation is demonstrated. Some explicit formulae determining the stability and the direction of the Hopf bifurcation periodic solutions bifurcating from Hopf bifurcations are obtained by using the normal form theory and center manifold theory. Some numerical simulation are also carried out to support our analytical findings. Finally, biological explanations and main conclusions are given.