• Geophysics and Geophysical Exploration
• /
• v.22 no.2
• /
• pp.56-61
• /
• 2019

In this paper, we studied efficient modeling strategies when we simulate the 3D time-domain acoustic wave propagation using a cell-based finite difference method which can handle the variations of both P-wave velocity and density. The standard finite difference method assigns physical properties such as velocities of elastic waves and density to grid points; on the other hand, the cell-based finite difference method assigns physical properties to cells between grid points. The cell-based finite difference method uses average physical properties of adjacent cells to calculate the finite difference equation centered at a grid point. This feature increases the computational cost of the cell-based finite difference method compared to the standard finite different method. In this study, we used additional memory to mitigate the computational overburden and thus reduced the calculation time by more than 30 %. Furthermore, we were able to enhance the performance of the modeling on several media with limited density variations by using the cell-based and standard finite difference methods together.

• Geophysics and Geophysical Exploration
• /
• v.6 no.2
• /
• pp.77-86
• /
• 2003

We designed a new time-domain, finite-difference, elastic wave modeling technique, based on a displacement formulation. which yields nearly correct solutions to Lamb's problem. Unlike the conventional, displacement-based, finite-difference method using a node-based grid set (where both displacements and material properties such as density and Lame constants are assigned to nodal points), in our new finite-difference method, we use a cell-based grid set (where displacements are still defined at nodal points but material properties within cells). In the case of using the cell-based grid set, stress-free conditions at the free surface are naturally described by the changes in the material properties without any additional free-surface boundary condition. Through numerical tests, we confirmed that the new second-order finite differences formulated in the cell-based grid let generate numerical solutions compatible with analytic solutions unlike the old second-order finite-differences formulated in the node-based grid set.

• Korean Journal of Metals and Materials
• /
• v.47 no.11
• /
• pp.699-706
• /
• 2009

Deformation behavior of pure aluminum during equal channel angular pressing (ECAP) was simulated using a three-dimensional version of the finite element method in conjunction with a constitutive model based on the dislocation density and cell evolution. The three-dimensional finite element analyses for the prediction of microstructural features, such as the variation of the dislocation density and the cell size with the number of ECAP, are reported. The calculated stress and strain and their distributions are also investigated for the route Bc ECAP processed pure aluminum. The results of finite element analyses are found to be in good agreement with experimental results for the dislocation cell size. Due to the accumulation of strain throughout the workpiece and an overall trend to saturation in cell size, a decrease of the difference in cell size with the number of passes (1~4) was predicted.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.214-219
• /
• 2007

Both the bubble rising in a fully filled container and the droplet splash are simulated by a solution code(PowerCFD). This code employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with interface capturing method (CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The present results are compared with other numerical solutions found in the literature. It is found that the present code simulate complex free surface flows such as multi phase flows due to large density difference efficiently and accurately.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.6
• /
• pp.443-453
• /
• 2009

The Rayleigh-Taylor instability, the bubble rising in both partially and fully filled containers and the droplet splash are simulated by an in-house solution code(PowerCFD), which are typical benchmark problems among multiphase flows with material interface due to density difference. The present method(code) employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with interface capturing method(CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The present results are compared with other numerical solutions found in the literature. It is found that the present method simulates efficiently and accurately complex free surface flows such as multiphase flows with material interface due to both density difference and instability.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.572-575
• /
• 2008

Two-dimensional multiphase flows due to density difference such as the Rayleigh-Taylor instability problem and the droplet splash are simulated by an in-house solution code(PowerCFD). This code employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with interface capturing method in a volume of fluid(VOF) scheme for phase interface capturing. The present results are compared with other numerical solutions found in the literature. It is found that the present code simulates complex free surface flows such as multiphase flows due to density difference efficiently and accurately.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.40 no.5
• /
• pp.405-412
• /
• 2012

A flow control on airfoil installed a vortex cell for high efficiency wind turbine blade in stationary and dynamic stall conditions have been numerically investigated by solving the compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with Roe's flux-difference splitting and an implicit time-integration method coupled with dual time step sub-iteration. The computed result for the airfoil in the stationary showed that lift-drag ratio increases due to low pressure by the vortex cell. The oscillating airfoil with the vortex cell showed that the magnitude of hysteresis loop is reduced due to the enhanced vortex in the cell.

• Nuclear Engineering and Technology
• /
• v.34 no.3
• /
• pp.250-258
• /
• 2002

This Paper is to apply the well-established coarse mesh finite difference(CMFD) method to the method of characteristics(MOC) transport calculation as an acceleration scheme. The CMFD problem is first formulated at the pin-cell level with the multi-group structure To solve the cell- based multi-group CMFD problem efficiently, a two-group CMFD formulation is also derived from the multi-group CMFD formulation. The performance of the CMFD acceleration is examined for three test problems with different sizes including a realistic quarter core PWR problem. The CMFD formulation provides a significant reduction in the number of ray tracings and thus only about 9 ray tracing iterations are enough for the realistic problem. In computing time, the CMFD accelerated case is about two or three times faster than the coarse-mesh rebalancing(CMR) accelerated case.

• Journal of Ocean Engineering and Technology
• /
• v.17 no.6
• /
• pp.7-15
• /
• 2003

A Digital Wave Tank simulation technique, based on a finite-difference method and a modified marker-and-cell (MAC) algorithm, is applied in order to investigate the characteristics of nonlinear Tsunami propagations and their interactions with a 2D sloping beach, Ohkushiri Island, and to predict maximum wove run-up around the island. The Navier-Stokes (NS) and continuity equation are governed in the computational domain, and the boundary values are updated at each time step, by a finite-difference time-marching scheme in the frame of the rectangular coordinate system. The fully nonlinear, kinematic, free-surface condition is satisfied by the modified marker-density function technique. The near shore Tsunami is assumed to be a solitary wave, and is generated from the numerical wave-maker in the developed Digital Wave Tank. The simulation results are compared with the experiments and other numerical methods, based on the shallow-water wave theory.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2003.05a
• /
• pp.231-239
• /
• 2003

A Digital Wave Tank simulation technique based on a finite-difference method and a modified marker-and-cell (MAC) algorithm is applied to investigate the characteristics of nonlinear Tsunami propagations and their interactions with a 2D sloping beach and Ohkushiri island, and to predict maximum wave run-up around the island. The Navier-Stokes (NS) and continuity equation are governed in the computational domain and the boundary values updated at each time step by a finite-difference time-marching scheme in the frame of rectangular coordinate system. The fully nonlinear kinematic free-surface condition is satisfied by the modified marker-density function technique. The Nearshore Tsunami is assumed to be a solitary wave and generated from the numerical wavemaker in the developed Digital Wave Tank. The simulation results are compared with the experiments and other numerical methods based on the shallow-water wave theory.