• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.84-90
• /
• 2001

Secondary instability of flow past a circular cylinder is examined using direct numerical simulation at Reynolds number 220 and 250. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme, and the pressure and viscosity terms are applied by the 4th order central difference. In spanwise, Navier-Stokes equation is distributed using Spectral Method. The critical Reynolds number for this instability is found to be about Re=190. The secondary instability leads re three-dimensionality with a spanwise wavelength about 4 cylinder diameters at onset (A-mode). Results of three-dimensional effect in wake of a circular cylinder are represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• Smart Structures and Systems
• /
• v.18 no.2
• /
• pp.355-374
• /
• 2016

This paper presents and compares a one-dimensional (1D) bending theory for piezoelectric thin beam-type structures with resistive-inductive electrodes to ANSYS$^{(R)}$ three-dimensional (3D) finite element (FE) analysis. In particular, the lateral deflections and vibrations of slender piezoelectric beams are considered. The peculiarity of the piezoelectric beam model is the modeling of electrodes in such a manner that is does not fulfill the equipotential area condition. The case of ideal, perfectly conductive electrodes is a special case of our 1D model. Two-coupled partial differential equations are obtained for the lateral deflection and for the voltage distribution along the electrodes: the first one is an extended Bernoulli-Euler beam equation (second-order in time, forth order in space) and the second one the so-called Telegrapher's equation (second-order in time and space). Analytical results of our theory are validated by 3D electromechanically coupled FE simulations with ANSYS$^{(R)}$. A clamped-hinged beam is considered with various types of electrodes for the piezoelectric layers, which can be either resistive and/or inductive. A natural frequency analysis as well as quasi-static and dynamic simulations are performed. A good agreement between the extended beam theory and the FE results is found. Finally, the practical relevance of this type of electrodes is shown. It is found that the damping capability of properly tuned resistive or resistive-inductive electrodes exceeds the damping performance of beams, where the electrodes are simply linked to an optimized impedance.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.3 no.1
• /
• pp.1-8
• /
• 1999

In the present paper an attempt has been made to simulate the secondary injection-primary flow interaction in the conical rocket nozzle and to derive the performance of secondary injection thrust vector control(SITVC) system. Complex three-dimensional flowfield induced by the secondary injection is numerically analyzed by solving unsteady three-dimensional Euler equation with Beam and Warming's implicit approximate factorization method. Emphasized in the present study is the effect of secondary injection such as secondary mass flow rates and the momentum of secondary/primary nozzle flow mass rates upon the gross system performance parameters such as thrust ratio, specific impulse ratio and deflection angle. The results obtained in terms of system performance parameters show that lower secondary mass flow rate is advantageous for to reduce secondary specific impulse loss. It is further found that the nozzle with secondary jet injected downstream and interacting with fast primary flow is preferable for efficient and stable SITVC over the wide range of use with the penalty of side specific impulse loss.

• Journal of the Society of Naval Architects of Korea
• /
• v.40 no.4
• /
• pp.16-21
• /
• 2003

The very large container ships have been built recently and those ships have very small structural rigidity compared with the other conventional ships. As a result, the destruction of ship hull is occurred by the springing including to warping phenomena due to encounter waves. In this study, the solutions of hydrodynamic coefficients are obtained by solving the three dimensional source distribution method and the forward speed Green function representing a translating and pulsating source potential for infinite water depth is used to calculating the integral equation. The vessel is longitudinally divided into various sections and the added mass, wave damping and wave exciting forces of each section is calculated by integrating the dynamic pressures over the mean wetted section surface. The equations for six degree freedom of motions is obtained for each section in the frequency domain and stiffness matrix is calculated by Euler beam theory. The computations are carried out for very large ship and effects of bending and torsional ridigity on the wave frequency and angle are investigated.

• Steel and Composite Structures
• /
• v.42 no.1
• /
• pp.75-90
• /
• 2022

The pointwise equilibrium polynomial (PEP) element considering local second-order effect has been widely used in direct analysis of many practical engineering structures. However, it was derived according to Euler-Bernoulli beam theory and therefore it cannot consider shear deformation, which may lead to inaccurate prediction for deep beams. In this paper, a novel beam-column element based on Timoshenko beam theory is proposed to overcome the drawback of PEP element. A fifth-order polynomial is adopted for the lateral deflection of the proposed element, while a quadric shear strain field based on equilibrium equation is assumed for transverse shear deformation. Further, an additional quadric function is adopted in this new element to account for member initial geometrical imperfection. In conjunction with a reliable and effective three-dimensional (3D) co-rotational technique, the proposed element can consider both member initial imperfection and transverse shear deformation for second-order direct analysis of frame structures. Some benchmark problems are provided to demonstrate the accuracy and high performance of the proposed element. The significant adverse influence on structural behaviors due to shear deformation and initial imperfection is also discussed.

• Coupled systems mechanics
• /
• v.13 no.3
• /
• pp.247-275
• /
• 2024

The paper studies the influence of the fluid flow velocity and flow direction in the initial state on the dispersion of the axisymmetric waves propagating in the inhomogeneously pre-stressed hollow cylinder containing this fluid. The corresponding eigenvalue problem is formulated within the scope of the three-dimensional linearized theory of elastic waves in bodies with initial stresses, and with linearized Euler equations for the inviscid compressible fluid. The discrete-analytical solution method is employed, and analytical expressions of the sought values are derived from the solution to the corresponding field equations by employing the discrete-analytical method. The dispersion equation is obtained using these expressions and boundary and related compatibility conditions. Numerical results related to the action of the fluid flow velocity and flow direction on the influence of the inhomogeneous initial stresses on the dispersion curves in the zeroth and first modes are presented and discussed. As a result of the analyses of the numerical results, it is established how the fluid flow velocity and flow direction act on the magnitude of the influence of the initial inhomogeneous stresses on the wave propagation velocity in the cylinder containing the fluid.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.38 no.3
• /
• pp.50-57
• /
• 1996

A numerical method is presented to obtain the natural frequencies and mode shapes of the arbitrary tapered beams with static deflection due to arbitrary distributed dead loads. The differential equation governing free vibration of such beams is derived and solved numerically. The double integration method using the trapezoidal rule is used to solve the static behaviour of beams loaded arbitrary distributed dead load. Also, the Improved Euler method and the determinant search method are used to integrate the differential equation subjected to the boundary conditions and to determine the natural frequencies of the beams, respectively. In the numerical examples, the various geometries of the beams are considered : (1) linearly tapered beams as the arbitrary variable cross-section, (2) the triangular, sinusoidal and uniform loads as the arbitrary distributed dead loads and (3) the hinged-hinged, clamped-clamped and hinged-clamped ends as the end constraints. All numerical results are shown as the non-dimensional forms of the system parameters. The lowest three natural frequencies versus load parameter, slenderness ratio and section ratio are reported in figures. And for the comparison purpose, the typical mode shapes with and without the effects of static deflection are presented in the figure. According to the numerical results obtained in this analysis, the following conclusions may be drawn : (1) the natural frequencies increase when the effects of static deflections are included, (2) the effects are larger at the lower modes than the higher ones and (3) it should be betteF to include the effect of static deflection for calculating the frequencies when the beams are supported by both hinged ends or one hinged end.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.3
• /
• pp.215-223
• /
• 2009

In this paper, the helicopter aerodynamics is simulated in hovering and forward flight. Also, an overlapped grid technique is applied in this simulation to consider the blade motion and moving effects. The Caradonna & Tung's rotor blade was selected to analyze the unsteady aerodynamics in hovering and non-lift forward flight. Also, the AH-1G rotor blade was selected in forward flight. In forward flight case, the numerical trim was applied to determine the cyclic pitching angles using Newton-Raphson method, and the numerical results were in good agreement with experimental data, especially, the BVI effects were well simulated in advancing side in comparison other numerical results. The governing equation is a three dimensional unsteady Euler equation, and the Riemann invariants condition is used for inflow and outflow at the boundary.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.12 no.6
• /
• pp.111-118
• /
• 2004

The work presented in this paper concerns the aerodynamic characteristics and compression wave generated in a tunnel when a high speed train enters it. A large number of solutions have been proposed to reduce the amplitude of the pressure gradient in tunnels and some of the most efficient solutions consist of (a) addition ofa blind hood, (b) addition of inclined part at the entrance, and (c) holes in the ceiling of the tunnel. These are numerically studied by using the three-dimensional unsteady compressible Euler equation solver with ALE, CFD code, based on FEM method. Computational results showed that the smaller inclined angle leads to the lower pressure gradient of compression wave front. This study indicated that the most efficient slant angle is in the range from $30^{\circ}$ to $50^{\circ}$. The maximum pressure gradient is reduced by $26.81\%$ for the inclined angle of $30^{\circ}$ as compared to vertical entry. Results also showed that maximum pressure gradient can be reduced by $15.94\%$ in blind hood entry as compared to $30^{\circ}$ inclined tunnel entry. Furthermore, the present analysis showed that inclined slant angle has little effect on aerodynamic drag. Comparison of the pressure gradient between the inclined tunnel hood and the vertical entry with air vent holes indicated that the optimum inclined tunnel hood is much more effective way in reducing pressure gradient and increasing the pressure rise time.

• Coupled systems mechanics
• /
• v.12 no.1
• /
• pp.41-68
• /
• 2023

The present work is concerned with the study of the influence of inhomogeneous initial stresses in a hollow cylinder containing a compressible inviscid fluid on the propagation of axisymmetric longitudinal waves propagating in this cylinder. The study is carried out using the so-called three-dimensional linearized theory of elastic waves in bodies with initial stresses to describe the motion of the cylinder and using the linearized Euler equations to describe the flow of the compressible inviscid fluid. It is assumed that the inhomogeneous initial stresses in the cylinder are caused by the internal pressure of the fluid. To solve the corresponding eigenvalue problem, the discrete-analytic solution method is applied and the corresponding dispersion equation is obtained, which is solved numerically, after which the corresponding dispersion curves are constructed and analyzed. To obtain these dispersion curves, parameters characterizing the magnitude of the internal pressure, the ratio of the sound velocities in the cylinder material and in the fluid, and the ratio of the material densities of the fluid and the cylinder are introduced. Based on these parameters, the influence of the inhomogeneous initial stresses in the cylinder on the dispersion of the above-mentioned waves in the considered hydro-elastic system is investigated. Moreover, based on these results, appropriate conclusions about this influence are drawn. In particular, it is found that the character of the influence depends on the wavelength. Accordingly, the inhomogeneous initial stresses before (after) a certain value of the wavelength lead to a decrease (increase) of the wave propagation velocity in the zeroth and first modes.