• Title/Summary/Keyword: 직접 수치해석

Search Result 555, Processing Time 0.031 seconds

An Introduction to Kinetic Monte Carlo Methods for Nano-scale Diffusion Process Modeling (나노 스케일 확산 공정 모사를 위한 동력학적 몬테칼로 소개)

  • Hwang, Chi-Ok;Seo, Ji-Hyun;Kwon, Oh-Seob;Kim, Ki-Dong;Won, Tae-Young
    • Journal of the Institute of Electronics Engineers of Korea SD
    • /
    • v.41 no.6
    • /
    • pp.25-31
    • /
    • 2004
  • In this paper, we introduce kinetic Monte Carlo (kMC) methods for simulating diffusion process in nano-scale device fabrication. At first, we review kMC theory and backgrounds and give a simple point defect diffusion process modeling in thermal annealing after ion (electron) implantation into Si crystalline substrate to help understand kinetic Monte Carlo methods. kMC is a kind of Monte Carlo but can simulate time evolution of diffusion process through Poisson probabilistic process. In kMC diffusion process, instead of. solving differential reaction-diffusion equations via conventional finite difference or element methods, it is based on a series of chemical reaction (between atoms and/or defects) or diffusion events according to event rates of all possible events. Every event has its own event rate and time evolution of semiconductor diffusion process is directly simulated. Those event rates can be derived either directly from molecular dynamics (MD) or first-principles (ab-initio) calculations, or from experimental data.

Bending Effect of Laminated Plates with a Circular Hole Repaired by Single-Sided Patch Based on p-Convergent Full Layerwise Model (p-수렴 완전층별모델에 의한 일면패치로 보강된 원공 적층판의 휨효과)

  • Woo, Kwang-Sung;Yang, Seung-Ho;Ahn, Jae-Seok;Shin, Young-Sik
    • Journal of the Computational Structural Engineering Institute of Korea
    • /
    • v.22 no.5
    • /
    • pp.463-474
    • /
    • 2009
  • Double symmetric patch repair of existing structures always causes membrane action only, however, in many cases this technique is not practical. On the other hand, the bending stiffness of the patch and the skin increases as tensile loading is increased and affects the bending deformation significantly in the case of single-sided patch repair. In this study, the p-convergent full layerwise model has been proposed to determine the stress concentration factor in the vicinity of a circular hole as well as across the thickness of plates with single-sided patch repair. In assumed displacement field, the strain-displacement relations and 3-D constitutive equations of a layer are obtained by the combination of 2-D and 3-D hierarchical shape functions. The transfinite mapping technique has been used to represent a circular boundary and Gauss-Lobatto numerical integration is implemented in order to directly obtain stresses occurred at the nodal points of each layer without other extrapolation techniques. The accuracy and simplicity of the present model are verified with comparison of the previous results in literatures using experiment and conventional 3-D finite element. Also, the bending effect has been investigated with various patch types like square, circular and annular shape.

용담댐 건설로 인한 금강 상류의 하천환경변화 분석 II - 수리변화분석 -

  • 정승권;정동양
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2004.05b
    • /
    • pp.971-976
    • /
    • 2004
  • 하천유지를 위한 하천내 최소유량의 충족은 사실상 환경문제에 제공되는 것으로서 기본적으로 하천변 수생식물의 생태계 서식처를 유지하거나 복원하기 위해 제공된다. 또한 하천유지유량은 하천의 고유종 동${\cdot}$식물을 보전하고, 아름다운 하천에 대한 심미적 질을 유지하며 여가활동 및 상업적 어획고를 최대로 늘릴 수 있거나 과학적 또는 문화적 관심사의 특성을 보호하는 기능을 위한 목적으로도 제공된다. 그러므로 생태학적인 하천유지유량 차원에서 무엇보다 중요한 것은 충분한 물을 확보하고 쾌적한 하천을 유지할 수 있는 맑은 물을 하천에 흐르도록 하는 것이다. 이와 더불어 하천의 생태학적인 측면에서는 서식지 보호를 고려한 계획이 뒷받침되어야 할 것이다. 기존의 개발은 하천정비, 골재채취, 수중 보와 하구 둑 등 하천시설물 건설, 그리고 댐 및 교량건설 등 인위적인 요인에 의하여 하천환경이 심각하게 변화는 경우가 많다. 특히 댐 건설로 인해 하류지역에 수리학적 변화와 하상 변동 등이 발생하고, 이로 인하여 생태계의 변화를 초래하기도 한다. 따라서 본 연구에서는 금강 최상류에 건설된 용담다목적댐에 의해 인위적으로 하천유입량이 감소하고 이로 인해 발생할 수 있는 수리학적 변화를 모의하여 댐 건설이 금강 상류에 미치는 영향을 분석하고자 한다. 이를 위해 2차원 유한요소모형을 이용하였으며, 댐 건설전 자연유출량 및 댐 건설후 하천유지유량에 의한 하천의 건조면적, 습지면적을 산정하였다. 또한 하천의 한계소류력과 하도를 형성하는 재료에 따라 바닥면 전단응력을 산정하여 그 차(差)를 활용하여 유량에 따른 유사이송현황을 모의하였다. 모의된 결과를 visual하게 표현하고자 GIS 기법을 이용하여 주요구간에 대한 유사이송현황을 도시하였다.유역의 고도차를 이용한 흐름특성 분석을 위해 수치고도자료를 이용하여 유역흐름특성을 분석할 수 있는 TOPAZ(Topographic PArameteri-Zation) 프로그램을 이용하였다. TOPAZ 프로그램을 통해 분석된 각 격자별 분포형 수문 매개변수는 적합한 관계식을 통해 분포형 유출량을 모의하는데 적용된다.다 정확한 유입량 예측이 가능할 것으로 사료된다.이 작은 오차를 발생하였으며, 전체적으로 퍼프 모형이 입자모형보다는 훨씬 적은 수의 계산을 통해서도 작은 오차를 나타낼 수 있다는 것을 알 수 있었다. 그러나 Gaussian 분포를 갖는 퍼프모형은 전단흐름에서의 긴 유선형 농도분포를 모의할 수 없었고, 이에 관한 오차는 전단계수가 증가함에 따라 비선형적으로 증가하였다. 향후, 보다 다양한 흐름영역에서 장${\cdot}$단점 분석 및 오차해석을 수행한 후에 각각의 Lagrangian 모형의 장점만을 갖는 모형결합 방법을 제시할 수 있을 것으로 판단된다.mm/$m^{2}$로 감소한 소견을 보였다. 승모판 성형술은 전 승모판엽 탈출증이 있는 두 환아에서 동시에 시행하였다. 수술 후 1년 내 시행한 심초음파에서 모든 환아에서 단지 경등도 이하의 승모판 폐쇄 부전 소견을 보였다. 수술 후 조기 사망은 없었으며, 합병증으로는 유미흉이 한 명에서 있었다. 술 후 10개월째 허혈성 확장성 심근증이 호전되지 않아 Dor 술식을 시행한 후 사망한 예를 제외한 나머지 6명은 특이 증상 없이 정상 생활 중이다 결론: 좌관상동맥 페동맥이상 기시증은 드물기는 하나, 영유아기에 심근경색 및 허혈성 심근증 또는 선천성 승모판 폐쇄 부전등을 초래하는 심각한 선천성 심질환이다. 그러나 진단 즉시 직접 좌관상동맥-대동맥 이식술로 수술

  • PDF

Seismic First Arrival Time Computation in 3D Inhomogeneous Tilted Transversely Isotropic Media (3차원 불균질 횡등방성 매질에 대한 탄성파 초동 주시 모델링)

  • Jeong, Chang-Ho;Suh, Jung-Hee
    • Geophysics and Geophysical Exploration
    • /
    • v.9 no.3
    • /
    • pp.241-249
    • /
    • 2006
  • 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 commonly used, however, may not give sufficiently precise computational results of traveltime data 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. Considering the complex geology of Korea, 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. The performance of the algorithm developed in this study is demonstrated by the comparison of the analytic and numerical solutions for the homogeneous anisotropic earth as well as through the numerical experiment for the two layer model whose anisotropic properties are greatly different each other. We expect that the developed modeling algorithm can be used in the development of processing and inversion schemes of seismic data acquired in strongly anisotropic environment, such as migration, velocity analysis, cross-well tomography and so on.

강지진동 분석의 최적화를 위한 고려요소

  • 이석태;조봉곤;이정모;조영삼
    • Proceedings of the International Union of Geodesy And Geophysics Korea Journal of Geophysical Research Conference
    • /
    • 2003.05a
    • /
    • pp.17-17
    • /
    • 2003
  • 한반도에 있어서의 지진의 영향을 분석하기 위해서는 강지진동 연구가 필수적이다. 강지진동 자료가 부족한 한반도의 특성상 모사를 통해 연구하고 있다. 강지진동 분석을 하기 위해서는 되도록 노이즈가 포함되어 있지 않은 지진파자료를 선택하여 그 지진자료의 스펙트럼 분석을 통해 감쇠상수 k, Q 등을 구한다. 이러한 감쇠상수 값을 통해 한반도의 진동 특성을 이해할 수 있다. 그러나 감쇠상수를 구하는 과정에서 감쇠상수 분석에 사용된 지진자료에 노이즈가 더해졌을 경우, 어떤 형태로 스펙트럼 영역에 영향을 미치고, 감쇠상수에는 어떤 영향을 미치는 지를 연구하여 노이즈효과를 제거할 수 있는 최적화된 분석에 관한 연구가 선행되어야 한다고 본다. 따라서 이번 연구에서는 강지진동 모사프로그램을 가지고 노이즈효과를 적용하면서 감쇠상수에 노이즈가 어떤 영향을 미치는 지에 대한 수치 해석적 연구를 실시하였다. 합성지진파에 이 합성지진파와 전혀 다른 주파수 형태를 보이는 노이즈를 강도를 달리하면서 합성해 본 결과, 노이즈효과를 고려할 수 있는 몇 가지 요소가 있음을 알 수 있었다. 감쇠상수 k값을 강지진동 모사프로그램으로부터 값을 달리하며 합성해 본 결과 노이즈효과를 보이는 것을 알 수 있었으며, 감쇠상수 k를 선형회귀를 통해 $k_{s}$$k_{q}$를 구할 때의 적용 주파수 범위를 변화시켰을 때도 일정한 양상의 노이즈 효과를 보였다. 또 지진자료와 노이즈를 중첩시킨 지진파 시계열 자료의 정부분만을 감쇠상수 k를 구하는 선형회귀에 이용했을 경우에도 노이즈 효과를 보였다. 또한 계산되어 나온 감쇠상수 값으로부터 특정지역의 지반운동의 특성을 이해할 수 있는 스펙트럼 가속도, 최대 가속도, 및 최대속도 값에 따른 감쇠식을 구하였다. 이것을 한반도와 같은 판 내부 환경인 ENA 값과 비교하였으며 기존의 연구와도 비교하였다.심으로부터 지오이드까지의 거리, 지오이드로부터 지표까지의 거리를 정의해주었으며, 각 격자점의 수직구조를 정의하기 위해 깊이에 따른 각 매질의 밀도, P파의 속도, S파의 속도, P파에 대한 Q값, S파에 대한 Q값을 정의 해주었다. S파의 속도를 구하기 위해서 지구 내부 물질을 포아송 매질이라는 가정 하에, 관계식을 $Vp{\;}={\;}SQRT(3){\;}{\times}{\;}Vs$ 이용하였다. 획득한 모델치들을 이용해 동해와 동해 인근 지역에 대한 초기모델을 구축하였다. 약 1 × 10/sup 6/ e/sup -//sec·n㎡ 의 전자선량에 해당되며 이를 기준으로 각각의 illumination angle에 대한 임계전자선량을 평가할 수 있었다. 실질적으로 Cibbsite와 같은 무기수화물의 직접가열실험 시 전자빔 조사에 의해 야기되는 상전이 영향을 배제하고 실험을 수행하려면 illumination angle 0.2mrad (Dose rate : 8000 e/sup -//sec·n㎡)이하로 관찰하고 기록되어야 함을 본 자료로부터 알 수 있었다.운동횟수에 의한 영향으로써 운동시간을 1일 6시간으로 설정하여, 운동횟수를 결정하기 위하여 오전, 오후에 각 3시간씩 운동시키는 방법과 오전부터 6시간동안 운동시키는 두 방법을 이용하여 품질을 비교하였다. 각 조건에 따라 운동시킨 참돔의 수분함량을 나타낸 것으로, 2회(오전 3시간, 오후 3시간)에 나누어서 운동시키기 위한 육의 수분함량은 73.37±2.02%를 나타냈으며, 1회(6시간 운동)운동시키기 위한 육은 71.74±1.66%을 나타내었다. 각각의 운동조건에서 양식된 참돔은 사육초기에는 큰 변화가 없었으나, 사육 5일 이후에는 수분함량이 증가하여 15일에는 76.40±0.14, 75.62±0.98%의 수분함량을 2회와 1회 운동시킨 참돔의 육에서 각각 나타났다. 운동횟수에 따른 지

  • PDF

Modeling and analysis of dynamic heat transfer in the cable penetration fire stop system by using a new hybrid algorithm (새로운 혼합알고리즘을 이용한 CPFS 내에서의 일어나는 동적 열전달의 수식화 및 해석)

  • Yoon En Sup;Yun Jongpil;Kwon Seong-Pil
    • Journal of the Korean Institute of Gas
    • /
    • v.7 no.4 s.21
    • /
    • pp.44-52
    • /
    • 2003
  • In this work dynamic heat transfer in a CPFS (cable penetration fire stop) system built in the firewall of nuclear power plants is three-dimensionally investigated to develop a test-simulator that can be used to verify effectiveness of the sealant. Dynamic heat transfer in the fire stop system is formulated in a parabolic PDE (partial differential equation) subjected to a set of initial and boundary conditions. First, the PDE model is divided into two parts; one corresponding to heat transfer in the axial direction and the other corresponding to heat transfer on the vertical planes. The first PDE is converted to a series of ODEs (ordinary differential equations) at finite discrete axial points for applying the numerical method of SOR (successive over-relaxation) to the problem. The ODEs are solved by using an ODE solver In such manner, the axial heat flux can be calculated at least at the finite discrete points. After that, all the planes are separated into finite elements, where the time and spatial functions are assumed to be of orthogonal collocation state at each element. The initial condition of each finite element can be obtained from the above solution. The heat fluxes on the vertical planes are calculated by the Galerkin FEM (finite element method). The CPFS system was modeled, simulated, and analyzed here. The simulation results were illustrated in three-dimensional graphics. Through simulation, it was shown clearly that the temperature distribution was influenced very much by the number, position, and temperature of the cable stream, and that dynamic heat transfer through the cable stream was one of the most dominant factors, and that the feature of heat conduction could be understood as an unsteady-state process.

  • PDF

A Study on the Conservation of Biodiversity by the Ecological Economic Numerical Model (생태경제수치모형에 의한 생물다양성 보존에 관한 연구)

  • Kim, Byung-Nam
    • The Journal of the Korea Contents Association
    • /
    • v.22 no.6
    • /
    • pp.629-637
    • /
    • 2022
  • It is at risk of depletion of biodiversity due to indiscriminate overfishing of ecosystems and destruction of habitats. Intensive fertilizers or development of related facilities to increase agricultural production in poor indigenous areas devastate the soil. Preservation of biodiversity is now emerging as an important issue of global human coexistence. After the Post-2020 GBF Declaration, all governance in agricultural development in indigenous agricultural areas should be supported and promoted as biodiversity conservation measures. A compromise plan to reduce ecosystem development and biodiversity loss can help establish public governance policies. In this paper, a viability kernel used for viable control feedback analysis is introduced to solve conflicting economic and ecological problems in ecosystem conservation, and a mathematical model on biodiversity conservation by the viability kernel is examined. Because all species in the ecosystem are interdependent, if the balance is broken, biodiversity is depleted, which is irreversible and eventually leads to extinction. For sustainable use and harmony of biological resources, a lot of policy consideration is required, such as creative governance that can efficiently protect all species. Subsidies or tax incentives have a direct impact on biodiversity conservation. The recovery of species in a state of decreasing biodiversity can be said to be of great economic value. Biodiversity will allow indigenous producers to be proud of their unique traditional knowledge and have a positive impact on local tourism, thereby enhancing regional identity and greatly contributing to the survival and prosperity of mankind.

Numerical Study on the Stress-distribution Ratio of Grouting Pile for Reinforced Ground (지반보강용 그라우팅 말뚝의 응력분담비에 대한 수치해석적 연구)

  • Yi, Gyeong-Ju;Lee, Joon-Kyu;Zhang Weiwei;Song, Ki-Il
    • Journal of the Korean Geotechnical Society
    • /
    • v.39 no.2
    • /
    • pp.19-30
    • /
    • 2023
  • Underground structures, such as compacted sand piles applied as soft ground countermeasures, are analyzed for settlement and stability by the composite ground design method. The basic principle of the composite ground design method is the arching effect. The reinforcing effect of the pile is evaluated as the stress-distribution ratio. When applying grouting piles with elastic properties using the ground reinforcement method, the existing stress-distribution ratio was only considered when the pile was installed. This study shows that the method of applying the stress-distribution ratio applied in previous studies should be changed when the ground reinforcement pile is installed at an arbitrary location in the ground without raising it to the ground surface. When high strength jet routing is applied, the stress-distribution ratio (n) to the in-situ ground generally ranges from 30 to 50. However, if the pile is located far from the surface and the depth goes down to the boundary depth of the stress sphere, the stress-distribution effect rapidly decreases, and the stress-distribution ratio converges to 1.5.

Numerical investigation on cavitation and non-cavitation flow noise on pumpjet propulsion (펌프젯 추진기의 공동 비공동 유동소음에 대한 수치적 연구)

  • Garam Ku;Cheolung Cheong;Hanshin Seol;Hongseok Jeong
    • The Journal of the Acoustical Society of Korea
    • /
    • v.42 no.3
    • /
    • pp.250-261
    • /
    • 2023
  • In this study, the noise contributions by the duct, stator and rotor, which are the propulsor components, are evaluated to identify the flow noise source in cavitation and non-cavitation conditions on pumpjet propulsion and the noise levels in both conditions are compared. The unsteady incompressible Reynolds averaged Navier-Stokes (RANS) equation based on the homogeneous mixture assumption is applied on the suboff submarine hull and pumpjet propeller in the cavitation tunnel, and the Volume of Fluid (VOF) method and Schnerr-Sauer cavitation model are used to describe the two-phase flow. Based on the flow simulation results, the acoustic analogy formulated by Ffowcs Williams and Hawkings (FW-H) equation is applied to predict the underwater radiated noise. The noise contributions are evaluated by using the three types of impermeable integral surface on the duct, stator and rotor, and the two types of permeable integral surface surrounding the propulsor. As a result of noise prediction, the contribution by the stator is insignificant, but it affects the generation of flow noise source due to flow separation in the duct and rotor, and the noise is predominantly radiated into the upward and right where the flow separations are. Also, the noise is radiated into the thrust direction due to pressure fluctuation between suction and pressure sides on the rotor blades, and the it can be seen that the cavitation effect into the noise can be considered through the permeable integral surface.

Development of Steel Composite Cable Stayed Bridge Weigh-in-Motion System using Artificial Neural Network (인공신경망을 이용한 강합성 사장교 차량하중분석시스템 개발)

  • Park, Min-Seok;Jo, Byung-Wan;Lee, Jungwhee;Kim, Sungkon
    • KSCE Journal of Civil and Environmental Engineering Research
    • /
    • v.28 no.6A
    • /
    • pp.799-808
    • /
    • 2008
  • The analysis of vehicular loads reflecting the domestic traffic circumstances is necessary for the development of adequate design live load models in the analysis and design of cable-supported bridges or the development of fatigue load models to predict the remaining lifespan of the bridges. This study intends to develop an ANN(artificial neural network)-based Bridge WIM system and Influence line-based Bridge WIM system for obtaining information concerning the loads conditions of vehicles crossing bridge structures by exploiting the signals measured by strain gauges installed at the bottom surface of the bridge superstructure. This study relies on experimental data corresponding to the travelling of hundreds of random vehicles rather than on theoretical data generated through numerical simulations to secure data sets for the training and test of the ANN. In addition, data acquired from 3 types of vehicles weighed statically at measurement station and then crossing the bridge repeatedly are also exploited to examine the accuracy of the trained ANN. The results obtained through the proposed ANN-based analysis method, the influence line analysis method considering the local behavior of the bridge are compared for an example cable-stayed bridge. In view of the results related to the cable-stayed bridge, the cross beam ANN analysis method appears to provide more remarkable load analysis results than the cross beam influence line method.