• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.22 no.10
• /
• pp.1005-1011
• /
• 2011

In this paper, we present a marching-on-in-degree(MOD) finite difference method(FDM) based on the Helmholtz wave equation for analyzing transient electromagnetic responses in a general dispersive media. The two issues related to the finite difference approximation of the time derivatives and the time consuming convolution operations are handled analytically using the properties of the Laguerre functions. The basic idea here is that we fit the transient nature of the fields, the flux densities, the permittivity with a finite sum of orthogonal Laguerre functions. Through this novel approach, not only the time variable can be decoupled analytically from the temporal variations but also the final computational form of the equations is transformed from finite difference time-domain(FDTD) to a finite difference formulation through a Galerkin testing. Representative numerical examples are presented for transient wave propagation in general Debye, Drude, and Lorentz dispersive medium.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.22 no.3
• /
• pp.201-215
• /
• 2018

In this paper we consider a class of singularly perturbed system of delay differential equations of convection diffusion type with integral boundary conditions. A finite difference scheme on an appropriate piecewise Shishkin type mesh is suggested to solve the problem. We prove that the method is of almost first order convergent. An error estimate is derived in the discrete maximum norm. Numerical experiments support our theoretical results.

• Journal of Magnetics
• /
• v.22 no.1
• /
• pp.104-108
• /
• 2017

Induction motors have wide applicability in many fields, both in industrial sectors and households, for their advantages of a high efficiency and robust structure. The introduction of power-source-containing harmonics into the induction motor winding lowers its efficiency and increases its temperature, greatly affecting its operation characteristics. In this study, we performed an electromagnetic field analysis using the time-difference finite-element method with the purpose of analyzing the steady-state current characteristics of an induction motor. Additionally, we calculated the steady-state current with a method combining an electromagnetic field equation and a circuit equation. In the electromagnetic field analysis, the nonlinearity was taken into account using the Newton-Raphson method, and a backward time-difference method was employed for the time derivative term. Then, we compared the steady-state current of the induction motor obtained by calculation with the experimentally measured values, thus validating the proposed algorithm. Furthermore, we analyzed the impacts of the shape and material of the rotor conductor bar of the induction motor on the steady-state current of the main winding.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.57-62
• /
• 2008

Finite difference method using not general SSG(standard staggered grid) but RSG(rotated staggered grid) was applied to simulation of elastic wave propagation. Special free surface boundary condition such as imaging method is needed in finite difference method using SSG in elastic wave propagation but free surface boundary condition in finite difference method using RSG is easily solved with adding air layer. Recently PML(Perfectly Matched layer) is widely used to eliminate artificial reflection waves from finite boundary because of its' greate efficiency. Absorbing ability of CPML(convolutional Perfectly Matched Layer) that is more efficient than that of PML was applied to FDM using RSG in this study. The results of CPML eliminated artificial boundary waves very effectively in FDM using RSG in being compared with that of Cerjan's absorbing method.

• 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 Mechanical Science and Technology
• /
• v.19 no.11
• /
• pp.2032-2039
• /
• 2005

The lattice Boltzman method (LBM) and the finite difference-based lattice Boltzmann method (FDLBM) are quite recent approaches for simulating fluid flow, which have been proven as valid and efficient tools in a variety of complex flow problems. They are considered attractive alternatives to conventional finite-difference schemes because they recover the Navier-Stokes equations and are computationally more stable, and easily parallelizable. However, most models of the LBM or FDLBM are for incompressible fluids because of the simplicity of the structure of the model. Although some models for compressible thermal fluids have been introduced, these models are for monatomic gases, and suffer from the instability in calculations. A lattice BGK model based on a finite difference scheme with an internal degree of freedom is employed and it is shown that a diatomic gas such as air is successfully simulated. In this research we present a 2-dimensional edge tone to predict the frequency characteristics of discrete oscillations of a jet-edge feedback cycle by the FDLBM in which any specific heat ratio $\gamma$ can be chosen freely. The jet is chosen long enough in order to guarantee the parabolic velocity profile of a jet at the outlet, and the edge is of an angle of $\alpha$=23$^{o}$. At a stand-off distance w, the edge is inserted along the centerline of the jet, and a sinuous instability wave with real frequency is assumed to be created in the vicinity of the nozzle exit and to propagate towards the downstream. We have succeeded in capturing very small pressure fluctuations resulting from periodic oscillation of the jet around the edge.

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

A parameter uniform numerical scheme is proposed for solving singularly perturbed parabolic partial differential-difference convection-diffusion equations with a small delay and advance parameters in reaction terms and spatial variable. Taylor's series expansion is applied to approximate problems with the delay and advance terms. The resulting singularly perturbed parabolic convection-diffusion equation is solved by utilizing the implicit Euler method for the temporal discretization and finite difference method for the spatial discretization on a uniform mesh. The proposed numerical scheme is shown to be an ε-uniformly convergent accurate of the first order in time and second-order in space directions. The efficiency of the scheme is proved by some numerical experiments and by comparing the results with other results. It has been found that the proposed numerical scheme gives a more accurate approximate solution than some available numerical methods in the literature.

• International Journal of Highway Engineering
• /
• v.17 no.5
• /
• pp.1-10
• /
• 2015

PURPOSES : A finite difference model considering snow melting process on porous asphalt pavement was derived on the basis of heat transfer and mass transfer theories. The derived model can be applied to predict the region where black-ice develops, as well as to predict temperature profile of pavement systems where a de-icing system is installed. In addition, the model can be used to determined the minimum energy required to melt the ice formed on the pavement. METHODS : The snow on the porous asphalt pavement, whose porosity must be considered in thermal analysis, is divided into several layers such as dry snow layer, saturated snow layer, water+pavement surface, pavement surface, and sublayer. The mass balance and heat balance equations are derived to describe conductive, convective, radiative, and latent transfer of heat and mass in each layer. The finite differential method is used to implement the derived equations, boundary conditions, and the testing method to determine the thermal properties are suggested for each layer. RESULTS: The finite differential equations that describe the icing and deicing on pavements are derived, and we have presented them in our work. The framework to develop a temperature-forecasting model is successfully created. CONCLUSIONS : We conclude by successfully creating framework for the finite difference model based on the heat and mass transfer theories. To complete implementation, laboratory tests required to be performed.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.15 no.3
• /
• pp.190-195
• /
• 2003

For many coastal problems, such as wave transformation, tidal circulation, sediment transports and diffusion phenomena, we resort to numerical techniques. The representative numerical techniques are the method of finite differences and finite elements. The approximate algebraic equations, referred to as finite difference equations(FDEs), are subsequently solved at discrete grid points within the domain of interests. Therefore, a set of grid points within the domain, as well as the boundaries of the domain, must be specified. The generation of grids for FDEs, with uniform spacing, is very simple compared to that of finite elements. However, within a very complex domain, there are few grid generation tools we can use conveniently. Unfortunately, most of the commercial grid generation programs are developed only for finite element method. In this paper, grid generation method using digitizer, with uniform rectangular spacing, are introduced in detail. Didger and Surfer programs by Golden Software are necessary to produce comparatively accurate and simple depth data.

• Journal of KSNVE
• /
• v.8 no.2
• /
• pp.267-273
• /
• 1998

This study aims at performing sensitivity analysis of piezoelectric smart structure for minimizing radiated noise from the structure, The structure consists of a flat plate on which disk shaped piezoelectric actuator is mounted, and finite element modeling is used for the structure. The finite element modeling uses a combination of three dimensional piezoelectric, flat shell and transition elements so thus it can take into account the coupling effects of the piezoelectric device precisely and it can also reduce the degrees of freedom of the finite element model. Electric potential on the piezoelectric actuator is taken as a design variable and total radiated power of the structure is chosen as an objective function. The objective function can be represented as Rayleigh's integral equation and is a function of normal displacements of the structure. For the convenience of computation, all degrees of freedom of the finite element equation is condensed out except the normal displacements of the structure. To perform the design sensitivity analysis, the derivative of the objective function with respect to the normal displacements is found, and the derivative of the norma displacements with respect to the design variable is calculated from the finite element equation by using so called the adjoint variable method. The analysis results are compared with those of the finite difference method, and shows a good agreement. This sensitivity analysis is faster and more accurate than the finite difference method. Once the sensitivity analysis program is used for gradient-based optimizations, one could achieve a better convergence rate than non-derivative methods for optimal design of piezoelectric smart structures.