• 한국전자파학회논문지
• /
• 제22권10호
• /
• pp.1005-1011
• /
• 2011

본 논문에서는 일반적인 분산 매질의 전자기 과도 응답을 해석하기 위하여 헬름홀츠 방정식에 근거한 MODFDM(Marching-on-in-Degree Finite Difference Method) 기법을 제안한다. 라게르 함수의 특성을 이용하여 시간에 대한 미분항과 상승 적분(convolution integral)의 근사를 해석적으로 처리하였다. 본 기법의 기본적인 독창성은 전장과 전속 밀도, 유전율 등을 모두 라게르 함수로 전개한 다음, 갤러킨 시험 과정을 적용하여 시간 변수를 완전히 제거하였을 뿐만 아니라, 기존의 FDTD(Finite Difference Time-Domain) 방법과 달리 최종 계산식에 공간적인 유한 차분만을 적용하는데 있다. 일반적인 분산 매질의 해석에 적용 가능함을 보이기 위하여 대표적인 드바이, 드루드 및 로렌츠 분산 매질에 대한 전자기 과도 응답을 수치예로 보인다.

• 한국소음진동공학회논문집
• /
• 제25권12호
• /
• pp.866-873
• /
• 2015

The transverse vibration of the deploying beam from rigid hub was analyzed by simulation and experiment. The linear governing equation of the deploying beam was obtained using the Euler-Bernoulli beam theory. To discretize the governing equation, the Galerkin method was used. After transforming the governing equation into the weak form, the weak form was discretized. The discretized equation was expressed by the matrix-vector form, and then the Newmark method was applied to simulate. To consider the damping effect of the beam, we conducted the modal test with various beam length. The mass proportional damping was selected by the relation of the first and second damping ratio. The proportional damping coefficient was calculated using the acquired natural frequency and damping ratio through the modal test. The experiment was set up to measure the transverse vibration of the deploying beam. The fixed beam at the carriage of the linear actuator was moved by moving the carriage. The transverse vibration of the deploying beam was observed by the Eulerian description near the hub. The deploying or retraction motion of the beam had the constant velocity and the velocity profile with acceleration and deceleration. We compared the transverse vibration results by the simulation and experiment. The observed response by the Eulerian description were analyzed.

• Wind and Structures
• /
• 제15권6호
• /
• pp.481-494
• /
• 2012

The greenhouse type metal structures are increasingly used in modern construction of livestock farms because they are less laborious to construct and they provide a more favorable microclimate for the growth of animals compared to conventional livestock structures. A key stress factor for metal structures is the wind. The external pressure coefficient ($c_{pe}$) is used for the calculation of the wind effect on the structures. A high pressure coefficient value leads to an increase of the construction weight and subsequently to an increase in the construction cost. The EC1 in conjunction with EN 13031-1:2001, which is specialized for greenhouses, gives values for this coefficient. This value must satisfy two requirements: the safety of the structure and a reduced construction cost. In this paper, the Navier - Stokes and continuity equations are solved numerically with the finite element method (Galerkin Method) in order to simulate the two dimensional, incompressible, viscous air flow over the vaulted roofs of single span and twin-span with eaves livestock greenhouses' structures, with a height of 4.5 meters and with length of span of 9.6 and 14 m. The simulation was carried out in a wind tunnel. The numerical results of pressure coefficients, as well as, the distribution of them are presented and compared with data from Eurocodes for wind actions (EC1, EN 13031-1:2001). The results of the numerical experiment were close to the values given by the Eurocodes mainly on the leeward area of the roof while on the windward area a further segmentation is suggested.

• 한국전자파학회논문지
• /
• 제20권9호
• /
• pp.936-942
• /
• 2009

본 논문에서는 함체 케이스의 내부에 손실 유전체를 충진한 부공동을 설치하고 함체 케이스 내부로 공급된 전력과 반사 계수를 계산하여 함체의 공진 특성을 검토하고 있다. 이론 해석으로는 내부 전자파원의 전류 분포 및 부공동 개구면에서의 전계 분포에 관한 연립 적분방정식을 유도하고 Galerkin의 모멘트법으로 해석하였다. 이론 해석 결과, carbon을 함유한 발포 폴리스티렌을 손실 유전체로 사용하여 부공동의 크기와 carbon 함유량을 조절함으로써 함체 케이스의 공진 특성을 제어하여 함체 내부의 전자파 방사 크기를 저감시킬 수 있음을 보이고 있으며, 공급 전력의 실험치와도 비교하여 이론 해석의 타당성을 확인하고 있다.

• 대한기계학회논문집A
• /
• 제24권12호
• /
• pp.2989-2994
• /
• 2000

A mechanical face seal is a tribo-element intended to control leakage of working fluid at the interface of a rotating shaft and its housing. Leakage of working fluid decreases drastically as the clearance between mating seal faces gets smaller. But the very small clearance may result in an increased reduction of seal life because of high wear and heat generation. Therefore, in the design of mechanical face seals a compromise between low leakage and acceptable seal life is important, ant it present a difficult and practical design problem. A fluid film or sealing dam geometry of the seal clearance affects seal lubrication performance very much, and thereby is optimization is one of the main design consideration. in this study the Reynolds equation for the sealing dam of mechanical face seals is numerically analyzed, using the Galerkin finite element method, which is readily applied to various seal geometries, to give lubrication performances, such as opening force, restoring moment, leakage, and axial and angular stiffness coefficients. Then, to improve the seal performance an optimization is performed, considering various design variables simultaneously. For the tested case the optimization ha successfully resulted in the optimal design values of outer and inner seal radii, coning, seal clearance, and balance radius while satisfying all the operation subjected constraints and design variable side-constraints, and improvements of axial and angular stiffness coefficients by 16.8% and 2.4% respectively and reduction of leakage by 38.4% have been achieved.

• 설비공학논문집
• /
• 제6권3호
• /
• pp.302-314
• /
• 1994

The effects of the inclination of enclosure and partition on natural convective flow and heat transfer were investigated numerically. The enclosure was composed of the lower hot and the upper cold horizontal walls and the adiabatic vertical walls, and a partition was positioned perpendicularly at the mid-height of one vertical insulated wall. The governing equations are solved by using the finite element method with Galerkin method. The computations were performed with the variations of the partition length and Rayleigh number based on the temperature difference between two horizontal walls and the enclosure height with water(Pr=4.95). The effects of the inclination angle of enclosure and partition on the heat transfer within an enclosure were also studied. As the results, the increase of the inclination angle of enclosure rapidly raised the heat transfer rate, while the inclination angle for the maximum Nusselt number was retarded with the increase of the partition length and the decrease of the heat transfer rate became larger in proportion to the increase of the partition length. The Nusselt number obtained by the inclination of partition was smaller than that of the inclination of enclosure. However, the difference of the heat transfer rates was considerably decreased at the longer partition lengths and the trends for the variation of the average Nusselt number were more similar with that of the inclination of enclosure. The upward oriented partition increases the convective heat transfer distinctly in contrast to that of the inclination of enclosure as the partition length increases.

• 대한기계학회논문집B
• /
• 제26권12호
• /
• pp.1738-1746
• /
• 2002

In the FSI (Fluid-Structure Interaction) problems, two different governing equations are to be solved together. One is fur the fluid and the other for the structure. Furthermore, a kinematic constraint should be imposed along the boundary between the fluid and the structure. We use the combined formulation, which incorporates both the fluid and structure equations of motion into a single coupled variational equation so that it is not necessary to calculate the fluid force on the surface of structure explicitly when solving the equations of motion of the structure. A two-dimensional channel flow divided by a Bernoulli-Euler beam is considered and the dynamic response of the beam under the influence of channel flow is studied. The Navier-Stokes equations are solved using a P2P1 Galerkin finite element method with ALE (Arbitrary Lagrangian-Eulerian) algorithm. The internal structural damping effect is not considered in this study and numerical results are compared with a previous work fer steady case. In addition to the Reynolds number, two non-dimensional parameters, which govern this fluid-structure system, are proposed. It is found that the larger the dynamic viscosity and density of the fluid are, the larger the damping of the beam is. Also, the added mass is found to be linearly proportional to the density of the fluid.

• Korean Chemical Engineering Research
• /
• 제47권1호
• /
• pp.72-78
• /
• 2009

전기전도성 이방성 복합재료의 방전가공에 대하여 비정상상태 수식모델을 세우고 갤러킨의 유한요소법으로 해를 구하였다. 피삭재의 온도 분포와 분화구의 모양 및 공작물 제거 속도를 공정 매개변수에 관하여 구득하였다. 계산의 정확도와 효율을 위하여 앞선 연구에서 최적치로 선정된 $12{\times}12$ 요소의 비규칙 체눈을 사용하였다. 알루미나/티타늄 카바이드 복합재료의 물성을 재료의 물성으로 선정하였고 51.4 V의 전압과 7 A의 전류를 갖는 전력을 적용하였으며 제거 효율을 10%로 전열 이방성 계수를 2와 3으로 가정하였다. 불꽃이 일어나면서 피삭재는 즉시 녹기 시작하였고 열적 손상 영역이 형성되었다. 또한 시간이 흘러감에 따라서 분화구의 경계가 이동하는 것이 확인되었다. 반경 방향과 축 방향의 열전도도가 독립적으로 커지면 온도분포와 분화구의 모양이 각각 반경 방향과 축 방향으로 이동하였다. 공작물 제거 속도는 축 방향의 열전도도보다 반경 방향의 열전도도가 증가할 때 더욱 커지는 것으로 나타났다.

• 한국지구과학회지
• /
• 제38권2호
• /
• pp.105-114
• /
• 2017

Two dimensional finite element method with quadrilateral basis functions was applied to the spherical high order filter on the spherical surface limited area domain. The basis function consists of four shape functions which are defined on separate four grid boxes sharing the same gridpoint. With the basis functions, the first order derivative was expressed as an algebraic equation associated with nine point stencil. As the theory depicts, the convergence rate of the error for the spherical Laplacian operator was found to be fourth order, while it was the second order for the spherical Laplacian operator. The accuracy of the new high order filter was shown to be almost the same as those of Fourier finite element high order filter. The two-dimension finite element high order filter was incorporated in the weather research and forecasting (WRF) model as a hyper viscosity. The effect of the high order filter was compared with the built-in viscosity scheme of the WRF model. It was revealed that the high order filter performed better than the built in viscosity scheme did in providing a sharper cutoff of small scale disturbances without affecting the large scale field. Simulation of the tropical cyclone track and intensity with the high order filter showed a forecast performance comparable to the built in viscosity scheme. However, the predicted amount and spatial distribution of the rainfall for the simulation with the high order filter was closer to the observed values than the case of built in viscosity scheme.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제10권6호
• /
• pp.692-710
• /
• 2018

The Three-dimensional (3-D) dynamical behaviors of a fluid-conveying pipe subjected to vortex-induced vibration are investigated with different internal flow velocity ${\nu}$. The values of the internal flow velocity are considered in both subcritical and supercritical regimes. During the study, the 3-D nonlinear equations are discretized by the Galerkin method and solved by a fourth-order Runge-Kutta method. The results indicate that for a constant internal flow velocity ${\nu}$ in the subcritical regime, the peak Cross-flow (CF) amplitude increases firstly and then decrease accompanied by amplitude jumps with the increase of the external reduced velocity. While two response bands are observed in the In-line (IL) direction. For the dynamics in the lock-in condition, 3-D periodic, quasi-periodic and chaotic vibrations are observed. A variety of CF and IL responses can be detected for different modes with the increase of ${\nu}$. For the cases studied in the supercritical regime, the dynamics shows a great diversity with that in the subcritical regime. Various dynamical responses, which include 3-D periodic, quasi-periodic as well as chaotic motions, are found while both CF and IL responses are coupled while ${\nu}$ is beyond the critical value. Besides, the responses corresponding to different couples of ${\mu}_1$ and ${\mu}_2$ are obviously distinct from each other.