• Journal of the Korean Society of Safety
• /
• v.21 no.1 s.73
• /
• pp.15-20
• /
• 2006

The finite element alternating method is a convenient and efficient method to analyze three-dimensional cracks embedded in an infinite or a finite body because the method has the property that the uncracked body and cracks can be modeled independently. In this paper the method was applied for fatigue crack growth simulation. A surface crack in a cylinder was considered as an initial crack and the crack configurations and stress intensity factors during the crack growth were obtained. In this paper the finite element alternating method proposed by Nikishkov, Park and Atluri was used after modification. In the method, as the required solution for a crack in an infinite body, the symmetric Galerkin boundary element method formulated by Li and Mear was used. And a crack was modeled as distribution of displacement discontinuities, and the governing equation was formulated as singularity-reduced integral equations.

• Journal of Mechanical Science and Technology
• /
• v.20 no.11
• /
• pp.1881-1890
• /
• 2006

In this paper, a Petrov-Galerkin natural element method (PG-NEM) based upon the natural neighbor concept is presented for the free vibration and dynamic response analyses of two-dimensional linear elastic structures. A problem domain is discretized with a finite number of nodes and the trial basis functions are defined with the help of the Voronoi diagram. Meanwhile, the test basis functions are supported by Delaunay triangles for the accurate and easy numerical integration with the conventional Gauss quadrature rule. The numerical accuracy and stability of the proposed method are verified through illustrative numerical tests.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.414-418
• /
• 2005

A finite element model is developed for the numerical simulation of bed elevation change in a curved channel. The SU/PG (Streamline-Upwind/Petrov-Galerkin) method is used to solve 2D shallow water equations and the BG (Bubnov-Galerkin) method is used for the Exner equation. For the time derivative terms, the Crank-Nicolson scheme is used. The developed model is a decoupled model in a sense that the bed elevation does not change simultaneously with the flow during the computational time step. The total load formula with is used for the sediment transport model. The slip conditions are described along the lateral boundaries. The effects of gravity force due to geometry change and the secondary flows in a curved channel are considered in the model. For the verification, the model is applied to two laboratory experiments. The first is $140^{\circ}$ bended channel data at Delft Hydraulics Laboratory and the second is $140^{\circ}$ bended channel data at Laboratory of Fluid Mechanics of the Delft University of Technology. The finite element grid is constructed with linear quadrilateral elements. It is found that the computed results are in good agreement with measured data, showing a point bar at the inner bank and a pool at the outer bank.

• Journal of applied mathematics & informatics
• /
• v.22 no.1_2
• /
• pp.223-235
• /
• 2006

By applying the Landau-type transformation, we transform a Stefan problem with nonlinear free boundary condition into a system consisting of a parabolic equation and the ordinary differential equations. Fully discrete finite element method is developed to approximate the solution of a system of a parabolic equation and the ordinary differential equations. We derive optimal orders of convergence of fully discrete approximations in $L_2,\;H^1$ and $H^2$ normed spaces.

• Bulletin of the Korean Mathematical Society
• /
• v.31 no.2
• /
• pp.309-318
• /
• 1994

The approach presented in this paper is based on the transformation of the Stefan problem in one space dimension to an initial-boundary value problem for the heat equation in a fixed domain. Of course, the problem is non-linear. The finite element approximation adopted here is the standared continuous Galerkin method in time. In this paper, only the regular case is discussed. This means the error analysis is based on the assumption that the solution is sufficiently smooth. The aim of this paper is the existence of the solution in a finite Galerkin system of ordinary equations.

• Steel and Composite Structures
• /
• v.49 no.5
• /
• pp.587-602
• /
• 2023

This article presents the numerical modelling of transient heat transfer in highly heterogeneous composite materials where the thermal conductivity, specific heat and density are assumed to be directional-dependent. This article uses a coupled finite element-finite difference scheme to perform the transient heat transfer analysis of unidirectional (1D) and multidirectional (2D/3D) functionally graded composite panels. Here, 1D/2D/3D functionally graded structures are subjected to nonuniform heat source and inhomogeneous boundary conditions. Here, the multidirectional functionally graded materials are modelled by varying material properties in individual or in-combination of spatial directions. Here, fully spatial-dependent material properties are evaluated using Voigt's micromechanics scheme via multivariable power-law functions. The weak form is obtained through the Galerkin method and solved further via the element-space and time-step discretisation through the 2D-isoparametric finite element and the implicit backward finite difference schemes, respectively. The present model is verified by comparing it with the previously reported results and the commercially available finite element tool. The numerous illustrations confirm the significance of boundary conditions and material heterogeneity on the transient temperature responses of 1D/2D/3D functionally graded panels.

• Communications of the Korean Mathematical Society
• /
• v.18 no.4
• /
• pp.767-779
• /
• 2003

We consider a model of population dynamics whose mortality function is unbounded. We approximate the solution of the model using a discontinuous Galerkin finite element for the age variable and a backward Euler for the time variable. We present several numerical examples. It is experimentally shown that the scheme converges at the rate of $h^{3/2}$ in the case of piecewise linear polynomial space.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.124-129
• /
• 2004

Finite element alternating method have been suggested and used for assessing the integrity of cracked structures. In the paper, in order to analyze arbitrarily shaped three dimensional cracks, the finite element alternating method is extended. The cracks are modeled as a distribution of displacement discontinuities by the displacement discontinuity method and the symmetric Galerkin boundary element method. Applied the proposed method to three dimensional crack included in the elbow, the efficiency and applicability of the method were demonstrated.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.409-413
• /
• 2005

Even though the relative importance of length scale of flow system allow us to simplify three dimensional flow problem to one or two dimensional representation, many systems still require three dimensional analysis. The objective of this study is to develop an efficient and accurate finite element model for analyzing and predicting three dimensional flow features in natural rivers and to offend to model spreading of pollutants and transport of sediments in the future. Firstly, three dimensional Reynolds averaged Navier-Stokes equations with the hydrostatic pressure assumption in generalized curvilinear coordinates were combined with the kinematic free-surface condition. Secondly. to simulate realistic high Reynolds number flow, the model employed the Streamline Upwind/Petrov-Galerkin(SU/PG) scheme as a weighting function for the finite element method in conjunction with an appropriate turbulence model(Smagorinsky scheme for the horizontal plain and Mellor-Yamada scheme for the vertical direction). Several tests is performed for the purpose of validation and verification of the developed model. A simple rectangular channel, 5-shaped and U-shaped channel are used for tests and comparisons are made with RMA-10 model. Runs for each case is converged stably without a oscillation and calculated water-surface deformation, longitudinal and transversal velocities, and velocity vector fields are in good agreement with the results of RMA-10 model.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.88-97
• /
• 1996

This paper treats an adaptive finite-element method for the viscous compressible flow governed by Navier-Stokes equations in two dimensions. The numerical algorithm is the two-step Taylor-Galerkin mettled using unstructured triangular grids. To increase accuracy and stability, combined moving node method and grid refinement method have been used for grid adaption. Validation of the present algorithm has been made by comparing the present computational results with the existing experimental data and other numerical solutions. Four benchmark problems are solved for demonstration of the present numerical approach. They include a subsonic flow over a flat plate, the Carter flat plate problem, a laminar shock-boundary layer interaction. and finally a laminar flow around NACA0012 airfoil at zero angle of attack and free stream Mach number of 0.85. The results indicates that the present adaptive triangular grid method is accurate and useful for laminar viscous flow calculations.