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

It is well known that preconditioning methods are efficient for convergence acceleration at compressible low Mach number flows. In this study, the original Euler equations and three preconditioners nondimensionalized differently are implemented in two dimensional inviscid bump flows using the 3rd order MUSCL and DADI schemes as flux discretization and time integration respectively. The multigrid and local time stepping methods are also used to accelerate the convergence. The test case indicates that a properly modified local preconditioning technique involving concepts of a global preconditioning one produces Mach number independent convergence. Besides, an asymptotic analysis for properties of preconditioning methods is added.

• Nuclear Engineering and Technology
• /
• 제49권6호
• /
• pp.1114-1124
• /
• 2017

Acceleration/preconditioning strategies available in the SCEPTRE radiation transport code are described. A flexible transport synthetic acceleration (TSA) algorithm that uses a low-order discrete-ordinates ($S_N$) or spherical-harmonics ($P_N$) solve to accelerate convergence of a high-order $S_N$ source-iteration (SI) solve is described. Convergence of the low-order solves can be further accelerated by applying off-the-shelf incomplete-factorization or algebraic-multigrid methods. Also available is an algorithm that uses a generalized minimum residual (GMRES) iterative method rather than SI for convergence, using a parallel sweep-based solver to build up a Krylov subspace. TSA has been applied as a preconditioner to accelerate the convergence of the GMRES iterations. The methods are applied to several problems involving electron transport and problems with artificial cross sections with large scattering ratios. These methods were compared and evaluated by considering material discontinuities and scattering anisotropy. Observed accelerations obtained are highly problem dependent, but speedup factors around 10 have been observed in typical applications.

• 한국가시화정보학회지
• /
• 제16권2호
• /
• pp.23-30
• /
• 2018

The present study compares flow fields reconstructed by data assimilation method with different combinations of parameters. As a data assimilation method, Vortex-in-Cell-sharp (VIC#), which supplements additional constraints and multigrid approximation to Vortex-in-Cell-plus (VIC+), is used to reconstruct flow fields from scattered particle tracks. Two parameters, standard deviation of Gaussian radial basis function (RBF) and grid spacing, are mainly tested using artificial data sets which contain few particle tracks. Consequent flow fields are analyzed in terms of flow structure sizes. It is demonstrated that sizes of the flow structures are proportional to an actual scale of the standard deviation of RBF. It implies that a combination of larger grid spacing and smaller standard deviation which preserves the actual standard deviation is able to save computational resources in case of a low track density. In addition, a simple comparison using an experimental data filled with dense particle tracks is conducted.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제17권3호
• /
• pp.197-207
• /
• 2013

The Cahn-Hilliard equation was proposed as a phenomenological model for describing the process of phase separation of a binary alloy. The equation has been applied to many physical applications such as amorphological instability caused by elastic non-equilibrium, image inpainting, two- and three-phase fluid flow, phase separation, flow visualization and the formation of the quantum dots. To solve the Cahn-Hillard equation, many numerical methods have been proposed such as the explicit Euler's, the implicit Euler's, the Crank-Nicolson, the semi-implicit Euler's, the linearly stabilized splitting and the non-linearly stabilized splitting schemes. In this paper, we investigate each scheme in finite-difference schemes by comparing their performances, especially stability and efficiency. Except the explicit Euler's method, we use the fast solver which is called a multigrid method. Our numerical investigation shows that the linearly stabilized stabilized splitting scheme is not unconditionally gradient stable in time unlike the known result. And the Crank-Nicolson scheme is accurate but unstable in time, whereas the non-linearly stabilized splitting scheme has advantage over other schemes on the time step restriction.

• 한국항공우주학회지
• /
• 제33권8호
• /
• pp.8-17
• /
• 2005

저속 압축성 유동에서 사용하는 예조건화 기법은 수렴성 증진에 효과적이다. 본 연구에서는 일반적인 오일러 지배 방정식과 각각 다르게 무차원화한 세 가지 종류의 예조건화 기법을 3차 공간 정확도의 MUSCL, DADI, 다중 격자, 국소 시간 전진 기법을 이용하여 2차원 비점성 bump 유동에 적용하였다. 결과적으로 국소 예조건화 기법에 전역 예조건화 기법의 압력 항 무차원화 방법을 적용하면, 마하수에 무관한 수렴 특성을 얻을 수 있다. 또한, 점근해석을 이용하여 각 예조건화 기법의 특성에 대해 언급하였다.

• 한국정보통신학회논문지
• /
• 제26권2호
• /
• pp.181-186
• /
• 2022

포아즌 볼츠만 방정식 (Poisson-Boltzmann equation, PBE)은 생물물리, 콜로이드 화학 등에서 등장하는 문제들을 모델링하는데 사용되는 방정식이다. 따라서 PBE의 수치해를 효율적으로 예측하는 것은 중요한 이슈이다. 저자들은 기존의 연구에서 PBE를 풀기위한 딥러닝 방법을 제안하였으나, 딥러닝을 훈련하기 위한 샘플을 생성하는 시간이 컸다는 어려움이 있었다. 본 논문에서는 FEM 수치해를 생성하는데 걸리는 시간을 줄이는 두가지 방안을 마련하였다. 첫째로 대수 방정식을 만들 때 bilinar form에 포함되는 penalty 파라메터를 실험적으로 조정하였다. 두 번째로, 대수적멀티그리드 기법을 활용하여 대수 방정식의 컨디션 넘버를 meshsize와 무관하게 만들었다. 따라서 PBE 방정식의 대수 방정식을 풀 때 계산 시간을 효과적으로 줄였다. 이러한 대수적 멀티그리드를 사용한 방법은 다양한 분야에서 딥러닝의 샘플을 생성하는데 효과적으로 활용될 수 있을 것으로 기대된다.