• Tunnel and Underground Space
• /
• v.30 no.6
• /
• pp.573-590
• /
• 2020

This study presents the current status of DECOVALEX-2023 project Task G and our research results so far. Task G, named 'Safety ImplicAtions of Fluid Flow, Shear, Thermal and Reaction Processes within Crystalline Rock Fracture NETworks (SAFENET)' aims at developing a numerical method to simulate the fracture creation and propagation, and the coupled thermohydro-mechanical processes in fracture in crystalline rocks. The first research step of Task G is a benchmark simulation, which is designed for research teams to make their modelling codes more robust and verify whether the models can represent an analytical solution for displacements of a single rock fracture. We reproduced the mechanical behavior of rock and embedded single fracture using a three-dimensional grain-based distinct element model for the simulations. In this method, the structure of the rock was represented by an assembly of rigid tetrahedral grains moving independently of each other, and the mechanical interactions at the grains and their contacts were calculated using 3DEC. The simulation results revealed that the stresses induced along the embedded fracture in the model were relatively low compared to those calculated by stress analysis due to stress redistribution and constrained fracture displacements. The fracture normal and shear displacements of the numerical model showed good agreement with the analytical solutions. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated using various experiments in a further study.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.1
• /
• pp.117-124
• /
• 2010

The bird strike to small aircraft has not been an issue because of its low speed and usage as a private aircraft. So, the compliance of the bird strike regulation is limited to large fixed-wing aircraft such as the commuter category in FAR Part 23 and the civil aircraft in FAR Part 25, generally. However, the forecast of dramatic increasing of VLJ(Very Light Jet), the usage of a composite material for an aircraft structure and flight time of general aviation due to Air-taxi for the point to point transportation, would rise up the need of bird strike regulations and a safety enhancement in normal and utility categorized aircraft. In this study, the safety of bird strike to small aircraft wing leading edge made of a metal and a composite material were compared using the explicit finite element analysis.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.37 no.1
• /
• pp.1-7
• /
• 2017

In general, the large loads which are applied from explosion, impact, earthquake and wind at a short time caused the materials of structures to large deformations, rotations and strains locally. If such phenomena will be analyzed, hydrocodes which can be considered fluid-structure interaction under computational continuum mechanics are inevitably needed. Also, the explosion mechanism is so complicated, it is reasonable that the behaviors of structure are predicted through explosion analyses and experiment at the same time. But, unfortunately, it is true that explosion experiments are limited to huge cost, large experiment facilities and safety problems. Therefore, in this study, it is shown that the results of explosion analyses using the AUTODYN are agreed with those of existing explosion experiments for reinforced concrete slabs within reasonable error limits. And the explosion damage of the same reinforced concrete slab are assessed for quite different reinforcement arrangement spacings, concrete cover depths, and vertical reinforcements. From the explosion analyses, it is known that the more the ratio of slab thickness to reinforcement arrangement spacing is increased, and small-diameter reinforcements are used than large-diameter reinforcements on the same reinforcement ratio, and vertical reinforcements are used, the more the anti-knock capacities are improved.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.3
• /
• pp.622-628
• /
• 2016

As the suction pad-supporting bike carrier attached to a car may be subject to an excessive dynamic load due to random vibrations and centrifugal forces during driving, its structural safety is of great concern. To examine this, the finite-element method with a fluid-structure interaction should be used because the pressure on the pad bottom is changed in real time according to the fluctuations of the force or the moment applied on the pad. This method, however, has high computing costs in terms of modeling efforts and software expense. Moreover, the accuracy of computation is not easily guaranteed. Therefore, a new method combining the experiment and computation is proposed in this paper: the bottom pressure and contact area of the pad under varying loads was measured in real time and the acquired data are then used in the nonlinear elastic finite-element calculations. The computational and experimental results obtained with the product under development showed that the safety margin of the pad under the axial loading is relatively sufficient, whereas with an excessive rotational loading, the pad is vulnerable to separation or a local surface damage; hence, the safety margin may not be secured. The predicted contact behavior under the variation of the magnitude and type of the loading were in good agreement with the one from the experiment. The proposed analysis method in this study could be used in the design of similar vacuum pad systems.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.9
• /
• pp.597-604
• /
• 2016

A parallel algorithm of bi-conjugate gradient method was developed based on CUDA for parallel computation of the incompressible Navier-Stokes equations. The governing equations were discretized using splitting P2P1 finite element method. Asymmetric stenotic flow problem was solved to validate the proposed algorithm, and then the parallel performance of the GPU was examined by measuring the elapsed times. Further, the GPU performance for sparse matrix-vector multiplication was also investigated with a matrix of fluid-structure interaction problem. A kernel was generated to simultaneously compute the inner product of each row of sparse matrix and a vector. In addition, the kernel was optimized to improve the performance by using both parallel reduction and memory coalescing. In the kernel construction, the effect of warp on the parallel performance of the present CUDA was also examined. The present GPU computation was more than 7 times faster than the single CPU by double precision.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.127-136
• /
• 2022

In this work a multi-fidelity non-intrusive polynomial chaos (MF-NIPC) has been applied to a structural wind engineering problem in architectural design for the first time. In architectural design it is important to design structures that are safe in a range of wind directions and speeds. For this reason, the computational models used to design buildings and bridges must account for the uncertainties associated with the interaction between the structure and wind. In order to use the numerical simulations for the design, the numerical models must be validated by experi-mental data, and uncertainties contained in the experiments should also be taken into account. Uncertainty Quantifi-cation has been increasingly used for CFD simulations to consider such uncertainties. Typically, CFD simulations are computationally expensive, motivating the increased interest in multi-fidelity methods due to their ability to lev-erage limited data sets of high-fidelity data with evaluations of more computationally inexpensive models. Previous-ly, the multi-fidelity framework has been applied to CFD simulations for the purposes of optimization, rather than for the statistical assessment of candidate design. In this paper MF-NIPC method is applied to flow around a rectan-gular 5:1 cylinder, which has been thoroughly investigated for architectural design. The purpose of UQ is validation of numerical simulation results with experimental data, therefore the radius of curvature of the rectangular cylinder corners and the angle of attack are considered to be random variables, which are known to contain uncertainties when wind tunnel tests are carried out. Computational Fluid Dynamics (CFD) simulations are solved by a solver that employs the Finite Element Method (FEM) for two turbulence modeling approaches of the incompressible Navier-Stokes equations: Unsteady Reynolds Averaged Navier Stokes (URANS) and the Large Eddy simulation (LES). The results of the uncertainty analysis with CFD are compared to experimental data in terms of time-averaged pressure coefficients and bulk parameters. In addition, the accuracy and efficiency of the multi-fidelity framework is demonstrated through a comparison with the results of the high-fidelity model.

• Advances in aircraft and spacecraft science
• /
• v.9 no.5
• /
• pp.415-431
• /
• 2022

Secondary flows have a huge impact on losses generation in modern low pressure gas turbines (LPTs). At design point, the interaction of the blade profile with the end-wall boundary layer is responsible for up to 40% of total losses. Therefore, predicting accurately the end-wall flow field in a LPT is extremely important in the industrial design phase. Since the inlet boundary layer profile is one of the factors which most affects the evolution of secondary flows, the first main objective of the present work is to investigate the impact of two different inlet conditions on the end-wall flow field of the T106A, a well known LPT cascade. The first condition, labeled in the paper as C1, is represented by uniform conditions at the inlet plane and the second, C2, by a flow characterized by a defined inlet boundary layer profile. The code used for the simulations is based on the Discontinuous Galerkin (DG) formulation and solves the Reynolds-averaged Navier-Stokes (RANS) equations coupled with the Spalart Allmaras turbulence model. Secondly, this work aims at estimating the influence of viscosity and turbulence on the T106A end-wall flow field. In order to do so, RANS results are compared with those obtained from an inviscid simulation with a prescribed inlet total pressure profile, which mimics a boundary layer. A comparison between C1 and C2 results highlights an influence of secondary flows on the flow field up to a significant distance from the end-wall. In particular, the C2 end-wall flow field appears to be characterized by greater over turning and under turning angles and higher total pressure losses. Furthermore, the C2 simulated flow field shows good agreement with experimental and numerical data available in literature. The C2 and inviscid Euler computed flow fields, although globally comparable, present evident differences. The cascade passage simulated with inviscid flow is mainly dominated by a single large and homogeneous vortex structure, less stretched in the spanwise direction and closer to the end-wall than vortical structures computed by compressible flow simulation. It is reasonable, then, asserting that for the chosen test case a great part of the secondary flows details is strongly dependent on viscous phenomena and turbulence.

• Economic and Environmental Geology
• /
• v.57 no.1
• /
• pp.1-15
• /
• 2024

The Melanesian region in the western Pacific is dominated by complex plate tectonics, with the largest oceanic plateau, the OntongJava plateau, and a hotspot, the Caroline Islands. To better understand the complex geodynamics of the region, we estimate P- and S-velocity models and 𝛿 (V_P/V_S) model by using relative teleseismic travel times measured at seismometers on land and the seafloor. Our results show high-velocity anomalies in the subduction zones of the Melanesian region to a depth of about 400 km, which is thought to be subducting Solomon Sea, Bismarck, and Australian plates along plate boundaries. Along subduction zones, positive 𝛿 (V_P/V_S) anomalies are found, which may be caused by partial melting due to dehydration. A broad high-velocity anomaly is observed at 600 km depth below the Ontong-Java plateau, with a negative 𝛿 (V_P/V_S) anomaly. This is thought to be a viscous and dry remnant of the Pacific plate that subducted at 45-25 Ma, with a low volume of fluids due to dehydration for a long period in the mantle transition zone. Beneath the Caroline Islands, a strong low-velocity anomaly is obseved to a depth of 800 km and appears to be connected to the underside of the remnant Pacific plate in the mantle transition zone. This suggests that the mantle plume originating in the lower mantle has been redirected due to the interaction with the remnant Pacific plate and has reached its current location. The mantle plume also has a positive 𝛿 (V_P/V_S) anomaly, which is thought to be due to the influence of embedded fluids or partial melting. A high-velocity anomaly, interpreted as an effect of the thick lithosphere beneath the Ontong-Java plateau, is observed down to 300 km depth with a negative 𝛿 (V_P/V_S) anomaly, which likely indicate that little fluid remains in the melt residue accumulated in the lithosphere.