• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• 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.

• Journal of electromagnetic engineering and science
• /
• 제17권2호
• /
• pp.91-97
• /
• 2017

An accurate and efficient evaluation for hypersingular integrals (HIs), strongly singular integrals (SSIs), and weakly singular integrals (WSIs) plays an essential role in the numerical solutions of 3D electromagnetic scattering problems. We derive analytic formulas for WSIs based on Stokes' theorem, which can be expressed in elementary functions. Several numerical examples are presented to validate these analytic formulas. Then, to show the feasibility of the proposed formulations for numerical methods, these formulations are used with the existing analytical expressions of HIs and SSIs to correct the near-field interaction in an iterative physical optics (IPO) scheme. Using IPO, the scattering caused by a dihedral reflector is analyzed and compared with the results of the method of moments and measurement data.

• 대한기계학회논문집A
• /
• 제20권11호
• /
• pp.3498-3506
• /
• 1996

A loading device to be used in fracture experiment is presented. It's loading angle can be adjusted from $-45^{\circ}$ to $105^{\circ}$ at intervals of $15^{\circ}$ for a CTS ( compact tension-shear) specimen, so that it is to be useful to measure mixed mode toughness. The equations to give the $K_ I$, $K_ I1$ and J-integral for the experiment are evluated though finite elemetn analysis in which the loading procedure is simulated and the behaviors of the specimen such as load-displacement curve are estimated. In the course of the evaluation the values $K_ I$, $K_ I1$ and J-integral calculated through recentrly released numerical methods are employed as the reference ones.

• Journal of applied mathematics & informatics
• /
• 제17권1_2_3호
• /
• pp.73-91
• /
• 2005

The application of Green's function in calculation of flow characteristics around submerged and floating bodies due to a regular wave is presented. It is assumed that the fluid is homogeneous, inviscid and incompressible, the flow is irrotational and all body motions are small. Two methods based on the boundary integral equation method (BIEM) are applied to solve associated problems. The first is a low order panel method with triangular flat patches and uniform distribution of velocity potential on each panel. The second method is a high order panel method in which the kernels of the integral equations are modified to make it nonsingular and amenable to solution by the Gaussian quadrature formula. The calculations are performed on a submerged sphere and some floating spheroids of different aspect ratios. The excellent level of agreement with the analytical solutions shows that the second method is more accurate and reliable.

• Structural Engineering and Mechanics
• /
• 제36권6호
• /
• pp.669-678
• /
• 2010

In general means, the stress concentration problem of elastic plate with a rectangular hole can be investigated by numerical methods, and only approximative results are derived. This paper deduces an analytical study of the stress concentration due to a rectangular hole in an elastic plate under bending loads. Base on classical elasticity theory and FEM applying the U-transformation technique, the uncoupled governing equations with 3-DOF are established, and the analytical displacement solutions of the finite element equations are derived in series form or double integral form. Therefore, the stress concentration factor can then be discussed easily and conveniently. For the plate subjected to unidirectional bending loads, the non-conforming plate bending element with four nodes and 12-DOF is taken as examples to demonstrate the application of the proposed method. The inner force distribution is obtained. The solutions are adequate for the condition when the hole is far away from the edges and the thin plate subjected to any transverse loadings.

• 대한설비공학회지:설비저널
• /
• 제3권3호
• /
• pp.185-190
• /
• 1974

Floor panel heating system is popular in Korea as dwelling house heating system. There are two methods for keeping floor surface warm. One method is delivering warm air under the floor such as Roman Hypocaust and Korean traditional Ondol. The other method is imbedding hot water pipes into the concrete floor slab. This paper gives basic equations for steady and transient solutions of heat conduction from hurried pipes. For steady-state solution, fin Efficiency Method and Sink and Source Method were introduced. Sink and Source Method is applied to transient state and basic solution is given in the form of Exponential Integral Function. Numerical solutions can be solved easily by digital computer from these equations.

• 지구물리
• /
• 제5권1호
• /
• pp.63-71
• /
• 2002

적분 방정식법은 매우 강력한 3차원 전자탐사 모델링 기법이다. 그러나 이 방법은 이상체내에의 전기장의 계산을 위하여 대형 선형 방정식의 해를 구해야 하므로 계산시간이 많이 소요된다는 단점이 있다. 특히 3차원 역산의 경우에는 이러한 적분방정식의 단점은 치명적이 될 수 밖에 없다. 이러한 단점을 극복하기 위하여 다양한 근사법이 개발되었지만, 이들 방법은 모두 Green 텐서의 적분을 수행해야 하며, 이 적분은 상당한 계산시간이 소요된다. 한편 전자탐사에서 Green 텐서는 거리가 증가함에 따라 급격히 감소하는 양상을 보이므로, 모델링을 위하여 이상체를 여러 개의 미소요소로 분할할 경우, 먼 거리에 있는 미소요소는 측정 미소요소에서의 전기장에 크게 영향을 미치지 못한다. 본 논문에서는 이점에 착안하여 먼 거리에 있는 미소요소에 의한 Green 텐서를 무시함으로써 획기적으로 계산시간을 단축할 수 있었다. 결과의 정확성 검토를 위하여 간단한 형태의 이상체에 대하여 본 방법을 사용하여 전기장을 계산한 결과 적분 방정식법, 확작 Born 근사 및 급수법의 결과와 거의 일치하는 정확한 결과를 얻을 수 있었으며, 계산시간을 단축할 수 있었다.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2002년도 춘계학술대회논문집
• /
• pp.1010-1015
• /
• 2002

Over the last few decades, noise has played a major role in the development of aircraft engines. The dominant noise is generated by the wake interactions of fan and downstream stator. Engine inlet and exhaust ducts are being fitted with liner materials that aid in damping fan related noise. In this paper, the radiation of duct internal noise from duct open ends with liners is studies via numerical methods. The linearized Euler's equations in generalized curvilinear coordinates are solved by the DRP scheme. The far field sound pressure levels are computed by the Kirchhoff integration method. Through comparison of sound directivity from bell-mouth duct with and without liners, it is shown that radiation from engine inlet is affected by liner effects or a soft wall boundary condition.

• 대한기계학회논문집
• /
• 제14권6호
• /
• pp.1408-1416
• /
• 1990

본 연구에서는 과도 연성 열탄성문제의 해를 구할때 사용되는 직접시간적분방 법과 Laplace 변환방법은 상호 장 단점을 가지고 있다. 각 방법들의 장단점은 서로 배타적이므로 서로의 장점을 살리는 수치방법이 필요하다. 그런데, 대부분의 과도 열탄성문제는 급격한 온도변화로 인한 물체의 변형에 관심이 있기때문에 이 형태의 문 제를 효율적으로 다루는 데 주안점을 두고 본 연구를 수행하였다. 유도된 유한요소 방정식은 결국 열탄성 지배 방정식 중 열전달방정식인 에너지보존식은 Gurtin의 범함 수로부터 유도된 원래의 형태를 사용하나 수치적 안정성(numerical stability)을 보장 하기 위하여 운동방정식은 시간에 대한 2차미분 형태로 수정하였다. 에너지보존식은 시간에 대한 합성적분(convolution)형태로 표현되므로 온도의 시간미분항이 소거되므 로 경계에서의 급격한 온도변화로 인한 수치 해석적 문제점은 간단히 해결된다. 그 러므로, 제안된 수치해법은 직접시간적분방법의 일종이나 결과식인 유한요소방정식은 기존의 문헌들과 상당한 차이가 있다. 과도 연성 열탄성해석을 위한 새로운 근사수 치해법의 장점을 이론적으로 설명하기보다 수치계산면에서의 안전성, 정확성 및 효율 성이 있음을 증명하기 위하여 이미 발표된 문헌들에서 다룬 예제를 선정하여 해석결과 를 비교하였다.

• Interaction and multiscale mechanics
• /
• 제4권2호
• /
• pp.85-105
• /
• 2011

The influence of surface elasticity and surface residual stress on the elastic field of an isotropic nanoscale elastic layer of finite thickness bonded to a rigid material base is considered by employing the Gurtin-Murdoch continuum theory of elastic material surfaces. The fundamental solutions corresponding to buried vertical and horizontal line loads are obtained by using Fourier integral transform techniques. Selected numerical results are presented for the cases of a finite elastic layer and a semi-infinite elastic medium to portray the influence of surface elasticity and residual surface stress on the bulk stress field. It is found that the bulk stress field depends significantly on both surface elastic constants and residual surface stress. The consideration of out-of-plane terms of the surface stress yields significantly different solutions compared to previous studies. The solutions presented in this study can be used to examine a variety of practical problems involving nanoscale/soft material systems and to develop boundary integral equations methods for such systems.