• Title/Summary/Keyword: Finite-difference method

Search Result 2,416, Processing Time 0.026 seconds

Numerical Simulation on Tidal Currents in a Bay - Application to Gamag Bay - (만내의 조류수치해석 - 가막만을 중심으로 -)

  • Lee, Gwan-Su;Lee, Yeong-Seok;Lee, Sam-No
    • Water for future
    • /
    • v.24 no.4
    • /
    • pp.41-48
    • /
    • 1991
  • This paper describes the characteristics of tidal currents in the Gamag Bay by using the two-dimensional nonlinear hydrodynamic euation. The basic equations are derived by Navier-Stokes momentum equation and continuity equation and its characteristics critically are reviewed, and they are analysed by the implicit finite difference method. The numerical model is constructed two-dimensional(depth-averaged) simple layer model. This paper investigates the stability of solution and convergence of solution in application of the method to Gamag Bay, and the reproducibility of the simulation is also discussed in comparison with the results of field survey. The following items are clarifide through the numerical investigation; i)the reproduc-ibility of tidal range and currents are quite acceptable, comparing with the results of model tests and field data, and ii) tidal cycle for convergent solution is four tidal cycle, also, iii) bottom friction is successfully represented as c=(1/n))$h^{1/6}$.

  • PDF

Comparisons between UPML and Liao's ABC in the FDTD method for 2D Cylindrical Coordinates (2D 원통형 좌표계를 위한 FDTD 방법에서 UPML과 Liao 흡수경계조건의 비교)

  • Hong, Ic-Pyo
    • Journal of the Korea Institute of Information and Communication Engineering
    • /
    • v.11 no.6
    • /
    • pp.1054-1061
    • /
    • 2007
  • In this paper, the comparison between UPML and Liao's absorbing boundary condition in the FDTD(Finite-Difference Time-Domain) method was performed for the analysis of the 2D cylindrical coordinate system. Generally, it is known as the absorbing characteristics of the UPML is bro than Liao's absorbing boundary condition in the 2D rectangular coordinate. The simulation results in this paper showed that Liao's original absorbing boundary condition is better than other two absorbing boundary conditions, Liao's modified condition and UPML. We concluded that more numerical, theoretical studies, simulations and verifications for various absorbing boundary conditions will be needed to get more accurate results for the design of useful 2D cylindrical microwave circuits.

Travel Time Calculation Using Mono-Chromatic Oneway Wave Equation (단일주파수 일방향파동방정식을 이용한 주시계산)

  • Shin, Chang-Soo;Shin, Sung-Ryul;Kim, Won-Sik;Ko, Seung-Won;Yoo, Hai-Soo
    • Geophysics and Geophysical Exploration
    • /
    • v.3 no.4
    • /
    • pp.119-124
    • /
    • 2000
  • A new fast algorithm for travel time calculation using mono-chromatic one-way wave equation was developed based on the delta function and the logarithms of the single frequency wavefield in the frequency domain. We found an empirical relation between grid spacing and frequency by trial and error method such that we can minimize travel time error. In comparison with other methods, travel time contours obtained by solving eikonal equation and the wave front edge of the snapshot by the finite difference modeling solution agree with our algorithm. Compared to the other two methods, this algorithm computes travel time of directly transmitted wave. We demonstrated our algorithm on migration so that we obtained good section showing good agreement with original model. our results show that this new algorithm is a faster travel time calculation method of the directly transmitted wave for imaging the subsurface and the transmission tomography.

  • PDF

Analysis on Design and Fabrication of High-diffraction-efficiency Multilayer Dielectric Gratings

  • Cho, Hyun-Ju;Lee, Kwang-Hyun;Kim, Sang-In;Lee, Jung-Hwan;Kim, Hyun-Tae;Kim, Won-Sik;Kim, Dong Hwan;Lee, Yong-Soo;Kim, Seoyoung;Kim, Tae Young;Hwangbo, Chang Kwon
    • Current Optics and Photonics
    • /
    • v.2 no.2
    • /
    • pp.125-133
    • /
    • 2018
  • We report an in-depth analysis of the design and fabrication of multilayer dielectric (MLD) diffraction gratings for spectral beam combining at a wavelength of 1055 nm. The design involves a near-Littrow grating and a modal analysis for high diffraction efficiency. A range of wavelengths, grating periods, and angles of incidence were examined for the near-Littrow grating, for the $0^{th}$ and $-1^{st}$ diffraction orders only. A modal method was then used to investigate the effect of the duty cycle on the effective indices of the grating modes, and the depth of the grating was determined for only the $-1^{st}$-order diffraction. The design parameters of the grating and the matching layer thickness between grating and MLD reflector were refined for high diffraction efficiency, using the finite-difference time-domain (FDTD) method. A high reflector was deposited by electron-beam evaporation, and a grating structure was fabricated by photolithography and reactive-ion etching. The diffraction efficiency and laser-induced damage threshold of the fabricated MLD diffraction gratings were measured, and the diffraction efficiency was compared with the design's value.

Numerical Study on Package Warpage as Structure Modeling Method of Materials for a PCB of Semiconductor Package (반도체 패키지용 PCB의 구조 모델링 방법에 따른 패키지의 warpage 수치적 연구)

  • Cho, Seunghyun;Ceon, Hyunchan
    • Journal of the Microelectronics and Packaging Society
    • /
    • v.25 no.4
    • /
    • pp.59-66
    • /
    • 2018
  • In this paper, we analyzed the usefulness of single-structured printed circuit board (PCB) modeling by using numerical analysis to model the PCB structure applied to a package for semiconductor purposes and applying modeling assuming a single structure. PCBs with circuit layer of 3rd and 4th were used for analysis. In addition, measurements were made on actual products to obtain material characteristics of a single structure PCB. The analysis results showed that if the PCB was modeled in a single structure compared to a multi-layered structure, the warpage analysis results resulting from modeling the PCB structure would increase and there would be a significant difference. In addition, as the circuit layer of the PCB increased, the mechanical properties of the PCB, the elastic coefficient and inertia moment of the PCB increased, decreasing the package's warpage.

Development of 2D Urban Inundation Analysis Model using Adaptive Mesh Refinement Method (메쉬 세분화 기법을 이용한 2차원 침수해석 모형의 개발)

  • Lee, Seung-soo
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2016.05a
    • /
    • pp.93-93
    • /
    • 2016
  • 최근 증가하고 있는 기후변화에 의해 설계빈도를 상회하는 강우의 발생빈도가 증가하고 있으며, 이로 인한 도시유역의 내수범람 피해가 증가하고 있다. 도시유역에서 발생하는 침수 피해의 경우 인적 물적 자원이 집중되어 있는 도시의 특성으로 인해 침수로 인한 직접적 피해 규모가 상당할 뿐만 아니라 침수 발생 후 세균 및 박테리아에 의해 발생하는 수인성 전염병의 유행 등과 같은 2차적 피해 또한 심각한 사회적 비용을 초래할 수 있어 도시유역의 침수 피해를 저감시키기 위한 대책이 절실히 요구되어지고 있다. 도시유역의 침수를 예방하기 위한 대책은 구조적 비구조적 대책으로 구분되어 질 수 있으며 구조적 대책의 경우 침수 피해 예방에 직접적인 효과를 낼 수 있다는 장점이 있으나 대규모 사업예산 및 사업 기간으로 인해 직접적 효과를 보기까지 상대적으로 긴 시간이 필요할 뿐만 아니라 사업 진행 중 대상지역 거주민들의 민원으로 인한 갈등 조정 등으로 인해 사업실행에 어려움을 겪고 있다. 이러한 측면에서 비구조적 대책의 일환인 수치해석을 통한 침수피해 재현 및 침수원인 파악을 통한 구조개선 제안은 구조적 대책의 단점을 보완할 수있는 좋은 대안이 될 수 있다. 도시유역의 경우 비도시유역과 대조적인 차이점으로는 높은 비율의 불투수층, 복잡한 지형, 다수의 인공 구조물 및 배수관망 시스템 등을 들 수 있으며, 침수해석 모형의 정확도를 높이기 위해서는 복잡한 지형의 효율적인 처리가 무엇보다 중요하다. 일반적으로 이용되는 2차원 침수해석 모형들은 직교구조 격자 또는 비구조 격자를 이용하여 지형을 묘사하고 있으며 DEM 자료를 직접 사용하는 직교구조 격자의 경우 지형 데이터 생성이 상대적으로 쉽다는 장점이 있으나 복잡한 지형을 표현하기 위해서는 불필요한 지역까지 높은 해상도를 이용해야 하며 이로 인하여 모의시간이 지나치게 길어지는 문제점이 발생한다. 비구조 격자의 경우 상대적으로 복잡한 도시 유역을 잘 묘사할 수 있다는 장점이 있으나 격자망 생성에 필요한 데이터가 많고 격자망 생성에 지나치게 많은 시간과 노력이 소요된다는 단점이 있다. 따라서 본 연구에서는 위에서 언급한 두 가지 방법의 장점만을 취할 수 있도록 메쉬 세분화 기법을 이용한 2차원 침수해석 모형을 개발 하여 복잡한 지형은 고해상도 메쉬를 이용하여 보다 자세히 묘사하고 상대적으로 복잡하지 않은 지형은 저해상도 메쉬를 이용하여 계산시간을 단축시킬 수 있도록 하였다. 수치해석 기법으로는 엇갈림 격자를 이용하는 Leap-Frog 기법과 유한차분 (Finite difference Method)기법을 이용하였다.

  • PDF

Prediction of Residual Settlement of Ground Improved by Vertical Drains Using the Elasto-Viscous Consolidation Model (I) - Verification of the Applicability of Theory - (탄-점성 압밀이론에 의한 버티칼 드레인 타설지반의 잔류침하 예측 (I) -이론의 적용성 검증)

  • Baek, Won-Jin
    • Journal of the Korean Geotechnical Society
    • /
    • v.23 no.4
    • /
    • pp.69-77
    • /
    • 2007
  • In this study, the consolidation behavior of clayey ground improved by vertical drain method was analyzed by the finite difference method based on the three-dimensional elasto-viscous consolidation theory, which can express the behavior of the secondary consolidation without considering the distinction of the normally consolidated and overconsolidated states. And the applicability of the elasto-viscous consolidation theory was discussed by comparing with the test results obtained from the model test of ground improved by vertical drain system. From these results, it is found that the amount of the settlement when the excess pore water pressure almost dissipated in the clay ground with vertical drains became smaller than that of the one-dimensional condition, and then the amount and rate of the residual settlement at secondary consolidation process became larger than those of the one-dimensional condition. finally, the effect of soil parameter on behavior of consolidation process was investigated by the results of a series of numerical analysis for the normalized and overconsoldiated states.

A Study on the Probabilistic Analysis Method Considering Spatial Variability of Soil Properties (지반의 공간적 변동성을 고려한 확률론적 해석기법에 관한 연구)

  • Cho, Sung-Eun;Park, Hyung-Choon
    • Journal of the Korean Geotechnical Society
    • /
    • v.24 no.8
    • /
    • pp.111-123
    • /
    • 2008
  • Geotechnical engineering problems are characterized by many sources of uncertainty. Some of these sources are connected to the uncertainties of soil properties involved in the analysis. In this paper, a numerical procedure for a probabilistic analysis that considers the spatial variability of soil properties is presented to study the response of spatially random soil. The approach integrates a commercial finite difference method and random field theory into the framework of a probabilistic analysis. Two-dimensional non-Gaussian random fields are generated based on a Karhunen-$Lo{\grave{e}}ve$ expansion in a fashion consistent with a specified marginal distribution function and an autocorrelation function. A Monte Carlo simulation is then used to determine the statistical response based on the random fields. A series of analyses were performed to study the effects of uncertainty due to the spatial heterogeneity on the settlement and bearing capacity of a rough strip footing. The simulations provide insight into the application of uncertainty treatment to the geotechnical problem and show the importance of the spatial variability of soil properties with regard to the outcome of a probabilistic assessment.

Development and validation of multiphysics PWR core simulator KANT

  • Taesuk Oh;Yunseok Jeong;Husam Khalefih;Yonghee Kim
    • Nuclear Engineering and Technology
    • /
    • v.55 no.6
    • /
    • pp.2230-2245
    • /
    • 2023
  • KANT (KAIST Advanced Nuclear Tachygraphy) is a PWR core simulator recently developed at Korea Advance Institute of Science and Technology, which solves three-dimensional steady-state and transient multigroup neutron diffusion equations under Cartesian geometries alongside the incorporation of thermal-hydraulics feedback effect for multi-physics calculation. It utilizes the standard Nodal Expansion Method (NEM) accelerated with various Coarse Mesh Finite Difference (CMFD) methods for neutronics calculation. For thermal-hydraulics (TH) calculation, a single-phase flow model and a one-dimensional cylindrical fuel rod heat conduction model are employed. The time-dependent neutronics and TH calculations are numerically solved through an implicit Euler scheme, where a detailed coupling strategy is presented in this paper alongside a description of nodal equivalence, macroscopic depletion, and pin power reconstruction. For validation of the steady, transient, and depletion calculation with pin power reconstruction capacity of KANT, solutions for various benchmark problems are presented. The IAEA 3-D PWR and 4-group KOEBERG problems were considered for the steady-state reactor benchmark problem. For transient calculations, LMW (Lagenbuch, Maurer and Werner) LWR and NEACRP 3-D PWR benchmarks were solved, where the latter problem includes thermal-hydraulics feedback. For macroscopic depletion with pin power reconstruction, a small PWR problem modified with KAIST benchmark model was solved. For validation of the multi-physics analysis capability of KANT concerning large-sized PWRs, the BEAVRS Cycle1 benchmark has been considered. It was found that KANT solutions are accurate and consistent compared to other published works.

Response of steel pipeline crossing strike-slip fault in clayey soils by nonlinear analysis method

  • Hadi Khanbabazadeh;Ahmet Can Mert
    • Geomechanics and Engineering
    • /
    • v.34 no.4
    • /
    • pp.409-424
    • /
    • 2023
  • Response of the pipeline crossing fault is considered as the large strain problem. Proper estimation of the pipeline response plays important role in mitigation studies. In this study, an advanced continuum modeling including material non-linearity in large strain deformations, hardening/softening soil behavior and soil-pipeline interaction is applied. Through the application of a fully nonlinear analysis based on an explicit finite difference method, the mechanics of the pipeline behavior and its interaction with soil under large strains is presented in more detail. To make the results useful in oil and gas engineering works, a continuous pipeline of two steel grades buried in two clayey soil types with four different crossing angles of 30°, 45°, 70° and 90° with respect to the pipeline axis have been considered. The results are presented as the fault movement corresponding to different damage limit states. It was seen that the maximum affected pipeline length is about 20 meters for the studied conditions. Also, the affected length around the fault cutting plane is asymmetric with about 35% and 65% at the fault moving and stationary block, respectively. Local buckling is the dominant damage state for greater crossing angle of 90° with the fault displacement varying from 0.4 m to 0.55 m. While the tensile strain limit is the main damage state at the crossing angles of 70° and 45°, the cross-sectional flattening limit becomes the main damage state at the smaller 30° crossing angles. Compared to the stiff clayey soil, the fault movement resulting 3% tensile strain limit reach up to 40% in soft clayey soil. Also, it was seen that the effect of the pipeline internal pressure reaches up to about 40% compared to non-pressurized condition for some cases.