• Geomechanics and Engineering
• /
• v.18 no.3
• /
• pp.225-234
• /
• 2019

A numerical stepwise approach for cavity expansion problem in strain-softening rock or soil mass is investigated, which is compatible with Mohr-Coulomb and generalized Hoek-Brown failure criteria. Based on finite difference method, plastic region is divided into a finite number of concentric rings whose thicknesses are determined internally to satisfy the equilibrium and compatibility equations, the material parameters of the rock or soil mass are assumed to be the same in each ring. For the strain-softening behavior, the strength parameters are assumed to be a linear function of deviatoric plastic strain (${\gamma}p^*$) for each ring. Increments of stress and strain for each ring are calculated with the finite difference method. Assumptions of large-strain for soil mass and small-strain for rock mass are adopted, respectively. A new numerical stepwise approach for limited pressure and plastic radius are obtained. Comparisons are conducted to validate the correctness of the proposed approach with Vesic's solution (1972). The results show that the perfectly elasto-plastic model may underestimate the displacement and stresses in cavity expansion than strain-softening coefficient considered. The results of limit expansion pressure based on the generalised H-B failure criterion are less than those obtained based on the M-C failure criterion.

• Advances in Computational Design
• /
• v.9 no.1
• /
• pp.39-52
• /
• 2024

If the governing differential equation arising from engineering problems is treated as an analytic, continuous and derivable function, it can be expanded by one point as a series of finite numbers. For the function to be zero for each value of its domain, the coefficients of each term of the same power must be zero. This results in a recursive relationship which, after applying the natural conditions or the boundary conditions, makes it possible to obtain the values of the derivatives of the function with acceptable accuracy. The elastoplastic analysis of an inhomogeneous thick sphere of metallic materials with linear variation of the modulus of elasticity, yield stress and Poisson's ratio as a function of radius subjected to internal pressure is presented. The Beltrami-Michell equation is established by combining equilibrium, compatibility and constitutive equations. Assuming axisymmetric conditions, the spherical coordinate parameters can be used as principal stress axes. Since there is no analytical solution, the natural boundary conditions are applied and the governing equations are solved using a proposed new method. The maximum effective stress of the von Mises yield criterion occurs at the inner surface; therefore, the negative sign of the linear yield stress gradation parameter should be considered to calculate the optimal yield pressure. The numerical examples are performed and the plots of the numerical results are presented. The validation of the numerical results is observed by modeling the elastoplastic heterogeneous thick sphere as a pressurized multilayer composite reservoir in Abaqus software. The subroutine USDFLD was additionally written to model the continuous gradation of the material.

• Journal of computational fluids engineering
• /
• v.13 no.2
• /
• pp.68-72
• /
• 2008

NUFLEX is a general purpose thermo/fluid flow analysis program which has various physical models including spray. In NUFLEX, spray models are composed of breakup and collision models of droplet. However, in case of diesel engine, interaction between wall-film and impingement model considering heat transfer is not coded in NUFLEX. In this study, Lee & Ryou impingement & wall-film model considering heat transfer is applied to NUFLEX. For the verification of this NUFLEX program, numerical results are compared with experimental data. Differences of film thickness and radius between numerical results and experimental data are within 10% error range. The results show that NUFLEX can be used for comprehensive analysis of spray phenomena.

• 한국연소학회:학술대회논문집
• /
• 2014.11a
• /
• pp.175-178
• /
• 2014

The flame cell dynamics in 2-D opposed nonpremixed tubular configuration was investigated using high-fidelity numerical simulations. The diffusive-thermal instability occurs as the $Damk{\ddot{o}}hler$ number, Da, approaches the 1-D extinction limit of the tubular flames and several flame cells are generated depending on Da, and flame radius. In general, the number of flame cells are found close to the largest wave number from the linear stability analysis. It was also found from the displacement speed analysis that during the local flame extinction and cell formation, negative edge flame speed is observed due to small gain from reaction compared to large loss from diffusion.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.56 no.9
• /
• pp.1571-1576
• /
• 2007

This report describes the analytical characteristic of sinewave air gap flux density for the brushless AC motor with surface permanent magnet. The analysis of sinewave air gap flux density is the key to expect the performance of back EMF for the design of brushless AC motor. The numerical analysis and FEM analysis are performed to adopt radial and parallel flux magnetization of magnet on the rotor. And it is also executed to vary the magnet arc angle and arc radius for the condition of constant and non constant air gap. This report is focused on the characteristic of sinewave air gap flux density for permanent magnet of surface brushless AC motor. This results also have more reliable data against the previous paper which had the representative numerical analysis of air gap flux density[1][2].

• Journal of computational fluids engineering
• /
• v.19 no.4
• /
• pp.45-51
• /
• 2014

We develop a numerical method for solving the radiation hydrodynamic equations in one-dimensional spherical coordinates. The present method is validated through simulations of shock tube, thermal radiative diffusion and point explosion problems. The transient growth of the fireball is investigated by varying explosion yields. The present study clearly captures well-known breakaway phenomena related to the shock separation between pressure waves and thermal shock front. The fireball radius at the breakaway point is roughly increased by the yield to power of 0.4.

• Journal of computational fluids engineering
• /
• v.8 no.4
• /
• pp.6-15
• /
• 2003

It is an aim of this study to perform extensive numerical study for analyzing the anisotropic turbulence effects on spatial and temporal behaviors of droplet for impinging sprays. The turbulence model of Durbin is used for comparisons with the k-ε model. The turbulence-induced dispersions of droplets are considered to describe the anisotropy of turbulence effectively and spray/wall interactions are simulated using the model of Lee and Ryou. Present study investigates the overall and the internal structures of impinging diesel sprays such as spray shapes, radius and height of wall sprays, Sauter mean diameter (SMD), local droplet velocity, and local gas velocity and compared the results with experimental data by two adopted turbulence models. When the anisotropy effect of turbulence is included, better predictions for both gas and droplet tangential velocities are obtained, compared to the k-ε model. It is concluded that anisotropic effect of turbulence should be considered for simulating impinging diesel sprays.

• Journal of Mechanical Science and Technology
• /
• v.18 no.4
• /
• pp.718-730
• /
• 2004

The solution for the natural convection in internally finned horizontal annuli is obtained by using a numerical simulation of time-dependent and two-dimensional governing equations. The fins existing in annuli influence the flow pattern, temperature distribution and heat transfer rate. The variations of the On configuration suppress or accelerate the free convective effects compared to those of the smooth tubes. The effects of fin configuration, number of fins and ratio of annulus gap width to the inner cylinder radius on the fluid flow and heat transfer in annuli are demonstrated by the distribution of the velocity vector, isotherms and streamlines. The governing equations are solved efficiently by using a parallel implementation. The technique is adopted for reduction of the computation cost. The parallelization is performed with the domain decomposition technique and message passing between sub-domains on the basis of the MPI library. The results from parallel computation reveal in consistency with those of the sequential program. Moreover, the speed-up ratio shows linearity with the number of processor.

• Tribology and Lubricants
• /
• v.14 no.4
• /
• pp.58-63
• /
• 1998

This paper presents a study of the elastohydrodynamic lubrication of elliptical contacts where lubricant entrainment coincides with the major axis of the Hertzian contact ellipse. A finite difference method and the Newton-Raphson method are applied to analyze the problem. Film contours and pressure distributions are compared with the results for lubricant entrainment coincides with the minor axis. Variations of the minimum and central film thicknesses with the radius ratio are also examined. The results showed that the present numerical scheme can be used generally in the analysis of the EHL of elliptical contacts where the lubricant entraining vector did not coincide with either of the principal axis of the conjunction.

• Journal of the Korean Society of Visualization
• /
• v.14 no.1
• /
• pp.46-50
• /
• 2016

Direct numerical simulations (DNSs) of turbulent pipe flows with temporal deceleration were performed to examine response of the turbulent flows to the deceleration. The simulations were started with a fully-developed turbulent pipe flow at the Reynolds number, $Re_D=24380$, based on the pipe radius and the laminar centerline velocity, and three different constant temporal decelerations were applied to the initial flow with varying dU/dt = -0.001274, -0.00625 and -0.025. It was shown that the mean flows were greatly affected by temporal decelerations with downward shift of log law, and turbulent intensities were increased in particular in the outer layer, compared to steady flows at a similar Reynolds number. The analysis of Reynolds shear stress showed that second- and fourth-quadrant Reynolds shear stresses were increased with the decelerations, and the increase of the turbulence was attributed to enhancement of outer turbulent vortical structures by the temporal decelerations.