• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.6
• /
• pp.1994-2004
• /
• 1996

Numerical calculation was applied to supersonic under-expanded jets, and compared with the results of a linear theory and other experiments. TVD difference scheme was employed to solve 2-dimensional and axisymmetric inviscid Euler equation. This paper aims to explore the effects of angle of divergence and design Mach number of nozzle on the structure of under-expanded jets. The angle of divergence was varied from 0 to 20 deg. The results show that the length of the first cell of the under-expanded jets decreases and Mach disk generates at lower nozzle pressure ratio, if the angle of divergence or design Mach number of nozzle increases. The distance from the nozzle exit to Mach disk in 2-dimensional jets becomes much larger than that of axisymmetric jets, and the widths of the jet boundary and the barrel shock wave are also larger than that of axisymmetric jets. Calculation results indicate that the configuration of the under-expanded jets is strongly dependent on the nozzle pressure ratio.

• Earthquakes and Structures
• /
• v.12 no.4
• /
• pp.457-468
• /
• 2017

When studying the vibration of a suspension bridge based on the traffic-bridge coupled system, most researchers ignored the contribution of the pavement response. For example, the pavement was simplified as a rigid base and the deformation of pavement was ignored. However, the action of deck pavement on the vibration of vehicles or bridges should not be neglected. This study is mainly focused on establishing a new methodology fully considering the effects of bridge deck pavement, probabilistic traffic flows, and varied road roughness conditions. The bridge deck pavement was modeled as a boundless Euler-Bernoulli beam supported on the Kelvin model; the typical traffic flows were simulated by the improved Cellular Automaton (CA) traffic flow model; and the traffic-pavement-bridge coupled equations were established by combining the equations of motion of the vehicles, pavement, and bridge using the displacement and interaction force relationship at the contact locations. The numerical studies show that the proposed method can more rationally simulate the effect of the pavement on the vibrations of bridge and vehicles.

• Steel and Composite Structures
• /
• v.42 no.3
• /
• pp.397-405
• /
• 2022

The existing researches on the dynamics of the fluid-conveying pipes only focus on stability and vibration problems, and there is no literature report on the wave propagation of the fluid-conveying pipes. Therefore, the purpose of this paper is to explore the propagation characteristics of longitudinal and flexural waves in the fluid-conveying pipes. First, it is assumed that the material properties of the fluid-conveying pipes vary based on a power function of the thickness. In addition, it is assumed that the material properties of both the fluid and the pipes are closely depended on temperature. Using the Euler-Bernoulli beam equation and based on the linear theory, the motion equations considering the thermal-mechanical-fluid coupling is derived. Then, the exact expressions of phase velocity and group velocity of longitudinal waves and bending waves in the fluid-conveying pipes are obtained by using the eigenvalue method. In addition, we also studied the free vibration frequency characteristics of the fluid-conveying pipes. In the numerical analysis, we successively studied the influence of temperature, functional gradient index and liquid velocity on the wave propagation and vibration problems. It is found that the temperature and functional gradient exponent decrease the phase and group velocities, on the contrary, the liquid flow velocity increases the phase and group velocities. However, for vibration problems, temperature, functional gradient exponent parameter, and fluid velocity all reduce the natural frequency.

• Advances in aircraft and spacecraft science
• /
• v.10 no.5
• /
• pp.419-437
• /
• 2023

This paper presents the study of the effects of rigid-body motion simultaneously with the presence of the effects of temporal variation due to the existence of morphing speed on the aeroelastic stability of the two-stage telescopic wings, and hence this is the main novelty of this study. To this aim, Euler-Bernoulli beam theory is used to model the bending-torsional dynamics of the wing. The aerodynamic loads on the wing in an incompressible flow regime are determined by using Peters' unsteady aerodynamic model. The governing aeroelastic equations are discretized employing a finite element method based on the beam-rod model. The effects of rigid-body motion on the length-based stability of the wing are determined by checking the eigenvalues of system. The obtained results are compared with those available in the literature, and a good agreement is observed. Furthermore, the effects of different parameters of rigid-body such as the mass, radius of gyration, fuselage center of gravity distance from wing elastic axis on the aeroelastic stability are discussed. It is found that some parameters can cause unpredictable changes in the critical length and frequency. Also, paying attention to the fuselage parameters and how they affect stability is very important and will play a significant role in the design.

• Wind and Structures
• /
• v.38 no.3
• /
• pp.161-170
• /
• 2024

In recent years, there have been an increasing number of experimental studies showing the need to include robustness criteria in the design process to develop complex active control designs for practical implementation. The paper investigates the crosswind aerodynamic parameters after the blocking phase of a two-dimensional square cross-section structure by measuring the response in wind tunnel tests under light wind flow conditions. To improve the accuracy of the results, the interpolation of the experimental curves in the time domain and the analytical responses were numerically optimized to finalize the results. Due to this combined effect, the three aerodynamic parameters decrease with increasing wind speed and asymptotically affect the upper branch constants. This means that the aerodynamic parameters along the density distribution are minimal. Taylor series are utilized to describe the fuzzy nonlinear plant and derive the stability analysis using polynomial function for analyzing the aerodynamic parameters and numerical simulations. Due to it will yield intricate terms to ensure stability criterion, therefore we aim to avoid kinds issues by proposing a polynomial homogeneous framework and utilizing Euler's functions for homogeneous systems. Finally, we solve the problem of stabilization under the consideration by SOS (sum of squares) and assign its fuzzy controller based on the feasibility of demonstration of a nonlinear system as an example.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.5
• /
• pp.414-420
• /
• 2010

This paper presents a validation of the accuracy of a reduced order model(ROM) and the efficiency of the design optimization using a Proper Orthogonal Decomposition(POD) to transonic wing/fuselage system. Three dimensional Euler equations are solved to extrude snapshot data of the full order aerodynamic analysis, and then a set of POD basis vectors reproducing the behavior of flow around the wing/fuselage system is calculated from these snapshots. In this study, reduced order model constructed through this procedure is applied to several validation cases, and then it is confirmed that the ROM has the capability of the prediction of flow field in the space of interest. Additionally, after the design optimization of the wing/fuselage system with the ROM is performed, results of the ROM are compared with results of the design optimization using response surface model(RSM). From these, it can be confirmed that the design optimization with the ROM is more efficient than RSM.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.10
• /
• pp.939-946
• /
• 2008

Euler codes or Navier-Stokes codes for compressible flows suffer severe degradation in convergence as Mach number approaches zero. The convergence problem arose from the wide disparity in characteristic speeds can be solved using preconditioning methods without large modifications. In this paper, a preconditioned RANS(Reynolds Averaged Navier-Stokes) solver is developed for analysis of low Mach number flows. In order to validate the method, computational examples are chosen and the results are compared with the experimental data and the existing computed results showing a good accuracy and convergence characteristics for steady inviscid, laminar and turbulent flows at low Mach number.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.169-172
• /
• 2002

When a plane shockwave reflects ken a concave wall, it is focused at a certain location, resulting in extremely high local pressure and temperature. This focusing is due to a nonlinear phenomenon of shock wave. The focusing phenomenon has been extensively applied to many diverse folds of engineering and medical treatment as well. In the current study, the focusing of shock wave over a reflector is numerically investigated using a CFD method. The Harten-Yee total variation diminishing (TVD) scheme is used to solve the unsteady, two-dimensional, compressible, Euler equations. The incident shock wave Mach number $M_{s}\;of\;1.1{\~}l.3$ is applied to the parabolic reflectors with several different depths. Detailed focusing characteristics of the shock wave are investigated in terms of peak pressure, gasdynamic and geometrical foci. The results obtained are compared with the previous experimental results. The results obtained show that the peak pressure of shock wave focusing and its location strongly depend on the magnitude of the incident shock wave and depth of parabolic reflector. It is also found that depending up on the depth of parabolic reflector, the weak shock wave focusing process can classified into three distinct patterns : the reflected shock waves do not intersect each other before and after focusing, the reflected shock waves do not intersect each other before focusing, but intersect after focusing, and the reflected shock waves intersect each other before and after focusing. The predicted Schlieren images represent the measured shock wave focusing with a good accuracy.

• Journal of computational fluids engineering
• /
• v.21 no.3
• /
• pp.15-23
• /
• 2016

In the present study, two phase flows around a projectile vertically launched from an underwater platform have been numerically investigated by using a three dimensional multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. The relative motion between the platform and projectile was described by six degrees of freedom equations of motion with Euler angles and a chimera technique. The propulsive power of the projectile was modeled as the fluid force acting on the lower surface of the body by the compressed air emitted from the underwater platform. Various flow conditions were considered to analyze the fluid-dynamics motion parameters of the projectile. The water level of platform and the current speed around the projectile were the main parametric variables. The numerical calculations were conducted up to 0.75sec in physical time scale. The dynamics tendency of the projectile was almost identical with respect to the water level variation due to the constant buoyancy term. The moving speed of the projectile along the vertical axis inside the platform decreased when the current speed increased. This is because the inflow from outside of the platform impeded development of the compressed air emitted from the floor surface of the launch platform. As a result, the fluid force acting on the lower surface of the projectile decreased, and injection time of the projectile from the platform was delayed.

• Journal of Korean Society of Environmental Engineers
• /
• v.33 no.7
• /
• pp.536-543
• /
• 2011

In order to find the best design of artificial stream for the riverbed filtration in multi-purpose filtration pond, a mathematical model was developed employing the energy line and the Manning's formula and was analyzed by the Euler's technique. Various design factors were investigated through scenario analyses of the artificial stream using the model. Results showed that the appropriate slope of the stream bottom was 2/10,000 and the appropriate infiltration rate at the streambed was $2.5m^3/m^2-day$ for the pond with the area of 100 ha, and that the Manning's roughness coefficient in this case was expected to be about 0.026 and the maximum water-depth was less than 1m. It was also shown that the longer the artificial stream the more advantageous it became for the riverbed filtration. Furthermore, results showed that it was not an efficient way to prevent clogging of the streambed by increasing the flow velocity of the stream and that the performance was higher near a weir with a large head drop.