• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2000.10a
• /
• pp.136-140
• /
• 2000

It is desirable to predict springback quantitatively and accurately for the tool and process design in sheet stamping operations, however, it is blown very difficult. The result of springback analysis by the finite element method is sensitively influenced by numerical factors such as blank element size, number of integration point, punch velocity, contact algorithm etc. In the present work, a parametric study by Taguchi method is performed in order to evaluate the influence of numerical factors on springback Quantitatively and to obtain the combination of numerical factors which yields the best approximation to experimental data. Since springback is determined by the residual stress after forming process, it is important to evaluate stress distribution accurately. The oscillation in the time history curve of stress obtained by explicit FEM says that the stress solution at termination time is in very unstable state. Therefore, a variability study is also carried out in this study in order to assess the stability of implicit springback analysis starting from the stress solution by explicit forming simulation. The 2D draw bending process, one of the NUMISHEET '93 benchmark problems, is adopted as an application model.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.24 no.2
• /
• pp.47-54
• /
• 1987

Linear and nonlinear hydrodynamic force, which acts on submerged circular and eilliptic cylinders in oscillations as well as in advancing motion, are investigated as an initial-boundary value problem using a numerical method, which makes use of the source distribution on the body surface and the spectral method for treating the free surface waves. In the numerical code developed here, the boundary condition at the body surface is linearized. Using the numerical code so attained, nonlinear effects for different forward speeds and of the large-amplitude motion are computed. One of the major findings is that, when the forward speed is large, the added mass has its minimum and the damping force change rapidly around the frequency corresponding to the speed-frequency parameter, $\tau$=0.25, Compared to the result of Grue's [10], who used linear theory to get abrupt changes in values of the added mass and the damping force at the frequency corresponding to $\tau$=0.25, the present study, which takes nonlinear effects into account, shows much smoother variations near the frequency.

• Wind and Structures
• /
• v.35 no.1
• /
• pp.35-54
• /
• 2022

A comparative study of aeroelastic vibrations of spring-mass cylinder and chimneys, with the help of a few wake oscillator models available in the literature, is presented. The models include those proposed by Facchinetti, Farshidian and Dolatabadi method-I, Farshidian and Dolatabadi method-II, de Langre, Skop and Griffin. Besides, the linear model proposed by Simiu and Scanlan is also incorporated in the study. For chimneys, the first mode oscillation is considered, and the top displacements of the chimneys are evaluated using the considered models. The results of the analytical model are compared with those obtained from the numerical solution of the wake-oscillator coupled equations. The response behavior of the cylinder and three chimneys of different heights are studied and compared with respect to critical parametric variations. The results of the study indicate that the numerical analysis is essential to capture the effect of non-linear aeroelastic phenomena in the solutions, especially for small damping. Further, except for the models proposed by Farshidian and Dolatabadi, other models predict nearly the same responses. The non-linear model predicts a much higher response as compared to the linear model.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.2B
• /
• pp.123-130
• /
• 2012

MLP (Multi dimensional Limiting Process) is implemented to simulate shallow water flows, and its performance over conventional TVD limiters in multidimensional flows is verified through several numerical simulations. MLP was developed to control oscillations for multi-dimensional compressible flows and proved to improve accuracy, efficiency and robustness in compressible flows. In this study, we applies MLP to modeling shallow water equations(SWEs) given that the SWEs are amenable to be solved using the large range of numerical methods developed to deal with compressible flows and MLP has been yet used for SWEs. Simulation results through the benchmark tests show that MLP has favorable features such as numerical oscillation control and convergence behaviors comparable to the conventional limiters. Both numerical accuracy and stability are improved in multi-dimensional discontinuous flows.

• The Journal of Engineering Geology
• /
• v.19 no.2
• /
• pp.205-215
• /
• 2009

The conventional numerical models to analyze flow in subsurface porous media under the transient state usually generate numerical oscillation and unstability due to local flux domain for critical cases such as infiltration into initially dry soil during rainfall period. In this case, it is required refined mesh and small time step, but it decrease efficiency of computation. In this study, numerical unstability in discontinuity domain is removed by applying particle tracking algorithm to simulate unsteady subsurface flow with inflow boundary condition. Finally the hybrid LE FEM improving numerical stability is proposed. The hypothetical domains with unsteady uniform and nonuniform flow field were used to demonstrated algorithm verification. In comparison with analytic solution, we obtained reasonable results and conducted simulation of hypothetical 3-D recharge/pumping area. The proposed algorithm can simulate saturated/unsaturated porous media with more practical problems and will greatly contribute to accuracy and stability of numerical computation.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.32-33
• /
• 2003

Mechanism of a periodic oscillation of shock-induced combustion over a two- dimensional wedges and axi-symmetric cones were investigated through a series of numerical simulations at off-attaching condition of oblique detonation waves(ODW). A same computational domain over 40 degree half-angle was considered for two-dimensional and axi-symmetric shock-induced combustion phenomena. For two-dimensional shock-induced combustion, a 2H2+02+17N2 mixture was considered at Mach number was 5.85with initial temperature 292 K and initial pressureof 12 KPa. The Rankine-Hugoniot relation has solution of attached waves at this condition. For axi-symmetric shock-induced combustion, a H2+2O2+2Ar mixture was considered at Mach number was 5.0 with initial temperature 288 K and initial pressure of 200 mmHg. The flow conditions were based on the conditions of similar experiments and numerical studies.[1, 3]Numerical simulation was carried out with a compressible fluid dynamics code with a detailed hydrogen-oxygen combustion mechanism.[4, 5] A series of calculations were carried out by changing the fluid dynamic time scale. The length wedge is varied as a simplest way of changing the fluid dynamic time scale. Result reveals that there is a chemical kinetic limit of the detached overdriven detonation wave, in addition to the theoretical limit predicted by Rankine-Hugoniot theory with equilibrium chemistry. At the off-attaching condition of ODW the shock and reaction waves still attach at a wedge as a periodically oscillating oblique shock-induced combustion, if the Rankine-Hugoniot limit of detachment isbut the chemical kinetic limit is not.Mechanism of the periodic oscillation is considered as interactions between shock and reaction waves coupled with chemical kinetic effects. There were various regimes of the periodicmotion depending on the fluid dynamic time scales. The difference between the two-dimensional and axi-symmetric simulations were distinct because the flow path is parallel and uniform behind the oblique shock waves, but is not behind the conical shock waves. The shock-induced combustion behind the conical shockwaves showed much more violent and irregular characteristics.From the investigation of characteristic chemical time, condition of the periodic instability is identified as follows; at the detaching condition of Rankine-Hugoniot theory, (1) flow residence time is smaller than the chemical characteristic time, behind the detached shock wave with heat addition, (2) flow residence time should be greater than the chemical characteristic time, behind an oblique shock wave without heat addition.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.6
• /
• pp.619-627
• /
• 2004

In most industrial applications, the geometrical complexity is combined with the moving boundaries. These problems considerably increase the computational difficulties since they require, respectively, regeneration and deformation of the grid. As a result, engineering flow simulation is restricted. In order to solve this kind of problems the immersed boundary method was developed. In this study, the immersed boundary method is applied to the numerical simulation of stationary, rotating and oscillating cylinders in the 2-dimensional square cavity. No-slip velocity boundary conditions are given by imposing feedback forcing term to the momentum equation. Besides, this technique is used with a second-order accurate interpolation scheme in order to improve the accuracy of flow near the immersed boundaries. The governing equations for the mass and momentum using the immersed boundary method are discretized on the non-staggered grid by using the finite volume method. The results agree well with previous numerical and experimental results. This study presents the possibility of the immersed boundary method to apply to the complex flow experienced in the industrial applications. The usefulness of this method will be confirmed when we solve the complex geometries and moving bodies.

• Journal of computational fluids engineering
• /
• v.21 no.4
• /
• pp.90-95
• /
• 2016

Leidenfrost drops with ellipsoidal shaping can control the bouncing height by adjusting the aspect ratio(AR) of the shape at the moment of impact. In this work, we focus on the effect of the AR and the impact Weber number(We) on the non-axisymmetrical spreading dynamics of the drop, which plays an important role in the control of bouncing. To understand the impact dynamics, the numerical simulation is conducted for the ellipsoidal drop impact upon the perfect non-wetting solid surface by using volume of fluid method, which shows the characteristics of the spreading behavior in each principal axis. As the AR increases, the drop has a high degree of the alignment into one principal axis, which leads to the consequent suppression of bouncing height with shape oscillation. As the We increases, the maximum spreading diameters in the principal axes both increase whereas the contact time on the solid surface rarely depends on the impact velocity at the same AR. The comprehensive understanding of the ellipsoidal drop impact upon non-wetting surface will provide the way to control of drop deposition in applications, such as surface cleaning and spray cooling.

• Structural Engineering and Mechanics
• /
• v.50 no.4
• /
• pp.487-500
• /
• 2014

An analytical solution for the nonlinear in-plane free oscillations of the suspended cable which contains the quadratic and cubic nonlinearities is investigated via the homotopy analysis method (HAM). Different from the existing analytical technique, the HAM is indeed independent of the small parameter assumption in the nonlinear vibration equation. The nonlinear equation is established by using the extended Hamilton's principle, which takes into account the effects of the geometric nonlinearity and quasi-static stretching. A non-zero equilibrium position term is introduced due to the quadratic nonlinearity in order to guarantee the rule of the solution expression. Therefore, the mth-order analytic solutions of the corresponding equation are explicitly obtained via the HAM. Numerical results show that the approximate solutions obtained by using the HAM are in good agreement with the numerical integrations (i.e., Runge-Kutta method). Moreover, the HAM provides a simple way to adjust and control the convergent regions of the series solutions by means of an auxiliary parameter. Finally, the effects of initial conditions on the linear and nonlinear frequency ratio are investigated.

• Geomechanics and Engineering
• /
• v.13 no.2
• /
• pp.269-284
• /
• 2017

A new technique is proposed to obtain more effective screening efficiency against the ground vibration using intermittent geofoam (IF) in-filled trench. The numerical analysis is performed by employing two-dimensional finite element method under dynamic condition. Vertically oscillated strip foundation is considered as the vibration source. In presence of the ground vibration, the vertical displacements at different locations (pick-up points) along the ground surface are captured to determine the amplitude reduction factor (ARF), which helps to assess the efficiency of the vibration screening technique. The efficiency of IF over continuous geofoam (CF) in-filled vibration barriers is assessed by varying the geofoam density, the location of trench and the frequency of excitation. The results from this study indicate that a significant reduction in ARF can be achieved by using intermittent geofoam as compared to continuous geofoam. Further, it is noticed that the efficiency of IF increases with an increase in the frequency of the vibrating source. These encouraging results put forward the potential of utilising intermittent geofoam as a vibration screening material.