• Journal of the Korean Society of Food Science and Nutrition
• /
• v.28 no.1
• /
• pp.178-190
• /
• 1999

This study was conducted to develop a simulation model for the growth dynamics of pigs and to describe quantitatively protein deposition depending on the amino acid composition of feed protein. In the model it is assumed that the essential processes that determine the utilization of feed protein in the whole body are protein synthesis, breakdown of protein, and oxidation of amino acid. Besides, it is also assumed that occurrence of protein deposition depends on genetic potential and amino acid composition of feed protein. The genetic potential for the protein deposition is the maximum capacity of protein synthesis, being dependent on the protein mass of the whole body. To describe the effect of amino acid composition of feed on the protein deposition, a factor, which consist of ten amino acid functions and lie between 0 and 1, is introduced. Accordingly a model was developed, which is described with 15 flux equations and 11 differential equations and is composed of two compartments. The model describes non linear structure of the protein utilization system of an organism, which is in non steady state. The objective function for the simulation was protein deposition(g/day) cal culated according to the empirical model, PAF(product of amino acid functions) of Menke. The mean of relative difference between the simulated protein deposition and PAF calculated values, lied in a range of 11.8%. The simulated protein synthesis and breakdown rates(g/day) in the whole body showed a parallel behavior in the course of growth.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.9 no.1
• /
• pp.1-7
• /
• 1998

The near field radiated from the monopole antenna of the cellular phone was calculated by using the modified finite difference time domain algorithm derived from the integral form of Maxwell's equations. Substituting the near field value into the differential form of Maxwell's equations, SAR's distribution in the human head was obtained. The human head was simulated by a model of 800,000 block cells with dielectric constant and conductivity. The cell size was taken to be 0.5 cm. the transmitted power of the cellular phone was assumed to be 0.6 watts at the frequency of 833 MHz. The distance between the head and the cellular phone was 2.0 cm, the maximum SAR induced in the human head was about 1.5 W/kg and was below the IEEE's upper safety limit of 1.6 W/kg.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.5 no.2
• /
• pp.63-69
• /
• 2001

In this study, we shall be concerned with computing in parallel a few of the smallest eigenvalues and their corresponding eigenvectors of the eigenvalue problem, $Ax={\lambda}Bx$, where A is symmetric, and B is symmetric positive definite. Both A and B are large and sparse. Recently iterative algorithms based on the optimization of the Rayleigh quotient have been developed, and CG scheme for the optimization of the Rayleigh quotient has been proven a very attractive and promising technique for large sparse eigenproblems for small extreme eigenvalues. As in the case of a system of linear equations, successful application of the CG scheme to eigenproblems depends also upon the preconditioning techniques. A proper choice of the preconditioner significantly improves the convergence of the CG scheme. The idea underlying the present work is a parallel computation of the Multi-Color Block SSOR preconditioning for the CG optimization of the Rayleigh quotient together with deflation techniques. Multi-Coloring is a simple technique to obatin the parallelism of order n, where n is the dimension of the matrix. Block SSOR is a symmetric preconditioner which is expected to minimize the interprocessor communication due to the blocking. We implemented the results on the CRAY-T3E with 128 nodes. The MPI(Message Passing Interface) library was adopted for the interprocessor communications. The test problems were drawn from the discretizations of partial differential equations by finite difference methods.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.7
• /
• pp.265-274
• /
• 2018

The increasing use of curved beams in buildings, vehicles, ships, and aircraft has prompted studies directed toward the development of an accurate method for analyzing the dynamic behavior of such structures. The stability behavior of elastic curved beams has been the subject of a large number of investigations. Solutions of the relevant differential equations have been obtained traditionally using standard finite difference or finite element methods. These techniques require a great deal of computer time as the number of discrete nodes becomes relatively large under the conditions of complex geometry and loading. One of the efficient procedures for the solution of partial differential equations is the method of differential quadrature. The differential quadrature method (DQM) has been applied to a large number of cases to overcome the difficulties of the complex algorithms of programming for the computer, as well as the excessive use of storage due to the conditions of complex geometry and loading. The in-plane buckling of curved beams considering the extensibility of the arch axis was analyzed under uniformly distributed radial loads using the DQM. The critical loads were calculated for the member with various parameter ratios, boundary conditions, and opening angles. The results were compared with the precise results by other methods for cases, in which they were available. The DQM, using only a limited number of grid points, provided results that agreed very well (less than 0.3%) with the exact ones. New results according to diverse variations were obtained, showing the important roles in the buckling behavior of curved beams, and can be used in comparisons with other numerical solutions or with experimental test data.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.20 no.5
• /
• pp.612-620
• /
• 2019

The increasing use of curved beams in buildings, vehicles, ships, and aircraft has results in considerable effort being directed toward developing an accurate method for analyzing the dynamic behavior of such structures. The stability behavior of elastic curved beams has been the subject of a large number of investigations. Solutions of the relevant differential equations have traditionally been obtained by the standard finite difference. These techniques require a great deal of computer time as the number of discrete nodes becomes relatively large under conditions of complex geometry and loading. One of the efficient procedures for the solution of partial differential equations is the method of differential quadrature. The differential quadrature method(DQM) has been applied to a large number of cases to overcome the difficulties of the complex algorithms of programming for the computer, as well as excessive use of storage due to conditions of complex geometry and loading. In this study, the in-plane extensional vibration for asymmetric curved beams with linearly varying cross-section is analyzed using the DQM. Fundamental frequency parameters are calculated for the member with various parameter ratios, boundary conditions, and opening angles. The results are compared with the result by other methods for cases in which they are available. According to the analysis of the solutions, the DQM, used only a limited number of grid points, gives results which agree very well with the exact ones.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.13 no.3
• /
• pp.174-180
• /
• 2010

The inhomogeneous Helmholtz equation is introduced for variable water depth and potential function and separation of variables are introduced for the derivation. Only harmonic wave motions are considered. The governing equation composed of the potential function for irrotational flow is directly applied to the still water level, and the inhomogeneous Helmholtz equation for variable water depth is obtained. By introducing the wave amplitude and wave phase gradient the governing equation with complex potential function is transformed into two equations of real variables. The transformed equations are the first and second-order ordinary differential equations, respectively, and can be solved in a forward marching manner when proper boundary values are supplied, i.e. the wave amplitude, the wave amplitude gradient, and the wave phase gradient at a side boundary. Simple spatially-centered finite difference numerical schemes are adopted to solve the present set of equations. The equation set is applied to two test cases, Booij’ inclined plane slope profile, and Bragg’ wavy bed profile. The present equations set is satisfactorily verified against other theories including the full linear equation, Massel's modified mild-slope equation, and Berkhoff's mild-slope equation etc.

• Journal of Korean Society of Steel Construction
• /
• v.12 no.3 s.46
• /
• pp.303-315
• /
• 2000

The main object of this study is to analyze mechanical behaviors as anisotropic three-dimensional body under various static loads. This paper presents the applicability of the finite difference method to three dimensional problem of anisotropic body. The finite difference method as applied here is generalized to anisotropic three-dimensional problem of elastic body where the governing differential equations of equilibrium of such bodies are expressed in terms of the displacement u, v, and w in the coordinates axes x, y and z, care being taken to modify the finite difference expressions to satisfy the appropriate boundary conditions. By adopting a new three dimensional finite difference modelling including elimination of pivotal difference points in the case of free boundary condition, the three dimensional problem of anisotropic body was successfully completed. Several numerical results show quick convergence and numerical validity of finite difference technique in three dimensional problem.

• Journal of applied mathematics & informatics
• /
• v.32 no.1_2
• /
• pp.185-201
• /
• 2014

In this paper, we construct a numerical method to solve singularly perturbed one-dimensional parabolic convection-diffusion problems. We use Euler method with uniform step size for temporal discretization and exponential-spline scheme on spatial uniform mesh of Shishkin type for full discretization. We show that the resulting method is uniformly convergent with respect to diffusion parameter. An extensive amount of analysis has been carried out to prove the uniform convergence with respect to the singular perturbation parameter. The obtained numerical results show that the method is efficient, stable and reliable for solving convection-diffusion problem accurately even involving diffusion parameter.

• Journal of Institute of Control, Robotics and Systems
• /
• v.5 no.4
• /
• pp.481-488
• /
• 1999

A thermal model based on the chamber geometry of the industry-standard AST SHS200MA rapid thermal processing system has been developed for the study of thermal uniformity and process repeatability thermal model combines radiation energy transfer directly from the tungsten-halogen lamps and the steady-state thermal conducting equations. Because of the difficulties of solving partial differential equation, calculation of wafer temperature was performed by using finite-difference approximation. The proposed thermal model was verified via titanium silicidation experiments. As a result, we can conclude that the thermal model show good estimation of wafer surface temperature distribution.

• 전기의세계
• /
• v.21 no.3
• /
• pp.48-52
• /
• 1972

The application of pontryagin's Maximum Principle to the optimal control eventually leads to the problem of solving the two point boundary value problem. Most of problems have been related to their own special factors, therfore it is very hard to recommend the best method of deriving their optimal solution among various methods, such as iterative Runge Kutta, analog computer, gradient method, finite difference and successive approximation by piece-wise linearization. The gradient method has been applied to the optimal control of two point boundary value problem in the power systems. The most important thing is to set up some objective function of which the initial value is the function of terminal point. The next procedure is to find out any global minimum value from the objective function which is approaching the zero by means of gradient projection. The algorithm required for this approach in the relevant differential equations by use of the Runge Kutta Method for the computation has been established. The usefulness of this approach is also verified by solving some examples in the paper.