• Journal of Korea Foundry Society
• /
• v.20 no.2
• /
• pp.129-140
• /
• 2000

This study aims at the phase-field modeling of the phase transformation in graphitization of the cast iron. As the first step, we constructed a phase-field model including the phases with negligible solubility. Under the dilute regular solution approximation, a simplified version of the phase-field model was obtained, which can be used for the phase transformation related with the stoichiometric phases. The results from the numerical calculation of the phase-field model was in good agreement with the exact analytic solution. The compositional shift due to Gibbs-Thomson effect can be reproduced within 0.5% error in the numerical calculation. The interface velocity, whereas, in numerical calculation of phase-field model appeared to be 15% larger than that from the analytic solution. This error is due to the shift of the interface position in phase-field model from the position with ${\phi}=0.5$.

• Journal of Mechanical Science and Technology
• /
• v.14 no.2
• /
• pp.251-258
• /
• 2000

Thermal deformation in an electronic package due to thermal strain mismatch is investigated. The warpage and the in-plane deformation of the package after encapsulation is analyzed using the laminated plate theory. An exact solution for the thermal deformation of an electronic package with circular shape is derived. Theoretical results are presented on the effects of the layer geometries and material properties on the thermal deformation. Several applications of the exact solution to electronic packaging product development are illustrated. The applications include lead on chip package, encapsulated chip on board and chip on substrate.

• Bulletin of the Korean Chemical Society
• /
• v.7 no.4
• /
• pp.271-275
• /
• 1986

The Smoluchowski equation with a step potential is solved in one-dimensional case and three-dimensional case with spherical symmetry. Exact analytic expressions for the solution and the remaining probability are obtained in one-dimensional case for the reflecting boundary condition and the long time behavior of the remaining probability is compared with the earlier work. In three-dimensional case, only the long time behavior is evaluated. More general case with the radiation boundary condition is also investigated and the results are shown to approach correct limits of the reflecting boundary condition.

• Journal of applied mathematics & informatics
• /
• v.40 no.3_4
• /
• pp.675-694
• /
• 2022

In this article, we investigate the solution of the inverse problem for one dimensional heat equation with Neumann and Robin boundary conditions, that is, we determine the temperature and source term with given initial and boundary conditions. Three radial basis functions(RBFs) have been used for numerical solution, and Homotopy perturbation method for analytic solution. Numerical solutions which are obtained by considering each of the three RBFs are compared to the exact solution. For appropriate value of shape parameter c, numerical solutions best approximates exact solutions. Furthermore, we have shown the impact of noisy data on the numerical solution of u and f.

• Nuclear Engineering and Technology
• /
• v.27 no.3
• /
• pp.374-384
• /
• 1995

The mathematical adjoint solution of the Analytic Function Expansion (AFEN) method is found by solving the transposed matrix equation of AFEN nodal equation with only minor modification to the forward solution code AFEN. The perturbation calculations are then performed to estimate the change of reactivity by using the mathematical adjoint The adjoint calculational scheme in this study does not require the knowledge of the physical adjoint or the eigenvalue of the forward equation. Using the adjoint solutions, the exact and first-order perturbation calculations are peformed for the well-known benchmark problems (i.e., IAEA-2D benchmark problem and EPRI-9R benchmark problem). The results show that the mathematical adjoint flux calculated in the code is the correct adjoint solution of the AFEN method.

• Coupled systems mechanics
• /
• v.11 no.2
• /
• pp.167-198
• /
• 2022

In this paper we deal with classical instability problems of heterogeneous Euler beam under conservative loading. It is chosen as the model problem to systematically present several possible solution methods from simplest deterministic to more complex stochastic approach, both of which that can handle more complex engineering problems. We first present classical analytic solution along with rigorous definition of the classical Euler buckling problem starting from homogeneous beam with either simplified linearized theory or the most general geometrically exact beam theory. We then present the numerical solution to this problem by using reduced model constructed by discrete approximation based upon the weak form of the instability problem featuring von Karman (virtual) strain combined with the finite element method. We explain how such numerical approach can easily be adapted to solving instability problems much more complex than classical Euler's beam and in particular for heterogeneous beam, where analytic solution is not readily available. We finally present the stochastic approach making use of the Duffing oscillator, as the corresponding reduced model for heterogeneous Euler's beam within the dynamics framework. We show that such an approach allows computing probability density function quantifying all possible solutions to this instability problem. We conclude that increased computational cost of the stochastic framework is more than compensated by its ability to take into account beam material heterogeneities described in terms of fast oscillating stochastic process, which is typical of time evolution of internal variables describing plasticity and damage.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.25 no.4A
• /
• pp.489-495
• /
• 2000

In this paper we propose an approximate method to estimate the cell loss probability(CLP) due to buffer overflow in ATM networks. The main idea is to relate the buffer capacity with the CLP target in explicit formula by using the approximate upper bound for the tail distribution of a queue. The significance of the proposition lies in the fact that we can obtain the expected CLP by using only the source traffic data represented by mean rate and its variance. To that purpose we consider the problem of estimating the cell loss measures form the statistical viewpoint such that the probability of cell loss due to buffer overflow does not exceed a target value. In obtaining the exact solution we use a typical matrix analytic method for GI/D/1B queue where B is the queue size. Finally, in order to investigate the accuracy of the result, we present both the approximate and exact results of the numerical computation and give some discussion.

• Proceedings of the KSR Conference
• /
• 2005.11a
• /
• pp.1165-1170
• /
• 2005

In order to perform the spatial buckling analysis of the curved beam element with nonsymmetric thin-walled cross section, exact static stiffness matrices are evaluated using equilibrium equations and force-deformation relations. Contrary to evaluation procedures of dynamic stiffness matrices, 14 displacement parameters are introduced when transforming the four order simultaneous differential equations to the first order differential equations and 2 displacement parameters among these displacements are integrated in advance. Thus non-homogeneous simultaneous differential equations are obtained with respect to the remaining 8 displacement parameters. For general solution of these equations, the method of undetermined parameters is applied and a generalized linear eigenvalue problem and a system of linear algebraic equations with complex matrices are solved with respect to 12 displacement parameters. Resultantly displacement functions are exactly derived and exact static stiffness matrices are determined using member force-displacement relations. The buckling loads are evaluated and compared with analytic solutions or results by ABAQUS's shell element.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.34 no.2
• /
• pp.99-103
• /
• 2010

A new analytic solution was derived for the diffusion into or from an immobile zone of a rectangular 2-D pore. For a long time, the new solution converges to a traditional mobile-immobile zone (MIM) model, but only if the latter is used with an apparent initial concentration that is smaller by almost 20% than the true one. This is the tradeoff for using a simple MIM model instead of an exact model based on the diffusion equation. The mass-transfer coefficient was found to be constant for a sufficiently long time; it was proportional to the molecular diffusion and inversely proportional to the square of the pore depth. The mass-transfer coefficient was time-dependent for a sufficiently short time and may be significantly larger than its asymptotic value.

• 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.