• 한국지구물리탐사학회:학술대회논문집
• /
• 2009.10a
• /
• pp.75-80
• /
• 2009

Despite its flexibility to complex geometry, three-dimensional (3D) electromagnetic(EM) modeling schemes using finite element method (FEM) have been faced to practical limitation due to the resulting large system of equations to be solved. An efficient 3D FEM modeling scheme has been developed, which can adopt either direct or iterative solver depending on the problems. The direct solver PARDISO can reduce the computing time remarkably by incorporating parallel computing on multi-core processor systems, which is appropriate for single frequency multi-source configurations. When limited memory, the iterative solver BiCGSTAB(1) can provide fast and stable convergence. Efficient 3D simulations can be performed by choosing an optimum solver depending on the computing environment and the problems to be solved. This modeling includes various types of controlled-sources and can be exploited as an efficient engine for 3D inversion.

• Journal of computational fluids engineering
• /
• v.20 no.1
• /
• pp.32-38
• /
• 2015

A 3-dimensional compressible turbulent boundary layer solver has been developed. A time marching method is used to integrate the turbulent boundary layer equations. While the direct integration of the boundary layer equations is performed for unseparated flow regions, the inverse integration is performed for separated flow regions. The program is verified for flows that have analytical solutions or other numerical results. The solver will be merged with an Euler solver for viscous-inviscid interaction.

• Structural Engineering and Mechanics
• /
• v.30 no.4
• /
• pp.467-480
• /
• 2008

An aggregation multigrid method (AMM) is a leading iterative solver in solid mechanics. Recently, AMM is applied for solving Schur Complement system in the FE analysis of shell structures. In this work, an extended application of AMM for solving Schur Complement system in the FE analysis of continuum elements is presented. Further, the performance of the proposed AMM in multiple load cases, which is a challenging problem for an iterative solver, is studied. The proposed method is developed by combining the substructuring and the multigrid methods. The substructuring method avoids factorizing the full-size matrix of an original system and the multigrid method gives near-optimal convergence. This method is demonstrated for the FE analysis of several elastostatic problems. The numerical results show better performance by the proposed method as compared to the preconditioned conjugate gradient method. The smaller computational cost for the iterative procedure of the proposed method gives a good alternative to a direct solver in large systems with multiple load cases.

• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.40-48
• /
• 2014

For rarefied gas flow regimes, physical phenomena such as velocity slip and temperature jump occur on the solid body surface. To predict these phenomena accurately, either the Navier-Stokes solver with a slip boundary condition or the direct simulation Monte Carlo method should be used. In the present study, flow simulations of a wedge were conducted in Mach-10 flow of argon gas for several different flow regimes using a two-dimensional Navier-Stokes solver with the Maxwell slip boundary condition. The results of the simulations were compared with those of the direct simulation Monte Carlo method to assess the present method. It was found that the values of the velocity slip and the temperature jump predicted increase as the Knudsen number increases. Also, the results are comparatively reasonable up to the Knudsen number of 0.05.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.6
• /
• pp.429-445
• /
• 2008

Two-phase compressible flow fields of air-water are investigated numerically in the fixed Eulerian grid framework. The phase interface is captured via volume fractions of each phase. A way to model two phase compressible flows as a single phase one is found based on an equivalent equation of states of Tait's type for a multiphase cell. The equivalent single phase field is discretized using the Roe‘s approximate Riemann solver. Two approaches are tried to suppress the pressure oscillation phenomena at the phase interface, a passive advection of volume fraction and a direct pressure relaxation with the compressible form of volume fraction equation. The direct pressure equalizing method suppresses pressure oscillation successfully and generates sharp discontinuities, transmitting and reflecting acoustic waves naturally at the phase interface. In discretizing the compressible form of volume fraction equation, phase interfaces are geometrically reconstructed to minimize the numerical diffusion of volume fraction and relevant variables. The motion of a projectile in a water-filled tube which is fired by the release of highly pressurized air is simulated presuming the flow field as a two dimensional one, and several design factors affecting the projectile movement are investigated.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1310-1317
• /
• 2017

GASFLOW-MPI is a widely used scalable computational fluid dynamics numerical tool to simulate the fluid turbulence behavior, combustion dynamics, and other related thermal-hydraulic phenomena in nuclear power plant containment. An efficient scalable linear solver for the large-scale pressure equation is one of the key issues to ensure the computational efficiency of GASFLOW-MPI. Several advanced Krylov subspace methods and scalable preconditioning methods are compared and analyzed to improve the computational performance. With the help of the powerful computational capability, the large eddy simulation turbulent model is used to resolve more detailed turbulent behaviors. A backward-facing step flow is performed to study the free shear layer, the recirculation region, and the boundary layer, which is widespread in many scientific and engineering applications. Numerical results are compared with the experimental data in the literature and the direct numerical simulation results by GASFLOW-MPI. Both time-averaged velocity profile and turbulent intensity are well consistent with the experimental data and direct numerical simulation result. Furthermore, the frequency spectrum is presented and a -5/3 energy decay is observed for a wide range of frequencies, satisfying the turbulent energy spectrum theory. Parallel scaling tests are also implemented on the KIT/IKET cluster and a linear scaling is realized for GASFLOW-MPI.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.10
• /
• pp.947-954
• /
• 2009

Starting from the Eu's GH(Generalized Hydrodynamic) theory, the multi-species GH numerical solver is developed in this research and its computatyional behaviors are examined for the hypersonic rarefied flow over an axisymmetric body. To improve the accuracy of the developed multi-species GH solver, various slip-wall boundary conditions are tested and the computed results are compared. Additionally, in order to validate the entropy characteristics of the GH equation, the entropy production and entropy generation rates of the GH equation are investigated in the 1-dimensional normal shock structure test at a high Knudsen number.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.3 no.2
• /
• pp.161-171
• /
• 1999

Aerodynamic sensitivity analysis codes are developed via the hand-differentiation using a direct differentiation method and an adjoint method respectively from discrete two-dimensional compressible Navier-Stokes equations. Unlike previous other researches, Baldwin-Lomax algebraic turbulence model is also differentiated by hand to obtain design sensitivities with respect to design variables of interest in turbulent flows. Discrete direct sensitivity equations and adjoint equations are efficiently solved by the same time integration scheme adopted in the flow solver routine. The required memory for the adjoint sensitivity code is greatly reduced at the cost of the computational time by allowing the large banded flux jacobian matrix unassembled. Direct sensitivity code results are found to be exactly coincident with sensitivity derivatives obtained by the finite difference. Adjoint code results of a turbulent flow case show slight deviations from the exact results due to the limitation of the algebraic turbulence model in implementing the adjoint formulation. However, current adjoint sensitivity code yields much more accurate sensitivity derivatives than the adjoint code with the turbulence eddy viscosity being kept constant, which is a usual assumption for the prior researches.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.11 no.3
• /
• pp.281-291
• /
• 1999

A computer code based on a node equation method has been developed for the analysis of pipe network. Both data structure and object-oriented programming technique are used for pipe and node modelling, in which simplification process is applied to complicated and large pipe network. The semi-direct solver, ILU-CGS, improves greatly both the accuracy and the rate of convergence. The computational result of high-pressure pipe network of city gas in Taegu shows the good agreement with the real data.

• International Journal of Fluid Machinery and Systems
• /
• v.4 no.1
• /
• pp.161-171
• /
• 2011

The availability of a high quality open source CFD simulation platform like OpenFOAM offers new R&D opportunities by providing direct access to models and solver implementation details. Efforts have been made by Hydro-Qu$\'{e}$bec to adapt OpenFOAM to hydroturbines for the development of steady-state capabilities. This paper describes the developments that have been made to implement new turbomachinery related capabilities: multiple frames of reference solver, domain coupling interfaces (GGI, cyclicGGI and mixing plane) and specialized boundary conditions. Practical use of the new turbomachinery capabilities are demonstrated for the analysis of a 195-MW Francis turbine.