• Journal of applied mathematics & informatics
• /
• v.33 no.5_6
• /
• pp.657-670
• /
• 2015

In this paper, we provide numerical analysis of so-called Legendre Gauss-Radau and Legendre-Gauss collocation methods for ordinary differential equations. After recasting these collocation methods as Runge-Kutta methods, we prove that the Legendre-Gauss collocation method is equivalent to the well-known Gauss method, while the Legendre-Gauss-Radau collocation method does not belong to the classes of Radau IA or Radau IIA methods in the Runge-Kutta literature. Making use of the well-established theory of Runge-Kutta methods, we study stability and accuracy of the Legendre-Gauss-Radau collocation method. Numerical experiments are conducted to confirm our theoretical results on the accuracy and numerical stability of the Legendre-Gauss-Radau collocation method, and compare Legendre-Gauss collocation method with the Gauss method.

• The Pure and Applied Mathematics
• /
• v.30 no.4
• /
• pp.347-355
• /
• 2023

For numerical evaluation of Cauchy principal value integrals, we present a simple rational function with a parameter satisfying some reasonable conditions. The proposed rational function is employed in coordinate transformation for accelerating the accuracy of the Gauss quadrature rule. The efficiency of the proposed rational transformation method is demonstrated by the numerical result of a selected test example.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.2
• /
• pp.264-277
• /
• 2015

This article investigates various transcription techniques for the Legendre pseudospectral (PS) method to compare the pros and cons of each approach. Eight combinations from four different types of collocation points and two discretization methods for dynamic constraints, which differentiate Legendre PS transcription techniques, are implemented to solve a carefully selected test set of nonlinear optimal control problems (NOCPs). The convergence property and prediction accuracy are compared to provide a useful guideline for selecting the best combination. The tested NOCPs consist of the minimum time, minimum energy, and problems with state and control constraints. Therefore, the results drawn from this comparative study apply to the solution of similar types of NOCPs and can mitigate much debate about the best combinations. Additionally, important findings from this study can be used to improve the numerical efficiency of the Legendre PS methods. Three PS applications to the aerospace engineering problems are demonstrated to prove this point.

• Journal of the Korean Data and Information Science Society
• /
• v.11 no.2
• /
• pp.139-155
• /
• 2000

An exclusive simulation study is conducted in computing means for order statistics in standard normal variate. Monte Carlo moments are used in Shapiro-Francia W' statistic computation. Finally, quantiles for Shapiro-Francia W' are generated. The study shows that in computing means for order statistics in standard normal variate, complicated distributions and intensive numerical integrations can be avoided by using Monte Carlo simulation. Lack of accuracy is minimal and computation simplicity is noteworthy.

• The KIPS Transactions:PartD
• /
• v.10D no.7
• /
• pp.1145-1148
• /
• 2003

We show the mixed finite element method which induces solutions that has the same order of errors for both the gradient of the solution and the solution itself. The technique to construct an efficient reusable matrix is suggested. Two families of mixed finite element methods are introduced with an automatic generating technique for matrix with my order of basis. The generated matrix by this technique has more accurate values and is a sparse matrix. This new technique is applied to solve a minimal surface problem.

• Journal of the Korean Mathematical Society
• /
• v.41 no.6
• /
• pp.957-976
• /
• 2004

It is well known that the application of the nonlinear coordinate transformations is useful for efficient numerical evaluation of weakly singular integrals. In this paper, we consider the trapezoidal rule combined with a nonlinear transformation $\Omega$$_{m}$(b;$\chi$), containing a parameter b, proposed first by Yun [14]. It is shown that the trapezoidal rule with the transformation $\Omega$$_{m}$(b;$\chi$), like the case of the Gauss-Legendre quadrature rule, can improve the asymptotic truncation error by using a moderately large b. By several examples, we compare the numerical results of the present method with those of some existing methods. This shows the superiority of the transformation $\Omega$$_{m}$(b;$\chi$).TEX>).

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.2
• /
• pp.87-95
• /
• 2018

This paper introduces a new explicit spectral element method for the simulation of SH-waves in semi-infinite domains. To simulate the wave motion in unbounded domains, it is necessary to reduce the infinite extent to a finite computational domain of interest. To prevent the wave reflection from the trunctated boundaries, perfectly matched layer(PML) wave-absorbing boundary is introduced. The forward problem for simulating SH-waves in PML-truncated domains can be formulated as second-order PDEs. The second-order semi-discrete form of the governing PDEs is constructed by using a mixed spectral elements with Legendre-gauss-Lobatto quadrature method, which results in a diagonalized mass matrix. Then the second-order semi-discrete form is transformed to a first-order, whose solutions are calculated by the fourth-order Runge-Kutta method. Numerical examples showed that solutions of SH-wave in the two-dimensional analysis domain resulted in stable and accurate, and reflections from truncated boundaries could be reduced by using PML boundaries. Elastic wave propagation analysis using explicit time integration method may be apt for solving larger domain problems such as three-dimensional elastic wave problem more efficiently.

• Structural Engineering and Mechanics
• /
• v.68 no.5
• /
• pp.603-619
• /
• 2018

In this research, an effective computational technique is carried out for nonlinear and post-buckling analyses of cracked imperfect composite plates. The laminated plates are assumed to be moderately thick so that the analysis can be carried out based on the first-order shear deformation theory. Geometric non-linearity is introduced in the way of von-Karman assumptions for the strain-displacement equations. The Ritz technique is applied using Legendre polynomials for the primary variable approximations. The crack is modeled by partitioning the entire domain of the plates into several sub-plates and therefore the plate decomposition technique is implemented in this research. The penalty technique is used for imposing the interface continuity between the sub-plates. Different out-of-plane essential boundary conditions such as clamp, simply support or free conditions will be assumed in this research by defining the relevant displacement functions. For in-plane boundary conditions, lateral expansions of the unloaded edges are completely free while the loaded edges are assumed to move straight but restricted to move laterally. With the formulation presented here, the plates can be subjected to biaxial compressive loads, therefore a sensitivity analysis is performed with respect to the applied load direction, along the parallel or perpendicular to the crack axis. The integrals of potential energy are numerically computed using Gauss-Lobatto quadrature formulas to get adequate accuracy. Then, the obtained non-linear system of equations is solved by the Newton-Raphson method. Finally, the results are presented to show the influence of crack length, various locations of crack, load direction, boundary conditions and different values of initial imperfection on nonlinear and post-buckling behavior of laminates.

• Journal of Radiation Protection and Research
• /
• v.17 no.1
• /
• pp.43-56
• /
• 1992

As one of the methods to ameliorate the ray effects which are the nature of anomalous computational effects due to the discretization of the angular variable in discrete ordinates approximations, a computational program, named TWODET (TWO dimensional Discrete Element Transport), has developed in 2 dimensional cartesian coordinates system using the discrete elements method, in which the discrete angle quadratures are steered by the spatially dependent angular fluxes. The results of the TWODET calculation with K-2, L-3 discrete angular quadratures, in the problem of a centrally located, isotropically emitting flat source in an absorbing square, are shown to be more accurate than that of the DOT 4.3 calculation with S-10 full symmetry angular quadratures, in remedy of the ray effect at the edge flux distributions of the square. But the computing time of the TWODET is about 4 times more than that of the DOT 4.3. In the problem of vacuum boundaries just outside of the source region in an absorbing square, the results of the TWODET calculation are shown severely anomalous ray effects, due to the sudden discontinuity between the source and the vacuum, like as the results of the DOT 4.3 calculation. In the probelm of an external source in an absorbing square in which a highly absorbing medium is added, the results of the TWODET calculation with K-3, L-4 show a good ones like as, somewhat more than, that of the DOT 4.3 calculation with S-10.