• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.49-54
• /
• 1998

3-dimensional Euler solver is parallelized. The spatial discretization method is the 2nd order TVD scheme and DADI method with multigrid is used as a time integration. In order to parallelize this solver, the domain decomposition method with overlapped grid and message passing techniques are used. The informations on the each inter-processor bound-aries are communicated with MPI library. Finally, the parallel performance repsented by calculating the ONERA M6 wing at transonic flow condition using CRAY T3E and C90.

• 한국통신학회논문지
• /
• 제21권10호
• /
• pp.2724-2730
• /
• 1996

In this paper, the dispersive characteristics of the Finite-Difference Frqequency-Domain method based on the Multi-Resolution Technique(MR-FDFD) are numerically analyzed. A dispersion analysis of the MR-FDFD ority of the MR-FDFD method to the spatial discretization is shown. We expect that the multi-resoluation technique will improve the disavantage of the finite difference techqnique which needs the large comutational memory for accurate electromagnetic analysis.

• Geomechanics and Engineering
• /
• 제24권2호
• /
• pp.167-178
• /
• 2021

Spatial variability is an inherent uncertainty of soil properties. Current reliability analyses generally incorporate random field theory and Monte Carlo simulation (MCS) when dealing with spatial variability, in which the computational efficiency is a significant challenge. This paper proposes a KL-FORM algorithm to improve the computational efficiency. In the proposed KL-FORM, Karhunen-Loeve (KL) expansion is used for discretizing random fields, and first-order reliability method (FORM) is employed for reliability analysis. The KL expansion and FORM can be used in conjunction, through adopting independent standard normal variables in the discretization of KL expansion as the basic variables in the FORM. To illustrate the effectiveness of this KL-FORM, it is applied to a case study of a strip footing in spatially variable unsaturated soil under rainfall, in which the bearing capacity of the footing is computed by numerical simulation. This case study shows that the KL-FORM is accurate and efficient. The parametric analyses suggest that ignoring the spatial variability of the soil may lead to an underestimation of the reliability index of the footing.

• 대한기계학회논문집B
• /
• 제24권9호
• /
• pp.1227-1238
• /
• 2000

This study examines the effects of the discretization schemes on numerical solutions of viscoelastic fluid flows. For this purpose, a temporally evolving mixing layer, a two-dimensional vortex pair interacting with a wall, and a turbulent channel flow are selected as the test cases. We adopt a fourth-order compact scheme (COM4) for polymeric stress derivatives in the momentum equations. For convective derivatives in the constitutive equations, the first-order upwind difference scheme (UD) and artificial diffusion scheme (AD), which are commonly used in the literature, show most stable and smooth solutions even for highly extensional flows. However, the stress fields are smeared too much and the flow fields are quite different from those obtained by higher-order upwind difference schemes for the same flow parameters. Among higher-order upwind difference schemes, a third-order compact upwind difference scheme (CUD3) shows most stable and accurate solutions. Therefore, a combination of CUD3 for the convective derivatives in the constitutive equations and COM4 for the polymeric stress derivatives in the momentum equations is recommended to be used for numerical simulation of highly extensional flows.

• 대한기계학회논문집B
• /
• 제29권9호
• /
• pp.1049-1056
• /
• 2005

A conservative pressure-based finite-volume numerical method has been developed for computing flow and heat transfer by using an unstructured grid system. The method admits arbitrary convex polyhedra. Care is taken in the discretization and solution procedures to avoid formulations that are cell-shape-specific. A collocated variable arrangement formulation is developed, i.e. all dependent variables such as pressure and velocity are stored at cell centers. Gradients required for the evaluation of diffusion fluxes and for second-order-accurate convective operators are found by a novel second-order accurate spatial discretization. Momentum interpolation is used to prevent pressure checkerboarding and the SIMPLE algorithm is used for pressure-velocity coupling. The resulting set of coupled nonlinear algebraic equations is solved by employing a segregated approach, leading to a decoupled set of linear algebraic equations fer each dependent variable, with a sparse diagonally dominant coefficient matrix. These equations are solved by an iterative preconditioned conjugate gradient solver which retains the sparsity of the coefficient matrix, thus achieving a very efficient use of computer resources.

• 한국항공우주학회지
• /
• 제35권8호
• /
• pp.677-684
• /
• 2007

본 연구에서는 삼차원 점성 유동을 효율적으로 해석하기 위해 사면체, 프리즘, 피라미드를 포함하는 비정렬 혼합격자계를 기반으로 하는 유동해석코드를 개발하였다. 유동의 지배방정식은 격자점 중심의 유한체적법을 사용하여 공간차분회었으며, 제어테적은 메디안 듀얼(median-dual)방법으로 구성하였다. 난류유동 해석은 Spalart-Allmaras 난류모형과 연계하여 계산되었다. 개발된 해석코드의 정상 유동 검증을 위해 삼차원 날개에 대한 층류, 난류유동을 해석하였으며, 비정상 유동 검증을 위해 조화운동에 의해 진동하는 삼차원 날개에 대한 유동해석을 수행하였다.

• Structural Engineering and Mechanics
• /
• 제52권4호
• /
• pp.787-813
• /
• 2014

Starting with Hamilton's variational principle, the governing field equations for the steady state response of thin-walled beams under harmonic forces are derived. The formulation captures shear deformation effects due to bending and warping, translational and rotary inertia effects and as well as torsional flexural coupling effects due to the cross section mono-symmetry. The equations of motion consist of four coupled differential equations in the unknown displacement field variables. A general closed form solution is then developed for the coupled system of equations. The solution is subsequently used to develop a family of shape functions which exactly satisfy the homogeneous form of the governing field equations. A super-convergent finite element is then formulated based on the exact shape functions. Key features of the element developed include its ability to (a) isolate the steady state response component of the response to make the solution amenable to fatigue design, (b) capture coupling effects arising as a result of section mono-symmetry, (c) eliminate spatial discretization arising in commonly used finite elements, (d) avoiding shear locking phenomena, and (e) eliminate the need for time discretization. The results based on the present solution are found to be in excellent agreement with those based on finite element solutions at a small fraction of the computational and modelling cost involved.

• Nuclear Engineering and Technology
• /
• 제54권7호
• /
• pp.2625-2634
• /
• 2022

Discretization errors are extremely challenging conundrums of discrete ordinates calculations for radiation transport problems with void regions. In previous work, we have presented a multi-collision source method (MCS) to overcome discretization errors, but the efficiency needs to be improved. This paper proposes a goal-oriented algorithm for the MCS method to adaptively determine the partitioning of the geometry and dynamically change the angular quadrature in remaining iterations. The importance factor based on the adjoint transport calculation obtains the response function to get a problem-dependent, goal-oriented spatial decomposition. The difference in the scalar fluxes from one high-order quadrature set to a lower one provides the error estimation as a driving force behind the dynamic quadrature. The goal-oriented algorithm allows optimizing by using ray-tracing technology or high-order quadrature sets in the first few iterations and arranging the integration order of the remaining iterations from high to low. The algorithm has been implemented in the 3D transport code ARES and was tested on the Kobayashi benchmarks. The numerical results show a reduction in computation time on these problems for the same desired level of accuracy as compared to the standard ARES code, and it has clear advantages over the traditional MCS method in solving radiation transport problems with reflective boundary conditions.

• Nuclear Engineering and Technology
• /
• 제55권6호
• /
• pp.2172-2194
• /
• 2023

A variational nodal method (VNM) with unstructured-mesh is presented for solving steady-state and dynamic neutron diffusion equations. Orthogonal polynomials are employed for spatial discretization, and the stiffness confinement method (SCM) is implemented for temporal discretization. Coordinate transformation relations are derived to map unstructured triangular nodes to a standard node. Methods for constructing triangular prism space trial functions and identifying unique nodes are elaborated. Additionally, the partitioned matrix (PM) and generalized partitioned matrix (GPM) methods are proposed to accelerate the within-group and power iterations. Neutron diffusion problems with different fuel assembly geometries validate the method. With less than 5 pcm eigenvalue (k_eff) error and 1% relative power error, the accuracy is comparable to reference methods. In addition, a test case based on the kilowatt heat pipe reactor, KRUSTY, is created, simulated, and evaluated to illustrate the method's precision and geometrical flexibility. The Dodds problem with a step transient perturbation proves that the SCM allows for sufficiently accurate power predictions even with a large time-step of approximately 0.1 s. In addition, combining the PM and GPM results in a speedup ratio of 2-3.

• Journal of applied mathematics & informatics
• /
• 제32권5_6호
• /
• pp.585-598
• /
• 2014

In this paper we consider the nonlinear parabolic problems with mixed boundary condition. Under comparatively mild conditions of the coefficients related to the problem, we construct the discontinuous Galerkin approximation of the solution to the nonlinear parabolic problem. We discretize spatial variables and construct the finite element spaces consisting of discontinuous piecewise polynomials of which the semidiscrete approximations are composed. We present the proof of the convergence of the semidiscrete approximations in $L^{\infty}(H^1)$ and $L^{\infty}(L^2)$ normed spaces.