• Coupled systems mechanics
• /
• v.8 no.2
• /
• pp.147-168
• /
• 2019

This work aims to simulate three-dimensional heavy oil flow in a reservoir with heater-wells. Mass, momentum and energy balances, as well as correlations for rock and fluid properties, are used to obtain non-linear partial differential equations for the fluid pressure and temperature, and for the rock temperature. Heat transfer is simulated using a two-equation model that is more appropriate when fluid and rock have very different thermal properties, and we also perform comparisons between one- and two-equation models. The governing equations are discretized using the Finite Volume Method. For the numerical solution, we apply a linearization and an operator splitting. As a consequence, three algebraic subsystems of linearized equations are solved using the Conjugate Gradient Method. The results obtained show the suitability of the numerical method and the technical feasibility of heating the reservoir with static equipment.

• Journal of the Korean Society of Visualization
• /
• v.18 no.3
• /
• pp.35-41
• /
• 2020

We conduct numerical simulations of the interaction of a deformable structure with two-phase compressible flow. The finite volume method (FVM) is used to simulate fluid phenomena including a shock wave, a gas bubble, and the deformation of free surface. The deformation of a floating structure is computed with the finite element method (FEM). The compressible two-phase volume of fluid (VOF) method is used for the generation and development of a cavitation bubble, and the immersed boundary method (IBM) is used to impose the effect of the structure on the fluid domain. The result of the simulation shows the generation of a shock wave, and the expansion of the bubble. Also, the deformation of the structure due to the hydrodynamic loading by the explosion is identified.

• Structural Engineering and Mechanics
• /
• v.77 no.5
• /
• pp.651-659
• /
• 2021

In this study, a new constitutive model has been developed to predict the elastic behavior of plain weave textile composites, using the finite element (FE) method. The geometric conditions and basic assumptions of this model are based on the basics of a continuum theory developed for the plane curved composites. In this model, the mechanical properties of the weave region and pure matrix region is calculated separately and then imported for the FE analysis. This new constitutive model is used to implement the mechanical properties of weave region in the representative volume element (RVE). The constitutive relations are implemented as user-material subroutine code (UMAT) in ABAQUS® FE software. The results of FE analysis have been compared with experimental results and other data available in the literature. These comparisons confirmed the capability of the presented model for the prediction of effective elastic properties of plain weave fabric composites.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.3
• /
• pp.59-70
• /
• 1992

A Finite volume method for the computation of the two-dimensional, incompressible, steady, laminar Navier-Stokes equation is developed using a non-staggered grid system in a general curvilinear coordinate. The numerical pressure fluctuations, usually encountered when the non-staggered grid system is used, is suppressed by the momentum interpolation method. Flows around a NACA0012 foil section have been computed by the present method and the results show good agreements with other experimental and numerical ones.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.8
• /
• pp.677-684
• /
• 2007

In the Present Study, a 3-D viscous flow solver, based on unstructured hybrid meshses containing tetrahedra, prisms and pyramids, has been developed. A finite-volume discretization scheme is used for solving the compressible Navier-Stokes equations. A cell-vertex median dual volume is used for spatial discretization. The one-equation Spalart-Allmaras turbulence model has been adopted to evaluate the eddy viscosity. Validation were made by computing laminar and turbulent flows around a 3-D wing for steady flows and turbulent flows around an oscillating 3-D wing in harmonic motion for unsteady flows.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.27 no.4
• /
• pp.197-201
• /
• 2015

Experimental and numerical analysis has been carried out in this paper to understand correlation between volume and yield of cellular bulkhead. It was firstly confirmed from these results that the conditions and parameters considered in the finite element analysis were reasonable and realistic due to the fact that the yield loads determined by the two different methods were equivalent in actual. Based on this results, a series of intensive numerical analysis has been further performed and revealed that the yield load varied in direct proportion to the change in the volume of the cellular bulkhead.

• Nuclear Engineering and Technology
• /
• v.41 no.8
• /
• pp.1045-1052
• /
• 2009

A calculation model is developed to predict the transient free surface flow on the containment floor following a loss-of-coolant accident (LOCA) of pressurized water reactors (PWR) for the use of debris transport evaluation. The model solves the two-dimensional Shallow Water Equation (SWE) using a finite volume method (FVM) with unstructured triangular meshes. The numerical scheme is based on a fully explicit predictor-corrector method to achieve a fast-running capability and numerical accuracy. The Harten-Lax-van Leer (HLL) scheme is used to reserve a shock-capturing capability in determining the convective flux term at the cell interface where the dry-to-wet changing proceeds. An experiment simulating a sudden break of a water reservoir with L-shape open channel is calculated for validation of the present model. It is shown that the present model agrees well with the experiment data, thus it can be justified for the free surface flow with accuracy. From the calculation of flow field over the simplified containment floor of APR1400, the important phenomena of free surface flow including propagations and interactions of waves generated by local water level distribution and reflection with a solid wall are found and the transient flow rates entering the Holdup Volume Tank (HVT) are obtained within a practical computational resource.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.5
• /
• pp.735-745
• /
• 1997

In this study, the compressibility of resin was considered in filling analysis to account for the possible packing type flow. A numerical simulation program employing a hybrid finite element/finite difference scheme was developed to solve Hele-Shaw flow of the compressible viscous fluid at non-isothermal conditions. To advance the melt front, a control volume approach was adopted. Thin complex 3-D shapes of cavities, runners, and sprues were discretized by employing triangular, cylindrical and/or rectangular strip elements. Mass conservation was applied to each control volume to solve for the pressure distribution. Directly applying a constant mass flow rate at the inlet removes calculation of the apparent pressure boundary conditions, resulting in better simulation condition. The Cross model was used to model viscosity and the Tait equation was employed to represent density as a function of temperature and pressure. The validity of the developed program was verified through comparisons with available data in the literature and the effect of compressibility on the pressure distribution was discussed. To reduce computation time, 1-D and 2-D elements were used instead of applying triangular elements and the numerical results were compared to each other.

• Journal of the Korean Institute of Telematics and Electronics D
• /
• v.34D no.4
• /
• pp.1-7
• /
• 1997

This paper presents an analysis of the wave propagation in cross-type tunnels using the finite volume time domain (FVTD) method. Because the FVTD method is based on the volume intergrations of themaxwell's equations with respected to arbitrary shaped small polyhedron cells and the fields at every center point of the cells ar eassigned in a average fashion, the method can handle arbitrary boundary problems with inhomogeneous media. In this paper, the wave propagation in cross-type tunnels has been analyzed using the fVTD method with the PML (perfectly matched layer) absorbing bundary conditons, and the numerical results are verified with a set of experimental data.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.65-73
• /
• 1989

Continuum formulations for the expressions of dynamic energy release rates and computational methods for dynamic stress intensity factors are developed for the analysis of dynamic fracture problems subjected to stress wave loading. Explicit volume integral expressions for instantaneous dynamic energy release rates are derived by modeling virtual crack extensions with the dynamic Eulerian-Lagrangian kinematic description. In the finite element applications a finite region around a crack-tip is modeled by using quarter-point singular isoparametric elements, and the volume integrals are evaluated for each crack-tip element during virtual crack extensions while the singularity is maintained. It is shown that the use of the present method is more reliable and accurate for the dynamic fracture analysis than that of other path-independent integral methods when the effects of stress waves are significant.