• 한국전산유체공학회지
• /
• 제13권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.

• 한국콘텐츠학회:학술대회논문집
• /
• 한국콘텐츠학회 2013년도 춘계 종합학술대회 논문집
• /
• pp.443-444
• /
• 2013

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

• 전산구조공학
• /
• 제11권4호
• /
• pp.187-196
• /
• 1998

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

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2016년도 학술발표회
• /
• 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
• /
• 제9권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
• /
• 제9권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.

• 한국음향학회지
• /
• 제15권2호
• /
• pp.72-78
• /
• 1996

본 논문에서는 2차원 양방향 포물선 방정식 법과 푸리에 변환을 이용하여 3차원 굴절현상 및 3차원 후방 산란파를 포함하는 $2\frac12$차원 문제를 푸는 방법에 대해 다루었다. 여기서 $2\frac12$ 차원 문제란 2차원적 해양환경 하에 3차원적 음원이 존재할 경우를 의미한다. 2차원 양방향 포물선 방정식법은 수치기법으로 깊이 방향과 수평거리 방향에 대해 각각 Galerkin법과 Crank-Nicolson법을 사용하며 수직 불연속 경계면에 의한 후방 산란파를 포함한 수평거리 의존 문제에 대해 유용한 해를 제공한다. 2차원 해양환경에서는 파수 k가 종 또는 횡 수평거리 방향과 깊이 방향에 대한 함수이므로 3차원 Helmholtz방정식 법을 이용해 스펙트럴 해를 구하여 다시 푸리에 역변환하면 최종 해를 구할 수 있다. 본 연구방법의 정확성을 시험하기 위해서 계단형 해저면을 갖는 간단한 해양환경에서 계산을 수행해 보았으며 대한해협의 특정지역에서의 3차원적 음파전달 특성을 살펴보았다.

• 한국전산구조공학회논문집
• /
• 제15권4호
• /
• pp.703-713
• /
• 2002

최근에 불연속 근사변위함수와 보조받침을 이용한 특이기저함수를 도입하여 균열의 불연속성과 특이성을 구현한 개선된 EFG(Element-Free Galerkin) 균열해석기법이 제안되었다. 개선된 EFG 균열해석기법은 균열의 성장에 따른 해석모형의 수정 없이도 높은 정확도로 균열전파해석을 수행할 수 있지만, 다른 무요소법과 마찬가지로 해석결과가 사용되는 해석계수에 의존하게 된다. 본 연구에서는 개선된 EFG 균열해석기법에서 사용하는 해석계수 즉, compact 받침 크기, 팽창계수, 선단주변에서의 형상함수의 평활화, 보조받침을 사용하는 절점개수가 수치해석 결과에 미치는 영향을 분석하였다. 균열문제에 대한 patch 시험을 통해 응력에 대한 L₂오차노름과 응력확대계수를 산정하여 해석계수의 영향을 분석하였으며, 그 결과는 해석계수의 선택에 대한 지침으로 제시된다.

• 한국자동차공학회논문집
• /
• 제15권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.

• 대한기계학회논문집B
• /
• 제23권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.