• Geophysics and Geophysical Exploration
• /
• v.15 no.2
• /
• pp.75-84
• /
• 2012

The sea layer in marine Controlled-Source Electromagnetic (mCSEM) survey changes the conventional definition of apparent resistivity which is used in the land CSEM survey. Thus, the development of a new algorithm, which computes apparent resistivity for mCSEM survey, can be an initiative of mCSEM data interpretation. First, we compared and analyzed electromagnetic responses of the 1D stratified gas hydrate model and the half-space model below the sea layer. Amplitude and phase components showed proper results for computing apparent resistivity than real and imaginary components. Next, the amplitude component is more sensitive to the subsurface resistivity than the phase component in far offset range and vice versa. We suggested the induction number as a selection criteria of amplitude or phase component to calculate apparent resistivity. Based on our study, we have developed a numerical algorithm, which computes appropriate apparent resistivity corresponding to measured mCSEM data using grid search method. In addition, we verified the validity of the developed algorithm by applying it to the stratified gas hydrate models with various model parameters. Finally, by constructing apparent resistivity pseudo-section from the mCSEM responses with 2D numerical models simulating gas hydrate deposits in the Ulleung Basin, we confirmed that the apparent resistivity can provide the information on the geometric distribution of the gas hydrate deposit.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.4
• /
• pp.114-131
• /
• 1992

In the present work, a two-dimensional boundary-value problem for a large amplitude motion is treated as an initial-value problem by satisfying the exact body-boundary and nonlinear free-surface boundary conditions. The present nonlinear numerical scheme is similar to that described by Vinje and Brevig(1981) who utilized the Cauchy's theorem and assumed the periodicity in the horizontal coordinate. In the present thesis, however, the periodicity in the horizontal coordinate is not assumed. Thus the present method can treat more realistic problems, which allow radiating waves to infinities. In the present method of solution, the original infinite fluid domain, is divided into two subdomains ; ie the inner and outer subdomains which are a local nonlinear subdomain and the truncated infinite linear subdomain, respectively. By imposing an appropriate matching condition, the computation is carried out only in the inner domain which includes the body. Here we adopt the nonlinear scheme of Vinje & Brevig only in the inner domain and respresent the solution in the truncated infinite subdomains by distributing the time-dependent Green function on the matching boundaries. The matching condition is that the velocity potential and stream function are required to be continuous across the matching boundary. In the computations we used, if necessary, a regriding algorithm on the free surface which could give converged stable solutions successfully even for the breaking waves. In harmonic oscillation problem, each harmonic component and time-mean force are obtained by the Fourier transform of the computed forces in the time domain. The numerical calculations are made for the following problems. $\cdot$ Forced harmonic large-amplitude oscillation(${\omega}{\neq}0,\;U=0$) $\cdot$ Translation with a uniform speed(${\omega}=0,\;U{\neq}0$) The computed results are compared with available experimental data and other analytical results.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.7 no.2
• /
• pp.176-197
• /
• 1995

Two-dimensional and three-dimensional structures of tidal currents on southeastern waters of Korea off Kyungnam coast were investigated via a series of numerical models based on dynamic principles. With a two-dimensional tidal model, tidal regimes of major eight tidal constituents (M$_2$, S$_2$, K$_1$, O$_1$, N$_2$, K$_2$, P$_1$, Q$_1$) were computed. Model results showed that the computed results were in good agreement with coastal observations. On the basis of these results the model was further improved to compute three-dimensional structure of tidal current in inner Jinhai and Masan Bay regions of the model area where severe pollutions occur due to red tide by combination of the previous two-dimensional model and inner three-dimensional model. For this work, three-dimensional Galerkin-Spectral model and two-dimensional depth-integrated model are dynamically combined by the method presented by Davies (1980). In addition to the previous work by Davies, the advective term and quadratic bottom friction term are included in present Three-dimensional numerical model. The computed results of M$_2$ tidal current ellipses with respect to depth showed general agreements with those of current observations by KORDI (1990).

• Journal of Korean Society of Environmental Engineers
• /
• v.32 no.9
• /
• pp.862-869
• /
• 2010

Numerical analysis was done to evaluate the gas flow distribution in small scale SCR system which has $2.4{\times}2.4{\times}3.1\;m^3$ in volume and 25,300 Sm3/hr in flue gas flow capacity. Various types of baffles proposed for controlling the flow uniformity were evaluated by the CFD analysis to find the optimal geometry of the baffle in the SCR system. By installing baffles in the SCR system, the RMS (%) value was raised up to 6.2% compared with the baffle-uninstalled state. The effect of baffle thicknesses on the RMS (%) value was not shown within 0 and 8 mm in thickness, but the RMS (%) value was raised by 2.5% in 10 mm of baffles thickness, which causes the unstability in flow. By comparison between the shape of baffles, it is known that the lattice type baffle has better performance in controlling the flow uniformity than the circular truncated cone type baffle or mixer type baffle. RMS (%) values have more that 10% difference according to the shape of baffle type.

• Journal of The Korean Astronomical Society
• /
• v.47 no.2
• /
• pp.77-86
• /
• 2014

Using a cosmological ${\Lambda}CDM$ simulation, we analyze the differences between the widely-used spin parameters suggested by Peebles and Bullock. The dimensionless spin parameter ${\lambda}$ proposed by Peebles is theoretically well-justified but includes an annoying term, the potential energy, which cannot be directly obtained from observations and is computationally expensive to calculate in numerical simulations. The Bullock's spin parameter ${\lambda}^{\prime}$ avoids this problem assuming the isothermal density profile of a virialized halo in the Newtonian potential model. However, we find that there exists a substantial discrepancy between ${\lambda}$ and ${\lambda}^{\prime}$ depending on the adopted potential model (Newtonian or Plummer) to calculate the halo total energy and that their redshift evolutions differ to each other significantly. Therefore, we introduce a new spin parameter, ${\lambda}^{\prime\prime}$, which is simply designed to roughly recover the value of ${\lambda}$ but to use the same halo quantities as used in ${\lambda}^{\prime}$. If the Plummer potential is adopted, the ${\lambda}^{\prime\prime}$ is related to the Bullock's definition as ${\lambda}^{\prime\prime}=0.80{\times}(1+z)^{-1/12}{\lambda}^{\prime}$. Hence, the new spin parameter ${\lambda}^{\prime\prime}$ distribution becomes consistent with a log-normal distribution frequently seen for the ${\lambda}^{\prime}$ while its mean value is much closer to that of ${\lambda}$. On the other hand, in case of the Newtonian potential model, we obtain the relation of ${\lambda}^{\prime\prime}=(1+z)^{-1/8}{\lambda}^{\prime}$; there is no significant difference at z = 0 as found by others but ${\lambda}^{\prime}$ becomes more overestimated than ${\lambda}$ or ${\lambda}^{\prime\prime}$ at higher redshifts. We also investigate the dependence of halo spin parameters on halo mass and redshift. We clearly show that although the ${\lambda}^{\prime}$ for small-mass halos with $M_h$ < $2{\times}10^{12}M_{\odot}$ seems redshift independent after z = 1, all the spin parameters explored, on the whole, show a stronger correlation with the increasing halo mass at higher redshifts.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.17 no.11 s.114
• /
• pp.1120-1131
• /
• 2006

This paper describes the numerical analysis techniques using the Combined Integration/Matrix Method to calculate ground potential rise which can be occurred in the various grounding systems. Combined Integration/Matrix Method is used to reduce the error and computation time with the analytical integration equation and the proper segmentaion of earth embedded conductor. To do it, optimal segmentaion method for the buried conductors is presented through error analysis which is capable of applying the practical scaled various grounding systems. The optimum length of segmented element is much co-related with the buried depth of grounding electrode and the maximum length of buried electrode. As a result, less 3 precent errors was obtained by proposed model. The proposed model is applied to verify an effect of multi-grounding problems which was aroused much controversy with separated or common grounding between the high power grounding system and low power grounding system such as signal and telecommunication grounding.

• KEPCO Journal on Electric Power and Energy
• /
• v.1 no.1
• /
• pp.157-162
• /
• 2015

Large scale superconducting magnetic energy storage (SMES) system requires very high magnetic energy density in its superconducting coils to enhance the energy capacity and efficiency of the system. The recent high temperature superconducting (HTS) conductors, so called 2G conductors, show very good performance under very high magnetic field so that they seem to be perfect materials for the large scale SMES coils. A general shape of the coil system with the 2G HTS conductor has been a tor oid, because the magnetic field applied perpendicularly to the surface of the 2G HTS conductor could be minimized in this shape of coil. However, a toroid coil requires a 3-dimensional computation to acquire the characteristics of its critical current density - magnetic field relations which needs very complicated numerical calculation, very high computer specification, and long calculation time. In this paper, we suggested an analytic and statistical calculation method to acquire the maximum magnetic flux density applied perpendicularly to the surface of the 2G HTS conductor and the stored energy in the toroid coil system. Although the result with this method includes some errors but we could reduce these errors within 5 percent to get a reasonable estimation of the important parameters for design process of the HTS toroid coil system. As a result, the calculation time by the suggested method could be reduced to 0.1 percent of that by the 3-dimensional numerical calculation.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.5
• /
• pp.287-296
• /
• 2019

In this study, the effect of ground boundary conditions on the evaluation of blast resistance performance of precast arch structures was evaluated by a numerical analysis method. Two types of boundary conditions, namely, fixed boundary conditions and a perfectly matched layer (PML) were applied to numerical models. Blast loads that were much higher than the design load of the target structure were applied to compare the effects of the boundary conditions. The distribution and path of the ground explosion pressure, structural displacement, fracture of concrete, stress of concrete, and reinforcing bars were compared according to the ground boundary condition settings. As a result, the reflecting pressure shock wave at the ground boundaries could be effectively eliminated using PML elements; furthermore, the displacement of the foundation was reduced. However, no distinct difference could be observed in the overall structural behavior including the fracture and stress of the concrete and rebar. Therefore, when blast simulations are performed in the design of protective structures, it is rational to apply the fixed boundary condition on the ground boundaries as conservative design results can be achieved with relatively short computation times.

• The Journal of the Acoustical Society of Korea
• /
• v.40 no.4
• /
• pp.261-269
• /
• 2021

Without any validation of the incompressible assumption, most of previous studies on cavitation flow and its noise have utilized numerical methods based on the incompressible Reynolds Average Navier-Stokes (RANS) equations because of advantage of its efficiency. In this study, to investigate the effects of the flow compressibility on the Tip Vortex Cavitation (TVC) flow and noise, both the incompressible and compressible simulations are performed to simulate the TVC flow, and the Ffowcs Williams and Hawkings (FW-H) integral equation is utilized to predict the TVC noise. The DARPA Suboff submarine body with an underwater propeller of a skew angle of 17 degree is targeted to account for the effects of upstream disturbance. The computation domain is set to be same as the test-section of the large cavitation tunnel in Korea Research Institute of Ships and Ocean Engineering to compare the prediction results with the measured ones. To predict the TVC accurately, the Delayed Detached Eddy Simulation (DDES) technique is used in combination with the adaptive grid techniques. The acoustic spectrum obtained using the compressible flow solver shows closer agreement with the measured one.

• The Journal of the Acoustical Society of Korea
• /
• v.42 no.1
• /
• pp.7-15
• /
• 2023

Recently, as interest in air quality in vehicles increases, the use of fine dust detection sensors for air quality measurement is becoming common. An axial-flow fan is inserted in the fine dust sensor installed in the air conditioning system in the vehicle to prevent dust from sinking directly on the sensor. When the sensor operates, the flow noise caused by the rotation of the axial-flow fan acts as a major noise source of the fine dust sensor. flow noise is recognized as one of the product competitiveness of fine dust sensors. In this study, the noise was gradually reduced at the same flow rate by improving the flow performance of the small axial flow fan. First, a virtual fan performance tester consisting of about 20 million grids was developed to analyze the aerodynamic performance of the target small axial-flow fan. In addition, the flow field was simulated by using compressible Large Eddy Simulation for direct computation of flow noise as well as high-accurate prediction of flow rate. The validity of numerical method are confirmed through the comparison of predicted results with experimental ones. After the effects of pitch angle on flow performance were analyzed using the verified numerical method, the pitch angle was determined to maximize the flow rate. It was found that the flow rate was increased by 8.1 % and noise was reduced by 0.8 dBA when the axial-flow fan with the optimum pitch angle was used.