• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2009.05a
• /
• pp.173-176
• /
• 2009

The Arbitrary Lagrangian-Eulerian(ALE, in short) method is the new description of continum motion, which combines the advantages of the classical kinematical descriptions, i.e. Lagrangian and Eulerian description, while minimizing their respective drawbacks. In this paper, the ALE description is adapted to simulate fluid-structure interaction problems. An automatic re-mesh algorithm and a fluid-structure coupling process are included to analyze the interaction and moving motion during the 2-D axisymmetric solid rocket interior FSI phenomena simulation.

• Coupled systems mechanics
• /
• v.7 no.6
• /
• pp.649-668
• /
• 2018

In this paper, we present a numerical model for fluid-structure interaction between structure built of porous media and acoustic fluid, which provides both pore pressure inside porous media and hydrodynamic pressures and hydrodynamic forces exerted on the upstream face of the structure in an unified manner and simplifies fluid-structure interaction problems. The first original feature of the proposed model concerns the structure built of saturated porous medium whose response is obtained with coupled discrete beam lattice model, which is based on Voronoi cell representation with cohesive links as linear elastic Timoshenko beam finite elements. The motion of the pore fluid is governed by Darcy's law, and the coupling between the solid phase and the pore fluid is introduced in the model through Biot's porous media theory. The pore pressure field is discretized with CST (Constant Strain Triangle) finite elements, which coincide with Delaunay triangles. By exploiting Hammer quadrature rule for numerical integration on CST elements, and duality property between Voronoi diagram and Delaunay triangulation, the numerical implementation of the coupling results with an additional pore pressure degree of freedom placed at each node of a Timoshenko beam finite element. The second original point of the model concerns the motion of the outside fluid which is modeled with mixed displacement/pressure based formulation. The chosen finite element representations of the structure response and the outside fluid motion ensures for the structure and fluid finite elements to be connected directly at the common nodes at the fluid-structure interface, because they share both the displacement and the pressure degrees of freedom. Numerical simulations presented in this paper show an excellent agreement between the numerically obtained results and the analytical solutions.

• Geomechanics and Engineering
• /
• v.28 no.6
• /
• pp.613-624
• /
• 2022

A study of the evolution of overburden fractures under the solid-fluid coupling state was conducted based on the geological and mining characteristics of the coal seam depth, weak strata cementation, and high-intensity mining in the mining areas of West China. These mining characteristics are key to achieving water conservation during mining or establishing groundwater reservoirs in coal mines. Based on the engineering background of the Daliuta Coal Mine, a non-hydrophilic simulation material suitable for simulating the weakly cemented rock masses in this area was developed, and a physical simulation test was carried out using a water-sand gushing test system. The study explored the spatial distribution and dynamic evolution of the fractured zone in the mining overburden under the coupling of stress and seepage. The experimental results show that the mining overburden can be vertically divided into the overall migration zone, the fracture extension zone and the collapse zone; additionally, in the horizontal direction, the mining overburden can be divided into the primary fracture zone, periodic fracture zone, and stop-fracture zone. The scope of groundwater flow in the overburden gradually expands with the mining of coal seams. When a stable water inrush channel is formed, other areas no longer generate new channels, and the unstable water inrush channels gradually close. Finally, the primary fracture area becomes the main water inrush channel for coal mines. The numerical simulation results indicate that the overlying rock breaking above the middle of the mined-out area allows the formation of the water-conducting channel. The water body will flow into the fracture extension zone with the shortest path, resulting in the occurrence of water bursting accidents in the mining face. The experimental research results provide a theoretical basis for the implementation of water conservation mining or the establishment of groundwater reservoirs in western mining areas, and this theoretical basis has considerable application and promotion value.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.709-714
• /
• 2000

A splitting method for the direct numerical simulation of solid-liquid mixtures is presented, where a symmetric pressure equation is newly proposed. Through numerical experiment, it is found that the newly proposed splitting method works well with a matrix-free formulation fer some bench mark problems avoiding an erroneous pressure field which appears when using the conventional pressure equation of a splitting method. When deriving a typical pressure equation of a splitting method, the motion of a solid particle has to be approximated by the 'intermediate velocity' instead of treating it as unknowns since it is necessary as a boundary condition. Therefore, the motion of a solid particle is treated in such an explicit way that a particle moves by the known form drag (pressure drag) that is calculated from the pressure equation in the previous step. From the numerical experiment, it was shown that this method gives an erroneous pressure field even for the very small time step size as a particle velocity increases. In this paper, coupling the unknowns of particle velocities in the pressure equation is proposed, where the resulting matrix is reduced to the symmetric one by applying the projector of the combined formulation. It has been tested over some bench mark problems and gives reasonable pressure fields.

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.396-401
• /
• 2008

The disc-pad brake system is an important part of automobile safety system. During braking, the kinetic energy and potential energies of a moving vehicle are converted into the thermal energy through frictional heat between the brake disc and the pads. Most of the thermal energy dissipated through the brake disc. The temperature could be exceed the critical value for a given material, which leads to undesirable effects, such as the brake fade, premature wear, brake fluid vaporization, bearing failure, thermal cracks, and thermallyexcited vibration. The object of the present study is to investigate temperature field and temperature variation of brake disc and pad during single brake. The brake disc is decelerated at the initial speed with constant acceleration, until the disc comes to stop. The pad-disc brake assembly is built by 3D model with the appropriate boundary condition. In the simulation process, the mechanical loads are applied to the thermomechanical coupling analysis in order to simulate the process of heat produced by friction.

• Wind and Structures
• /
• v.37 no.4
• /
• pp.275-287
• /
• 2023

Subgrade differential settlement of high-speed railways was a pivotal issue that could increase the risk of trains operation. The risk will be further increased when trains in the subsidence zone are affected by crosswinds. In this paper, the computational fluid dynamics (CFD) model and finite element (FE) model were established, and the data transmission interface of the two models was established by fluid-solid interaction (FSI) method to form a systematic crosswind-train-track-subgrade dynamic model. The risk of high-speed train encountering crosswind in settlement area was analyzed. The results showed that the aerodynamic force of the trains increased significantly with the increase in crosswind speed. The aerodynamic force of the trains could reach 125.14 kN, significantly increasing the risk of derailment and overturning. Considering the influence of crosswind, the risk of train operation could be greatly increased. The safety indices and the wheel-rail force both increased with the increase of the wind speed. For the high-speed train running at 350 km/h, the warning value of wind speed was 10.2 m /s under the condition of subgrade settlement with wavelength of 20 m and amplitude of 15 mm.

• Journal of Ocean Engineering and Technology
• /
• v.27 no.2
• /
• pp.53-58
• /
• 2013

Recently, some fluid-structure interaction (FSI) problems involving the fluid impact loads interacting with structures, such as sloshing, slamming, green-water, etc., have been considered, especially in the ocean engineering field. The governing equations for both an elastic solid model and flow model were originally derived from similar continuum mechanics principles. In this study, an elastic model based on a particle method, the MPS method, was developed for simulating the FSI problems. The developed model was first applied to a simple cantilever deflection problem for verification. Then, the model was coupled with the fluid flow model, the PNU (Pusan National University modified)-MPS method, and applied to the numerical investigation of the coupling effects between a cantilever and a mass of water, which has variable density, free-falling to the end of the cantilever.

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

The aim of this study was to develop a new efficient strategy that uses the Vector form Intrinsic Finite-element (VFIFE) method to conduct the static and dynamic analyses of marine pipes. Nonlinear problems, such as large displacement, small strain, and contact and collision, can be analyzed using a unified calculation process in the VFIFE method according to the fundamental theories of point value description, path element, and reverse motion. This method enables analysis without the need to integrate the stiffness matrix of the structure, because only motion equations of particles established according to Newton's second law are required. These characteristics of the VFIFE facilitate the modeling and computation efficiencies in analyzing the nonlinear dynamic problem of flexible pipe with large deflections. In this study, a three-dimensional (3-D) dynamical model based on 3-D beam element was established according to the VFIFE method. The deep-sea flexible pipe was described by a set of spatial mass particles linked by 3-D beam element. The motion and configuration of the pipe are determined by these spatial particles. Based on this model, a simulation procedure to predict the 3-D dynamical behavior of flexible pipe was developed and verified. It was found that the spatial configuration and static internal force of the mining pipe can be obtained by calculating the stationary state of pipe motion. Using this simulation procedure, an analysis was conducted on the static and dynamic behaviors of the flexible mining pipe based on a 1000-m sea trial system. The results of the analysis proved that the VFIFE method can be efficiently applied to the static and dynamic analyses of marine pipes.

• Journal of Ocean Engineering and Technology
• /
• v.31 no.1
• /
• pp.14-21
• /
• 2017

The erosion of solid particles in a pipe elbow was numerically investigated. A numerical procedure to estimate the sand erosion rate, as well as the particle motion, in the pipe elbow flow was introduced. This procedure was performed based on the combined empirical erosion model and computational fluid dynamics (CFD) analysis to consider the interaction between the particle motion and the eroded surface. The underlying turbulent flow on an Eulerian frame is described by the Reynolds averaged Navier-Stokes (RANS) equations with a $k-{\epsilon}$ turbulent model. The one-way coupled Eulerian-Lagrangian motion of the air flow and sand particles is employed to simulate the particle trajectories and particle-wall interactions on the pipe surfaces. The predicted CFD erosion magnitudes are compared with experimental data from pipe elbows. The erosion rate results do not reveal a good accordance between the simulation and experimental results. It seems that the CFD shows a slightly over-predicted erosion ratio.

• Nuclear Engineering and Technology
• /
• v.54 no.3
• /
• pp.849-859
• /
• 2022

Feed-water valves (FWVs) are used to regulate the flow rate of water entering steam generators, which are very important devices in nuclear power plants. Due to the working environment of relatively high pressure and temperature, there is strength failure problem of valve body in some cases. Based on the thermo-fluid-solid coupling model, the valve body stress of the feed-water valve in the opening process is investigated. The flow field characteristics inside the valve and temperature change of the valve body with time are studied. The stress analysis of the valve body is carried out considering mechanical stress and thermal stress comprehensively. The results show that the area with relatively high-velocity area moves gradually from the bottom of the cross section to the top of the cross section with the increase of the opening degree. The whole valve body reaches the same temperature of 250 ℃ at the time of 1894 s. The maximum stress of the valve body meets the design requirements by stress assessment. This work can be referred for the design of FWVs and other similar valves.