• 제목/요약/키워드: matrix solver

검색결과 86건 처리시간 0.025초

유한요소법에 의한 인체내 전계분포 해석 용 소프트웨어의 개발 (Finite Element Software Package for Analysis of Electric Field Distribution in Human Body)

  • 우응제
    • 대한의용생체공학회:학술대회논문집
    • /
    • 대한의용생체공학회 1993년도 춘계학술대회
    • /
    • pp.66-69
    • /
    • 1993
  • We have developed a software package for the analysis of electric field distribution in human body. It includes the graphical finite element mesh generator, linear system of equations solver using sparse matrix and vector technique, and post-processor for the display of the results. This software package can be used in various research areas of biomedical engineering where we inject current or apply voltage to human body. The software package was developed on Macintosh II computer and the size of the model is only limited by the main memory.

  • PDF

모델 축소법을 이용한 타이어 인터페이싱 개발 (A Development of the Tire Interfacing Using the Reduction Method)

  • 임문수;김영배;조규종;정광용
    • 한국CDE학회논문집
    • /
    • 제8권2호
    • /
    • pp.109-114
    • /
    • 2003
  • In order to develop the reduced tire modal model for analyzing a full tire model, the Craig-Bampton method is utilized in this paper. When the tire contacts the road, the Abaqus solver extracts the condensed stiffness, coupled mass and mode shape matrix about the node, which contacts the road. The Abaqus full tire model is reduced using the substructure method utilizing Craig-Bampton algorithm. Then, the extracted matrices are interfaced with the superelement, which is fed to the Nastran reduction algorithm. Eventually, the reduced tire model is verified from experiment and various reduction parameters (i.e. modal number, reduction point, etc.) are studied for the effectiveness of the proposed paper.

레이어가 있는 하프스페이스에서 페라이트코아가 있는 와류탐침에 대한 연구 (A Study on Eddy-current Probe with Ferrite Cores over a Layered Half-Space)

  • 김태원;변기량;최재훈;강유석;황호정
    • 대한전자공학회:학술대회논문집
    • /
    • 대한전자공학회 1998년도 추계종합학술대회 논문집
    • /
    • pp.613-616
    • /
    • 1998
  • In this paper, a model of a Eddy-current probe coil with a ferrite core in the presence of a half-space with a layer is developed. The half-space with a layer is accounted for by computing the appropriate Green's function by using Bessel transforms. Upon introducing equivalent Amperian currents within a core to explain effect to a impedance change in the coil due to a (ferrite) core, we derive a volume integral equation, The integral equation is transformed via the method of moments into a vector-matrix equation, which is then solved using a linear equation solver. Through the above processing, we computed impedance value in a Eddy-current probe coil due to a conductivity change of layer.

  • PDF

Treatment of non-resonant spatial self-shielding effect of double heterogeneous region

  • Tae Young Han;Hyun Chul Lee
    • Nuclear Engineering and Technology
    • /
    • 제55권2호
    • /
    • pp.749-755
    • /
    • 2023
  • A new approximation method was proposed for treating the non-resonant spatial self-shielding effects of double heterogeneous region such as the double heterogeneous effect of VHTR fuel compact in the thermal energy range and that of BP compact with BISO. The method was developed based on the effective homogenization method and a spherical unit cell model with explicit coated layers and a matrix layer. The self-shielding factor was derived from the relation between the collision probabilities for a double heterogeneous compact and the effective cross section for the homogenized compact. First, the collision probabilities and transmission probabilities for all layers of the spherical model were calculated using conventional collision probability solver. Then, the effective cross section for the homogenized sphere cell representing the homogenized compact was obtained from the transmission probability calculated using the probability density function of a chord length. The verification calculations revealed that the proposed method can predict the self-shielding factor with a maximum error of 2.3% and the double heterogeneous effect with a maximum error of 200 pcm in the typical VHTR problems with various packing fractions and BP compact sizes.

대규모 자유도 문제의 구조해석을 위한 병렬 알고리즘 (A Parallel Algorithm for Large DOF Structural Analysis Problems)

  • 김민석;이지호
    • 한국전산구조공학회논문집
    • /
    • 제23권5호
    • /
    • pp.475-482
    • /
    • 2010
  • 본 논문에서는 대규모 자유도 시스템의 병렬처리를 위하여 2단계로 이루어진 영역분할법(Domain Decomposition Method) 기반의 병렬 알고리즘을 제안하였다. 분할된 영역의 내부 및 외부 경계를 상위영역문제로 정의하고 국부영역문제는 변위 경계조건이 모두 주어지는 분할영역에서의 Dirichlet 문제로 구성한다. 상위영역에서는 전체 상위영역에 대한 강성 행렬의 어셈블이 필요없는 반복법을 통하여 변위를 구하고, 이를 바탕으로 국부영역에서 Multi-Frontal Sparse Solver (MFSS)를 이용하여 변위를 계산한다. 상위영역문제의 연산에서 프로세서 간의 데이터 교환을 최소화하여 계산효율을 유지하며, 동시에 해석 가능한 자유도를 증대시키는 병렬 PCG(Preconditioned Conjugate Gradient)법 기반의 알고리즘을 개발하였다. 제안된 알고리즘을 적용하여 수치해석을 수행한 결과, 프로세서 수가 증가할수록 계산성능의 손실없이 해석 가능한 자유도가 비례하여 증가하는 선형 확장성을 관찰할 수 있었으며, 대규모 자유도 문제에 효과적으로 사용 가능함을 확인하였다.

전달 루틴의 병렬화를 통한 SAT 알고리즘의 GPGPU 가속화 (GPGPU Acceleration of SAT Algorithm with Propagation Routine Parallelization)

  • 강형주
    • 한국정보통신학회논문지
    • /
    • 제20권10호
    • /
    • pp.1919-1926
    • /
    • 2016
  • 대량의 데이터를 병렬적으로 처리할 수 있는 General-Purpose Graphics Processing Unit(GPGPU)가 최근 많은 분야에서 적용되고 있으며, 이는 전자 설계 자동화 분야에서도 예외가 아니다. SAT 알고리즘은 다양한 전자 설계 자동화 문제에 적용되는 대표적인 알고리즘 중 하나이다. GPGPU를 이용해서 SAT 알고리즘을 가속화하기 위해 노력이 이루어져 왔으나, SAT 알고리즘 자체의 특성으로 인해 병렬화에 어려움이 있어왔다. 이 논문에서는 SAT 알고리즘의 내부 과정 중 비교적 병렬화가 용이한 전달 루틴을 병렬화함으로써 GPGPU 가속화를 적용하였다. 전달 루틴이 희소 행렬의 곱셈과 유사한 점에 착안하여 데이터 구조를 구성하고 이에 맞추어서 병렬적인 전달 루틴을 작성하였다. 병렬적으로 동작하는 쓰레드들 사이의 데이터 손실을 방지하기 위해 아토믹(atomic) 연산을 이용하였다. 벤치마크 SAT 문제들에 대해 기존의 GPGPU 기반 SAT solver에 비해 성능이 10배 이상 향상되었음을 확인하였다.

Topological Design Sensitivity on the Air Bearing Surface of Head Slider

  • Yoon, Sang-Joon;Kim, Min-Soo;Park, Dong-Hoon
    • Journal of Mechanical Science and Technology
    • /
    • 제16권8호
    • /
    • pp.1102-1108
    • /
    • 2002
  • In this study, a topological design sensitivity of the ai. bearing surface (ABS) is suggested by using an adjoint variable method. The discrete form of the generalized lubrication equation based on a control volume formulation is used as a compatible condition. A residual function of the slider is considered as an equality constraint function, which represents the slider in equilibrium. The slider thickness parameters at all grid cells are chosen as design variables since they are the topological parameters determining the ABS shape. Then, a complicated adjoint variable equation is formulated to directly handle the highly nonlinear and asymmetric coefficient matrix and vector in the discrete system equation of air-lubricated slider bearings. An alternating direction implicit (ADI) scheme is utilized for the numerical calculation. This is an efficient iterative solver to solve large-scale problem in special band storage. Then, a computer program is developed and applied to a slider model of a sophisticated shape. The simulation results of design sensitivity analysis (DSA) are directly compared with those of FDM at the randomly selected grid cells to show the effectiveness of the proposed approach. The overall distribution of DSA results are reported, clearly showing the region on the ABS where special attention should be given during the manufacturing process.

수치해를 이용한 선박의 점성저항 해석 (Visous resistance analysis of a ship using numerical solutions)

  • 곽영기
    • 한국해양공학회지
    • /
    • 제11권2호
    • /
    • pp.100-106
    • /
    • 1997
  • Viscous flow around an actual ship is calculated by an use of RANS(Reynolds-averaged Navier-Stokes) solver. Reynolds stress is modelled by using k-$\varepsilon$ turbulence model and the law of wall is applied near the body. Body fitted coordinates are introduced for the treatment of the complex boundary of the ship hull form. The transformed equations in the computational domain are numerically solved by an employment of FVM(Finite Volume Method). SIMPLE(Semi-Implcit Pressure Linked Equation) method is adopted in the calculation of pressure and the solution of the disssssssscretized equation is obtained by the line-by-line method with the use of TDMA(Tri-Diagonal Matrix Algorithme). The subject ship model of actual calculation is 4,410 TEU class container carrier. For 4 geosim models the calculated viscous resistancce values are compared with the model test results and analyzed on their componentss. The resistance performance of an actual ship is predicted very resonably, so this mothod may be utilized as a design tool of hull form.

  • PDF

연속적인 최대-최소 연결비율 문제: 회선망에서의 공정성 및 효율성을 보장하는 경로설정 (Successive Max-min Connection-Ratio Preoblem:Routing with Fairness and Efficiency in Circuit Telecommunication Networks)

  • 박구현;우재현
    • 한국경영과학회지
    • /
    • 제22권2호
    • /
    • pp.13-29
    • /
    • 1997
  • This paper considers a new routing problem, successive max-min connection ratio problem (SMCRP), arised in circuit telecommunication networks such as SONET and WDM optical transport network. An optimization model for SMCRP is established based on link-flow formulation. It's first optimization process is an integral version of maximum concurrent flow problem. Integer condition does not give the same connection-ratio of each node-pair at an optimal solution any more. It is also an integral multi-commodity flow problem with fairness restriction. In order to guarantee fairness to every node-pair the minimum of connection ratios to demand is maximized. NP- hardness of SMCRP is proved and a heuristic algorithm with polynomial-time bound is developed for the problem. Augmenting path and rerouting flow are used for the algorithm. The heuristic algorithm is implemented and tested for networks of different sizes. The results are compared with those given by GAMS/OSL, a popular commercial solver for integer programming problem.n among ferrite-pearlite matrix, the increase in spheroidal ratio with increasing fatigue limitation, 90% had the highest, 14.3% increasing more then 70%, distribution range of fatigue.ife was small in same stress level. (2) $\sqrt{area}_{max}$ of graphite can be used to predict fatigue limit of Ductile Cast Iron. The Statistical distribution of extreme values of $\sqrt{area}_{max}$ may be used as a guideline for the control of inclusion size in the steelmaking.

  • PDF

다중프론트 해법의 공유메모리 병렬화 (Parallelization of Multifrontal Solution Method for Shared Memory Architecture)

  • 김민기;김정호;박찬익;김승조
    • 한국항공우주학회지
    • /
    • 제40권11호
    • /
    • pp.972-978
    • /
    • 2012
  • 본 논문은 유한요소 구조해석의 선형해법으로 널리 사용되는 다중프론트 해법의 공유메모리 환경하의 병렬화 방법을 논의한다. 다중프론트 해법은 병렬성이 내재되어 있어서 여타 해법보다 상대적으로 병렬화가 용이한 방법이다. 다중프론트 해법의 공유메모리 컴퓨터에서 최적의 성능을 내도록 병렬 계산을 수행하기 위한 기법들이 제시되었다. 주로 독립적인 계산 작업 시에 필요한 주 메모리 용량을 줄이는 데 초점을 맞춘 방법들로서 프론트 행렬 연성화와 행렬 분리로 명명된 두 기법에 대해 자세히 설명한다. 개발된 방법으로 기존의 알고리즘과의 성능 비교를 수행하여 본지에 제안한 방법이 현대의 다중코어 컴퓨터에서 훨씬 더 효율적인 기법임을 입증하였다.