• Journal of Korea Water Resources Association
• /
• v.50 no.7
• /
• pp.455-465
• /
• 2017

Understanding of the two-dimensional velocity field is crucial in terms of analyzing various hydrodynamic and fluvial processes in the riverine environments. Until recently, many numerical models have played major roles of providing such velocity field instead of in-situ flow measurements, because there were limitations in instruments and methodologies suitable for efficiently measuring in the broad range of river reaches. In the last decades, however, the advent of modernized instrumentations started to revolutionize the flow measurements. Among others, acoustic Doppler current profilers (ADCPs) became very promising especially for accurately assessing streamflow discharge, and they are also able to provide the detailed velocity field very efficiently. Thus it became possible to capture the velocity field only with field observations. Since most of ADCPs measurements have been mostly conducted in the cross-sectional lines despite their capabilities, it is still required to apply appropriate interpolation methods to obtain dense velocity field as likely as results from numerical simulations. However, anisotropic nature of the meandering river channel could have brought in the difficulties for applying simple spatial interpolation methods for handling dynamic flow velocity vector, since the flow direction continuously changes over the curvature of the channel shape. Without considering anisotropic characteristics in terms of the meandering, therefore, conventional interpolation methods such as IDW and Kriging possibly lead to erroneous results, when they dealt with velocity vectors in the meandering channel. Based on the consecutive ADCP cross-sectional measurements in the meandering river channel. For this purpose, the geographic coordinate with the measured ADCP velocity was converted from the conventional Cartesian coordinate (x, y) to a curvilinear coordinate (s, n). The results from application of A-VIM showed significant improvement in accuracy as much as 41.5% in RMSE.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.2
• /
• pp.135-146
• /
• 2020

To define an interface in a conventional level-set method, an analytical function must be revealed for an interfacial geometry. However, it is not always possible to define a functional form of level sets when interfaces become complex in a Cartesian coordinate system. To overcome this difficulty, we have developed a new level-set formalism that discriminates the interior from the exterior of a CAD modeled interface by parameterizing the stereolithography (STL) file format. The work outlined here confirms the accuracy and scalability of the hydrodynamic reactive solver that utilizes the new level set concept through a series of tests. In particular, the complex interaction between shock and geometrical confinements towards deflagration-to-detonation transition is numerically investigated. Also, a process of flame spreading and damages caused by point source detonation in a real-sized plant facility have been simulated to confirm the validity of the new method built for reactive hydrodynamic simulation in any complex three-dimensional geometries.

• Journal of Astronomy and Space Sciences
• /
• v.20 no.1
• /
• pp.21-28
• /
• 2003

In case that we independently obtain orbital informations about the low earth satellites of foreign countries using radar systems, we develop the orbit determination algorithm for this purpose using a SGP4 model with an analytical orbit model and the extended Kalman filter with a real-time processing method. When the state vector is Keplerian orbital elements, singularity problems happen to compute partial derivative with respect to inclination and eccentricity orbit elements. To cope with this problem, we set state vector osculating to mean equinox and true equator cartesian elements with coordinate transformation. The state transition matrix and the covariance matrix are numerically computed using a SGP4 model. Observational measurements are the type of azimuth, elevation and range, filter process to each measurement in a lump. After analyzing performance of the developed orbit determination algorithm using TOPEX/POSEIDON POE(precision 0.bit Ephemeris), its position error has about 1 km. To be similar to performance of NORAD system that has up to 3km position accuracy during 7 days need to radar system performance that have accuracy within 0.1 degree for azimuth and elevation and 50m for range.

• Economic and Environmental Geology
• /
• v.53 no.3
• /
• pp.271-285
• /
• 2020

In this study, a numerical experiment related to the stress-strain analysis was performed on 3-D discrete fracture network(DFN) systems based on the distinct element method to evaluate the effect of joint geometry on deformability of jointed rock masses. Using one or two joint sets with deterministic orientation, a total of 12 3-D DFN blocks having 10m cube domain were generated with different joint density and size distribution. Directional deformation modulus of the DFN cube blocks were estimated along the axis directions of 3-D cartesian coordinate. In addition, deviatoric stress directions were chosen at every 30° of trend and plunge in 3-D for some DFN blocks to examine the variability of directional deformation modulus with respect to joint geometry. The directional deformation modulus of the DFN block were found to reduce with the increase of joint size distribution. The increase in joint density was less likely to have a significant effect on directional deformation modulus of the DFN block in case of the effect of rock bridges was relatively large because of short joint size distribution. It, however, was evaluated that the longer the joint size, the increase in the joint density had a more significant effect on the anisotropic deformation modulus of the DFN block. The variation of the anisotropic deformation modulus according to the variations in joint density and size distribution was highly dependent on the number of joint sets and their orientation in the DFN block. Finally, this study addressed a numerical procedure for stress-strain analysis of jointed rock masses considering joint geometry and discussed a methodology for practical application at the field scale.

• Journal of Korean Society of Steel Construction
• /
• v.10 no.1 s.34
• /
• pp.47-61
• /
• 1998

The behavior of horizontally curved I-girder bridges is complex because the flexural and torsional behavior of curved girders are coupled due to their initial curvature. Also, the behavior is affected by cross beams. To investigate the behavior of horizontally curved I-girder bridges, it is necessary to consider curved girders with cross beams. In order to perform free vibration analyses of horizontally curved I-girder bridges, a finite element formulation is presented here and a finite element analysis program is developed. The formulation that is presented here consists of curved and straight beam elements, including the warping degree of freedom. Based on the theory of thin-walled curved beams, the shape functions of the curved beam elements are derived from homogeneous solutions of the static equilibrium equations. Third-order hermits polynomials are used to form the shape functions of the straight beam elements. In the finite element analysis program, global stiffness and mass matrix are composed, based on the Cartesian coordinate system. The Gupta method is used to efficiently solve the eigenvalue problem. Comparing the results of several examples here with those of previous studies, the formulation presented is verified. The validity of the program developed is shown by comparing results with those analyzed by the shell element.

• Investigative Magnetic Resonance Imaging
• /
• v.18 no.2
• /
• pp.107-119
• /
• 2014

Purpose : To efficiently evaluate phased array coil performance using a software tool box with which we can make visual comparison of the sensitivity of every coil element between the real experiment and EM simulation. Materials and Methods: We have developed a $C^{{+}{+}}$- and MATLAB-based software tool called Phased Array Coil Evaluator (PACE). PACE has the following functions: Building 3D models of the coil elements, importing the FDTD simulation results, and visualizing the coil sensitivity of each coil element on the ordinary Cartesian coordinate and the relative coil position coordinate. To build a 3D model of the phased array coil, we used an electromagnetic 3D tracker in a stylus form. After making the 3D model, we imported the 3D model into the FDTD electromagnetic field simulation tool. Results: An accurate comparison between the coil sensitivity simulation and real experiment on the tool box platform has been made through fine matching of the simulation and real experiment with aids of the 3D tracker. In the simulation and experiment, we used a 36-channel helmet-style phased array coil. At the 3D MRI data acquisition using the spoiled gradient echo sequence, we used the uniform cylindrical phantom that had the same geometry as the one in the FDTD simulation. In the tool box, we can conveniently choose the coil element of interest and we can compare the coil sensitivities element-by-element of the phased array coil. Conclusion: We expect the tool box can be greatly used for developing phased array coils of new geometry or for periodic maintenance of phased array coils in a more accurate and consistent manner.