• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 2000년도 유체기계 연구개발 발표회 논문집
• /
• pp.247-252
• /
• 2000

Flow through turbine flow meter is simulated by solving the incompressible Navier-Stockes equations. The solution method is based on the pseudocompressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel Line relaxation method. The equations are solved steadily in rotating reference frames and the centrifugal force and the Coriolis force are added to the equation of motion. The standard k-$\epsilon$ model is employed to evaluate turbulent viscosity.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 추계학술대회논문집B
• /
• pp.573-578
• /
• 2000

Flow through turbine flow meter is simulated by solving the incompressible Navier-Stockes equations. The solution method is based on the pseudocompressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. The equations are solved steadily in rotating reference frames and the centrifugal force and tile Coriolis force are added to the equation of motion. The standard $k-{\varepsilon}$ model is employed to evaluate turbulent viscosity. At first the stability and accuracy of the program is verified with the flow through a square duct with a $90^{\circ}$ bend and on the flat plate.

• 대한수학회지
• /
• 제50권1호
• /
• pp.81-93
• /
• 2013

We consider a distributed optimal flux control problem: finding the potential of which gradient approximates the target vector field under an elliptic constraint. Introducing the Lagrange multiplier and a change of variables the Euler-Lagrange equation turns into a coupled equation of an elliptic equation and a reaction diffusion equation. The change of variables reduces iteration steps dramatically when the Gauss-Seidel iteration is considered as a solution method. For the elliptic equation solver we consider the Cell Boundary Element (CBE) method, which is the finite element type flux preserving methods.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1999년도 추계 학술대회논문집
• /
• pp.59-64
• /
• 1999

The Equations of Chemical kinetics ate very stiff, which forces the use of an implicit scheme. The problem of implicit scheme, however, is that the jacobian must be solved at each time step. In this paper, we examined the approximate chemical jacobian methods such as Gauss-Seidel, Jacobi partial jacobian and diagonalized jacobian that can be stable without full jacobian, We show that Gauss-Seidel jacobian method is stable and accurate as well as full jacobian and that this is more efficient in supersonic combustion problem about $20\%$ than the full jacobian method with same accuracy,

• 한국전산유체공학회지
• /
• 제21권1호
• /
• pp.10-18
• /
• 2016

Numerical boundary conditions are as important as the governing equations when analyzing the fluid flows numerically. An explicit boundary condition method updates the solutions at the boundaries with extrapolation from the interior of the computational domain, while the implicit boundary condition method in conjunction with an implicit time integration method solves the solutions of the entire computational domain including the boundaries simultaneously. The implicit boundary condition method, therefore, is more robust than the explicit boundary condition method. In this paper, steady compressible 2-Dimensional Navier-Stokes solver is developed. We present the implicit boundary condition method coupled with LU-SGS(Lower Upper Symmetric Gauss Seidel) method. Also, the explicit boundary condition method is implemented for comparison. The preconditioning Navier-Stokes equations are solved on unstructured meshes. The numerical computations for a number of flows show that the implicit boundary condition method can give accurate solutions.

• Nuclear Engineering and Technology
• /
• 제55권7호
• /
• pp.2516-2533
• /
• 2023

Several practical methods for accelerating the depletion calculation in a GPU-based Monte Carlo (MC) code PRAGMA are presented including the multilevel spectral collapse method and the vectorized Chebyshev rational approximation method (CRAM). Since the generation of microscopic reaction rates for each nuclide needed for the construction of the depletion matrix of the Bateman equation requires either enormous memory access or tremendous physical memory, both of which are quite burdensome on GPUs, a new method called multilevel spectral collapse is proposed which combines two types of spectra to generate microscopic reaction rates: an ultrafine spectrum for an entire fuel pin and coarser spectra for each depletion region. Errors in reaction rates introduced by this method are mitigated by a hybrid usage of direct online reaction rate tallies for several important fissile nuclides. The linear system to appear in the solution process adopting the CRAM is solved by the Gauss-Seidel method which can be easily vectorized on GPUs. With the accelerated depletion methods, only about 10% of MC calculation time is consumed for depletion, so an accurate full core cycle depletion calculation for a commercial power reactor (BEAVRS) can be done in 16 h with 24 consumer-grade GPUs.

• Tribology and Lubricants
• /
• 제28권4호
• /
• pp.160-166
• /
• 2012

This paper presents a Reynolds equation solver for hydrostatic gas bearings, implemented to run on graphics processing units (GPUs). The original analysis code for the central processing unit (CPU) was modified for the GPU by using the compute unified device architecture (CUDA). The red-black Gauss-Seidel (RBGS) algorithm was employed instead of the original Gauss-Seidel algorithm for the iterative pressure solver, because the latter has data dependency between neighboring nodes. The implemented GPU program was tested on the nVidia GTX580 system and compared to the original CPU program on the AMD Llano system. In the iterative pressure calculation, the implemented GPU program showed 20-100 times faster performance than the original CPU codes. Comparison of the wall-clock times including all of pre/post processing codes showed that the GPU codes still delivered 4-12 times faster performance than the CPU code for our target problem.

• 한국자기학회지
• /
• 제6권1호
• /
• pp.1-6
• /
• 1996

MR/TbCo 박막 계면의 자기교환 결합의 특성을 시뮬레이션 하기 위해서 박막내에 존재 하는 자기 에너지를 최소화 하여 자화의 방향을 구하였다. 수치해석 방법으로는 Newton법과 Gauss-Seidel 방법이 이용되었다. TbCo의 이방성 상수에 비례적으로 자기 이력곡선의 폭이 변하고, TbCo 천이영역의 두께가 얇을수록 자기이력곡선이 벌어지지 않았다. 특히 TbCo 천이영역의 두께가 $400\;\AA$ 일때 단일 곡선(보자력 10 Oe 이하)을 보였다. 또한 계면 교환 결합상수가 작을수록 천이영역의 두께가 얇을 때와 비슷한 곡선을 보였고, 완전히 결합된 상태의 1/100 정도 일때, 단일 곡선을 보였다. 제작한 박막의 자기이력곡선과 비교하여 계면교환 결합상수를 추정할 수 있었다.

• 대한기계학회논문집B
• /
• 제27권8호
• /
• pp.1122-1133
• /
• 2003

Numerical characteristics of implicit upwind schemes, such as upwind ADI, line Gauss-Seidel (LGS) and point Gauss-Seidel (LU) algorithms, for Navier-Stokes equations have been investigated. Time-derivative preconditioning method was applied for efficient convergence at low Mach/Reynolds number regime as well as at large grid aspect ratios. All the algorithms were expressed in approximate factorization form and von Neumann stability analysis was performed to identify stability characteristics of the above algorithms in the presence of high grid aspect ratios. Stability analysis showed that for high aspect ratio computations, the ADI and LGS algorithms showed efficient damping effect up to moderate aspect ratio if we adopt viscous preconditioning based on min-CFL/max-VNN time-step definition. The LU algorithm, on the other hand, showed serious deterioration in stability characteristics as the grid aspect ratio increases. Computations for several practical applications also verified these results.

• 공학논문집
• /
• 제3권1호
• /
• pp.85-95
• /
• 1998

본 논문에서 blocking이 발생되고 다중 customer를 취급하는 유한용량의 Coxian 대기 네트워크를 분석하기 위한 회귀적인 방법을 연구한다. 다중 customer를 나타내기 위한 대용량의 state space로 인해 발생되는 메모리 문제와 rate matrix의 부행렬들에 대한 inverse를 구할 수 없는 관계로 발생될 수 있는 singular 문제를 해결하는 방법을 제안한다. 2가지 blocking 방식(FBFE, HPFE)에 대한 분석을 Gauss-Seidel 방법과 비교하여 본다.