• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.38-47
• /
• 2005

Substructuring methods are often used in finite element structural analyses. In this study a multi-level substructuring(MLSS) algorithm is developed and proposed as a possible candidate for finite element fluid solvers. The present algorithm consists of four stages such as a gathering, a condensing, a solving and a scattering stage. At each level, a predetermined number of elements are gathered and condensed to form an element of higher level. At the highest level, each sub-domain consists of only one super-element. Thus, the inversion process of a stiffness matrix associated with internal degrees of freedom of each sub-domain has been replaced by a sequential static condensation of gathered element matrices. The global algebraic system arising from the assembly of each sub-domain matrices is solved using a well-known iterative solver such as the conjugare gradient(CG) or the conjugate gradient squared(CGS) method. A time comparison with CG has been performed on a 2-D Poisson problem. With one domain the computing time by MLSS is comparable with that by CG up to about 260,000 d.o.f. For 263,169 d.o.f using 8 x 8 sub-domains, the time by MLSS is reduced to a value less than $30\%$ of that by CG. The lid-driven cavity problem has been solved for Re = 3200 using the element interpolation degree(Deg.) up to cubic. in this case, preconditioning techniques usually accompanied by iterative solvers are not needed. Finite element formulation for the incompressible flow has been stabilized by a modified residual procedure proposed by Ilinca et al.[9].

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• v.15 no.12
• /
• pp.4567-4583
• /
• 2021

This study proposes an analytical approximation algorithm based on extreme value theory (EVT) for the inverse of the power of the incomplete Gamma function. First, the Gumbel function is used to approximate the power of the incomplete Gamma function, and the corresponding inverse problem is transformed into the inversion of an exponential function. Then, using the tail equivalence theorem, the normalized coefficient of the general Weibull distribution function is employed to replace the normalized coefficient of the random variable following a Gamma distribution, and the approximate closed form solution is obtained. The effects of equation parameters on the algorithm performance are evaluated through simulation analysis under various conditions, and the performance of this algorithm is compared to those of the Newton iterative algorithm and other existing approximate analytical algorithms. The proposed algorithm exhibits good approximation performance under appropriate parameter settings. Finally, the performance of this method is evaluated by calculating the thresholds of space-time block coding and space-frequency block coding pattern recognition in multiple-input and multiple-output orthogonal frequency division multiplexing. The analytical approximation method can be applied to other related situations involving the maximum statistics of independent and identically distributed random variables following Gamma distributions.

• Tunnel and Underground Space
• /
• v.19 no.5
• /
• pp.407-420
• /
• 2009

In this study geostatistical technique using indicator kriging was performed to evaluate the optimal rock mass classification by integrating the various geophysical information such as borehole data and geophysical data. To get the optimal kriging result, it is necessary to devise the suitable technique to integrate the hard (borehole) and soft (geophysical) data effectively. Also, the model parameters of the variogram must be determined as a priori procedure. Iterative non-linear inversion method was implemented to determine the model parameters of theoretical variogram. To verify the algorithm, behaviour of object function and precision of convergence were investigated, revealing that gradient of the range is extremely small. This algorithm for the field data was applied to a mountainous area planned for a large-scale tunneling construction. As for a soft data, resistivity information from AMT survey is incorporated with RMR information from borehole data, a sort of hard data. Finally, RMR profiles were constructed and attempted to be interpreted at the tunnel elevation and the upper 1D level.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2009.10a
• /
• pp.75-80
• /
• 2009

Despite its flexibility to complex geometry, three-dimensional (3D) electromagnetic(EM) modeling schemes using finite element method (FEM) have been faced to practical limitation due to the resulting large system of equations to be solved. An efficient 3D FEM modeling scheme has been developed, which can adopt either direct or iterative solver depending on the problems. The direct solver PARDISO can reduce the computing time remarkably by incorporating parallel computing on multi-core processor systems, which is appropriate for single frequency multi-source configurations. When limited memory, the iterative solver BiCGSTAB(1) can provide fast and stable convergence. Efficient 3D simulations can be performed by choosing an optimum solver depending on the computing environment and the problems to be solved. This modeling includes various types of controlled-sources and can be exploited as an efficient engine for 3D inversion.

• Geophysics and Geophysical Exploration
• /
• v.18 no.1
• /
• pp.1-8
• /
• 2015

To estimate the shear-wave velocity (${\nu}_s$ beneath the OTVZ seismic station on Ngauruhoe volcano in New Zealand, we calculated receiver functions (RFs) using 127 teleseismic data ($Mw{\geq}5.5$) with high signal-to-noise ratios recorded during November 11, 2011 to September 11, 2013. The genetic inversion algorithms was applied to 21 RFs calculated by the iterative time-domain deconvolution method. In the 1-D ${\nu}_s$ model derived by the inversion, the Moho is observed at a 14 km depth, marked by a ${\nu}_s$ transition from 3.7 km/s to 4.7 km/s. The average ${\nu}_s$ of the overlying crust is 3.4 km/s, and the average ${\nu}_s$ of a greater than 9-km thick low-velocity layer (LVL) in the lower crust is 3.1 km/s. The LVL becomes thinner with increasing distance from the station. Another LVL thicker than 10 km with ${\nu}_s$ less than 4.3 km/s is found in the upper mantle. Those LVLs in the lower crust and the upper mantle and the relatively thin crust might be related to the magma activity caused by the subducting Pacific plate.

• Geophysics and Geophysical Exploration
• /
• v.17 no.2
• /
• pp.66-73
• /
• 2014

To estimate the shear-velocity ($v_s$) structure beneath the WIZ station on White Island in New Zealand, we applied receiver function (RF) inversion and H-${\kappa}$ stacking methods to 362 teleseismic events (Mw > 5.5) recorded during April 20, 2007 to September 6, 2013. Using 71 RFs with errors less than 20% after 200 iterative computations, we determined that the depth to Moho of $v_s$ = 4.35 km/s is $24{\pm}1km$ within a 15 km radius of the station. In an 1-d $v_s$ model derived by RF inversions, a 4-km thick low-velocity layer (LVL) at depths of 18 ~ 22 km was identified in the lower crust. This LVL, which is 0.15 km/s slower than the rocks above and below it, may indicate the presence of a deep magma reservoir. The H-${\kappa}$ stacking method yielded an estimate of the depth to the Moho of 24.5 km, which agrees well with the depth determined by RF inversions. The low $v_p/v_s$ ratio of 1.64 may be due to the presence of gas-filled rock or hot crystallizing magma.

• The Journal of the Acoustical Society of Korea
• /
• v.36 no.5
• /
• pp.293-299
• /
• 2017

Small unmanned platforms maneuvering underwater are the key naval future forces, utilized as the asymmetric power in war. As a method of detecting and identifying such platforms, we introduce a property estimation technique based on an iterative numerical analysis. The property estimation technique can estimate not only the position of a target but also its physical properties. Moreover, it will have a potential in detecting and classifying still target or multiple targets. In this study, we have conducted the property estimation of an small underwater target using the data acquired from the lake experiment. As a result, it shows that the properties of a small platform may be roughly estimated from the in site data even using one channel.

• Geophysics and Geophysical Exploration
• /
• v.14 no.2
• /
• pp.140-145
• /
• 2011

In order to analyze the influence of the noise on a cross-well traveltime tomography to detect a small scale low velocity body in a homogeneous medium, the first arrival travel times were computed one a tunnel model by a finite-difference ray tracing scheme. Three different types and four different intensity levels of white noises were added to the computed first arrival travel times, and velocity tomograms were constructed using an iterative inversion method (SIRT). Tomograms with the noise intensity up to 10% of the maximum traveltime delay in the tunnel model, showed the exact location of the tunnel. However, the velocity shown at the tunnel location was not close to air velocity but only slightly less than the velocity of the background medium. The additive random noise showed significantly less degree of influence on the resulting tomogram than the source- and receiver consistent noise.

• Geophysics and Geophysical Exploration
• /
• v.7 no.2
• /
• pp.148-154
• /
• 2004

This article reviews the development of three-dimensional (3-D) magnetotelluric (MT) modeling. The 3-D modeling of electromagnetic fields is essential in understanding the physics of MT soundings, and in implementing an inversion method to reconstruct a 3-D resistivity image. Although various numerical schemes have been developed over the last two decades, practical methods have been quite limited. However, the recent rapid improvement in computer speed and memory, as well as the advance in iterative solution algorithms for a large system of equations, makes it possible to model the MT responses of complex 3-D structures, which have been very difficult to simulate before. The use of staggered grids in finite difference method has become popular, conserving a magnetic flux and an electric current and allowing for realistic discontinuous fields. The convergence of numerical solutions has been greatly accelerated by adopting Krylov subspace methods, proper preconditioning techniques, and static divergence corrections. The vector finite-element method using edge elements is also free from the discontinuity problem, and seems a natural choice for modeling complex structures including irregular topography because its flexibility allows one to capture full geometric complexity.

• Journal of the Korean earth science society
• /
• v.30 no.2
• /
• pp.153-164
• /
• 2009

We examined the regional 1-D deep resistivity structure of the Korean Peninsula using MT data acquired at seven sites located in the Kyongsang Basin and Kyonggi Massif. At the sites located in the Kyongsang Basin, surrounding sea distorts observed MT response and hence this distortion, so called "sea effect", is corrected using an iterative tensor stripping method. The 1-D layered inversion results for the seven MT sites reveal 4 layered structure, which is composed of 1) near surface layer, 2) upper crust, 3) lower crust and upper mantle, and 4) asthenosphere from the surface downward. Conrad interface, which is a boundary between upper and lower crust, is distinctly identified beneath all the MT sites. Conrad interface depth is estimated to about be 17km in the Kyongsang Basin and about 12km in the Kyonggi Massif, while the upper crust of the Kyongsang Basin is about 5 times more resistive than that of the Kyonggi Massif. Finally, asthenosphere is inferred to exist below a depth of approximately 100km with a resistivity of 200-300 ohm-m.