• Journal of The Korean Society of Agricultural Engineers
• /
• v.60 no.1
• /
• pp.59-65
• /
• 2018

The object of this paper was to analyze combined load acting on agricultural reservoir. This study was carried out to 3-D numerical modeling for displacement characteristic and seismic acceleration characteristic. The results of study were analyzed and summarized as follow. It was found that the displacement caused by combined load acting on railway and agricultural reservoir did not reflect the effect of load and the seismic wave consistently. The ground accelerations that occur in railway and dam were amplified because 3-D numerical analysis program interprets ground as an elastic body. Actual ground shows characteristics of elasticity and plasticity, so measured values will show different tendency. As a result of analyzing displacement characteristics, it is considered to be related to stiffness. The Ofunato seismic wave, the displacement (77.1 mm) of the body satisfied the allowable displacement (220 mm), but The Hachinohe seismic wave (282.8 mm) did not. It is considered that displacement caused by combined load is affected not only by acceleration but also by characteristics of materials.

• Geophysics and Geophysical Exploration
• /
• v.4 no.3
• /
• pp.70-79
• /
• 2001

The integral equation method is a powerful tool for numerical electromagnetic modeling. But the difficulty of this technique is the size of the linear equations, which demands excessive memory and calculation time to invert. This limitation of the integral equation method becomes critical in inverse problem. The conventional Born approximation, where the electric field in the anomalous body is approximated by the background field, is very rapid and easy to compute. However, the technique is inaccurate when the conductivity contrast between the body and the background medium is large. Quasi-linear, quasi-analytical and extended Born approximations are novel approaches to 3-D EM modeling based on the linearization of the integral equations for scattered EM field. These approximation methods are much less time consuming than full integral equation method and more accurate than conventional Born approximation. They we, however, still approximate methods for 3-D EM modeling. Iterative series methods such as modified Born, quasi-linear and quasi-analytical can be used to increase the accuracy of various approximation methods. Comparisons of numerical performance against a full integral equation and various approximation codes show that the iterative series methods are very accurate and almost always converge. Furthermore, they are very fast and easy to implement on a computer. In this study, extended Born series method is developed and it shows more accurate result than that of other series methods. Therefore, Iterative series methods, including extended Born series, open principally new possibilities for fast and accurate 3-D EM modeling and inversion.

• Computers and Concrete
• /
• v.18 no.3
• /
• pp.319-336
• /
• 2016

In this study, both two- and three-dimensional (2D and 3D) finite-volume-based models were developed to analyze the heat transfer mechanisms through the porous structures of cellular concretes under steady-state heat transfer conditions and to investigate the differences between the 2D and 3D modeling results. The 2D and 3D reconstructed pore networks were generated from the microstructural information measured by 3D images captured by X-ray computerized tomography (X-CT). The computed effective thermal conductivities based on the 2D and 3D calculations performed on the reconstructed porous structures were found to be nearly identical to those evaluated from the 2D cross-sectional images and the 3D X-CT images, respectively. In addition, the 3D computed effective thermal conductivity was found to agree better with the measured values, in comparison with the 2D reconstruction and real cross-sectional images. Finally, the thermal conductivities computed for different reconstructed porous 3D structures of cellular concretes were compared with those obtained from 2D computations performed on 2D reconstructed structures. This comparison revealed the differences between 2D and 3D image-based modeling. A correlation was thus derived between the results of the 3D and 2D models.

• Structural Engineering and Mechanics
• /
• v.26 no.3
• /
• pp.343-361
• /
• 2007

This article presents a numerical tool for dimensioning two-threaded fasteners connecting prismatic parts subjected to fatigue tension loads that are coplanar with the screw axis. A simplified numerical model is developed from unidirectional finite elements, modeling the connected parts and screws with bent elements and the elastic contact layer between the parts with springs. An algorithm updating the contact stiffness matrix, calculating forces and displacements at each node of the structure and thus normal stresses in the screws in both static and fatigue is further developed using C language. An experimental study is also conducted in parallel with the numerical approach to validate the developed model assumptions, the numerical model and the 3D finite element results. Since stiffness values for the compressive zones in the parts are analytically difficult to determine, a statistical software method is used, from which a tuning factor is derived for identifying these stiffness values. The method is also applied to set out the influence of each parameter on the fatigue behaviour of each screw. Finally, the developed model will be used to establish a new, sophisticated, fast and accurate tool for dimensioning bolted mechanical structures.

• Geophysics and Geophysical Exploration
• /
• v.20 no.3
• /
• pp.129-136
• /
• 2017

3D seismic data processing methods such as full waveform inversion or reverse-time migration require 3D wave propagation modeling and heavy calculations. We compared efficiency and accuracy of a Xeon Phi coprocessor to those of a high-end server CPU using 3D frequency-domain wave propagation modeling. We adopted the OpenMP parallel programming to the time-domain finite difference algorithm by considering the characteristics of the Xeon Phi coprocessors. We applied the Fourier transform using a running-integration to obtain the frequency-domain wavefield. A numerical test on frequency-domain wavefield modeling was performed using the 3D SEG/EAGE salt velocity model. Consequently, we could obtain an accurate frequency-domain wavefield and attain a 1.44x speedup using the Xeon Phi coprocessor compared to the CPU.

• International Journal of Fluid Machinery and Systems
• /
• v.10 no.3
• /
• pp.254-263
• /
• 2017

The cavitation erosion remains an industrial issue for many applications. This paper deals with the cavitation intensity, which can be described as the fluid mechanical loading leading to cavitation damage. The estimation of this quantity is a challenging problem both in terms of modeling the cavitating flow and predicting the erosion due to cavitation. For this purpose, a numerical methodology was proposed to estimate cavitation intensity from 3D unsteady cavitating flow simulations. CFD calculations were carried out using Code_Saturne, which enables U-RANS equations resolution for a homogeneous fluid mixture using the Merkle's model, coupled to a $k-{\varepsilon}$ turbulence model with the Reboud's correction. A post-process cavitation intensity prediction model was developed based on pressure and void fraction derivatives. This model is applied on a flow around a hydrofoil using different physical (inlet velocities) and numerical (meshes and time steps) parameters. The article presents the cavitation intensity model as well as the comparison of this model with experimental results. The numerical predictions of cavitation damage are in good agreement with experimental results obtained by pitting test.

• Structural Engineering and Mechanics
• /
• v.57 no.3
• /
• pp.425-439
• /
• 2016

This paper deals with the behavior of fresh concrete that is under vibration using mass-spring model (MSM). To this end, behaviors of two different full scale precast concrete molds were investigated experimentally and theoretically. Experiments were performed under vibration with the use of a computer-based data acquisition system. Transducers were used to measure time-dependent lateral displacements at some points on mold while mold is empty and full of fresh concrete. Analytical modeling of molds used in experiments were prepared by three dimensional finite element method (3D FEM) using software. Modeling of full mold, using MSM, was made to solve the problem of dynamic interaction between fresh concrete and mold. Numerical displacement histories obtained from time history analysis were compared with experimental results. The comparisons show that the measured and computed results are compatible.

• Proceedings of the KIEE Conference
• /
• 1998.07a
• /
• pp.113-115
• /
• 1998

This paper presents the indirect boundary integral equation method(BIEM) to analyze 3D moving conductor problem. Instead of an artificial upwind algothm, the proposed method uses a fundamental Green's function which is a particular solution of diffusion equation. Therefore, this method yields a stable and accurate solution regardless of the Peclet number. The indirect BIEM is compared with 3D upwind FEM for a numerical model which has analytic solutions.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.06a
• /
• pp.203-208
• /
• 2007

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms for the travel time computation commonly used, however, may not give sufficiently precise results particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. We assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution.

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.222-227
• /
• 2008

This study is focused on an integrated numerical modeling enabling one to investigate the dynamic behavior and failure of 2-D textile composite and 3-D orthogonal woven composite structures weakened by micro-cracks and subjected to an impact load. The integrated numerical modeling is based on: I) determination of governing equations via a three-level hierarchy: micro-mechanical unit cell analysis, layer-wise analysis accounting for transverse strains and stresses, and structural analysis based on anisotropic plate layers, II) development of an efficient computational approach enabling one to perform transient response analyses of 2-D plain woven and 3-D orthogonal woven composite structures featuring the matrix cracking and exposed to time-dependent loads, III) determination of the structural characteristics of the textile-layered composites and their degraded features under various geometrical yarn shapes, and finally, IV) assessment of the implications of stiffness degradation on dynamic response to impact loads.