• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.6
• /
• pp.1465-1477
• /
• 1993

Powder Injection Molding(PM) is an advanced and complicated technology for manufacturing ceramic or metal products making use of a conventional injection molding process, which is generally used for plastic products. Among many technologies involved in the successful PIM, injection molding process is one of the key steps to form a desired shape out of powder/binder mixtures. Thus, it is of great importance to have a numerical tool to predict the powder injection molding filling process. In this regard, a finite element analysis system has been developed for numerical simulations of filling process of powder injection molding. Powder/polymer mixtures during the filling pro cess of injection molding can be rheologically characterized as Non-Newtonian fluids with a so called yield phenomena and have a peculiar feature of apparent slip phenomena on the wall boundaries surrounding mold cavity. Therefore, in the present study, a physical modeling of the filling process of powder/polymer mixtures was developed to take into account both the yield stress and slip phenomena and a finite element formulation was developed accordingly. The numerical analysis scheme for filling simulation is accomplished by combining a finite element method with control volume technique to simulate the movement of flow front and a finite difference method to calculate the temperature distribution. The present study presents the modeling, numerical scheme and some numerical analysis results showing the effect of the yield stress and slip phenomena.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.5
• /
• pp.1145-1158
• /
• 1994

A three-dimensional, finite difference, numerical model with free surface was developed on ${\sigma}$-coordinate. A semi-implicit numerical scheme in time has been adopted for computational efficiency. The scheme is essentially independent of the stringent stability criteria (CFL condition) for explicit schemes of external surface gravity wave. Implicit algorithm was applied for vertical shear stress, Coriolis force and pressure gradient terms. The reliability of the model with vertically variable grid system was checked by the comparison of simulation results with analytic solution of wind-driven currents in a one-dimensional channel. Sensitivity analysis of differencing parameters was carried out by applying the model to the calculation of wind-driven currents in a square lake.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.17 no.3
• /
• pp.20-27
• /
• 2013

Simple Iteration Method for calculating the natural frequency is presented in this paper. This method is simple but exact method of calculating natural frequencies corresponding to the modes of vibration of beams and tower structures with irregular cross sections and arbitrary boundary conditions. This method consists of determining the deflected mode shape of the member due to the inertia force under resonance condition. Finite difference method is used for this purpose. The influence of the $D_{22}$ stiffness on the natural frequency is rigorously investigated. In this paper, the influence of the loading sizes, different cross section on the natural frequency of vibration of some structural elements is presented. This method extends to two dimensional problems including advanced composite material structures.

• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.12 no.7
• /
• pp.1165-1172
• /
• 2008

In this paper, the scattering parameters of the metamaterial slab are obtained using the 2D FDTD(Finite-Difference Time-Domain) method. FDTD method is one of strongest electromagnetic numerical method which is widely used to analyze the metamaterial structure because of its simplicity. But it is very difficult to obtain frequency response of metamaterial itself because frequency dispersive model such as Lorentz, Drude model are used in FDTD. We used the well-known m-n-m cycle sine pulse to obtain the frequency response of the metamaterials. Comparisons between the wideband Gaussian input pulse and band-limited m-n-m cycle sine pulse are performed in this paper also. From the results, we concluded that the scattering parameters in frequency domain can be obtained using specific input pulse in FDTD even if the response has valid only for limited bandwidth.

• Water for future
• /
• v.21 no.2
• /
• pp.173-182
• /
• 1988

Characteristics of finite difference schemes for St. Venant equation were compared with two input cases. One is the monoclinal wave which has large friction slope without discontinuity and the other is the shock wave with discontinuity. For monoclinal wave, Keller Box scheme is the best in terms of accuracy, efficiency and stability when two parameters were selected with a rele : $0.5{\leq}{\theta}{\leq}1.0$, ${\theta}+{\psi}$=1, But for shock wave only the Preissmann type of parameter ${\psi}$(=0.5) showed stable results. Numerical experiments of monoclinal wave showed that Lax-Wendroff and Richtmyer schemes were more stable than leap Frog and more accurate than Lax-Fredrich scheme. For shock wave Lax-Fredrich showed larger numerical dissipation than other explicit schemes and Leap Frog produced slower mass transport.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.6
• /
• pp.1382-1390
• /
• 1990

Accurate load-carrying capacities and moment distributions of thin circular plates are obtained for clamped or simply-supported boundary condition under various concentrated circular loadings. The material is regarded as perfectly-plastic based on an arbitrary yield function such as the Tresca yield function, the Johansen yield function, and the farmily of .betha.-norms which possesses the von Mises yield function and the Frobenius norm. To obtain the lower bound solutions, a maximization formulation is derived and implemented for efficient numerical calculation with a finite difference method and the modified Newton's method. The numerical results demonstrate plastic collapse behavior of circular plates and provide their design criteria.

• Journal of Korea Water Resources Association
• /
• v.37 no.11
• /
• pp.889-896
• /
• 2004

The fractal advection-diffusion equation (ADE) is a generalization of the classical AdE in which the second-order derivative is replaced with a fractal order derivative. While the fractal ADE have been analyzed with a stochastic process In the Fourier and Laplace space so far, in this study a fractal ADE for describing solute transport in rivers is derived with a finite difference scheme in the real space. This derivation with a finite difference scheme gives the hint how the fractal derivative order and fractal diffusion coefficient can be estimated physically In contrast to the classical ADE, the fractal ADE is expected to be able to provide solutions that resemble the highly skewed and heavy-tailed time-concentration distribution curves of contaminant plumes observed in rivers.

• Journal of Soil and Groundwater Environment
• /
• v.7 no.4
• /
• pp.10-16
• /
• 2002

The multi-aquifer system separated by semipervious leaky beds was analyzed. The finite difference scheme of the Crank-Nicolson method is applied to obtain the solution for this system. The solution of this scheme was compared with the analytical solution for two-layer aquifer systems with one-dimensional steady state. The results showed a good agreement between analytical and numerical solution for two-layer aquifer systems. So, the numerical scheme can be extended to multi-aquifer system. When the pumping is tried for single or multi aquifer, the computation of the groundwater heads was possible for each aquifer in the multi-aquifer with two-dimensional system. So, it might be helpful for the effective groundwater management.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.23 no.4
• /
• pp.846-856
• /
• 1998

In this paper, we have applied Berenger's perfectly matched layer(PML) absorbing boundary conditition(ABC) with three dimension at the propagation direction of field, side boundary and upper boundary in which evanescent and radiation field is exist. Even though the same computationaal domain is used, PML is superior to other ABC in the absorbing ability that is obtained reflective coefficient about $10^{-4}$ unit for the propagation wave, but PML is ineffective in absorbing evanescent field. Also we have compared dispersive characteristics of the microstrip with Kobayashi's emprical method, and obtaind very similar result, but 0.8% error is generated at the dc or near.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.28 no.2A
• /
• pp.63-69
• /
• 2003

In this paper, the scattered field from a perfectly conducting fractal surface by Finite-Difference Time-Domain(FDTD) method was computed. A one-dimensional fractal surface was generated by using the fractional Brownian motion model. Back scattering coefficients are calculated with different values of the spectral parameter(S0), fractal dimension(D) which determine characteristics of the fractal surface. The number of surface realization for the computed field, the point number, and the width of surface realization are set to be 80, 1024, 16λ, respectively. In order to verify the computed results these results are compared with those of small perturbation methods, which show good agreement between them.