• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.21-26
• /
• 2006

A large-scale finite element simulation and modeling method is presented for environmental flows in urban area. Parallel stabilized finite element method based on domain decomposition method is employed for the numerical simulation. Several GIS and CAD data are used for the preparation of the shape model for landform and urban structures. The present method Is applied to the simulation of flood flow and wind flow In urban area. The present method is shown to be a useful planning and design tool for the natural disasters and the change of environments in urban area.

• Journal of applied mathematics & informatics
• /
• v.15 no.1_2
• /
• pp.1-28
• /
• 2004

We first apply a first order splitting to a semilinear reaction-diffusion equation and then discretize the resulting system by an $H^1$-Galerkin mixed finite element method in space. This semidiscrete method yields a system of differential algebraic equations (DAEs) of index one. A priori error estimates for semidiscrete scheme are derived for both differ-ential as well as algebraic components. For fully discretization, an implicit Runge-Kutta (IRK) methods is applied to the temporal direction and the error estimates are discussed for both components. Finally, we conclude the paper with a numerical example.

• Journal of applied mathematics & informatics
• /
• v.15 no.1_2
• /
• pp.29-51
• /
• 2004

In this paper, we propose a new mixed finite element method, called the characteristics-mixed method, for approximating the solution to Burgers' equation. This method is based upon a space-time variational form of Burgers' equation. The hyperbolic part of the equation is approximated along the characteristics in time and the diffusion part is approximated by a mixed finite element method of lowest order. The scheme is locally conservative since fluid is transported along the approximate characteristics on the discrete level and the test function can be piecewise constant. Our analysis show the new method approximate the scalar unknown and the vector flux optimally and simultaneously. We also show this scheme has much smaller time-truncation errors than those of standard methods. Numerical example is presented to show that the new scheme is easily implemented, shocks and boundary layers are handled with almost no oscillations. One of the contributions of the paper is to show how the optimal error estimates in $L^2(\Omega)$ are obtained which are much more difficult than in the standard finite element methods. These results seem to be new in the literature of finite element methods.

• Wind and Structures
• /
• v.26 no.1
• /
• pp.25-34
• /
• 2018

This paper presents a beam finite element model of a vibrate wind blade in large elastic deformation subjected to the aerodynamic, centrifugal, gyroscopic and gravity loads. The gyroscopic loads applied to the blade are induced by her simultaneous vibration and rotation. The proposed beam finite element model is based on a simplex interpolation method and it is mainly intended to the numerical analysis of wind blades vibration in large elastic deformation. For this purpose, the theory of the sheared beams and the finite element method are combined to develop the algebraic equations system governing the three-dimensional motion of blade vibration. The applicability of the theoretical approach is elucidated through an original case study. Also, the static deformation of the used wind blade is assessed by appropriate software using a solid finite element model in order to show the effectiveness of the obtained results. To simulate the nonlinear dynamic response of wind blade, the predictor-corrector Newmark scheme is applied and the stability of numerical process is approved during a large time of blade functioning. Finally, the influence of the modified geometrical stiffness on the amplitudes and frequencies of the wind blade vibration induced by the sinusoidal excitation of gravity is analyzed.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.3
• /
• pp.243-259
• /
• 2016

The Richards equation for water movement in unsaturated soil is highly nonlinear partial differential equations which are not solvable analytically unless unrealistic and oversimplifying assumptions are made regarding the attributes, dynamics, and properties of the physical systems. Therefore, conventionally, numerical solutions are the only feasible procedures to model flow in partially saturated porous media. The standard Finite element numerical technique is usually coupled with an Euler time discretizations scheme. Except for the fully explicit forward method, any other Euler time-marching algorithm generates nonlinear algebraic equations which should be solved using iterative procedures such as Newton and Picard iterations. In this study, lumped mass and distributed mass in the frame of Picard and Newton iterative techniques were evaluated to determine the most efficient method to solve the Richards equation with finite element model. The accuracy and computational efficiency of the scheme and of the Picard and Newton models are assessed for three test problems simulating one-dimensional flow processes in unsaturated porous media. Results demonstrated that, the conventional mass distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Even though small mesh sizes are applied for all the test problems, it is shown that the traditional mass-distributed scheme can still generate an incorrect response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. On the other hand, non oscillatory solutions are obtained and non-physics solutions for these problems are evaded by using the mass-lumped finite element method.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.11b
• /
• pp.85-89
• /
• 2005

For the analysis of a vibrating two dimensional structure such as the simply supported rectangular plate, Spectral Finite Element Method (SFEM) has been studied. Under the condition that two parallel edges are simply supported at least and the other two edges can be arbitrary, Spectral Finite Element has been developed. Using this element SFEM is applied to the vibrating rectangular plate which all edges are simply supported, and obtain the frequency response function in frequency domain and the dynamic response in time domain. To evaluate these results normal mode method and finite element method (FEM) are also accomplished and compared. It is seen that SFEM is more powerful analysis tool than FEM in high frequency range.

• Journal of Power System Engineering
• /
• v.16 no.1
• /
• pp.78-83
• /
• 2012

The authors suggest the algorithm for the static analysis of a three dimensional solid structure by using the finite element-transfer stiffness coefficient method (FE-TSCM) and the hexahedral element of the finite element method (FEM). MATLAB codes were made by both FE-TSCM and FEM for the static analysis of three dimensional solid structure. They were applied to the static analyses of a very thick plate structure and a three dimensional solid structure. In this paper, as we compare the results of FE-TSCM with those of FEM, we confirm that FE-TSCM introducing the hexahedral element for the static analysis of a three dimensional solid structure is very effective from the viewpoint of the computational accuracy, speed, and storage.

• Journal of Mechanical Science and Technology
• /
• v.15 no.11
• /
• pp.1499-1506
• /
• 2001

Various transition elements are used in general for the effective finite element analysis of complicated mechanical structures. In this paper, a solid-to-beam transition finite element, which can b e used for connecting a C1-continuity beam element to a continuum solid element, is proposed. The shape functions of the transition finite element are derived to meet the compatibility condition, and a transition element equation is formulated by the conventional finite element procedure. In order to show the effectiveness and convergence characteristics of the proposed transition element, numerical tests are performed for various examples. As a result of this study, following conclusions are obtained. (1) The proposed transition element, which meets the compatibility of the primary variables, exhibits excellent accuracy. (2) In case of using the proposed transition element, the number of nodes in the finite element model may be considerably reduced and the model construction becomes more convenient. (3) This formulation method can be applied to the usage of higher order elements.

• 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.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.24 no.2
• /
• pp.215-227
• /
• 2020

A new multiscale finite element method for elliptic problems with highly oscillating coefficients are introduced. A hybridization yields a locally flux-conserving numerical scheme for multiscale problems. Our approach naturally induces a homogenized equation which facilitates error analysis. Complete convergence analysis is given and numerical examples are presented to validate our analysis.