• Communications of the Korean Mathematical Society
• /
• v.10 no.1
• /
• pp.39-48
• /
• 1995

The ADI method was introduced by Peaceman and Rachford [6] in 1955, to solve the discretized boundary value problems for elliptic and parabolic PDEs. The finite difference discretization of the model elliptic problem $$ (1) -\Delta u = f, \Omega = [0, 1] \times [0, 1] $$ $$ u = 0 on \delta \Omega $$ with 5-point centered finite difference discretization, with n +2 mesh-points in the x - direction and m + 2 points in the y direction, leads to the solution of a linear system of equations of the form $$ (2) Au = b $$ where A is a matrix of dimension $N = n \times m$. Without loss of generality and for the sake of simplicity, we will assume for the remainder of this paper that m = n, so that $N = n^2$.

• Korean Journal of Hydrosciences
• /
• v.3
• /
• pp.61-79
• /
• 1992

In oder to predict the dispersion of a thermal discharge with strong turbulent and buoyant effects, the development of a numerical model using turbulence model and its application are significantly increased. In this study, a three-dimensional steady-state model for the surface discharge heated water into quiescent water body is developed. For the model closure of turbulent terms the four-equation turbulence model is used. For economic mumerical simulation, the elliptic governing equations are transformed to the partially parabolic equations. In general, the simulated results by the present model agree well with the experimental results by Pande and Rajratnam (1977). The model characteristics are presented in comparison with the predicted results from the two-eqauation turbulence model by McGuirk and Rodi (1979).

• Structural Engineering and Mechanics
• /
• v.13 no.6
• /
• pp.713-730
• /
• 2002

This paper presents a finite element approach for determining the natural frequencies for planar inclined arches of various shapes vibrating in three-dimensional space. The profile of inclined arches, represented by undeformed centriodal axis of cross-section, is defined by the equation of plane curves expressed in the rectangular coordinates which are : circular, parabolic, sine, elliptic, and catenary shapes. In free vibration state, the arch is slightly displaced from its undeformed position. The linear relationship between curvature-torsion and axial strain is expressed in terms of the displacements in three-dimensional space. The finite element discretization along the span length is used rather than the total are length. Numerical results for arches of various shapes are given and they are in good agreement with those reported in literature. The natural frequency parameters and mode shapes are reported as functions of two nondimensional parameters: the span to cord length ratio (e) and the rise to cord length ratio (f).

• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• 1998.10a
• /
• pp.201-209
• /
• 1998

The propagation of water waves over irregular bottom bathymetry and around islands involves many process. In this study of numerical model is developed current in water of varying depth. The method used is splitting method and minimax approximation. This numerical method used is Crank-Nicolson scheme. This model is applied to Vincent shoal and compared with laboratory data. The results agreed well with laboratory data. The results agreed well with laboratory data. Current effect is considered in this study. So, the model is used for the estimation of rip current in the slowly varying topography.

• The Journal of the Acoustical Society of Korea
• /
• v.9 no.6
• /
• pp.62-73
• /
• 1990

한국 포항 연해에서 Range-dependent 환경을 수용할 수 있는 Parabolic Equation에 기초한 수치 해석법을 도입하여 수층의 음속구조, 음원 깊이 및 음원의 주파수 변화가 전달손실에 미치는 영향을 규 명하였다. 수층의 음속 구조에 따라 전달손실 분포는 크게 달라지는데 혼합층 구조와 성층화된 구조 및 수온전선 구조간의 전달손실 변화는 25~40 dB의 큰 차이를 보인다. 성층화된 음속 구조에서는 음원 깊 이 변화에 관계없이 거리에 따라 큰 전달손실을 나타낸다. 그러나 같음 음속구조하에서 주파수를 변화 시켰을 때에는 전달손실 차이가 15~30dB에 이른다. 결국 포항 여해에서 전달손실 분포는 특정 주파수에 대해서 음속구조의 시공간적 변화에 따라 크게 영향을 받는다. 또한 주어진 음속구조하에서 주파수 변 화에 따라서도 전달손실 분포가 크게 변하는데, 이는 Range-dependent 환경하에서도 최적 전달 주파수 가 결정될 수 있는 가능성을 제시한다.

• Smart Structures and Systems
• /
• v.19 no.4
• /
• pp.403-413
• /
• 2017

The present article encompasses a nonlinear finite element (FE) and genetic algorithm (GA) based optimal vibration energy harvesting from nonprismatic piezo-laminated cantilever beams. Three cases of cross section profiles (such as linear, parabolic and cubic) are modelled to analyse the geometric nonlinear effects on the output responses such as displacement, voltage, and power. The simultaneous effects of taper ratios (such as breadth and height taper) on the output power are also studied. The FE based nonlinear dynamic equation of motion has been solved by an implicit integration method (i.e., Newmark method in conjunction with the Newton-Raphson method). Besides this, a real coded GA based constrained optimization scheme has also been proposed to determine the best set of design variables for optimal harvesting of power within the safe limits of beam stress and PZT breakdown voltage.

• Structural Engineering and Mechanics
• /
• v.60 no.2
• /
• pp.271-299
• /
• 2016

In this article, based on the higher-order shear deformation plate theory, buckling analysis of a rectangular plate made of functionally graded piezoelectric materials and its effective parameters are investigated. Assuming the transverse distribution of electric potential to be a combination of a parabolic and a linear function of thickness coordinate, the equilibrium equations for the buckling analysis of an FGP rectangular plate are established. In addition to the Maxwell equation, all boundary conditions including the conditions on the top and bottom surfaces of the plate for closed and open circuited are satisfied. Considering double sine solution (Navier solution) for displacement field and electric potential, an analytical solution is obtained for full simply supported boundary conditions. The accurate buckling load of FGP plate is presented for both open and closed circuit conditions. It is found that the critical buckling load for open circuit is more than that of closed circuit in all loading conditions. Furthermore, it is observed that the influence of dielectric constants on the critical buckling load is more than those of others.

• Proceedings of the Acoustical Society of Korea Conference
• /
• 1996.06a
• /
• pp.64-68
• /
• 1996

수중음파전달 모델은 benchmark 시험을 통해 정확도, 적용범위, 계산시간 등의 성능을 평가받는다. 본 논문에서는 analytic 모델, 정상 모드 모델(normal mode model), 포물선 방정식 모델(parabolic equation model), 가우시안 빔 모델(Gaussian beam model), 스펙트럼 모델(spectral model) 등 거리의존 모델에 대해 benchmark 시험을 수행하였으며, benchmark 시험은 다음과 같은 세 가지 거리의존 해양환경으로 나누어 실시했다 : 1) 해수면과 해저면이 Dirichlet 경계조건인 이상 쐐기 문제(ideal wedge problem), 2) 해수면은 앞서 말한 Dirichlet 경계조건이나 해저면은 전달 손실이 있는 손실 통과 해저면 쐐기 문제(penetrable lossy bottom wedge problem), 3) 해수면은 앞서 말한 Dirichlet 경계조건이고 해저면은 Neumann 경계조건으로 서로 평행이면 음파전달 속도가 거리방향 의존인 경우, 경우 1은 anaytic 모델을 사용하고 경우 2는 정상 모드 모델, 포물선 방정식 모델, 스펙트럼 모델을 사용하였으며, 경우 3에 대해서는 가우시안 빔 모델과 포물선 방정식 모델을 사용하였다.

• Proceedings of the Acoustical Society of Korea Conference
• /
• 1996.06a
• /
• pp.33-36
• /
• 1996

해양에서 내부파 (internal wave)는 수직적으로 밀도가 크게 다른 두 층의 경계면에서 생성되는데, 일반적으로 경계면은 수온약층 (thermocline)에 해당된다. 본 연구에서는 간단한 단주기(하나의 주기만을 갖는)내부파를 가정하고 이 내부파에 의해서 야기될 수 있는 저주파 (100-400 HZ) 전파손실의 변동 정도를 모델링을 통해 추정하였다. 내부파의 파장은 1Km이고 진폭은 25m이며 수심 50-100m 사이에 존재하는 것으로 가정하였다. 전파손실 계산은 PE (parabolic equation) 기법을 도입한 모델을 이용하였다.수층의 두께는 2000m로 균일하고 퇴적층 및 기반암층의 두께는 각각 500m로 가정하였다. 모델링 결과 단주기 내부파의 1 파장이 진행하는 동안 고정된 수신 깊이별로 20km 거리에서 10 dB 이상의 큰 변동이 야기 될 수 있음을 확인하였다. 이는 내부파 내에서 큰 굴절을 변화로 인한 음파의 산란 때문으로 볼 수 있다. 실제 해양에서는 내부파가 단주기 파동 형태로 나타나기 보다는 여러 주기의 파동이 중첩되어 나타나므로 정확한 내부파 영향을 추정하기 위해서는 보다 실제에 가까운 내부파를 모델링할 필요가 있다.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1997.10a
• /
• pp.171-178
• /
• 1997

This paper describes the application of discretized continuum-type optimality criteria (DCOC) for the multispan partially prestressed concrete beams. The cost of construction as objective function which includes the costs of concrete, prestressing steel, non-prestressing steel and formwork is minimized. The design constraints include limits on the maximum deflection, flexural and shear strengths, in addition to ductility requirements, and upper and lower bounds on design variables as stipulated by the design code. Based on Kuhn-Tucker necessary conditions, the optimality criteria are explicitly derived in terms of the design variables-effective depth, eccentricity of prestressing steel and non-prestressing steel ratio. The prestressing profile is prescribed by parabolic functions. The self-weight of the structure is included in the equilibrium equation of the real system, as is the secondary effect resulting from the prestressing force. Two numerical examples of multispan PPC beams with rectangular cross-section are solved to show the applicability and efficiency fo the DCOC-based technique.