• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.3
• /
• pp.467-476
• /
• 2000

An unsteady numerical simulation was performed for locally-forced separated and reattaching flow over a backward-facing step. The local forcing was given to the separated and reattaching flow by means of a sinusoidally oscillating jet from a separation line. A version of the $k-{\varepsilon}-f_{\mu}$ model was employed, in which the near-wall behavior without reference to distance and the nonequilibrium effect in the recirculation region were incorporated. The Reynolds number based on the step height (H) was fixed at $Re_H=33000$, and the forcing frequency was varied in the range $0{\leq}St_H{\leq}2$. The predicted results were compared and validated with the experimental data of Chun and Sung. It was shown that the unsteady locally-forced separated and reattaching flows are predicted reasonably well with the $k-{\varepsilon}-f_{\mu}$ model. To characterize the large-scale vortex evolution due to the local forcing, numerical flow visualizations were carried out.

• Wind and Structures
• /
• v.25 no.6
• /
• pp.589-608
• /
• 2017

A novel probabilistic approach is presented for estimating the equivalent static wind loads that produce a static response of the structure, which is "equivalent" in a probabilistic sense, to the extreme dynamic responses due to the unsteady pressure random field induced by the wind. This approach has especially been developed for complex structures (such as stadium roofs) for which the unsteady pressure field is measured in a boundary layer wind tunnel with a turbulent incident flow. The proposed method deals with the non-Gaussian nature of the unsteady pressure random field and presents a model that yields a good representation of both the quasi-static part and the dynamical part of the structural responses. The proposed approach is experimentally validated with a relatively simple application and is then applied to a stadium roof structure for which experimental measurements of unsteady pressures have been performed in boundary layer wind tunnel.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.4
• /
• pp.839-846
• /
• 1994

This paper presents a simultaneous solution algorithm for one-dimensional unsteady flow routing through a dendritic channel system. This simulations solution algorithm is based on the double-sweep method and utilizes separate recursion equations for continuity, momentum and energy equations for each of the individual components of a dendritic channel system. Through separate recursion equations for each of the components. the new algorithm converts a dendritic channel network problem into a single-channel problem. The new algorithm is utilized in conjunction with a linearized unsteady flow model using full dynamic flow equations. The required computer storage for the coefficient matrix of the whole system is reduced significantly from the $2N{\times}2N$ matrix to a $2N{\times}4$ matrix, where N is the number of cross sections used in the computation of flow variables in a dendritic channel system. The algorithm presented in this paper provides an efficient and accurate modeling of unsteady flow events through a dendritic channel system.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.6 no.1
• /
• pp.1-11
• /
• 2002

The suitable turbulence model is found to be required in the course of establishing a proper analysis methodology for thermal stripping phenomena which are shown in strong temperature variation area such as reactors and propulsion devices. Three different turbulence models of $\kappa$-$\varepsilon$ model, modified $\kappa$-$\varepsilon$ model, and full Reynolds stress(FRS) model, are applied to analyze unsteady turbulent flows with temperature variation. Three test cases are selected for verification. These are vertical jet flows with water and sodium, and parallel jet flow with sodium. Analysis yields the conclusion that 3-D computation with FRS betters others. However, modified modeling is required to improve its heat transfer characteristic analysis. Further analysis is performed to find momentum variation effects on temperature distribution. It is found that the momentum increase results increase of fluid mixing and magnitude of temperature variation.

• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.7-13
• /
• 2010

Flight vehicles such as wheel wells and bomb bays have many cavities. The flow around a cavity is characterized as an unsteady flow because of the formation and dissipation of vortices brought by the interaction between the free stream shear layer and the internal flow of the cavity. The resonance phenomena can damage the structures around the cavity and negatively affect the aerodynamic performance and stability of the vehicle. In this study, a numerical analysis was performed for the cavity flows using the unsteady compressible three-dimensional Reynolds-Averaged Navier-Stokes (RANS) equation with Wilcox's turbulence model. The Message Passing Interface (MPI) parallelized code was used for the calculations by PC-cluster. The cavity has aspect ratios (L/D) of 5.5 ~ 7.5 with width ratios (W/D) of 2 ~ 4. The Mach and Reynolds numbers are 0.4 ~ 0.6 and $1.6{\times}10^6$, respectively. The occurrence of oscillation is observed in the "shear layer and transient mode" with a feedback mechanism. Based on the Sound Pressure Level (SPL) analysis of the pressure variation at the cavity trailing edge, the dominant frequencies are analyzed and compared with the results of Rossiter's formula. The dominant frequencies are very similar to the result of Rossiter's formula and other experimental datum in the low aspect ratio cavity (L/D = ~4.5). In the high aspect ratio cavity, however, there are other low dominant frequencies of the leading edge shear layer with the dominant frequencies of the feedback mechanism.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.13-18
• /
• 2009

Flight vehicles such as wheel wells and bomb bays have many cavities. The flow around a cavity is characterized as an unsteady flow because of the formation and dissipation of vortices brought about by the interaction between the free stream shear layer and the internal flow of the cavity. The resonance phenomena can damage the structures around the cavity and negatively affect the aerodynamic performance and stability of the vehicle. In this study, a numerical analysis was performed for the cavity flows using the unsteady compressible three-dimensional Reynolds-Averaged Navier-Stokes (RANS) equation with Wilcox's turbulence model. The Message Passing Interface (MPI) parallelized code was used for the calculations by PC-cluster. The cavity has aspect ratios (L/D) of 2.5 ~ 7.5 with width ratios (W/D) of 2 ~ 4. The Mach and Reynolds numbers are 0.4 ~ 0.6 and $1.6{\times}106$, respectively. The occurrence of oscillation is observed in the "shear layer and transient mode" with a feedback mechanism. Based on the Sound Pressure Level (SPL) analysis of the pressure variation at the cavity trailing edge, the dominant frequencies are analyzed and compared with the results of Rossiter's formula. The dominant frequencies are very similar to the result of Rossiter's formula and other experimental data in the low aspect ratio cavity (L/D = ~ 4.5). In the large aspect ratio cavity, however, there are other low dominant frequencies due to the leading edge shear layer with the dominant frequencies of the feedback mechanism. The characteristics of the acoustic wave propagation are analyzed using the Correlation of Pressure Distribution (CPD).

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.12 no.1
• /
• pp.42-47
• /
• 2002

Sirocco fans are widely used in HVAC and air conditioning systems, and the noise generated by these machines causes one of the most serious problems. In general, the sirocco fan noise is often dominated by tones at BPF(blade passage frequency) and broadband noise. However, only a few researches have been carried out on predicting the aeroacoustic noise because of the difficulty in obtaining detailed information about the flow field and casing effects on noise radiation. The objective of this study is to develop a prediction method for the unsteady flow field and the acoustic pressure field of a sirocco fan. We assume that the impeller rotates with a constant angular velocity and the flow field around the impeller is incompressible and inviscid. So, a discrete vortex method (DVM) is used to model the centrifugal fan and to calculate the flow field. The force of each element on the blade is calculated by the unsteady Bernoulli equation. Lowson\`s method is used to predict the acoustic source. Reasonable results are obtained not only fur the tonal noise but also far the amplitudes of the broadband noise. Acoustic pressure is proportional to (Ω)2.3, which is the similar value with the measured data.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.8 s.251
• /
• pp.796-805
• /
• 2006

In this study, the problem of using the mixing index as a measure of the mixing performance for a certain flow field has been discussed. The flow model subjected to this study is the two-dimensional unsteady lid-driven cavity flow. The transport equation for the concentration within the cavity was solved by using the finite volume method where the convective terms are discretized with the central difference scheme. It was shown that both the concentration dispersion and the mixing index depend highly on the initial distribution of the concentration, and therefore the mixing index obtained from the concentration dispersion equation loses its universal applicability.

• 유체기계공업학회:학술대회논문집
• /
• 1999.12a
• /
• pp.208-217
• /
• 1999

The present work describes the prediction method for the unsteady flow field and the acoustic pressure field of a ducted axial fan. The prediction method is comprised of time-marching free-wake method, acoustic analogy, and the Helmholtz-Kirchhoff BEM. The predicted sound signal of a rotor is similar to the experiment one. We assume that the rotor rotates with a constant angular velocity and the flow field around the rotor is incompressible and inviscid. Then, a time-marching free-wake method is used to model the fan and to calculate the flow field. The force of each element on the blade is calculated by the unsteady Bernoulli equation. Lowson's method is used to predict the acoustic source. The newly developed Helmholtz-Kirchhoff BEM for thin body is used to calculate the sound field of the ducted fan. The ducted fan with 6 blades is analysed and the sound field around the duct is calculated.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1997.04a
• /
• pp.43-51
• /
• 1997

Steady and unsteady numerical simulations are conducted for the experiments performed to investigate the ram accelerator flow field by using the expansion tube facility in Stanford University. Navier-Stokes equations for chemically reacting flows are analyzed by fully implicit and time accurate numerical methods with Jachimowski's detailed chemistry model for hydrogen-air combustion involving 9 species and 19 reaction steps. Although the steady state assumption shows a good agreement with the experimental schlieren and OH PLIF images for the case of $2H_2$+$O_2$+$17N_2$, it fails in reproducing the combustion region behind the shock intersection point shown in the case of $2H_2$+$O_2$+$12N_2$, mixture. Therefore, an unsteady numerical simulation is conducted for this case and the result shows all the detailed flow stabilization process. The experimental result is revealed to be an instantaneous result during the flow stabilization process. The combustion behind the shock intersection point is the result of a normal detonation formed by the intersection of strong oblique shocks that exist at early stage of the stabilization process. At final stage, the combustion region behind the shock intersection point disappears and the steady state result is retained. The time required for stabilization of the reacting flow in the model ram accelerator is found to be very long in comparison with the experimental test time.