• Membrane Journal
• /
• v.16 no.4
• /
• pp.252-258
• /
• 2006

A numerical analysis was performed for a separation process of carbon dioxide from a flue gas stream using polyethersulfone hollow fiber membranes. Countercurrent flow governing equations were regarded to be two point boundary-value problem and the nonlinear ordinary differential equation were simultaneously solved using the finite- difference method. A computer program was developed using the Compaq Visual Fortran 6.6 software. The carbon dioxide permeate driving force and the fred gas residence time at the inside of membrane were found to be very important factors affecting the permeation characteristics of carbon dioxide. The carbon dioxide concentration in the permeate and the flow rate of the permeate were found to be slightly larger by a few percent with a countercurrent flow analysis than those with a cocurrent flow analysis.

• Journal of the Korean Society of Food Science and Nutrition
• /
• v.26 no.4
• /
• pp.704-713
• /
• 1997

This study was conducted to describe qualitatively the protein metabolism of pigs during growth depending on the feed protein quality and to describe quantitatively amino acids requirements, using a simulation model. The used model has a non-linear structure. In the used model, the protein utilization system of a pig, which is in the non-steady-state, is described with 15 flux equations and 11 differential equations and is composed with two compartments. Protein deposition(g/day) of pigs on the 30th, 60th, 90th, and 120th day of feeding duration with three-quality protein, beginning with body weight 20kg, were calculated according to the empirical model, PAF(the product of amino acid functions) of Menke, and was used as object function for the simulation. The mean of relative difference between the simulated protein deposition and PAF calculated values, lied in a range of 8.8%. The simulated protein deposition showed different behavior according to feed protein quality. In the high-quality protein, it showed paraboloidal form with extending growth simulation up to 150eh day. So the maximum of protein deposition was acquired on the 105th day of simulate growth time and then it decreased fast. In the low-quality protein, this form of protein deposition in the course of simulated growth did not appear until 150th day. The simulated protein mass also showed a difference in accordance with feed protein quality. The difference was small on the 30th day of simulated growth, but with duration of the simulated growth it was larger. On the 150th day the simulated protein deposition of high quality protein was 1.5 times higher as compared to the low-quality protein. The simulated protein synthesis and break-down rates(g/day) in the whole body showed a parallel behavior in the course of growth, according to feed protein quality. It was found that the improvement of feed protein quality increased protein deposition in the whole body through a increase of both protein synthesis and breakdown during growth. Also protein deposition efficiency, which was calculated from simulated protein deposition and protein synthesis, showed a difference in dependence on the protein qualify of feed protein. The protein deposition efficiency was higher in pigs fed with high quality protein, especially at the simulation time 30th day. But this phenomena disappeared with growth, so on the 150th day of growth, the protein deposition of the high feed protein quality was lowest among the three different quality of feed protein. The simulated total requirement of the 10 essential amino acids for the growth of pigs was 28.1(g/100g protein), similar to NRC. The requirement of lysine was 4.2(g/100g protein).

• 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 Biomedical Engineering Research
• /
• v.20 no.2
• /
• pp.173-182
• /
• 1999

Contrast-enhanced CT has an important role in assessing liver lesions, the optimal protocol to get most effective result is not clear. The mein goal when deciding injention protocol is to optimize lesion detectability with rapid scanning when lesion to liver contrast is maximum. For this purpose, we developed a physiological model of the contrast medium enhancement based on the compartment modeling and pharmacokinetics. Blood supply to liver is achieved in two paths. This dual supply characteristic distinguishes the CT enhancement of liver from that of the other organs. The first path is by hepatic artery and to second, by portal vein. However, it is assumed that only gepatic artery can supply blood to hepatocellular carcinoma(HCC) compartment, thus, the difference of contrast enhancement is resulted between normal liver tissue and hepatic tumor. By solving differential equations for each compartment simultaneously using the computer program Matlab, CT contrast-enhancement curves were simulated. The simulated enhancement curves for aortic, hepatic, portal vein, and HCC compartments were compared with the mean enhancement curves from 24 patients exposed to the same protocols as the simulation. These enhancement curves showed a good agreement. Furthermore, we simulated lesion-to-liver curves for various injection protocols, and the effects were analyzed. The variables to be considered in the injection protocol were injection rate, dose, and concentration of contrast material. These data may help to optimize scanning protocols for better diagnosis.

• Journal of Korean Ophthalmic Optics Society
• /
• v.13 no.1
• /
• pp.5-13
• /
• 2008

Purpose: This review article was written to investigate what kind of forces are acting on the contact lens fitted on the cornea and its subsequent motion. Methods: A capillary action-induced force develops in the tear layer between the lens and cornea, which leads to the restoring force due to difference in layer thickness according to lens rotation. The characteristics of the lens movement can be determined by the various factors such as friction between eyelid and lens, acceleration force based on blinking and the restoring force incorporated with the viscous damping force. A mathematical model which consists of the differential equations and their numerical solution was proposed to analyze the damped motion of lenses. The model predicts the time dependence of lenses during and after the blink varying the BC, blink period and eyelid pressure. Results: It was found that both the blink period and lid pressure increases the movement increases because of the enhanced lid friction. As the BC increases the viscous damping reduces due to the lacrimal layer's increase which resulted in the enhanced lens motion. After blink the lens illustrates the damped oscillation because of the restoring force by the increased lacrimal layer thickness and reduced viscous resistance. The time for the lens to return to the equilibrium shortens as the BC increase because of the resistance reduction. Conclusions: The movement of the contact lens is governed by the characteristics of the lacrimal layer between the lens and cornea as well as the lid blink.

• Journal of the Institute of Electronics Engineers of Korea SD
• /
• v.41 no.6
• /
• pp.25-31
• /
• 2004

In this paper, we introduce kinetic Monte Carlo (kMC) methods for simulating diffusion process in nano-scale device fabrication. At first, we review kMC theory and backgrounds and give a simple point defect diffusion process modeling in thermal annealing after ion (electron) implantation into Si crystalline substrate to help understand kinetic Monte Carlo methods. kMC is a kind of Monte Carlo but can simulate time evolution of diffusion process through Poisson probabilistic process. In kMC diffusion process, instead of. solving differential reaction-diffusion equations via conventional finite difference or element methods, it is based on a series of chemical reaction (between atoms and/or defects) or diffusion events according to event rates of all possible events. Every event has its own event rate and time evolution of semiconductor diffusion process is directly simulated. Those event rates can be derived either directly from molecular dynamics (MD) or first-principles (ab-initio) calculations, or from experimental data.

• Analytical Science and Technology
• /
• v.10 no.1
• /
• pp.66-74
• /
• 1997

Influence of the axial dispersion on immobilized enzyme catalytic bed was investigated in order to examine the kinetic behavior of the biocatalysis. The enzyme employed in this study was the tyrosinase(EC 1.14.18.1) immobilized on carbon support : this system requires two substrates of phenol and oxygen. This enzyme has potential application for phenol degradation in waste water. A simulated reactor was a packed-bed reactor of 2.54cm in diameter and 10cm long, loaded with the immobilized carbon particle with an average diameter of $550{\mu}m$. A phenol feed in the strength of 55.5mM(5220ppm) was used to observe the behavior of the immobilized enzyme column at three different dissolved oxygen levels of 0.08445mM(2.7ppm), 0.1689mM(5.4ppm) and 0.3378mM(9.5ppm) with the flow rates in the range of 60(1mL/s) to 180mL/min(3mL/s). Examination of the Biot number and Damkolher numbers of the immobilized system enables us to eliminate the contribution of external mass transfer to set of differential equations derived from the dispersion model. Solution of the equation was finally obtained numerically with the application of the Danckwert boundary conditions and the assumed zero-and first order rates on the non-linear two substrate enzyme kinetics. Higher conversion of phenol was observed at the low flow rates and at the higher oxygen concentration. Comparison of axial dispersion and plug flow model showed that no detectable difference was observed in the column outlet conversion between the axial and the plug flow models which was in complete agreement with the previous studies.