• Journal of applied mathematics & informatics
• /
• v.9 no.1
• /
• pp.27-43
• /
• 2002

Finite element Galerkin solutions for the strongly damped extensible beam equations are considered. The semidiscrete scheme and a fully discrete time Galerkin method are studied and the corresponding stability and error estimates are obtained. Ratios of numerical convergence are given.

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

• Bulletin of the Korean Mathematical Society
• /
• v.57 no.1
• /
• pp.219-244
• /
• 2020

In this paper, we approximate the solution of the coupled nonlinear Schrödinger equations by using a fully discrete finite element scheme based on the standard Galerkin method in space and implicit midpoint discretization in time. The proposed scheme guarantees the conservation of the total mass and the energy. First, a priori error estimates for the fully discrete Galerkin method is derived. Second, the existence of the approximated solution is proved by virtue of the Brouwer fixed point theorem. Moreover, the uniqueness of the solution is shown. Finally, convergence orders of the fully discrete Crank-Nicolson scheme are discussed. The end of the paper is devoted to some numerical experiments.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.12 no.3
• /
• pp.133-141
• /
• 2007

A finite element model is developed for the extended Boussinesq equations that is capable of simulating the dynamics of long and short waves. Galerkin weighted residual method and the introduction of auxiliary variables for 3rd spatial derivative terms in the governing equations are used for the model development. The Adams-Bashforth-Moulton Predictor Corrector scheme is used as a time integration scheme for the extended Boussinesq finite element model so that the truncation error would not produce any non-physical dispersion or dissipation. This developed model is applied to the problems of solitary wave propagation. Predicted results is compared to available analytical solutions and laboratory measurements. A good agreement is observed.

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

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

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

• Journal of the Korean Mathematical Society
• /
• v.61 no.4
• /
• pp.621-657
• /
• 2024

In this study, the numerical solution of a space-fractional Burgers' equation with initial and boundary conditions is considered. This equation is the simplest nonlinear model for diffusive waves in fluid dynamics. It occurs in a variety of physical phenomena, including viscous sound waves, waves in fluid-filled viscous elastic pipes, magneto-hydrodynamic waves in a medium with finite electrical conductivity, and one-dimensional turbulence. The proposed QBS/CNG technique consists of the Galerkin method with a function basis of quadratic B-splines for the spatial discretization of the space-fractional Burgers' equation. This is then followed by the Crank-Nicolson approach for time-stepping. A linearized scheme is fully constructed to reduce computational costs. Stability analysis, error estimates, and convergence rates are studied. Finally, some test problems are used to confirm the theoretical results and the proposed method's effectiveness, with the results displayed in tables, 2D, and 3D graphs.

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