• Journal of Mechanical Science and Technology
• /
• v.15 no.5
• /
• pp.671-680
• /
• 2001

Boundary condition is one of the major factors to influence the numerical stability and solution accuracy in numerical analysis. One of the most important physical boundary conditions in the flowfield analysis is the wall boundary condition imposed on the body surface. To solve a two-dimensional Euler equation, totally four numerical wall boundary conditions should be prescribed. Two of them are supplied by the flow tangency condition. The other two conditions, therefore, should be prepared additionally in a suitable way. In this paper, four different sets of wall boundary conditions are proposed and then applied to solve high-speed flowfields around a quarter circle geometry. A two-dimensional compressible Euler solver is prepared based on the finite volume method. This solver hires three different upwind schemes; Steger-Warmings flux vector splitting, Roes flux difference splitting, and Lious advection upstream splitting method. It is found that the way to specify the additional numerical wall boundary conditions strongly affects the overall stability and accuracy of the upwind schemes in high-speed flow calculation. The optimal wall boundary conditions should be also chosen very carefully depending on the numerical schemes used to solve the problem.

• Journal of the Korean Institute of Intelligent Systems
• /
• v.21 no.6
• /
• pp.673-678
• /
• 2011

This paper analyzes the standing balance of four-legged robots which are useful for delivering objects or investigating of information. For this, we specify an effective model of general four-legged robots and propose a boundary condition based on the standing stability of the four-legged walking. To verify such a standing balance, we consider some exemplary free motions at the standing mode of the robot and discuss on the robot's balance margin. The analysis specified in this paper will be applicable for effective balancing control of various quadruped robotic walking.

• Journal of the Korean Mathematical Society
• /
• v.58 no.5
• /
• pp.1109-1129
• /
• 2021

This paper considers an anisotropic polytropic infiltration equation with a source term $$u_t={\sum\limits_{i=1}^{N}}{\frac{{\partial}}{{\partial}x_i}}\(a_1(x){\mid}u{\mid}^{{\alpha}_i}{\mid}u_{x_i}{\mid}^{p_i-2}u_{x_i}\)+f(x,t,u)$$, where p_i > 1, α_i > 0, a_i(x) ≥ 0. The existence of weak solution is proved by parabolically regularized method. Based on local integrability $u_{x_i}{\in}W_{loc}^{1,p_i}(\Omega)$, the stability of weak solutions is proved without boundary value condition by the weak characteristic function method. One of the essential characteristics of an anisotropic equation different from an isotropic equation is found originally.

• Communications of the Korean Mathematical Society
• /
• v.30 no.3
• /
• pp.297-311
• /
• 2015

We present a robust and accurate boundary condition for pricing financial options that is a hybrid combination of the payoff-consistent extrapolation and the Dirichlet boundary conditions. The payoff-consistent extrapolation is an extrapolation which is based on the payoff profile. We apply the new hybrid boundary condition to the multi-dimensional Black-Scholes equations with a high correlation. Correlation terms in mixed derivatives make it more difficult to get stable numerical solutions. However, the proposed new boundary treatments guarantee the stability of the numerical solution with high correlation. To verify the excellence of the new boundary condition, we have several numerical tests such as higher dimensional problem and exotic option with nonlinear payoff. The numerical results demonstrate the robustness and accuracy of the proposed numerical scheme.

• Steel and Composite Structures
• /
• v.11 no.5
• /
• pp.359-374
• /
• 2011

In this paper, the influence of fibre orientation and aspect ratio on stability analysis of simply supported skew plates subjected to in plane loading is studied by using a four noded hybrid plate finite element. The formulation of the element is based on Hellinger-Reissner variational principle. The element is developed by combining a hybrid plane stress element and a hybrid plate element. Some numerical problems are solved and the effects of skew angle, aspect ratio, fibre orientation and loading type on the critical buckling loads are highlighted.

• International Journal of Aeronautical and Space Sciences
• /
• v.13 no.1
• /
• pp.43-57
• /
• 2012

In a modern aircraft, there are many variations in its mass, stiffness, and aerodynamic characteristics. Recently, an analytical approach was proposed, and this approach uses the idea of uncertainty to find out the most critical flight flutter boundary due to the variations in such aerodynamic characteristics. An analytical method that has been suggested to predict robust stability is the mu method. We previously analyzed the robust flutter boundary by using the mu method, and in that study, aerodynamic variations in the Mach number, atmospheric density, and flight speed were taken into consideration. The authors' previous attempt and the results are currently quoted as varying Mach number mu analysis. In the author's previous method, when the initial flight conditions were located far from the nominal flutter boundary, conservative predictions were obtained. However, relationships among those aerodynamic parameters were not applied. Thus, the varying Mach number mu analysis results required validation. Using an optimization approach, the varying Mach number mu analysis was found out to be capable of capturing a reasonable robust flutter boundary, i.e., with a low percentage difference from boundaries that were obtained by optimization. Regarding the optimization approach, a discrete nominal flutter boundary is to be obtained in advance, and based on that boundary, an interpolated function was established. Thus, the optimization approach required more computational effort for a larger number of uncertainty variables. And, this produced results similar to those from the mu method which had lower computational complexity. Thus, during the estimation of robust aeroelastic stability, the mu method was regarded as more efficient than the optimization method was. The mu method predicts reasonable results when an initial condition is located near the nominal flutter boundary, but it does not consider the relationships that are among the aerodynamic parameters, and its predictions are not very accurate when the initial condition is located far from the nominal flutter boundary. In order to provide predictions that are more accurate, the relationships among the uncertainties should also be included in the mu method.

• Fire Science and Engineering
• /
• v.29 no.4
• /
• pp.33-38
• /
• 2015

Loads applied on the floor are transferred through beams to columns. The beams can be designed as both end fixed or simple beams. The load bearing capacity of a beam depends on each boundary condition. However, when the load bearing capacity of a beam is evaluated in fire tests, all kinds of beams are tested using simple beam conditions. In this study, an analytical method performed using heat transfer theory and heat stress analysis based on the mechanical and thermal properties of SS-400 steel at high temperature. This method was used to clarify the differences between the two types of boundary conditions at normal and high temperature. The results show that the load bearing capacity of a both-end fixed beam at high temperature is superior to that of a simple beam. Therefore, the application of simple beam conditions in fire tests for evaluation of load bearing capacity is conservatively safe compared to fixed boundary conditions.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.3
• /
• pp.179-202
• /
• 2019

This paper presents second derivative generalized extended backward differentiation formulas (SDGEBDFs) based on the second derivative linear multi-step formulas of Cash [1]. This class of second derivative linear multistep formulas is implemented as boundary value methods on stiff problems. The order, error constant and the linear stability properties of the new methods are discussed.

• Tunnel and Underground Space
• /
• v.19 no.6
• /
• pp.498-506
• /
• 2009

The stability assessment of a concrete liner of a compressed air storage tunnel should be performed by an approach which is different from that commonly used for the liners of road tunnels, since the liner is exposed to high air pressure. In this study, the stability analysis method for the liner of compressed air storage tunnel is proposed based on the elastic and elasto-plastic solutions of the thick-walled cylinder problem. In case of elastic analysis, the yield initiation condition at the inner boundary is considered as the failure condition of the liner, while the condition which results in the extension of yielding zone to a certain depth is taken as a failure indicator of the liner in the elasto-plastic analysis taking Mohr-Coulomb criterion. The application of the proposed method revealed that the influence of the relative magnitude of boundary loads on the stability of liner is considerable. In particular, noting that the estimation of the outer boundary load may be relatively difficult, it is thought that the precise prediction of outer boundary load is very important in the analysis. Accordingly, the emphasis is put on the selection of the liner installation time, which may govern the magnitude of outer boundary load.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.10a
• /
• pp.244-251
• /
• 2002

The purpose of this paper is to investigate the stability of tapered columns with general boundary condition(translational and rotational elastic support) at one end and carrying a tip mass of rotatory inertia with translational elastic support at the other end. The column model is based on the classical Bernoulli-Euler beam theory which neglects the effects of rotatory inertia and shear deformation. The governing differential equation for the free vibrations of linearly tapered columns subjected to a subtangential follower force is solved numerically using the corresponding boundary conditions. And the bisection method is used to calculate the critical divergence/flutter load. After having verified the results of the present study, the frequency and critical divergence/flutter load are presented as functions of various nondimensional system parameters.