• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.12 no.6
• /
• pp.931-939
• /
• 2001

The slot coupling characteristics was analysed in a radial-line slot antenna for its design. The previously proposed waveguide model with a periodic boundary condition on its narrow walls and periodically arranged slots on its wide wall was used. The magnetic field integral equation and two dyadic Green\`s functions for respective regions was derived and the method of moments was used. To maximize the efficiency of numerical analysis and to extract singularities, two different kinds of basis functions, the entire domain basis function and the sub-domain one, are used. In addition, the Ewald sum technique for the rectangular waveguide and the Shanks transform for the half space were used to accelerate the computation of the slowly convergent potential Green\`s functions. Simulation results expressed the effects of the various design parameters on the slot coupling.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.6
• /
• pp.549-558
• /
• 2012

A topology optimization method for phononic crystals is developed for the design of sound barriers, using the level set approach. Given a frequency and an incident wave to the phononic crystals, an optimal shape of periodic inclusions is found by minimizing the norm of transmittance. In a sound field including scattering bodies, an acoustic wave can be refracted on the obstacle boundaries, which enables to control acoustic performance by taking the shape of inclusions as the design variables. In this research, we consider a layered structure which is composed of inclusions arranged periodically in horizontal direction while finite inclusions are distributed in vertical direction. Due to the periodicity of inclusions, a unit cell can be considered to analyze the wave propagation together with proper boundary conditions which are imposed on the left and right edges of the unit cell using the Bloch theorem. The boundary conditions for the lower and the upper boundaries of unit cell are described by impedance matrices, which represent the transmission of waves between the layered structure and the semi-infinite external media. A level set method is employed to describe the topology and the shape of inclusions. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set function, which facilitates to handle complicated topological shape changes. Through several numerical examples, the applicability of the proposed method is demonstrated.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.84-84
• /
• 2011

터널화재의 위험요소에 대한 해석을 위해서는 실제 상황을 재현한 실대형 실험이 가장 유용하겠지만 현실적으로 시간적, 공간적, 경제적인 제약이 따르기 때문에 CFD Modeling 기술의 이용 및 검증이 필요하고, 실제 상황에 가까운 현상의 재현을 위해서는 시뮬레이션의 정확도에 대한 향상이 필수적이다. 또한, CFD Modeling을 터널화재에 적용할 때 시뮬레이션의 질에 영향을 미칠 수 있는 요소들에 대한 결정이 선행되어야 한다. 우선, 터널의 기하학적 구조와 경계조건의 확립이 필요한데 필요한 정보를 얻기 위해서 어느정도 길이의 터널이 적절한지에 대해 생각할 필요가 있으며, 단면변화에 대한 결정을 통해 모델링을 수행하여야 한다. 모델링 작업이 선행된 후에 화재의 위치, 성장률, 최대 크기, 환기시스템 사항 등의 고려가 필요한데 이러한 조건들은 CFD Modeling의 결과에 직접적인 영향을 주기 때문에 충분한 사전조사가 이루어져야 하고, 각 사항들의 변수를 고려하여 다양한 화재시나리오의 도출이 가능할 수 있다. 마지막으로, 화재에서 발생된 열중 약 30%가 복사에 의해 주위 벽으로 전달될 수 있고 열은 연기가 가득찬 영역내에서 재분배될 수 있는데, 열전달 및 연기의 유동 등에 관한 자료를 기초로 화재현상에 대한 분석이 가능하다. 이러한 과정들을 통해 실제 상황에 가까운 설계화재 시나리오를 예측할 수 있다. 본 연구에서는 우리나라 최장대터널인 죽령터널에 대해 합리적인 가정을 통한 설계화재 시나리오를 기초로 화재시뮬레이션은 FDS(Fire Dynamics Simulator) 프로그램을 사용하여 화재 및 연기의 이동 양상을 분석하고, 피난시뮬레이션은 SIMULEX 프로그램을 사용하여 피난시간을 예측 함으로써 터널화재의 CFD Modeling에 의한 피난안전성을 검토하고자 한다.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1418-1423
• /
• 2004

Optimization of a heat exchanger with internally finned circular tubes has been performed for three-dimensional periodically fully developed turbulent flow and heat transfer. The design variables of fin number N, fin width ($d_1,d_2$) and fin height(H), are numerically optimized for the limiting conditions of $N=22{\sim}37$, $d_1=0.5{\sim}1.5$ mm, $d_2=0.5{\sim}1.5$ mm, $H=0.1{\sim}1.5$. Due to the periodic boundary conditions along main flow direction, the three layers of meshes are considered. The CFD and the mathematical optimization are coupled to optimize the heat exchanger. The flow and thermal fields are predicted using the finite volume method and the optimization is carried out by using the sequential quadratic programming (SQP) method which is widely used in the constrained nonlinear optimization problem.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.4 s.235
• /
• pp.460-468
• /
• 2005

Analyses of flow and heat transfer characteristics and shape optimization of internally finned circular tubes have been performed for three-dimensional periodically fully developed turbulent flow and heat transfer. CFD and mathematical optimization are coupled in order to optimize the shape of heat exchanger. The design variables such as fin widths $(d_{1},\;d_{2})$ and fin height (h) are numerically optimized by minimizing the pressure loss and maximizing the heat transfer rate for limiting conditions of $d_{1}=0.2\~1.5\;mm,\;d_{2}=0.2\~1.5\;mm,$ and $h=0.2\~1.5mm$. Due to the periodic boundary conditions along main flow direction, the three layers of meshes are considered. The flow and thermal fields are predicted using the finite volume method and the optimization is carried out by means of the sequential quadratic programming (SQP) method which is widely used in the constrained nonlinear optimization problem.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.3
• /
• pp.185-192
• /
• 1993

A numerical model has been used for the prediction of wave agitations in a harbor which are induced by the intrusion and transformation of incident waves. Based on linear wave theory a mild-slope equation has been used. A partial absorbing boundary condition has been used on solid boundary. Functional has been derived following Chen and Mei(l974)'s technique based on Hybrid Element Method which uses finite discretisation in the inner region and analytical solution of Helmholtz equation in the outer region. Final simultaneous equation has been solved using the Gaussian Elimination Method. Helmholtz natural period and second peak period of seiche in Donghae Harbor coincide very well with the results from numerical calculation. Computed amplification factors show good agreement, especially when the reflection coefficient on solid boundary is 0.99, with those of measurements.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.2
• /
• pp.197-206
• /
• 2003

A three-dimensional (3-D) method of analysis is presented for determining the free vibration frequencies and mode shapes of solid and hollow hemispherical shells of revolution of arbitrary wall thickness having arbitrary constraints on their boundaries. Unlike conventional shell theories, which are mathematically two-dimensional (2-D), the present method is based upon the 3-D dynamic equations of elasticity. Displacement components μ/sub Φ/, μ/sub z/, and μ/sub θ/ in the meridional, normal, and circumferential directions, respectively, are taken to be sinusoidal in time, periodic in θ, and algebraic polynomials in the Φ and z directions. Potential (strain) and kinetic energies of the hemispherical shells are formulated, and the Ritz method is used to solve the eigenvalue problem, thus yielding upper bound values of the frequencies obtained by minimizing the frequencies. As the degree of the polynomials is increased, frequencies converge to the exact values. Novel numerical results are presented for solid and hollow hemispheres with linear thickness variation. The effect on frequencies of a small axial conical hole is also discussed. Comparisons are made for the frequencies of completely free, thick hemispherical shells with uniform thickness from the present 3-D Ritz solutions and other 3-D finite element ones.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1995.10a
• /
• pp.221-226
• /
• 1995

터빈 익렬, 펌프 익차, 원형 냉각탑, 치차 등과 같이 동일한 형상이 원주 방향으로 반복되어 있는 순환 대칭 구조물의 진동특성을 유한 요소법을 사용하여 해석하는 경우에 전체구조를 모델링하는 대신에 구조물을 동일한 형상의 부분구조로 분할하여 부분구조 한개만을 모델링하고 분할된 경계에서 적절한 경계조건을 부과하여 진동해석을 수행함으로서 컴퓨터 기억용량을 절감시키고 계산시간을 단축할 수 있는 방법이 널리 사용되고 있다. Orris and Petyt[1]는 부분구조의 양쪽 분할 경계면, 즉 연결 경계상에 있는 절점변위의 상관관계를 복소파동전파식을 이용해서 구하여 부분구조의 감소된 복소강성행렬 및 질량행렬을 만들고 실수부와 허수부를 분리하여 유한요소해석을 수행하는 방법을 제안하였다. 유한요소 프로그램 ANSYS[2]에서는 이와 같은 방법을 사용하고 있다. Thomas[3]는 순회 정규모드를 이용하였고, 참고문헌[4]에서는 순회행렬을 이용하였다. 또한 유한요소 프로그램 MSC/NASTRAN[5]에서는 푸리에 급수를 이용하고 유한요소 절점의 위치 및 변위를 원통 좌표계를 표현하여 순환대칭구조물의 유한요소해석을 수행할 수 있도록 되어있다. 본 논문에서는 순환 대칭구조물의 형상의 주기성과 순환성을 고려하여 이산퓨리에 변환을 이용함으로써 순환대칭구조물의 유한요소진동해석을 체계적으로 저용량의 컴퓨터에서 신속하고 정확하게 수행할 수 있는 방법을 제안하고자 한다.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.31 no.3A
• /
• pp.321-327
• /
• 2006

In this paper, when a H-polarized plane wave is incident on the grating consisting of uniform resistive strips, electromagnetic scattering is analyzed using the moment of methods (MoM). The current density of each resistive strip on a grounded dielectric plane is fixed by zero at both edges. To satisfy the condition at both ends of each resistive strip, the induced surface current density is expanded in a series of cosine and sine functions. The scattered electromagnetic fields are expanded in a series of floquet mode functions. The boundary conditions are applied to obtain the unknown current coefficients. According to the variation of the involving parameters such as strip width and spacing and angle of the incident field, numerical simulations are performed by applying the Fourier-Galerkin moment method. The numerical results of the normalized reflected power for resistive strips case for zero and several resistivities are obtained.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.1
• /
• pp.47-54
• /
• 1993

For the numerical analysis of free oscillation characteristics in a harbor with general boundary and bottom topography, finite element method is applied. The governing Helmholtz equation is transformed into a generalized matrix eigenvalue problem using the standard finite element procedure. A computer code is developed for the numerical evaluation of natural frequencies and free oscillation modes. In the eigensolution process, a shifting strategy is devised for the treatment of numerical singularity. Scaling of coefficient matrix is also found to be effective for the alleviation of numerical ill-conditioning. For the test problems, firstly, analytical and numerical solutions are compared and validity of the code is obtained. Hence the method is successfully applicable for the real-world problems with general geometric boundaries and bottom topography.