• International Journal of Aerospace System Engineering
• /
• v.7 no.2
• /
• pp.6-12
• /
• 2020

The precise prediction of future position of satellite depends on the accurate determination of orbit, which is also helpful in performing orbit maneuvers and trajectory correction maneuvers. For estimating the orbit of satellite many methods are being used. Some of the conventional methods are based on (i) Differential Correction (DC) (ii) Extended Kalman Filter (EKF). In this paper, Differential Evolution (DE) is used to determine the orbit. Orbit Determination using DC and EKF requires some initial guess of the state vector to initiate the algorithm, whereas DE does not require an initial guess since a wide range of bounds for the design unknown variables (orbital elements) is sufficient. This technique is uniformly valid for all orbits viz. circular, elliptic or hyperbolic. Simulated observations have been used to demonstrate the performance of the method. The observations are generated by including random noise. The simulation model that generates the observations includes the perturbation due to non-spherical earth up to second zonal harmonic term.

• Journal of the Korean Society of Safety
• /
• v.37 no.2
• /
• pp.1-9
• /
• 2022

Most factories deal with toxic or flammable chemicals in their industrial processes. These hazardous substances pose a risk of leakage due to accidents, such as fire and explosion. In the event of chemical release, massive casualties and property damage can result; hence, quantitative risk prediction and assessment are necessary. Several methods are available for evaluating chemical dispersion in the atmosphere, and most analyses are considered neutral in dispersion models and under far-field wind condition. The foregoing assumption renders a model valid only after a considerable time has elapsed from the moment chemicals are released or dispersed from a source. Hence, an initial dispersion model is required to assess risk quantitatively and predict the extent of damage because the most dangerous locations are those near a leak source. In this study, the dispersion model for initial consequence analysis was developed with three-dimensional unsteady advective diffusion equation. In this expression, instantaneous leakage is assumed as a puff, and wind velocity is considered as a coordinate transform in the solution. To minimize the buoyant force, ethane is used as leaked fuel, and two different diffusion coefficients are introduced. The calculated concentration field with a molecular diffusion coefficient shows a moving circular iso-line in the horizontal plane. The maximum concentration decreases as time progresses and distance increases. In the case of using a coefficient for turbulent diffusion, the dispersion along the wind velocity direction is enhanced, and an elliptic iso-contour line is found. The result yielded by a widely used commercial program, ALOHA, was compared with the end point of the lower explosion limit. In the future, we plan to build a more accurate and general initial risk assessment model by considering the turbulence diffusion and buoyancy effect on dispersion.

• Proceedings of the KSRS Conference
• /
• 2005.10a
• /
• pp.741-744
• /
• 2005

A numerical scattering model for artificial metal structure based on physical optics approximation is developed to identify the height of phase center, and the result is compared with interferometric SAR DEM. The interferometric SAR data were gathered by AIRSAR during PACRIM- II campaign on Jeju Island. Power transmission towers on piedmont pasture along the slopes of Mt. Halla look like elliptic risings in TOPSAR DEM. The heights of risings are quantitatively analyzed using a scattering model in the way of achieving the height of phase centers of power transmission towers. A numerical algorithm is developed on the basis of physical optics approximation. The structure of power transmission tower was decomposed into hundreds of rectangular metal plates, of which the scattering matrix is known in analytic form, and the calculated scattering fields were summed coherently. The effect of direct backscattering component, ground-scatterer component and scatterer-ground component are decomposed and computed individually for each rectangular metal plate. The $\Deltak-radar$ equivalent was used to calculate height of phase center of the scatterer. The heights of a selected power transmission tower and scattering algorithm results give existence and location of the transmission towers but not actual tower heights.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.24 no.2
• /
• pp.84-88
• /
• 2012

A numerical skill for effective treatment of inclined boundaries in a finite element method is introduced. A finite element method has been frequently used to simulate hydraulic phenomena in a coastal zone since it can be applied to irregular and complex geometry. In case elliptic partial equations are governing equations for a finite element model, however, there is a difficulty in treating boundary conditions properly for cases in which boundaries are vertically inclined. In this study, a method to treat such inclined boundaries using Bessel functions for a finite element method is introduced and compared with analytical solutions.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.10 s.253
• /
• pp.966-972
• /
• 2006

A numerical simulation is performed fur developing turbulent flow in a strongly curved 180 deg pipe and its downstream tangent by a new solution code(PowerCFD) which adopts an unstructured cell-centered method. The governing equations are discretized as the full elliptic from of the equations of motion. Three typical two-equation turbulence models of low-Reynolds-number form are used to approximate the turbulent stress field. Solutions fur both streamwise and circumferential velocity components are compared with the experimental data by Azzola et at.(1986). The ${\kappa}-{\omega}$ model by Wilcox(1988) is found to give better prediction performance than the other two. Predicted secondary velocities and streamwise velocity component contours at sequential longitudinal stations are also presented in order to enable a detailed description of the complete flow. It is also found that, in the bend both mean streamwise and secondary velocities never achieve a fully-developed state and the code is capable of producing very well the complex nature of steady flow in a strongly curved pipe.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.3
• /
• pp.716-724
• /
• 1990

Oblique impinging plane jets were investigated experimentally and numerically at Reynolds number 21000. The inclination angle was varied from 90.deg.(normal to the impinging plate) to 60.deg.. The distance H between the nozzle exit and the stagnation point on the impinging plate was fixed at H/D=8. The working fluid was air. The mean velocity components and turbulent quantities were measured by a hot-wire anemometer. And the static pressure distributions on the impinging plate were measured by a Pitot tube. In numerical computation, the governing partial differential equations of elliptic type were solved with conventional k-.epsilon. turbulence model. The measurements show that, after impingement, the jet half width alone the wall increases in both directions, and that similarity for each turbulent quantity such as Reynolds shear stress or turbulent kinetic energy is revealed in the wall jet region. The computed results show some deviation from experimental data in the impingement region, where streamline curvature is significant. However, the computed results agree qualitatively well with measurements.

• Journal of the Korea Institute of Information Security & Cryptology
• /
• v.21 no.3
• /
• pp.11-25
• /
• 2011

The Miller algorithm is employed in the typical pairing computation such as Weil, Tate and Ate for implementing ID based cryptosystem. By analyzing the Mrabet's attack that is one of fault attacks against the Miller algorithm, this paper presents au efficient fault attack in Affine coordinate system, it is the most basic coordinates for construction of elliptic curve. The proposed attack is the effective model of a count check fault attack, it is verified to work well by practical fault injection experiments and can omit the probabilistic analysis that is required in the previous counter fault model.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.2
• /
• pp.204-213
• /
• 2000

This study deals with jet impingement, which is extensively used in the process industries to achieve intense heating, cooling or drying rates and also widely employed as a test flow for turbulent models due to its complex flow configuration, on a flat plate by numerical methods. In this calculation, the finite volume method was employed to solve the Navier-stokes equation based on the non-orthogonal coordinate with non-staggered variable arrangement. To get a better understanding for the fluid flow and heat transfer characteristics of the turbulent jet impingements, $k-{\varepsilon}-{\overline{v^{'2}}}$ turbulent model was adapted and compared with the experimental data and the result of standard $k-{\varepsilon}$ turbulent model. Numerical calculations were carried out with various flow rates, nozzle to plate distances. In the case of the axisymmetric jet impingement on a flat plate, $k-{\varepsilon}-{\overline{v^{'2}}}$ turbulent model showed better agreement with the experimental data than the standard $k-{\varepsilon}$ turbulent model in the prediction of the mean velocity profiles, the turbulent velocity profiles. the turbulent shear stress and the heat transfer rate. The highest heat transfer rate can be obtained when the impingement occurs within the potential core..

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.1 no.1
• /
• pp.22-30
• /
• 1989

Wave propagation is changed by the effect of shoaling, current-depth refraction and shelter-ing etc. To solve these problems. numerous models have been developed. In the present study, a coordinate system is proposed based on the wave ray equation with the wave number equation including diffraction effects . The governing equation for the study was derived from the mild slope wave equation in non-steady state, including current effects (Kirby, 1986a) and trans-formed into an orthogonal curvilinear coordinate system on the basis of the wave ray equation. To obtain a numerical solution, an explicit finite difference scheme was used, and solved by the relaxation method. This model was tested for various cases: Firstly a submersed circular shoal and a constant unit depth. Secondly a submerged elliptic shoal on a slope, and finally a breakwater harbour with obliquely incident waves on a slope. The model was found to simulate the experimental results and other theoretical results in wave height and wave angle fairy well, and the applicability of the model around an arbitrary shaped coastal structure was also verified. To demonstrate the general usefullness of the present approach , the model is to be applied to a field situation with a complex bed topography.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.8 no.2
• /
• pp.146-150
• /
• 1996

The waves are most important dynamical factors for the analyses of structural stability and topographical changes on coastal engineering field. However, wind-generated waves are very irregular in shape and transformed through refraction, diffraction and shoaling when they propagate into shallow water where bottom topography and water depth vary significantly. Recently, Vincent and Briggs (1989) reported hydraulic model experiments for the transformation of monochromatic and directionally-spread irregular waves passing over a submerged elliptical mound. They concluded that for the case of combined refraction-diffraction of waves by a shoal, the propagation characteristics of the irregular and equivalent regular wave conditions can be vastly different. On the irregular wave transformation have been made theoretical and numerical studies for several years. Although theoretical and laboratory studies on wave transformation have progressed considerably, field measurement and comparison of numerical results with related theories are still necessary for the prediction of the phenomena in reality. In this study, field measurement of both incident and transformed waves in Youngil Bay were made using various kinds of equipments, and numerical computations were made on the transformed frequency spectra of large waves propagating over the shoal using Chae and Jeong's (1992) elliptic model. It is shown that this model results agree very well with field data, and thus the applicability of the model is now validated.