• 한국항공우주학회지
• /
• 제36권10호
• /
• pp.939-946
• /
• 2008

압축성 유동장 해석을 위해 개발된 전산유체역학 프로그램의 응용범위를 확장하여 비압축성 유동장을 해석하고자 할 경우 해의 수렴성에 문제가 발생할 수 있다. 이러한 수렴성 저하 문제는 지배방정식의 수학적 특성에 기인하며 예조건화 방법을 이용하면 지배방정식의 큰 수정 없이 해결할 수 있다. 본 논문은 예조건화 기법을 적용한 압축성 유동해석 프로그램을 개발하여 수렴성의 문제없이 저속의 비압축성 유동장 해석의 적용에 관한 내용을 다루고 있다. 개발된 프로그램의 정확성과 수렴특성을 확인하기 위해 정상상태의 비점성 및 층류 그리고 난류 문제에 대하여 해석하였으며 이를 실험치와 비압축성 계산치와 비교하였다

• Journal of Mechanical Science and Technology
• /
• 제20권12호
• /
• pp.2292-2314
• /
• 2006

Efficient numerical method is developed for the prediction of aerodynamic noise generation and propagation in low Mach number flows such as aeolian tone noise. The proposed numerical method is based on acoustic/viscous splitting techniques of which acoustic solvers use simplified linearised Euler equations, full linearised Euler equations and nonlinear perturbation equations as acoustic governing equations. All of acoustic equations are forced with immersed surface dipole model which is developed for the efficient computation of aerodynamic noise generation and propagation in low Mach number flows in which dipole source, originating from unsteady pressure fluctuation on a solid surface, is known to be more efficient than quadrupole sources. Multi-scale overset grid technique is also utilized to resolve the complex geometries. Initially, aeolian tone from single cylinder is considered to examine the effects that the immersed surface dipole models combined with the different acoustic governing equations have on the overall accuracy of the method. Then, the current numerical method is applied to the simulation of the aeolian tones from twin cylinders aligned perpendicularly to the mean flow and separated 3 diameters between their centers. In this configuration, symmetric vortices are shed from twin cylinders, which leads to the anti-phase of the lift dipoles and the in-phase of the drag dipoles. Due to these phase differences, the directivity of the fluctuating pressure from the lift dipoles shows the comparable magnitude with that from the drag dipoles at 10 diameters apart from the origin. However, the directivity at 100 diameters shows that the lift-dipole originated noise has larger magnitude than, but still comparable to, that of the drag-dipole one. Comparison of the numerical results with and without mean flow effects on the acoustic wave emphasizes the effects of the sheared background flows around the cylinders on the propagating acoustic waves, which is not generally considered by the classic acoustic analogy methods. Through the comparison of the results using the immersed surface dipole models with those using point sources, it is demonstrated that the current methods can allow for the complex interactions between the acoustic wave and the solid wall and the effects of the mean flow on the acoustic waves.

• 한국항공우주학회지
• /
• 제47권12호
• /
• pp.831-840
• /
• 2019

플러시 대기자료 측정장치는 비행체 표면에서 측정되는 압력 데이터를 이용하여 대기자료를 예측한다. FADS는 돌출된 프로브가 없으므로 고성능 항공기, 스텔스 비행체 및 극초음속 비행체에 적합하다. 본 논문에서는 구-원추 형상을 갖는 비행체에 대해서 아음속부터 초음속 비행까지 대기자료를 예측할 수 있는 FADS의 교정 절차와 계산 알고리즘을 제시한다. 표면 압력 데이터 측정을 위해 노즈부 표면에 5개 플러시 압력공들을 마련하였다. 유동각 예측과 압력 관련 변수의 예측을 분리하는 개념이며, 아음속 유동의 포텐셜 유동해와 극초음속 유동의 수정 뉴톤식을 결합한 압력모델을 사용한다. 교정 압력 데이터는 Euler 방정식을 푸는 전산유체역학 코드를 만들어서 마흐수 0.5 ~ 3.0의 범위에서 구축하였다. 비행 마흐 수 0.6~3.0, 받음각과 옆미끄럼각은 각각 -10° ~ +10°의 범위에서 여러 비행조건에 대해서 테스트를 수행하였다. 예측된 대기자료는 받음각, 옆미끄럼각, 마흐수, 자유류 정압이며 참고 데이터와 비교하여 정확도를 분석하였다.

• Smart Structures and Systems
• /
• 제25권2호
• /
• pp.219-228
• /
• 2020

This work presents a modified continuum model to explore and investigate static and vibration behaviors of perforated piezoelectric NEMS structure. The perforated nanostructure is modeled as a thin perforated nanobeam element with Euler-Bernoulli kinematic assumptions. A size scale effect is considered by included a nonlocal constitutive equation of Eringen in differential form. Modifications of geometrical parameters of perforated nanobeams are presented in simplified forms. To satisfy the Maxwell's equation, the distribution of electric potential for the piezoelectric nanobeam model is assumed to be varied as a combination of a cosine and linear functions. Hamilton's principle is exploited to develop mathematical governing equations. Modified numerical finite model is adopted to solve the equation of motion and equilibrium equation. The proposed model is validated with previous respectable work. Numerical investigations are presented to illustrate effects of the number of perforated holes, perforation size, nonlocal parameter, boundary conditions, and external electric voltage on the electro-mechanical behaviors of piezoelectric nanobeams.

• Advances in nano research
• /
• 제6권1호
• /
• pp.21-37
• /
• 2018

In the present article, wave dispersion behavior of a temperature-dependent functionally graded (FG) nanobeam undergoing rotation subjected to thermal loading is investigated according to nonlocal strain gradient theory, in which the stress numerates for both nonlocal stress field and the strain gradient stress field. The small size effects are taken into account by using the nonlocal strain gradient theory which contains two scale parameters. Mori-Tanaka distribution model is considered to express the gradually variation of material properties across the thickness. The governing equations are derived as a function of axial force due to centrifugal stiffening and displacements by applying Hamilton's principle according to Euler-Bernoulli beam theory. By applying an analytical solution, the dispersion relations of rotating FG nanobeam are obtained by solving an eigenvalue problem. Obviously, numerical results indicate that various parameters such as angular velocity, gradient index, temperature change, wave number and nonlocality parameter have significant influences on the wave characteristics of rotating FG nanobeams. Hence, the results of this research can provide useful information for the next generation studies and accurate deigns of nanomachines including nanoscale molecular bearings and nanogears, etc.

• Advances in nano research
• /
• 제6권1호
• /
• pp.57-80
• /
• 2018

In this study, the thermal loading effect on free vibration characteristics of carbon nanotubes (CNTs) with multiple cracks is studied. Various boundary conditions for nanotube are taken in to account. In order to take the small scale effect, the nonlocal elasticity of Eringen is employed in the framework of Euler-Bernoulli beam theory. This theory states that the stress at a reference point is a function of strains at all points in the continuum. A cracked nanotube is assumed to be consisted of two segments that are connected by a rotational spring which is located in the position of the cracked section. Hamilton's principle is used to achieve the governing equations. Influences of the nonlocal parameter, crack severity, temperature change and the number of cracks on the system frequencies are investigated. Also, it is found that at room or lower temperature the natural frequency for CNT decreases as the value of temperature change increases, while at temperature higher than room temperature the natural frequency of CNT increases as the value of temperature change increases. Various boundary conditions have been applied to the nanotube.

• International Journal of Fluid Machinery and Systems
• /
• 제8권3호
• /
• pp.169-182
• /
• 2015

When simulating the dynamic behaviour of a hydro power plant, it is essential to have a good representation of the turbine behaviour. The pressure transients in the system occurs because the flow changes, which the turbine defines. The flow through the turbine is a function of the pressure, the speed of rotation and the wicket gate opening and is, most often described in a performance diagram or Hill diagram. In the Hill diagram, the efficiency is drawn like contour lines, hence the name. A turbines Hill diagram is obtained by performance tests on scaled model in a laboratory. However, system dynamic simulations have to be performed in the early stage of a project, before the turbine manufacturer has been chosen and the Hill diagram is known. Therefore one have to rely on diagrams for a turbine with similar speed number. The Hill diagram is drawn through measured points, so for using the diagram in a simulation program, one have to iterate in the diagram based on curve fitting of the measured points. This paper describes an alternative method. By means of the Euler turbine equation, it is possible to set up two differential equations which represents the turbine performance with good enough accuracy for the dynamic simulations. The only input is the turbine's main geometry, the runner blade in- and outlet angle and the guide vane angle at best efficiency point of operation (BEP). In the paper, simulated turbine characteristics for a high head Francis turbine, and for a reversible pump turbine are compared with laboratory measured characteristics.

• 한국전산유체공학회지
• /
• 제6권4호
• /
• pp.26-34
• /
• 2001

The flow past a circular cylinder forced to vibrate transversely is numerically simulated by solving the two-dimensional Navier-Stokes equations modified by the vibration velocity of a circular cylinder at a Reynolds number of 164. The higher-order finite difference scheme is employed for the spatial discretization along with the second order Adams-Bashforth and the first order backward-Euler time integration. The calculated cylinder vibration frequency is between 0.60 and 1.30 times of the natural vortex-shedding frequency. The calculated oscillation amplitude extends to 25% of the cylinder diameter and in the case of the lock-in region it is 60%. It is made clear that the cylinder oscillation has influence on the wake pattern, the time histories of the drag and lift forces, power spectral density and phase diagrams, etc. It is found that these results include both the periodic (lock-in) and the quasi-periodic (non-lock-in) state. The vortex shedding frequency equals the driving frequency in the lock-in region but is independent in the non-lock-in region. The mean drag and the maximum lift coefficient increase with the increase of the forcing amplitude in the lock-in state. The lock-in boundaries are also established from the present direct numerical simulation.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2001년도 춘계 학술대회논문집
• /
• pp.181-188
• /
• 2001

The flow past a circular cylinder forced to vibrate transversely is numerically simulated by solving the two-dimensional Wavier-Stokes equations modified by the vibration velocity of a circular cylinder at a Reynolds number of 164. The higher-order finite difference scheme is employed for the spatial discretization along with the second order Adams-Bashforth and the first order backward-Euler time integration. The calculated cylinder vibration frequency is between 0.60 and 1.30 times of the natural vortex-shedding frequency. The calculated oscillation amplitude extends to $25\%$ of the cylinder diameter and in the case of the lock-in region it is $60\%$. It is made clear that the cylinder oscillation has influence on the wake pattern, the time histories of the drag and lift forces, power spectral density and phase diagrams, etc. It is found that these results include both the periodic (lock-in) and the quasi-periodic (non-lock-in) state. The vortex shedding frequency equals the driving frequency in the lock-in region but is independent in the non-lock-in region. The mean drag and the maximum lift coefficient increase with the increase of the forcing amplitude in the lock-in state. The lock-in boundaries are also established from the present direct numerical simulation.

• Structural Engineering and Mechanics
• /
• 제43권6호
• /
• pp.739-759
• /
• 2012

The tailoring optimization technique recently developed by the author for improving structural response and energy absorption of composites is extended to sandwich shells using a previously developed zig-zag shell model with hierarchic representation of displacements. The in-plane variation of the stiffness properties of plies and the through-the thickness variation of the core properties are determined solving the Euler-Lagrange equations of an extremal problem in which the strain energy due to out-of-plane strains and stresses is minimised, while that due to their in-plane counterparts is maximised. In this way, the energy stored by unwanted out-of-plane modes involving weak properties is transferred to acceptable in-plane modes. As shown by the numerical applications, the critical interlaminar stress concentrations at the interfaces with the core are consistently reduced without any bending stiffness loss and the strength to debonding of faces from the core is improved. The structural model was recently developed by the author to accurately describe strain energy and interlaminar stresses from the constitutive equations. It a priori fulfills the displacement and stress contact conditions at the interfaces, considers a second order expansion of Lame's coefficients and a hierarchic representation that adapts to the variation of solutions. Its functional d.o.f. are the traditional mid-plane displacements and the shear rotations, so refinement implies no increase of the number of functional d.o.f. Sandwich shells are represented as multilayered shells made of layers with different thickness and material properties, the core being treated as a thick intermediate layer.