• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.359-366
• /
• 2008

Current state and perspective of DNS of turbulence and turbulent combustion are discussed with feature trend of the fastest supercomputer in the world. Based on the perspective of DNS of turbulent combustion, possibility of perfect simulations of IC engine is shown. In 2020, the perfect simulation will be realized with 30 billion grid points by 1EXAFlops supercomputer, which requires 4 months CPU time. The CPU time will be reduced to about 4 days if several developments were achieved in the current fundamental researches. To shorten CPU time required for DNS of turbulent combustion, two numerical methods are introduced to full-explicit full-compressible DNS code. One is compact finite difference filter to reduce spatial resolution requirements and numerical oscillations in small scales, and another is well-known point-implicit scheme to avoid quite small time integration of the order of nanosecond for fully explicit DNS. Availability and accuracy of these numerical methods have been confirmed carefully for auto-ignition, planar laminar flame and turbulent premixed flames. To realize DNS of IC engine with realistic kinetic mechanism, several DNS of elemental combustion process in IC engines has been conducted.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.359-366
• /
• 2008

Current state and perspective of DNS of turbulence and turbulent combustion are discussed with feature trend of the fastest supercomputer in the world. Based on the perspective of DNS of turbulent combustion, possibility of perfect simulations of IC engine is shown. In 2020, the perfect simulation will be realized with 30 billion grid points by 1EXAFlops supercomputer, which requires 4 months CPU time. The CPU time will be reduced to about 4 days if several developments were achieved in the current fundamental researches. To shorten CPU time required for DNS of turbulent combustion, two numerical methods are introduced to full-explicit full-compressible DNS code. One is compact finite difference filter to reduce spatial resolution requirements and numerical oscillations in small scales, and another is well-known point-implicit scheme to avoid quite small time integration of the order of nanosecond for fully explicit DNS. Availability and accuracy of these numerical methods have been confirmed carefully for auto-ignition, planar laminar flame and turbulent premixed flames. To realize DNS of IC engine with realistic kinetic mechanism, several DNS of elemental combustion process in IC engines has been conducted.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.1
• /
• pp.1-14
• /
• 1998

A two-dimensional numerical model for tidal currents based on the depth averaged equation is developed. The mode1 employs a rectangular grid system for its simplicity in the application of complicate coastal shore lines. To raise computing efficiency, implicit approximate factorization scheme is implemented in solving governing equations. An upwind-differencing is used to discretize convective terms, which provides a numerical dissipation automatically and suppresses any oscillations caused by nonlinear instabilities. Some numerical tests are made against the analytic solutions of a linearized shallow water equation to validate the developed numerical scheme, and comparisons of the model prediction with the analytic solution are satisfactory. As a real application, the tidal currents are computed on the Inchon area where the tidal currents are important for the design of new canal which is under construction.

• Journal of the Korean Institute of Navigation
• /
• v.15 no.2
• /
• pp.87-96
• /
• 1991

The change of water level at Mokpo Harbour and its adjacent coastal area due to the construction of the Youngsan Estuary Barrage and the Third Land Reclamation Work of estuary barren had been roughly expected. Periodical floods, which occur 2 times per month, are also being observed at the low lying commercial areas near the Mokpo Old Harbor. Although it is said that the highest tidal current component among the tidal current records at the approaching channel to Mokpo Harbor is reduced to 6 kts, because of the esturary barrage, they do not give any precise statement or a deep analysis for the flooding and periodical water level change under certain environmental conditions. Moreover, they never tried the analysis of development plan considering the natural disaster such as typhoon or other extreme conditions. Thus, it is necessary to collect and analyze the data related to floodings, harbor oscillations, currents, and water quality , etc. because of the development considering the extreme condition. Thus, it is necessary to collect and analyze the data related to floodings, harbor oscillations, currents, and water quality, etc. because of the development considering the extreme condition and to evaluate the field observation and measurement, including the numerical model simulation based on the scientific approaches. This study deals the problem of the water level change among the integrated analyses of the coastal area changes. The result can be used for the integrated planning to give a strong foundation and it will contribute to the development of local area.

• Journal of computational fluids engineering
• /
• v.13 no.1
• /
• pp.21-27
• /
• 2008

A numerical study for laminar flow in the entrance region of helical tubes for uniform inlet velocity conditions is carried out by means of the finite volume method to investigate the effects of Reynolds number, pitch and curvature ratio on the flow development. This results cover a curvature ratio range of 1/10$\sim$1/320, a pitch range of 0.0$\sim$3.2, and a Reynolds number range of 125$\sim$2000. It has been found that the curvature ratio does significantly effect on the angle of flow development, but the pitch and Reynolds number do not. The characteristic angle $\phi_c(=\phi/\sqrt{\delta})$, or the non-dimensional length $\overline{l}(=l\sqrt{\delta}cos(atan\lambda)/d)$ can be used to represent the flow development for uniform inlet velocity conditions. In uniform inlet velocity conditions, the growth of boundary layer delays the flow development attributed to centrifugal force, and in which conditions the amplitude of flow oscillations is smaller than that in parabolic inlet velocity conditions. If the pitch increases or if the curvature ratio or Reynolds number decreases, the minimum friction factor and the fully developed average friction factor normalized with the friction factor of a straight tube and the flow oscillations decrease.

• Journal of computational fluids engineering
• /
• v.21 no.2
• /
• pp.1-11
• /
• 2016

Considerable efforts are required to develop a monotone, robust and stable high-order numerical scheme for solving the hyperbolic system. The discontinuous Galerkin(DG) method is a natural choice, but elimination of the spurious oscillations from the high-order solutions demands a new development of proper limiters for the DG method. There are several available limiters for controlling or removing unphysical oscillations from the high-order approximate solution; however, very few studies were directed to analyze the exact role of the limiters in the hyperbolic systems. In this study, the performance of the several well-known limiters is examined by comparing the high-order($p^1$, $p^2$, and $p^3$) approximate solutions with the exact solutions. It is shown that the accuracy of the limiter is in general problem-dependent, although the Hermite WENO limiter and maximum principle limiter perform better than the TVD and generalized moment limiters for most of the test cases. It is also shown that application of the troubled cell indicators may improve the accuracy of the limiters under some specific conditions.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.10a
• /
• pp.157-163
• /
• 2003

A systematic numerical experiment has been conducted to study droplet gasification in high pressure environments with pressure oscillations. The general frame of previous rigorous model[1] is retained but tailored for flash equilibrium calculation of vapor-liquid interfacial thermodynamics. Time-dependent conservation equations of mass, momentum, energy, and species concentrations are formulated in axisymmetric coordinate system for both the droplet interior and ambient gases. In addition, a unified property evaluation scheme based on the fundamental equation of state and empirical methods are used to find fluid thermophysical properties over the entire thermodynamic domain of interest. The governing equations with appropriate physical boundary conditions are numerically time integrated using an implicit finite-difference method with a dual time-stepping technique. A series of calculation have been carried out to investigate the gasification of an isolated n-pentane droplet in a nitrogen gas environment over a wide range of ambient pressures and frequencies. Results show that the mean pressures and frequencies of the ambient gas have strong influences on the characteristics of the droplet gasification. The amplitude of the response increases with increasing pressure, and the magnitude of the vaporization response increases with the frequency.

• International Journal of Ocean System Engineering
• /
• v.4 no.1
• /
• pp.19-28
• /
• 2014

The tension leg platform (TLP) is a vertically moored structure with excess buoyancy. The TLP is regarded as moored structure in horizontal plan, while inherit stiffness of fixed platform in vertical plane. In this paper, a numerical study using modified Morison equation was carried out in the time domain to investigate the influence of nonlinearities due to hydrodynamic forces and the coupling effect between surge, sway, heave, roll, pitch and yaw degrees of freedom on the dynamic behavior of TLP's. The stiffness of the TLP was derived from a combination of hydrostatic restoring forces and restoring forces due to cables and the nonlinear equations of motion were solved utilizing Newmark's beta integration scheme. The effect of tethers length and wave characteristics such as wave period and wave height on the response of TLP's was evaluated. Only uni-directional waves in the surge direction was considered in the analysis. It was found that for short wave periods (i.e. 10 sec.), the surge response consisted of small amplitude oscillations about a displaced position that is significantly dependent on tether length, wave height; whereas for longer wave periods, the surge response showed high amplitude oscillations about that is significantly dependent on tether length.

• Journal of Ocean Engineering and Technology
• /
• v.18 no.5
• /
• pp.15-21
• /
• 2004

An edge tone is the discrete tone or narrow-band sound produced by an oscillating free shear layer, impinging on a rigid surface. In this paper, we present a 2-D edge tone to predict the frequency characteristics of the discrete oscillations of a jet-edge feedback cycle, using the finite difference lattice Boltzmann method (FDLBM). We use a modified version of the lattice BGK compressible fluid model, adding an additional term and allowing for longer time increments, compared to a conventional FDLBM, and also use a boundary fitted coordinates system. The jet is chosen long enough in order to guarantee the parabolic velocity profile of the jet at the outlet, and the edge consists of a wedge with an angle of ${\alpha}$ = 23. At a stand-off distance, the edge is inserted along the centerline of the jet, and a sinuous instability wave, with real frequency, is assumed to be created in the vicinity of the nozzle and propagates towards the downstream. We have succeeded in capturing very small pressure fluctuations, resulting from periodical oscillations of a jet around the edge. The pressure fluctuations propagate with the speed of sound. Its interaction with the wedge produces an non-rotational feedback field, which, near the nozzle exit, is a periodic transverse flow, producing the singularities at the nozzle lips.

• International Journal of Aeronautical and Space Sciences
• /
• v.11 no.4
• /
• pp.285-302
• /
• 2010

Active aeroelastic control is an emerging technology aimed at providing solutions to structural systems that under the action of aerodynamic loads are prone to instability and catastrophic failures, and to oscillations that can yield structural failure by fatigue. The purpose of the aeroelastic control among others is to alleviate and even suppress the vibrations appearing in the flight vehicle subcritical flight regimes, to expand its flight envelope by increasing the flutter speed, and to enhance the post-flutter behavior usually characterized by the presence of limit cycle oscillations. Recently adaptive and robust control strategies have demonstrated their superiority to classical feedback strategies. This review paper discusses the latest development on the topic by the authors. First, the available control techniques with focus on adaptive control schemes are reviewed, then the attention is focused on the advanced single-input and multi-input multi-output adaptive feedback control strategies developed for lifting surfaces operating at subsonic and supersonic flight speeds. A number of concepts involving various adaptive control methodologies, as well as results obtained with such controls are presented. Emphasis is placed on theoretical and numerical results obtained with the various control strategies.