• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.3
• /
• pp.341-349
• /
• 1997

The characteristics of defect correction method are discussed in a sample heat conduction problem showing the numerical solution of the error correction equation can predict the error of the numerical solution of the original governing equation. A way of using defect correction method combined with the existing algorithm for the incompressible fluid flow, is proposed and subsequently tested for the driven square cavity problem. The error correction equations for the continuity equation and the momentum equations are considered to estimate the errors of the numerical solutions of the original governing equations. With this new approach, better velocity and pressure fields can be obtained by correcting the original numerical solutions using the estimated errors. These calculated errors also can be used to estimate the orders of magnitude of the errors of the original numerical solutions.

• Journal of KSNVE
• /
• v.9 no.4
• /
• pp.687-692
• /
• 1999

Digital holographic reconstruction method [W. A. Veronesi, et al, 1989, JASA, Vol. 85, pp. 588∼598] reconstructs a whole sound field by using measured sound pressure. This paper reports numerical errors of the method that occur at specific frequencies. The errors occur due to the truncation errors included in the calculation of transfer matrixes. The frequencies of the errors depend on the size of boundary element surfaces. Moreover, a modified calculation technique is proposed in the paper. The technique prevents the truncation errors by employing an indirect calculation procedure.

• Atmosphere
• /
• v.30 no.2
• /
• pp.115-124
• /
• 2020

The Northern Hemisphere extratropical prediction skill of the Korea Meteorological Administration (KMA) Global Data Assimilation and Prediction System (GDAPS) is examined for January 2019. The real-time prediction skill, evaluated with mean squared skill score (MSSS) of 30-90°N geopotential height field at 500 hPa (Z500), is ~8 days in the troposphere. The MSSS of Z500 considerably decreases after 3 days mainly due to the increasing eddy errors. The eddy errors are largely explained by the eddy-phased errors with minor contribution of amplitude errors. In particular, planetary-scale eddy errors are considered as a main reason of rapidly increasing errors. It turns out that such errors are associated with the blocking highs over North Pacific (NP) and Euro-Atlantic (EA) regions. The model overestimates the blocking highs over NP and EA regions in time, showing dependence of blocking predictability on blocking initializations. This result suggests that the extratropical prediction skill could be improved by better representing blocking in the model.

• 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 Theoretical and Applied Mechanics
• /
• v.3 no.1
• /
• pp.16-33
• /
• 2002

A concept of hierarchical modeling, the newest modeling technology. has been introduced early In 1990. This nu technology has a goat potential to advance the capabilities of current computational mechanics. A first step to Implement this concept is to construct hierarchical models, a family of mathematical models which are sequentially connected by a key parameter of the problem under consideration and have different levels in modeling accuracy, and to investigate characteristics In their numerical simulation aspects. Among representative model problems to explore this concept are elastic structures such as beam-, arch-. plate- and shell-like structures because the mechanical behavior through the thickness can be approximated with sequential accuracy by varying the order of thickness polynomials in the displacement or stress fields. But, in the numerical analysis of hierarchical models, two kinds of errors prevail: the modeling error and the numerical approximation errors. To ensure numerical simulation quality, an accurate estimation of these two errors Is definitely essential. Here, a local a posteriori error estimator for elastic structures with thin domain such as plate- and shell-like structures Is derived using element residuals and flux balancing technique. This method guarantees upper bounds for the global error, and also provides accurate local error Indicators for two types of errors, in the energy norm. Comparing to the classical error estimators using flux averaging technique, this shows considerably reliable and accurate effectivity indices. To illustrate the theoretical results and to verify the validity of the proposed error estimator, representative numerical examples are provided.

• 제어로봇시스템학회:학술대회논문집
• /
• 1992.10a
• /
• pp.923-926
• /
• 1992

ARMA fast Transversal filter(FTF) algorithm solves the extended least squres estimation problems in a very efficient way. But unfortunately, it exhibits a very unstable behavior, due to the accumulation of round-off errors. So, in this paper, two effective method to stabilize ARMA FTF algorithm is proposed. They are based on the analysis of the propagation of the numerical errors according to a first order linear model. The proposed methods modify the numerical properties of the variables responsible for the numerical instability, while proeserving the theoretical form of the algorithm. The proposed algorithms still have the nice complexity properties of the original algorithm, but have a much more stable brhavior.

• Journal of the Optical Society of Korea
• /
• v.16 no.4
• /
• pp.372-380
• /
• 2012

The optical method Digital Speckle Correlation Measurement (DSCM) has been extensively applied due its capability to measure the entire displacement field over a body surface. A formula of displacement measurement errors by the gradient-based DSCM method was derived. The errors were found to explicitly relate to the image grayscale errors consisting of sub-pixel interpolation algorithm errors, image noise, and subset deformation mismatch at each point of the subset. A power-law dependence of the standard deviation of displacement measurement errors on the subset size was established when the subset deformation was rigid body translation and random image noise was dominant and it was confirmed by both the numerical and experimental results. In a gradient-based algorithm the basic assumption is rigid body translation of the interrogated subsets, however, this is in contradiction to the real circumstances where strains exist. Numerical and experimental results also indicated that, subset shape function mismatch was dominant when the order of the assumed subset shape function was lower than that of the actual subset deformation field and the power-law dependence clearly broke down. The power-law relationship further leads to a simple criterion for choosing a suitable subset size, image quality, sub-pixel algorithm, and subset shape function for DSCM.

• Structural Engineering and Mechanics
• /
• v.8 no.5
• /
• pp.513-529
• /
• 1999

A concept of hierarchical modeling, the newest modeling technology, has been introduced in early 1990's. This new technology has a great potential to advance the capabilities of current computational mechanics. A first step to implement this concept is to construct hierarchical models, a family of mathematical models sequentially connected by a key parameter of the problem under consideration and have different levels in modeling accuracy, and to investigate characteristics in their numerical simulation aspects. Among representative model problems to explore this concept are elastic structures such as beam-, arch-, plate- and shell-like structures because the mechanical behavior through the thickness can be approximated with sequential accuracy by varying the order of thickness polynomials in the displacement or stress fields. But, in the numerical, analysis of hierarchical models, two kinds of errors prevail, the modeling error and the numerical approximation error. To ensure numerical simulation quality, an accurate estimation of these two errors is definitely essential. Here, a local a posteriori error estimator for elastic structures with thin domain such as plate- and shell-like structures is derived using the element residuals and the flux balancing technique. This method guarantees upper bounds for the global error, and also provides accurate local error indicators for two types of errors, in the energy norm. Compared to the classical error estimators using the flux averaging technique, this shows considerably reliable and accurate effectivity indices. To illustrate the theoretical results and to verify the validity of the proposed error estimator, representative numerical examples are provided.

• Atmosphere
• /
• v.33 no.1
• /
• pp.33-47
• /
• 2023

In this study, for the application of observation errors to the Korean Integrated Model (KIM) to utilize the Constellation Observing System for Meteorology, Ionosphere & Climate-2 (COSMIC-2) new satellites, the observation errors were diagnosed based on the Desroziers method using the cost function in the process of variational data assimilation. We calculated observation errors for all observational species being utilized for KIM and compared with their relative values. The observation error of the calculated the Global Navigation Satellite System Radio Occultation (GNSS RO) was about six times smaller than that of other satellites. In order to balance with other satellites, we conducted two experiments in which the GNSS RO data expanded by about twice the observation error. The performance of the analysis field was significantly improved in the tropics, where the COSMIC-2 data are more available, and in the Southern Hemisphere, where the influence of GNSS RO data is significantly greater. In particular, the prediction performance of the Southern Hemisphere was improved by doubling the observation error in global region, rather than doubling the COSMIC-2 data only in areas with high density, which seems to have been balanced with other observations.

• Journal of the Korean Society for Aviation and Aeronautics
• /
• v.16 no.2
• /
• pp.21-27
• /
• 2008

An inverse method is introduced to construct the problem for the error analysis of the numerical solution of initial value problem. These problems constructed through this method have a known exact solution, even though analytical solutions are generally not obtainable. The process leading to the exact solution makes use of an initially available approximate numerical solution. A smooth interpolation of the approximate solution is forced to exactly satisfy the differential equation by analytically deriving a small forcing function to absorb all of the errors in the interpolated approximate solution. Using this special case exact solution, it is possible to investigate the relationship between global errors of a candidate numerical solution process and the associated tuning parameters for a given problem. Under the assumption that the original differential equation is well-posed with respect to the small perturbations, we thereby obtain valuable information about the optimal choice of the tuning parameters and the achievable accuracy of the numerical solution.