• Advances in Computational Design
• /
• v.5 no.4
• /
• pp.363-380
• /
• 2020

In this paper, a cylindrical shell is immersed in a non-viscous fluid using first order shell theory of Sander. These equations are partial differential equations which are solved by approximate technique. Robust and efficient techniques are favored to get precise results. Employment of the Rayleigh-Ritz procedure gives birth to the shell frequency equation. Use of acoustic wave equation is done to incorporate the sound pressure produced in a fluid. Hankel's functions of second kind designate the fluid influence. Mathematically the integral form of the Lagrange energy functional is converted into a set of three partial differential equations. Throughout the computation, simply supported edge condition is used. Expressions for modal displacement functions, the three unknown functions are supposed in such way that the axial, circumferential and time variables are separated by the product method. Comparison is made for empty and fluid-filled cylindrical shell with circumferential wave number, length- and height-radius ratios, it is found that the fluid-filled frequencies are lower than that of without fluid. To generate the fundamental natural frequencies and for better accuracy and effectiveness, the computer software MATLAB is used.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.24 no.4
• /
• pp.354-361
• /
• 2004

This Paper describes a crack scattering model for SH wave based on the boundary integral equation(BIE) method, where the fundamental unknown is crack opening displacement(COD). When a time harmonic plane wave was incident on a 2-D isolated crack (slit) of width 2a, the COD distributions were numerically calculated as a function of ka. The calculated COD agreed well with results obtained with other methods. The far-field scattering amplitude, which completely characterizes the flaw response, was calculated in two ways. The Kirchhoff approximation and the BIE-COD exact formulation were compared in terms of incidence angle and frequency ka in a pulse-echo mode. Maximum response was obtained for both methods at the specular reflection direction. Away from the specular direction, the Kirchhoff approximation becomes less accurate. The time domain crack response was also calculated using a band-limited spectrum of center frequency 10 MHz. At oblique incidence to the crack both methods show the existence of an antisymmetric flash points occurring from the crack edge. The Kirchhoff approximation provides an exact time interval between flash points, although it unrealistically gives the same amplitude.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.22 no.10
• /
• pp.1005-1011
• /
• 2011

In this paper, we present a marching-on-in-degree(MOD) finite difference method(FDM) based on the Helmholtz wave equation for analyzing transient electromagnetic responses in a general dispersive media. The two issues related to the finite difference approximation of the time derivatives and the time consuming convolution operations are handled analytically using the properties of the Laguerre functions. The basic idea here is that we fit the transient nature of the fields, the flux densities, the permittivity with a finite sum of orthogonal Laguerre functions. Through this novel approach, not only the time variable can be decoupled analytically from the temporal variations but also the final computational form of the equations is transformed from finite difference time-domain(FDTD) to a finite difference formulation through a Galerkin testing. Representative numerical examples are presented for transient wave propagation in general Debye, Drude, and Lorentz dispersive medium.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.4
• /
• pp.15-26
• /
• 1990

This paper describes a potential-baoed panel method formulated for the analysis cf a supercavitating two-dimensional symmetri strut. The method employs normal dipoles and sources distributed on the foil and cavity surfaces to represent the potential flow around the cavitating hydrofoil. The kinematic boundary condition on the wetted portion of the foil surface is satisfied by requiring that the total potential vanish in the fictitious inner flow region of the foil, and the dynamic boundary condition on the cavity surface is satisfied by requiring that the potential vary linearly, i.e., the tangential velocity be constant. Green's theorem then results in a potential-based integral equation rather than the usual velocity-based formulation of Hess & Smith type, With the singularities distributed on the exact hydrofoil surface, the pressure distributions are predicted with improved accuracy compared to those of the linearized lifting surface theory, especially near the leading edge. The theory then predicts the cavity shape and cavitation number for an assumed cavity length. To improve the accuracy, the sources and dipoles on the cavity surface are moved to the newly computed cavity surface, where the boundary conditions are satisfied again. This iteration process is repeated until the results are converged.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.2
• /
• pp.159-173
• /
• 1991

This paper describes a potential-based panel method formulated for the analysis of a super-cavitating two-dimensional hydrofoil. The method employs normal dipoles and sources distributed on the foil and cavity surfaces to represent the potential flow around the cavitating hydrofoil. The kinematic boundary condition on the wetted portion of the foil surface is satisfied by requiring that the total potential vanish in the fictitious inner flow region of the foil, and the dynamic boundary condition on the cavity surface is satisfied by requiring thats the potential vary linearly, i.e., the tangential velocity be constant. Green's theorem then results in a potential-based integral equation rather than the usual velocity-based formulation of Hess & Smith type. With the singularities distributed on the exact hydrofoil surface, the pressure distributions are predicted with improved accuracy compared to those of the linearized lilting surface theory, especially near the leading edge. The theory then predicts the cavity shape and cavitation number for an assumed cavity length. To improve the accuracy, the sources and dipoles on the cavity surface are moved to the newly computed cavity surface, where the boundary conditions are satisfied again. This iteration process is repeated until the results are converged. Characteristics of iteration and discretization of the present numerical method are much faster and more stable than the existing nonlinear theories. The theory shows good correlations with the existing theories and experimental results for the super-cavitating flow. In the region of small angles of attack, the present prediction shows and excellent comparison with the Geurst's linear theory. For the long cavity, the method recovers the trends of the Wu's nonlinear theory. In the intermediate regions of the short super-cavitation, the method compares very well with the experimental results of Parkin and also those of Silberman.