• 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.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2004.04a
• /
• pp.57-62
• /
• 2004

Streamline simulation researches have been extensively accomplished due to the swiftness of computation and the reduction of numerical dispersion. In this study, we developed a streamline simulation model using a semianalytical solution of ID transport equation. To validate accuracy of the developed model, we compared simulation results of contaminant transport, which were acquired by streamline simulation models using an analytical solution, a numerical solution, and a semianalytical solution. The developed model using the semianalytical solution matched well with the model using an analytical solution. However, streamline simulation model using a numerical solution showed numerical dispersion. For an advection-dominant flow, there was little difference in the simulation results between the developed model and tile analytical model, but the differences between the analytical model and the numerical model were cleary shown. From the comparison of computing time we know that the streamline simulation using the semianalytical solution is 2-60 times as fast as the streamline simulation using the numerical solution.

• International Journal of Aeronautical and Space Sciences
• /
• v.3 no.2
• /
• pp.46-58
• /
• 2002

An inverse method is introduced to construct benchmark problems for the numerical solution of initial value problems. Benchmark problems constructed through this method have a known exact solution, even though analytical solutions are generally not obtainable. The solution is constructed such that it lies near a given approximate numerical solution, and therefore the special case solution can be generated in a versatile and physically meaningful fashion and can serve as a benchmark problem to validate approximate solution methods. 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. A multi-variable orthogonal function expansion method and computer symbol manipulation are successfully used for this process. Using this special case exact solution, it is possible to directly investigate the relationship between global errors of a candidate numerical solution process and the associated tuning parameters for a given code and 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. Illustrative examples show the utility of this method not only for the ordinary differential equations (ODEs) but for the partial differential equations (PDEs).

• 제어로봇시스템학회:학술대회논문집
• /
• 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.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.05a
• /
• pp.73-80
• /
• 2006

In order to examine the validity of an asymptotic solution for nonlinear interaction in asymmetric vibration modes of a perfect circular plate, we obtain the numerical solution. The motion of the plate is governed by nonlinear partial differential equation. The initial and boundary value problem is solved by using the finite difference method. The numerical solution is compared with the asymptotic solution. It is found that traveling waves relating clockwise and counterclockwise as well as standing wave are depicted by the numerical solution.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.8 no.5
• /
• pp.434-441
• /
• 2016

Numerical simulations of 1D Burgers equation and 2D sloshing problem were carried out to study numerical discrepancy between serial and parallel computations. The numerical domain was decomposed into 2 and 4 subdomains for parallel computations with message passing interface. The numerical solution of Burgers equation disclosed that fully explicit boundary conditions used on subdomains of parallel computation was responsible for the numerical discrepancy of transient solution between serial and parallel computations. Two dimensional sloshing problems in a rectangular domain were solved using OpenFOAM. After a lapse of initial transient time sloshing patterns of water were significantly different in serial and parallel computations although the same numerical conditions were given. Based on the histograms of pressure measured at two points near the wall the statistical characteristics of numerical solution was not affected by the number of subdomains as much as the transient solution was dependent on the number of subdomains.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.5
• /
• pp.461-471
• /
• 2008

An analytical solution of one dimensional mild slope equation is derived by use of the WKB method, which has a form similar to Porter's solution(2003). The present solution is so general in the sense of application that it is comparable to the corresponding numerical solutions. In the derivation we also presented the solution of refraction equation in terms of surface displacement. Some numerical results of the present solution by use of Bremmer's method are presented which agree with existing numerical solutions.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.26 no.7
• /
• pp.795-805
• /
• 2016

Nonlinear normal modes(NNMs) is a branch of periodic solution of nonlinear dynamic systems. Determination of stable periodic solution is very important in many engineering applications since the stable periodic solution can be an attractor of such nonlinear systems. Periodic solutions of nonlinear system are usually calculated by perturbation methods and numerical methods. In this study, numerical method is used in order to calculate the NNMs. Iteration of the solution is presented by multiple shooting method and continuation of solution is presented by pseudo-arclength continuation method. The stability of the NNMs is analyzed using Floquet multipliers, and bifurcation points are calculated using indirect method. Proposed analyses are applied to two nonlinear numerical models. In the first numerical model nonlinear spring-mass system is analyzed. In the second numerical model Jeffcott rotor system which has unstable equilibria is analyzed. Numerical simulation results show that the multiple shooting method can be applied to self excited system as well as the typical nonlinear system with stable equilibria.

• Water Engineering Research
• /
• v.2 no.3
• /
• pp.187-196
• /
• 2001

An approximate analytical solution of a nonlinear hydro-thermo coupled diffusion equation is derived using the dimensionless form of the equation and transformation method. To derive an analytical solution, it is drastically assumed that the product of first order derivatives in the non-dimensionalized governing equation has little influence on the solution of heat and moisture behavior problem. The validity of this drastic assumption is demonstrated. Some numerical simulation is performed to investigate the applicability of a derived approximate analytical solution. The results show a good agreement between analytical and numerical solutions. The proposed solution may provide a useful tool in the verification process of the numerical models. Also, the solution can be used for the analysis of one-dimensional coupled heat and moisture movements in unsaturated porous media.

• Structural Engineering and Mechanics
• /
• v.34 no.1
• /
• pp.85-95
• /
• 2010

In this paper, numerical solution of the singular integral equation for the multiple curved branch-cracks is investigated. If some quadrature rule is used, one difficult point in the problem is to balance the number of unknowns and equations in the solution. This difficult point was overcome by taking the following steps: (a) to place a point dislocation at the intersecting point of branches, (b) to use the curve length method to covert the integral on the curve to an integral on the real axis, (c) to use the semi-open quadrature rule in the integration. After taking these steps, the number of the unknowns is equal to the number of the resulting algebraic equations. This is a particular advantage of the suggested method. In addition, accurate results for the stress intensity factors (SIFs) at crack tips have been found in a numerical example. Finally, several numerical examples are given to illustrate the efficiency of the method presented.