• Journal of computational fluids engineering
• /
• v.13 no.3
• /
• pp.8-13
• /
• 2008

We present a numerical simulation technique and some preliminary results of the impact and spreading of a droplet containing particles on the solid substrate in 2D. We used the 2nd-order Adams-Bashforth / Crank-Nicholson method to solve the Navier-Stokes equation and employed the level-set method with the continuous surface stress for description of droplet spreading with interfacial tension. The impact velocity has been generated by the instantaneous gravity. The distributed Lagrangian-multipliers method has been combined for the implicit treatment of rigid particles and the discontinuous Galerkin method has been used for the stabilization of the interface advection equation. We investigated the droplet spreading by the inertial force and discussed effects of the presence of particles on the spreading behavior using an example problem. We observed reduced oscillation and spread for the particulate droplet.

• Proceedings of the Korea Contents Association Conference
• /
• 2013.05a
• /
• pp.443-444
• /
• 2013

빈발하고 있는 대규모 홍수와 자연재해는 정확도가 높은 하천 흐름 수치해석 모델에 대한 관심의 증대로 이어지고 있다. 현재 하천에서 발생하는 일반적인 흐름은 기존에 개발된 여러 형태의 천수방정식을 지배방정식으로 하는 수치기법에 의해 해석되고 있으나, 연속적이지 않은 형태의 흐름을 해석하거나 매우 정확한 해석을 필요로 하는 경우에는 기존의 수치해석기법은 많은 한계를 보여 주고 있다. 본 연구에서는 불연속 갤러킨 기법 기반의 흐름 모델을 개발하고, 이를 이용하여 천이류로 분류되는, 댐 붕괴파, 둔덕위 흐름과 2차원 사류의 모의에 적용하여 기존의 수치해와 잘 일치함을 확인하였다.

• Computational Structural Engineering
• /
• v.11 no.4
• /
• pp.187-196
• /
• 1998

불연속 갤러킨 정식화에 기초를 둔 시간적분법에 대하여 시간을 변수로 한 유한요소적 접근법을 시도하였다. 단일 형상함수와 두 형상함수 정식화에 대해 각각 선형, 이차 형상함수를 적용하여 모두 네 종류의 시간적분법을 유도하였으며, 각 방법에 대하여 시간시텝의 증가에 따른 변위와 속도의 관계를 나타내는 증폭행렬을 계산하였다. 유도된 방법들의 성능을 평가하기 위하여 부하가 갑자기 변화는 진동 문제를 해석하고 변위의 오차를 비교하였다. 네 가지의 방법에 대하여 국부 오차 추정치를 개발하였으며, 오차 추정치의 정확도를 수치예를 이용하여 평가하였다. 단일 형상함수 정식화에서 이차 형상함수를 이용한 오차 추정치가 실제 국부오차를 잘 나타내었으며 유도된 오차 추정치는 시간간격제어 기법에서 시간간격의 크기를 결정하는 척도로 이용 가능하다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.470-474
• /
• 2016

Madsen et al. (2002)이 제안한 일차원 고차 Boussinesq 방정식에 대하여 불연속갤러킨 유한요소법(Discontinuous Galerkin Finite Element Method)을 적용하였다. 연속적인 Boussinesq 방정식에서 각 요소경계에 불연속을 허용할 수 있도록 공간차분하고, 시간방향으로 4차 Runge-Kutta 시간적분법, 각 요소사이에는 Lax-Friedrichs 수치흐름률을 사용하였다. 계산영역의 양쪽에 불필요한 파랑의 반사를 억제하도록 흡수층을 설치하였으며, 영역 내부에서 조파할 수 있도록 하였다. Luth et al.(1994)의 수중잠제 실험에 적용하여 관측값과 잘 일치함을 확인하였다.

• Advances in aircraft and spacecraft science
• /
• v.9 no.5
• /
• pp.463-477
• /
• 2022

This study compares various ways of calculating flows for the problems with the presence of shock waves by first-order schemes and higher-order DG method on the tests from the Quirk list, namely: Quirk's problem and its modifications, shock wave diffraction at a 90 degree corner, the problem of double Mach reflection. It is shown that the use of HLLC and Godunov's numerical schemes flows in calculations can lead to instability, the Rusanov-Lax-Friedrichs scheme flow can lead to high dissipation of the solution. The most universal in heavy production calculations are hybrid schemes flows, which allow the suppression of the development of instability and conserve the accuracy of the method.

• Advances in aircraft and spacecraft science
• /
• v.9 no.5
• /
• pp.415-431
• /
• 2022

Secondary flows have a huge impact on losses generation in modern low pressure gas turbines (LPTs). At design point, the interaction of the blade profile with the end-wall boundary layer is responsible for up to 40% of total losses. Therefore, predicting accurately the end-wall flow field in a LPT is extremely important in the industrial design phase. Since the inlet boundary layer profile is one of the factors which most affects the evolution of secondary flows, the first main objective of the present work is to investigate the impact of two different inlet conditions on the end-wall flow field of the T106A, a well known LPT cascade. The first condition, labeled in the paper as C1, is represented by uniform conditions at the inlet plane and the second, C2, by a flow characterized by a defined inlet boundary layer profile. The code used for the simulations is based on the Discontinuous Galerkin (DG) formulation and solves the Reynolds-averaged Navier-Stokes (RANS) equations coupled with the Spalart Allmaras turbulence model. Secondly, this work aims at estimating the influence of viscosity and turbulence on the T106A end-wall flow field. In order to do so, RANS results are compared with those obtained from an inviscid simulation with a prescribed inlet total pressure profile, which mimics a boundary layer. A comparison between C1 and C2 results highlights an influence of secondary flows on the flow field up to a significant distance from the end-wall. In particular, the C2 end-wall flow field appears to be characterized by greater over turning and under turning angles and higher total pressure losses. Furthermore, the C2 simulated flow field shows good agreement with experimental and numerical data available in literature. The C2 and inviscid Euler computed flow fields, although globally comparable, present evident differences. The cascade passage simulated with inviscid flow is mainly dominated by a single large and homogeneous vortex structure, less stretched in the spanwise direction and closer to the end-wall than vortical structures computed by compressible flow simulation. It is reasonable, then, asserting that for the chosen test case a great part of the secondary flows details is strongly dependent on viscous phenomena and turbulence.

• The Journal of the Acoustical Society of Korea
• /
• v.15 no.2
• /
• pp.72-78
• /
• 1996

This thesis deals with a method to solve a two-and-one-half-dimensional ($2\frac12$ D) problem, which means that the ocean environment is two-dimensional whereas the source is fully three-dimensionally propagating, including three-dimensional refraction phenomena and three-dimensional back-scattering, using two-dimensional two-way parabolic equation method combined with Fourier synthesis. Two dimensional two-way parabolic equation method uses Galerkin's method for depth and Crank-Nicolson method and alternating direction for range and provides a solution available to range-dependent problem with wave-field back-scattered from discontinuous interface. Since wavenumber, k, is the function of depth and vertical or horizontal range, we can reduce a dimension of three-dimensional Helmholtz equation by Fourier transforming in the range direction. Thus transformed two-dimensional Helmholtz equation is solved through two-way parabolic equation method. Finally, we can have the $2\frac12$ D solution by inverse Fourier transformation of the spectral solution gained from in the last step. Numerical simulation has been carried out for a canonical ocean environment with stair-step bottom in order to test its accuracy using the present analysis. With this spectral parabolic equation method, we have examined three-dimensional acoustic propagation properties in a specified site in the Korean Straits.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.4
• /
• pp.703-713
• /
• 2002

Recently, an improved EFG(Element-Free Galerkin) crack analysis technique, which includes a discontinuous approximation and a singular basis function on the auxiliary supports, was developed. The technique is able to accurately analyze the crack propagation problem without any modification of the analysis model; however, it shows some dependency on the analysis parameters used. In this study, the effect of analysis parameters such as the size of compact support, dilation parameter, the smoothness of shape function around the crack tip, and the number of node using auxiliary supports on the accuracy of solution has been investigated. Through a patch test with a crack, relative L₂ error norm of stresses and the stress intensity factor were computed and compared for various analysis parameters and the results were presented as guidelines for adequate choice of analysis parameters.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.15 no.2
• /
• pp.101-108
• /
• 2007

This paper presents an efficient meshless method for analyzing cracked piezoelectric structures subjected to mechanical and electrical loading. The method employs an element free Galerkin (EFG) formulation and an enriched basic function as well as special shape functions that contain discontinuous derivatives. Based on the moving least squares (MLS) interpolation approach, The EFG method is one of the promising methods for dealing with problems involving progressive crack growth. Since the method is meshless and no element connectivity data are needed, the burdensome remeshing procedure required in the conventional finite element method (FEM) is avoided. The numerical results show that the proposed method yields an accurate near-tip stress field in an infinite piezoelectric plate containing an interior hole. Another example is to study a ceramic multilayer actuator. The proposed model was found to be accurate in the simulation of stress and electric field concentrations due to the abrupt end of an internal electrode.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.1
• /
• pp.43-57
• /
• 1999

The fluid flow and heat transfer in a thin liquid film are investigated numerically. The flow Is assumed to be two-dimensional laminar and surface tension is considered. The most important characteristics of this flow is the existence of a hydraulic jump through which the flow undergoes very sharp and discontinuous change. Arbitrary Lagrangian-Eulerian(ALE) method is used to describe moving free boundary and a modified SIMPLE algorithm based on streamline upwind Petrov-Galerkin(SUPG) finite element method is used for time marching iterative solution. The numerical results obtained by solving unsteady full Navier-Stokes equations are presented for planar and radial flows subject to constant wall temperature or constant wall heat flux, and compared with available experimental data. It Is discussed systematically how the inlet Reynolds and Froude numbers and surface tension affect the formation of a hydraulic jump. In particular, the effect of temperature dependent fluid properties is also discussed.