• 한국수자원학회논문집
• /
• 제31권5호
• /
• pp.599-612
• /
• 1998

수송방정식의 양면적은 특성으로 인하여 이송항이 지배적인 흐름에 있어서 수송방정식의 수채해석은 매우 난해하다. 특히 유한요소법을 사용하여 수치해석할 때, 상류방향으로 더 많은 가중치를 두기 위하여 변화된 시행함수를 사용한다. 이러한 방법을 Petrov-Galerkin 방법이라고 한다. 본 논문에서는 N+1 과 N+2 Petrov-Galerkin 방법을 격자 Peclet 수가 큰 수송문제에 적용하였다. 주파수맞춤 기법을 사용하여 N+2 Petrov- Galerkin 방법을 격자 Peclet 수가 큰 소송문제에 적용하였다. 주파수맞춤 기법을 사용하여 N+2 Petrov-Galerkin 방법의 적정 풍상정도를 찾아내었고, 그 결과를 토의하였다. 이 기법의 시행함수는 이송항과 확산항의 상대적 크기에 따라 그 모양이 변화된다. 수치실험을 통하여 제시된 새로운 수치해석기법의 우수성을 설명하였다.

• 대한기계학회논문집B
• /
• 제29권6호
• /
• pp.703-710
• /
• 2005

A finite element discretization of the advection and redistancing equations of level set method has been studied. It has been shown that Galerkin spatial discretization combined with Crank-Nicolson temporal discretization of the advection equation of level set yields a good result and that consistent streamline upwind Petrov-Galerkin(CSUPG) discretization of the redistancing equation gives satisfactory solutions for two test problems while the solutions of streamline upwind Petrov-Galerkin(SUPG) discretization are dissipated by the numerical diffusion added for the stability of a hyperbolic system. Furthermore, it has been found that the solutions obtained by CSUPG method are comparable to those by second order ENO method.

• 대한기계학회논문집
• /
• 제19권10호
• /
• pp.2667-2677
• /
• 1995

A new numerical algorithm using equal order linear finite element and fractional step method has been developed which is capable of analyzing unsteady fluid flow and heat transfer problems. Streamline Upwind Petrov-Galerkin (SUPG) method is used for the weighted residual formulation of the Navier-Stokes equations. It is shown that fractional step method, in which pressure term is splitted from the momentum equation, reduces computer memory and computing time. In addition, since pressure equation is derived without any approximation procedure unlike in the previously developed SIMPLE algorithm based FEM codes, the present numerical algorithm gives more accurate results than them. The present algorithm has been applied preferentially to the well known bench mark problems associated with steady flow and heat transfer, and proves to be more efficient and accurate.

• Journal of Mechanical Science and Technology
• /
• 제20권10호
• /
• pp.1741-1752
• /
• 2006

A combined Streamline Upwind Petrov-Galerkin method (SUPG) and segregated finite element algorithm for solving conjugate heat transfer problems where heat conduction in a solid is coupled with heat convection in viscous fluid flow is presented. The Streamline Upwind Petrov-Galerkin method is used for the analysis of viscous thermal flow in the fluid region, while the analysis of heat conduction in solid region is performed by the Galerkin method. The method uses the three-node triangular element with equal-order interpolation functions for all the variables of the velocity components, the pressure and the temperature. The main advantage of the presented method is to consistently couple heat transfer along the fluid-solid interface. Four test cases, which are the conjugate Couette flow problem in parallel plate channel, the counter-flow in heat exchanger, the conjugate natural convection in a square cavity with a conducting wall, and the conjugate natural convection and conduction from heated cylinder in square cavity, are selected to evaluate efficiency of the presented method.

• Journal of applied mathematics & informatics
• /
• 제22권1_2호
• /
• pp.317-329
• /
• 2006

In this paper, a singularly perturbed ordinary differential equation with non-smooth data is considered. The numerical method is generated by means of a Petrov-Galerkin finite element method with the piecewise-exponential test function and the piecewise-linear trial function. At the discontinuous point of the coefficient, a special technique is used. The method is shown to be first-order accurate and singular perturbation parameter uniform convergence. Finally, numerical results are presented, which are in agreement with theoretical results.

• 물과 미래
• /
• 제29권6호
• /
• pp.155-166
• /
• 1996

개수로내의 검변 및 급변 부정류 해석을 위해서 dynamic wave식을 기본방정식으로 하고 이를 불연속 보간함수와 upstream weighting 을 도입한 Petrov-Galerkin 기법에 의해 해석하는 유한요소모형을 개발하였다. 매트릭스 안정성 해석 결과 Petrov-Galerkin기법은 단파장에서의 선택적 감쇠능력과 위상오차에 있어 우수한 것으로 나타났다. 반면에 Preissmann기법은 단파장에서의 선택적 감쇠능력과 위상오차에 있어 열등한 것으로 나타났고, Bubnov-Galerkin 기법은 비감쇠특성을 나타내고 있어 단파장 영역에서 발산해를 일으키는 주요원인임을 확인할 수 있었다. Petrov-Galerkin 방법은 Courant수의 넓은 범위에서 높은 주파수를 가진 진행파에 대한 선택적인 감쇠와 작은 Courant 수의 범위에서 양호한 위상정도를 가지는 이상적인 조합을 나타내고 있어 점변 및 급변 부정류 해석에 있어 이상적인 기법으로 활용될 수 있을 것으로 판단되었다.

• Journal of Mechanical Science and Technology
• /
• 제20권11호
• /
• pp.1881-1890
• /
• 2006

In this paper, a Petrov-Galerkin natural element method (PG-NEM) based upon the natural neighbor concept is presented for the free vibration and dynamic response analyses of two-dimensional linear elastic structures. A problem domain is discretized with a finite number of nodes and the trial basis functions are defined with the help of the Voronoi diagram. Meanwhile, the test basis functions are supported by Delaunay triangles for the accurate and easy numerical integration with the conventional Gauss quadrature rule. The numerical accuracy and stability of the proposed method are verified through illustrative numerical tests.

• 한국전산구조공학회논문집
• /
• 제18권2호
• /
• pp.113-121
• /
• 2005

무요소기법이 공통적으로 내재하고 있는 수치적분의 부정확성을 해결하기 위해, 페트로프-갤러킨 자연요소법이라 불리는 향상된 자연요소법을 제안한다. 제안된 방법은 라플라스 기저함수를 시도 형상함수로 사용하는 반면, 시험 형상함수로서 델라우니 삼각형이 지지영역이 되는 함수를 새롭게 정의한다. 이러한 접근은 통상적인 적분영역과 적분함수 지지영역간의 불일치를 제거하게 하며, 이는 적용이 편리할 뿐만 아니라 수치적분의 정확성을 보장한다 본 논문에서는 2차윈 선형 탄성의 대표적인 검증문제를 통하여 제안된 방법의 타당성을 검증한다. 비교를 위해 기존의 부브노프-갤러킨 자연요소법과 일정 변형률 유한요소법을 이용한 해석을 동시에 수행한다. 조각 시험과 수렴율 평가를 통해 제안된 기법의 우수성을 확인할 수 있다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.289-294
• /
• 2010

Numerical stud of an oscillating body in incompressible fluid is performed. Stabilized finite element method comprising of Streamline-Upwind/Petrov-Galerkin (SUPG) and Pressure-Stabilizing/Petrov-Galerkin (PSPG) formulations of linear triangular elements was employed to solve 2D incompressible Navier-Stokes equations whereas the motion of the body was considered by incorporating the arbitrary Langrangian-Eulerian(ALE) formulation. An algebraic moving mesh strategy is utilized for obtaining body conforming mesh deformation at each time step. Two tests cases, namely motion of a circular cylinder and of an airfoil in incompressible flow were analyzed. The model is first validated against the stationary cases and then the capability to handle moving boundaries is demonstrated.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집E
• /
• pp.386-391
• /
• 2001

A finite element code based on P2P1 tetra element has been developed for the large eddy simulation (LES) of turbulent flows around a complex geometry. Fractional 4-step algorithm is employed to obtain time accurate solution since it is less expensive than the integrated formulation, in which the velocity and pressure fields are solved at the same time. Crank-Nicolson method is used for second order temporal discretization and Galerkin method is adopted for spatial discretization. For very high Reynolds number flows, which would require a formidable number of nodes to resolve the flow field, SUPG (Streamline Upwind Petrov-Galerkin) method is applied to the quadratic interpolation function for velocity variables, Noting that the calculation of intrinsic time scale is very complicated when using SUPG for quadratic tetra element of velocity variables, the present study uses a unique intrinsic time scale proposed by Codina et al. since it makes the present three-dimensional unstructured code much simpler in terms of implementing SUPG. In order to see the effect of numerical diffusion caused by using an upwind scheme (SUPG), those obtained from P2P1 Galerkin method and P2P1 Petrov-Galerkin approach are compared for the flow around a sphere at some Reynolds number. Smagorinsky model is adopted as subgrid scale models in the context of P2P1 finite element method. As a benchmark problem for code validation, turbulent flows around a sphere and a MIRA model have been studied at various Reynolds numbers.