• The Transactions of the Korean Institute of Electrical Engineers C
• /
• v.51 no.11
• /
• pp.559-566
• /
• 2002

In this paper, we present an accurate and stable method for the solution of the transient electromagnetic response from the conducting wire structures using the time domain integral equation. By using an implicit scheme with the central finite difference approximation for the time domain electric field integral equation, we obtain the transient response from a wire scatterer illuminated by a plane wave and a conducting wire antenna with an impressed voltage source. Also, we consider a wire above a 3-dimensional conducting object. Numerical results are presented, which show the validity of the presented methodology, and compared with a conventional method using backward finite difference approximation.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.5A
• /
• pp.457-465
• /
• 2009

This study presents a moving least squares (MLS) finite difference method for solving elastic crack problems with stress singularity at the crack tip. Near-tip functions are intrinsically employed in the MLS approximation to model near-tip field inducing singularity in stress field. employment of the functions does not lose the merit of the MLS Taylor polynomial approximation which approximates the derivatives of a function without actual differentiating process. In the formulation of crack problem, computational efficiency is considerably improved by taking the strong formulation instead of weak formulation involving time consuming numerical quadrature Difference equations are constructed on the nodes distributed in computational domain. Numerical experiments for crack problems show that the intrinsically enriched MLS finite difference method can sharply capture the singular behavior of near-tip stress and accurately evaluate stress intensity factors.

• ETRI Journal
• /
• v.45 no.4
• /
• pp.704-712
• /
• 2023

This paper presents a flexible finite-difference technique for analyzing the nonuniform guiding structures. Because the voltage and current variations along the nonuniform structure differ for each segment, this work considers the adaptable discretization steps. This technique increases the accuracy of the final response. Moreover, by applying the singular value decomposition and discarding the nonprincipal singular values, an optimal lower rank approximation of the discretization matrix is obtained. The computational cost of the introduced method is significantly reduced using the optimal discretization matrix. Also, the proposed method can be extended to the nonuniform waveguides. The technique is verified by analyzing several practical transmission lines and waveguides with nonuniform profiles.

• Tribology and Lubricants
• /
• v.5 no.2
• /
• pp.48-54
• /
• 1989

The effects of angular misalignment, coning and the temperature difference between the primary seal ring and the seal seat on the pressure distribution in mechanical face seals are analyzed. The modified Reynolds equation for the temperature dependent viscosity was solved by a finite difference approximation and Gauss-Seidel method. It is shown that the amplitude of pressure is significantly affected by the misalignment of the seals and a large temperature difference between the rotor and the stator.

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

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.10 no.2
• /
• pp.180-186
• /
• 1999

Conformal methods such as staircase approximation and effective medium theory have been used to model arbitrary dielectric surfaces in Cartesian FDTD(finite-difference time-domain) grids. However, the staircase approximation doesn't guarantee accurate results and the effective medium theory cannot be applied for frequency dispersive medium. In this paper, a simple conformal method is introduced for analyzing arbitrary dielectric surfaces. These arbitrarily shaped dielectric surfaces make nonuniformly filled FDTD cells. E(H)-field in a nonuniformly filled FDTD cell is calculated by considering the cell as the combination of two kinds of uniformly filled cells whose material boundaries are normal and parallel to the E(H)-field in the cell. TM scattering from a rotated 2-D dielectric and a ferrite square cylinder is analyzed to show the accuracy of the proposed method.

• Journal of applied mathematics & informatics
• /
• v.35 no.3_4
• /
• pp.277-302
• /
• 2017

In this paper, we have proposed a numerical method for Singularly Perturbed Boundary Value Problems (SPBVPs) of reaction-diffusion type of third order Ordinary Differential Equations (ODEs). The SPBVP is reduced into a weakly coupled system of one first order and one second order ODEs, one without the parameter and the other with the parameter ${\varepsilon}$ multiplying the highest derivative subject to suitable initial and boundary conditions, respectively. The numerical method combines boundary value technique, asymptotic expansion approximation, shooting method and finite difference scheme. The weakly coupled system is decoupled by replacing one of the unknowns by its zero-order asymptotic expansion. Finally the present numerical method is applied to the decoupled system. In order to get a numerical solution for the derivative of the solution, the domain is divided into three regions namely two inner regions and one outer region. The Shooting method is applied to two inner regions whereas for the outer region, standard finite difference (FD) scheme is applied. Necessary error estimates are derived for the method. Computational efficiency and accuracy are verified through numerical examples. The method is easy to implement and suitable for parallel computing. The main advantage of this method is that due to decoupling the system, the computation time is very much reduced.

• Journal of the korean Society of Automotive Engineers
• /
• v.14 no.2
• /
• pp.54-64
• /
• 1992

During combustion process in a diesel engine radiation heat transfer is the same order of magnitude as the convection heat transfer. An approximation of heat and momentum source distributions is applied at a level consistent with those used in modelling the soot distribution and the turbulence instead of modelling the fuel spray and the chemical kinetics. This paper illustrates a use of the third order spherical harmonics approximation to the radiative transfer equation and delta-Eddington approximation to the scattering phase function for droplets in the flow. Results are obtained numerically by a time marching finite difference scheme. This study aims to compare heat transfer with convection heat transfer and to investigate the importance of scattering by fuel droplets and of accounting for spatial variations in the extinction coefficient on the radiative heat flux distributions at the walls of a disc shaped diesel engine.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.181-182
• /
• 2003

A preconditioned numerical method for gas-liquid to-phase flow is applied to solve cavitating flow. The present method employs a density based finite-difference method of dual time-stepping integration procedure and Roe's flux difference splitting approximation with MUSCL-TVD scheme. A homogeneous equilibrium cavitation model is used. The method permits simple treatment of the whole gas-liquid two-phase flow field including wave propagation, large density changes and incompressible flow characteristics at low Mach number. By this method, two-dimensional internal flows through a venturi tuve and decelerating cascades are computed and discussed.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.2 no.3
• /
• pp.134-142
• /
• 1990

A wide-angle approximation in the parabolic equation method is presented to calculate wave transformation in the shallow water. The parabolic approximation to the mild-slope equation is obtain-ed by the use of a splitting matrix, which leads to a generalized equation in form. A numerical model based on a finite difference scheme is presented and computational results are provided to test the model against the laboratory measurements of circular and elliptical shoals. The numerical results are in good agreement with most of experimental data. Therefore it can be concluded that the model shows greater capability to reproduce the characteristics of waves in the refractive focus.