• Bulletin of the Korean Space Science Society
• /
• 2008.10a
• /
• pp.36.1-36.1
• /
• 2008

An Unscented Kalman Filter(UKF) for estimation of attitude and rate of a spacecraft using only magnetometer vector measurement is presented. The dynamics used in the filter is nonlinear rotational equation which is augmented by the quaternion kinematics to construct a process model. The filter is designed for low Earth orbit satellite, so the disturbance torques include gravity-gradient torque, magnetic disturbance torque, and aerodynamic drag. The magnetometer measurements are simulated based on time-varying position of the spacecraft. The filter has been tested not only in the standby mode but also in the detumbling mode. To stabilize the attitude, linear PD controller is applied and the actuator is assumed to be thruster. A Monte-Carlo simulation has been done to guarantee the stability of the filter performance to the various initial conditions. The UKF performance is compared to that of EKF and it reveals that UKF outperforms EKF.

• The Journal of the Petrological Society of Korea
• /
• v.23 no.4
• /
• pp.351-366
• /
• 2014

We assumed the situation where an eruption column had been formed by the explosive Plinian eruption from Mt. Baekdu and that the collapse of eruption column had caused pyroclastic density currents to occur. Based on this assumption, we simulated by using a Titan2D model. To find out about the range of the impacts of pyroclastic density currents by volcanic eruption scenarios, we studied the distance for the range of the impacts by VEIs. To compare the results by each volcanic eruption scenario, we set the location of the vent on the 8-direction flank of the outer rim and on the center of the caldera, the internal friction angle of the pyroclastic density currents as $35^{\circ}$, the bed friction angle as $16^{\circ}$. We set the pile height of column collapse and the vent diameter with various VEIs. We properly assumed the height of the column collapse, the diameter of the vent, the initial rates of the column collapse and the simulation period, based on the VEIs, gravity and the volume of the collapsed volcanic ash. According to the comparative analysis of the simulation results based on the increase of the eruption, the higher VEI by the increase of eruption products, the farther the pyroclastic density currents disperse. To the northwest from the vent on the northeast slope of the outer rim of the caldera, the impact range was 3.3 km, 4.6 km, 13.2 km, 24.0 km, 50.2 km, 83.4 km or more from VEI=2 to VEI=7, respectively. Once the database has been fully constructed, it can be used as a very important material in terms of disaster prevention and emergency management, which aim to minimize human and material damages in the vicinity of Mt. Baekdu when its eruption causes the pyroclastic density currents to occur.

• Economic and Environmental Geology
• /
• v.40 no.2 s.183
• /
• pp.153-170
• /
• 2007

Loading time and loading environment of the Dokdo seamounts were studied from flexure model and VGP(Virtual Geomagnetic Pole) determined by gravity and magnetic data. In spite of their similarity in size. a large difference about 50 mGal between gravity anomaly peaks of Dokdo and the Isabu Tablemount suggests different compensation degrees. Flexural modeling results show that the flexural rigidity(effective elastic thickness) of lithosphere for Dokdo is stronger(thicker) than that for the Isabu Tablemount. Also, it implies that the age of lithosphere at the time of loading of the Isabu Tablemount may be younger than that of Dokdo. Magnetic anomalies occur complicated over the Dokdo seamounts. Paleomagnetism was studied from VGP estimated by the least square and the seminorm magnetization methods with 1500 m upward continued magnetic anomalies. Age dating of Dokdo from previous study, flexural modeling, VGP, and geomagnetic polarity time scale suggest that after the cease of spreading in the Ulleung Basin, the Isabu Tablemount was formed first in normal polarity interval and followed by Dokdo. Also, they indicate that the fist large eruption of Dokdo was in normal polarity interval and the second large eruption in reversed polarity interval. The Simheungtaek Tablemount was formed in normal polarity interval between the formations of the Isabu Tablemount and Dokdo. These loading times for the Dokdo seamounts show a good coherence with the compressive stress period after the end of the opening of the East Sea. The Dokdo seamounts probably was caused by volcanism associated with the compressive stress.

• Journal of Astronomy and Space Sciences
• /
• v.22 no.4
• /
• pp.463-472
• /
• 2005

In this paper, the dynamic model development for interplanetary navigation has been discussed. The Cowell method for special perturbation theories was employed to develop an interplanetary trajectory propagator including the perturbations due to geopotential, the Earth's dynamic polar motion, the gravity of the Sun, the Moon and the other planets in the solar system, the relativistic effect of the Sun, solar radiation pressure, and atmospheric drag. The equations of motion in dynamic model were numerically integrated using Adams-Cowell 11th order predictor-corrector method. To compare the influences of each perturbation, trajectory propagation was performed using initial transfer orbit elements of the Mars Express mission launched in 2003, because it can be the criterion to choose proper perturbation models for navigation upon required accuracy. To investigate the performance of dynamic model developed, it was tested whether the spacecraft can reach the Mars. The interplanetary navigation tool developed in this study demonstrated the spacecraft entering the Mars SOI(Sphere of Influence) and its velocity .elative to the Mars was less than the escape velocity of the Mars, hence, the spacecraft can arrive at the target planet. The obtained results were also verified by using the AGI Satellite Tool Kit. It is concluded that the developed program is suitable for supporting interplanetary spacecraft mission for a future Korean Mars mission.

• Journal of Astronomy and Space Sciences
• /
• v.26 no.1
• /
• pp.31-46
• /
• 2009

An Unscented Kalman Filter (UKF) for estimation of the attitude and rate of a spacecraft using only magnetometer vector measurement is developed. The attitude dynamics used in the estimation is the nonlinear Euler's rotational equation which is augmented with the quaternion kinematics to construct a process model. The filter is designed for small satellite in low Earth orbit, so the disturbance torques include gravity-gradient torque, magnetic disturbance torque, and aerodynamic drag torque. The magnetometer measurements are simulated based on time-varying position of the spacecraft. The filter has been tested not only in the standby mode but also in the detumbling mode. Two types of actuators have been modeled and applied in the simulation. The PD controller is used for the two types of actuators (reaction wheels and thrusters) to detumble the spacecraft. The estimation error converged to within 5 deg for attitude and 0.1 deg/s for rate respectively when the two types of actuators were used. A joint state parameter estimation has been tested and the effect of the process noise covariance on the parameter estimation has been indicated. Also, Monte-Carlo simulations have been performed to test the capability of the filter to converge with the initial conditions sampled from a uniform distribution. Finally, the UKF performance has been compared to that of the EKF and it demonstrates that UKF slightly outperforms EKF. The developed algorithm can be applied to any type of small satellites that are actuated by magnetic torquers, reaction wheels or thrusters with a capability of magnetometer vector measurements for attitude and rate estimation.

• Aerospace Engineering and Technology
• /
• v.12 no.1
• /
• pp.87-93
• /
• 2013

The conceptual design lunar lander demonstrator has been developed to use as a test bed for advanced spacecraft technologies and to test a prototype planetary lander capable of vertical takeoff and landing. Size of the lunar lander demonstrator is the same as that of lunar lander conceptually designed, however, the weight of lunar lander demonstrator is designed in 1/6 scale in consideration of gravity difference between moon and earth. The thruster clustering and virtual flight test were performed in the demonstrator fixed on the ground. The demonstrator ground test has been conducted for two months in the test site for the solid motor combustion of the Goheung Flight Center. The purposes of ground test of demonstrator are to demonstrate and verify essential electronics, propulsion system, control algorithm, embedded software, structure and system operation technologies before developing the flight model lander. This paper is described about the virtual flight test including test configuration, test aims and test facilities

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.330-339
• /
• 2008

Ground-based interceptors(GBI) comprise a major element of the strategic defense against hostile targets like Intercontinental Ballistic Missiles(ICBM) and reentry vehicles(RV) dispersed from them. An optimum design of the subsystems is required to increase the performance and reliability of these GBI. Propulsion subsystem design and optimization is the motivation for this effort. This paper describes an effort in which an entire GBI missile system, including a multi-stage solid rocket booster, is considered simultaneously in a Genetic Algorithm(GA) performance optimization process. Single goal, constrained optimization is performed. For specified payload and miss distance, time of flight, the most important component in the optimization process is the booster, for its takeoff weight, time of flight, or a combination of the two. The GBI is assumed to be a multistage missile that uses target location data provided by two ground based RF radar sensors and two low earth orbit(LEO) IR sensors. 3Dimensional model is developed for a multistage target with a boost phase acceleration profile that depends on total mass, propellant mass and the specific impulse in the gravity field. The monostatic radar cross section (RCS) data of a three stage ICBM is used. For preliminary design, GBI is assumed to have a fixed initial position from the target launch point and zero launch delay. GBI carries the Kill Vehicle(KV) to an optimal position in space to allow it to complete the intercept. The objective is to design and optimize the propulsion system for the GBI that will fulfill mission requirements and objectives. The KV weight and volume requirements are specified in the problem definition before the optimization is computed. We have considered only continuous design variables, while considering discrete variables as input. Though the number of stages should also be one of the design variables, however, in this paper it is fixed as three. The elite solution from GA is passed on to(Sequential Quadratic Programming) SQP as near optimal guess. The SQP then performs local convergence to identify the minimum mass of the GBI. The performance of the three staged GBI is validated using a ballistic missile intercept scenario modeled in Matlab/SIMULINK.

• Proceedings of the Korean Geotechical Society Conference
• /
• 1999.10a
• /
• pp.449-456
• /
• 1999

The concrete wall is the most useful of retaining structure which can obtain the engineering stability, but has problems that is not friendly with nature environment in a fine view, such as poor rear drainage, and shrinkage crack by temperature difference, etc. Because of this problems, the research for a segmental crib retaining wall has been performed. A segmental crib retaining wall is quickly and easily erected because is possible to be erected as the individual members, and is not sensitive to differential settlement and earthquakes. Also, it shows effective drainage and has a friendly advantage with nature environment because of being able to be planted with vines and shrubs in retaining walls The design of crib retaining walls has traditionally been based on classical soil mechanics theories. These theories, originally derived by Rankine(1857) and Coulomb(1776), assume that the wall acts as a rigid body. This assumption results in failure being predicted by either monolithic overturning or base sliding mechanisms. However, the wall consists of individual members which have been created a three dimensional grid. This grid confines an fill mass which becomes part of the wall. The filled wall resists the earth pressure with the same mechanism of classical gravity walls. Because of the flexibility of the individual segment, it allows relative movement between the individual members within the wall. The three dimensional flexible grid leads to stress redistribution when the wall is subjected to external or fill loads. Due to the flexibility and the stress redistribution, the failure of segmental crib wall consists of not only overturing and base sliding but the local deformation and the failure between the segmental members. It has been researched in the field that due to this flexibility and load redistribution, serviceability failure of segmental crib walls is unlikely to be due to overturning or base sliding. Therefore, in this study, the relative displacement appearance of retaining wall due to variation of inclination is measured to examine this behavior characteristics. Also, the behavior characteristics of retaining walls by surcharge load, and location of acting point of retaining wall rear, and the displacement characteristics and deflections are estimated about the existence and nonexistence of Rear Stretcher performing an role in transmitting earth pressure of Header and Stretcher organizing retaining walls. This research focuses on the characteristics due to the behavior of retaining walls. This research focuses on the characteristics due to the behavior of retaining walls.

• Ocean and Polar Research
• /
• v.44 no.4
• /
• pp.269-285
• /
• 2022

The interaction between ocean and ice shelf is a critical physical process in relation to water mass transformations and ice shelf melting/freezing at the ocean-ice interface. However, it remains challenging to thoroughly understand the process due to a lack of observational data with respect to ice shelf cavities. This is the first study to simulate the variability and circulation of water mass both overlying the continental shelf and underneath an ice shelf and an ice tongue in the Terra Nova Bay (TNB), East Antarctica. To explore the properties of water mass and circulation patterns in the TNB and the corresponding effects on sub ice shelf basal melting, we explicitly incorporate the dynamic-thermodynamic processes acting on the ice shelf in the Regional Ocean Modeling System. The simulated water mass formation and circulation in the TNB region agree well with previous studies. The model results show that the TNB circulation is dominated by the geostrophic currents driven by lateral density gradients induced by the releasing of brine or freshwater at the polynya of the TNB. Meanwhile, the circulation dynamics in the cavity under the Nansen Ice shelf (NIS) are different from those in the TNB. The gravity-driven bottom current induced by High Salinity Shelf Water (HSSW) formed at the TNB polynya flows towards the grounding line, and the buoyance-driven flow associated with glacial meltwater generated by the HSSW emerges from the cavity along the ice base. Both current systems compose the thermohaline overturning circulation in the NIS cavity. This study estimates the NIS basal melting rate to be 0.98 m/a, which is comparable to the previously observed melt rate. However, the melting rate shows a significant variation in space and time.

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.99-108
• /
• 2008

To extend our detailed knowledge for the Hwasan caldera, we carried out magnetotelluric (MT) survey, which is pretty sensitive to electrical property variation in both horizontal and vertical direction of subsurface, across the Hwasan caldera with the direction of EW. The 2-D inversion results of observed MT data lead to following conclusions. Firstly, the depth of the basin basement inferred by the MT inversion results matches well with that suggested by previous potential studies, but the basement resistivity seems fairly low when compared to that of general case. This feature might be related with the large-scaled, highly conductive layer beneath the Euisung Sub-basin suggested by the previous MT study. Secondly, the high resistivity zones reaching to 4000 $\Omega{\cdot}m$ are imaged around two external ring fault boundaries. These zones are thought of as the response of the rhyolitic dykes intruding along the ring fault, and in the previous gravity data correspond to relatively high density anomalies. Thirdly, low resistivity zone reaching to 200 $\Omega{\cdot}m$ is detected around a depth of 1km beneath the central part of the caldera, which has not been yet reported in korean geophysical literatures. If we take account of the evolution model of the Hwasan caldera, this zone is regarded as the past sedimentary layer that subsided during the period of forming external ring fault system. In addition, the relatively low density anomaly observed in the central part of the caldera may be attributed to this sedimentary layer.