• 제목/요약/키워드: Galerkin methods

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

페트로프-갤러킨 자연요소법 : II. 선형 정탄성 해석 (The Petrov-Galerkin Natural Element Method : II. Linear Elastostatic Analysis)

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

전자파 수치 해석을 위해 갤러킨 기법과 보간법을 혼용하여 개선시킨 모멘트법 (Improved Method of Moments Using Hybrid Technique of Galerkin's and Interpolation Methods for Numerical Analysis of Electromagnetic Waves)

  • 황지환;권순구;오이석
    • 한국전자파학회논문지
    • /
    • 제23권4호
    • /
    • pp.541-550
    • /
    • 2012
  • 본 논문에서는 3차원 공간의 전자파 수치 해석을 위한 모멘트법(method of moments)의 개선된 해석 기법을 선보인다. 전자파 산란 특성을 해석하기 위해 기본적으로 EFIE(Electric Field Integral Equation)와 RWG(Rao-Wilton-Glisson) 기저 함수를 이용하였으며, 계산 효율을 높이기 위해 기존의 갤러킨(Galerkin) 기법과 중심점 보간(interpolation)법을 혼용하여 해석 시간을 단축시켰다. 이때, 계산 정확도 유지를 위해 임피던스 행렬의 각 원소간 거리를 상대 거리 지수로 정의하여 보간법 적용이 가능한 먼 거리 원소를 구분하였다. 제안된 해석 기법의 성능 검증은 금속구의 Mie-series 해법을 이용한 이론적 RCS(Radar Cross Section)를 비교/분석하였다. 또한, 본 연구 결과를 삼면-/전방향- 전파반사기와 같은 산란체에 적용하여 레이더 후방 산란 특성을 분석하였다.

고차 정확도 수치기법의 GPU 계산을 통한 효율적인 압축성 유동 해석 (EFFICIENT COMPUTATION OF COMPRESSIBLE FLOW BY HIGHER-ORDER METHOD ACCELERATED USING GPU)

  • 장태규;박진석;김종암
    • 한국전산유체공학회지
    • /
    • 제19권3호
    • /
    • pp.52-61
    • /
    • 2014
  • The present paper deals with the efficient computation of higher-order CFD methods for compressible flow using graphics processing units (GPU). The higher-order CFD methods, such as discontinuous Galerkin (DG) methods and correction procedure via reconstruction (CPR) methods, can realize arbitrary higher-order accuracy with compact stencil on unstructured mesh. However, they require much more computational costs compared to the widely used finite volume methods (FVM). Graphics processing unit, consisting of hundreds or thousands small cores, is apt to massive parallel computations of compressible flow based on the higher-order CFD methods and can reduce computational time greatly. Higher-order multi-dimensional limiting process (MLP) is applied for the robust control of numerical oscillations around shock discontinuity and implemented efficiently on GPU. The program is written and optimized in CUDA library offered from NVIDIA. The whole algorithms are implemented to guarantee accurate and efficient computations for parallel programming on shared-memory model of GPU. The extensive numerical experiments validates that the GPU successfully accelerates computing compressible flow using higher-order method.

최소 제곱 무요소법과 적분 오차 (Least-Squares Meshfree Method and Integration Error)

  • 박상훈;윤성기
    • 대한기계학회논문집A
    • /
    • 제25권10호
    • /
    • pp.1605-1612
    • /
    • 2001
  • Least-squares meshfree method is presented. Conventional meshfree methods based on the Galerkin formulation suffer from inaccurate numerical integration. Least-squares formulation exhibits rather different integration-related characteristics. It is demonstrated through numerical examples that least-squares formulation is much more robust to integration errors than the Galerkin's. Therefore efficient meshfree methods can be devised by combining very simple integration algorithms and least-squares formulation.

Prediction of initiation time of corrosion in RC using meshless methods

  • Yao, Ling;Zhang, Lingling;Zhang, Ling;Li, Xiaolu
    • Computers and Concrete
    • /
    • 제16권5호
    • /
    • pp.669-682
    • /
    • 2015
  • Degradation of reinforced concrete (RC) structures due to chloride penetration followed by reinforcement corrosion has been a serious problem in civil engineering for many years. The numerical simulation methods at present are mainly finite element method (FEM) and finite difference method (FDM), which are based on mesh. Mesh generation in engineering takes a long time. In the present article, the numerical solution of chloride transport in concrete is analyzed using radial point interpolation method (RPIM) and element-free Galerkin (EFG). They are all meshless methods. RPIM utilizes radial polynomial basis, whereas EFG uses the moving least-square approximation. A Galerkin weak form on global is used to attain the discrete equation, and four different numerical examples are presented. MQ function and appropriate parameters have been proposed in RPIM. Numerical simulation results are compared with those obtained from the finite element method (FEM) and analytical solutions. Two case of chloride transport in full saturated and unsaturated concrete are analyzed to test the practical applicability and performance of the RPIM and EFG. A good agreement is obtained among RPIM, EFG, and the experimental data. It indicates that RPIM and EFG are reliable meshless methods for prediction of chloride concentration in concrete structures.

1-D 오일러 방정식에 관한 Modal 불연속 갤러킨 기법에서의 Limiter 성능 비교 (PERFORMANCE OF LIMITERS IN MODAL DISCONTINUOUS GALERKIN METHODS FOR 1-D EULER EQUATIONS)

  • 아볼파즐 카르차니;명노신
    • 한국전산유체공학회지
    • /
    • 제21권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.

Computation of 2-D mixed-mode stress intensity factors by Petrov-Galerkin natural element method

  • Cho, Jin-Rae
    • Structural Engineering and Mechanics
    • /
    • 제56권4호
    • /
    • pp.589-603
    • /
    • 2015
  • The mixed-mode stress intensity factors of 2-D angled cracks are evaluated by Petrov-Galerkin natural element (PG-NE) method in which Voronoi polygon-based Laplace interpolation functions and CS-FE basis functions are used for the trial and test functions respectively. The interaction integral is implemented in a frame of PG-NE method in which the weighting function defined over a crack-tip integral domain is interpolated by Laplace interpolation functions. Two Cartesian coordinate systems are employed and the displacement, strains and stresses which are solved in the grid-oriented coordinate system are transformed to the other coordinate system aligned to the angled crack. The present method is validated through the numerical experiments with the angled edge and center cracks, and the numerical accuracy is examined with respect to the grid density, crack length and angle. Also, the stress intensity factors obtained by the present method are compared with other numerical methods and the exact solution. It is observed from the numerical results that the present method successfully and accurately evaluates the mixed-mode stress intensity factors of 2-D angled cracks for various crack lengths and crack angles.

Energy extraction from the motion of an oscillating water column

  • Wang, Hao;Falzarano, Jeffrey M.
    • Ocean Systems Engineering
    • /
    • 제3권4호
    • /
    • pp.327-348
    • /
    • 2013
  • An Oscillating Water Column (OWC) is a relatively practical and convenient device that converts wave energy to a usable form, which is electricity. The OWC is kept inside a fixed truncated vertical cylinder, which is a hollow structure with one open end submerged in the water and with an air turbine at the top. This research adopts potential theory and Galerkin methods to solve the fluid motion inside the OWC. Using an air-water interaction model, OWC design for energy extraction from regular wave is also explored. The hydrodynamic coefficients of the scattering and radiation potentials are solved for using the Galerkin approximation. The numerical results for the free surface elevation have been verified by a series of experiments conducted in the University of New Orleans towing tank. The effect of varying geometric parameters on the response amplitude operator (RAO) of the OWC is studied and modification of the equation for evaluating the natural frequency of the OWC is made. Using the model of air-water interaction under certain wave parameters and OWC geometric parameters, a computer program is developed to calculate the energy output from the system.

Topology Optimization for Large-displacement Compliant Mechanisms Using Element Free Galerkin Method

  • Du, Yixian;Chen, Liping
    • International Journal of CAD/CAM
    • /
    • 제8권1호
    • /
    • pp.1-10
    • /
    • 2009
  • This paper presents a topology optimization approach using element-free Galerkin method (EFGM) for the optimal design of compliant mechanisms with geometrically non-linearity. Meshless method has an advantage over the finite element method(FEM) because it is more capable of handling large deformation resulted from geometrical nonlinearity. Therefore, in this paper, EFGM is employed to discretize the governing equations and the bulk density field. The sensitivity analysis of the optimization problem is performed by incorporating the adjoint approach with the meshless method. The Lagrange multipliers method adjusted for imposition of both the concentrated and continuous essential boundary conditions in the EFGM is proposed in details. The optimization mathematical formulation is developed to convert the multi-criteria problem to an equivalent single-objective problem. The popularly applied interpolation scheme, solid isotropic material with penalization (SIMP), is used to indicate the dependence of material property upon on pseudo densities discretized to the integration points. A well studied numerical example has been applied to demonstrate the proposed approach works very well and the non-linear EFGM can obtain the better topologies than the linear EFGM to design large-displacement compliant mechanisms.