• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.5
• /
• pp.1158-1174
• /
• 1988

A numerical technique is developed for the solution of fully developed turbulent recirculating flow in the passage of variable area using the non-orthogonal coordinate transformation. In the numerical analysis, primitive pressure-velocity finite difference equations were solved by SIMPLER algorithm with 2-equation turbulence model and algebraic stress model (ASM). QUICK scheme on the differencing of convective terms which is free from the inaccuracies of numerical diffusion has been applied to the variable grids and the results compared with those from HYBRID scheme. In order to test the effect of streamline curvatures on turbulent diffusion Lee and Choi streamline curvature correction model which has been obtained by modifying the Leschziner and Rodi's model is testes. The ASM was also employed and the results are compared to those from another turbulence model. The results show that difference of convective differencing schemes and turbulence models give significant differences in the prediction of velocity fields in the expansion region and outlet region of the combustor, however show little differences in the parallel flow region.

• Advances in aircraft and spacecraft science
• /
• v.5 no.3
• /
• pp.363-383
• /
• 2018

One-dimensional (1D) models of incompressible flows, can be of interest for many applications in which fast resolution times are demanded, such as fluid-structure interaction of flows in compliant pipes and hemodynamics. This work proposes a higher-order 1D theory for the flow-field analysis of incompressible, laminar, and viscous fluids in rigid pipes. This methodology is developed in the domain of the Carrera Unified Formulation (CUF), which was first employed in structural mechanics. In the framework of 1D modelling, CUF allows to express the primary variables (i.e., velocity and pressure fields in the case of incompressible flows) as arbitrary expansions of the generalized unknowns, which are functions of the 1D computational domain coordinate. As a consequence, the governing equations can be expressed in terms of fundamental nuclei, which are invariant of the theory approximation order. Several numerical examples are considered for validating this novel methodology, including simple Poiseuille flows in circular pipes and more complex velocity/pressure profiles of Stokes fluids into non-conventional computational domains. The attention is mainly focused on the use of hierarchical McLaurin polynomials as well as piece-wise nonlocal Lagrange expansions of the generalized unknowns across the pipe section. The preliminary results show the great advantages in terms of computational costs of the proposed method. Furthermore, they provide enough confidence for future extensions to more complex fluid-dynamics problems and fluid-structure interaction analysis.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.743-750
• /
• 2008

Asymmetric tip clearance in an axial compressor induces pressure and velocity redistributions along the circumferential direction in an axial compressor. This paper presents the mechanism of the flow redistribution due to the asymmetric tip clearance with a simple numerical modeling. The flow field of a rotor of an axial compressor is predicted when an asymmetric tip clearance occurs along the circumferential direction. The modeling results are supported by CFD results not only to validate the present modeling but also to investigate more detailed flow fields. Asymmetric tip clearance makes local flow area and resultant axial velocity vary along the circumferential direction. This flow redistribution 'seed' results in a different flow patterns according to the flow coefficient. Flow field redistribution patterns are largely dependent on the local tip clearance performance at low flow coefficients. However, the contribution of the main flow region becomes dominant while the tip clearance effect becomes weak as the flow coefficient increases. The flow field redistribution pattern becomes noticeably strong if a blockage effect is involved when the flow coefficient increases. The relative flow angle at the small clearance region decreases which result in a negative incidence angle at the high flow coefficient. It causes a recirculation region at the blade pressure surface which results in the flow blockage. It promotes the strength of the flow field redistribution at the rotor outlet. These flow pattern changes have an effect on the blade loading perturbations. The integration of blade loading perturbation from control volume analysis of the circumferential momentum leads to well-known Alford's force. Alford's force is always negative when the flow blockage effects are excluded. However when the flow blockage effect is incorporated into the modeling, main flow effects on the flow redistribution is also reflected on the Alford's force at the high flow coefficient. Alford's force steeply increases as the flow coefficient increases, because of the tip leakage suppression and strong flow redistribution. The predicted results are well agreed to CFD results by Kang and Kang(2006).

• Wind and Structures
• /
• v.19 no.4
• /
• pp.421-441
• /
• 2014

With the assistance of typhoon field data at aerial elevation level observed by meteorological satellites and wind velocity and direction records nearby the ground gathered in Guangzhou Weather Station between 1985 and 2001, some key wind field parameters under typhoon climate in Guangzhou region were calibrated based on Monte-Carlo stochastic algorithm and Meng's typhoon numerical model. By using Peak Over Threshold method (POT) and Generalized Pareto Distribution (GPD), Wind field characteristics during typhoons for various return periods in several typical engineering fields were predicted, showing that some distribution rules in relation to gradient height of atmosphere boundary layer, power-law component of wind profile, gust factor and extreme wind velocity at 1-3s time interval are obviously different from corresponding items in Chinese wind load Codes. In order to evaluate the influence of typhoon field parameters on long-span flexible bridges, 1:100 reduced-scale wind field of type B terrain was reillustrated under typhoon and normal conditions utilizing passive turbulence generators in TJ-3 wind tunnel, and wind-induced performance tests of aero-elastic model of long-span Guangzhou Xinguang arch bridge were carried out as well. Furthermore, aerodynamic admittance function about lattice cross section in mid-span arch lib under the condition of higher turbulence intensity of typhoon field was identified via using high-frequency force-measured balance. Based on identified aerodynamic admittance expressions, Wind-induced stochastic vibration of Xinguang arch bridge under typhoon and normal climates was calculated and compared, considering structural geometrical non-linearity, stochastic wind attack angle effects, etc. Thus, the aerodynamic response characteristics under typhoon and normal conditions can be illustrated and checked, which are of satisfactory response results for different oncoming wind velocities with resemblance to those wind tunnel testing data under the two types of climate modes.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.15-27
• /
• 2022

This work comes within the framework of the "Benchmark on the Aerodynamics of a Rectangular Cylinder" that investigates a rectangular cylinder of length-to-depth ratio equal to 5. The present study reports and discusses velocity fields acquired using planar Particle Image Velocitmetry for several angles of attack and Reynolds numbers. In particular, for a cylinder depth-based Reynolds number of 2 × 10⁴ and zero incidence angle, the flow features along the lateral (parallel to the freestream) upper and lower surfaces of the cylinder are reported. Using first and second order statistics of the velocity field, the main flow features are discussed, especially the size and location of the time-averaged flow structures and the distribution of the Reynolds stresses. The variation of the flow features with the incidence is also studied considering angles of attack up to 6°. It is shown that the time-averaged flow is fully detached for incidence higher than 2°. For an angle of attack of 0°, the effects of the Reynolds number varying between 5 × 10³ and 2 × 10⁴ are investigated looking at flow statistics. It is shown that the time-averaged location of the reattachment point and the shape and position of the time-averaged main vortex are mostly constant with the Reynolds number. However, the size of the inner region located below the time-averaged shear layer and just downstream the leading edge corner appears to be strongly dependent on the Reynolds number.

• Journal of Ocean Engineering and Technology
• /
• v.36 no.3
• /
• pp.143-152
• /
• 2022

To reach a port, a ship must pass through a shallow water zone where seabed effects alter the hydrodynamics acting on the ship. This study examined the maneuvering characteristics of an autonomous surface ship at 3-DOF (Degree of freedom) motion in deep water and shallow water based on the in-port speed of 1.54 m/s. The CFD (Computational fluid dynamics) method was used as a specialized tool in naval hydrodynamics based on the RANS (Reynolds-averaged Navier-Stoke) solver for maneuvering prediction. A virtual captive model test in CFD with various constrained motions, such as static drift, circular motion, and combined circular motion with drift, was performed to determine the hydrodynamic forces and moments of the ship. In addition, a model test was performed in a square tank for a static drift test in deep water to verify the accuracy of the CFD method by comparing the hydrodynamic forces and moments. The results showed changes in hydrodynamic forces and moments in deep and shallow water, with the latter increasing dramatically in very shallow water. The velocity fields demonstrated an increasing change in velocity as water became shallower. The least-squares method was applied to obtain the hydrodynamic coefficients by distinguishing a linear and non-linear model of the hydrodynamic force models. The course stability, maneuverability, and collision avoidance ability were evaluated from the estimated hydrodynamic coefficients. The hydrodynamic characteristics showed that the course stability improved in extremely shallow water. The maneuverability was satisfied with IMO (2002) except for extremely shallow water, and collision avoidance ability was a good performance in deep and shallow water.

• Journal of Astronomy and Space Sciences
• /
• v.40 no.3
• /
• pp.101-111
• /
• 2023

We conducted a statistical study of polar mesospheric summer echoes (PMSEs) in relation to magnetic local time (MLT), considering the geomagnetic conditions using the K-index (or K). Additionally, we performed a case study to examine the velocity profile, specifically for high velocities (≥ ~100 m/s) varying with high temporal resolution at high K-index values. This study utilized the PMSE data obtained from the mesosphere-stratosphere-troposphere radar located in Esrange, Sweden (63.7°N, 21°E). The change in K-index in terms of MLT was high (K ≥ 4) from 23 to 04 MLT, estimated for the time PMSE was present. During the near-midnight period (0-4 MLT), both PMSE occurrence and signal-to-noise ratio (SNR) displayed an asymmetric structure with upper curves for K ≥ 3 and lower curves for K < 3. Furthermore, the occurrence of high velocities peaked at 3-4 MLT for K ≥ 3. From case studies focusing on the 0-3 MLT period, we observed persistent eastward-biased high velocities (≥ 200 m/s) prevailing for ~18 min. These high velocities were accompanied with the systematic motion of profiles at 85-88 km, including large shear formation. Importantly, the rapid variations observed in velocity could not be attributed to neutral wind effects. The present findings suggest a strong substorm influence on PMSE, especially in the midnight and early dawn sectors. The large zonal drift observed in PMSE were potentially energized by local electromagnetic fields or the global convection field induced by the electron precipitation during substorms.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.26 no.2
• /
• pp.94-101
• /
• 1989

Transport properties of 2DEG at AlGaAs/GaAs interface such as average electron energy, flight distance, each valley occupancy ratio, average electron velocity for various fields are investigated by MONTE-CARLO method. As the electric field increases, more electrons transit drastically from (000) valley to (000) upper valley. This phenomenon shows the nonstationary effect such as velocity overshoot. The duration of the transient decreases from about 1.4 psec for electric field E = 7KV/cm to about 0.7 psec for 12KV/cm. The average electron velocity during transient transport in 2DEG is about 8 times the steady-state velocity for E = 12KV/cm at room temperature. In comparison with bulk GaAs the peak velocity in the 2DEG is higher than that in even pure bulk GaAs at electric field E = 7 KV/cm. On the basis of the fact that the electrons in the 2DEG have larger peak velocity and shorter transient time of velocity than those in the bulk GaAs, it is suggested that the device with 2DEG may obtain higher mobility than that with bulk GaAs.

• Journal of Energy Engineering
• /
• v.6 no.1
• /
• pp.58-67
• /
• 1997

In this study numerical computations were performed to predict reacting flow fields of gasification processes of pulverized subbituminous coal in a cylindrical coal gasifier. To check the size effects of particles on flow fields in the gasifier, simulations were performed for five cases with four sizes of particles such as 40 $\mu\textrm{m}$, 60 $\mu\textrm{m}$, 80 $\mu\textrm{m}$ and 100 $\mu\textrm{m}$. Each case has a unique size of particles with one more case that has evenly mixed four sizes of particles. Predictions showed that the gasification which uses coals of the mixed sizes reveals more preferable gas velocity and temperature distributions than that uses coals of a unique size. Predicted gas temperature at the exit of the gasifier ranged 1,400 to 1580$^{\circ}C$, 1,480 to 1,700$^{\circ}C$, 1,600 to 1740$^{\circ}C$, 1630 to 1790$^{\circ}C$ and 1500 to 1680$^{\circ}C$ for particle sizes of 40 $\mu\textrm{m}$, 60 $\mu\textrm{m}$, 80 $\mu\textrm{m}$, 100 $\mu\textrm{m}$ and the evenly mixed four sizes, respectively.

• Korean Journal of Veterinary Research
• /
• v.38 no.4
• /
• pp.898-908
• /
• 1998

The aim of this experiments were to assess the time-interval change of motional characteristics in frozen-thawed semen of Korean native cattle (KNC) by using computer aided semen analysis (CASA) technology. Twenty-six KNC frozen semen straws were obtained from Korean KNC improvement department, livestock improvement main division, national livestock cooperatives federation in Korea. Specimens were allowed to thaw at $37^{\circ}C$ for 30 sec in water bath. Semen analysis was performed on semen image analysis system (SIAS, Medical supply, Korea) adjusted to the gate settings and used the semen droplet ($5{\mu}l$) placed on Makler counting chamber (Sefi medical instrument, Israel) prewarmed at $37^{\circ}C$. The same person used the same micropipette to fill the Makler counting chamber. A total of 150 or more of sperms were analysed in each specimen by a single trained person by scanning at least 5 to 10 fields. The measurement parameters in SIAS were as follows ; frame rate = 30 frames per sec, image capture = 1 sec, minimum motile speed = $10{\mu}m/s$, maximum countable sperm number = 400. Statistical analysis was done by Student t-test with use of the Sigma plot program on a IBM personal computer. The dancemean(DNM) and hyperactivated sperm(HYP) of frozen-thawed KNC semen kinematics were significantly decreased(p < 0.05) after 10 min of incubation at $37^{\circ}C$ water bath. But, wobble(WOB) of same sample semen was significantly increased(p < 0.05) after 10 min of incubation and significantly decrease(p < 0.05) after 60 min of same incubation. And, after 30 mim of incubation, significantly differences were found most of motion kinematics, motifity(MOT), curvilinear velocity(VCL), straight line velocity(VSL), average path velocity(VAP), amplitude of lateral head displacement(ALH), beat cross frequency(BCF), mean angular displacement(MAD), dance(DNC), on same sample semen. The DNM of KNC semen sample was variable kinematics after 30 min of incubation. Also, the linearity(LIN) and straightness(STR) was significantly decreased(p < 0.05) from 60 min of incubation. In conclusion, the AI within 30 min after thawing of frozen semen can be an effective method for obtaining high fertility rate in KNC reproductive program.