• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.63-72
• /
• 1998

The numerical simulations of flowfield and pollutant dispersion over two-dimensional hills of various shapes are described. The Reynolds-averaged Wavier-Stokes equations and concentration diffusion equation based on the gradient diffusion theory have been applied to the atmospheric shear flow over the bell-shaped hills which are basic components of the complex terrain. The flow characteristics such as velocity profiles of the geophysical boundary layer, speed-up phenomena, mean pollutant concentration profiles are compared with experimental data to validate the present numerical procedure and it has been found that the present numerical results agree well with experiments and other numerical data. It has been also found that the distributions of ground level concentration are strongly influenced by the source location and height.

• Progress in Medical Physics
• /
• v.18 no.3
• /
• pp.161-166
• /
• 2007

A model theory of passive-targeted drug delivery in sphere-shaped solid tumors is introduced on the basis of Fick's law of diffusion, with appropriate boundary and initial conditions. For a uniform initial concentration inside the tumor, the concentration is obtained as a function of time and radial position. The concentration is shown to approach the equilibrium distribution as the time elapses, as is expected by the Gedanken Experiment. The time-evolution rate is found to be determined by the diffusion coefficient of the drug in the tissue, the size of the tumor, the volume of the drug-injected region, and the concentration gradient at the boundary.

• Journal of Korean Society for Atmospheric Environment
• /
• v.13 no.5
• /
• pp.383-396
• /
• 1997

Numerical prediction of the pollutant dispersion over a two-dimensional hilly terrain is presented. The dispersion model used in the present work is based on the gradient diffusion theory and the finite-volume method on a non-orthogonal boundary-fitted grid system. The numerical model is validated by comparing the results with the available experimental data for the flat-floor dispersion within a turbulent boundary-layer. The numerical error analysis is performed based on the guideline of Kasibhatla et al.(1988) for the elevated-source dispersion in the flat-floor boundary layer having a power-law velocity and linear eddy-diffusivity profile. The influences of the two-dimensional hilly terrain on the dispersion from a continuously released source are numerically investigated by changing the emission locations and heights. It is found that the distributions of ground-level concentration are strongly influenced by the source location and the emission height. Hence, the terrain amplification factor is greatly enhanced when the pollutant source is located within a flow separation region. Dispersion from a source of short duration is also simulated and the duration time of the pollutant is compared at several downstream locations on a hilly terrain. The results of the numerical prediction are applied to the evaluation of environmental impacts due to the automobile exhausts at the seashore highway with a parallel mountain range.

• Bulletin of the Korean Chemical Society
• /
• v.13 no.4
• /
• pp.398-404
• /
• 1992

We have investigated the kinetics of diffusion-influenced bimolecular reactions in which one reactant has an internal mode, called the gating mode, that activates or deactivates its reactivity intermittently. The rate law and an expression for the time-dependent rate coefficient have been obtained from the general formalism based on the hierarchy of kinetic equations involving reactant distribution functions. The analytic expression obtained for the steady-state reaction rate constant coincides with the one obtained by Szabo et al., who derived the expression by employing the conventional concentration-gradient approach. For the time-dependent reaction rate coefficient, we obtained for the first time an exact analytic expression in the Laplace domain which was then inverted numerically to give the time-domain results.

• Journal of Preventive Medicine and Public Health
• /
• v.42 no.2
• /
• pp.143-150
• /
• 2009

Objectives : South Korea has experienced unprecedented ups and downs in the sex ratio at birth(SRB), which has been a unique phenomenon in the last two decades. However, little is known about socioeconomic factors that influence the SRB. Employing the diffusion theory by Rogers, this study was undertaken to examine the trends in social variations in the SRB from 1981 to 2004 in Korea. Methods : The data was taken from Vital Birth Statistics for the period from 1981-2004. We computed the annual male proportion of live births according to the parental education(university, middle/high school, primary) and occupation(non-manual, manual, others). Logistic regression analysis was employed to estimate the odds ratios of male birth according to social position for the equidistant three time periods(1981-1984, 1991-1994, and 2001-2004). Results : An increased SRB was detected among parents with higher social position before the mid 1980s. Since then, however, a greater SRB was found for the less educated and manual jobholders. The inverse social gradient for the SRB was most prominent in early 1990s, but the gap has narrowed since the late 1990s. The mother's socioeconomic position could be a sensitive indicator of the social variations in the sex ratio at birth. Conclusions : Changes in the relationship of parental social position with the SRB were detected during the 1980-2004 in Korea. This Korean experience may well be explained by diffusion theory, suggesting there have been socioeconomic differences in the adoption and spread of sex-detection technology.

• Journal of computational fluids engineering
• /
• v.6 no.4
• /
• pp.1-7
• /
• 2001

Wind flow perturbations, recirculations and turbulence generated by buildings often dominate air pollutant distributions around buildings. This paper describes dispersion of contaminants in the vicinity of a building by solving the concentration equation based on previously simulated wind flow field. Turbulence closure is achieved by using the standard k-ε two-equation model. The paper shows application of the CIP method for solving a species concentration equation of contaminant gas around a rectangular building for two different sources under conditions of neutral atmospheric stratification. Results have been compared to the experimental data and the previous numerical results by hybrid scheme. The computational results of concentration profiles by the CIP method agree well with experimental data.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.723-728
• /
• 2001

Wind flow perturbations, recirculations and turbulence generated by buildings often dominate air pollutant distributions around buildings. This paper describes dispersion of contaminants in the vicinity of a building by solving the concentration equation based on previously simulated wind flow field. Turbulence closure is achieved by using the standard k-e two-equation model. The paper shows application of the CIP method for solving a species concentration equation of contaminant gas around a rectangular building for two different sources under conditions of neutral atmospheric stratification. Results have been compared to the experimental data and the previous numerical results by hybrid scheme. The computational results of concentration profiles by the CIP method agree well with experimental data.

• Journal of KIISE:Computer Systems and Theory
• /
• v.32 no.11_12
• /
• pp.692-704
• /
• 2005

With the improvement of computer hardware, GPUs(Graphics Processor Units) have tremendous memory bandwidth and computation power. This leads GPUs to use in general purpose computation. Especially, GPU implementation of compute-intensive physics based simulations is actively studied. In the solution of differential equations which are base of physics simulations, tridiagonal matrix systems occur repeatedly by finite-difference approximation. From the point of view of physics based simulations, fast solution of tridiagonal matrix system is important research field. We propose a fast GPU implementation for the solution of tridiagonal matrix systems. In this paper, we implement the cyclic reduction(also known as odd-even reduction) algorithm which is a popular choice for vector processors. We obtained a considerable performance improvement for solving tridiagonal matrix systems over Thomas method and conjugate gradient method. Thomas method is well known as a method for solving tridiagonal matrix systems on CPU and conjugate gradient method has shown good results on GPU. We experimented our proposed method by applying it to heat conduction, advection-diffusion, and shallow water simulations. The results of these simulations have shown a remarkable performance of over 35 frame-per-second on the 1024x1024 grid.

• Journal of Mechanical Science and Technology
• /
• v.17 no.3
• /
• pp.440-448
• /
• 2003

In the present paper, numerical simulations of buoyant plume dispersion in a neutral and stable atmospheric boundary layer have been carride out. A Lagrangian Stochastic Model (LSM) with a Non-Linear Eddy Viscosity Model (NLEVM) for turbulence is used to generate a Reynolds stress field as an input condition of dispersion simulation. A modified plume-rise equation is included in dispersion simulation in order to consider momentum effect in an initial stage of plume rise resulting in an improved prediction by comparing with the experimental data. The LSM is validated by comparing with the prediction of an Eulerian Dispersion Model (EDM) and by the measured results of vertical profiles of mean concentration in the downstream of an elevated source in an atmospheric boundary layer. The LSM predicts accurate results especially in the vicinity of the source where the EDM underestimates the peak concentration by 40% due to inherent limitations of gradient diffusion theory. As a verification study, the LSM simulation of buoyant plume dispersions under a neutral and stable atmospheric condition is compared with a wind-tunnel experiment, which shows good qualitative agreements.

• The KIPS Transactions:PartB
• /
• v.10B no.5
• /
• pp.521-526
• /
• 2003

In this paper, we propose a geometric active contour model based on intensity information and curve evolution for detecting region boundaries. We put boundary extraction problem as the minimization of the difference between the average intensity of the region and the intensity of the expanding closed curves. We used level set theory to implement the curve evolution for optimal solution. It offered much more freedom in the initial curve position than a general active contour model. Our methods could detect regions whose boundaries are not necessarily defiened by gradient compared to general edge based methods and detect multiple boundaries at the same time. We could improve the result by using anisotropic diffusion filter in image preprocessing. The performance of our model was demonstrated on several data sets like CT and MRI medical images.