• 대한기계학회:학술대회논문집
• /
• 대한기계학회 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.

• 한국음향학회지
• /
• 제8권5호
• /
• pp.51-58
• /
• 1989

유한요소법은 일반적으로 변분원리를 이용해 정식화를 하고 있으나, 본 연구에서는 웨이티드 잔차법으로 아주 좋은 근사해를 얻을 수 있다는 Galerkin법에 의해 Helmholtz방정식으로부터 직접 유산요소 정식화하는 방법을 소개하고, 정식화한 수치계산법을 2, 3차원 음장의 고유모드 및 음향방사상태해석에 응용하였다. 또한 수치 계산결과를 확인하기 위하여 간단한 모형을 제작, 실내음향 모드와 음압분포 등의 측정도 병행하였으며 그 결과, 유한요소법에 의한 수치해석결과의 측정치가 잘 맞는 것을 알았다.

• 한국광학회:학술대회논문집
• /
• 한국광학회 1989년도 제4회 파동 및 레이저 학술발표회 4th Conference on Waves and lasers 논문집 - 한국광학회
• /
• pp.157-160
• /
• 1989

The most serious difficulty in using the finite element method is the appearance of the so-called spurious, nonphysical modes. We have proposed the finite element formulation of the variational expression in the three-component magnetic field based on Galerkin's method. In this approach, the divergence relation H is satisfied and spurious modes does not appear and finite-element solutions agree with the exact solutions.

• Journal of applied mathematics & informatics
• /
• 제40권3_4호
• /
• pp.603-618
• /
• 2022

In this paper, a numerical solution of the generalized diffusion equation with a delay has been obtained by a numerical technique based on the Galerkin finite element method by applying the cubic B-spline basis functions. The time discretization process is carried out using the forward Euler method. The numerical scheme is required to preserve the delay-independent asymptotic stability with an additional restriction on time and spatial step sizes. Both the theoretical and computational rates of convergence of the numerical method have been examined and found to be in agreement. As it can be observed from the numerical results given in tables and graphs, the proposed method approximates the exact solution very well. The accuracy of the numerical scheme is confirmed by computing L₂ and L_∞ error norms.

• International Journal of Precision Engineering and Manufacturing
• /
• 제2권2호
• /
• pp.5-10
• /
• 2001

This paper attempt to explore the shape of stress wave propagation of 3-dimensional stress field which in made in the process of the time increment. A finite element program about 3-dimensional stress wave propagation is developed for investigating the changing shape of the stress by the impact load. The finite element program, which is the solution for the 3-dimensional stress wave analysis, based on Galerkin and Newmark-${\beta}$ method at time increment step. The tensile stress and compressive stress become larger with the order of the middle , the upper and the opposite layers when the impact load is applied. In a while the shear stress become larger according to the order of the upper, the middle and the opposite layers when impact load applied.

• 한국항해항만학회지
• /
• 제31권9호
• /
• pp.793-799
• /
• 2007

The effect of seasonal wind on the tidal circulation in Jeju harbor was examined by using a numerical shallow water model. A finite element for analyzing shallow water flow is presented. The Galerkin method is employed for spatial discretization. Two step explicit finite element scheme is used to discretize the time function, which has advantage in problems treating large numbers of elements and unsteady state. The numerical simulation is compared with three cases; Case 1 does not consider the effect of wind, Case 2 and Case 3 consider the effect of summer and winter seasonal wind, respectively. According to result considering effect of seasonal wind, velocity of current vector shows slightly stronger than that of case 1 in the flow field. It can be concluded that the present method is a useful and effective tool in tidal current analysis.

• 물과 미래
• /
• 제14권2호
• /
• pp.53-62
• /
• 1981

댐이나 제방과 같이 흙으로 축조된 수리구조물에 있어서, 자유지하수면의 최상부 침투선을 해석하였다. 자유지하수면에 작용하는 압력은 대기압이고, 침투선은 유선이라는 원리에 의하여 연구를 수행하였다. 미지의 침투선을 해석하기 위하여 Galerkin 원리에 기초를 둔 유한요소법에 의하여 다공체속을 흐르는 정류상태의 침투수를 해석하여 실험치와 이론치를 비교하였고 그 결과 이론치와 실험치가 거의 일치함을 알았다. 결론적으로 침투선해석에 있어서 유한요소법이 실험적인 방법보다 더 실용적이라는 것을 알았다.

• Nuclear Engineering and Technology
• /
• 제48권1호
• /
• pp.43-54
• /
• 2016

In the present paper, development of the three-dimensional (3D) computational code based on Galerkin finite element method (GFEM) for solving the multigroup forward/adjoint diffusion equation in both rectangular and hexagonal geometries is reported. Linear approximation of shape functions in the GFEM with unstructured tetrahedron elements is used in the calculation. Both criticality and fixed source calculations may be performed using the developed GFEM-3D computational code. An acceptable level of accuracy at a low computational cost is the main advantage of applying the unstructured tetrahedron elements. The unstructured tetrahedron elements generated with Gambit software are used in the GFEM-3D computational code through a developed interface. The forward/adjoint multiplication factor, forward/adjoint flux distribution, and power distribution in the reactor core are calculated using the power iteration method. Criticality calculations are benchmarked against the valid solution of the neutron diffusion equation for International Atomic Energy Agency (IAEA)-3D and Water-Water Energetic Reactor (VVER)-1000 reactor cores. In addition, validation of the calculations against the $P_1$ approximation of the transport theory is investigated in relation to the liquid metal fast breeder reactor benchmark problem. The neutron fixed source calculations are benchmarked through a comparison with the results obtained from similar computational codes. Finally, an analysis of the sensitivity of calculations to the number of elements is performed.

• Journal of Mechanical Science and Technology
• /
• 제17권1호
• /
• pp.64-76
• /
• 2003

In this work, a new penalty formulation is proposed for the analysis of Mindlin-Reissner plates by using the element-free Galerkin method. A penalized weak form for the Mindlin-Reissner Plates is constructed through the exterior penalty method to enforce the essential boundary conditions of rotations as well as transverse displacements. In the numerical examples, some typical problems of Mindlin-Reissner plates are analyzed, and parametric studies on the order of integration and the size of influence domain are also carried out. The effect of the types of background cells on the accuracy of numerical solutions is observed and a proper type of background cell for obtaining optimal accuracy is suggested. Further, optimal order of integration and basis order of Moving Least Squares approximation are suggested to efficiently handle the irregularly distributed nodes through the triangular type of background cells. From the numerical tests, it is identified that unlike the finite element method, the proposed element-free Galerkin method with penalty technique gives highly accurate solution without shear locking in dealing with Mindlin-Reissner plates.

• 한국안전학회지
• /
• 제8권4호
• /
• pp.201-206
• /
• 1993

동적 문제에 대한 공간-시간 유한요소법을 제시하였다. 이 방법은 공간과 시간을 동일한 변수로 취급하였으며 공간-시간 영역에서의 유한요소 전개에 있어서는 연속적 갤러킨 방법에 근거하여 가중여분법을 이용하였다. 이 방법은 조건부 안정을 주는 고차원적 정확성을 주는 해법인 것이다.