• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.494-498
• /
• 2001

When supersonic jet impinges on wall from the nozzle, complex flow pattern appears such as Mach disc, expansion fan, and jet boundary. The numerical computation of this supersonic jet is important on flame deflecctor design for launch space especially. In this paper, we analyzed supersonic jet structure impinging on deflector wall using three dimensional steady and unsteady compressible equation and showed temperature and pressure distribution on the wall surface. As a result, some dominant factors of jet flows are discussed for conceptual design of flame deflector.

• Journal of computational fluids engineering
• /
• v.17 no.3
• /
• pp.25-32
• /
• 2012

This work was devoted to compare two different cavitation models to study the dependency of model constants. The cavitation model of Merkle et al.(2006) and Kunz et al.(2000) were used for the present computational study. The cavitation models were coupled with the incompressible unsteady Reynolds-Averaged Navier-Stokes solver to indicate the vaporization and condensation processes. For this purpose, a preconditioning method was added as the pseudo-time term to solve the unsteady stiffness problems. For the validation of the numerical simulation, the computation was performed for the cavitating flow in a converging-diverging channel. The present results show that Merkle's cavitation model is independent to the model constants, and the higher numerical accuracy over Kunz's cavitation model.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.3
• /
• pp.716-724
• /
• 1990

Oblique impinging plane jets were investigated experimentally and numerically at Reynolds number 21000. The inclination angle was varied from 90.deg.(normal to the impinging plate) to 60.deg.. The distance H between the nozzle exit and the stagnation point on the impinging plate was fixed at H/D=8. The working fluid was air. The mean velocity components and turbulent quantities were measured by a hot-wire anemometer. And the static pressure distributions on the impinging plate were measured by a Pitot tube. In numerical computation, the governing partial differential equations of elliptic type were solved with conventional k-.epsilon. turbulence model. The measurements show that, after impingement, the jet half width alone the wall increases in both directions, and that similarity for each turbulent quantity such as Reynolds shear stress or turbulent kinetic energy is revealed in the wall jet region. The computed results show some deviation from experimental data in the impingement region, where streamline curvature is significant. However, the computed results agree qualitatively well with measurements.

• Earthquakes and Structures
• /
• v.10 no.3
• /
• pp.495-521
• /
• 2016

Ancient monuments of Greek and Roman classical architecture usually consist of multi-drum columns that are constructed of stone blocks placed on top of each other. Several research studies deal with the seismic behaviour of such structures, since earthquakes are common causes of destruction of such monuments. This paper investigates the effect of multiple earthquakes on the seismic performance of multi-drum columns, through numerical simulations and parametric analyses. The Discrete Element Method and an appropriate contact model have been implemented in a specially developed software application that is able to efficiently perform the necessary simulations in two dimensions. Specifically, various strong ground excitations are used in series for the computation of the collective final deformation of multi-drum columns. In order to calculate this cumulative deformation for a series of ground motions, the individual deformation of the column for each excitation is computed and then used as initial conditions for the next earthquake excitation. Various multi-drum columns with different dimensions are also considered in the analyses in order to examine how the geometric characteristics of columns can affect their seismic sequence behaviour, in combination with the excitation frequency content.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.98-104
• /
• 2011

Numerical simulation in exploration geophysics provides important insights into subsurface wave propagation phenomena. Although elastic wave simulations take longer to compute than acoustic simulations, an elastic simulator can construct more realistic wavefields including shear components. Therefore, it is suitable for exploration of the responses of elastic bodies. To overcome the long duration of the calculations, we use a Graphic Processing Unit (GPU) to accelerate the elastic wave simulation. Because a GPU has many processors and a wide memory bandwidth, we can use it in a parallelised computing architecture. The GPU board used in this study is an NVIDIA Tesla C1060, which has 240 processors and a 102 GB/s memory bandwidth. Despite the availability of a parallel computing architecture (CUDA), developed by NVIDIA, we must optimise the usage of the different types of memory on the GPU device, and the sequence of calculations, to obtain a significant speedup of the computation. In this study, we simulate two- (2D) and threedimensional (3D) elastic wave propagation using the Finite-Difference Time-Domain (FDTD) method on GPUs. In the wave propagation simulation, we adopt the staggered-grid method, which is one of the conventional FD schemes, since this method can achieve sufficient accuracy for use in numerical modelling in geophysics. Our simulator optimises the usage of memory on the GPU device to reduce data access times, and uses faster memory as much as possible. This is a key factor in GPU computing. By using one GPU device and optimising its memory usage, we improved the computation time by more than 14 times in the 2D simulation, and over six times in the 3D simulation, compared with one CPU. Furthermore, by using three GPUs, we succeeded in accelerating the 3D simulation 10 times.

• Journal of the Korean Geotechnical Society
• /
• v.35 no.12
• /
• pp.35-44
• /
• 2019

In this study, the axial stiffness of reinforcing piles (K_vr) for the vertical extension remodeling structures was estimated through 3D finite element analysis. In the computation of the minimum required axial stiffness of reinforcing piles, proposed maximum axial stiffness of old and deteriorated existing piles (K_ve) based on theoretical and experimental approaches will be applied. Through this, the required increase rate of axial stiffness of reinforcing piles in order to support the increased structural loading was proposed for end-bearing and friction piles by different slenderness ratio (L/D). The numerical model was validated by comparing the computed results with actual field measurements. Based on the computed results, it was concluded that the end-bearing reinforcing pile needs 44% - 67% increase in axial stiffness to deal with the deterioration of existing piles and support the additional structural load due to vertical extension remodeling.

• Journal of computational fluids engineering
• /
• v.13 no.2
• /
• pp.27-41
• /
• 2008

An unstructured hybrid mesh flow solver has been developed for the simulation of three-dimensional steady and unsteady incompressible flow fields. The incompressible Navier-Stokes equations with an artificial compressibility method were discretized by using a node-based finite-volume method. For the unsteady time-accurate computation, a dual-time stepping method was adopted to satisfy a divergence-free flow field at each physical time step. An implicit time integration method with local time stepping was implemented to accelerate the convergence in the pseudo-time sub-iteration procedure. The one-equation Spalart-Allmaras turbulence model has been adopted to solve high-Reynolds number flow fields. The flow solver was parallelized to minimize the CPU time and to overcome the computational overhead. This method has been applied to calculate steady and unsteady flow fields around submarine configurations and a 3-D infinite cylinder. Validations were made by comparing the predicted results with those of experiments or other numerical results. It was demonstrated that the present method is efficient and robust for the prediction of steady and unsteady incompressible flow fields.

• Wind and Structures
• /
• v.10 no.6
• /
• pp.511-532
• /
• 2007

The present study aims to generate turbulent inflow data to more accurately represent the turbulent flow around a square cylinder when the inflow turbulence level is significant. The modified random flow generation (RFG) technique in conjunction with a previously developed LES code is successfully adopted into a finite element based fluid flow solver to generate the required inflow turbulence boundary conditions for the three-dimensional (3-D) LES computations of transitional turbulent flow around a square cylinder at Reynolds number of 22,000. The near wall region is modelled without using wall approximate conditions and a wall damping coefficient is introduced into the calculation of sub-grid length scale in the boundary layer of the cylinder wall. The numerical results obtained from simulations are compared with each other and with the experimental data for different inflow turbulence boundary conditions in order to discuss the issues such as the synthetic inflow turbulence effects on the 3-D transitional flow behaviour in the near wake and the free shear layer, the basic mechanism by which stream turbulence interacts with the mean flow over the cylinder body and the prediction of integral flow parameters. The comparison among the LES results with and without inflow turbulence and the experimental data emphasizes that the turbulent inflow data generated by the present RFG technique for the LES computation can be a viable approach in accurately predicting the effects of inflow turbulence on the near wake turbulent flow characteristics around a bluff body.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2003.04a
• /
• pp.184-187
• /
• 2003

Complete prediction of second order permeability tensor for three dimensional circular braided preform is critical to understand the resin transfer molding process of composites. The permeability can be predicted by considering resin flow through the multi-axial fiber structure. In this study, permeability tensor for a 3-D circular braided preform is calculated by solving a boundary problem of a periodic unit cell. Flow field through the unit cell is obtained by using a 3-D finite volume method (FVM) and Darcy's law is utilized to obtain permeability tensor. Flow analysis for two cases that a fiber tow is regarded as impermeable solid and permeable porous medium is carried out respectively. It is found that the flow within the intra-tow region of the braided preform is negligible if inter-tow porosity is relatively high but the flow through the tow must be considered when the porosity is low. To avoid checkerboard pressure field and improve the efficiency of numerical computation, a new interpolation function for velocity variation is proposed on the basis of analytic solutions. Permeability of the braided preform is measured through a radial flow experiment and compared with the permeability predicted numerically.

• Journal of Environmental Science International
• /
• v.3 no.4
• /
• pp.367-372
• /
• 1994

A predictive model is demonstrated for gas removal rates from the aklosphere by dw deposition. Typical deposition velocities are complex functions of surface types, atmospheric stabilities, friction velocities, air pollutants, and so on. In this paper we simulated the calculation of dry deposition velocities near the earth surfaces, simultaneously we estimated real dry deposition velocities using the previous simulation. The measurement taken over a deciduous forest by Padro et d.(1988) were used to verify this model. In the comparison of the value of deposition velocity between numerical computation and observation, there are partially overestimations and underestimations between them, but we can speak that they are in a good accordance.