• 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 the Earthquake Engineering Society of Korea
• /
• v.26 no.6
• /
• pp.227-235
• /
• 2022

For a member model in nonlinear structural analysis, a lumped plastic model that idealizes its flexural bending, shear, and axial behaviors by springs with the nonlinear hysteretic model is widely adopted because of its simplicity and transparency compared to the other rigorous finite element methods. On the other hand, a challenging task in its numerical solution is to satisfy the equilibrium condition between nonlinear flexural bending and shear springs connected in series. Since the local forces between flexural and shear springs are not balanced when one or both springs experience stiffness changes (e.g., cracking, yielding, and unloading), the additional unbalanced force due to overshooting or undershooting each spring force is also generated. This paper introduces an iterative scheme for numerical solutions satisfying the equilibrium conditions between flexural bending and shear springs. The effect of equilibrium iteration on analysis results is shown by comparing the results obtained from the proposed method to those from the conventional scheme, where the equilibrium condition is not perfectly satisfied.

• Journal of electromagnetic engineering and science
• /
• v.10 no.2
• /
• pp.45-49
• /
• 2010

This paper presents an analytical formulation of the numerical reflection coefficient of the ID-FDTD scheme at the planar dielectric boundary for a TM wave incidence. The reflection coefficient is formulated in an approximate manner, and the accuracy of this method is numerically verified. The effective dielectric constant for a grid on the interface is obtained, and then reduced to that of the Yee scheme for a small cell size.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2006.05a
• /
• pp.253-256
• /
• 2006

A numerical study for spray flow of fuel and oxidizer droplets in the combustion chamber has been conducted prior to the analysis of spray combustion of the liquid rocket engine. As the spray combustion model, DSF model and Euler-Lagrange scheme have been used. While the coupling effects of the droplets between gas phase and evaporated vapor have been calculated using PSIC model, SIMPLER algorithm and QUICK scheme have been used as numerical schemes. As the results, the calculations have shown velocity and temperature distribution in combustion chamber as well as mole fraction of fuel and oxidizer.

• Bulletin of the Korean Mathematical Society
• /
• v.50 no.5
• /
• pp.1737-1751
• /
• 2013

An efficient and stable point collocation scheme based on a meshfree method is studied for the stationary incompressible Navier-Stokes equations. We describe the diffuse derivatives associated with the moving least square method. Using these diffuse derivatives, we propose a point collocation method to fit in solving the Navier-Stokes equations which improves the stability of the direct point collocation scheme. The convergence of the numerical solution is investigated from numerical examples. The driven cavity ow and the backward facing step ow are implemented for the reliability of the scheme. Also, the viscous ow on complicated geometry is successfully calculated such as the ow past a circular cylinder in duct.

• 유체기계공업학회:학술대회논문집
• /
• 2001.11a
• /
• pp.23-29
• /
• 2001

Three-dimensional flow through a tubular centrifugal fan with airfoil type blades are analyzed using CFX-TASCflow. Standard k - $\epsilon$ model and k - $\omega$ model are used as turbulence closures. The numerical schemes for convetion terms, i.e., Upwind Differencing Scheme(UDS), Mass Weighted Skewed Upstream Differencing Scheme(MWS), Linear Profile Skewed Upstream Differencing Scheme(LPS), and Modified Linear Profile Skewed Upstream Differencing Scheme(MLPS) are also tested. And, the performance of these schemes coupled with two turbulence models are evaluated. Computational static pressure distributions are compared with experimental data obtained in this work.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.27-30
• /
• 2002

As computer capacity has been progressed continuously, the studies of the flow characteristics have been performing by the numerical methods actively. In this study, 3-dimensional unsteady incompressible Wavier-Stokes equation was solved by numerical method using the fractional step method with the fourth order compact pade' scheme to achieve high accuracy To validate the present code and algorithm, 3D flow-field around a cylinder was simulated. The drag coefficient and lift coefficient were computed and, then, compared with experiment. The present code will be tailored to LES simulation for more accurate turbulent flow analysis.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.14 no.3
• /
• pp.532-537
• /
• 2011

This paper describes the development of SPH(Smooth Particle Hydrodynamics) scheme and integration into the multi-material shock physics code(ExLO) for the purpose of the application to the extreme large deformation problems. SPH numerical scheme has been extended into the fluid dynamics and the high-speed impact events, such as space structure protection against space debris and meteorite catering. Like other hydrocodes, SPH scheme also solves the conservation equations with the constitutive equation including equation of state. The benchmark problem, Taylor-Impact test, was simulated and the predictions show good agreements with both the published numerical data and experimental data. Currently, the contact treatment between materials is under development.

• Journal of applied mathematics & informatics
• /
• v.39 no.5_6
• /
• pp.623-641
• /
• 2021

In this paper, numerical treatment of singularly perturbed parabolic delay differential equations is considered. The considered problem have small delay on the spatial variable of the reaction term. To treat the delay term, Taylor series approximation is applied. The resulting singularly perturbed parabolic PDEs is solved using Crank Nicolson method in temporal direction with non-standard finite difference method in spatial direction. A detail stability and convergence analysis of the scheme is given. We proved the uniform convergence of the scheme with order of convergence O(N^-1 + (∆t)²), where N is the number of mesh points in spatial discretization and ∆t is mesh length in temporal discretization. Two test examples are used to validate the theoretical results of the scheme.

• The Journal of the Acoustical Society of Korea
• /
• v.25 no.1E
• /
• pp.32-35
• /
• 2006

The screech tone of underexpanded jet is numerically calculated without any specific modeling for the screech tone itself. Fourth-order optimized compact scheme and fourth-order Runge-Kutta method are used to solve the 2D axisymmetric Euler equation. Adaptive nonlinear artificial dissipation model and generalized characteristic boundary condition are also used. The screech tone, generated by a closed loop between instability waves and quasi-periodic shock cells at the near field, is reasonably analyzed with present numerical methods for the underexpanded jet having Mach number 1.13. First of all, the centerline mean pressure distribution is calculated and compared with experimental and other numerical results. The instantaneous density contour plot shows Mach waves due to mixing layer convecting supersonically, which propagate downstream. The pressure signal and its Fourier transform at upstream and downstream shows the directivity pattern of screech tone very clearly. Most of all, we can simulate the axisymmetric mode change of screech tone very precisely with present method. It can be concluded that the basic phenomenon of screech tone including the frequency can be calculated by using high-order and high-resolution schemes without any specific numerical modeling for screech tone feedback loop.