• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.4
• /
• pp.18-27
• /
• 2002

Aerodynamic shape optimization of two-dimensional airfoils in inviscid compressible flows is performed using a continuous adjoint formulation on unstructured meshes. Accurate evaluation of the gradient is achieved by using a reconstruction scheme based on the Laplacian averaging. A least-square method with extended stencil is used for flow gradient calculations. Proper convergence criterion is studied on Euler and adjoint equations for efficient design. The present method has been applied to RAE2822 and NACA0012 airfoils such that wave drag can be minimized by removing the shock wave. An inverse design is also performed to recover the shock wave on the designed RAE2822 airfoil.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.221-224
• /
• 2003

A numerical study is carried out for the detonation wave propagation through a T-branch. The T-branch is a crucial part of the combustion wave igniter, a novel concept of rocket ignition system aimed for the simultaneous ignition of multiple combustion chambers by delivering detonation waves. Euler equation and induction parameter equation are used as governing equations with a reaction term modeled from the chemical kinetics database obtained from a detailed chemistry mechanism. Second-order accurate implicit time integration and third-order space accurate TVD algorithm were used for solution of the coupled equations. Over two-million grid points enabled the capture of the dynamics of the detonation wave propagation including the degeneration and re-initiation phenomena, and some of the design factors were be obtained for the CWI flame tubes.

• Journal of Korea Water Resources Association
• /
• v.32 no.6
• /
• pp.607-616
• /
• 1999

A fractional step finite difference model for the longitudinal dispersion of nonconservative contaminants is developed. It is based on splitting the longitudinal dispersion equation into a set of three equations each to be solved over a one-third time step. The fourth-order Holly-Preissmann scheme, an analytic solution, and the Crank-Nicholson scheme are used to solve the equations for the pure advection, the first-order decay, and the diffusion, respectively. To test the model, it is applied to simulate the longitudinal dispersion of continuous source released into a nonuniform flow field as well as the dispersion of an instantaneous source in a uniform flow field. The results are compared with the exact solution and those computed by an existing model. Compared to the existing model which uses Euler method for the first-order decay equation, the present model yield more accurate results as the decay coefficient increases.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.78-84
• /
• 2001

A dual-time stepping algorithm combined with a parallelized multigrid DADI method is presented to predict the dynamic damping coefficients. The Basic Finner model is chosen to validate the prediction capability of the present unsteady Euler method. The linearity of the pitch- and roll-damping coefficients is shown in the low angular rates and the interesting large drop and stiff increment in transonic region for roll-damping coefficients are explained in detail. Through the analysis for the pressure distributions at Mach number 1.0 to 1.2, the sudden drop results from the normal shock and the stiff increment of roll-damping reflects the transition of the normal shock to the oblique shock. The results also show that the Euler equations can give the damping coefficients with a comparable accuracy.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.175-181
• /
• 1995

Three-Dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for the various contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreements.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.21 no.4
• /
• pp.315-324
• /
• 2011

This paper presents a historical study on the derivation of the differential equation of motion for the single-degree-of-freedom m-k system with the harmonic excitation. It was Euler for the first time in the history of vibration theory who tackled the equation of motion for that system analytically, then gave the solution of the free vibration and described the resonance phenomena of the forced vibration in his famous paper E126 of 1739. As a result of the chronological progress in mechanics like pendulum condition from Galileo to Euler, the author asserts two conjectures that Euler could apply to obtain the equation of motion at that time.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.3
• /
• pp.27-36
• /
• 2002

Numerical stability is studied based on numerics and mathematics. It is frequently observed in the region where velocity is zero. In that region, the Euler equation have numerous solutions and, thus, it is impossible to determine a unique solution with only governing equations. However, a unique solution can be determined by additional outer flow conditions or outer numerical discontinuity calculation since the information or a unique solution under undisturbed conditions is lost by disturbances. In this reason, the numerical scheme comsistent with Euler equations cannot remove shock instability completely.

• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.58-65
• /
• 2014

The adaptive wavelet method is studied for the enhancement of computational efficiency of three-dimensional flows. For implementation of the method for three-dimensional Euler equation, wavelet decomposition process is introduced based on the previous two-dimensional adaptive wavelet method. The order of numerical accuracy of an original solver is preserved by applying modified thresholding value. In order to assess the efficiency of the proposed algorithm, the method is applied to the computation of flow field around ONERA-M6 wing in transonic regime with 4th and 6th order interpolating polynomial respectively. Through the application, it is confirmed that the three-dimensional adaptive wavelet method can reduce the computational time while conserving the numerical accuracy of an original solver.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2013.04a
• /
• pp.528-532
• /
• 2013

In this study, an analytical model of micro-beam structure including thermoelastic damping with irregularly distributed masses is investigated. The significance of thermoelastic damping for micro-scale mechanical resonators is evaluated to design -with high quality factor(Q-factor). The beam model of this work is based on Euler-Bernoulli beam theory. In order to determine the natural frequency of the model, energy method is applied. Also, the thermoelatic damping effects are considered by using heat conduction equations, and the Q-factor can be determined. To derive the equation of motion, non-dimensionalization is employed for systematic analysis. Results of the model are verified, and present mode shapes and Q-factors for the micro-beam with thermoelastic damping including random point masses.

• Computational Structural Engineering
• /
• v.10 no.4
• /
• pp.177-184
• /
• 1997

In this paper, numerical methods are developed for solving the elastica of simple beams with variable-arc-length subjected to a point loading. The beam model is based on Bernoulli-Euler beam theory. The Runge-Kutta and Regula-Falsi methods, respectively, are used to solve the governing differential equations and to compute the beam's rotation at the left end of the beams. Extensive numerical results of the elastica responses, including deflected shapes, rotations of cross-section and bending moments, are presented in non-dimensional forms. The possible maximum values of the end rotation, deflection and bending moment are determined by analyzing the numerical data obtained in this study.