• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.25 no.10
• /
• pp.1087-1094
• /
• 2014

An numerical technique of impedance boundary condition to improve an efficiency in the process of moment method with CFIE(Combined Field Integral Equation), which is widely used to analyze the scattering property of dielectric scatterers, and results of its cross-validations are presented in this study. Application of the impedance boundary allows to represent the equivalent surface currents of dielectric scatterer depicted by both kinds of electric/magnetic surface currents(Js, Ms) to the single surface current by Js or Ms only. Accuracy of this technique is validated by the existing CFIE and theoretical values such as Mie-series solution and small perturbation scattering model. The computational difference of less than 1 dB was verified within an imaginary part of dielectric constant more than 12, as well.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.6
• /
• pp.447-457
• /
• 2018

In this study, a two-dimensional fully nonlinear transient wave numerical tank was developed using a desingularized indirect boundary integral equation method. The desingularized indirect boundary integral equation method is simpler and faster than the conventional boundary element method because special treatment is not required to compute the boundary integral. Numerical simulations were carried out in the time domain using the fourth order Runge-Kutta method. A mixed Eulerian-Lagrangian approach was adapted to reconstruct the free surface at each time step. A numerical damping zone was used to minimize the reflective wave in the downstream region. The interpolating method of a Gaussian radial basis function-type artificial neural network was used to calculate the gradient of the free surface elevation without element connectivity. The desingularized indirect boundary integral equation using an isolated point source and radial basis function has no need for information about the element connectivity and is a meshless method that is numerically more flexible. In order to validate the accuracy of the numerical wave tank based on the desingularized indirect boundary integral equation method and meshless technique, several numerical simulations were carried out. First, a comparison with numerical results according to the type of desingularized source was carried out and confirmed that continuous line sources can be replaced by simply isolated sources. In addition, a propagation simulation of a $2^{nd}$-order Stokes wave was carried out and compared with an analytical solution. Finally, simulations of propagating waves in shallow water and propagating waves over a submerged bar were also carried and compared with published data.

• The Pure and Applied Mathematics
• /
• v.17 no.1
• /
• pp.39-49
• /
• 2010

This paper deals with the second-order differential equation (p(t)x'(t))' + g(t)f(t, x(t), x'(t)) = 0, a.e. in (0, $\infty$) with the boundary conditions $$x(0)={\int}^{\infty}_0g(s)x(s)ds,\;{lim}\limits_{t{\rightarrow}{\infty}}p(t)x'(t)=0,$$ where $g\;{\in}\;L^1[0,{\infty})$ with g(t) > 0 on [0, $\infty$) and ${\int}^{\infty}_0g(s)ds\;=\;1$, f is a g-Carath$\acute{e}$odory function. By applying the coincidence degree theory, the existence of at least one solution is obtained.

• The Bulletin of The Korean Astronomical Society
• /
• v.44 no.1
• /
• pp.46.1-46.1
• /
• 2019

We present an accurate and efficient method to calculate the gravitational potential of an isolated system in three-dimensional Cartesian and cylindrical coordinates subject to vacuum (open) boundary conditions. Our method consists of two parts: an interior solver and a boundary solver. The interior solver adopts an eigenfunction expansion method together with a tridiagonal matrix solver to solve the Poisson equation subject to the zero boundary condition. The boundary solver employs James's method to calculate the boundary potential due to the screening charges required to keep the zero boundary condition for the interior solver. A full computation of gravitational potential requires running the interior solver twice and the boundary solver once. We develop a method to compute the discrete Green's function in cylindrical coordinates, which is an integral part of the James algorithm to maintain second-order accuracy. We implement our method in the {\tt Athena++} magnetohydrodynamics code, and perform various tests to check that our solver is second-order accurate and exhibits good parallel performance.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.15-28
• /
• 1998

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

• Kyungpook Mathematical Journal
• /
• v.63 no.3
• /
• pp.437-450
• /
• 2023

In this paper, we present a sufficient condition for the unique existence of solutions for a coupled system of nonlinear fractional Langevin equations with a new class of multipoint and nonlocal integral boundary conditions. We define a 𝓩^*_λ-contraction mapping and present the sufficient condition by identifying the problem with an equivalent fixed point problem in the context of b-metric spaces. Finally, some numerical examples are given to validate our main results.

• Journal of applied mathematics & informatics
• /
• v.30 no.1_2
• /
• pp.27-47
• /
• 2012

It has become apparent from the recent work by Choi et al. [3] on the nonlinear beam deflection problem, that analysis of the integral operator $\mathcal{K}$ arising from the beam deflection equation on linear elastic foundation is important. Motivated by this observation, we perform investigations on the eigenstructure of the linear integral operator $\mathcal{K}_l$ which is a restriction of $\mathcal{K}$ on the finite interval [$-l,l$]. We derive a linear fourth-order boundary value problem which is a necessary and sufficient condition for being an eigenfunction of $\mathcal{K}_l$. Using this equivalent condition, we show that all the nontrivial eigenvalues of $\mathcal{K}l$ are in the interval (0, 1/$k$), where $k$ is the spring constant of the given elastic foundation. This implies that, as a linear operator from $L^2[-l,l]$ to $L^2[-l,l]$, $\mathcal{K}_l$ is positive and contractive in dimension-free context.

• Journal of Mechanical Science and Technology
• /
• v.18 no.9
• /
• pp.1549-1558
• /
• 2004

This study is concerned with the general solution of the field intensity factors and energy release rate for a Griffith crack in a piezoelectric ceramic of finite radius under combined anti-plane mechanical and in-plane electrical loading. Both electrically continuous and impermeable crack surface conditions are considered. Employing Mellin transforms and Fourier series, the problem is reduced to dual integral forms. The solution to the resulting expressions is expressed in terms of Fredholm integral equation of the second kind. The solutions are provided to study the influence of the crack length, the crack surface boundary conditions on the intensity factors and the energy release rate.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.4
• /
• pp.775-786
• /
• 2002

A Mixed Volume and Boundary Integral Equation Method is applied for the effective analysis of elastic wave scattering problems and plane elastostatic problems in unbounded solids containing general anisotropic inclusions and voids or isotropic inclusions. It should be noted that this newly developed numerical method does not require the Green's function for anisotropic inclusions to solve this class of problems since only Green's function for the unbounded isotropic matrix is involved in their formulation for the analysis. This new method can also be applied to general two-dimensional elastodynamic and elastostatic problems with arbitrary shapes and number of anisotropic inclusions and voids or isotropic inclusions. In the formulation of this method, the continuity condition at each interface is automatically satisfied, and in contrast to finite element methods, where the full domain needs to be discretized, this method requires discretization of the inclusions only. Finally, this method takes full advantage of the pre- and post-processing capabilities developed in FEM and BIEM. Through the analysis of plane elastostatic problems in unbounded isotropic matrix with orthotropic inclusions and voids or isotropic inclusions, and the analysis of plane wave scattering problems in unbounded isotropic matrix with isotropic inclusions and voids, it will be established that this new method is very accurate and effective for solving plane wave scattering problems and plane elastic problems in unbounded solids containing general anisotropic inclusions and voids/cracks or isotropic inclusions.

• Journal of the Korean Mathematical Society
• /
• v.41 no.2
• /
• pp.345-368
• /
• 2004

A preconditioning algorithm is developed in this paper for the iterative solution of the linear system of equations resulting from the p-version boundary element approximation of the three dimensional integral equation with hypersingular operators. The preconditioner is derived by first making the nodal and side basis functions locally orthogonal to the element internal bases, and then by decoupling the nodal and side bases from the internal bases. Its implementation consists of solving a global problem on the wire-basket and a series of local problems defined on a single element. Moreover, the condition number of the preconditioned system is shown to be of order $O((1＋ln/p)^{7})$. This technique can be applied to discretization with triangular elements and with general basis functions.