• Journal of Korean Society for Atmospheric Environment
• /
• v.16 no.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
• /
• v.4 no.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
• /
• v.17 no.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
• /
• v.52 no.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
• /
• v.24 no.4
• /
• pp.315-326
• /
• 2007

Recently introduced, several explicit solutions of relative motion between neighboring elliptic satellite orbits are reviewed. The performance of these solutions is compared with an analytic solution of the general linearized equation of motion. The inversion solution by the Hill-Clohessy-Wiltshire equations is used to produce the initial condition of numerical results. Despite the difference of the reference orbit, the relative motion with the relatively small eccentricity shows the similar results on elliptic case and circular case. In case of the 'chief' satellite with the relatively large eccentricity, HCW equation with the circular reference orbit has relatively larger error than other elliptic equation of motion does.

• Journal of the Korean Statistical Society
• /
• v.25 no.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
• /
• v.8 no.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.

• Journal of KSNVE
• /
• v.7 no.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.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1996.10a
• /
• 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
• /
• v.61 no.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.