• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.177-185
• /
• 1999

The airfoil design optimization procedures based on the Navier-Stokes equations were developed, This procedure enables more realistic and practical transonic airfoil designs. The modified Hicks-Henne functions were used to generate the shape of airfoils. Five Hick-Henne functions were used to design upper surface of airfoil only. To enhance the ability of Hick-Henne function to generate various airfoil shape with limited number of functions, the positions of control points were adjusted through optimization procedure. The design procedure was applied to the single-point design for the drag minimization problem with lift and area constraints. The result shows the capability of the procedure to generate much realistic airfoils with very small drag-creep in the low transonic regime. This is mainly due to the viscosity effect of Navier-Stokes flow analysis. However, in the higher transonic range tile drag-creep appears. The multi-point design is shown to be an effective way to avoid the drag-creep and improve off-design performance which is very similar in the Euler design.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.565-572
• /
• 2004

Numerical stability analysis of one-dimensional axial flow in solid rocket motors is performed based on the Euler equation coupled with an unsteady combustion equation of solid propellant. In order to check the numerical scheme, behavior of a standing wave in a closed tube is examined. A standing wave in solid rocket motor decays or grows depending on the total effect of propellant combustion, nozzle flow, and so on. The stability boundary of the fundamental mode standing wave is determined by changing one of the combustion parameters. In addition growth rates of the wave are calculated numerically in relatively low Mach number flow region for the motors with different port and nozzle throat diameters. The results obtained here agree well with the approximate solution. The same scheme is applied to a motor with shorter length and L＊-instability is observed.

• Proceedings of the KIEE Conference
• /
• 1993.07b
• /
• pp.683-685
• /
• 1993

Semiconductor switches are modelled as binary inductors, as a very low value of inductance during conduction, and infinite value of inductance otherwise. The system matrix is partitioned in such a manner as to permit efficient handling of switch status. The backward Euler method of integration is used for the solution of equation to ensure convergence. The application of thristor switching to induction motor speed control has resulted in a number of unconventional supply systems. In this paper, an analytic method for predicting the steady-state of performance of system is presented.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2001.05a
• /
• pp.52-57
• /
• 2001

Laminar two-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for the low Reynolds number (Re=164∼280). The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. The convection term is applied by the 7th order up wind scheme and the pressure and viscosity terms are applied by the 4th order central difference. The grid system makes use of the regular grid system and it is generated by an equation. The calculated results of drag coefficients, lift coefficients, pressure distributions, and vorticity contours and other information are compared with experimental and numerical ones. These results obtained by the present DNS show good agreement with the previous studies.

• 한국가시화정보학회:학술대회논문집
• /
• 2002.11a
• /
• pp.113-116
• /
• 2002

The shock wave discharged from an annular duct leads to very complicated flow features, such as Mach stem, spherical waves, and vortex rings. In the current study, the merging phenomenon and propagation characteristics of the shock wave are numerically investigated using a CFD method. The Harten-Yee's total variation diminishing (TVD) scheme is used to the unsteady, axisymmetric, two-dimensional, compressible Euler equations. The Mach number of incident shock wave $M_s$ is varied in the range below 2.0. The computational results are visualized to observe the major features of the annular shock waves discharged from the tube. On the symmetric axis, the peak pressure produced by the shock wave and its location depend upon strongly the radius of the annular tubes. A Mach stem is generated along the symmetric axis of the annular tubes.

• East Asian mathematical journal
• /
• v.35 no.1
• /
• pp.23-30
• /
• 2019

Since Mittag-Leffler introduced the so-called Mittag-Leffler function, a number of its extensions have been investigated due mainly to their applications in a variety of research subjects. Shukla and Prajapati presented a lot of formulas involving a generalized Mittag-Leffler function in a systematic manner. Motivated mainly by Shukla and Prajapati's work, we aim to investigate a generalized multi-index Mittag-Leffler function and, among possible numerous formulas, choose to present several formulas involving this generalized multi-index Mittag-Leffler function such as a recurrence formula, derivative formula, three integral transformation formulas. The results presented here, being general, are pointed out to reduce to yield relatively simple formulas including known ones.

• Communications of the Korean Mathematical Society
• /
• v.34 no.4
• /
• pp.1105-1115
• /
• 2019

For positive integers n and k, let $S_k(n)$ and $S^{\prime}_k(n)$ be the sums of the elements in the finite sets {$x^k:1{\leq}x{\leq}n$, (x, n) = 1} and {$x^k:1{\leq}x{\leq}n/2$, (x, n) = 1}respectively. The formulae for both $S_k(n)$ and $S^{\prime}_k(n)$ are established. The explicit formulae when k = 1, 2, 3 are also given.

• Geomechanics and Engineering
• /
• v.32 no.5
• /
• pp.541-551
• /
• 2023

Snap-buckling is one of the main failure modes of structures, because it will lead to the reduction of structural bearing capacity, durability loss and even structural damage. Boundary condition plays an important role in the research of engineering mechanics. Further discussion on the boundary conditions problems will help to analyze the dynamic and static behavior of structures more accurately. Therefore, in order to understand the dynamic and static behavior of curved beams more comprehensively, this paper mainly studies the nonlinear snap-through buckling and forced vibration characteristics of functionally graded graphene reinforced composites (FG-GPLRCs) curved beams with two different boundary conditions (including clamped-hinged and hinged-hinged) using Euler-Bernoulli beam theory (E-BBT). In addition, the effects of the curved beam radius, the GLPs distributions, number of GLPs layers, the mass fraction of GLPs and elastic foundation parameters on the nonlinear snap-through buckling and forced vibration behavior are discussed respectively.

• Journal of Ocean Engineering and Technology
• /
• v.37 no.6
• /
• pp.256-265
• /
• 2023

Many numerical methods have been developed since 1961, but unresolved issues remain. This study developed a numerical method to address these issues and determine the coefficients and properties of rotational waves with a shear current in a finite water depth. The number of unknown constants was reduced significantly by introducing a wavelength-independent coordinate system. The reference depth was calculated independently using the shooting method. Therefore, there was no need for partial derivatives with respect to the wavelength and the reference depth, which simplified the numerical formulation. This method had less than half of the unknown constants of the other method because Newton's method only determines the coefficients. The breaking limit was calculated for verification, and the result agreed with the Miche formula. The water particle velocities were calculated, and the results were consistent with the experimental data. Dispersion relations were calculated, and the results are consistent with other numerical findings. The convergence of this method was examined. Although the required series order was reduced significantly, the total error was smaller, with a faster convergence speed.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.10
• /
• pp.939-946
• /
• 2008

Euler codes or Navier-Stokes codes for compressible flows suffer severe degradation in convergence as Mach number approaches zero. The convergence problem arose from the wide disparity in characteristic speeds can be solved using preconditioning methods without large modifications. In this paper, a preconditioned RANS(Reynolds Averaged Navier-Stokes) solver is developed for analysis of low Mach number flows. In order to validate the method, computational examples are chosen and the results are compared with the experimental data and the existing computed results showing a good accuracy and convergence characteristics for steady inviscid, laminar and turbulent flows at low Mach number.