• Coupled systems mechanics
• /
• v.4 no.1
• /
• pp.1-18
• /
• 2015

The state-based peridynamics is considered a nonlocal method in which the equations of motion utilize integral form as opposed to the partial differential equations in the classical continuum mechanics. As a result, the enforcement of boundary conditions in solid mechanics analyses cannot follow the standard way as in a classical continuum theory. In this paper, a new approach for the boundary condition enforcement in the state-based peridynamic formulation is presented. The new method is first formulated based on a convex kernel approximation to restore the Kronecker-delta property on the boundary in 1-D case. The convex kernel approximation is further localized near the boundary to meet the condition that recovers the correct boundary particle forces. The new formulation is extended to the two-dimensional problem and is shown to reserve the conservation of linear momentum and angular momentum. Three numerical benchmarks are provided to demonstrate the effectiveness and accuracy of the proposed approach.

• Interaction and multiscale mechanics
• /
• v.7 no.1
• /
• pp.525-542
• /
• 2014

The state-based peridynamics is considered a nonlocal method in which the equations of motion utilize integral form as opposed to the partial differential equations in the classical continuum mechanics. As a result, the enforcement of boundary conditions in solid mechanics analyses cannot follow the standard way as in a classical continuum theory. In this paper, a new approach for the boundary condition enforcement in the state-based peridynamic formulation is presented. The new method is first formulated based on a convex kernel approximation to restore the Kronecker-delta property on the boundary in 1-D case. The convex kernel approximation is further localized near the boundary to meet the condition that recovers the correct boundary particle forces. The new formulation is extended to the two-dimensional problem and is shown to reserve the conservation of linear momentum and angular momentum. Three numerical benchmarks are provided to demonstrate the effectiveness and accuracy of the proposed approach.

• Journal of Mechanical Science and Technology
• /
• v.14 no.1
• /
• pp.84-92
• /
• 2000

A new bounce back boundary method of the second order in error is proposed for the lattice Boltzmann fluid simulation. This new method can be used for the arbitrarily irregular lattice geometry of a non-slip boundary. The traditional bounce back boundary condition for the lattice Boltzmann simulation is of the first order in error. Since the lattice Boltzmann method is the second order scheme by itself, a boundary technique of the second order has been desired to replace the first order bounce back method. This study shows that, contrary to the common belief that the bounce back boundary condition is unilaterally of the first order, the second order bounce back boundary condition can be realized. This study also shows that there exists a generalized bounce back technique that can be characterized by a single interpolation parameter. The second order bounce back method can be obtained by proper selection of this parameter in accordance with the detailed lattice geometry of the boundary. For an illustrative purpose, the transient Couette and the plane Poiseuille flows are solved by the lattice Boltzmann simulation with various boundary conditions. The results show that the generalized bounce back method yields the second order behavior in the error of the solution, provided that the interpolation parameter is properly selected. Coupled with its intuitive nature and the ease of implementation, the bounce back method can be as good as any second order boundary method.

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

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.2
• /
• pp.113-120
• /
• 2021

The wave-propagation phenomenon in an infinite medium has been used to describe the physics in many fields of engineering and natural science. Analytical or numerical methods have been developed to obtain solutions to problems related to the wave-propagation phenomenon. Energy radiation into infinite regions must be accurately considered for accurate solutions to these problems; hence, various numerical and mechanical models as well as boundary conditions have been developed. This paper proposes a new boundary condition that can be applied to scalar-wave or horizontally-polarized shear-wave (or SH-wave) propagation problems in layered waveguides. A governing equation is obtained for the SH waves by applying finite-element discretization in the vertical direction of the waveguide and subsequently modified to derive the boundary condition for the infinite region of the waveguide. Using the orthogonality of the eigenmodes for the SH waves in a layered waveguide, the new boundary condition is shown to be equivalent to the existing root-finding absorbing boundary condition; further, the accuracy is shown to increase with the degree of the new boundary condition, and its stability can be proven. The accuracy and stability are then demonstrated by applying the proposed boundary condition to wave-propagation problems in layered waveguides.

• Structural Engineering and Mechanics
• /
• v.78 no.4
• /
• pp.497-518
• /
• 2021

In this study, a new method for treating the wall boundary in smoothed particle hydrodynamics (SPH) is proposed to simulate free surface flows effectively. Unlike conventional methods of wall boundary treatment through boundary particles, in the proposed method, the wall boundary condition is directly imposed by adding boundary truncation terms to the mass and momentum conservation equations. Thus, boundary particles are not used in boundary modeling. Doing so, the wall boundary condition is accurately imposed, boundary modeling is simplified, and computation is made efficient without losing stability in SPH. Performance of the proposed method is demonstrated through several numerical examples: dam break, dam break with a wedge, sloshing, inclined bed, cross-lever rotation, pulsating tank and sloshing with a flexible baffle. These results are compared with available experimental results, analytical solutions, and results obtained using the boundary particle method.

• Journal of the Korean Society for Precision Engineering
• /
• v.28 no.8
• /
• pp.894-903
• /
• 2011

In this paper, the capability index is introduced in order to improve the reliability of new concept machine tools and the method to improve the machine accuracy from the analysis of cutting process, statistical methodology and influence factors are proposed. In addition, the rib structure of bed and column in machine tools is analyzed by using the thermal impact method in order to analogize the rib pattern which has the small thermal deformation under thermal boundary condition. In the analysis of column rib structure, thermal boundary condition is separated to heat conduction and heat transfer to appropriate real boundary condition. Finally, performance chart of bed and column rib structure is provided for designer to estimate each rib pattern and select rib structure appropriating to thermal condition.

• Earthquakes and Structures
• /
• v.9 no.1
• /
• pp.173-194
• /
• 2015

Site response analysis is an important topic in earthquake engineering. A time-domain numerical method called as one-dimensional (1D) finite element artificial boundary method is proposed to simulate the homogeneous plane elastic wave propagation in a layered half space subjected to the obliquely incident plane body wave. In this method, an exact artificial boundary condition combining the absorbing boundary condition with the inputting boundary condition is developed to model the wave absorption and input effects of the truncated half space under layer system. The spatially two-dimensional (2D) problem consisting of the layer system with the artificial boundary condition is transformed equivalently into a 1D one along the vertical direction according to Snell's law. The resulting 1D problem is solved by the finite element method with a new explicit time integration algorithm. The 1D finite element artificial boundary method is verified by analyzing two engineering sites in time domain and by comparing with the frequency-domain transfer matrix method with fast Fourier transform.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2797-2802
• /
• 2007

It has been confirmed that implementation of the no-slip boundary conditions for the lattice-Boltzmann method play an important role in the overall accuracy of the numerical solutions as well as the stability of the solution procedure. We in this paper propose a new algorithm, i.e. the method of the dynamic boundary condition for no-slip boundary condition. The distribution functions on the wall along each of the links across the physical boundary are assumed to be composed of equilibrium and nonequilibrium parts which inherit the idea of Guo's extrapolation method. In the proposed algorithm, we apply a dynamic equation to reflect the computational slip velocity error occurred on the actual wall boundary to the correction; the calculated slip velocity error dynamically corrects the fictitious velocity on the wall nodes which are subsequently employed to the computation of equilibrium distribution functions on the wall nodes. Along with the dynamic selfcorrecting process, the calculation efficiently approaches the steady state. Numerical results show that the dynamic boundary method is featured with high accuracy and simplicity.

• Journal of applied mathematics & informatics
• /
• v.27 no.5_6
• /
• pp.1109-1118
• /
• 2009

In this paper, we are concerned with the following four point boundary value problem with one-dimensional p-Laplacian, $\{({\phi}_p(x'(t)))'+h(t)f(t,x(t),|x'(t)|)=0$, 0< t<1, $x'(0)-{\delta}x(\xi)=0,\;x'(1)+{\delta}x(\eta)=0$, where $\phi_p$ (s) = |s|$^{p-2}$, p > $\delta$ > 0, 1 > $\eta$ > $\xi$ > 0, ${\xi}+{\eta}$ = 1. By using a fixed point theorem in a cone, we obtain the existence of at least three symmetric positive solutions. The interesting point is that the boundary condition is a new Sturm-Liouville-like boundary condition, which has rarely been treated up to now.