• 한국가시화정보학회지
• /
• 제5권2호
• /
• pp.20-25
• /
• 2007

The geometry of a commercial passenger airplane is realized based on a Boeing 747-400 model through the photographic scanning and reverse engineering. The each element consisting of the plane such as fuselage, wing, vertical fin, stabilizer and engines, is individually generated and then the whole body is assembled by the photomodeler. The maximum error in the realized airplane is about 1.4% comparing with the real one. The three-dimensional inviscid steady compressible governing equations are solved in the unstructured tetrahedron grid system, and in a finite volume method using STAR-CD when the airplane flies at the cruise condition. The pressure distribution on the surface and the wing-tip vortices are visualized, and in addition to the aerodynamics coefficients, lift and drag are estimated.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 춘계학술대회
• /
• pp.1957-1962
• /
• 2004

The objective of this study is to perform the numerical investigation of flow characteristics in static mixers. Simulations are carried out for mixers consisting of up to six Kenics and PPM elements placed end-to-end at an angle of $90^{\circ}$and for a range of Reynolds number($1{\leq}Re{\leq}100$). The pressure drop across a six-element Kenics mixer is computed and compared with the previous experimental correlations. The results are in good agreement with the previous correlations. The simulated flow field of Kenics mixer is extremely complex and contains regions of transverse flow that is dominated by the interaction of vortices produced by the mixer elements.

• International Journal of Aeronautical and Space Sciences
• /
• 제15권4호
• /
• pp.356-365
• /
• 2014

This paper is the second in a series and aims to build a high-fidelity mathematical model for a propeller-driven airplane using the propeller's aerodynamics and inertial models, as developed in the first paper. It focuses on aerodynamic models for the fuselage, the main wing, and the stabilizers under the influence of the wake trailed from the propeller. For this, application of the vortex lattice method is proposed to reflect the propeller's wake effect on those aerodynamic surfaces. By considering the maneuvering flight states and the flow field generated by the propeller wake, the induced velocity at any point on the aerodynamic surfaces can be computed for general flight conditions. Thus, strip theory is well suited to predict the distribution of air loads over wing components and the viscous flow effect can be duly considered using the 2D aerodynamic coefficients for the airfoils used in each wing. These approaches are implemented in building a high-fidelity mathematical model for a propeller-driven airplane. Flight dynamic analysis modules for the trim, linearization, and simulation analyses were developed using the proposed techniques. The flight test results for a series of maneuvering flights with a scaled model were used for comparison with those obtained using the flight dynamics analysis modules to validate the usefulness of the present approaches. The resulting good correlations between the two data sets demonstrate that the flight characteristics of the propeller-driven airplane can be analyzed effectively through the integrated framework with the propeller and airframe aerodynamic models proposed in this study.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.151-158
• /
• 2009

To implement the insects' flapping flight for developing flapping MAVs(micro air vehicles), the unsteady flow characteristics of the insects' forward flight is investigated. In this paper, two-dimensional FSI(Fluid-Structure Interaction) simulations are conducted to examine realistic flow features of insects' flapping flight and to examine the flexibility effects of the insect's wing. The unsteady incompressible Navier-Stokes equations with an artificial compressibility method are implemented as the fluid module while the dynamic finite element equations using a direct integration method are employed as the solid module. In order to exchange physical information to each module, the common refinement method is employed as the data transfer method. Also, a simple and efficient dynamic grid deformation technique based on Delaunay graph mapping is used to deform computational grids. Compared to the earlier researches of two-dimensional rigid wing simulations, key physical phenomena and flow patterns such as vortex pairing and vortex staying can still be observed. For example, lift is mainly generated during downstroke motion by high effective angle of attack caused by translation and lagging motion. A large amount of thrust is generated abruptly at the end of upstroke motion. However, the quantitative aspect of flow field is somewhat different. A flexible wing generates more thrust but less lift than a rigid wing. This is because the net force acting on wing surface is split into two directions due to structural flexibility. As a consequence, thrust and propulsive efficiency was enhanced considerably compared to a rigid wing. From these numerical simulations, it is seen that the wing flexibility yields a significant impact on aerodynamic characteristics.

• 대한기계학회논문집B
• /
• 제25권9호
• /
• pp.1157-1164
• /
• 2001

Centrifugal fans are widely used and the noise generated by these machines causes one of the most serious problems. In general, the centrifugal fan noise is often dominated by tones at BPF(blade passage frequency) and its higher harmonics. This is a consequence of the strong interaction between the flow discharged from the impeller and the cutoff in the casing. However, only a few research have been carried out on predicting the noise because of the difficulty in obtaining detailed information about the flow field and casing effects on noise radiation. The objective of this study is to understand the generation mechanism of sound and to develop a prediction method for the unsteady flow field and the acoustic pressure field of a centrifugal fan. We assume that the impeller rotates with a constant angular velocity and the flow field of the impeller is incompressible and inviscid. So, a discrete vortex method(DVM) is used to model the centrifugal fan and to calculate the flow field. The force of each element on the blade is calculated by the unsteady Bernoulli equation. Lowsons method is used to predict the acoustic source. In order to compare the experimental data, a centrifugal impeller and wedge introduced by Weidemann are used in the numerical calculation and the results are compared with the experimental data. Reasonable results are obtained not only for the peak frequencies but also for the amplitudes of the tonal.

• Smart Structures and Systems
• /
• 제25권1호
• /
• pp.81-96
• /
• 2020

Wind measurements were made on the Canton Tower at a height of 461 m above ground during the Typhoon Vincente, the wind-induced accelerations and displacements of the tower were recorded as well. Comparisons of measured wind parameters at upper level of atmospheric boundary layer with those adopted in wind tunnel testing were presented. The measured turbulence intensity can be smaller than the design value, indicating that the wind tunnel testing may underestimate the crosswind structural responses for certain lock-in velocity range of vortex shedding. Analyses of peak factors and power spectral density for acceleration response shows that the crosswind responses are a combination of gust-induced buffeting and vortex-induced vibrations in the certain range of wind directions. The identified modal frequencies and mode shapes from acceleration data are found to be in good agreement with existing experimental results and the prediction from the finite element model. The damping ratios increase with amplitude of vibration or equivalently wind velocity which may be attributed to aerodynamic damping. In addition, the natural frequencies determined from the measured displacement are very close to those determined from the acceleration data for the first two modes. Finally, the relation between displacement responses and wind speed/direction was investigated.

• 한국해안·해양공학회논문집
• /
• 제24권4호
• /
• pp.235-246
• /
• 2012

연안방재계획 수립을 위한 기초자료 생성의 일환으로 태풍에 의한 해일 범람모의를 수행하였다. 통상적인 침수-노출 처리가 불가한 직립안벽을 갖는 항만역의 특징이 반영되도록 가상의 내부제방 개념을 도입하여 해석하였다. 목포항에 대한 실제해역 적용에 앞서 단순해역에서 다수의 적용성 평가를 통해 내부제방의 폭은 1 m가 적정한 것으로 평가되었다. 해일범람으로 인한 침수면적과 침수심 산정에 민감하게 작용하는 육역의 침수-노출 처리기법에서의 최소침수심은 격자에 의존적이나 10 m 내외의 유한요소 상세격자에서는 0.01 m가 적정한 것으로 판단된다. 기존 태풍경로에 준거한 100년 빈도급 중심기압과 최대풍속을 가지며 비대칭 경도풍과 파랑효과가 고려된 가상태풍의 내습에 따른 목포해역의 침수예상역이 모의되었다. 침수면적은 기존 연구결과와 유사성을 나타내지만 항만권역 일부에서 확대 가능성이 있는 것으로 평가된다.

• 대한토목학회논문집
• /
• 제33권6호
• /
• pp.2267-2275
• /
• 2013

기존의 천수흐름 해석 상용모형에서는 내부 경계조건을 단순히 완전활동조건으로 가정하여 유체의 흐름을 해석함으로써 구조물 주위에서의 유속, 와도, 수위, 전단력의 분포, 항력 및 양력의 시간에 따른 변화 등을 올바르게 해석하지 못하였다. 본 연구에서는 구조물 주위에서의 흐름특성을 정확하게 예측할 수 있는 유한요소모형을 개발하고, 구조물에서의 경계조건을 활동길이를 이용한 부분활동조건으로 묘사하여 내부경계조건에 따른 원형 실린더 후면에서의 층류 흐름특성을 분석하였다. 종횡방향 유속 및 와도의 시간에 따른 변화, 후류길이, 활동길이에 따른 와류열의 변화와 질량보존율을 비교한 결과 완전활동조건을 부여한 경우에는 와류열이 전혀 형성되지 않고 완전한 층류흐름이 발생하였다. 부분활동조건을 입력한 경우 실린더 표면에서의 유속분포가 변화되어 전단력의 크기와 와도의 발생에 영향을 미치므로 무활조건을 부여한 경우에 비해 와류열의 발생 주기가 짧아졌다. 최대 질량보존 오차는 무활조건을 적용한 경우 0.73%로 나타났으며, 무활조건에 비해 부분활동조건을 부여한 경우의 오차율이 최대 0.21% 감소하였다.

• 한국자동차공학회논문집
• /
• 제6권2호
• /
• pp.12-20
• /
• 1998

Exter ballistics of a typical high-speed projectile is studied through a flow-visualization experiment and an unstructured grid Navier-Srokes computation. Experiment produced a schlieren photograph that adequately shows the characteristic features of this complex flow, namely two kinds of oblique cone shocks and turbulent wake developing into the downstream. A hybrid scheme of finite volume-element method is used to simulate the compressible Reynolds-Averaged Navier-Stok- es solution on unstructured grids. Osher's approximate Riemann solver is used to discretize the cinvection term. Higher-order spatial accuracy is obtained by MUSCL extension and van Albada ty- pe flux limiter is used to stabilize the numerical oscillation near the solution discontinuity. Accurate Gakerkin method is used to discretize the viscous term. Explict fourth-order Runge-Kutta method is used for the time-stepping, which simplifies the application of MUSCL extension. A two-layer k-$\varepsilon$ turbulence model is used to simulate the turbulent wakes accurately. Axisymmetric folw and two-dimensional flow with an angle of attack have been computed. Grid-dependency is also checked by carrying out the computation with doubled meshes. 2-D calculation shows that effect of angle of attack on the flow field is negligible. Axi-symmetric results of the computation agrees well with the flow visualization. Primary oblique shock is represented within 2-3 meshes in numerical results, and the varicose mode of the vortex shedding is clearly captured in the turbulent wake region.

• Korea-Australia Rheology Journal
• /
• 제15권3호
• /
• pp.131-150
• /
• 2003

A new numerical algorithm of finite element methods is presented to solve high Deborah number flow problems with geometric singularities. The steady inertialess planar 4 : 1 contraction flow is chosen for its test. As a viscoelastic constitutive equation, we have applied the globally stable (dissipative and Hadamard stable) Leonov model that can also properly accommodate important nonlinear viscoelastic phenomena. The streamline upwinding method with discrete elastic-viscous stress splitting is incorporated. New interpolation functions classified as rational interpolation, an alternative formalism to enhance numerical convergence at high Deborah number, are implemented not for the whole set of finite elements but for a few elements attached to the entrance comer, where stress singularity seems to exist. The rational interpolation scheme contains one arbitrary parameter b that controls the singular behavior of the rational functions, and its value is specified to yield the best stabilization effect. The new interpolation method raises the limit of Deborah number by 2∼5 times. Therefore on average, we can obtain convergent solution up to the Deborah number of 200 for which the comer vortex size reaches 1.6 times of the half width of the upstream reservoir. Examining spatial violation of the positive definiteness of the elastic strain tensor, we conjecture that the stabilization effect results from the peculiar behavior of rational functions identified as steep gradient on one domain boundary and linear slope on the other. Whereas the rational interpolation of both elastic strain and velocity distorts solutions significantly, it is shown that the variation of solutions incurred by rational interpolation only of the elastic strain is almost negligible. It is also verified that the rational interpolation deteriorates speed of convergence with respect to mesh refinement.