• Journal of electromagnetic engineering and science
• /
• 제11권1호
• /
• pp.11-15
• /
• 2011

Preserving high-order accuracy in high-order FDTD solutions across dielectric interfaces is very important for practical time-domain electromagnetic simulations. This paper presents a fourth-order accurate numerical boundary scheme for the planar dielectric interface to be used in the fourth-order FDTD method proposed earlier by the author. The interface scheme for the two-dimensional (2-D) transverse magnetic (TM) polarization case is derived and validated by monitoring the $L_2$ norm errors in the numerical solutions of a partially-filled cavity demonstrating its fourth-order convergence and long-time numerical stability in the presence of the planar dielectric interface.

• 대한조선학회지
• /
• 제19권4호
• /
• pp.19-29
• /
• 1982

Galerkin's finite element method is applied to a two-dimensional heat convection-diffusion problem arising in the hydrodynamic lubrication of thrust bearings used in naval vessels. A parabolized thermal energy equation for the lubricant, and thermal diffusion equations for both bearing pad and the collar are treated together, with proper juncture conditions on the interface boundaries. it has been known that a numerical instability arises when the classical Galerkin's method, which is equivalent to a centered difference approximation, is applied to a parabolic-type partial differential equation. Probably the simplest remedy for this instability is to use a one-sided finite difference formula for the first derivative term in the finite difference method. However, in the present coupled heat convection-diffusion problem in which the governing equation is parabolized in a subdomain(Lubricant), uniformly stable numerical solutions for a wide range of the Peclet number are obtained in the numerical test based on Galerkin's classical finite element method. In the present numerical convergence errors in several error norms are presented in the first model problem. Additional numerical results for a more realistic bearing lubrication problem are presented for a second numerical model.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제1권1호
• /
• pp.20-28
• /
• 2009

In thus paper we validate a numerical model for wave-structure interaction by comparing numerical results with laboratory data. The numerical model is based on the Navier-Stokes (N-S) equations for an incompressible fluid. The N-S equations are solved by a two-step projection finite volume scheme and the free surface displacements are tracked by the volume of fluid (VOF) method The numerical model is used to simulate solitary waves and their interaction with a group of slender vertical piles. Numerical results are compared with the laboratory data and very good agreement is observed for the time history of free surface displacement, fluid particle velocity and wave force. The agreement for dynamic pressure on the cylinder is less satisfactory, which is primarily caused by instrument errors.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2002년도 추계학술대회논문집
• /
• pp.623-629
• /
• 2002

The spectrum of impulse response signal from an impulse hammer testing is widely used to obtain frequency response function(FRF) of the structure. However the FRFs obtained from impact hammer testing have not only leakage errors but also finite record length errors when the record length for the signal processing is not sufficiently long. The errors cannot be removed with the conventional signal analyzer which treats the signals as if they are always steady and periodic. Since the response signals generated by the impact hammer are transient and have damping, they are undoubtedly non-periodic. It is inevitable that the signals be acquired for limited recording time, which causes the finite record length error and the leakage error. In this paper, the errors in the frequency response function of multi degree of freedom system are formulated theoretically. And the method to remove these errors is also suggested. This method is based on the optimization technique. A numerical example of 3-dof model shows the validity of the proposed method.

• 한국소음진동공학회논문집
• /
• 제12권10호
• /
• pp.792-798
• /
• 2002

The spectrum of impulse response signal from an impulse hammer testing is widely used to obtain frequency response function (FRF) of the structure. However the FRFs obtained from impact hammer testing have not only leakage errors but also finite record length errors when the record length for the signal processing is not sufficiently long. The errors cannot be removed with the conventional signal analyzer which treats the signals as if they are always steady and periodic. Since the response signals generated by the impact hammer are transient and have damping, they are undoubtedly non-periodic. It is inevitable that the signals be acquired for limited recording time, which causes the finite record length error and the leakage error. In this paper, the errors in the frequency response function of multi degree of freedom system are formulated theoretically. And the method to remove these errors is also suggested. This method is based on the optimization technique. A numerical example of 3-dof model shows the validity of the proposed method.

• Structural Engineering and Mechanics
• /
• 제16권2호
• /
• pp.219-240
• /
• 2003

One of the most versatile approaches for analyzing the dynamic behavior of structural systems is direct time integration of semi-discrete equations of motion. However responses computed by time integration are generally inexact and hence the corresponding errors would rather be studied in advance. In spite of the various error estimation formulations that exist in the literature, it is accepted practice to repeat the analyses with smaller time steps, followed by a comparison between the results. In this paper, after a review of this simple method and disregarding the round-off errors, a more efficient, reliable and yet simple method for estimating errors and enhancing the accuracy is proposed. The main objectives of this research are more realistic error estimation based on the concept of convergence, approximately controlling the reliability by comparing the actual rate of convergence with the integration method's order of accuracy, and enhancement of reliability by applying Richardson's extrapolation. Starting from the errors at specific time instants, the study is then generalized to cases in which the errors should be estimated and decreased at specific events e.g. peak responses. Numerical study illustrates the efficacy of the proposed method.

• 대기
• /
• 제25권1호
• /
• pp.1-18
• /
• 2015

The most objective way to overcome the limitation of numerical weather prediction model is to represent the uncertainty of prediction by introducing probabilistic forecast. The uncertainty of the numerical weather prediction system developed due to the parameterization of unresolved scale motions and the energy losses from the sub-scale physical processes. In this study, we focused on the growth of model errors. We performed ensemble forecast to represent model uncertainty. By employing the multi-physics scheme (PHYS) and the stochastic kinetic energy backscatter scheme (SKEBS) in simulating typhoon Rusa (2002), we assessed the performance level of the two schemes. The both schemes produced better results than the control run did in the ensemble mean forecast of the track. The results using PHYS improved by 28% and those based on SKEBS did by 7%. Both of the ensemble mean errors of the both schemes increased rapidly at the forecast time 84 hrs. The both ensemble spreads increased gradually during integration. The results based on SKEBS represented model errors very well during the forecast time of 96 hrs. After the period, it produced an under-dispersive pattern. The simulation based on PHYS overestimated the ensemble mean error during integration and represented the real situation well at the forecast time of 120 hrs. The displacement speed of the typhoon based on PHYS was closest to the best track, especially after landfall. In the sensitivity tests of the model uncertainty of SKEBS, ensemble mean forecast was sensitive to the physics parameterization. By adjusting the forcing parameter of SKEBS, the default experiment improved in the ensemble spread, ensemble mean errors, and moving speed.

• Kyungpook Mathematical Journal
• /
• 제62권4호
• /
• pp.699-713
• /
• 2022

This paper proposes a hybrid successive over-relaxation iterative method for the numerical solution of a nonlinear diffusion in a one-dimensional porous medium. The considered mathematical model is discretized using a computational complexity reduction scheme called half-sweep finite differences. The local truncation error and the analysis of the stability of the scheme are discussed. The proposed iterative method, which uses explicit group technique and modified successive over-relaxation, is formulated systematically. This method improves the efficiency of obtaining the solution in terms of total iterations and program elapsed time. The accuracy of the proposed method, which is measured using the magnitude of absolute errors, is promising. Numerical convergence tests of the proposed method are also provided. Some numerical experiments are delivered using initial-boundary value problems to show the superiority of the proposed method against some existing numerical methods.

• Journal of applied mathematics & informatics
• /
• 제42권2호
• /
• pp.321-331
• /
• 2024

We consider the numerical treatment of blow-up problems having non-monotonic singular solutions that tend to infinity at some point in the domain. The use of standard numerical methods for solving problems with blow-up solutions can lead to significant errors. The reason being that solutions of such problems have singularities whose positions are unknown in advance. To be able to integrate such non-monotonic blow-up problems, we describe and use a method of non-local transformations. To show the efficiency of the method, we present a comparison of exact and numerical solutions in addition to some comparison with the work of other authors.

• Structural Engineering and Mechanics
• /
• 제33권1호
• /
• pp.15-30
• /
• 2009

In this study, the nonlinear vibrations of stepped beams having different boundary conditions were investigated. The equations of motions were obtained by using Hamilton's principle and made non dimensional. The stretching effect induced non-linear terms to the equations. Natural frequencies are calculated for different boundary conditions, stepped ratios and stepped locations by Newton-Raphson Method. The corresponding nonlinear correction coefficients are also calculated for the fundamental mode. At the second part, an alternative method is produced for the analysis. The calculated natural frequencies and nonlinear corrections are used for training an artificial neural network (ANN) program which has a multi-layer, feed-forward, back-propagation algorithm. The results of the algorithm produce errors less than 2.5% for linear case and 10.12% for nonlinear case. The errors are much lower for most cases except clamped-clamped end condition. By employing the ANN algorithm, the natural frequencies and nonlinear corrections are easily calculated by little errors, and the computational time is drastically reduced compared with the conventional numerical techniques.