• Journal of Institute of Control, Robotics and Systems
• /
• v.7 no.6
• /
• pp.465-470
• /
• 2001

This paper considers a new weighed model reduction method using block diagonal solutions of Lyapunov inequalities. With the input and/or output weighting function, the stability of the reduced order system is guaranteed and an a priori error bound is proposed. to achieve this after finding the solutions of two Lyapunov inequalities and balancing the full order system, we find the reduced order systems using the direct truncation and the singular perturbation approximation. The proposed method is compared with other existing methods using numerical examples.

• Journal of the Korean Mathematical Society
• /
• v.58 no.3
• /
• pp.553-569
• /
• 2021

In this paper, we present and analyze a fully discrete numerical method for solving the time-fractional diffusion wave equation: ∂^β_tu - div(a∇u) = f, 1 < β < 2. We first construct a difference formula to approximate ∂^β_tu by using an interpolation of derivative type. The truncation error of this formula is of O(△t^2+δ-β)-order if function u(t) ∈ C^2,δ[0, T] where 0 ≤ δ ≤ 1 is the Hölder continuity index. This error order can come up to O(△t^3-β) if u(t) ∈ C³ [0, T]. Then, in combinination with the linear finite volume discretization on spatial domain, we give a fully discrete scheme for the fractional wave equation. We prove that the fully discrete scheme is unconditionally stable and the discrete solution admits the optimal error estimates in the H¹-norm and L₂-norm, respectively. Numerical examples are provided to verify the effectiveness of the proposed numerical method.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.23 no.9A
• /
• pp.2365-2371
• /
• 1998

Discrete cosine transform (DCT) has wide applications in speech and image coding. In this paper, we propose a novel fast dCT scheme with the property of reduced multiplication stages and the smaller number of additions and multiplications. This exploits the symmetry property of the DCT kernel to decompose the N-point dCT to N/2 point, and can be generally applied recursively to $2^{m}$-point. The proposed algorithm has a structure that most of multiplications tend to be performed at final stage, and this reduces propagation of truncation error which could occur in the fixed-point computation. Also the minimization of the multiplication stages further decreases the error.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.38 no.9
• /
• pp.933-941
• /
• 2014

This paper provides a comparison between the Krylov subspace method (KSM) and modal truncation method (MTM), which are typical projection-based model order reduction methods. The frequency responses are compared to determine the numerical accuracies and efficiencies. In order to compare the numerical accuracies of the KSM and MTM, the frequency responses and relative errors according to the order of the reduced model and frequency of interest are studied. Subsequently, a numerical examination shows whether a reduced order can be determined automatically with the help of an error convergence indicator. As for the numerical efficiency, the computation time needed to generate the projection matrix and the solution time to perform a frequency response analysis are compared according to the reduced order. A finite element model for a car suspension is considered as an application example of the numerical comparison.

• Nuclear Engineering and Technology
• /
• v.54 no.1
• /
• pp.110-116
• /
• 2022

Human failure event (HFE) dependency analysis is a part of human reliability analysis (HRA). For efficient HFE dependency analysis, a maximum number of minimal cut sets (MCSs) that have HFE combinations are generated from the fault trees for the probabilistic safety assessment (PSA) of nuclear power plants (NPPs). After collecting potential HFE combinations, dependency levels of subsequent HFEs on the preceding HFEs in each MCS are analyzed and assigned as conditional probabilities. Then, HFE recovery is performed to reflect these conditional probabilities in MCSs by modifying MCSs. Inappropriate HFE dependency analysis and HFE recovery might lead to an inaccurate core damage frequency (CDF). Using the above process, HFE recovery is performed on MCSs that are generated with a non-zero truncation limit, where many MCSs that have HFE combinations are truncated. As a result, the resultant CDF might be underestimated. In this paper, a new method is suggested to incorporate HFE recovery into the MCS generation stage. Compared to the current approach with a separate HFE recovery after MCS generation, this new method can (1) reduce the total time and burden for MCS generation and HFE recovery, (2) prevent the truncation of MCSs that have dependent HFEs, and (3) avoid CDF underestimation. This new method is a simple but very effective means of performing MCS generation and HFE recovery simultaneously and improving CDF accuracy. The effectiveness and strength of the new method are clearly demonstrated and discussed with fault trees and HFE combinations that have joint probabilities.

• Journal of Computing Science and Engineering
• /
• v.3 no.2
• /
• pp.73-87
• /
• 2009

A hybrid system is a dynamical system in which states can be changed continuously and discretely. Simulation based on numerical methods is the widely used technique for analyzing complicated hybrid systems. Numerical simulation of hybrid systems, however, is subject to two types of numerical errors: truncation error and round-off error. The effect of such errors can make an impossible transition step to become possible during simulation, and thus, to generate a simulation behavior that is not allowed by the model. The possibility of an incorrect simulation behavior reduces con.dence in simulation-based analysis since it is impossible to know whether a particular simulation trace is allowed by the model or not. To address this problem, we define the notion of Instrumented Hybrid Automata (IHA), which considers the effect of accumulated numerical errors on discrete transition steps. We then show how to convert Hybrid Automata (HA) to IRA and prove that every simulation behavior of IHA preserves the discrete transition steps of some behavior in HA; that is, simulation of IHA is sound with respect to HA.

• Journal of applied mathematics & informatics
• /
• v.15 no.1_2
• /
• pp.29-51
• /
• 2004

In this paper, we propose a new mixed finite element method, called the characteristics-mixed method, for approximating the solution to Burgers' equation. This method is based upon a space-time variational form of Burgers' equation. The hyperbolic part of the equation is approximated along the characteristics in time and the diffusion part is approximated by a mixed finite element method of lowest order. The scheme is locally conservative since fluid is transported along the approximate characteristics on the discrete level and the test function can be piecewise constant. Our analysis show the new method approximate the scalar unknown and the vector flux optimally and simultaneously. We also show this scheme has much smaller time-truncation errors than those of standard methods. Numerical example is presented to show that the new scheme is easily implemented, shocks and boundary layers are handled with almost no oscillations. One of the contributions of the paper is to show how the optimal error estimates in $L^2(\Omega)$ are obtained which are much more difficult than in the standard finite element methods. These results seem to be new in the literature of finite element methods.

• The Journal of the Korea Contents Association
• /
• v.6 no.2
• /
• pp.162-168
• /
• 2006

This paper presents a high performance CORDIC circuit based on redundant numbers yielding a minimal number of iteration stages. The minimal number of iteration stages reflects the iteration number yielding a smaller computation error than the truncation error. The minimal number of iterations is found n-4 for $n\geq16$, where n is the number of input angle bits. The CORDIC circuit is based on a redundant number system with a constant scale factor The circuit performs sine and cosine calculations with a delay of {5 (n-4)+ 2[$log_{2}n$]}$\DeltaT$.

• Earthquakes and Structures
• /
• v.14 no.1
• /
• pp.45-53
• /
• 2018

A family of the structure-dependent methods seems very promising for time integration since it can simultaneously have desired numerical properties, such as unconditional stability, second-order accuracy, explicit formulation and numerical dissipation. However, an unusual overshoot, which is essentially different from that found by Goudreau and Taylor in the transient response, has been experienced in the steady-state response of a high frequency mode. The root cause of this unusual overshoot is analytically explored and then a remedy is successfully developed to eliminate it. As a result, an improved formulation of this family method can be achieved.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.11a
• /
• pp.39-44
• /
• 1996

A nodal transport method based on the F$_{N}$ method is developed for the transport calculation in x- y geometry and tested for benchmark problems. Using transverse integration, the two-dimensional transport equation is converted to one-dimensional equations for x, y-directions and the one-dimensional equations are integrated over azimuthal angle. With proper approximations for the transverse leakage, the one-dimensional equations are discretized by using the F$_{N}$ method without truncation error. At present, isotropic approximation of the transverse with a quadratic or flat shape in spatial variable is tested.ted.