• 제목/요약/키워드: N-S equation solver

검색결과 12건 처리시간 0.017초

우수한 수렴특성을 갖는 3차원 포아송 방정식의 이산화 방법 (A discretization method of the three-dimensional poisson's equation with excellent convergence characteristics)

  • 김태한;이은구;김철성
    • 전자공학회논문지D
    • /
    • 제34D권8호
    • /
    • pp.15-25
    • /
    • 1997
  • The integration method of carier concentrations to redcue the discretization error of th box integratio method used in the discretization of the three-dimensional poisson's equation is presented. The carrier concentration is approximated in the closed form as an exponential function of the linearly varying potential in the element. The presented method is implemented in the three-dimensional poisson's equation solver running under the windows 95. The accuracy and the convergence chaacteristics of the three-dimensional poisson's equation solver are compared with those of DAVINCI for the PN junction diode and the n-MOSFET under the thermal equilibrium and the DC reverse bias. The potential distributions of the simulatied devices from the three-dimensional poisson's equation solver, compared with those of DAVINCI, has a relative error within 2.8%. The average number of iterations needed to obtain the solution of the PN junction diode and the n-MOSFET using the presented method are 11.47 and 11.16 while the those of DAVINCI are 21.73 and 23.0 respectively.

  • PDF

다화학종 GH 방정식의 정확성 향상을 위한 벽면 경계조건 연구 및 GH 방정식의 엔트로피 특성 고찰 (The Slip-Wall Boundary Conditions Effects and the Entropy Characteristics of the Multi-Species GH Solver)

  • 안재완;김종암
    • 한국항공우주학회지
    • /
    • 제37권10호
    • /
    • pp.947-954
    • /
    • 2009
  • 본 연구에서는 다화학종 희박유동의 해석을 위해, GH(Generalized Hydrodynamic) 방정식을 기반으로 한 축대칭 유동 해석이 가능한 다화학종 GH 수치해석 기법을 전산유체역학 수치해석자로서 개발하였다. 최초로 구현된 다화학종 GH 수치기법은 축대칭 형상의 물체 주위의 극초음속 희박유동을 대상으로 하여, DSMC(Direct Simulation Monte Carlo) 및 N-S(Navier-Stokes) 방정식의 결과와의 비교를 통해 정확도를 검증하고자 하였다. GH 해석자의 정확도 향상을 위해 고체 벽면에서의 여러 가지 slip-wall 경계조건을 적용하고 각각의 결과를 비교하였다. 또한, 높은 Knudsen 수의 1차원 수직 충격파 구조 문제를 통해 GH 방정식의 엔트로피 특성 및 정확성을 고찰하였다.

Flexible 효과를 고려한 다물체 시스템의 동역학적 해석에 관한 연구

  • 최승렬;한창수
    • 한국정밀공학회:학술대회논문집
    • /
    • 한국정밀공학회 1992년도 추계학술대회 논문집
    • /
    • pp.349-353
    • /
    • 1992
  • The purpose of this paper is to develop methods for the dynamic analysis of multibody system that consist of interconnected rigid and deformable component. The equations of motion are derived by using the Lagrange's equation and finite element theory for the elastic mechanism systems. The type of equation of motion is the differential algebraic equation included kinematic nonlinear algebraic equation. The generalized coordinate partitioning method is used for solving this equation. To show the validity of this analysis solver, couple of models were canalized and those results were compared with the commercial package(ADAMS).

정상상태에서 드리프트-확산 방정식의 소신호 해석 프로그램 개발 (A Development of the Small Signal Analyzer for the Stationary Drift-Diffusion Equation)

  • 임웅진;이은구;김태한;김철성
    • 전자공학회논문지D
    • /
    • 제36D권11호
    • /
    • pp.45-55
    • /
    • 1999
  • 정상상태에서 반도체 소자의 전기적 특성을 분석할 수 있는 소신호 해석용 시뮬레이터를 개발하였다. 전위와 전자 및 정공 농도의 Slotboom변수에 대한 소신호 응답을 정의한 후 $S^3A$방법을 적용하여 ,DC 동작점에서 반도체 방정식을 선형화하였다. 행렬풀이를 위해 전진해법을 사용하여 메모리 소비량을 최소화하며 고주파 영역에서 소신호 해의 정확성을 향상시켰다. 구현된 알고리즘의 검증을 위해 3차원 구조를 갖는 N'P 다이오드 및 2차원 구조를 갖는 n-MOSFET에 대해 소신호 해석을 수행하여 MEDICI와 비교한 결과, 인가 전압에 따른 컨덕턴스와 캐패시턴스의 평균 상대 오차는 N'P 다이오드에서는 0.87%와 2.6%이고, n-MOSFET의 경우 7.75%와 2.24%로 나타났다. n-MOSFET에 대하여 입력신호의 주파수 변화에 따른 컨덕턴스와 캐패시턴스를 비교한 결과 MEDICI의 경우 10GHz까지 예측한 반면 본 논문이 제시한 방법을 이용한 모의실험은 100GHz까지 예측이 가능하여 고주파 영역에서 모의실험의 정확도가 향상됨을 확인하였다.

  • PDF

이중시간적분법을 이용한 순차적 유동해석 기법 (DELTA-FORMULATION OF A SEGREGATED NAVIER-STOKES SOLVER WITH A DUAL-TIME INTEGRATION)

  • 김종태;탁남일;김상백;김민환;이원재
    • 한국전산유체공학회:학술대회논문집
    • /
    • 한국전산유체공학회 2006년도 추계 학술대회논문집
    • /
    • pp.31-35
    • /
    • 2006
  • The delta-formulation of the Navier-Stokes equations has been popularly used in the aerodynamics area. Implicit algorithm can be easily implemented in that by using Taylor series expansion. This formulation is extended for an unsteady analysis by using a dual-time integration. In the meanwhile, the incompressible flows with heat transfers which occur in the area of thermo-hydraulics have been solved by a segregated algorithm such as the SIMPLE method, where each equation is discretised by using an under-relaxed deferred correction method and solved sequentially. In this study, the dual-time delta formulation is implemented in the segregated Navier-Stokes solver which is based on the collocated cell-centerd scheme with un unstructured mesh FVM. The pressure correction equation is derived by the SIMPLE method. From this study, it was found that the Euler dual-time method in the delta formulation can be combined with the SIMPLE method.

  • PDF

An IE-FFT Algorithm to Analyze PEC Objects for MFIE Formulation

  • Seo, Seung Mo
    • Journal of electromagnetic engineering and science
    • /
    • 제19권1호
    • /
    • pp.6-12
    • /
    • 2019
  • An IE-FFT algorithm is implemented and applied to the electromagnetic (EM) solution of perfect electric conducting (PEC) scattering problems. The solution of the method of moments (MoM), based on the magnetic field integral equation (MFIE), is obtained for PEC objects with closed surfaces. The IE-FFT algorithm uses a uniform Cartesian grid to apply a global fast Fourier transform (FFT), which leads to significantly reduce memory requirement and speed up CPU with an iterative solver. The IE-FFT algorithm utilizes two discretizations, one for the unknown induced surface current on the planar triangular patches of 3D arbitrary geometries and the other on a uniform Cartesian grid for interpolating the free-space Green's function. The uniform interpolation of the Green's functions allows for a global FFT for far-field interaction terms, and the near-field interaction terms should be adequately corrected. A 3D block-Toeplitz structure for the Lagrangian interpolation of the Green's function is proposed. The MFIE formulation with the IE-FFT algorithm, without the help of a preconditioner, is converged in certain iterations with a generalized minimal residual (GMRES) method. The complexity of the IE-FFT is found to be approximately $O(N^{1.5})$and $O(N^{1.5}logN)$ for memory requirements and CPU time, respectively.

THE BUCHSTAB'S FUNCTION AND THE OPERATIONAL TAU METHOD

  • Aliabadi, M.Hosseini
    • Journal of applied mathematics & informatics
    • /
    • 제7권3호
    • /
    • pp.905-915
    • /
    • 2000
  • In this article we discuss some aspects of operational Tau Method on delay differential equations and then we apply this method on the differential delay equation defined by $\omega(u)\;=\frac{1}{u}\;for\;1\lequ\leq2$ and $(u\omega(u))'\;=\omega(u-1)\;foru\geq2$, which was introduced by Buchstab. As Khajah et al.[1] applied the Recursive Tau Method on this problem, they had to apply that Method under the Mathematica software to get reasonable accuracy. We present very good results obtained just by applying the Operational Tau Method using a Fortran code. The results show that we can obtain as much accuracy as is allowed by the Fortran compiler and the machine-limitations. The easy applications and reported results concerning the Operational Tau are again confirming the numerical capabilities of this Method to handle problems in different applications.

Beckward Facing Step의 층류 유동 수치계산 (Numerical Computation of Laminar Flow over a Backward Facing Step)

  • ;반석호;김형태
    • 한국해양공학회지
    • /
    • 제7권2호
    • /
    • pp.150-161
    • /
    • 1993
  • 원초변수를 이용한 Navier-Stokes 방정식의 수치계산기법을 개발하고, 이를 응용하여 backward facing step의 층류 유동을 계산하였다. 직교좌표계에서의 비압축성 Navier-Stokes방정식을 풀기위해 시간과 공간항을 2차 정도의 유한 차분을 사용하여 이산화하였고 비교차격자계를 사용하여 양해법으로 수치 계산하였다. 운동량방정식과 연속방정식으로 부터 유도된 압력방정식(pressure-poisson equation)을 이용하여 무발산 조건을 만족시켰ㄲ다. Backward facing step의 층류 유동을 100.$\leq$R$_e$$\leq$1000 범위에 대해서 수치 계산하였으며 실험결과와 잘 일치하는 결과를 구할 수 있었다. 특히 step뒤에서 생기는 박리구간의 길이는 다른 계산결과들보다 실험치에 가까운 값을 얻을 수 있었으며, Re가 600보다 클때는 위쪽 벽에 또 다른 박리 유동이 발생되는 현상이 예측되었다.

  • PDF

극초음속 희박유동 해석을 위한 축대칭 다화학종 GH 방정식의 개발 (EVELOPMENT OF AXISYMMETRIC MULTI-SPECIES GH EQUATION FOR HYPERSONIC RAREFIED FLOW ANALYSES)

  • 안재완;김종암
    • 한국전산유체공학회:학술대회논문집
    • /
    • 한국전산유체공학회 2008년도 학술대회
    • /
    • pp.84-91
    • /
    • 2008
  • Generalized hydrodynamic (GH) theory for multi-species gas and the computational models are developed for the numerical simulation of hypersonic rarefied gas flow on the basis of Eu's GH theory. The rotational non-equilibrium effect of diatomic molecules is taken into account by introducing excess normal stress associated with the bulk viscosity. The numerical model for the diatomic GH theory is developed and tested. Moreover, with the experience of developing the dia-tomic GH computational model, the GH theory is extended to a multi-species gas including 5 species; O$_2$, N$_2$, NO, O, N. The multi-species GH model includes diffusion relation due to the molecular collision and thermal phenomena. Two kinds of GH models are developed for an axisymmetric flow solver. By compar-ing the computed results of diatomic and multi-species GH theories with those of the Navier-Stokes equations and the DSMC results, the accuracy and physical consistency of the GH computational models are examined.

  • PDF

극초음속 희박유동 해석을 위한 축대칭 다화학종 GH 방정식의 개발 (EVELOPMENT OF AXISYMMETRIC MULTI-SPECIES GH EQUATION FOR HYPERSONIC RAREFIED FLOW ANALYSES)

  • 안재완;김종암
    • 한국전산유체공학회:학술대회논문집
    • /
    • 한국전산유체공학회 2008년 추계학술대회논문집
    • /
    • pp.84-91
    • /
    • 2008
  • Generalized hydrodynamic (GH) theory for multi-species gas and the computational models are developed for the numerical simulation of hypersonic rarefied gas flow on the basis of Eu's GH theory. The rotational non-equilibrium effect of diatomic molecules is taken into account by introducing excess normal stress associated with the bulk viscosity. The numerical model for the diatomic GH theory is developed and tested. Moreover, with the experience of developing the dia-tomic GH computational model, the GH theory is extended to a multi-species gas including 5 species; $O_2,\;N_2$, NO, O, N. The multi-species GH model includes diffusion relation due to the molecular collision and thermal phenomena. Two kinds of GH models are developed for an axisymmetric flow solver. By compar-ing the computed results of diatomic and multi-species GH theories with those of the Navier-Stokes equations and the DSMC results, the accuracy and physical consistency of the GH computational models are examined.

  • PDF