• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.40 no.11
• /
• pp.972-978
• /
• 2012

This paper discusses the parallelization of multifrontal solution method, widely used for finite element structural analyses, for a shared memory architecture. Multifrontal method is easier than other linear solution methods because the solution procedure implies that unknowns can be eliminated simultaneously. Two innovative ideas are introduced to achieve optimal solver performance on a shared memory computer. Those are pairing two frontal matrices and splitting the frontal matrix in order to reduce the temporal memory space required by independent computing tasks. Performance comparisons between original algorithm and proposed one prove that proposed method is more computationally efficient on current multicore machines.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.5
• /
• pp.755-761
• /
• 1987

Two-Dimensional incompressible laminar boundary layer with the reversed flow region is computed using the parially parabolized Navier-Stokes equations in primitive variables. The velocities and the pressure are explicity coupled in the difference equation and the resulting penta-diagonal matrix equations are solved by a streamwise marching technique. The test calculations for the trailing edge region of a finite flat plate and Howarth's linearly retarding flows demonstrate that the method is accurate, efficient and capable of predicting the reversed flow region.

• Journal of ILASS-Korea
• /
• v.5 no.1
• /
• pp.13-22
• /
• 2000

The present study has numerically analyzed the vaporization characteristics of fuel droplets in the high temperature convective flow field. The axisymmetric governing equations for mass, momentum, energy, and species are solved by an iterative and implicite unstructured finite-volume method. The moving boundary due to vaporization is handled by the deformable unstructured grid technique. The pressure-velocity coupling in the density-variable flows is treated by the SIMPLEC algorithm. In terms of the matrix solver, Bi-CGSTAB is employed for the numerically efficient and stable convergence. The n-decane is used as a liquid fuel and the initial droplet temperature is 300K. Computations are performed for the nonevaporating and evaporating droplets with the relative interphase velocity(25m/s). The unsteady vaporization process has been simulated up to the nondimensional time, 25. Numerical results indicate that the mathematical model developed in this study succesfully simulates the main features of the droplet vaporization process in the convective environment.

• Journal of Ocean Engineering and Technology
• /
• v.13 no.2 s.32
• /
• pp.138-146
• /
• 1999

The viscous flow around a ship hull is calculated by the use of RANS(Reynolds-averaged Navier-Stokes) solver. Reynolds stresses are midelled by using the k-${epsilon}$ turbulence model and the law is applied near the body. Body fitted corrdinates are introduced for the treatment of the complex boundary of the ship hull form and the governing equations in the physical domain transformed into ones in the computational domain. The transformed equations are numerically solved by an employment of FVM(Finite Volume Method). SIMPLE(Semi-Implicit Pressure Linked Equation) method is adopted in the calculation of pressure and the solution of the sidcretized equation is obtained by the line-by-line method with the use of TDMA(Tri-Diagonal Matrix Algorithme). To assure the proprietty of this computing method, HSVA tanker and Dyne hull are calculated ar both model and ship scale Reynolds number. Their reaults of pressure distributions on fore and aft body, axial velocity contours and transverse velocity velocity vectors and viscous resistance coefficients are compared with other's experiments and calculations.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.12 no.7
• /
• pp.1076-1085
• /
• 2001

In this paper, the FMM(Fast Multipole Method) combined with GMRES(Generalized Minimal RESidual Method) matrix solver is used to extract the parasitic impedance for complicated 3-D structures in uniform dielectric materials which limit the use of MoM(Method of Moment) due to its large computation time and memory requirement. This algorithm is a fast multipole-accelerated method based on quasistatic analysis and is very efficient for computing impedance between conductors. This paper proved the accuracy and efficiency of the FMM by comparing with MoM in simple examples. Finally the parasitic impedance of FC-PGA(Flip Chip Pin Grid Array) Package pins has been extracted by this algorithm and we have considered the possibility of the EMI/EMC problem caused by the signal interference.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.247-247
• /
• 2021

국내외적으로 하도 내의 흐름을 해석하기 위해 다양한 2차원 흐름해석 모형이 적용되고 있다. 2차원 흐름해석 모형은 기존의 1차원 흐름 해석 모형에서 해석하기 어려운 확산형 홍수파 해석에 강점을 가지고 있어 도심 하천의 외수 범람 예측 등에도 사용되고 있다. 하지만 복잡한 지형 형상을 어떻게 격자로 구성하는가에 따라 해석의 효율성과 정확성이 크게 좌우된다. 초기의 2차원 흐름해석 모형은 주로 정형격자 기반의 단순한 셀을 제작하여 구동되었다. 하지만 매우 빠른 유속과 복잡한 형상을 반영하기 위해서는 전체 격자를 조밀하게 구성할 필요가 있으므로 계산 효율이 떨어지는 문제점이 있다. 그렇기 때문에 대안으로 삼각망과 혼합망 등 비정형 격자를 사용하여 필요한 구역만 격자를 조밀하게 구성하는 방법을 사용하고 있지만 이 방법 또한 추가적인 계산 과정에 따른 계산 시간의 증가가 필연적이다. 따라서 최근에는 정형격자와 비정형격자에 대하여 wet-dry front matrix 최적화, 절점제거법 등 다양한 기법을 통하여 계산 효율을 향상시키고 있는 실정이다. 따라서 본 연구에서는 HLLC Rimann solver와 2차 정확도 기법인 MUSCL-Hancock Method를 적용한 유한체적기반 천수방정식을 기반으로 다양한 격자 구성에 따른 2차원 흐름해석 모형의 효율성 분석을 수행하고, 이를 통해 최적의 흐름해석 방안을 제시하고자 한다.

• KIEAE Journal
• /
• v.16 no.1
• /
• pp.89-94
• /
• 2016

This study is part of three-dimensional(3D) heat transfer analysis program developmental process. The program is being developed without it's own built in 3D-modeller. So 3D-model must be created from another 3D-modeller such as generic CAD programs and imported to the developed program. After that, according to the 3D-geometric data form imported model, 3D-mesh created for numerical calculation. But the 3D-model created from another 3D-modeller is likely to have errors in it's geometric data such as mismatch of position between vertexes or surfaces. these errors make it difficult to create 3D-mesh for calculation. These errors are must be detected and cured in the pre-process before creating 3D-mesh. So, in this study four kinds of filters and functions are developed and tested. Firstly, 'vertex error filter' is developed for detecting and curing for position data errors between vertexes. Secondly, 'normal vector error filter' is developed for errors of surface's normal vector in 3D-model. Thirdly, 'intersection filter' is developed for extracting and creating intersection surface between adjacent objects. fourthly, 'polygon-line filter' is developed for indicating outlines of object in 3D-model. the developed filters and functions were tested on several shapes of 3D-models. and confirmed applicability. these developed filters and functions will be applied to the developed program and tested and modified continuously for less errors and more accuracy.

• Journal of the Korean Institute of Telematics and Electronics D
• /
• v.36D no.11
• /
• pp.45-55
• /
• 1999

The small signal analyzer for the stationary drift-diffusion equation is developed. The slotboom variables of the potential, electron and hole concentrations for the response of applied small signal are defined and the stationary drift-diffusion equation is linearlized on DC operation point by $S^3A$ method. Frontal solver, which is used to solve the global matrix, progresses the accuracy of the solution in high frequency and minimizes the requirement of the memory. The simulations are executed on the structure of 3 dimensional N'P junction diode and 2 dimensional n-MOSFET to verify the proposed algorithm. The average relative errors of the conductance and the capacitance compared with MEDICI are about 26% and 0.67 for N'P junction diode and 7.75% and 2.24% for n-MOSFET. The simulation by the proposed algorithm can analyze the stationary drift-diffusion equation for applied small signal in high frequency region about 100GHz.

• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.38-47
• /
• 2005

Substructuring methods are often used in finite element structural analyses. In this study a multi-level substructuring(MLSS) algorithm is developed and proposed as a possible candidate for finite element fluid solvers. The present algorithm consists of four stages such as a gathering, a condensing, a solving and a scattering stage. At each level, a predetermined number of elements are gathered and condensed to form an element of higher level. At the highest level, each sub-domain consists of only one super-element. Thus, the inversion process of a stiffness matrix associated with internal degrees of freedom of each sub-domain has been replaced by a sequential static condensation of gathered element matrices. The global algebraic system arising from the assembly of each sub-domain matrices is solved using a well-known iterative solver such as the conjugare gradient(CG) or the conjugate gradient squared(CGS) method. A time comparison with CG has been performed on a 2-D Poisson problem. With one domain the computing time by MLSS is comparable with that by CG up to about 260,000 d.o.f. For 263,169 d.o.f using 8 x 8 sub-domains, the time by MLSS is reduced to a value less than $30\%$ of that by CG. The lid-driven cavity problem has been solved for Re = 3200 using the element interpolation degree(Deg.) up to cubic. in this case, preconditioning techniques usually accompanied by iterative solvers are not needed. Finite element formulation for the incompressible flow has been stabilized by a modified residual procedure proposed by Ilinca et al.[9].

• Nuclear Engineering and Technology
• /
• v.42 no.5
• /
• pp.474-499
• /
• 2010

This paper presents the mostortant developments implemented in the APOLLO2 spectral code since its last general presentation at the 1999 M&C conference in Madrid. APOLLO2 has been provided with new capabilities in the domain of cross section self-shielding, including mixture effects and transfer matrix self-shielding, new or improved flux solvers (CPM for RZ geometry, heterogeneous cells for short MOC and the linear-surface scheme for long MOC), improved acceleration techniques ($DP_1$), that are also applied to thermal and external iterations, and a number of sophisticated modules and tools to help user calculations. The method of characteristics, which took over the collision probability method as the main flux solver of the code, allows for whole core two-dimensional heterogeneous calculations. A flux reconstruction technique leads to fast albeit accurate solutions used for industrial applications. The APOLLO2 code has been integrated (APOLLO2-A) within the $ARCADIA^{(R)}$ reactor code system of AREVA as cross section generator for PWR and BWR fuel assemblies. APOLLO2 is also extensively used by Electricite de France within its reactor calculation chain. A number of numerical examples are presented to illustrate APOLLO2 accuracy by comparison to Monte Carlo reference calculations. Results of the validation program are compared to the measured values on power plants and critical experiments.