• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.21 no.8
• /
• pp.698-707
• /
• 2011

This paper is to analyze the stability of the thick plate on inhomogeneous Pasternak foundation, with linearly varying thickness and concentrated mass by finite element method. To verify this finite element method, the results of natural frequencies and buckling stresses by the proposed method are compared with the existing solutions. The dynamic instability regions are decided by the dynamic stability analysis of the thick plate on inhomogeneous Pasternak foundation, with linearly varying thickness and concentrated mass. The non-dimensional Winkler foundation parameter is applied as 100, 1000 and non-dimensional shear foundation parameter is applied as 5. The tapered ratios are applied as 0.25 and 1.0, the ratios of concentrated mass to plate mass as 0.25 and 1.0, and the ratio of in-plane force to critical load as 0.4. As the result of numerical analysis of the thick plate on inhomogeneous Pasternak foundation for $u{\times}v=300cm{\times}300cm$ and $a{\times}b=600cm{\times}600cm$, instability areas of the thick plate which has the larger rigidity of inner area are farther from ${\beta}$-axis and narrower than those which has the larger rigidity of outer area.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.5
• /
• pp.695-702
• /
• 2002

Nonlinear dynamics of striped diffusion flames, by the diffusional-thermal instability with Lewis numbers sufficiently less than unity, is numerically investigated by examining various two-dimensional flame-structure solutions. The Lewis numbers for fuel and oxidizer are assumed to be identical and an overall single-step Arrhenius-type chemical reaction rate is employed in the model. Particular attention is focused on identifying the flame-stripe solution branches corresponding to each distinct stripe pattern and hysteresis encountered during the transition. At a Damkohler number slightly greater than the extinction Damkohler number, eight-stripe solution first emerges from one dimensional solution. The eight-stripe solution survives Damkohler numbers much smaller than the extinction Damkohler number until the transition to four-stripe solution occurs at the first forward transition Damkohler number. At the second forward transition Damkohler number, somewhat smaller than the first transition Damkohler number, the transition to two-stripe solution occurs. However, anu further transition from two-stripe solution to one-stripe solution is not always possible even if one-stripe solution can be independently accessed for particular initial conditions. The Damkohler number ranges for two-stripe and one-stripe solutions are found to be virtually identical because each stripe is an independent structure if distance between stripes is sufficiently large. By increasing the Damkohler number, the backward transition can be observed. In comparison with the forward transition Damkohler numbers, the corresponding backward transition Damkohler numbers are always much greater, thereby indicating significant hysteresis between the stripe patterns of strained diffusion flames.

• Structural Engineering and Mechanics
• /
• v.61 no.4
• /
• pp.497-510
• /
• 2017

A framed structure may be composed of two sub-structures, which are linked by a hinged joint. One sub-structure is the primary system and the other is the secondary system. The primary system, which is subjected to the periodic external load, can give rise to an auto-parametric resonance of the second system. Considering the geometric-stiffness effect produced by the axially internal force, the element equation of motion is derived by the extended Hamilton's principle. The element equations are then assembled into the global non-homogeneous Mathieu-Hill equations. The Newmark's method is introduced to solve the time-history responses of the non-homogeneous Mathieu-Hill equations. The energy-growth exponent/coefficient (EGE/EGC) and a finite-time Lyapunov exponent (FLE) are proposed for determining the auto-parametric instability boundaries of the structural system. The auto-parametric instabilities are numerically analyzed for the two frames. The influence of relative stiffness between the primary and secondary systems on the auto-parametric instability boundaries is investigated. A phenomenon of the "auto-parametric internal resonance" (the auto-parametric resonance of the second system induced by a normal resonance of the primary system) is predicted through the two numerical examples. The risk of auto-parametric internal resonance is emphasized. An auto-parametric resonance experiment of a ${\Gamma}$-shaped frame is conducted for verifying the theoretical predictions and present calculation method.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.39 no.3
• /
• pp.232-241
• /
• 2011

The hybrid rocket shows interesting characteristics of complicated mixing layer developed by the interaction between turbulent oxidizer flow and injected surface mass flow from fuel vaporization. In this study, the compressible LES was conducted to explore the physical phenomena of surface oscillatory flow induced by the flow interferences in a duct domain. From the numerical results, the wall injection generates the stronger streamwise vorticites and the negative components of axial velocity accompanied with the azimuthal vorticity near the surface. And the vortex shedding with a certain time scale was found to be developed by hydrodynamic instability in the mixing layer. The pressure fluctuations in this calculation exhibit a peculiar peak at a specific angular frequency($\omega$=8.8) representing intrinsic oscillation due to the injection.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.8 no.3
• /
• pp.437-450
• /
• 1996

A set of stability equations is formulated for natural convection flows adjacent to a vertical isothermal surface melting in cold pure water. It takes account of the nonparallelism of the base flows. The melting rate is regarded as a blowing velocity at the ice surface. The numerical solutions of the linear stability equations which constitute a two-point boundary value problem are accurately obtained for various values of the density extremum parameter $R=(T_m-T_{\infty})/(T_0-T_{\infty})$ in the range $0.3{\leq}R{\leq}0.6$, by using a computer code COLNEW. The blowing effects on the base flow becomes more significant as ambient temperature ($T_{\infty}$) increases to $T_{\infty}=10^{\circ}C$. The maximum decrease of heat transfer rate is about 6.4 percent. The stability results show that the melting at surface causes the critical Grashof number $G^*$ and the maximum frequency of disturbances to decrease. In comparision with the results for the conventional parallel flow model, the nonparallel flow model has a higher critical Grashof number but has lower amplification rates of disturbances than does the parallel flow model. The spatial amplification contours exhibit that the selective frequency $B_0$ of the nonparallel flow model is higher than that of the parallel flow model and that the effects of melting are rather small. The present study also indicates that the selective frequency $B_0$ can be easily predicted by the value of the frequency parameter $B^*$ at $G^*$, which comes from the neutral stability results of the nonparallel flow model.

• Wind and Structures
• /
• v.25 no.6
• /
• pp.507-535
• /
• 2017

In the strong wind shutdown state, the blade position significantly affects the streaming behavior and stability performance of wind turbine towers. By selecting the 3M horizontal axis wind turbine independently developed by Nanjing University of Aeronautics and Astronautics as the research object, the CFD method was adopted to simulate the flow field of the tower-blade system at eight shutdown positions within a single rotation period of blades. The effectiveness of the simulation method was validated by comparing the simulation results with standard curves. In addition, the dynamic property, aerostatic response, buckling stability and ultimate bearing capacity of the wind turbine system at different shutdown positions were calculated by using the finite element method. On this basis, the influence regularity of blade shutdown position on the wind-induced response and stability performance of wind turbine systems was derived, with the most unfavorable working conditions of wind-induced buckling failure of this type of wind turbines concluded. The research results implied that within a rotation period of the wind turbine blade, when the blade completely overlaps the tower (Working condition 1), the aerodynamic performance of the system is the poorest while the aerostatic response is relatively small. Since the influence of the structure's geometrical nonlinearity on the system wind-induced response is small, the maximum displacement only has a discrepancy of 0.04. With the blade rotating clockwise, its wind-induced stability performance presents a variation tendency of first-increase-then-decrease. Under Working condition 3, the critical instability wind speed reaches its maximum value, while the critical instability wind speed under Working condition 6 is the smallest. At the same time, the coupling effect between tower and blade leads to a reverse effect which can significantly improve the ultimate bearing capacity of the system. With the reduction of the area of tower shielded by blades, this reverse effect becomes more obvious.

• Nuclear Engineering and Technology
• /
• v.16 no.3
• /
• pp.125-130
• /
• 1984

In the application of the RETRAN code to the analysis of very small LOCA one of main concerns is placed on use of the bubble rise model in the upper plenum, because the bubble rise model nay cause a numerical divergence problem and coefficients used to describe it are based on experimental results of large LOCA. In order to solve this problem, a method, which enables us to predict the mixture level in the upper plenum without use of the bubble rise model, was proposed. For this method the local void distribution in the core and upper plenum was derived by using a simplified slip model. It was shown that results predicted from the derived equation are in excellent agreement with experimental data. Additionally it was found that local void in the upper plenum has a uniform distribution unlike a linear distribution in large LOCA. Communication between the upper plenum and upper head was investigated. By introducing the concept of Taylor instability, it was proved that counter-current Hon between them is possible.

• Tribology and Lubricants
• /
• v.39 no.4
• /
• pp.139-147
• /
• 2023

Cavitation phenomena observed during the operation of a submerged plain journal bearing (PJB) can affect bearing performance parameters such as dynamic coefficients, whirl frequency ratio, and critical mass. This study presents numerical solutions of the Reynolds equation for steadily and dynamically loaded submerged PJBs with half-Sommerfeld (HS), Reynolds, and Jakobsson-Floberg-Olsson (JFO) cavitation models when the supply pressure is larger or equal to the cavitation pressure. The loads at various eccentricity ratios are identical; however, the attitude angle is approximately 6% smaller when the eccentricity ratio is between 0.2 and 0.7 and the JFO model is used, compared to that when the Reynolds model is used. Dynamic coefficients obtained with the HS and Reynolds model show good agreement with each other, except for k_xz, which is sensitive to changes in the force normal to the rotor weight, and is attributed to the difference in the attitude angle obtained with each cavitation model. Stiffness coefficients are determined using the pressure distribution in the film, and therefore, when the JFO model is used, the direct stiffness coefficients are affected and show opposite signs for most eccentricity ratios. The mass-conservative JFO model can predict at least a 30% smaller critical mass compared to that using the HS and Reynolds models. Thus, the instability analysis results can change based on the cavitation model used in a submerged PJB. The results of this research indicate that the JFO model should be used when designing a rotor system supported by submerged PJBs.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.4
• /
• pp.343-351
• /
• 2016

Cryogenic spray characteristics of a nitrogen swirl injector operating in supercritical environment have been numerically investigated. By comparing the equation of states(EOS) used for supercritical condition, SRK EOS was applied to predict the nitrogen thermodynamic property under supercritical environment. A Chung's method was implemented for the calculation of viscosity and conductivity and Takahashi's correlation based on Fuller's Theorem was implemented for the calculation of diffusion coefficient. By injecting the nitrogen with 5 bar differential pressure into 50 bar chamber filled with nitrogen, numerical simulation has been conducted. The dynamic Smagorinsky sub-grid scale (SGS) model has been compared with the algebraic Smagorinsky SGS model using FFT frequency analysis. The instability at the liquid film and gas core inside injector and the propagation of pressure oscillation into the injector has been investigated. The spreading angle of swirl injector obtained by numerical calculation has been validated with experimental result.

• Wind and Structures
• /
• v.10 no.4
• /
• pp.367-382
• /
• 2007

A nonlinear numerical method was developed to assess the stability of suspension bridge catwalks under a wind load. A section model wind tunnel test was used to obtain a catwalk's aerostatic coefficients, from which the displacement-dependent wind loads were subsequently derived. The stability of a suspension bridge catwalk was analyzed on the basis of the geometric nonlinear behavior of the structure. In addition, a full model test was conducted on the catwalk, which spanned 960 m. A comparison of the displacement values between the test and the numerical simulation shows that a numerical method based on a section model test can be used to effectively and accurately evaluate the stability of a catwalk. A case study features the stability of the catwalk of the Runyang Yangtze suspension bridge, the main span of which is 1490 m. Wind can generally attack the structure from any direction. Whenever the wind comes at a yaw angle, there are six wind load components that act on the catwalk. If the yaw angle is equal to zero, the wind is normal to the catwalk (called normal wind) and the six load components are reduced to three components. Three aerostatic coefficients of the catwalk can be obtained through a section model test with traditional test equipment. However, six aerostatic coefficients of the catwalk must be acquired with the aid of special section model test equipment. A nonlinear numerical method was used study the stability of a catwalk under a yaw wind, while taking into account the six components of the displacement-dependent wind load and the geometric nonlinearity of the catwalk. The results show that when wind attacks with a slight yaw angle, the critical velocity that induces static instability of the catwalk may be lower than the critical velocity of normal wind. However, as the yaw angle of the wind becomes larger, the critical velocity increases. In the atmospheric boundary layer, the wind is turbulent and the velocity history is a random time history. The effects of turbulent wind on the stability of a catwalk are also assessed. The wind velocity fields are regarded as stationary Gaussian stochastic processes, which can be simulated by a spectral representation method. A nonlinear finite-element model set forepart and the Newmark integration method was used to calculate the wind-induced buffeting responses. The results confirm that the turbulent character of wind has little influence on the stability of the catwalk.