• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.8-15
• /
• 2012

A code developed using the flux-difference splitting scheme on the hybrid Cartesian/immersed boundary method is applied to simulate three-dimensional internal waves. The material interface is regarded as a moving contact discontinuity and is captured on the basis of mass conservation without any additional treatment across the interface. Inviscid fluxes are estimated using the flux-difference splitting scheme for incompressible fluids of different density. The hybrid Cartesian/immersed boundary method is used to enforce the boundary condition for a moving three-dimensional body. Immersed boundary nodes are identified within an instantaneous fluid domain on the basis of edges crossing a boundary. The dependent variables are reconstructed at the immersed boundary nodes along local normal lines to provide the boundary condition for a discretized flow problem. The internal waves are simulated, which are generated by an pitching ellipsoid near an material interface. The effects of density ratio and location of the ellipsoid on internal waves are compared.

• Steel and Composite Structures
• /
• v.28 no.6
• /
• pp.655-670
• /
• 2018

A coupled method, that combines the Ritz method and the finite element (FE) method, is proposed to solve the vibration problem of rectangular thin and thick plates with general boundary conditions. The eigenvalue partial differential equation(s) of the plate is (are) first reduced to a set of eigenvalue ordinary differential equations by the application of the Ritz method. The resulting eigenvalue differential equations are then reduced to an eigenvalue algebraic equation system using the finite element method. The natural boundary conditions of the plate problem including the free edge and free corner boundary conditions are also implemented in a simple and accurate manner. Various boundary conditions including simply supported, clamped and free boundary conditions are considered. Comparisons with existing numerical and analytical solutions show that the proposed mixed method can produce highly accurate results for the problems considered using a small number of Ritz terms and finite elements. The proposed mixed Ritz-FE formulation is also compared with the mixed FE-Ritz formulation which has been recently proposed by the present author and his co-author. It is found that the proposed mixed Ritz-FE formulation is more efficient than the mixed FE-Ritz formulation for free vibration analysis of rectangular plates with Levy-type boundary conditions.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.2
• /
• pp.127-136
• /
• 2007

In this paper, the equations of the scaled boundary finite element method are derived for non-homogeneous half plane and analyzed numerically In the scaled boundary finite element method, partial differential equations are weaken in the circumferential direction by approximation scheme such as the finite element method, and the radial direction of equations remain in analytical form. The scaled boundary equations of non-homogeneous half plane, its elastic modulus varies as power function, are newly derived by the virtual work theory. It is shown that the governing equation of this problem is the Euler-Cauchy equation, therefore, the logarithm mode used in the half plane problem is not valid in this problem. Two numerical examples are analysed for the verification and the feasibility.

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

In this paper, we consider the multipoint boundary value problem for the one-dimensional p-Laplacian $({\phi}_p(u'))'$(t)+q(t)f(t,u(t),u'(t))=0, t $\in$ (0, 1), subject to the boundary conditions: $u(0)=\sum\limits_{i=1}^{n-2}{\alpha}_iu({\xi}_i),\;u(1)=\sum\limits_{i=1}^{n-2}{\beta}_iu({\xi}_i)$ where $\phi_p$(s) = $|s|^{n-2}s$, p > 1, $\xi_i$ $\in$ (0, 1) with 0 < $\xi_1$ < $\xi_2$ < $\cdots$ < $\xi{n-2}$ < 1 and ${\alpha}_i,\beta_i{\in}[0,1)$, 0< $\sum{\array}{{n=2}\\{i=1}}{\alpha}_i,\sum{\array}{{n=2}\\{i=1}}{\beta}_i$<1. Using a fixed point theorem due to Bai and Ge, we study the existence of at least three positive solutions to the above boundary value problem. The important point is that the nonlinear term f explicitly involves a first-order derivative.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.30 no.2A
• /
• pp.131-136
• /
• 2010

In this study, most of famous algorithms for the corner problem are listed. By comparing these with implemented codes and theoretical dissections, new algorithms are developed. These algorithms are combined by the existing auxiliary equations. All relating algorithms are numerically tested with 3 problems. Two problems have well-known analytical solutions and the result of another example is compared with the one of the published paper. The conducted research reveals the characteristics of existing algorithms and demonstrates newly developed algorithms can produce a reasonable solution by reflecting various type of boundary conditions.

• Structural Engineering and Mechanics
• /
• v.12 no.1
• /
• pp.101-118
• /
• 2001

The reliable prediction of elastic vibrations of wetted complex structures, as ships, tanks, offshore structures, propulsion components etc. represent a theoretical and numerical demanding task due to fluid-structure interaction. The paper presented is addressed to the vibration analysis by a combined FE-BE-procedure based on the added mass concept utilizing a direct boundary integral formulation of the potential fluid problem in interior and exterior domains. The discretization is realized by boundary element collocation method using conventional as well as infinite boundary element formulation with analytical integration scheme. Particular attention is devoted to modelling of interior problems with both several separate or communicating fluid domains as well as thin-walled structures wetted on both sides. To deal with this specific kind of interaction problems so-called "virtual" boundary elements in areas of cut outs are placed to satisfy the kinematical conditions in partial connected fluid domains existing in realistic tank systems. Numerical results of various theoretical and practical examples demonstrate the performance of the BE-methodology presented.

• Journal of the Korean Mathematical Society
• /
• v.37 no.1
• /
• pp.31-44
• /
• 2000

We study a free boundary problem satisfying Bernoulli type boundary condition along which the gradient of a piecewise harmonic solution jumps zero to a given constant value. In such problem, the free boundary splits the domain into two regions, the zero set and the harmonic region. Our main interest is to identify the global shape and the location of the zero set. In this paper, we find the lower and the upper bound of the zero set. In a convex domain, easier estimation of the upper bound and faster disk test technique are given to find a rough shape of the zero set. Also a simple proof on the convexity of zero set is given for a connected zero set in a convex domain.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.6
• /
• pp.829-836
• /
• 2001

A dependable boundary reconstruction technique is proposed. The finite element method is used for the analysis of the direct heat conduction problem to realize the deformable grid system. An appropriate strategy for grid update is suggested. A complete sensitivity analysis is performed to obtain the derivatives required for restoration of the optimal boundary. With the result of the sensitivity analysis, the unknown boundary is sought using the sequential quadratic programming. The method is applied to reconstruction of boundaries with sinusoidal, step, and cavity form. The overall performance of the proposed method is examined by comparison between the estimated the exact boundaries.

• Communications of the Korean Mathematical Society
• /
• v.16 no.3
• /
• pp.405-413
• /
• 2001

The inverse boundary value problem for the steady state heat equation with convection term is considered in a simply connected bounded domain with smooth boundary. Taking the Dirichlet to Neumann map which maps the temperature on the boundary to the that flux on the boundary as an observation data, the global uniqueness for identifying the convection term from the Dirichlet to Neumann map is proved.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.6
• /
• pp.1001-1004
• /
• 1987

In this study mathematical properties of variational solution and solution of the boundary integral equation of the linear elasticity problem are studied. It is first reviewed that a variational solution for the three-dimensional linear elasticity problem exists in the Sobolev space [ $H^{1}$(.OMEGA.)]$^{3}$ and, then, it is shown that a unique solution of the boundary integral equation is identical to the variational solution in [ $H^{1}$(.OMEGA.)]$^{3}$. To represent the boundary integral equation, the Green's formula in the Sobolev space is utilized on the solution domain excluding a ball, with small radius .rho., centered at the point where the point load is applied. By letting .rho. tend to zero, it is shown that, for the linear elasticity problem, boundary integral equation is valid for the variational solution. From this fact, one can obtain a numerical approximatiion of the variational solution by the boundary element method even when the classical solution does not exist.exist.