• Nuclear Engineering and Technology
• /
• v.31 no.2
• /
• pp.172-181
• /
• 1999

Series expansion methods to compute the exponential of a matrix have been compared by applying them to fuel depletion calculations. Specifically, Taylor, Pade, Chebyshev, and rational Chebyshev approximations have been investigated by approximating the exponentials of bum matrices by truncated series of each method with the scaling and squaring algorithm. The accuracy and efficiency of these methods have been tested by performing various numerical tests using one thermal reactor and two fast reactor depletion problems. The results indicate that all the four series methods are accurate enough to be used for fuel depletion calculations although the rational Chebyshev approximation is relatively less accurate. They also show that the rational approximations are more efficient than the polynomial approximations. Considering the computational accuracy and efficiency, the Pade approximation appears to be better than the other methods. Its accuracy is better than the rational Chebyshev approximation, while being comparable to the polynomial approximations. On the other hand, its efficiency is better than the polynomial approximations and is similar to the rational Chebyshev approximation. In particular, for fast reactor depletion calculations, it is faster than the polynomial approximations by a factor of ∼ 1.7.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.12
• /
• pp.2631-2635
• /
• 2002

The semi-analytic sensitivity analysis using Pade approximation is presented for linear elastic structures. Although the semi-analytic method has several advantages, accuracy of the method prevents it from practical application. One of promising remedies is the use of geometric series for the matrix inversion. Though series expansion of order three has been successfully applied to the calculation of the structural sensitivity in the most range of the design perturbation, it is prone to have a slow convergence for large perturbation. To overcome this shortage, Pade approximation is introduced so that it can broaden the trust region of the perturbation without adding expansion terms. Numerical results show that the confident sensitivity can be obtained with tiny expenses of computation effort.

• Interaction and multiscale mechanics
• /
• v.2 no.3
• /
• pp.295-319
• /
• 2009

In this work, we discuss a reproducing kernel collocation method (RKCM) for solving $2^{nd}$ order PDE based on strong formulation, where the reproducing kernel shape functions with compact support are used as approximation functions. The method based on strong form collocation avoids the domain integration, and leads to well-conditioned discrete system of equations. We investigate the convergence and the computational complexity for this proposed method. An important result obtained from the analysis is that the degree of basis in the reproducing kernel approximation has to be greater than one for the method to converge. Some numerical experiments are provided to validate the error analysis. The complexity of RKCM is also analyzed, and the complexity comparison with the weak formulation using reproducing kernel approximation is presented.

• Journal of Biomedical Engineering Research
• /
• v.41 no.1
• /
• pp.28-34
• /
• 2020

In this paper, we introduce motion artifact reduction algorithm for interleaved MRI using an advanced 3D approximation algorithm. The motion artifact framework of this paper is data corrected by post-processing with a new 3-D approximation algorithm which uses data structure for each voxel. In this study, we simulate and evaluate our algorithm using Shepp-Logan phantom and T1-MRI template for both scattered dataset and uniform dataset. We generated motion artifact using random generated motion parameters for the interleaved MRI. In simulation, we use image coregistration by SPM12 (https://www.fil.ion.ucl.ac.uk/spm/) to estimate the motion parameters. The motion artifact correction is done with using full dataset with estimated motion parameters, as well as use only one half of the full data which is the case when the half volume is corrupted by severe movement. We evaluate using numerical metrics and visualize error images.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.6
• /
• pp.553-563
• /
• 2008

If an arbitrary topography is approximated to a number of vertical steps, both variational approximation and eigenfunction expansion method can be used to compute linear wave transformation over the bottom. In this study a scatterer method associated with variational approximation is proposed to calculate reflection and transmission coefficients. Present method may be shown to be more simple and direct than the successive-application-matrix method by O'Hare and Davies. And Several numerical examples are given which are in good agreement with existing results.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.27 no.3
• /
• pp.180-193
• /
• 2023

Approximating the implied volatilities and estimating the model parameters are important topics in quantitative finance. This study proposes an approximation formula for short-maturity near-the-money implied volatilities in stochastic volatility models. A general second-order nonlinear PDE for implied volatility is derived in terms of time-to-maturity and log-moneyness from the Feyman-Kac formula. Using regularity conditions and the Taylor expansion, an approximation formula for implied volatility is obtained for short-maturity nearthe-money call options in two stochastic volatility models: Heston model and SABR model. In addition, we proposed a novel numerical method to estimate model parameters. This method reduces the number of model parameters that should be estimated. Generating sample data on log-moneyness, time-to-maturity, and implied volatility, we estimate the model parameters fitting the sample data in the above two models. Our method provides parameter estimates that are close to true values.

• The Pure and Applied Mathematics
• /
• v.30 no.1
• /
• pp.67-81
• /
• 2023

The current article manages with new generalization of Post-Widder operators preserving constant function and other test functions in Bohmann-Korovkin sense and studies the approximation properties via different estimation tools like modulus of continuity and approximation in weighted spaces. The viability of the recently modified operators as per classical Post-Widder operators is introduced in specific faculties also. Numerical examples are additionally introduced to verify our theortical results. In second last section we introduce Grüss-Voronovskaya results and in last section, we show the better approximation our new modified operators via graphical exmaples using Mathematica.

• 제어로봇시스템학회:학술대회논문집
• /
• 1998.10a
• /
• pp.440-444
• /
• 1998

In this paper, we consider numerical solution of a H-J-B (Hamilton-Jacobi-Bellman) equation of elliptic type arising from the stochastic control problem. For the numerical solution of the equation, we take an approach involving contraction mapping and finite difference approximation. We choose the It(equation omitted) type stochastic differential equation as the dynamic system concerned. The numerical method of solution is validated computationally by using the constructed test case. Map of optimal controls is obtained through the numerical solution process of the equation. We also show how the method applies by taking a simple example of nonlinear spacecraft control.

• Kyungpook Mathematical Journal
• /
• v.61 no.2
• /
• pp.383-393
• /
• 2021

In this paper, numerical techniques are presented for solving initial value problems of fractional differential equations with variable coefficients. The method is derived by applying a Taylor vector approximation. Moreover, the operational matrix of fractional integration of a Taylor vector is provided in order to transform the continuous equations into a system of algebraic equations. Furthermore, numerical examples demonstrate that this method is applicable and accurate.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.5
• /
• pp.446-456
• /
• 2007

Numerical analyses of stem waves, the interaction between incident and reflected waves of obliquely incident regular waves along a vertical wall in a constant water depth, are presented. For the numerical model of the analysis, the two-layer Boussinesq equations developed by Lynett and Liu(2004a,b) are employed. Numerical results are compared with both laboratory measurements and those obtained using parabolic approximation model. The overall comparisons between the results from the two numerical models and the experiments are good. However, the two-layer Boussinesq model is more accurate than the parabolic approximation model as the angle of incident waves increases. In particular, the higher harmonic generation due to the wave nonlinearity is captured only in the Boussinesq model.