• 제목/요약/키워드: Petrov-Galerkin Method

검색결과 48건 처리시간 0.021초

An Adaptive Mesh-Independent Numerical Integration for Meshless Local Petrov-Galerkin Method

  • Cho, Jin-Yeon;Jee, Young-Burm
    • Journal of Mechanical Science and Technology
    • /
    • 제17권7호
    • /
    • pp.986-998
    • /
    • 2003
  • In this paper, an adaptive numerical integration scheme, which does not need non-overlapping and contiguous integration meshes, is proposed for the MLPG (Meshless Local Petrov-Galerkin) method. In the proposed algorithm, the integration points are located between the neighboring nodes to properly consider the irregular nodal distribution, and the nodal points are also included as integration points. For numerical integration without well-defined meshes, the Shepard shape function is adopted to approximate the integrand in the local symmetric weak form, by the values of the integrand at the integration points. This procedure makes it possible to integrate the local symmetric weak form without any integration meshes (non-overlapping and contiguous integration domains). The convergence tests are performed, to investigate the present scheme and several numerical examples are analyzed by using the proposed scheme.

GIS와 연계한 한강 하류부에서의 동수역학적 수치모형의 개발 (Development of Hydrodynamic Model on the Downstream of Han River by Using Geographical Information System)

  • 한건연;이을래;박재홍
    • 한국수자원학회논문집
    • /
    • 제34권2호
    • /
    • pp.107-118
    • /
    • 2001
  • 본 연구에서는 2차원 동수역학적 수치해석모형을 개발하고, 한강하류부에 적용하였다. 과학적인 천수흐름 거동해석을 위해서 GIS에 기초한 등수역학모형을 구성하였으며, 개발된 모형의 검증을 위해서 모의수행에 의한 결과 값과 현장관측값과의 비교를 수행하였다. 본 연구에서 제시된 수치해석 모형은 Petrol-Galerkin 유한요소법이다. 또한 GIS를 이용한 하천흐름해석기법을 관측 하였으며, 이는 다양한 해석을 통해 직접적으로 사용자의 결정에 도움을 주게 된다. 본 연구의 결과는 다양한 하천, 저수지 그리고 하구의 해석에 관한 지침이 될 수 있을 것으로 사료된다.

  • PDF

안정화된 유한요소법을 이용한 진동하는 2차원 물체 주의 유동해석 (A STABILZED FINITE ELEMENT COMPUTATION OF FLOW AROUND OSCILLATING 2D BODIES)

  • 안형택;라술 라힐
    • 한국전산유체공학회:학술대회논문집
    • /
    • 한국전산유체공학회 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.

  • PDF

Elastic analysis of arbitrary shape plates using Meshless local Petrov-Galerkin method

  • Edalati, H.;Soltani, B.
    • Wind and Structures
    • /
    • 제27권4호
    • /
    • pp.235-245
    • /
    • 2018
  • In this study the stress analysis of orthotropic thin plate with arbitrary shapes for different boundary conditionsis investigated. Meshfreemethod is applied to static analysis of thin plates with various geometries based on the Kirchhoff classical plate theory. According to the meshfree method the domain of the plates are expressed through a set of nodes without using mesh. In this method, a set of nodes are defined in a standard rectangular domain, then via a third order map, these nodes are transferred to the main domain of the original geometry; therefore the analysis of the plates can be done. Herein, Meshless local Petrov-Galerkin (MLPG) as a meshfree numerical method is utilized. The MLS function in MLPG does not satisfy essential boundary conditions using Delta Kronecker. In the MLPG method, direct interpolation of the boundary conditions can be applied due to constructing node by node of the system equations. The detailed parametric study is conducted, focusing on the arbitrary geometries of the thin plates. Results show that the meshfree method provides better accuracy rather than finite element method. Also, it is found that trend of the figures have good agreement with relevant published papers.

P2P1 유한요소를 이용한 LES (Large Eddy simulation using P2P1 finite element formulation)

  • 최형권;남영석;유정열
    • 대한기계학회:학술대회논문집
    • /
    • 대한기계학회 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.

  • PDF

부정류 해석을 위한 Hermite-Collocation 기법의 적용 (Application of Hermite-Collocation Method for Unsteady Flow Analysis)

  • 한건연;이을래;박재홍
    • 한국수자원학회논문집
    • /
    • 제32권3호
    • /
    • pp.237-246
    • /
    • 1999
  • 자유수면 흐름의 모의를 위한 유한요소모형이 동수역학적 흐름방정식과 collocation 유한요소법에 의해 모의하였다. collocation 기법은 Hermite 다항식을 가진 접합점에서 적용이 되며, 메크릭스 방정식은 skyline 기법에 의해 해석하였다. 본 연구 모형은 마찰이 없는 수평수로에서의 정상도수, 비선형 표면전파 그리고 댐 파괴해석에 적용하였다. 계산결과 Bubnov-Galerkin 과 Petrov-Galerkin 기법과 비교하였다. 실제하천에 대한 적용성을 검토하기 위해서 북한강 유역에 적용하여 해석하였는데, 계산결과는 유량수문곡선에 있어서 기존의 DWOPER 모형의 결과와 일치하였다. Collocation 기법은 개수로 흐름에서의 점변 및 급변 부정류흐름을 모의하기 위해서 적절한 기법임을 확인할 수 있었다.

  • PDF

Meshless Local Petrov-Galerkin (MLPG) method for dynamic analysis of non-symmetric nanocomposite cylindrical shell

  • Ferezghi, Yaser Sadeghi;Sohrabi, Mohamadreza;Nezhad, Seyed Mojtaba Mosavi
    • Structural Engineering and Mechanics
    • /
    • 제74권5호
    • /
    • pp.679-698
    • /
    • 2020
  • In this paper, the meshless local Petrov-Galerkin (MLPG) method is developed for dynamic analysis of non-symmetric nanocomposite cylindrical shell equations of elastic wave motion with nonlinear grading patterns under shock loading. The mechanical properties of the nanocomposite cylinder are obtained based on a micro-mechanical model. In this study, four kinds of grading patterns are assumed for carbon nanotube mechanical properties. The displacements can be approximated using shape function so, the multiquadrics (MQ) Radial Basis Functions (RBF) are used as the shape function. In order to discretize the derived equations in time domains, the Newmark time approximation scheme with suitable time step is used. To demonstrate the accuracy of the present method for dynamic analysis, at the first a problem verifies with analytical solution and then the present method compares with the finite element method (FEM), finally, the present method verifies by using the element free Galerkin (EFG) method. The comparison shows the high capacity and accuracy of the present method in the dynamic analysis of cylindrical shells. The capability of the present method to dynamic analysis of non-symmetric nanocomposite cylindrical shell is demonstrated by dynamic analysis of the cylinder with different kinds of grading patterns and angle of nanocomposite reinforcements. The present method shows high accuracy, efficiency and capability to dynamic analysis of non-symmetric nanocomposite cylindrical shell, which it furnishes a ground for a more flexible design.

NUMERICAL MODELING OF TWO-DIMENSIONAL ADVECTION-DISPERSION IN OPEN CHANNEL

  • Lee, Myung-Eun;Kim, Young-Han;Seo, Il-Won
    • Water Engineering Research
    • /
    • 제4권1호
    • /
    • pp.45-58
    • /
    • 2003
  • Two-dimensional depth-averaged advection-dispersion equation was simulated using FEM. In the straight rectangular channel, the advection-dispersion processes are simulated so that these results can be compared with analyti-cal solutions for the transverse line injection and the point injection. In the straight domain the standard Galerkin method with the linear basis function is found to be inadequate to the advection-dispersion analysis compared to the upwind finite element scheme. The experimental data in the S-curved channel were compared with the result by the numerical model using SUPG(Streamline upwind Petrov-Galerkin) method.

  • PDF

Stabilized finite element technique and its application for turbulent flow with high Reynolds number

  • Huang, Cheng;Yan, Bao;Zhou, Dai;Xu, Jinquan
    • Wind and Structures
    • /
    • 제14권5호
    • /
    • pp.465-480
    • /
    • 2011
  • In this paper, a stabilized large eddy simulation technique is developed to predict turbulent flow with high Reynolds number. Streamline Upwind Petrov-Galerkin (SUPG) stabilized method and three-step technique are both implemented for the finite element formulation of Smagorinsky sub-grid scale (SGS) model. Temporal discretization is performed using three-step technique with viscous term treated implicitly. And the pressure is computed from Poisson equation derived from the incompressible condition. Then two numerical examples of turbulent flow with high Reynolds number are discussed. One is lid driven flow at Re = $10^5$ in a triangular cavity, the other is turbulent flow past a square cylinder at Re = 22000. Results show that the present technique can effectively suppress the instabilities of turbulent flow caused by traditional FEM and well predict the unsteady flow even with coarse mesh.