• 지구물리와물리탐사
• /
• 제9권1호
• /
• pp.121-128
• /
• 2006

헬리콥터를 이용한 항공자력탐사는 정해진 고도를 따라 지표면에 평행하게 이루어지지만, 고해상도 탐사에서는 특히 측정이 이루어지는 고도가 너무 변화하여 평탄면으로 간주할 수 없는 경우가 있다. 이 연구에서는 모서리 효과를 조절할 수 있도록 주변 자력원이 포함되는 등가원 방법을 이용하여 이러한 자료를 변환하는 방법을 개발하였고, 3차원적으로 무작위하게 분포하는 점의 자료를 직접적으로 모델화하였다. 이 문제는 일반적으로 under-determined 이지만 CG 법은 최소 norm 해를 찾을 수 있으며, 자력이상을 자력원과 연관시키는 조화함수를 선택할 자유가 있는데, 상향연속 함수 연산자가 선택되면 등가원 자체가 자력이상이 된다. 기본자기장의 방향으로의 자기 쌍극자분포를 자력원으로 선택하면, 자기 쌍극자의 방향을 수직으로 돌려줌으로써 쉽게 자극화 변환 이상을 유도할 수 있다.

• 대한기계학회논문집B
• /
• 제29권9호
• /
• pp.1049-1056
• /
• 2005

A conservative pressure-based finite-volume numerical method has been developed for computing flow and heat transfer by using an unstructured grid system. The method admits arbitrary convex polyhedra. Care is taken in the discretization and solution procedures to avoid formulations that are cell-shape-specific. A collocated variable arrangement formulation is developed, i.e. all dependent variables such as pressure and velocity are stored at cell centers. Gradients required for the evaluation of diffusion fluxes and for second-order-accurate convective operators are found by a novel second-order accurate spatial discretization. Momentum interpolation is used to prevent pressure checkerboarding and the SIMPLE algorithm is used for pressure-velocity coupling. The resulting set of coupled nonlinear algebraic equations is solved by employing a segregated approach, leading to a decoupled set of linear algebraic equations fer each dependent variable, with a sparse diagonally dominant coefficient matrix. These equations are solved by an iterative preconditioned conjugate gradient solver which retains the sparsity of the coefficient matrix, thus achieving a very efficient use of computer resources.

• Nuclear Engineering and Technology
• /
• 제41권9호
• /
• pp.1171-1180
• /
• 2009

The internal core baffle structure of a French Pressurized Water Reactor (PWR) consists of a collection of baffles and formers that are attached to the barrel. The connections are done thanks to a large number of bolts (about 1500). After inspection, some of the bolts have been found cracked. This has been attributed to the Irradiation Assisted Stress Corrosion Cracking (IASCC). The $Electricit\acute{e}$ De France (EDF) has set up a research program to gain better knowledge of the temperature distribution, which may affect the bolts and the whole structure. The temperature distribution in the structure was calculated thanks to the thermal code SYRTHES that used a finite element approach. The heat transfer between the by-pass flow inside the cavities of the core baffle and the structure was accounted for thanks to a strong thermal coupling between the thermal code SYRTHES and the CFD code named Code_Saturne. The results for the CP0 plant design show that both the high temperature and strong temperature gradients could potentially induce mechanical stresses. The CPY design, where each bolt is individually cooled, had led to a reduction of temperatures inside the structures. A new parallel version of SYRTHES, for calculations on very large meshes and based on MPI, has been developed. A demonstration test on the complete structure that has led to about 1.1 billion linear tetraedra has been calculated on 2048 processors of the EDF Blue Gene computer.

• Journal of Mechanical Science and Technology
• /
• 제20권12호
• /
• pp.2218-2229
• /
• 2006

A conservative finite-volume numerical method for unstructured grids with the cell-centered method has been developed for computing flow and heat transfer by combining the attractive features of the existing pressure-based procedures with the advances made in unstructured grid techniques. This method uses an integral form of governing equations for arbitrary convex polyhedra. Care is taken in the discretization and solution procedure to avoid formulations that are cell-shape-specific. A collocated variable arrangement formulation is developed, i.e. all dependent variables such as pressure and velocity are stored at cell centers. For both convective and diffusive fluxes the forms superior to both accuracy and stability are particularly adopted and formulated through a systematic study on the existing approximation ones. Gradients required for the evaluation of diffusion fluxes and for second-order-accurate convective operators are computed by using a linear reconstruction based on the divergence theorem. Momentum interpolation is used to prevent the pressure checkerboarding and a segregated solution strategy is adopted to minimize the storage requirements with the pressure-velocity coupling by the SIMPLE algorithm. An algebraic solver using iterative preconditioned conjugate gradient method is used for the solution of linearized equations. The flow analysis code (PowerCFD) developed by the present method is evaluated for its application to several 2-D structured-mesh benchmark problems using a variety of unstructured quadrilateral and triangular meshes. The present flow analysis code by using unstructured grids with the cell-centered method clearly demonstrate the same accuracy and robustness as that for a typical structured mesh.

• 지구물리와물리탐사
• /
• 제7권3호
• /
• pp.207-212
• /
• 2004

자기지전류(MT) 자료의 3차원 역산에 대해 소개한다. MT 자료의 역산 문제는 기본적으로 악조건이므로 유일한 해가 존재하지 않는다. 이러한 비유일성을 줄이고 정확한 역산해를 구하기 위해서는 역산 시 사전정보를 추가하는 제약조건을 가해야 한다. 물리탐사 분야에서 비선형 역산에 사용되는 가장 일반적인 방법은 일련의 선형화된 역산문제를 푸는 Gauss-Newton법이다. 이 알고리듬은 수렴 시, 모델 공간에서 역산문제에 대한 목적함수를 최소화하는 최적해를 준다. 그러나 이러한 반복적 선형화기법은 3차원 MT 역산의 경우 Jacobian 행렬을 구하기 힘들기 때문에 그 유용성에 한계가 있다. 이러한 어려움은 CG법에 의해 완화할 수 있다. 선형 CG법은 Gauss-Newton 반복의 각 단계를 근사적으로 풀기 위해서 사용된다. 한편 비선형 CG법은 목적함수의 최소화에 직접적으로 적용된다. 이들 CG법은 Jacobian 행렬의 계산 및 대형 선형방정식의 해를 반복 당 세 번의 모델링으로 대치할 수 있어서 3차원 역산에 적합하다.

• Journal of Animal Science and Technology
• /
• 제46권2호
• /
• pp.137-144
• /
• 2004

불연속 범주형 자료에 대한 잠재변수가 존재한다는 가정하에 임계값을 추정하고 잠재변수를 생성하며 생성된 잠재변수 및 기타 연속변량에 대한 관측치를 포함하는 다변량 임계개체모형을 설정하고 유전능력을 예측하기 위한 방법을 제시하였다. 각각의 범주형 조사 자료의 특성을 갖는 형질에 있어서 임계점의 추정은 추정 가능한 임계점에 대한 1차 미분값(gradient)과 2차 미분값(Hessian)을 이용한 Newton 방법을 이용하면 추정가능하며 지역모수인 육종가의 추정은 PCG 방법으로 구현 가능하다. 이러한 이론은 Quaas(2001)가 제시한 하나의 이산형 자료와 하나의 연속형 자료의 2변량 동시 분석방법을 확장하여 전개한 것이며 이때 잠재변수 및 임계점의 추정은 기타 형질의 잔차 회귀계수 및 상관을 고려해야 한다. 본 연구를 위한 모의실험은 2개의 연속변량으로 체중과 유량을 고려하였고 또 다른 2개의 불연속 변량인 분만난이도와 출생시 생존유무를 고려하여 4형질 동시 분석을 실시하였다. 임계모형에 의한 육종가 추정치의 정확도는 4개의 구간으로 분류되어 기록된 분만난이도의 경우에 91${\sim}$92%의 정확도를 보였고 이항분포인 분만시 생존유무에 대하여는 87~89%의 정확도를 보였다. 반면에 이들 범주형 자료를 선형으로 간주하고 분석한 선형 동물개체 혼합모형에서는 72${\sim}$84% 및 59${\sim}$70%으로 비교적 낮은 추정의 정확도를 보였다. 따라서 범주형 자료의 유전분석은 선형 혼합모형 보다 임계형 혼합모형이 크게 타당할 것으로 사료되었다.

• Investigative Magnetic Resonance Imaging
• /
• 제5권1호
• /
• pp.57-65
• /
• 2001

목적 : 초고속 Echo Planar Imaging (EPI) 의 데이터 측정시간을 반으로 단축시킨 half courier EPI와 임상에서 널리 사용하는 $T_2$ 강조 영상을 결합한 $T_2$-weighted half courier EPI(T2HEPl) 영상기법을 제안하였다. 제안한 방법으로 강한 $T_2$ 대조도를 갖는 $128{\times}128$의 고해상도 EPI 영상을 single scan으로 얻고자 한다. 대상 및 방법 : 초고속 EPI 기법의 데이터 측정시간을 줄이기 위하여 k-space의 절반만을 측정하고, 나머지 절반은 conjugate symmetry 성질을 이용하여 재구성한다. 따라서 64개의 에코로 $128{\times}128$의 고해상도 single shot 영상을 얻을 수 있다. 또한 k-space 데이터의 순서를 시간 축에서 조절하여 강한 $T_2$ 대조도를 갖는 영상을 얻을 수 있다 Eddy current의 영향으로 phase encoding 방향으로 잔류 경사자계가 있을 경우 측정된 데이터들은 k-space에서 이동이 되며 , 이것은 재구성 영상에 심각한 문제점을 초래하게 된다. 본 논문에서는 pre-scan에서 얻은 기준 데이터와 dc에 해당하는 에코간의 상관도를 측정하여 이동을 추정하며, 초기의 phase encoding gradient의 크기를 조정하여 이러한 이동을 제거한다. 결과 : 제안한 $T_2$-weighted half courier EPI 영상기법을 1.0 Tesla 전신 MRI 시스템에 적용하였다. 실험 조건은 single shot으로 TR은 무한대이고, TE는 72ms와 96ms로 설정하였다. 제안한 영상기법을 이용하여 single scan으로 강한 $T_2$ 대조도를 갖는 128x128의 고해상도 EPI 영상을 얻을 수 있었다. 결론 : 제안한 Half Fourier기법을 이용하여 기존의 $64{\times}64$ 영상보다 해상도를 높인 $128{\times}128$ EPI영상을 single scan으로 얻을수 있었으며, 에코의 적절한 배치를 통해 임상에서 널리 사용되는 $T_2$ 대조도를 얻을 수 있었다. 제안한 방법은 특별한 하드웨어의 추가 없이, 펄스 시퀀스와 tuning 및 재구성 알고리즘 등의 소프트웨어적인 방법만으로 구현이 가능하여 많은 임상 응용이 기대된다.