• Title/Summary/Keyword: finite difference technique

Search Result 330, Processing Time 0.025 seconds

Application of 3D magnetotelluric investigation for geothermal exploration - Examples in Japan and Korea

  • Uchida Toshihiro;Song Yoonho;Mitsuhata Yuji;Lee Seong Kon
    • 한국지구물리탐사학회:학술대회논문집
    • /
    • 2003.11a
    • /
    • pp.390-397
    • /
    • 2003
  • A three-dimensional (3D) inversion technique has been developed for interpretation of magnetotelluric (MT) data. The inversion method is based on the linearized least-squares (Gauss-Newton) method with smoothness regularization. In addition to the underground 3D resistivity distribution, static shifts are also treated as unknown parameters in the inversion. The forward modeling is by the staggered-grid finite difference method. A Bayesian criterion ABle is applied to search the optimum trade-off among the minimization of the data misfit, model roughness and static shifts. The method has been applied to several MT datasets obtained at geothermal fields in Japan and other Asian countries. In this paper, two examples will be discussed: one is the data at the Ogiri geothermal area, southwestern Japan, and the other is at the Pohang low-enthalpy geothermal field, southeastern Korea. The inversion of the Ogiri data has been performed stably, resulting in a good fitting between the observed and computed apparent resistivities and phases. The recovered 3D resistivity structure is generally similar to the two-dimensional (2D) inversion models, although the deeper portion of the 3D model seems to be more realistic than that of the 2D model. The 3D model is also in a good agreement with the geological model of the geothermal reservoirs. 3D interpretation of the Pohang MT data is still preliminary. Although the fitting to the observed data is very good, the preliminary 3D model is not reliable enough because the station coverage is not sufficient for a 3D inversion.

  • PDF

Reactor core analysis through the SP3-ACMFD approach Part II: Transient solution

  • Mirzaee, Morteza Khosravi;Zolfaghari, A.;Minuchehr, A.
    • Nuclear Engineering and Technology
    • /
    • v.52 no.2
    • /
    • pp.230-237
    • /
    • 2020
  • In this part, an implicit time dependent solution is presented for the Boltzmann transport equation discretized by the analytic coarse mesh finite difference method (ACMFD) over the spatial domain as well as the simplified P3 (SP3) for the angular variable. In the first part of this work we proposed a SP3-ACMFD approach to solve the static eigenvalue equations which provide the initial conditions for temp dependent equations. Having solved the 3D multi-group SP3-ACMFD static equations, an implicit approach is resorted to ensure stability of time steps. An exponential behavior is assumed in transverse integrated equations to establish a relationship between flux moments and currents. Also, analytic integration is benefited for the time-dependent solution of precursor concentration equations. Finally, a multi-channel one-phase thermal hydraulic model is coupled to the proposed methodology. Transient equations are then solved at each step using the GMRES technique. To show the sufficiency of proposed transient SP3-ACMFD approximation for a full core analysis, a comparison is made using transport peers as the reference. To further demonstrate superiority, results are compared with a 3D multi-group transient diffusion solver developed as a byproduct of this work. Outcomes confirm that the idea can be considered as an economic interim approach which is superior to the diffusion approximation, and comparable with transport in results.

The Aging Effect of Dredging Clayey Soil on the Consolidation Characteristics (준설점성토의 압밀특성에 미치는 시간효과)

  • 김형주
    • Geotechnical Engineering
    • /
    • v.10 no.1
    • /
    • pp.71-82
    • /
    • 1994
  • According to the field measurement of dredging-reclaimed land, the actual self-weight consolidation settlement has been frequently reported to be less than the predicted values based on the laboratory tests results. The author estimates that one of the reasons is the reduction of the compressibility due to the sedimentation of the dredging material, Furthemore, the aging effect is ignored in the consolidation characteristics of the very low stress range as a pump dredging-reclaimed land. In this paper, a series of seepage consotidation tests has been carried out by applying the seepage force to the specimen prepared by sedimentation in consolidmeter in order to clarflfy the aging-effect on the compressibility of dredging clayey soil, Also, with a view to overcome unstable consolidation solution occurring in the case where the initial water content is higher as pumpdredging reclaimed land, the finite difference analysis technique using predictorforrector method is suggested that it gets good agreement with ezperimental results. Finally, the compressibility of the dredging clayey soil is depended on self-weight consolidation time.

  • PDF

A Back-Analysis of Tunnels in Multi-Layered Underground Structures (다층구조계내 터널 거동의 역해석)

  • 전병승;이상도;나경웅;김문겸
    • Tunnel and Underground Space
    • /
    • v.4 no.1
    • /
    • pp.17-23
    • /
    • 1994
  • This study consists of two procedures on back analysis and forward analysis which is a basic tool of the former. For a safe and economical construction of underground structures, it is required to identify the structural parameters and analyze the structural behavior as exactly as possible. In this paper, a boundary element method to analyze the behavior of multi-alyered underground structures is studied, in which body forces and initial stresses are considered. That is, each layer is discritized into subregions using infinite fundamental solutions, and terms of body forces and initial stresses are transformed into boundary integral where the applied direct integral method is used. And the system of equations containing body forces and initial stresses are considered. That is, each layer is discritized into subregions using infinite fundamental solutions, and terms of body forces and initial stresses are transformed into boundary integral where the applied direct integral method is used. And the system of equations containing body forces and initial stresses are composed, then the method to solve unknowns is used with applying compatibility and equilibrium conditions between interfaces. As well, the direct search method is applied in back analysis problems. By Powell's method as a technique to search unknown parameters, assuming displacements calculated from boundary element analysis as in-situ displacements, elastic moduli and initial stresses are presumed. As consequences of this study, the results of boundary element analysis of the behavior of multilayered structure considering body forces and initial stresses are agreed with those of finite element analysis. And results of back analysis of elastic moduli and initial stresses in each layers are agreed with exact values with a little difference. Therefore, it is known that this study can be efficiently applied for analyzing the behavior of underground structures including back analysis problems.

  • PDF

Dispersion of Nonconservative Contaminants Accidentally Released into Natural Streams (사고에 의하여 자연하천으로의 방류된 비보존성 오염물질의 종확산)

  • Jo, Seong-U;Jeon, Gyeong-Su
    • Journal of Korea Water Resources Association
    • /
    • v.34 no.4
    • /
    • pp.289-301
    • /
    • 2001
  • A fractional step finite difference model for the longitudinal dispersion of nonconservative pollutants is applied to the Nakdong River to simulate the phenol spill accident which occurred on March, 1971. Prior to the dispersion calculation, the flow conditions are simulated to provide inputs to the dispersion model. An unsteady flow model based on Preissmann's four-point scheme is used for this purpose. Sensitivities of the dispersion calculation to empirical equations for dispersion coefficient and to the first-order decay coefficient are analyzed. The time to peak concentration at a downstream location is significantly different depending on the formula for the dispersion coefficient. Although the decay coefficient does not affect the shape of the temporal concentration distribution, the concentration values depend on the decay coefficient verb significantly. An optimization technique is used to calibrate the dispersion model as well as the flow model. The time to the peak concentration is simulated for major positions of water intake along the Nakdong River.

  • PDF

Comparative Study on Sloshing Impact Flows between PIV and CFD (슬로싱 충격현상 해석을 위한 모형실험과 수치해석 적용에 관한 비교 연구: PIV vs. CFD)

  • Yang, Kyung-Kyu;Kim, Jieung;Kim, Sang-Yeob;Kim, Yonghwan
    • Journal of Ocean Engineering and Technology
    • /
    • v.29 no.2
    • /
    • pp.154-162
    • /
    • 2015
  • In this study, experimental and numerical methods were applied to observe sloshing impact phenomena. A two-dimensional rectangular tank filled with water and air was considered with a specific excitation condition that induced a hydrodynamic impact without an air pocket at the top corner of the tank. High-speed cameras and a pressure measurement system were synchronized, and a particle image velocimetry (PIV) technique was applied to measure the velocity field and corresponding pressure. The experimental condition was implemented in a numerical computation to solve incompressible two-phase flows using a Cartesian-grid method. The discretized solution was obtained using the finite difference and constraint-interpolation-profile (CIP) methods, which adopt a fractional step scheme for coupling the pressure and velocity. The tangent of the hyperbola for interface capturing (THINC) scheme was used with the weighed line interface calculation (WLIC) method to capture the interface between the air and water. The calculated impact pressures and velocity fields were compared with experimental data, and the relationship between the local velocity and pressure was investigated based on the computational results.

COMPARISONS OF PARALLEL PRECONDITIONERS FOR THE COMPUTATION OF SMALLEST GENERALIZED EIGENVALUE

  • Ma, Sang-Back;Jang, Ho-Jong;Cho, Jae-Young
    • Journal of applied mathematics & informatics
    • /
    • v.11 no.1_2
    • /
    • pp.305-316
    • /
    • 2003
  • Recently, an iterative algorithm for finding the interior eigenvalues of a definite matrix by CG-type method has been proposed. This method compares to the inverse power method. The given matrices A, and B are assumed to be large and sparse, and SPD( Symmetric Positive Definite) The CG scheme for the optimization of the Rayleigh quotient has been proven a very attractive and promising technique for large sparse eigenproblems for smallest eigenvalue. Also, it is very amenable to parallel computations, like the CG method for the linear systems. A proper choice of the preconditioner significantly improves the convergence of the CG scheme. But for parallel computations we need to find an efficient parallel preconditioner. Our candidates we ILU(0) in the wave-front order, ILU(0) in the multi-coloring order, Point-SSOR(Symmetric Successive Overrelaxation), and Multi-Color Block SSOR preconditioner. Wavefront order is a simple way to increase parallelism in the natural order, and Multi-coloring realizes a parallelism of order(N), where N is the order of the matrix. Another choice is the Multi-Color Block SSOR(Symmetric Successive OverRelaxation) preconditioning. Block SSOR is a symmetric preconditioner which is expected to minimize the interprocessor communication due to the blocking. We implemented the results on the CRAY-T3E with 128 nodes. The MPI (Message Passing Interface) library was adopted for the interprocessor communications. The test problem was drawn from the discretizations of partial differential equations by finite difference methods. The results show that for small number of processors Multi-Color ILU(0) has the best performance, while for large number of processors Multi-Color Block SSOR performs the best.

A New Evaluation Method for Interfacial Properties of Composites using the Gradual Multi-Fiber Fragmentation Test (단계적 다섬유 Fragmentation 시험법을 이용한 복합재료의 계면적 특성에 대한 새로운 평가방법)

  • ;;Koichi Goda
    • Composites Research
    • /
    • v.12 no.2
    • /
    • pp.10-25
    • /
    • 1999
  • A new evaluation method for the interfacial properties of fibrous composites based on a fragmentation technique is proposed by using the gradual multi-fiber composite, in which the inter-fiber spacing is gradually changed. The results showed that as the inter-fiber distance increased, the aspect ratio of broken fibers decreased while the interfacial shear strength between the fiber and matrix increased. When the reciprocal of the inter-fiber destance was taken for the above relations, both the aspect ratio and interfacial shear strength showed a saturated value. This means that the gradual multi-fiber composite indicates an upper bound in aspect ratio and an upper bound in interfacial shear strength. It was concluded that this fragmentation test could be a new method for composite evaluation, since reducing a difference between these two bounds is effective for composite strengthening. In addition an elastoplastic finite element analysis was carried out to relate the above results with fiber stress a distribution around fiber breaks. It was proved that the bound obtained in the gradual multi-fiber composite test is closely related to stress concentrations caused by a group of multi-fiber breaks.

  • PDF

Analysis of A New Crossbar Embedded Structure for Improved Attenuation Characteristics on the Various Lossy Media (다양한 손실매질내의 손실특성 개선을 위한 새로운 크로스바 구조의 해석)

  • Kim, Yoon-Suk
    • Journal of the Institute of Electronics Engineers of Korea TC
    • /
    • v.43 no.12 s.354
    • /
    • pp.83-88
    • /
    • 2006
  • In this paper, we propose a new cross bar embedded structure for improvement of attenuation characteristics along the different lossy media. A general characterization procedure based on the extraction of the characteristic impedance and propagation constant for analyzing a single MIS(Metal-Insulator-Semiconductor) transmission line used and an analysis for a new substrate shielding MIS structure consisting of grounded crossbars at the interface between Si and Sio2 layer using the Finite-Difference Time-Domain(FDTD) technique is used. In order to reduce the substrate effects on the transmission line characteristics, a shielding structure consisting of grounded cross bar lines over time-domain signal has been examined. The extracted, distributed frequency-dependent transmission line parameters as well as the line voltages and currents, and also corresponding equivalent circuit parameters have been examined as function of frequency. It is shown that the quality factor of the transmission line can be improved without significant changes in the characteristic impedance and effective dielectric constant.

Simulation of the Combined Effects of Dipole Emitter Orientation, Mie Scatterers, and Pillow Lenses on the Outcoupling Efficiency of an OLED (쌍극자 광원의 진동방향, Mie 산란자, 그리고 Pillow 렌즈가 OLED의 광추출효율에 미치는 영향에 대한 시뮬레이션 연구)

  • Lee, Ju Seob;Lee, Jong Wan;Park, Jaehoon;Ko, Jae-Hyeon
    • Korean Journal of Optics and Photonics
    • /
    • v.25 no.4
    • /
    • pp.193-199
    • /
    • 2014
  • The net effect of the emitter orientation, Mie scatters, and pillow lenses on the outcoupling efficiency (OCE) of a bottom-emitting OLED having an internal photonic crystal layer was investigated by a combined optical simulation based on the finite-difference time-domain method (FDTD) and the ray-tracing technique. The simulation showed that when the emitter orientation was horizontal with respect to the OLED surface, the OCE could be increased by 54% when a photonic crystal layer was employed, while it could be improved by 86% under optimized conditions of Mie scatters and pillow lenses applied to the glass substrate. The peculiar intensity distribution of the OLED, caused by the periodic lattice structure of the photonic crystal layer, could be ameliorated by inserting Mie scatters into the glass substrate. This study suggests that conventional outcoupling structures combined with control of the emitter orientation could improve the OCE substantially.