• 한국항공우주학회지
• /
• 제30권4호
• /
• pp.18-27
• /
• 2002

비정렬 격자계에서 continuous adjoint 방정식을 사용하여 비점성 압축성 유동장에서의 이차원 날개꼴에 대한 공력 형상 최적설계를 수행하였다. 정확한 민감도를 구하기 위한 재구성 기법으로는 Laplacian averaging을 사용하였으며, extended stencil을 사용한 최소자승법으로 유동변수의 미분을 구하였다. 또한 유동장과 adjoint 방정식의 적절한 수렴 조건에 대해서 알아보았으며, 이를 이용하여 천음속 비점성 유동장에서 RAE2822 및 NACA0012 날개꼴에 대해서 조파항력 최소화 설계를 수행하였다. 설계된 날개꼴에서 다시 충격파를 복원하는 역설계문제를 적용하여 본 연구에서 사용된 최적설계 기법을 검증하였다.

• International Journal of Aeronautical and Space Sciences
• /
• 제3권1호
• /
• pp.74-85
• /
• 2002

The numerical experiment has been conducted to investigate the unsteady shock wave reflecting phenomena. The cell-vertex finite-volume, Roe's upwind flux difference splitting method with unstructured grid is implemented to solve unsteady Euler equations. The $4^{th}$-order Runge-Kutta method is applied for time integration. A linear reconstruction of the flux vector using the least-square method is applied to obtain the $2^{nd}$-order accuracy for the spatial derivatives. For a better resolution of the shock wave and slipline, the dynamic grid adaptation technique is adopted. The new concept of grid adaptation technique, which is much simpler than that of conventional techniques, is introduced for the current study. Three error indicators (divergence and curl of velocity, and gradient of density) are used for the grid adaptation procedure. Considering the quality of the solution and the numerical efficiency, the grid adaptation procedure was updated up to $2^{nd}$ level at every 20 time steps. For the convenience of comparison with other experimental and analytical results, the case of interaction between the straight incoming shock wave and a sharp wedge is simulated for various flow conditions. The numerical results show good agreement with other experimental and analytical results, in the shock wave reflecting structure, slipline, and the trajectory of the triple points. Some critical cases show disagreement with the analytical results, but these cases also have been proven to show hysteresis phenomena.

• 대한기계학회논문집B
• /
• 제40권3호
• /
• pp.181-189
• /
• 2016

수평 관군에 대하여 수직이고 상향으로 흐르는 2상 유동에 의한 감쇠비를 예측하기 위한 해석모델이 Sim에 의하여 개발되었다. 이 모델에서 평가된 2상 유동의 기공률, 압력손실 등의 유동변수는 기존의 실험식을 사용하여 계산하였다. 그러나 관군의 경우에 사용하기에는 약간의 개선이 요구된다. 따라서 관군 내에 흐르는 2상 유동의 유동 변수에 대한 더 많은 정보를 획득하기 위하여 실험적으로 연구할 필요가 있다. 실험은 공기 - 물의 2상 유동이 흐르는 정사각형 배열 관군에서의 압력계수와 2상 유동 마찰승수를 계산하기 위하여 수행되었다. 피치 직경 비는 1.35이었고, 실린더의 직경은 18 mm이다. 압력센서와 신호처리 장치를 이용하여 관군에서의 압력차를 측정하였다. 2상 유동 마찰승수와 오일러수를 계산하기 위하여 관군에 적용되는 비균질 유동의 기공률은 Feenstra 등의 실험식을 사용하여 계산하였다. 균질과 비균질 2상 유동의 마찰승수와 오일러의 수를 실험적으로 구하고 Sim의 어림적 모델에 근거한 이론적 해석 결과와 비교 분석하였다.

• Smart Structures and Systems
• /
• 제25권4호
• /
• pp.501-514
• /
• 2020

This article presents a comprehensive model to investigate a free vibration and resonance frequencies of nanostructure perforated beam element as nano-resonator. Nano-scale size dependency of regular square perforated beam is considered by using nonlocal differential form of Eringen constitutive equation. Equivalent mass, inertia, bending and shear rigidities of perforated beam structure are developed. Kinematic displacement assumptions of both Timoshenko and Euler-Bernoulli are assumed to consider thick and thin beams, respectively. So, this model considers the effect of shear on natural frequencies of perforated nanobeams. Equations of motion for local and nonlocal elastic beam are derived. After that, analytical solutions of frequency equations are deduced as function of nonlocal and perforation parameters. The proposed model is validated and verified with previous works. Parametric studies are performed to illustrate the influence of a long-range atomic interaction, hole perforation size, number of rows of holes and boundary conditions on fundamental frequencies of perforated nanobeams. The proposed model is supportive in designing and production of nanobeam resonator used in nanoelectromechanical systems NEMS.

• Journal of information and communication convergence engineering
• /
• 제7권4호
• /
• pp.539-544
• /
• 2009

An estimation of 3D velocity field including occluded parts without maxing tracer to the fluid had not only never been proposed but also impossible by the conventional computer vision algorithm. In this paper, we propose a new method of three dimensional optical flow of the fluid based on physical model, where some boundary conditions are given from a priori knowledge of the flow configuration. Optical flow is obtained by minimizing the mean square errors of a basic constraint and the matching error terms with visual data using Euler equations. Here, Navier-Stokes motion equations and the differences between occluded data and observable data are employed as the basic constrains. we verify the effectiveness of our proposed method by applying our algorithm to simulated data with partly artificially deleted and recovering the lacking data. Next, applying our method to the fluid of observable surface data and the knowledge of boundary conditions, we demonstrate that 3D optical flow are obtained by proposed algorithm.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 추계학술대회
• /
• pp.550-555
• /
• 2003

When a shock wave arrives at an open end of tube, an impulse wave is discharged from the tube exit and complicated flow is formed near tube exit. The flow field is influenced by the cross-sectional geometry of tube exit, such as circular, square, rectangular, trapezoid and etc. In the current study, three-dimensional propagation characteristics of impulse wave discharged from the tube exit with non-circular cross section are numerically investigated using a CFD method. Total variation diminishing (TVD) scheme is used to solve the three-dimensional, unsteady, compressible Euler equations. Computations are performed for the Mach numbers of the incident shock wave $M_{s}$ below 1.5. The results obtained show that the peak pressure of the impulse wave and propagation directivity depends on the cross-sectional geometry of tube exit and the Mach number of incident shock wave.

• Smart Structures and Systems
• /
• 제16권3호
• /
• pp.401-414
• /
• 2015

A developed hybrid method for crack identification of beams is presented. Based on the Euler-Bernouli beam theory and concepts of fracture mechanics, governing equation of the cracked beams is reformulated. Finite element (FE) method as a powerful numerical tool is used to discritize the equation in space domain. After transferring the equations from time domain to frequency domain, frequencies and mode shapes of the beam are obtained. Efficiency of the governed equation for free vibration analysis of the beams is shown by comparing the results with those available in literature and via ANSYS software. The used equation yields to move the influence of cracks from the stiffness matrix to the mass matrix. For crack identification measured data are produced by applying random error to the calculated frequencies and mode shapes. An objective function is prepared as root mean square error between measured and calculated data. To minimize the function, hybrid genetic algorithms (GAs) and particle swarm optimization (PSO) technique is introduced. Efficiency, Robustness, applicability and usefulness of the mixed optimization numerical tool in conjunction with the finite element method for identification of cracks locations and depths are shown via solving different examples.

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

The Pressure-based Unstructured-grid Navier-Stokes Solver PUNS-2/3D for incompressible steady and unsteady viscous flows has been developed. It is based on nonstaggered cell-centered finite volume method. Second-order upwind scheme with least-square reconstruction is used for convective fluxes. The SIMPLE method is implemented to couple the pressure and velocity fields. And the time derivatives in the momentum equations are discretised using a second-order Euler backward-differencing scheme. The discretised linear equations are solved by the preconditioned Biconjugate Gradient Stabilized method(Bi-CGSTAB). The developed solver is applied to validation problems using hybrid meshes.

• 제어로봇시스템학회:학술대회논문집
• /
• 제어로봇시스템학회 1987년도 한국자동제어학술회의논문집(한일합동학술편); 한국과학기술대학, 충남; 16-17 Oct. 1987
• /
• pp.932-938
• /
• 1987

The dynamic equations of robotic manipulators can be derived from either Newton-Euler equation or Lagrangian equation. Model parameters which appear in the resulting dynamic equation are the nonlinear functions of both the inertial parameters and the geometric parameters of robotic manipulators. The identification of the model parameters is important for advanced robot control. In the previous methods for the identification of the model parameters, the geometric parameters are required to be predetermined, or the robotic manipulators are required to follow some special motions. In this paper, we propose an approach to the identification of the model parameters, in which prior knowledge of the geometric parameters is not necessary. We show that the estimation equation for the model parameters can be formulated in an upper block triangular form. Utilizing the special structures, we obtain a simplified least-square estimation algorithm for the model parameter identification. To illustrate the practical use of our method, a 4DOF SCARA robot is examined.

• 대한조선학회지
• /
• 제14권1호
• /
• pp.1-10
• /
• 1977

The major objects of this report are to supplement data of natural frequencies of thin elastic rectangular plates to the available data, and to give an experimental verification for natural frequencies obtained by Rayleigh-Ritz method, the generation set of which are eigenfunctions of Euler beams. For the first object the following five models, for which data only for the fundamental mode or data only for square plates are available, are adopted; (1) two opposed edges are clamped and the other two opposed edges simply supported (C-C, S-S), (2) one edge is simply supported and the other three edges clamped (C-C, C-S), (3) one edge is free and the other three edges clamped (C-C, C-F), (4) two adjacent edges are clamped and the other two adjacent edges free (C-F, C-F). For the (C-C, S-S) model the frequency equation obtained with the mode shapes assumed as of a single trigonometric series is solved. And for the other four models Rayleigh-Ritz method taking eigenfunctions of Euler beams as the generating set is applied. The numerical examples are obtained up to the fourth, the fifth or the sixth order depending on the range of the aspect ratio (0.1-10.0). The number of terms in the generating set for Rayleigh-Ritz method is fifteen for all models. For the experiment three models made of 3.2mm thickness mild steel plate for general structure use were prepared in following size; $300mm{\times}600mm,\;600mm{\times}600mm\;and\;900mm{\times}600mm$. Their boundary conditions are made to fit (C-C, C-F) condition. From the experiment mechanical impedance curves based on the frequency response method were obtained together with phase relation diagrams. The experimental data are resulted in good conformity to calculated values.