• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년 추계학술대회논문집
• /
• pp.57-70
• /
• 2008

A high-order accurate Euler flow solver based on a discontinuous Galerkin finite-element method has been developed for the numerical simulations of blade-vortex interaction phenomena on unstructured meshes. A free vortex in freestream was investigated to assess the vortex-preserving property and the accuracy of the present flow solver. Blade-vortex interaction problems in subsonic and transonic freestreams were simulated by adopting a multi-level solution-adaptive dynamic mesh refinement/coarsening technique. The results were compared with those of other numerical and experimental methods. It was shown that the present discontinuous Galerkin flow solver can preserve the vortex structure for significantly longer vortex convection time and can accurately capture the complex unsteady blade-vortex interaction flows, including generation and propagation of acoustic waves.

• 한국전산유체공학회지
• /
• 제3권2호
• /
• pp.17-26
• /
• 1998

The three-dimensional incompressible Navier-Stokes equations have been solved by a node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method with Jacobi matrix solver is used for the time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetragedra, prisms, pyramids, hexahedra, or mixed-element grid. Inviscid bump flow is solved to check the accuracy of high order convective flux discretisation. And viscous flows around a circular cylinder and a sphere are studied to show the efficiency and accuracy of the solver.

• 전자공학회논문지D
• /
• 제34D권8호
• /
• pp.41-51
• /
• 1997

The device simulator (BANDIS) which can analyze efficiently the electrical characteristics of the semiconductor devices under the three dimensional stationary conditions on the IBM PC was developed. Poisson, electon and hole continuity equations are discretized y te galerkin method using a tetrahedron as af finite element. The frontal solver which has exquisite data structures and advanced input/output functions is dused for the matrix solver which needs the highest cost in the three dimensional device simulation. The discretization method of the continuity equations used in BANDIS are compared with that of the scharfetter-gummel method used in the commercial three-dimensional device. To verify an accuracy and the efficiency of the discretization method, the simulation results of the PN junction diode and the BJT from BANDIS are compared with those of the commercial three-dimensiional device simulator such as DAVINCI. The maximum relative error within 2% and the average number of iterations needed for the convergence is decreased by more than 20%. The total simulation time of the BJT with 25542 nodes is decreased to about 60% compared with that of DAVINCI.

• 소성∙가공
• /
• 제9권3호
• /
• pp.233-241
• /
• 2000

It is well known that the quality and efficiency of the design of metal forming processes can be significantly improved with the aid of effective numerical simulations. In the present study, a two-and three-dimensional finite element simulation system, CAMP form, was developed for the analysis of metal forming processes in the PC environment. It is composed of a solver based on the thermo-rigid-viscoplastic approach and graphic user interface (GUI) based pre-and post-processors to be used for the effective description of forming conditions and graphic display of simulation results, respectively. In particular, in the case of CAMPform 2D (two-dimensional), as the solver contains an automatic remeshing module which determines the deformation step when remeshing is required and reconstructs the new mesh system, it is possible to carry out simulations automatically without any user intervention. Also, the forming analysis considers ductile fracture of the workpiece and wear of dies for better usage of the system. In the case of CAMPform 3D, general three-dimensional problems that involve complex die geometries and require remeshing can be analyzed, but full automation of simulations has yet to be achieved. In this paper, the overall structure and computational background of CAMPform will be briefly explained and analysis results of several forming processes will be shown. From the current results, it is construed that CAMPform can be used in providing useful information to assist the design of forming processes.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2004년도 봄 학술발표회 논문집
• /
• pp.389-395
• /
• 2004

We consider the interface between the parallel distributed memory multifrontal solver and the finite element method. We give in detail the requirement and the data structure of parallel FEM interface which includes the element data and the node array. The full procedures of solving a large scale structural problem are assumed to have pre-post processors, of which algorithm is not considered in this paper. The main advantage of implementing the parallel FEM interface is shown up in the case that we use a distributed memory system with a large number of processors to solve a very large scale problem. The memory efficiency and the performance effect are examined by analyzing some examples on the Pegasus cluster system.

• 전산구조공학
• /
• 제6권4호
• /
• pp.83-88
• /
• 1993

시간변수에 대하여 불연속성을 주는 시간불연속 Galerkin 방법을 유한요소법으로 해석하였다. 이 방법은 미분방정식 관점에서 지금까지 요소간에 연속성을 준 일반적 유한요소법과 다르게 임의의 시간요소를 선택, 매 시간단계에서 요소경계에 불연속을 허락함으로서 해의 정확성을 높이고 무조건의 안정을 주는 상미분방정식의 해법인 것이다.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제3권4호
• /
• pp.274-285
• /
• 2011

The supercavitating vehicle is an underwater vehicle that is surrounded almost completely by a supercavity to reduce hydrodynamic drag substantially. Since the cruise speed of the vehicle is much higher than that of conventional submarines, the drag force is huge and a buckling may occur. The buckling phenomenon is analyzed in this study through static and dynamic approaches. Critical buckling load and pressure as well as buckling mode shapes are calculated using static buckling analysis and a stability map is obtained from dynamic buckling analysis. When the finite element method (FEM) is used for the buckling analysis, the solver requires a linear static solver and an eigenvalue solver. In this study, these two solvers are integrated and a consolidated buckling analysis module is constructed. Furthermore, Particle Swarm Optimization (PSO) algorithm is combined in the buckling analysis module to perform a design optimization computation of a simplified supercavitating vehicle. The simplified configuration includes cylindrical shell structure with three stiffeners. The target for the design optimization process is to minimize total weight while maintaining the given structure buckling-free.

• 한국산학기술학회논문지
• /
• 제9권2호
• /
• pp.286-290
• /
• 2008

본 논문은 고체 추진 로켓의 연소 중에 발생하는 고체추진체의 동적 파괴 현상 및 유체-구조 상호작용을 시뮬레이션 하기 위한 프로그램 개발에 대한 것이다. 개발된 프로그램은 구조해석을 위한 CVFE (cohesive Volumetric Finite Element) 방법과 외재적 ALE (Arbitrary Lagrangian Eulerian) 방법을 응용한 유한요소법 코드와 유동해석을 위한 외재적 비정렬 유한 체적 오일러 코드(Explicit Unstructured Finite Volume Euler code)로 구성된다. 개발된 프로그램의 또 다른 중요한 특징은 균열의 전파와 고체추진체의 변형에 따라 생기는 추진제 형상의 대변형이 발생할 때, 새로 생긴 유체 영역에서의 격자의 확장과 복구되는 능력이다.

• 대한기계학회논문집 C: 기술과 교육
• /
• 제3권2호
• /
• pp.107-115
• /
• 2015

유연 다물체 동역학은 실제 시스템을 가능한 유사하게 수치화하여 해석할 수 있기 때문에 일반 동역학 연구에 대한 차세대 주제로 각광을 받고 있다. 이러한 유연 다물체 동역학에 대한 해석 기술은 리커다인이라는 상용 소프트웨어에 효과적으로 적용되어 있는데, 특히 강체와 유연체를 통합하여 하나의 솔버에서 해석을 할 수 있는 특징을 가지고 있다. 본 논문에서는 이러한 리커다인의 유연 다물체 동역학 솔버의 기술들을 살펴보고자 한다. 기본적으로 리커다인의 유연 다물체 동역학 해법은 동시회전 기법을 사용하는 증분 유한요소 정식화를 기존의 순환공식을 이용한 동역학 정식화에 결합함으로써 구현되어 진다. 이 과정에서 강체와 유연체 사이의 조인트나 힘 요소 등의 효율적 처리를 위해 가상 바디와 유연체 조인트 개념이 사용된다.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제6권4호
• /
• pp.1024-1040
• /
• 2014

The theoretical background and a numerical solution procedure for a time domain hydroelastic code are presented in this paper. The code combines a VOF-based free surface flow solver with a flexible body motion solver where the body linear elastic deformation is described by a modal superposition of dry mode shapes expressed in a local floating frame of reference. These mode shapes can be obtained from any finite element code. The floating frame undergoes a pseudo rigid-body motion which allows for a large rigid body translation and rotation and fully preserves the coupling with the local structural deformation. The formulation relies on the ability of the flow solver to provide the total fluid action on the body including e.g. the viscous forces, hydrostatic and hydrodynamic forces, slamming forces and the fluid damping. A numerical simulation of a flexible barge is provided and compared to experiments to show that the VOF-based flow solver has this ability and the code has the potential to predict the global hydroelastic responses accurately.