• Journal of applied mathematics & informatics
• /
• v.29 no.3_4
• /
• pp.813-829
• /
• 2011

In this article, we present a numerical scheme for solving singularly perturbed (i.e. highest -order derivative term multiplied by small parameter) Burgers-Huxley equation with appropriate initial and boundary conditions. Most of the traditional methods fail to capture the effect of layer behavior when small parameter tends to zero. The presence of perturbation parameter and nonlinearity in the problem leads to severe difficulties in the solution approximation. To overcome such difficulties the present numerical scheme is constructed. In construction of the numerical scheme, the first step is the dicretization of the time variable using forward difference formula with constant step length. Then, the resulting non linear singularly perturbed semidiscrete problem is linearized using quasi-linearization process. Finally, differential quadrature method is used for space discretization. The error estimate and convergence of the numerical scheme is discussed. A set of numerical experiment is carried out in support of the developed scheme.

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

• Proceedings of the KIEE Conference
• /
• 1990.11a
• /
• pp.413-416
• /
• 1990

This paper proposes a design procedure based on the LQG/LTR method for a launch vehicle autopilot. Continuous-discrete type LQG/LTR compensators are designed using the $\delta$-transformation [1] in order to overcome numerical problems occurring in the process of discretization. The $\delta$-LQG/LTR compensator using the $\delta$-transformation is compared with the $\delta$-LQG/LTR compensator using the $\delta$-transformation. The performance of the overall system controlled by the $\delta$-LQG/LTR compensator is evaluated via simulations, which show that the discretization error problem is resolved and the control performances are satisfied in the proposed compensator.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1994.04a
• /
• pp.1-8
• /
• 1994

An automated procedure which allows adaptation of spatial and time discretization simultaneously in finite element analysis of linear elastodynamic problems is presented. For dynamic problems having responses dominated by high frequency modes, such as those with impact, explosive, traveling and earthquake loads high gradient stress regions change their locations from time to time. And the time step size may need to vary in order to deal with whole process ranging from transient phase to steady state phase. As the sizes of elements in space vary in different regions, the procedure also permits different time stepping. In such a way, the best performance attainable by the finite element method can be achieved. In this study, we estimate both of the kinetic energy error and stran energy error induced by spatial and time discretization in a consistent manner. Numerical examples are used to demonstrate the performance of the procedure.

• Journal of applied mathematics & informatics
• /
• v.29 no.5_6
• /
• pp.1129-1141
• /
• 2011

In this paper, a class of large time-stepping method based on the finite element discretization for the Cahn-Hilliard equation with the Neumann boundary conditions is developed. The equation is discretized by finite element method in space and semi-implicit schemes in time. For the first order fully discrete scheme, convergence property is investigated by using finite element analysis. Numerical experiment is presented, which demonstrates the effectiveness of the large time-stepping approaches.

• Journal of computational fluids engineering
• /
• v.8 no.3
• /
• pp.50-55
• /
• 2003

Recent numerical simulation has a tendency to require the higher-order accuracy in time, as well as in space. This tendency is more true in LES and acoustic noise simulation. In the present work, the accuracy of a Fractional step method, which is widely used in LES simulation, has been increased to the fourth-order accurate compact Pade discretization. To validate the present code, the flow-field past a cylinder was simulated and compared with experiment. A good agreement with experiment was achieved.

• Journal of KSNVE
• /
• v.3 no.4
• /
• pp.353-359
• /
• 1993

In this study, the Taylor-Galerkin method is modified to take into consideration the third order term in the Taylor series of the fundamental variable. In the Taylor-Galerkin method, after expressing the governing equation of motion in conservation form, the temporal discretization is done first and then spatial discretization follows in contrast to the conventional approaches. A predictor-corrector type algorithm has been developed previously by the same author. A new computationally efficient direct algorithm is proposed in this study. A study on convergency and accuracy of the solution is carried out. Numerical examples show that this new algorithm exhibits the same order of accuracy with less computational effort.

• Structural Engineering and Mechanics
• /
• v.75 no.6
• /
• pp.685-699
• /
• 2020

Polygonal finite element provides a great flexibility in mesh generation of crack propagation problems where the topology of the domain changes significantly. However, the control of the discretization error in such problems is a main concern. In this paper, a polygonal-FEM is presented in modeling of crack propagation problems via an automatic adaptive mesh refinement procedure. The adaptive mesh refinement is accomplished based on the Zienkiewicz-Zhu error estimator in conjunction with a weighted SPR technique. Adaptive mesh refinement is employed in some steps for reduction of the discretization error and not for tracking the crack. In the steps that no adaptive mesh refinement is required, local modifications are applied on the mesh to prevent poor polygonal element shapes. Finally, several numerical examples are analyzed to demonstrate the efficiency, accuracy and robustness of the proposed computational algorithm in crack propagation problems.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2558-2561
• /
• 2008

In the present study, a parallel Taylor-Galerkin/level set based two-phase flow code was developed using finite element discretization and domain decomposition method based on MPI (Message Passing Interface). The proposed method can be utilized for the analysis of a large scale free surface problem in a complex geometry due to the feature of FEM and domain decomposition method. Four-step fractional step method was used for the solution of the incompressible Navier-Stokes equations and Taylor-Galerkin method was adopted for the discretization of hyperbolic type redistancing and advection equations. A Parallel ILU(0) type preconditioner was chosen to accelerate the convergence of a conjugate gradient type iterative solvers. From the present parallel numerical experiments, it has been shown that the proposed method is applicable to the simulation of large scale free surface flows.

• Nuclear Engineering and Technology
• /
• v.53 no.12
• /
• pp.3861-3878
• /
• 2021

The steady-state simplified spherical harmonics equations (SP_N equations) are a higher order approximation to the neutron transport equations than the neutron diffusion equation that also have reasonable computational demands. This work extends these results for the analysis of transients by comparing of two formulations of time-dependent SP_N equations considering different treatments for the time derivatives of the field moments. The first is the full system of equations and the second is a diffusive approximation of these equations that neglects the time derivatives of the odd moments. The spatial discretization of these methodologies is made by using a high order finite element method. For the time discretization, a semi-implicit Euler method is used. Numerical results show that the diffusive formulation for the time-dependent simplified spherical harmonics equations does not present a relevant loss of accuracy while being more computationally efficient than the full system.