• 한국대기환경학회지
• /
• 제16권4호
• /
• pp.339-350
• /
• 2000

The three-dimensional new corrdinate system over a single hill double hills and complex terrain with a single hill and a rectangular obstacle was generated using a body-fitted coordinate system. Control of the coordinate line distribution in the field was executed by generalizing the elliptic generating system to Poisson equation. ▽2ξ=P. The new coordinate system was well fitted to the surface boundary of single hill and double hills. But in the case of complex terrain with hill and rectangular obstacle there was smoothing tendency around the rectangular obstacle. In order to show the validity of the body-fitted coordinate system the heat diffusion equation was transformed and the temperature distribution was calculated over the various terrain. The results showed the temperature distribution was very symmetrical and stable around hills and obstacle. As a result the couple of a body-fitted coordinate system and the heat diffusion equation were executed successfully. Wind field over complex terrain with hill and rectangular obstacle which represent urban area was simulated stably in body-fitted coordinate system. The qualitative result show the enhancement of wind speed at the upwind direction of a hill and a rectangular obstacle and the recirculation zone at the downwind direction.

• International Journal of Aeronautical and Space Sciences
• /
• 제4권2호
• /
• pp.1-8
• /
• 2003

Satellite formation flying is the placing micro-satellites with the same mission into nearby orbits to form a cluster. Clohessy-Wiltshire equations are used to describe the relative motion and control strategies between satellites within a cluster, which are known as Hill's equations. Even though Hill's equations are powerful in determining initial conditions for the satellite formation flying, they can not accurately express the relative motion under J2 perturbation. Some methods have been developed for the determination of initial conditions to avoid limits of Hill's equation. This paper gives a new method of determining initial conditions using mean elements. For this research mean elements were transformed to osculating elements using Brouwer's theory and the orbit was propaeated with the consideration of J2-J8 to get a relative position. The results show that satellites within a cluster are maintained in the desired boundary for long period and the method is effective on a fuel saving for satellite formation flying.

• Journal of Mechanical Science and Technology
• /
• 제17권7호
• /
• pp.1054-1062
• /
• 2003

The numerical simulations of flow and pollutant particle dispersion are described for two-dimensional bell shaped hills with various aspect ratios. The Reynolds-averaged incompressible Navier-Stokes equations with low Reynolds number $\kappa$-$\varepsilon$ turbulent model are used to simulate the flowfield. The gradient diffusion equation is used to solve the pollutant dispersion field. The code was validated by comparison of velocity, turbulent kinetic energy, Reynolds shear stress, speed-up ratio, and ground level concentration with experimental and numerical data. Good agreement has been achieved and it has been found that the pollutant dispersion pattern and ground level concentration have been strongly influenced by the hill shape and aspect ratio, as well as the location and height of the source.

• Structural Engineering and Mechanics
• /
• 제52권3호
• /
• pp.525-540
• /
• 2014

The Lyapunov exponent and moment Lyapunov exponents of Hill's equation with frequency and damping coefficient fluctuated by correlated wideband random processes are studied in this paper. The method of stochastic averaging, both the first-order and the second-order, is applied. The averaged $It\hat{o}$ differential equation governing the pth norm is established and the pth moment Lyapunov exponents and Lyapunov exponent are then obtained. This method is applied to the study of the almost-sure and the moment stability of the stationary solution of the thin simply supported beam subjected to time-varying axial compressions and damping which are small intensity correlated stochastic excitations. The validity of the approximate results is checked by the numerical Monte Carlo simulation method for this stochastic system.

• Journal of Astronomy and Space Sciences
• /
• 제24권4호
• /
• pp.315-326
• /
• 2007

최근에 제안된, 이웃하는 타원 궤도의 상대운동에 대한 몇 가지 양함수형 해를 분석하였다. 이 해를 이용한 상대운동 결과를 일반 선형화 운동 방정식의 해석적 해와 비교했다. 수치계산 결과를 위한 초기 조건은 Hill-Clohessy-Wiltshire(HCW) 운동방정식에 의한 해의 역함수로 구했다. 기준 궤도의 차이에도 불구하고 상대적으로 작은 이심률의 궤도 일 경우에는 타원 상대운동 궤도와 원 상대운동 궤도의 결과는 근접했다. 주위성의 궤도가 상대적으로 큰 이심률을 가질 경우에는, 기본 궤도로 원 궤도를 이용하는 HCW 운동방정식은 다른 타원 상대운동 궤도 방정식의 해보다 상대적으로 큰 오차를 갖는다.

• Journal of the Korean Statistical Society
• /
• 제25권2호
• /
• pp.161-173
• /
• 1996

When we deal with a Hilbert space-valued Stochastic Differential Equation (SDE) (or Stochastic Partial Differential Equation (SPDE)), depending on some unknown parameters, the solution usually has a Fourier series expansion. In this situation we consider the maximum likelihood method for the statistical estimation problem and derive the asymptotic properties (consistency and normality) of the Maximum Likelihood Estimator (MLE).

• 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.

• 소음진동
• /
• 제7권3호
• /
• pp.489-498
• /
• 1997

In this paper, chaotic motions of a curved pipe conveying oscillatory flow are theoretically investigated. The nonliear partial differential equation of motion is derived by Newton's method. The transformed nonlinear ordinary differential equation is a type of Hill's equation, which has the external and parametric excitation with a same frequency. Bifurcation curves of chaotic motion of the piping systems are obtained by applying Melnikov's method. Numerical simulations are performed to demonstrate theoretical results and show the strange attractor of the chaotic motion.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 1996년도 추계학술대회논문집; 한국과학기술회관, 8 Nov. 1996
• /
• pp.288-294
• /
• 1996

In this paper, Chaotic motions of a curved pipe conveying oscillatory flow are theoretically investigated. The nonlinear partial differential equation of motion is derived by Newton's method. The transformed nonlinear ordinary differential equation is a type of Hill's equation, which have the parametric and external excitation. Bifurcation curves of chaotic motion of the piping systems are obtained by applying Melnikov's method. Poincare maps numerically demonstrate theoretical results and show transverse homoclinic orbit of the chaotic motion.

• Structural Engineering and Mechanics
• /
• 제61권4호
• /
• pp.497-510
• /
• 2017

A framed structure may be composed of two sub-structures, which are linked by a hinged joint. One sub-structure is the primary system and the other is the secondary system. The primary system, which is subjected to the periodic external load, can give rise to an auto-parametric resonance of the second system. Considering the geometric-stiffness effect produced by the axially internal force, the element equation of motion is derived by the extended Hamilton's principle. The element equations are then assembled into the global non-homogeneous Mathieu-Hill equations. The Newmark's method is introduced to solve the time-history responses of the non-homogeneous Mathieu-Hill equations. The energy-growth exponent/coefficient (EGE/EGC) and a finite-time Lyapunov exponent (FLE) are proposed for determining the auto-parametric instability boundaries of the structural system. The auto-parametric instabilities are numerically analyzed for the two frames. The influence of relative stiffness between the primary and secondary systems on the auto-parametric instability boundaries is investigated. A phenomenon of the "auto-parametric internal resonance" (the auto-parametric resonance of the second system induced by a normal resonance of the primary system) is predicted through the two numerical examples. The risk of auto-parametric internal resonance is emphasized. An auto-parametric resonance experiment of a ${\Gamma}$-shaped frame is conducted for verifying the theoretical predictions and present calculation method.