• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 1999.11a
• /
• pp.146-153
• /
• 1999

A numerical dynamic simulation is necessary to investigate the capacity of the HDD. The slider surface become more and more complicated to make the magnetized area smaller and readback signal stronger. So a numerical dynamic simulation must be preceded to develop a new slider in HDD. The dynamic simulations of air-lubricated slider bearing have been peformed using FIFD(Factored Implicit Finite Difference) method. The governing equation, Reynolds equation Is modified with Fukui and Kaneko model(FK model) which includes the first and the second-order slip. The equations of motion for the slider bearing are solved simultaneously with the modified Reynolds equation for the case of three degrees of freedom. The slider transient response for disk step bump and slider impulse force is given for various case and for iteration algorithm and new algorithm.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.48-54
• /
• 1998

Three-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method is used for time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetrahedra, prisms, pyramids, hexahedra, or mixed-element grid. The numerical efficiency and accuracy of the present method is critically evaluated for several example problems.

• Journal of Korea Foundry Society
• /
• v.13 no.4
• /
• pp.323-332
• /
• 1993

An implicit finite difference formulation with three methods of latent heat treatment, such as equivalent specific heat method, temperature recovery method and enthalpy method, was applied to solidification analysis. The Neumann problem was solved to compare the numerical results with the exact solution. The implicit solutions with the equivalent specific heat method and the temperature recovery method were comparatively consistent with the Neumann exact solution for smaller time steps, but its error increased with increasing time step, especially in predicting the solidification beginning time. Although the computing time to solve energy equation using temperature recovery method was shorter than using enthalpy method, the method of releasing latent heat is not realistic and causes error. The implicit formulation of phase change problem requires enthalpy method to treat the release of latent heat reasonably. We have modified the enthalpy formulation in such a way that the enthalpy gradient term is not needed, and as a result of this modification, the computation stability and the computing time were improved.

• Journal of electromagnetic engineering and science
• /
• v.9 no.1
• /
• pp.12-16
• /
• 2009

This paper presents a new iterative locally one-dimensional [mite-difference time-domain(LOD-FDTD) method that has a simpler formula than the original iterative LOD-FDTD formula[l]. There are fewer arithmetic operations than in the original LOD-FDTD scheme. This leads to a reduction of CPU time compared to the original LOD-FDTD method while the new method exhibits the same numerical accuracy as the iterative ADI-FDTD scheme. The number of arithmetic operations shows that the efficiency of this method has been improved approximately 20 % over the original iterative LOD-FDTD method.

• Transactions of Materials Processing
• /
• v.9 no.4
• /
• pp.395-403
• /
• 2000

Formation of Shear Band under the adiabatic condition is widely observed In the engineering materials during rapidly forming process lot a thermally rate-dependent material. The shear band stems from evolution of a narrow region in which an intensive plastic flow occurs. The shear band often plays a role of a precursor of the ductile fracture during a forming process. The objective of this study is to investigate the localization behavior using numerical method. In this work, the implicit finite difference scheme is employed due to the ease of convergence and the numerical stability It is noted that physical and mechanical properties of materials determine how the shear band is formed and then localized. Material properties can be characterized with inertia number dissipation number and diffusion number. It is observed that the dimensionless numbers effect on localization. Using a parametric study, comparison was made between CRS-1018 steel with WHA (tungsten heavy alloy). The deformation behavior of material in this study include an isotropic hardening as well as thermal softening. Moreover, this study suggests that a kinematic hardening constitutive relation be required to predict a more accurate strain level at a shear band.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.3
• /
• pp.197-207
• /
• 2013

The Cahn-Hilliard equation was proposed as a phenomenological model for describing the process of phase separation of a binary alloy. The equation has been applied to many physical applications such as amorphological instability caused by elastic non-equilibrium, image inpainting, two- and three-phase fluid flow, phase separation, flow visualization and the formation of the quantum dots. To solve the Cahn-Hillard equation, many numerical methods have been proposed such as the explicit Euler's, the implicit Euler's, the Crank-Nicolson, the semi-implicit Euler's, the linearly stabilized splitting and the non-linearly stabilized splitting schemes. In this paper, we investigate each scheme in finite-difference schemes by comparing their performances, especially stability and efficiency. Except the explicit Euler's method, we use the fast solver which is called a multigrid method. Our numerical investigation shows that the linearly stabilized stabilized splitting scheme is not unconditionally gradient stable in time unlike the known result. And the Crank-Nicolson scheme is accurate but unstable in time, whereas the non-linearly stabilized splitting scheme has advantage over other schemes on the time step restriction.

• Journal of computational fluids engineering
• /
• v.3 no.2
• /
• pp.17-26
• /
• 1998

The three-dimensional incompressible Navier-Stokes equations have been solved by a node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method with Jacobi matrix solver is used for the time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetragedra, prisms, pyramids, hexahedra, or mixed-element grid. Inviscid bump flow is solved to check the accuracy of high order convective flux discretisation. And viscous flows around a circular cylinder and a sphere are studied to show the efficiency and accuracy of the solver.

• Journal of applied mathematics & informatics
• /
• v.40 no.1_2
• /
• pp.25-45
• /
• 2022

In this study, the non-standard finite difference method for the numerical solution of singularly perturbed parabolic reaction-diffusion subject to Robin boundary conditions has presented. To discretize temporal and spatial variables, we use the implicit Euler and non-standard finite difference method on a uniform mesh, respectively. We proved that the proposed scheme shows uniform convergence in time with first-order and in space with second-order irrespective of the perturbation parameter. We compute three numerical examples to confirm the theoretical findings.

• Journal of Korea Water Resources Association
• /
• v.45 no.2
• /
• pp.179-189
• /
• 2012

To simulate the transport of pollutant, a numeric model for the advection-diffusion equation with the decay term is developed. This is finite-difference model using the implicit method (with the weight factor ${\alpha}$) and Gauss-Seidel SOR(successive over-relaxation). This model is compared to the analytical solutions (of simpler dimensional or boundary conditions), and in the condition of Peclet number < 5~20, the result shows stable condition, and Crank-Nicolson method (${\alpha}$=0.5) shows the more accurate results than fully-implicit method (${\alpha}$=1). The mass of advection, diffusion and decay is calculated and the error of mass balance is less than 3%. This model can evaluate the 3-D concentrations of the advection-diffusion and decay problems, but this model uses only the finite-difference method with the fixd grid system, so it can be effectively used in the problems with small Peclet numbers like the pollutant transport in groundwater.

• Kyungpook Mathematical Journal
• /
• v.60 no.3
• /
• pp.629-645
• /
• 2020

A uniformly convergent numerical method is developed for solving singularly perturbed 1-D parabolic convection-diffusion problems. The developed method applies a non-standard finite difference method for the spatial derivative discretization and uses the implicit Runge-Kutta method for the semi-discrete scheme. The convergence of the method is analyzed, and it is shown to be first order convergent. To validate the applicability of the proposed method two model examples are considered and solved for different perturbation parameters and mesh sizes. The numerical and experimental results agree well with the theoretical findings.