• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.18 no.11
• /
• pp.690-696
• /
• 2017

Various types of alternative energy sources to petroleum are being developed both domestically and internationally as clean energy that does not emit greenhouse gases. In particular, offshore wind power has been studied because the wind resources are relatively limitless and the wind power is relatively smaller than onshore. In this study, to analyze the scour phenomenon around offshore wind foundations, mono pile and tripod pile foundations were simulated using a FLOW-3D model. The scour phenomenon was evaluated for mono piles: one is a pile with a 5 m diameter and d=1.69 m and the other is a pile with a 5 m diameter. Numerical analysis showed that in the latter, the falling-flow increased and the maximum scour depth occurred more than 1.7 times. For a tripod pile foundation, the measured velocity and the maximum wave condition were applied to the upstream boundary condition, respectively, and the scour phenomenon was evaluated. When the maximum wave condition was applied, the maximum scour depth occurred more than about 1.3 times. When the LES model was applied, the scour depth reached equilibrium, whereas the numerical results of the RNG model show that the scour phenomenon occurred in the entire boundary area and the scour depth did not reach equilibrium. To evaluate the scour phenomenon around offshore wind foundations, it is reasonable to apply the wave condition and the LES turbulence model to numerical model applications.

• Geomechanics and Engineering
• /
• v.14 no.2
• /
• pp.141-149
• /
• 2018

In this paper, the creep behavior of bolted rock was analyzed by using the unconfined creep tests and the numerical results. Based on the test results, the Bolted Burgers creep model (B-B model) was proposed to clarify the creep mechanism of rock mass due to rock bolts. As to the simulation of the creep behaviour of bolted rock, a new user-defined incremental iterative format of the B-B model was established and the open-source $FLAC^{3D}$ code was written by using the object-oriented language (C++). To check the reliability of the present B-B creep constitutive model program, a numerical model of a tunnel with buried depth of 1000 m was established to analyze the creep response of the tunnel with the B-B model support, the non-support and the bolt element support. The simulation results show that the present B-B model is consistent with the calculated results of the inherent bolt element in $FLAC^{3D}$, and the convergence deformation can be more effectively controlled when the proposed B-B model is used in the $FLAC^{3D}$ software. The big advantage of the present B-B creep model secondarily developed in the $FLAC^{3D}$ software is the high computational efficiency.

• Computers and Concrete
• /
• v.4 no.3
• /
• pp.221-241
• /
• 2007

Based on the experiment results, the damage and fracture behavior of concrete at the ages of 1d, 2d, 7d and 28d, in three-point bending and uniaxial tensile tests, were simulated with a finite element program, ABAQUS. The critical stress intensity factor $K_{IC}^s$ and the critical crack tip opening displacement ($CTOD_C$) of concrete were calculated with effective-elastic crack approach for the three-point bending test of grade C30 concrete. Based on the crack band model, a bilinear strain-softening curve was derived to simulate the LOAD-CMOD curves and LOAD-Displacement curves. In numerical analysis of the uniaxial tension test of concrete of grade C40, the damage and fracture mechanics were combined. The smeared cracking model coupling with damaged variable was adopted to evaluate the onset and development of microcracking of uniaxial tensile specimen. The uniaxial tension test was simulated by invoking the damage plastic model which took both damage and plasticity as inner variables with user subroutines. All the numerical simulated results show good agreement with the experimental results.

• Structural Engineering and Mechanics
• /
• v.40 no.2
• /
• pp.257-277
• /
• 2011

In this paper the geometrically nonlinear continuum plate finite element model, hitherto not reported in the literature, is developed using the total Lagrange formulation. With the layerwise displacement field of Reddy, nonlinear Green-Lagrange small strain large displacements relations (in the von Karman sense) and linear elastic orthotropic material properties for each lamina, the 3D elasticity equations are reduced to 2D problem and the nonlinear equilibrium integral form is obtained. By performing the linearization on nonlinear integral form and then the discretization on linearized integral form, tangent stiffness matrix is obtained with less manipulation and in more consistent form, compared to the one obtained using laminated element approach. Symmetric tangent stiffness matrixes, together with internal force vector are then utilized in Newton Raphson's method for the numerical solution of nonlinear incremental finite element equilibrium equations. Despite of its complex layer dependent numerical nature, the present model has no shear locking problems, compared to ESL (Equivalent Single Layer) models, or aspect ratio problems, as the 3D finite element may have when analyzing thin plate behavior. The originally coded MATLAB computer program for the finite element solution is used to verify the accuracy of the numerical model, by calculating nonlinear response of plates with different mechanical properties, which are isotropic, orthotropic and anisotropic (cross ply and angle ply), different plate thickness, different boundary conditions and different load direction (unloading/loading). The obtained results are compared with available results from the literature and the linear solutions from the author's previous papers.

• Journal of computational fluids engineering
• /
• v.20 no.2
• /
• pp.46-51
• /
• 2015

In this paper, we have developed numerical model for particulated flow through narrow slit using Eulerian-Granular method. Commercial software (FLUENT) was utilized as simulation tool and main focus was to identify the effect from various numerical options for modeling of solid particles as continuos phase in granular flow. Gidaspow model was chosen as basic model for solid viscosity and drag model. And lun-et-al model was used as solid pressure and radial distribution model, respectively. Several other model options in FLUENT were tested considering the cross related effect. Mass flow rate of the particulate through the slit was measured to compare. Due to the high volume density of the stacked particulates above the slit, effect from various numerical options were not significant. The numerical results from basic model were also compared with experimental results and showed very good agreement. The effects from the characteristics of particles such as diameter, angle of internal friction, and collision coefficient were also analyzed for future design of velocity resistance layer in solar thermal absorber. Angle of internal friction was found to be the dominat variable for the particle mass flow rate considerably. More defined 3D model along with energy equation for complete solar thermal absorber design is currently underway.

• Wind and Structures
• /
• v.17 no.6
• /
• pp.647-670
• /
• 2013

The paper presents a numerical approach to study of fluid flow characteristics and to predict performance of wind turbines. The numerical model is based on Finite-volume method (FVM) discretization of unsteady Reynolds-averaged Navier-Stokes (URANS) equations. The movement of turbine blades is modeled using moving mesh technique. The turbulence is modeled using commonly used turbulence models: Renormalization Group (RNG) k-${\varepsilon}$ turbulence model and the standard k-${\varepsilon}$ and k-${\omega}$ turbulence models. The model is validated with the experimental data over a large range of tip-speed to wind ratio (TSR) and blade pitch angles. In order to demonstrate the use of numerical method as a tool for designing wind turbines, two dimensional (2-D) and three-dimensional (3-D) simulations are carried out to study the flow through a small scale Darrieus type H-rotor Vertical Axis Wind Turbine (VAWT). The flows predictions are used to determine the performance of the turbine. The turbine consists of 3-symmetrical NACA0022 blades. A number of simulations are performed for a range of approaching angles and wind speeds. This numerical study highlights the concerns with the self-starting capabilities of the present VAWT turbine. However results also indicate that self-starting capabilities of the turbine can be increased when the mounted angle of attack of the blades is increased. The 2-D simulations using the presented model can successfully be used at preliminary stage of turbine design to compare performance of the turbine for different design and operating parameters, whereas 3-D studies are preferred for the final design.

• Structural Engineering and Mechanics
• /
• v.26 no.3
• /
• pp.343-361
• /
• 2007

This article presents a numerical tool for dimensioning two-threaded fasteners connecting prismatic parts subjected to fatigue tension loads that are coplanar with the screw axis. A simplified numerical model is developed from unidirectional finite elements, modeling the connected parts and screws with bent elements and the elastic contact layer between the parts with springs. An algorithm updating the contact stiffness matrix, calculating forces and displacements at each node of the structure and thus normal stresses in the screws in both static and fatigue is further developed using C language. An experimental study is also conducted in parallel with the numerical approach to validate the developed model assumptions, the numerical model and the 3D finite element results. Since stiffness values for the compressive zones in the parts are analytically difficult to determine, a statistical software method is used, from which a tuning factor is derived for identifying these stiffness values. The method is also applied to set out the influence of each parameter on the fatigue behaviour of each screw. Finally, the developed model will be used to establish a new, sophisticated, fast and accurate tool for dimensioning bolted mechanical structures.

• Structural Engineering and Mechanics
• /
• v.61 no.3
• /
• pp.325-334
• /
• 2017

In the present study, a numerical model for probability analysis of optimal design of fatigue non-uniform crack growth behaviour of a cracked aluminium 2024 T3 plate repaired with a bonded composite patch is investigated. The proposed 3D numerical model has advanced in literatures, which gathers in a unique study: problems of reliability, optimization, fatigue, cracks and repair of plates subjected to tensile loadings. To achieve this aim, a finite element modelling is carried out to determine the evolution of the stress intensity factor at the crack tip Paris law is used to predict the fatigue life for a give n crack. To have an optimal volume of our patch satisfied the practical fatigue life, a procedure of optimization is proposed. Finally, the probabilistic analysis is performed in order to a show that optimized patch design is influenced by uncertainties related to mechanical and geometrical properties during the manufacturing process.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.181-184
• /
• 2006

High-speed flight vehicle have various cavities. The supersonic cavity flow is complicated due to vortices, flow separation and reattachment, shock and expansion waves. The general cavity flow phenomena include the formation and dissipation of vortices, which induce oscillation and noise. The oscillation and noise greatly affect flow control, chemical reaction, and heat transfer processes. The supersonic cavity' flow with high Reynolds number is characterized by the pressure oscillation due to turbulent shear layer, cavity geometry, and resonance phenomenon based on external flow conditions, The resonance phenomena can damage the structures around the cavity and negatively affect aerodynamic performance and stability. In the present study, we performed numerical analysis of cavities by applying the unsteady, compressible three dimensional Reynolds-Averaged Navier-Stokes(RANS) equations with the ${\kappa}-{\omega}$ turbulence model. The cavity model used for numerical calculation had a depth(D) of 15mm cavity aspect ratio(L/D) of 3, width to spanwise ratio(W/D) of 1.0 to 5.0. Based on the PSD(Power Spectral Density) and CSD(Cross Spectral Density) analysis of the pressure variation, the dominant frequency was analyized and compared with the results of Rossiter's Eq.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.62-66
• /
• 2006

High-speed flight vehicle have various cavities. The supersonic cavity flow is complicated due to vortices, flow separation, reattachment, shock waves and expansion waves. The general cavity flow phenomena includes the formation and dissipation of vortices, which induce oscillation and noise. The oscillation and noise greatly affect flow control, chemical reaction, and heat transfer processes. The supersonic cavity flow with high Reynolds number is characterized by the pressure oscillation due to turbulent shear layer, cavity geometry, and resonance phenomenon based on external flow conditions. The resonance phenomena can damage the structures around the cavity and negatively affect aerodynamic performance and stability. In the present study, we performed numerical analysis of cavities by applying the unsteady, compressible three dimensional Reynolds-Averaged Navier-Stokes(RANS) equations with the ${\kappa}-{\omega}$ turbulence model. The cavity model used for numerical calculation had a depth(D) of 15mm cavity aspect ratio (L/D) of 3, width to spanwise ratio(W/D) of 1.0 to 5.0. Based on the PSD(Power Spectral Density) and CSD(Cross Spectral Density) analysis of the pressure variation, the dominant frequency was analyzed and compared with the results of Rossiter's Eq.