• Geophysics and Geophysical Exploration
• /
• v.7 no.1
• /
• pp.80-87
• /
• 2004

During excavation work for the construction of a highway bridge in a limestone area in Korea, several cavities were found, and construction work was stopped temporarily. Cavities under the bridge piers might seriously threaten the safety of the planned bridge, because they could lead to excessive subsidence and differential settlement of the pier foundations. In order to establish a method for reinforcement of the pier foundations, borehole radar reflection and tomography surveys were carried out, to locate cavities under the planned pier locations and to determine their sizes where they exist. Since travel time data from the crosshole radar survey showed anisotropy, we applied an anisotropic tomography inversion algorithm assuming heterogeneous elliptic anisotropy, in order to reconstruct three kinds of tomograms: tomograms of maximum and minimum velocities, and of the direction of the symmetry axis. The distribution of maximum velocity matched core logging results better than that of the minimum velocity. The degree of anisotropy, defined by the normalized difference between maximum and minimum velocities, was helpful in deciding whether an anomalous zone in a tomogram was a cavity or not. By careful examination of borehole radar reflection and tomography images, the spatial distributions of cavities were delineated, and most of them were interpreted as being filled with clay and/or water. All the interpretation results implied that two faults imaged clearly by a DC resistivity survey were among the most important factors controlling the groundwater movement in the survey area, and therefore were closely related to the development of cavities. The method of reinforcement of the pier foundations was based on the interpretation results, and the results were confirmed when construction work was resumed.

• Journal of Astronomy and Space Sciences
• /
• v.14 no.1
• /
• pp.136-141
• /
• 1997

To calculate the position and velocity of the artificial satellite precisely, one has to build a mathematical model concerning the perturbations by understanding and analysing the space environment correctly and then quantifying. Due to these space environment model, the total acceleration of the artificial satellite can be expressed as the 2nd order differential equation and we build an orbit propagation algorithm by integrating twice this equation by using the Cowell's method which gives the position and velocity of the artificial satellite at any given time. Perturbations important for the orbits of geostationary spacecraft are the Earth's gravitational potential, the gravitational influences of the sun and moon, and the solar radiation pressure. For precise orbit propagation in Cowell' method, 40 x 40 spherical harmonic coefficients can be applied and the JPL DE403 ephemeris files were used to generate the range from earth to sun and moon and 8th order Runge-Kutta single step method with variable step-size control is used to integrate the the orbit propagation equations.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.62 no.6
• /
• pp.833-839
• /
• 2013

Cardiopulmonary resuscitation (CPR) is performed by thoracic compression and artificial ventilation for the patient under emergent situation to maintain at least the minimum level of respiration and blood circulation for life survival. Good quality CPR requires monitoring respiration, however, traditional respiratory air flow transducers cannot be used because the transducer elements are facing the whole area perpendicular to the flow axis. The present study developed a new air flow transducer conveniently applicable to CPR. Specially designed "sensing rod" samples the air velocity at 3 different locations of the flow cross-section, then transforms into average dynamic pressure by the Bernoulli's law. The symmetric structure of the sensing holes of the sensing rod enables bi-directional measurement simply by taking the difference in pressure by a commercial differential pressure transducer. Both inspiratory and expiratory flows were obtained with symmetric measurement characteristics. Quadratic curve fitting provided excellent calibration formula with a correlation coefficient>0.999 (P<0.0001) and the mean relative error<1%. The present results can be usefully applied to accurately monitor the air flow rate during CPR.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.21 no.5
• /
• pp.577-582
• /
• 2015

This paper is the fundamental study for the application of cold storage system to the transportation equipment by sea and land. This numerical study presents the solid-liquid phase change phenomenon of calcium chloride solution of 30wt %. The governing equations are 1-dimensional unsteady state heat transfer equations of $1^{st}$ order partial differential equations. This type of latent heat storage material is often usable in fishery vessel for controlling the temperature of container with constant condition. The governing equation was discretized with finite difference method and the program was composed with Mathcad program. The main parameters of this solution were the initial temperature of heat storage material, ambient temperature of cold air and the velocity of cold air. The data of boundary layer thickness becomes thin with the increasing of cold air flowing velocity and also the heat storage completion time become shorten.

• Ocean Systems Engineering
• /
• v.7 no.4
• /
• pp.435-460
• /
• 2017

The main objective of this study is to investigate the turning and zig-zag maneuvering performance of the well-known naval surface combatant DTMB (David Taylor Model Basin) 5415 hull with URANS (Unsteady Reynolds-averaged Navier-Stokes) method. Numerical simulations of static drift tests have been performed by a commercial RANS solver based on a finite volume method (FVM) in an unsteady manner. The fluid flow is considered as 3-D, incompressible and fully turbulent. Hydrodynamic analyses have been carried out for a fixed Froude number 0.28. During the analyses, the free surface effects have been taken into account using VOF (Volume of Fluid) method and the hull is considered as fixed. First, the code has been validated with the available experimental data in literature. After validation, static drift, static rudder and drift and rudder tests have been simulated. The forces and moments acting on the hull have been computed with URANS approach. Numerical results have been applied to determine the hydrodynamic maneuvering coefficients, such as, velocity terms and rudder terms. The acceleration, angular velocity and cross-coupled terms have been taken from the available experimental data. A computer program has been developed to apply a fast maneuvering simulation technique. Abkowitz's non-linear mathematical model has been used to calculate the forces and moment acting on the hull during the maneuvering motion. Euler method on the other hand has been applied to solve the simultaneous differential equations. Turning and zig-zag maneuvering simulations have been carried out and the maneuvering characteristics have been determined and the numerical simulation results have been compared with the available data in literature. In addition, viscous effects have been investigated using Eulerian approach for several static drift cases.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.123-130
• /
• 2002

A general formulation for shape design sensitivity analysis over three dimensional beam structure is developed based on a variational formulation of the beam in linear elasticity. Sensitivity formula is derived based on variational equations in cartesian coordinates using the material derivative concept and adjoint variable method for the displacement and Von-Mises stress functionals. Shape variation is considered for the beam shape in general 3-dimensional direction as well as for the orientation angle of the beam cross section. In the sensitivity expression, the end points evaluation at each beam segment is added to the integral formula, which are summed over the entire structure. The sensitivity formula can be evaluated with generality and ease even by employing piecewise linear design velocity field despite the bending model is fourth order differential equation. For the numerical implementation, commercial software ANSYS is used as analysis tool for the primal and adjoint analysis. Once the design variable set is defined using ANSYS language, shape and orientation variation vector at each node is generated by making finite difference to the shape with respect to each design parameter, and is used for the computation of sensitivity formula. Several numerical examples are taken to show the advantage of the method, in which the accuracy of the sensitivity is evaluated. The results are found excellent even by employing a simple linear function for the design velocity evaluation. Shape optimization is carried out for the geometric design of an archgrid and tilted bridge, which is to minimize maximum stress over the structure while maintaining constant weight. In conclusion, the proposed formulation is a useful and easy tool in finding optimum shape in a variety of the spatial frame structures.

• Journal of Radiation Protection and Research
• /
• v.28 no.2
• /
• pp.97-108
• /
• 2003

A three-dimensional wind field model based on the variational technique has been developed for estimating the overall wind patterns over a complex terrain. The three-dimensional elliptic partial differential equations on Cartesian and terrain-following coordinates have been established to obtain the Lagrangian multiplier and the adjusted wind velocity. The simulations were performed to evaluate the variations of the velocity vectors on the hemisphere, half-cylinder, and saddle type obstacles. Also, the wind field model in the terrain-following coordinate has been applied for evaluating the characteristics of wind patterns according to the variations of Gauss precision moduli on the hemispheric topography. The results showed that the horizontal and vertical wind components were strongly governed by the selection of the values of Gauss precision moduli.

• Journal of Navigation and Port Research
• /
• v.26 no.2
• /
• pp.221-226
• /
• 2002

To provide reliable data for drift prediction models, field experiments were carried out in the coastal region off Busan port. Four different size of vessels(10, 30, 50, 90G/T ton) were deployed for the experiment. Among them G/T 50ton class vessel was equipped with instruments measuring the currents, winds, headings and trajectory data. In the rest of vessels only the position data were recorded for the purpose of target divergence study. The trajectories of each vessel were measured by DGPS(Differential Global Positioning System) and collected by APRS(Automatic Position Reporting System). The experiment was done in wind of 2~10m/s and current of 0.5~1.5m/s. The leeway was derived by subtracting surface current velocity from target drifting velocity. The leeway rate of G/T 50ton vessel was found to be about 3.6% and the computed leeway speed equation was $U_L$=0.042 W - 0.034. The processed leeway angle data were deflected by $-30^{\circ}$~$40^{\circ}$ from the direction of ship drift.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.9 no.2
• /
• pp.97-104
• /
• 2005

A simple but efficient grid generation technique by using the modified compressible form of stream function has been formulated. Transformation of a physical plane to a streamline plane, the Von Mises Transformation, has been widely used to solve the differential equations governing flow phenomena, however, limitation arises in low velocity region of boundary layer, mixing layer and wake region where the relatively large grid spacing is inevitable. Modified Von Mises Transformation with simple mathematical adjustment for the stream function is suggested and applied to solve the confined coaxial turbulent jet mixing with simple $\kappa-\epsilon$ turbulence model. Comparison with several experimental data of axial mean velocity, turbulent kinetic energy, and Reynolds shear stress distribution shows quite good agreement in the mixing layer except in the centerline where the turbulent kinetic energy distributions were somewhat under estimated. This formulation is strongly suggested to be utilized specially for free turbulent mixing layers in axisymmetric flow conditions such as the investigation of mixing behavior, jet noise production and reduction for Turbofan engines.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.1A
• /
• pp.95-104
• /
• 2008

Dynamic behaviors of the two-span continuous bridge which is one of prototypes on Gyoung-Bu high-speed railway are analyzed and some methods for reducing the resonance of the bridge are proposed. The bridge is modeled as a two-span continuous beam and the load is a vehicle of TGV-K (2p+18T) with length of 380.15 meter traveling on the railway bridge at some constant velocity. The equations governing the dynamic behaviors of the bridge are partial differential equations produced by using a system with distributed mass and elasticity. The analysis of the governing equations is performed by the mode superposition method which has modal coordinates solved by Duhamel's integral. Without considering the train velocity the dynamic reponses can be greatly reduced at some special lengths of bridge. It is different from the results of simple bridges researched so far. When the dynamic responses increase rapidly to make a resonance phenomenon depending on the train velocities, the several methods are proposed to deduce the resonance.