• Communications of the Korean Mathematical Society
• /
• v.18 no.3
• /
• pp.569-580
• /
• 2003

The Lotka-McKendrick equation which describes the evolution of a single population under the phenomenological conditions is developed from the well-known Malthus’model. In this paper, we introduce the Lotka-McKendrick equation for the description of the dynamics of a population. We apply a discontinuous Galerkin finite element method in age-time domain to approximate the solution of the system. We provide some numerical results. It is experimentally shown that, when the mortality function is bounded, the scheme converges at the rate of $h^2$ in the case of piecewise linear polynomial space. It is also shown that the scheme converges at the rate of $h^{3/2}$ when the mortality function is unbounded.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.15 no.3
• /
• pp.159-166
• /
• 2003

In this study, we develop a finite element model to directly solve the Laplace equation while keeping the same computational efficiency as the models based on the extended mild-slope equation which has been widely used for calculation of wave transformation in shallow water. For this, the computational domain is discretized into finite elements with a single layer in the vertical direction. The velocity potential in the element is then expressed in terms of the potentials at the nodes located at water surface, and the Galerkin method is used to construct the numerical model. A common shape function is adopted in horizontal direction, and the cosine hyperbolic function in vertical direction, which describes the vertical behavior of progressive waves. The model was developed for vertical two-dimensional problems. In order to verify the developed model, it is applied to vertical two-dimensional problems of wave reflection and transmission. It is shown that the present finite element model is comparable to the models based on extended mild-slope equations in both computational efficiency and accuracy.

• East Asian mathematical journal
• /
• v.26 no.5
• /
• pp.615-626
• /
• 2010

In this paper, we analyze discontinuous Galerkin methods with penalty terms, namly symmetric interior penalty Galerkin methods, to solve nonlinear parabolic equations. We construct finite element spaces on which we develop fully discrete approximations using extrapolated Crank-Nicolson method. We adopt an appropriate elliptic-type projection, which leads to optimal ${\ell}^{\infty}$ ($L^2$) error estimates of discontinuous Galerkin approximations in both spatial direction and temporal direction.

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

In this paper, a singularly perturbed ordinary differential equation with non-smooth data is considered. The numerical method is generated by means of a Petrov-Galerkin finite element method with the piecewise-exponential test function and the piecewise-linear trial function. At the discontinuous point of the coefficient, a special technique is used. The method is shown to be first-order accurate and singular perturbation parameter uniform convergence. Finally, numerical results are presented, which are in agreement with theoretical results.

• Journal of applied mathematics & informatics
• /
• v.28 no.3_4
• /
• pp.953-966
• /
• 2010

In this paper, we develop discontinuous Galerkin methods with penalty terms, namaly symmetric interior penalty Galerkin methods to solve nonlinear parabolic equations. By introducing an appropriate projection of u onto finite element spaces, we prove the optimal convergence of the fully discrete discontinuous Galerkin approximations in ${\ell}^2(L^2)$ normed space.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.2
• /
• pp.25-32
• /
• 1988

A formulation of an isoparametric quadrilateral higher-order plate bending finite element is presented. The 8-noded 28-d.o.f. plate element has been degenerated from the three-dimensional continuum by introducing the plate assumptions and considering higher-order in-plane displacement profile. The element characteristics have been derived by the Galerkin's weighted residual method and computed by using the selective reduced integration technique to avoid shear-locking phenomenon. Several numerical examples are given to demonstrate the accuracy and versatility of the proposed quadrilateral higher-order plate bending element over the other existing plate finite elements in both static and dynamic analyses.

• Nuclear Engineering and Technology
• /
• v.41 no.5
• /
• pp.649-654
• /
• 2009

The finite element based lattice Boltzmann method (FELBM) has been developed to model complex fluid domain shapes, which is essential for studying fluid-structure interaction problems in commercial nuclear power systems, for example. The present study addresses a new finite element formulation of the lattice Boltzmann equation using a general weighted residual technique. Among the weighted residual formulations, the collocation method, Galerkin method, and method of moments are used for finite element based Lattice Boltzmann solutions. Different finite element geometries, such as triangular, quadrilateral, and general six-sided solids, were used in this work. Some examples using the FELBM are studied. The results were compared with both analytical and computational fluid dynamics solutions.

• Proceedings of the KIEE Conference
• /
• 1999.07a
• /
• pp.76-78
• /
• 1999

This paper presents Galerkin- and Upwind-finite element analyses solution in the travelling magnetic filed problem. The travelling magnetic field problem is subject to convective- diffusion equation. Therefore, the solution derived from Galerkin-FEM with linear interpolation function may oscillate between the adjacent nodes. A simple model with Derichlet, Noumann and periodic boundary condition respectively, have been analyzed to investigate stabilities of solutions. It is concluded that the solution of Galerkin-FEM may oscillate according to boundary condition and element type, but that of Upwind-FEM is stable regardless boundary condition.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.2
• /
• pp.145-152
• /
• 2009

SGBEM(Symmetric Galerkin Boundary Element Method)-FEM alternating method has been proposed by Nikishkov, Park and Atluri. In the proposed method, arbitrarily shaped three-dimensional crack problems can be solved by alternating between the crack solution in an infinite body and the finite element solution without a crack. In the previous study, the SGBEM-FEM alternating method was extended further in order to solve elastic-plastic crack problems and to obtain elastic-plastic stress fields. For the elastic-plastic analysis the algorithm developed by Nikishkov et al. is used after modification. In the algorithm, the initial stress method is used to obtain elastic-plastic stress and strain fields. In this paper, elastic-plastic J integrals for three-dimensional cracks are obtained using the method. For that purpose, accurate values of displacement gradients and stresses are necessary on an integration path. In order to improve the accuracy of stress near crack surfaces, coordinate transformation and partitioning of integration domain are used. The coordinate transformation produces a transformation Jacobian, which cancels the singularity of the integrand. Using the developed program, simple three-dimensional crack problems are solved and elastic and elastic-plastic J integrals are obtained. The obtained J integrals are compared with the values obtained using a handbook solution. It is noted that J integrals obtained from the alternating method are close to the values from the handbook.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.475-483
• /
• 2008

2D finite element model for analysis of transport of accidentally released pollutants in the flow was developed by SUPG method, and the mass balance of this model was checked though two example problems: line source and point source problem in the straight channel and unidirectional 2D flow field, respectively. All the test cases were simulated with both SUPG and conventional Galerkin method to compare the accuraccy of the numerical mass balance. Test results show that the model with SUPG can adequately conserve the released mass though simulation than the model using Galerkin method, so the developed model verified to be appropriate to solve this accidental mass release problem.