• 제목/요약/키워드: Governing Equation

검색결과 1,286건 처리시간 0.023초

Markov Envelope를 이용한 지진동의 위상차 확률분포와 전파지연시간의 추정 (Inference of the Probability Distribution of Phase Difference and the Path Duration of Ground Motion from Markov Envelope)

  • 최항;윤병익
    • 한국지진공학회논문집
    • /
    • 제26권5호
    • /
    • pp.191-202
    • /
    • 2022
  • Markov envelope as a theoretical solution of the parabolic wave equation with Markov approximation for the von Kármán type random medium is studied and approximated with the convolution of two probability density functions (pdf) of normal and gamma distributions considering the previous studies on the applications of Radiative Transfer Theory (RTT) and the analysis results of earthquake records. Through the approximation with gamma pdf, the constant shape parameter of 2 was determined regardless of the source distance ro. This finding means that the scattering process has the property of an inhomogeneous single-scattering Poisson process, unlike the previous studies, which resulted in a homogeneous multiple-scattering Poisson process. Approximated Markov envelope can be treated as the normalized mean square (MS) envelope for ground acceleration because of the flat source Fourier spectrum. Based on such characteristics, the path duration is estimated from the approximated MS envelope and compared to the empirical formula derived by Boore and Thompson. The results clearly show that the path duration increases proportionately to ro1/2-ro2, and the peak value of the RMS envelope is attenuated by exp (-0.0033ro), excluding the geometrical attenuation. The attenuation slope for ro≤100 km is quite similar to that of effective attenuation for shallow crustal earthquakes, and it may be difficult to distinguish the contribution of intrinsic attenuation from effective attenuation. Slowly varying dispersive delay, also called the medium effect, represented by regular pdf, governs the path duration for the source distance shorter than 100 km. Moreover, the diffraction term, also called the distance effect because of scattering, fully controls the path duration beyond the source distance of 300 km and has a steep gradient compared to the medium effect. Source distance 100-300 km is a transition range of the path duration governing effect from random medium to distance. This means that the scattering may not be the prime cause of peak attenuation and envelope broadening for the source distance of less than 200 km. Furthermore, it is also shown that normal distribution is appropriate for the probability distribution of phase difference, as asserted in the previous studies.

유한 요소 해석을 이용한 고주파 간 종양 절제술의 입력 파형 최적화를 위한 연구 (A Study For Optimizing Input Waveforms In Radiofrequency Liver Tumor Ablation Using Finite Element Analysis)

  • 임도형;남궁범석;이태우;최진승;탁계래;김한성
    • 대한의용생체공학회:의공학회지
    • /
    • 제28권2호
    • /
    • pp.235-243
    • /
    • 2007
  • Hepatocellular carcinoma is significant worldwide public health problem with an estimated annually mortality of 1,000,000 people. Radiofrequency (RF) ablation is an interventional technique that in recent years has come to be used for treatment of the hepatocellualr carcinoma, by destructing tumor tissues in high temperatures. Numerous studies have been attempted to prove excellence of RF ablation and to improve its efficiency by various methods. However, the attempts are sometimes paradox to advantages of a minimum invasive characteristic and an operative simplicity in RF ablation. The aim of the current study is, therefore, to suggest an improved RF ablation technique by identifying an optimum RF pattern, which is one of important factors capable of controlling the extent of high temperature region in lossless of the advantages of RF ablation. Three-dimensional finite element (FE) model was developed and validated comparing with the results reported by literature. Four representative Rf patterns (sine, square, exponential, and simulated RF waves), which were corresponding to currents fed during simulated RF ablation, were investigated. Following parameters for each RF pattern were analyzed to identify which is the most optimum in eliminating effectively tumor tissues. 1) maximum temperature, 2) a degree of alteration of maximum temperature in a constant time range (30-40 second), 3) a domain of temperature over $47^{\circ}C$ isothermal temperature (IT), and 4) a domain inducing over 63% cell damage. Here, heat transfer characteristics within the tissues were determined by Bioheat Governing Equation. Developed FE model showed 90-95% accuracy approximately in prediction of maximum temperature and domain of interests achieved during RF ablation. Maximum temperatures for sine, square, exponential, and simulated RF waves were $69.0^{\circ}C,\;66.9^{\circ}C,\;65.4^{\circ}C,\;and\;51.8^{\circ}C$, respectively. While the maximum temperatures were decreased in the constant time range, average time intervals for sine, square, exponential, and simulated RE waves were $0.49{\pm}0.14,\;1.00{\pm}0.00,\;1.65{\pm}0.02,\;and\;1.66{\pm}0.02$ seconds, respectively. Average magnitudes of the decreased maximum temperatures in the time range were $0.45{\pm}0.15^{\circ}C$ for sine wave, $1.93{\pm}0.02^{\circ}C$ for square wave, $2.94{\pm}0.05^{\circ}C$ for exponential wave, and $1.53{\pm}0.06^{\circ}C$ for simulated RF wave. Volumes of temperature domain over $47^{\circ}C$ IT for sine, square, exponential, and simulated RF waves were 1480mm3, 1440mm3, 1380mm3, and 395mm3, respectively. Volumes inducing over 63% cell damage for sine, square, exponential, and simulated RF waves were 114mm3, 62mm3, 17mm3, and 0mm3, respectively. These results support that applying sine wave during RF ablation may be generally the most optimum in destructing effectively tumor tissues, compared with other RF patterns.

정형 사각 격자 기반의 2차원 지표면 침수해석 모형 개발 및 평가 (Development and evaluation of a 2-dimensional land surface flood analysis model using uniform square grid)

  • 최윤석;김주훈;최천규;김경탁
    • 한국수자원학회논문집
    • /
    • 제52권5호
    • /
    • pp.361-372
    • /
    • 2019
  • 본 연구의 목적은 운동량방정식에서 이송가속도항을 제외한 지배방정식을 이용하여 정형 사각 격자 기반의 2차원 지표면 침수해석 모형을 개발하는 것이다. 공간적 이산화는 유한체적법을 이용하였으며, 시간적 이산화는 음해법을 적용하였다. 모형의 실행시간을 단축하기 위해서 CPU를 이용한 병렬계산 기법을 적용하였다. 개발된 모형의 검증을 위해서 해석해와 비교하고, 가상 도메인에서 수치실험을 통해 모형의 거동을 평가하였다. 또한 국내의 장호원 지역과 모로코의 Sebou 강 지역에 대해서 각기 다른 공간해상도로 침수해석을 수행하고, 그 결과를 CAESER-LISFLOOD (CLF) 모형을 이용한 해석 결과와 비교하였다. 모형의 검증 결과 해석해와 잘 일치된 모의 결과를 나타내었고, 가상 도메인에서의 흐름 해석도 타당한 것으로 평가되었다. 장호원 지역과 Sebou 강 지역에 대한 본 연구와 CLF 모형의 침수모의 결과는 침수심과 침수범위에서 서로 유사하게 나타났으며, 장호원 지역의 경우 홍수위험지도의 침수범위와도 유사한 값을 보였다. 본 연구와 CLF 모형의 모의결과에서 상이한 부분에 대해서는 각각의 모의결과를 비교 평가하였다. 연구결과 본 연구에서 제시된 모형은 홍수터에서의 침수 양상을 잘 모의할 수 있는 것으로 평가되었다. 그러나 본 연구에서 제시된 모형을 이용하여 침수해석을 할 경우에는 도메인 구성 방법과 지배방정식 및 해석 방법에 의한 모형의 특징과 한계점을 충분히 고려해야 할 것이다.

유체-구조물-지반 상호작용을 고려한 해상풍력발전기의 지진응답해석 (Earthquake Response Analysis of an Offshore Wind Turbine Considering Fluid-Structure-Soil Interaction)

  • 이진호;이상봉;김재관
    • 한국지진공학회논문집
    • /
    • 제16권3호
    • /
    • pp.1-12
    • /
    • 2012
  • 이 논문에서는 유체-구조물-지반의 상호작용을 고려한 해상풍력발전기의 지진응답해석법을 제시하였다. 풍력발전기는 tower와 그 정점에 집중된 질량으로 모델링 되었다. 이 tower는 유연한 해저지반에 기초하고 있는 튜브형 cantilever로 이상화하였다. Tower와 해수 간의 동적 상호작용, 기초와 지반간의 동적 상호작용이 고려된 유체-구조물-지반 연성계의 지배방정식은 부분구조법과 Rayleigh-Ritz방법에 의해서 유도되었다. 해수는 압축성 비점성 이상 유체로 이상화하였다. 해수로 포화된 층상지반에 놓인 footing의 동적 강성은 Thin Layer법에 의해서 계산하여 상부구조물 모델과 결합시켰다. 이 해석법을 해상풍력발전기 모델의 지진응답해석에 적용하였다. 해석 결과를 준거해와 비교해서 제안한 해석법의 타당성을 검증하였다. Tower의 유연성, 지반의 강성이 해상풍력발전기 지진거동에 미치는 영향을 분석하였다. 유체-구조물 상호작용과 지반-구조물 상호작용의 지진응답에 대한 상대적인 중요도를 비교 평가하였다.

L자형 이동상수로에서 댐 붕괴파의 수치해석 (Numerical Analysis of Dam-break Waves in an L-shaped Channel with a Movable Bed)

  • 김대근;황건
    • 한국수자원학회논문집
    • /
    • 제45권3호
    • /
    • pp.291-300
    • /
    • 2012
  • 댐 붕괴 직후에 댐 붕괴 근역에서의 댐 붕괴파 및 댐 붕괴파로 인한 하상변동을 추적하기 위하여 RANS를 지배방정식으로 하는 FLOW-3D를 이용한 3차원 수치모의를 수행하였다. 특히 이동상에서의 하상변동과 그 하상변동으로 인한 홍수파의 변동을 고정상에서의 해석결과와 비교하였다. 그 해석 결과를 정리하면 다음과 같다. 첫째, L자형 실험수로에서의 홍수파 해석 및 이동상 실험수로에서의 홍수파와 하상변동 해석 결과는 해당 수리실험을 만족스럽게 재현하고 있다. 둘째, 부유사의 농도는 홍수파의 전면에서 가장 높은 값을 보이며, 하상침식은 흐름이 급변하는 댐 직하류 지점에서 가장 크게 발생한다. 수로의 상류에서 발생하는 부유사로 인해 하류에서는 침식과 퇴적이 번갈아서 발생한다. 홍수파가 도달하는 초기에는 L자형 만곡부의 내측에서 침식이 우세하게 발생하나, 시간이 지나며 침식은 점차 만곡부의 외측으로 이동하는 양상을 보인다. 셋째, L자형 이동상에서의 홍수파는 하상의 침식 및 퇴적으로 인해 고정상에서의 홍수파에 비해 그 전파가 지체되며 홍수위가 크게 상승한다.

나노 스케일 판의 좌굴해석을 위한 비국소 탄성 이론의 적용 (Application of nonlocal elasticity theory for buckling analysis of nano-scale plates)

  • 이원홍;한성천;박원태
    • 한국산학기술학회논문지
    • /
    • 제13권11호
    • /
    • pp.5542-5550
    • /
    • 2012
  • Eringen의 비국소 탄성이론을 이용한 3차 전단변형이론을 정식화 하였고 비국소 탄성이론이 적용된 평형방정식을 유도하였다. 비국소 탄성 이론은 미소 규모 효과를 고려할 수 있고 3차원 전단변형이론은 나노 판의 두께방향으로의 전단변형률과 전단응력의 곡선변화 효과를 고려할 수 있다. 모든 변이 단순지지된 나노-스케일 판의 지배방정식을 풀기 위해 Navier 방법을 사용하였다. 비국소 변수의 효과를 나타내기 위한 나노-스케일 판의 해석적 좌굴하중을 제시하였다. 국소 탄성이론과의 관계를 수치해석 결과를 통하여 고찰하였다. 또한 (i) 나노-스케일 판의 크기, (ii) 비국소 계수, (iii) 형상비 그리고 (iv) 모드 수 등이 나노-스케일 판의 무차원 좌굴하중에 미치는 효과에 대하여 관찰하였다. 본 연구의 결과를 검증하기 위해 참고문헌의 결과들과 비교 분석하였으며 해석결과는 참고문헌의 결과들과 잘 일치함을 알 수 있었다. 비국소 이론에 의한 나노-스케일 판의 좌굴해석에 관한 연구는 향후 관련연구에 비교자료로 활용될 수 있을 것이다.

CFT 기둥의 비선형 유한요소해석을 위한 개선된 강관-콘크리트 간 부착 모델 개발 (An Improved Bond Slip Model of CFT Columns for Nonlinear Finite Element Analysis)

  • 권양수;곽효경;황주영;김진국;김종민
    • 한국전산구조공학회논문집
    • /
    • 제28권2호
    • /
    • pp.213-220
    • /
    • 2015
  • 본 논문에서는 CFT 구조의 강관과 내부 충전 콘크리트 간 복합거동을 유한요소해석 시 적절하게 반영하기 위해 강관과 콘크리트 간 부착 슬립관계 묘사를 위한 알고리즘을 제시하였다. 내부 충전 콘크리트에 축방향 하중 발생 시, 강관과 콘크리트 간 마찰로 인해 강관으로 하중이 전달되며, 이에 따른 강관 슬립량과 힘의 평형관계를 통해 등가강성을 통해 부착관계를 파악할 수 있다. 실제 원형 CFT 부재의 부착응력 실험을 통해 측정된 수직 및 수평 방향 응력 분포 결과와 제안된 해석 기법을 통해 산정된 응력 분포의 비교를 통해 제안된 해석 기법의 타당성을 검증하였다. 또한 비선형 유한요소해석 시 강관과 콘크리트의 부착 거동 묘사에 따라 CFT 기둥의 거동 특성에 영향을 미치게 되므로 축방향 하중이 작용하는 CFT 부재 실험 결과와 제안된 부착-슬립 모델을 반영한 유한요소해석 결과의 하중-변위 곡선 관계 비교를 통해 제안된 기법의 적합성을 검증하였다.

가스 하이드레이트 부존층의 하부 경계면을 규명하기 위한 심도영역 탄성파 구간속도 분석 (Seismic interval velocity analysis on prestack depth domain for detecting the bottom simulating reflector of gas-hydrate)

  • 고승원;정부흥
    • 한국신재생에너지학회:학술대회논문집
    • /
    • 한국신재생에너지학회 2005년도 춘계학술대회
    • /
    • pp.638-642
    • /
    • 2005
  • For gas hydrate exploration, long offset multichannel seismic data acquired using by the 4km streamer length in Ulleung basin of the East Sea. The dataset was processed to define the BSRs (Bottom Simulating Reflectors) and to estimate the amount of gas hydrates. Confirmation of the presence of Bottom Simulating reflectors (BSR) and investigation of its physical properties from seismic section are important for gas hydrate detection. Specially, faster interval velocity overlying slower interval velocity indicates the likely presences of gas hydrate above BSR and free gas underneath BSR. In consequence, estimation of correct interval velocities and analysis of their spatial variations are critical processes for gas hydrate detection using seismic reflection data. Using Dix's equation, Root Mean Square (RMS) velocities can be converted into interval velocities. However, it is not a proper way to investigate interval velocities above and below BSR considering the fact that RMS velocities have poor resolution and correctness and the assumption that interval velocities increase along the depth. Therefore, we incorporated Migration Velocity Analysis (MVA) software produced by Landmark CO. to estimate correct interval velocities in detail. MVA is a process to yield velocities of sediments between layers using Common Mid Point (CMP) gathered seismic data. The CMP gathered data for MVA should be produced after basic processing steps to enhance the signal to noise ratio of the first reflections. Prestack depth migrated section is produced using interval velocities and interval velocities are key parameters governing qualities of prestack depth migration section. Correctness of interval velocities can be examined by the presence of Residual Move Out (RMO) on CMP gathered data. If there is no RMO, peaks of primary reflection events are flat in horizontal direction for all offsets of Common Reflection Point (CRP) gathers and it proves that prestack depth migration is done with correct velocity field. Used method in this study, Tomographic inversion needs two initial input data. One is the dataset obtained from the results of preprocessing by removing multiples and noise and stacked partially. The other is the depth domain velocity model build by smoothing and editing the interval velocity converted from RMS velocity. After the three times iteration of tomography inversion, Optimum interval velocity field can be fixed. The conclusion of this study as follow, the final Interval velocity around the BSR decreased to 1400 m/s from 2500 m/s abruptly. BSR is showed about 200m depth under the seabottom

  • PDF

T/R코일프로브를 이용한 원전 SG세관 검사의 와전류탐상 신호해석 (Signal Analysis of Eddy Current Test Using T/R Coil Probe for Inspection of Steam Generator Tube in NPP)

  • 임건규;이향범
    • 조명전기설비학회논문지
    • /
    • 제22권4호
    • /
    • pp.159-165
    • /
    • 2008
  • 본 논문에서는 전자기 유한요소 해석을 통하여 배열와전류프로브의 T/R 코일의 와전류탐상 특성을 해석하였다. 신호해석을 위해 사용된 결함은 Notch 결함이며, 결함의 깊이는 관두께를 기준으로 40[%]로 하였으며, 결함의 위치는 관의 내부 및 외부에 있는 것으로 하였다. Transmit 코일을 중심으로 Receive 코일의 위치를 원주방향으로 $0[^{\circ}]$, $30[^{\circ}]$, $60[^{\circ}]$, $90[^{\circ}]$에 위치시키면서 신호해석을 수행하였다. 프로브의 전자기적 특성을 해석하기 위하여 맥스웰 방정식을 이용하여 지배방정식을 유도하였고, 이를 3차원 유한요소법을 이용하여 수치 해석을 수행하였다. 두 종류 결함의 수치해석 비교 결과 내부결함의 신호가 외부결함보다 크게 발생하였고, Transmit 코일에 대한 Receive 코일의 각도 및 위치 변화시 결함신호의 차이를 확인할 수 있었다. ASME 표준 시험편을 이용한 배열와전류 프로브의 와전류탐상 실험신호와 비교결과 유사한 신호를 확인할 수 있었다. 본 논문의 결과는 배열와전류 프로브의 와전류 탐상 신호 평가시 도움이 될 것이다.

하이브리드 복합재 플라이휠 로터에 작용하는 내압의 효과 (Effect of the Inner Pressure on a Hybrid Composite Flywheel Retor)

  • 오제훈;한상철;김명훈;하성규
    • Composites Research
    • /
    • 제18권1호
    • /
    • pp.45-54
    • /
    • 2005
  • 필라멘트 와인딩 공정으로 제작된 고속 회전용 복합재 플라이휠 로터는 층간분리 현상에 의해 에너지 저장용량이 저하된다. 그리고 기존의 링 타입 허브는 복합재 로터 내측면에 인장력을 가하게 되고. 이는 로터내의 반경방향 인장응력을 가중시켜 로터의 한계 회전수를 저하시킨다. 복합재 로터의 응력해석을 위해서 2차원 평형방정식과 경계조건이 사용되었고, 이를 근거로 강도비를 최소화시키는 최적의 내압이 존재함을 수치적으로 제시하였다. 이러한 최적의 내압을 발생시키기 위해서 원주방향으로 분할된 스플릿 타입 허브를 제안하고, 링 타입과 스플릿 타입 허브의 두께변화에 따른 내압분포의 영향을 제시하였다. 스플릿 타입 허브의 유효성을 검증하기 위해 허브를 포함한 복합재 로터를 제작한 다음, 최대 회전수 40,000rpm까지 파손 없이 스핀 테스트를 수행하였다. 동시에 로터 표면에 4개의 원주방향 및 반경방향 스트레인게이지를 부착하여 변형률을 무선으로 측정하였다. 측정된 변형률은 해석결과와 매우 잘 일치하였다. 특히 반경방향의 응력을 크게 낮출 수 있었고, 반경방향으로 모두 압축 변형률이 발생함을 확인하였다. 결국 스플릿 타입 허브는 플라이휠 로터의 단점인 반경방향의 낮은 강도를 보안하는 효과를 나타내어, 저장에너지 밀도를 증가시킴으로써 대형 고출력 플라이횔 에너지 저장 시스템의 개발 가능성을 제시하였다.