• 대한수학회지
• /
• 제49권4호
• /
• pp.765-778
• /
• 2012

We construct a symmetric finite volume element (SFVE) scheme for a self-adjoint elliptic problem on tetrahedron grids and prove that our new scheme has optimal convergent order for the solution and has superconvergent order for the flux when grids are quasi-uniform and regular. The symmetry of our scheme is helpful to solve efficiently the corresponding discrete system. Numerical experiments are carried out to confirm the theoretical results.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제12권4호
• /
• pp.249-260
• /
• 2008

In this work we propose a mixed finite volume method for the Signorini problem which are based on the idea of Keller's finite volume box method. The triangulation may consist of both triangles and quadrilaterals. We choose the first-order nonconforming space for the scalar approximation and the lowest-order Raviart-Thomas vector space for the vector approximation. It will be shown that our mixed finite volume method is equivalent to the standard nonconforming finite element method for the scalar variable with a slightly modified right-hand side, which are crucially used in a priori error analysis.

• 대한의용생체공학회:의공학회지
• /
• 제28권2호
• /
• pp.174-177
• /
• 2007

Trabecular bone can be accurately represented using image-based finite element modeling and analysis of these bone models is widely used to predict their mechanical properties. However, the choice of thresholding technique, a necessary step in converting grayscale images to finite element models which can thus significantly influence the structure of the resulting finite element model, is often overlooked. Therefore, we investigated the effects of thresholding techniques on micro-computed tomography (micro-CT) based finite element models of trabecular bone. Three types of thresholding techniques were applied to micro-CT images of trabecular bone which resulted in three unique finite element models for each specimen. Bone volume fractions and apparent moduli were predicted for each model and compared to experimental results. Our findings suggest that predictions of apparent properties agree well with experimental measurements regardless of the choice of thresholding technique in micro CT images of trabecular bone.

• 한국복합재료학회:학술대회논문집
• /
• 한국복합재료학회 2005년도 추계학술발표대회 논문집
• /
• pp.66-69
• /
• 2005

Mechanical strength of functionally graded composite plates that composed of ceramic, functionally graded material and metal layers is investigated using 3-D finite element method. In FGM layer, material properties are assumed to be varied continuously in the thickness direction according to a simple power law distribution in terms of the volume fraction of a ceramic and metal. The 3-D finite element model is adopted by using an IS-node solid element to analyze more accurately the variation of material properties in the thickness direction. Numerical results are compared with those of the previous works. In addition, the displacements, the tensile stresses and the compressive stresses are analyzed for the variation of FGM thickness ratio and volume fraction distribution.Mechanical strength of functionally graded composite plates that composed of ceramic, functionally graded material and metal layers is investigated using 3-D finite element method. In FGM layer, material properties are assumed to be varied continuously in the thickness direction according to a simple power law distribution in terms of the volume fraction of a ceramic and metal. The 3-D finite element model is adopted by using an IS-node solid element to analyze more accurately the variation of material properties in the thickness direction. Numerical results are compared with those of the previous works. In addition, the displacements, the tensile stresses and the compressive stresses are analyzed for the variation of FGM thickness ratio and volume fraction distribution.

• Structural Engineering and Mechanics
• /
• 제80권5호
• /
• pp.553-562
• /
• 2021

In this paper, a hybrid scaled boundary finite element and finite volume method (SBFEM-FVM) is proposed for simulating hydraulic-fracture propagation in brittle concrete materials. As a semi-analytical method, the scaled boundary finite element method is introduced for modelling concrete crack propagation under both an external force and water pressure. The finite volume method is employed to model the water within the crack and consider the relationship between the water pressure and the crack opening distance. The cohesive crack model is used to analyse the non-linear fracture process zone. The numerical results are compared with experimental data, indicating that the F-CMOD curves and water pressure changes under different loading conditions are approximately the same. Different types of water pressure distributions are also studied with the proposed coupled model, and the results show that the internal water pressure distribution has an important influence on crack propagation.

• 대한용접접합학회:학술대회논문집
• /
• 대한용접접합학회 2003년도 추계학술발표대회 개요집
• /
• pp.36-38
• /
• 2003

In order to understand the basic knowledge on the model of heat source in hybrid welding, authors have conducted finite element analysis to calculate heat distribution using three heat source models of non-split type and split type(Volume, Volume-Volume, Volume-Surface). From the research result, we can confirm that Volume-Volume heat source of split type is suitable for the analysis of heat distribution.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 추계학술대회논문집B
• /
• pp.551-556
• /
• 2001

The filling pattern technique based on the finite element method and Eulerian mesh advancement approach has been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The penalty and predictor-corrector methods are used effectively for finite element formulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among four filling patterns at each triangular control volume. Using the proposed numerical technique, the collapse of a dam has been analyzed to predict flow phenomenon of fluid and the predicted front positions versus time have been compared with the reported experimental result.

• 대한기계학회논문집B
• /
• 제21권9호
• /
• pp.1230-1243
• /
• 1997

A numerical technique for simulating incompressible viscous flow with free surface is presented. The flow field is obtained by penalty finite element formulation. In this work, a modified volume of fluid (VOF) method which is compatible with 4-node element is proposed to track the moving free surface. This scheme can be applied to irregular mesh system, and can be easily extended to three dimensional geometries. Numerical analyses were done for two benchmark examples, namely the broken dam problem and the solitary wave propagation problem. The numerical results were in close agreement with the existing data. Illustrative examples were studied to show the effectiveness of the proposed numerical scheme.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 춘계학술대회논문집B
• /
• pp.51-56
• /
• 2000

The Finite Element Solutions Is reported on solid-liquid phase change in porous media with natural convection including freezing. The model is based on volume averaged transport equations, while phase change is assumed to occur over a small temperature range. The FEM (Finite Element Method) algorithm used in this study is 3-step time-splitting method which requires much less execution time and computer storage the velocity-pressure integrated method and the penalty method. And the explicit Lax-Wendroff scheme is applied to nonlinear convective term in the energy equation. For natural convection including melting and solidification the numerical results show reasonable agreement with FDM (Finite Difference Method) results.

• 대한수학회지
• /
• 제53권6호
• /
• pp.1411-1429
• /
• 2016

We consider the nite volume element method for elliptic problems using isoparametric bilinear elements on quadrilateral meshes. A gradient recovery method is presented by using the patch interpolation technique. Based on some superclose estimates, we prove that the recovered gradient $R({\nabla}u_h)$ possesses the superconvergence: ${\parallel}{\nabla}u-R({\nabla}u_h){\parallel}=O(h^2){\parallel}u{\parallel}_3$. Finally, some numerical examples are provided to illustrate our theoretical analysis.