• East Asian mathematical journal
• /
• v.35 no.5
• /
• pp.597-612
• /
• 2019

In this study, we propose a new logarithmic barrier approach to solve linear programming problem using the projective method of Karmarkar. We are interested in computation of the direction by Newton's method and of the step-size using minorant functions instead of line search methods in order to reduce the computation cost. Our new approach is even more beneficial than classical line search methods. We reinforce our purpose by many interesting numerical simulations proved the effectiveness of the algorithm developed in this work.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.1
• /
• pp.104-112
• /
• 1998

The mold filling process has been a central issue in the development of numerical methods to solve the casting processes. A mold filling which is inherently transient free surface fluid flow, is important because the quality of casting highly depends on such phenomenon, Most of the existing numerical schemes to solve mold filling process have severe limitations in time step restrictions or Courant criteria since explicit time integration is used. Therefore, a large computation time is required to analyze casting processes. In this study, the well known SOLA-VOF method has been modified implicitly to simulate the mold filling process. Solutions to example filling problems show that the proposed method is more efficient in computation time than the original SOLA -VOF method.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.15 no.12
• /
• pp.1095-1102
• /
• 2003

Flow characteristics of turbulent steady fluid flow past a cylinder in rectangular duct are measured by 5 W laser doppler velocity meter. The fluid flow is also computed by commercial software of STAR-CD for comparison between the measurement and computation. The turbulent models applied in the computations are standard K-epsilon model, RNG K-epsilon model and Chen K-epsilon model. Acurracy of standard K-epsilon model is a little bit better than acurracies of other models even though those models have almost the same order of error compared to measured data. The computations predict satisfactorily the measured velocity profiles at middle section of the circular cylinder before the fluid flow diverges. However, there are some disagreements between them at down stream from the circular cylinder.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.7 no.5
• /
• pp.440-446
• /
• 1996

The multiple-ridged circular waveguides is analyzed using Fourier series and the mode matching technique. The enforcement of the boundary conditions yields the simultaneous equations for the field coefficient inside the waveguides. The simultaneous equations are solved to represent a dispersion relation in an analytic series form. The numerical computation is performed to illustrate the behavior of the cutoff wavenumbers in terms of number, length and angle of ridges. The presented series solution is exact and rapidly-convergent so that it is efficient for numerical computation. A simple dispersion relation based on the dominant mode analysis is obtained and is shown to be very accurate for most practical applications.

• Journal of applied mathematics & informatics
• /
• v.10 no.1_2
• /
• pp.63-73
• /
• 2002

The governing equation locating component centers in the degree-n bifurcation set is a polynomial with a very high degree and its root-finding lacks numerical accuracy. The equation is transformed to have its degree reduced by a factor(n-1). Newton's method applied to the transformed equation improves the accuracy with properly chosen initial values. The numerical implementation is done with Maple V using a large number of computational precision digits. Many cases are studied for 2 $\leq$ n $\leq$ 25 and show a remarkably improved computation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.2
• /
• pp.202-212
• /
• 1997

In this study, performance of the multi-stage explicit methods for numerical computation of the unsteady Navier-Stokes equations is investigated. Three methods under consideration are 1 st-, 2 nd-, and 4 th-order Runge-Kutta (R-K) methods. Compared in this estimation is stability, accuracy, and CPU time of each method. The computational codes developed are applied to the two-dimensional flow in a square cavity driven by an oscillating lid. It turned out that at Reynolds number 400, the 1 st-order R-K method is the best, while at 3200 the 2 nd-order R-K is recommended. At higher Reynolds numbers, it is conjectured that the 4 th-order R-K method will be the best algorithm among three due to its highest stability.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1996.04a
• /
• pp.316-321
• /
• 1996

The authors has developed the transfer influence coefficient method(TICM) which is an algorithm for the analysis of vibration suitable for a personal computer and applied to many structures such as straight-line beam, plate and shell structures. But TICM can't be applied to the closed loop system and is required extremely much time on the computation of complicate and large structures. Therefore, we now suggest the transfer stiffness coefficient method(TSCM) overcoming these two problems, which consists on the concept of the substructure synthesis method and transfer influence coefficient method. TSCM is formulated for the free vibration analyses of a straight-line distributed beam structure and confirmed by the results of numerical computation to the numerical high accuracy and the high speed.

• Journal of Mechanical Science and Technology
• /
• v.18 no.4
• /
• pp.718-730
• /
• 2004

The solution for the natural convection in internally finned horizontal annuli is obtained by using a numerical simulation of time-dependent and two-dimensional governing equations. The fins existing in annuli influence the flow pattern, temperature distribution and heat transfer rate. The variations of the On configuration suppress or accelerate the free convective effects compared to those of the smooth tubes. The effects of fin configuration, number of fins and ratio of annulus gap width to the inner cylinder radius on the fluid flow and heat transfer in annuli are demonstrated by the distribution of the velocity vector, isotherms and streamlines. The governing equations are solved efficiently by using a parallel implementation. The technique is adopted for reduction of the computation cost. The parallelization is performed with the domain decomposition technique and message passing between sub-domains on the basis of the MPI library. The results from parallel computation reveal in consistency with those of the sequential program. Moreover, the speed-up ratio shows linearity with the number of processor.

• The Journal of the Acoustical Society of Korea
• /
• v.19 no.7
• /
• pp.85-89
• /
• 2000

In this paper, by extending the previously developed self- and mutual-radiation impedance calculation method for a regular-square vibrating surface by using numerical series, we proposed a method to obtain the self-radiation impedance of a rectangular transducer with an arbitrary integer ratio of the length to width. The proposed method exhibits high accuracy and a short computation time. After investigating the accuracy and computation time as the number of elements changes, we have calculated the self-radiation impedance of several rectangular transducers, and compared the results with those in the literature.

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

Unsteady compressible Euler equation is solved and the high-order, high-resolution numerical solver, physical boundary condition, adaptive nonlinear artificial dissipation model and conformal mapping are applied to computation of steady transonic flow and unsteady acoustics. The acoustic characteristics of axi-symmetric duct and two dimensional straight/S channel are studied and the computation results shows good agreements with linear analysis. In transonic case, local time stepping and canceling-the-residual techniques are used for convergence acceleration. The aspect of flow and acoustics in S-channel and the Pattern of noise radiation is changed by inflow Mach no. and static pressure at fan-face.