• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• 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.

• The Journal of the Acoustical Society of Korea
• /
• v.8 no.5
• /
• pp.51-58
• /
• 1989

The finite element method is usually formulated by utilizing the variation principle. In this paper, we introduce the approximate equation of finite element from Helmholtz eqation by means of the Galerkin method, which provides the best approximation of those methods known as the method of weighted residuals, and a numerical simulation based of the finite element method is applied to analysing the acoustic modes and the pattern of sound radiation in two and three dimensional sound fields. Beside the numerical calculations, the acoustic modes and the sound pressure level are mesured by scale model experiments. The finite element analysis of the model shows very good agreement with the mesured results.

• Proceedings of the Optical Society of Korea Conference
• /
• 1989.02a
• /
• 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
• /
• v.40 no.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
• /
• v.2 no.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.

• Journal of Navigation and Port Research
• /
• v.31 no.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.

• Water for future
• /
• v.14 no.2
• /
• pp.53-62
• /
• 1981

In the Hydraulic Structure, Such as dam body or levee of river that is constructed with soil, We analyzed a top line of free ground water table. This study is based on the logical reason that the pressure on the free surface is atmospheric and the seepage line is a stream line. In order to research for the unknown seepage line. We analyzed seepage water of steady flow through parous media by Finite Element method based on Galerkin Principle, and compared the comluted value with experimental value. The results show that the computed value was nearly equal to the experimental value. Finally, it noticed that finite Element method was more practical than Experimental Method for Seepage line analysis.

• Nuclear Engineering and Technology
• /
• v.48 no.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
• /
• v.17 no.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.

• Journal of the Korean Society of Safety
• /
• v.8 no.4
• /
• pp.201-206
• /
• 1993

A space-time finite element method was presented for time dependent problem. The method which treat both the space and time unformly were proposed and numerically tested. The weighted residual process was used to formulate a finite element method in a space-time domain based upon continuous Galerkin method. This method leads to a conditional stabie high-order accurate solver.