• Korean Journal of Mathematics
• /
• v.26 no.3
• /
• pp.467-481
• /
• 2018

In this paper we present a postprocessing scheme for the Raviart-Thomas mixed finite element approximation of the second order elliptic eigenvalue problem. This scheme is carried out by solving a primal source problem on a higher order space, and thereby can improve the convergence rate of the eigenfunction and eigenvalue approximations. It is also used to compute a posteriori error estimates which are asymptotically exact for the $L^2$ errors of the eigenfunctions. Some numerical results are provided to confirm the theoretical results.

• Journal of applied mathematics & informatics
• /
• v.36 no.1_2
• /
• pp.135-154
• /
• 2018

In this paper, an almost second order overlapping Schwarz method for singularly perturbed third order convection-diffusion type problem is constructed. The method splits the original domain into two overlapping subdomains. A hybrid difference scheme is proposed in which on the boundary layer region we use the combination of classical finite difference scheme and central finite difference scheme on a uniform mesh while on the non-layer region we use the midpoint difference scheme on a uniform mesh. It is shown that the numerical approximations which converge in the maximum norm to the exact solution. We proved that, when appropriate subdomains are used, the method produces convergence of second order. Furthermore, it is shown that, two iterations are sufficient to achieve the expected accuracy. Numerical examples are presented to support the theoretical results. The main advantages of this method used with the proposed scheme are it reduce iteration counts very much and easily identifies in which iteration the Schwarz iterate terminates.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.303-309
• /
• 2007

For the propagation of elastic waves in unbounded domains, absorbing boundary conditions at the fictitious numerical boundaries have been proposed. In this paper we focus on both first and second order paraxial boundary conditions(PBCs) in the framework of variational approximations which are based on paraxial approximations of the scalar and elastic wave equations. We propose a penalty function method for the treatment of PBCs and apply these into finite element analysis. The numerical verification of the efficiency is carried out through comparing PBCs with Lysmer-Kuhlemeyer's boundary conditions.

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

A robust numerical method for a singularly perturbed second-order ordinary differential equation having two parameters with a discontinuous source term is presented in this article. Theoretical bounds are derived for the derivatives of the solution and its smooth and singular components. An appropriate piecewise uniform mesh is constructed, and classical upwind finite difference schemes are used on this mesh to obtain the discrete system of equations. Parameter-uniform error bounds for the numerical approximations are established. Numerical results are provided to illustrate the convergence of the numerical approximations.

• Journal of applied mathematics & informatics
• /
• v.30 no.1_2
• /
• pp.235-244
• /
• 2012

In this paper, we solve the three-dimensional biharmonic equation with Dirichlet boundary conditions of second kind using the full multigrid (FMG) algorithm. We derive a finite difference approximations for the biharmonic equation on a 18 point compact stencil. The unknown solution and its second derivatives are carried as unknowns at grid points. In the multigrid methods, we use a fourth order interpolation to producing a new intermediate unknown functions values on a finer grid, and the full weighting restriction operators to calculating the residuals at coarse grid points. A set of test problems gives excellent results.

• Journal of computational fluids engineering
• /
• v.11 no.1 s.32
• /
• pp.8-15
• /
• 2006

The existing approximations of diffusion flux in unstructured cell-centered finite volume methods are examined in detail with each other and clarified to have indefinite expressions in several respects. A new numerical approximation of diffusion flux at cell face center is then proposed, which is second-order accurate even on irregular grids and may be easily implemented in CFD code using cell-centered finite volume method with unstructured grids composed of arbitrary convex polyhedral shape.

• Communications for Statistical Applications and Methods
• /
• v.7 no.1
• /
• pp.233-244
• /
• 2000

Confidence intervals for survival function give useful information about the lifetime distribution. In this paper we develop Edgeworkth expansions as approximation to the true and bootstrap distributions of normalized nonparametric maximum likelihood estimator of survival function in the Koziol-Green model and then use these results to show that the bootstrap approximations have second order accuracy.

• Communications for Statistical Applications and Methods
• /
• v.26 no.6
• /
• pp.527-538
• /
• 2019

This paper considers the issue of obtaining the optimal design in Poisson regression model when the sample size is small. Poisson regression model is widely used for the analysis of count data. Asymptotic theory provides the basis for making inference on the parameters in this model. However, for small size experiments, asymptotic approximations, such as unbiasedness, may not be valid. Therefore, first, we employ the second order expansion of the bias of the maximum likelihood estimator (MLE) and derive the mean square error (MSE) of MLE to measure the quality of an estimator. We then define DM-optimality criterion, which is based on a function of the MSE. This criterion is applied to obtain locally optimal designs for small size experiments. The effect of sample size on the obtained designs are shown. We also obtain locally DM-optimal designs for some special cases of the model.

• Nuclear Engineering and Technology
• /
• v.49 no.1
• /
• pp.29-42
• /
• 2017

The purpose of the present study is the presentation of the appropriate element and shape function in the solution of the neutron diffusion equation in two-dimensional (2D) geometries. To this end, the multigroup neutron diffusion equation is solved using the Galerkin finite element method in both rectangular and hexagonal reactor cores. The spatial discretization of the equation is performed using unstructured triangular and quadrilateral finite elements. Calculations are performed using both linear and quadratic approximations of shape function in the Galerkin finite element method, based on which results are compared. Using the power iteration method, the neutron flux distributions with the corresponding eigenvalue are obtained. The results are then validated against the valid results for IAEA-2D and BIBLIS-2D benchmark problems. To investigate the dependency of the results to the type and number of the elements, and shape function order, a sensitivity analysis of the calculations to the mentioned parameters is performed. It is shown that the triangular elements and second order of the shape function in each element give the best results in comparison to the other states.

• Journal of Energy Engineering
• /
• v.16 no.1 s.49
• /
• pp.46-52
• /
• 2007

The piecewise constant angular approximation is developed to replace the conventional angular quadrature sets in the solution of the second-order, multi-dimensional $S_{N}$ neutron transport equations. The newly generated quadrature sets by this method substantially mitigate ray effects and can be used in the same manner as the conventional quadrature sets are used. The discrete-ordinates and the piecewise-constant approximations are applied to both the first-order Boltzmann and the second-order form of neutron transport equations in treating angular variables. The result is that the mitigation of ray effects is only achieved by the piecewise-constant method, in which new angular quadratures are generated by integrating angle variables over the specified region. In other sense, the newly generated angular quadratures turn out to decrease the contribution of mixed-derivative terms in the even-parity equation that is one of the second-order neutron transport equation. This result can be interpreted as the entire elimination or substantial mitigation of ray effect are possible in the simplified even-parity equation which has no mixed-derivative terms.