• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.237-255
• /
• 2013

In this paper, a meshfree-enriched finite element method (ME-FEM) is introduced for the large deformation analysis of nonlinear path-dependent problems involving contact. In linear ME-FEM, the element formulation is established by introducing a meshfree convex approximation into the linear triangular element in 2D and linear tetrahedron element in 3D along with an enriched meshfree node. In nonlinear formulation, the area-weighted smoothing scheme for deformation gradient is then developed in conjunction with the meshfree-enriched element interpolation functions to yield a discrete divergence-free property at the integration points, which is essential to enhance the stress calculation in the stage of plastic deformation. A modified variational formulation using the smoothed deformation gradient is developed for path-dependent material analysis. In the industrial metal forming problems, the mortar contact algorithm is implemented in the explicit formulation. Since the meshfree-enriched element shape functions are constructed using the meshfree convex approximation, they pose the desired Kronecker-delta property at the element edge thus requires no special treatments in the enforcement of essential boundary condition as well as the contact conditions. As a result, this approach can be easily incorporated into a conventional displacement-based finite element code. Two elasto-plastic problems are studied and the numerical results indicated that ME-FEM is capable of delivering a volumetric locking-free and pressure oscillation-free solutions for the large deformation problems in metal forming analysis.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.18 no.4
• /
• pp.26-33
• /
• 2019

This study is a thermal and flow analysis of Organic Rankine Cycle (ORC) pipe line for 250 kW grade waste gas heat recovery. We attempted to obtain the boundary condition data through the process design of the ORC, which can produce an electric power of 250 kW through the recovery of waste heat. Then, we conducted a simulation by using STAR-CCM+ to verify the model for the pipe line stream of the 250 kW class waste heat recovery system. Based on the results of the thermal and flow analyses of each pipe line applied to the ORC system, we gained the following conclusion. The pressure was relatively increased at the pipe outside the refracted part due to the pipe shape. Moreover, the heat transfer amount of the refrigerant gas line is relatively higher than that of the liquid line.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.12 no.1
• /
• pp.343-353
• /
• 2020

Based on the Precise Transfer Matrix Method (PTMM), the dynamic model is constructed to observe the vibration behaviour of cylindrical shells with variable thickness by solving a set of first-order differential equations. The free vibration of stiffened cylindrical shells with variable thickness can be obtained to compare with the exact solution and FEM results. The reliability of the present method of free vibration is well proved. Furthermore, the effect of thickness on the vibration responses of the cylindrical shell is also discussed. The acoustic response of immersed cylindrical shells is analyzed by a Pluralized Wave Superposition Method (PWSM). The sound pressure coefficient can be gained by collocating points along the meridian line to satisfy the Neumann boundary condition. The mode convergence analysis of the cylindrical shell is carried out to guarantee calculation precision. Also, the reliability of the present method on sound radiation is verified by comparing with experimental results and numerical results.

• Journal of the Korea Computer Graphics Society
• /
• v.20 no.4
• /
• pp.9-16
• /
• 2014

This paper proposes a novel method that couples fluid and deformation/fracture. Our method considers two interaction types: fluid-object interaction and fluid-fluid interaction. In fluid-fluid interaction, we simulate water and smoke separately and blend their velocities in the intersecting region depend on their densities. Our method separates projection process into two steps for each of water and smoke. This reduces the number of grid cells required for projection in order to optimize the number of iterations for convergence and improve stability of the simulation. In water projection step, smoke region regarded as the cells with Dirichlet boundary condition. The smoke projection step solves water region with Neumann boundary condition. To take care of fluid-object interaction, we make use of the fluid pressure to update velocities of the each of the mass points so that the object can deform or fracture. Although our method doesn't provide physically accurate results, the various examples show that our method generate appealing visuals with good performance.

• Journal of Korea Water Resources Association
• /
• v.38 no.7 s.156
• /
• pp.537-544
• /
• 2005

For the proper management of high pressurized gas storage caverns, analysis of groundwater flow field and inflow quantity according to the groundwater head, gas storage pressure and water curtain head should be performed. The finite element method has been widely used for the groundwater flow analysis surrounding underground storage cavern because it can reflect the exact shape of cavern. But the various simulations according to the change of design factors such as the width of water curtain, shape of cavern etc. are not easy when elements were set up. To overcome these limitations, two dimensional groundwater flow model is established based on the boundary element method which compute the unknown variable by using only the boundary shape and condition. For the exact computation of drainage rate into cavern, the model test is performed by using the exact solution and pre-developed finite element model. The test result shows that the model could be used as an alternative to finite element model when various flow simulations are needed to determine the optimizing cavern shape and arrangement of water curtain holes and so forth.

• Journal of the Korean Geotechnical Society
• /
• v.24 no.8
• /
• pp.149-160
• /
• 2008

Offshore pile foundations are prone to lateral soil pressures resulting from embankment construction for the reclamation on deepwater soft clay. Since the 1990s, offshore reclamation has actively progressed in Korea, connecting with the development of Songdo newtown, Incheon newport, and Busan newport representatively. Special attention has been given to lateral soil-structure interaction problems related to passively-loaded offshore pile foundations. Based on a plane strain large deformation finite element (LDFE) approach, this paper presents the results of investigation into undrained (short-term) and drained (long-term) behavior of passively-loaded offshore pile foundations. This study examines the effects of major factors, such as soil profile, pile head boundary condition, magnitude of embankment load, and average degree of consolidation. The results allowed quantification of differences in the magnitude of lateral soil pressure acting on the piles between undrained and drained phases.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.2
• /
• pp.197-215
• /
• 2015

Modeling water flow in variably saturated, porous media is important in many branches of science and engineering. Highly nonlinear relationships between water content and hydraulic conductivity and soil-water pressure result in very steep wetting fronts causing numerical problems. These include poor efficiency when modeling water infiltration into very dry porous media, and numerical oscillation near a steep wetting front. A one-dimensional finite element formulation is developed for the numerical simulation of variably saturated flow systems. First order backward Euler implicit and second order Crank-Nicolson time discretization schemes are adopted as a solution strategy in this formulation based on Picard and Newton iterative techniques. Five examples are used to investigate the numerical performance of two approaches and the different factors are highlighted that can affect their convergence and efficiency. The first test case deals with sharp moisture front that infiltrates into the soil column. It shows the capability of providing a mass-conservative behavior. Saturated conditions are not developed in the second test case. Involving of dry initial condition and steep wetting front are the main numerical complexity of the third test example. Fourth test case is a rapid infiltration of water from the surface, followed by a period of redistribution of the water due to the dynamic boundary condition. The last one-dimensional test case involves flow into a layered soil with variable initial conditions. The numerical results indicate that the Crank-Nicolson scheme is inefficient compared to fully implicit backward Euler scheme for the layered soil problem but offers same accuracy for the other homogeneous soil cases.

• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.4
• /
• pp.55-63
• /
• 1995

This paper describes the procedure to analyze the performance of the end-plated propeller(EPP) by a boundary integral method. The screw blade(SB) and end-plate(EP) are represented by a set of quadrilateral panels, where the source and normal dipole of uniform strength are distributed. The perturbation velocity potential, being the only unknown via the potential-based formulation, is determined by satisfying the flow tangency condition on the blade and the end-plate at the same time. The Kutta condition is satisfied through an iterative process by requiring the null pressure jump across the upper and lower sides of the trailing edges of both the SH and the EP. Sample calculations indicate that the EP increases the loading near the tip of the SB while spreading the trailing vortices along the trailing edge of the EP, thus avoiding the strong tip-vortex formation. Predicted performance of the EPP shows good correlations with the experimental results. The method is therefore considered applicable in designing and analyzing the EPP which may be an alternative for energy-saving propulsive devices.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.2
• /
• pp.79-87
• /
• 2013

In this paper, an Euler equation solver based on a Cartesian-grid method and non-uniform staggered grid system is applied to predict the ship motion response and added resistance in waves. Water, air, and solid domains are identified by a volume-fraction function for each phase and in each cell. For capturing the interface between air and water, the tangent of hyperbola for interface capturing (THINC) scheme is used with a weighed line interface calculation (WLIC) method. The volume fraction of solid body embedded in a Cartesian-grid system is calculated by a level-set based algorithm, and the body boundary condition is imposed by volume weighted formula. Added resistance is calculated by direct pressure integration on the ship surface. Numerical simulations for a Wigley III hull and an S175 containership in regular waves have been carried out to validate the newly developed code, and the ship motion responses and added resistances are compared with experimental data. For S175 containership, grid convergence test has been conducted to investigate the sensitivity of grid spacing on the motion responses and added resistances.

• Nuclear Engineering and Technology
• /
• v.56 no.1
• /
• pp.9-18
• /
• 2024

In-vessel retention through external reactor vessel cooling (IVR-ERVC) is a severe accident management (SAM) strategy that has been adopted and used in many nuclear reactors such as AP1000, APR1400, and light water reactor etc. Some reactor accidents have raised concerns about nuclear reactors among residents, leading to a decrease in residents' acceptability and many studies on SAM are being conducted. Experiments on IVR-ERVC are almost impossible due to its specificity, so fluid characteristics are analyzed through BALI experiments with similar condition. In this study, computational fluid dynamics (CFD) via Reynolds-averaged Navier-Stokes (RANS) and large eddy simulation (LES) for BALI experiments were performed. Steady-state CFD analysis was performed on three turbulence models, and SST k-ω model was in good agreement with the experimental measurement temperature within the maximum error range of 1.9%. LES CFD analysis was performed based on the RANS analysis results and it was confirmed that the temperature and wall heat flux for depth was consistent within an error range of 1.0% with BALI experiment. The LES CFD analysis results were compared with those of the Lagrangian-based solver. LES matched the temperature distribution better than SOPHIA, but SOPHIA calculated the position of boundary between stratified layer and convective layer more accurately. On the other hand, Lagrangian-based solver predicted several small eddy behaviors of the convective layer and LES predicted large vortex behavior. The vibration characteristics near the cooling part of the BALI experimental device were confirmed through Fast Fourier Transform (FFT) investigation. It was found that the power spectral density for pressure at least 10 times higher near the side cooling than near the top cooling.