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

The existing numerical approximations of convection flux, especially the spatial higher-order difference schemes, in unstructured cell-centered finite volume methods are examined in detail with each other and evaluated with respect to the accuracy through their application to a 2-D benchmark problem. Six higher-order schemes are examined, which include two second-order upwind schemes, two central difference schemes and two hybrid schemes. It is found that the 2nd-order upwind scheme by Mathur and Murthy(1997) and the central difference scheme by Demirdzic and Muzaferija(1995) have more accurate prediction performance than the other higher-order schemes used in unstructured cell-centered finite volume methods.

• Journal of the Korean Mathematical Society
• /
• v.49 no.3
• /
• pp.549-569
• /
• 2012

In this paper asymptotic error expansions for mixed finite element approximations to a class of second order elliptic optimal control problems are derived under rectangular meshes, and the Richardson extrapolation of two different schemes and interpolation defect correction can be applied to increase the accuracy of the approximations. As a by-product, we illustrate that all the approximations of higher accuracy can be used to form a class of a posteriori error estimators of the mixed finite element method for optimal control problems.

• Communications for Statistical Applications and Methods
• /
• v.29 no.4
• /
• pp.471-486
• /
• 2022

Autoregressive moving average (ARMA) models involve nonlinearity in the model coefficients because of unobserved lagged errors, which complicates the likelihood function and makes the posterior density analytically intractable. In order to overcome this problem of posterior analysis, some approximation methods have been proposed in literature. In this paper we first review the main analytic approximations proposed to approximate the posterior density of ARMA models to be analytically tractable, which include Newbold, Zellner-Reynolds, and Broemeling-Shaarawy approximations. We then use the Kullback-Leibler divergence to study the relation between these three analytic approximations and to measure the distance between their derived approximate posteriors for ARMA models. In addition, we evaluate the impact of the approximate posteriors distance in Bayesian estimates of mean and precision of the model coefficients by generating a large number of Monte Carlo simulations from the approximate posteriors. Simulation study results show that the approximate posteriors of Newbold and Zellner-Reynolds are very close to each other, and their estimates have higher precision compared to those of Broemeling-Shaarawy approximation. Same results are obtained from the application to real-world time series datasets.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.18 no.4
• /
• pp.337-350
• /
• 2014

In this paper, we consider discontinuous Galerkin finite element methods with interior penalty term to approximate the solution of nonlinear parabolic problems with mixed boundary conditions. We construct the finite element spaces of the piecewise polynomials on which we define fully discrete discontinuous Galerkin approximations using the Crank-Nicolson method. To analyze the error estimates, we construct an appropriate projection which allows us to obtain the optimal order of a priori ${\ell}^{\infty}(L^2)$ error estimates of discontinuous Galerkin approximations in both spatial and temporal directions.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.4
• /
• pp.301-327
• /
• 2019

The proper orthogonal decomposition (POD) method for time-dependent problems significantly reduces the computational time as it reduces the original problem to the lower dimensional space. Even a higher degree of reduction can be reached if the solution is smooth in space and time. However, if the solution is discontinuous and the discontinuity is parameterized e.g. with time, the POD approximations are not accurate in the reduced space due to the lack of ability to represent the discontinuous solution as a finite linear combination of smooth bases. In this paper, we propose to post-process the sample solutions and re-initialize the POD approximations to deal with discontinuous solutions and provide accurate approximations while the computational time is reduced. For the post-processing, we use the Gegenbauer reconstruction method. Then we regularize the Gegenbauer reconstruction for the construction of POD bases. With the constructed POD bases, we solve the given PDE in the reduced space. For the POD approximation, we re-initialize the POD solution so that the post-processed sample solution is used as the initial condition at each sampling time. As a proof-of-concept, we solve both one-dimensional linear and nonlinear hyperbolic problems. The numerical results show that the proposed method is efficient and accurate.

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

• Structural Engineering and Mechanics
• /
• v.19 no.6
• /
• pp.749-755
• /
• 2005

In this work, thermo-elastic stability behavior of laminated cross-ply elliptical cylindrical shells subjected to uniform temperature rise is studied employing the finite element approach based on higher-order theory that accounts for the transverse shear and transverse normal deformations, and nonlinear in-plane displacement approximations through the thickness with slope discontinuity at the layer interfaces. The combined influence of higher-order shear deformation, shell geometry and non-circularity on the prebuckling thermal stress distribution and critical temperature parameter of laminated elliptical cylindrical shells is examined.

• Bulletin of the Korean Mathematical Society
• /
• v.52 no.2
• /
• pp.593-601
• /
• 2015

Let $f{\in}C^r$ ([-1,1]), $r{\geq}0$ and let $L^*$ be a linear left fractional differential operator such that $L^*$ $(f){\geq}0$ throughout [0, 1]. We can find a sequence of polynomials $Q_n$ of degree ${\leq}n$ such that $L^*$ $(Q_n){\geq}0$ over [0, 1], furthermore f is approximated left fractionally and simulta-neously by $Q_n$ on [-1, 1]. The degree of these restricted approximations is given via inequalities using a higher order modulus of smoothness for $f^{(r)}$.

• Journal of the Korean Mathematical Society
• /
• v.37 no.6
• /
• pp.1031-1042
• /
• 2000

The equations prescribed in Ω⊂R(sup)n are called loaded, if they contain some operations of the traces of desired solution on manifolds (of dimension which is strongly less than n) from closure Ω. These equations result from approximations of nonlinear equations by linear ones, in the problems of optimal control when the control when the control actions depends on a part of independent variables, in investigations of the inverse problems and so on. In present work we study the nonlocal boundary value problems for first-order loaded differential operator equations. Criterion of unique solvability is established. We illustrate the obtained results by examples.

• Steel and Composite Structures
• /
• v.8 no.5
• /
• pp.343-359
• /
• 2008

In the present study, an efficient method for the optimum design of three-dimensional (3D) steel framed structures is proposed. In this method, in addition to choosing the best position of columns based on architectural requirements, the optimum cross-sectional dimensions of elements are determined. The preliminary design variables are considered as the number of columns in structural plan, which are determined by a direct optimization method suitable for discrete variables, without requiring the evaluation of derivatives. After forming the geometry of structure, the main variables of the cross-sectional dimensions are evaluated, which satisfy the design constraints and also achieve the least-weight of the structure. To reduce the number of finite element analyses and the overall computational time, a new third order approximate function is introduced which employs only the diagonal elements of the higher order derivatives matrices. This function produces a high quality approximation and also, a robust optimization process. The main feature of the proposed techniques that the higher order derivatives are established by the first order exact derivatives. Several examples are solved and efficiency of the new approximation method and also, the proposed method for the best position of columns in 3D steel framed structures is discussed.