• Honam Mathematical Journal
• /
• v.41 no.2
• /
• pp.403-420
• /
• 2019

The Helmholtz equation represents acoustic or electromagnetic scattering phenomena. The Method of Lines are known to have many advantages in simulation of forward and inverse scattering problems due to the usage of angle rays and Bessel functions. However, the method does not account for the jump phenomena on obstacle boundary and the approximation includes many high order Bessel functions. The high order Bessel functions have extreme blow-up or die-out features in resonance region obstacle boundary. Therefore, in particular, when we consider shape reconstruction problems, the method is suffered from severe instabilities due to the logical confliction and the severe singularities of high order Bessel functions. In this paper, two approximation formulas for the Helmholtz equation are introduced. The formulas are new and powerful. The derivation is based on Method of Lines, Huygen's principle, boundary jump relations, Addition Formula, and the orthogonality of the trigonometric functions. The formulas reduce the approximation dimension significantly so that only lower order Bessel functions are required. They overcome the severe instability near the obstacle boundary and reduce the computational time significantly. The convergence is exponential. The formulas adopt the scattering jump phenomena on the boundary, and separate the boundary information from the measured scattered fields. Thus, the sensitivities of the scattered fields caused by the boundary changes can be analyzed easily. Several numerical experiments are performed. The results show the superiority of the proposed formulas in accuracy, efficiency, and stability.

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

• 전기의세계
• /
• v.22 no.1
• /
• pp.28-34
• /
• 1973

This paper treats with the sub-optimal control problem of noninear systems by approximation method. This method involves the approximation by linearization which provides the sub-optimal solution of non-linear control problems. The result of this work shows that, in the problem in which the controlled plant is characterized by an ordinary differential equation of first order, the solution obtained by this method coincides with the exact solution of problem. In of case of the second or higher order systems, it is proved analytically that this method of linearization produces the sub-optimal solution of the given problem. It is also shown that the sub-optimality of solution by the method can be evaluated by introducing the upper and lower bounded performance indices. Discussion is made on the procedure with some illustrative examples whose performance indices are given in the quadratic forms.

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.36 no.8
• /
• pp.584-593
• /
• 1987

A new combined approximation method using Routh stability array and mean-square error (MSE) method is proposed for deriving reduced-order z-transter functions for discrete time systems. The Routh stability array is used to obtain the reduced-order denominator polynomial, and the numerator polynomial is obtained by minimizing the mean-square error between the unit step responses of the original system and reduced model. The advantages of the new combined approximation method are that the reduced model is always stable provided the original model is stable and the initial and steady-state characteristics of the original model can be preserved in the reduced model.

• Journal of the Korean Mathematical Society
• /
• v.51 no.2
• /
• pp.267-288
• /
• 2014

We propose and analyze two-scale product approximation for semilinear heat equations in the mixed finite element method. In order to efficiently resolve nonlinear algebraic equations resulting from the mixed method for semilinear parabolic problems, we treat the nonlinear terms using some interpolation operator and exploit a two-scale grid algorithm. With this scheme, the nonlinear problem is reduced to a linear problem on a fine scale mesh without losing overall accuracy of the final system. We derive optimal order $L^{\infty}((0, T];L^2({\Omega}))$-error estimates for the relevant variables. Numerical results are presented to support the theory developed in this paper.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.41 no.12
• /
• pp.1736-1738
• /
• 2016

In modern digital communication systems adapting high-order modulation and high performance channel code, log-likelihood ratios involving the repeated calculations of the logarithm of sum of exponential functions are necessary for demodulation and decoding. In this paper, the approximation methods called Min and MinC are applied to demodulation and decoding together and their complexity and joint performance are analyzed.

• Journal of the Korean Society of Safety
• /
• v.26 no.4
• /
• pp.80-86
• /
• 2011

Structural optimization algorithm of steel box girder bridges using improved higher-order approximate reanalysis technique is proposed in this paper. The proposed approximation method is a generalization of the convex approximation method. The order of the approximate reanalysis for each function is analytically adjusted in the optimization process. This self-adjusted capability makes the approximate structural analysis values conservative enough to maintain the optimum design point of the approximate problem. The efficiency of proposed optimazation algorithm, compared with conventional algorithm, is successfully demonstrated in the steel box girder bridges. The efficiency and robustness of proposed algorithm is also demonstrated in practical steel box girder bridges.

• The Journal of the Acoustical Society of Korea
• /
• v.33 no.4
• /
• pp.232-237
• /
• 2014

A passive sonar of warship is composed of several directional or omni-directional sensors. In order to model the acoustic signal received into a warship sonar, the wave propagation modeling is usually required from arbitrary noise source to all sensors equipped to the sonar. However, the full calculation for all sensors is time-consuming and the performance of sonar simulator deteriorates. In this study, we suggest an asymptotic method to estimate the sonar signal arrived to sensors adjacent to the reference sensor, where it is assumed that all information of eigenrays is known. This method is developed using Taylor series for the time delay of eigenray and similar to Fraunhofer and Fresnel approximation for sonar aperture. To validate the proposed method, some numerical experiments are performed for the passive sonar. The approximation when the second-order term is kept is vastly superior. In addition, the error criterion for each approximation is provided with a practical example.

• International Journal of CAD/CAM
• /
• v.13 no.1
• /
• pp.1-11
• /
• 2013

Just by adjusting the control points iteratively, progressive iterative approximation (PIA) presents an intuitive and straightforward scheme such that the resulting limit curve (surface) can interpolate the original data points. In order to obtain more flexibility, adjusting only a subset of the control points, a new method called local progressive iterative approximation (LPIA) has also been proposed. But to this day, there are two problems about PIA and LPIA: (1) Only an approximation process is discussed, but the accurate convergence curves (surfaces) are not given. (2) In order to obtain an interpolating curve (surface) with high accuracy, recursion computations are needed time after time, which result in a large workload. To overcome these limitations, this paper gives an explicit matrix expression of the control points of the limit curve (surface) by the PIA or LPIA method, and proves that the column vector consisting of the control points of the PIA's limit curve (or surface) can be obtained by multiplying the column vector consisting of the original data points on the left by the inverse matrix of the collocation matrix (or the Kronecker product of the collocation matrices in two direction) of the blending basis at the parametric values chosen by the original data points. Analogously, the control points of the LPIA's limit curve (or surface) can also be calculated by one-step. Furthermore, the $G^1$ joining conditions between two adjacent limit curves obtained from two neighboring data points sets are derived. Finally, a simple LPIA method is given to make the given tangential conditions at the endpoints can be satisfied by the limit curve.

• Journal of applied mathematics & informatics
• /
• v.27 no.1_2
• /
• pp.441-452
• /
• 2009

In this paper, a fifth order numerical method is presented for solving singularly perturbed differential-difference equations with negative shift. In recent papers the term negative shift has been using for delay. Similar boundary value problems are associated with expected first exit time problem of the membrane, potential in models for neuron and in variational problems in control theory. In the numerical treatment for such type of boundary value problems, first we use Taylor approximation to tackle terms containing small shifts which converts it to a boundary value problem for singularly perturbed differential equation. The two point boundary value problem is transformed into general first order ordinary differential equation system. A discrete approximation of a fifth order compact difference scheme is presented for the first order system and is solved using the boundary conditions. Several numerical examples are solved and compared with exact solution. It is observed that present method approximates the exact solution very well.