• Journal of computational fluids engineering
• /
• v.11 no.1 s.32
• /
• pp.57-66
• /
• 2006

An optimal shape design approach is presented for a subsonic S-shaped intake using aerodynamic sensitivity analysis. Two-equation turbulence model is employed to capture strong counter vortices in the S-shaped duct more precisely. Sensitivity analysis is performed for the three-dimensional Navier-Stokes equations coupled with two-equation turbulence models using a discrete adjoint method For code validation, the result of the flow solver is compared with experiment data and other computational results of bench marking test. To study the influence oj turbulence models and grid refinement on the duct flow analysis, the results from several turbulence models are compared with one another and the minimum number of grid points, which can yield an accurate solution is investigated The adjoint variable code is validated by comparing the complex step derivative results. To realize a sufficient and flexible design space, NURBS equations are introduced as a geometric representation and a new grid modification technique, Least Square NURBS Grid Approximation is applied With the verified flow solver, the sensitivity analysis code and the geometric modification technique, the optimization of S-shaped intake is carried out and the enhancement of overall intake performance is achieved The designed S-shaped duct is tested in several off-design conditions to confirm the robustness of the current design approach. As a result, the capability and the efficiency of the present design tools are successfully demonstrated in three-dimensional highly turbulent internal flow design and off-design conditions.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.21-24
• /
• 2008

The present study numerically investigates the glottal airflow characteristics as well as acoustic features of phonation fully coupled with dynamic behavior of vocal folds. The vocal folds are described by a low-dimensional body-covered model characterized by bio-mechanical parameters such as glottal width, vocal folds stiffness, and subglottal pressure. The flow in the vocal tract is modeled as an incompressible, axisymmetric form of the Navier-Stokes equations (INS), while the acoustic field is predicted by the linearized perturbed compressible equations (LPCE). The computed result shows that a two-mass model of vocal folds is sufficient to reproduce temporal variations in oral airflow and glottis motion produced by female speakers. It is also found that i) the glottal width has a significant effect on the amplitude of glottal flow, and thus on the amplitude of acoustic wave in the vocal tract, ii) the vocal fold tension is the main control parameter for the fundamental frequency of phonation, iii) the subglottal pressure plays an appreciable role on reproduction of the self-sustained oscillation of vocal folds, and iv) the strength of pulsating airflow and vortical structures are primarily affected by glottal width and subglottal pressure, and are closely related to pitch, loudness, and voice quality. Finally, more comprehensive explanation about the difference between one- and two-mass models is presented with discussion of effectiveness of vocal folds oscillation and voice quality.

• 한국신재생에너지학회:학술대회논문집
• /
• 2011.11a
• /
• pp.82.2-82.2
• /
• 2011

A channel design which is closely related with the mass transport overpotential is one of the most important procedures to optimize the whole fuel cell performance. In this study, three dimensional results of a numerical study for gas distribution in channels of a molten carbonate fuel cell (MCFC) unit cell for a 1kW class stack was presented. The relationship between the fuel and air distribution in the anode and cathode channels of the unit cell and the electric performance was observed. A charge balance model in the electrodes and the electrolyte coupled with a heat transfer model and a fluid flow model in the porous electrodes and the channels was solved for the mass, momentum, energy, species and charge conservation. The electronic and ionic charge balance in the anode and cathode current feeders, the electrolyte and GDEs were solved for using Ohm's law, while Butler-Volmer charge transfer kinetics described the charge transfer current density. The material transport was described by the diffusion and convection equations and Navier-Stokes equations govern the flow in the open channel. It was assumed that heat is produced by the electrochemical reactions and joule heating due to the electrical currents.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.19 no.6
• /
• pp.551-559
• /
• 2009

In this study, flow-induced vibration(FIV) analyses have been conducted for a 3D compressor blade model. Advanced computational analysis system based on computational fluid dynamics(CFD) and computational structural dynamics(CSD) has been developed in order to investigate detailed dynamic responses of designed compressor blades. Fluid domains are modeled using the computational grid system with local grid deforming and remeshing techniques. Reynolds-averaged Navier-Stokes equations with $\kappa-\epsilon$ turbulence model are solved for unsteady flow problems of the rotating compressor model. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of the 3D compressor blade for fluid-structure interaction(FSI) problems. Detailed dynamic responses and instantaneous pressure contours on the blade surfaces considering flow-separation effects are presented to show the multi-physical phenomenon of the rotating compressor blade.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.1
• /
• pp.104-114
• /
• 1993

In this paper, effects of free surface on viscous flow is investigated. Continuity equation coupled with Navier-Stokes equations are solved numerically by using an artificial compressibility method[1, 2]. The body-fitted generalized curvilinear coordinate system is employed to deal with arbitrary body shape. The IAF scheme with finite difference method is used to solve the equations, and a diagonal algorithm is applied to time-varying Jacobian matrices for the computational economics. Free surface shape is obtained by applying zero pressure condition to still water surface at each time step. A numerical test is made for larminar flow around a circular cylinder vertically piercing the free surface. Computed flow patterns are largely affected by the existance of free surface in low Reynolds number flows treated in this paper. Free surface causes viscous pressure drag to vary much in depth direction in accordance with the variations of flow pattern.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.2
• /
• pp.835-850
• /
• 2019

The numerical simulations for the Vortex-induced Vibration (VIV) of the cylinders with different combinations of mass ratio and frequency ratio were performed under the Reynolds (Re) number ranges of 1450-10200, 5800-40800 and 13050-91800 by using the embedded programs in OpenFoam. By combining with the modified SST k-ω turbulence model, the coupled Unsteady Reynolds-Average Navier-Stokes equations and double-degree-of-freedom vibration equations were solved. After analyzing the results, it is found that the some characteristics of the VIV have changed with the increase of the range of Re number, and the effects of Re number on vibration characteristics are also different under different combinations of mass ratio and frequency ratio. On this basis, the influence law of Re number on the characteristics of VIV of the cylinders is summarized, which can provide a reference for the research of VIV under higher Re number.

• Coupled systems mechanics
• /
• v.9 no.3
• /
• pp.289-309
• /
• 2020

The paper studies the fluid flow profile contained between the orthotropic plate and rigid wall under the action of the moving load on the plate and main attention is focused on the fluid velocity profile in the load moving direction. It is assumed that the plate material is orthotropic one and the fluid is viscous and barotropic compressible. The plane-strain state in the plate and the plane flow of the fluid is considered. The motion of the plate is described by utilizing the exact equations of elastodynamics for anisotropic bodies, however, the flow of the fluid by utilizing the linearized Navier-Stokes equations. For the solution of the corresponding boundary value problem, the moving coordinate system associated with the moving load is introduced, after which the exponential Fourier transformation is employed with respect to the coordinate which indicates the distance of the material points from the moving load. The exact analytical expressions for the Fourier transforms of the sought values are obtained, the originals of which are determined numerically. Presented numerical results and their analyses are focused on the question of how the moving load acting on the face plane of the plate which is not in the contact with the fluid can cause the fluid flow and what type profile has this flow along the thickness direction of the strip filled by the fluid and, finally, how this profile changes ahead and behind with the distance of the moving load.

• Wind and Structures
• /
• v.16 no.5
• /
• pp.411-431
• /
• 2013

A computational study of vortex-induced transverse vibrations of a cylinder with low mass-damping is presented. An Arbitrary Lagrangian-Eulerian (ALE) formulation of the Unsteady Reynolds-Averaged Navier-Stokes equations (URANS), along with the Spalart-Allmaras (SA) one-equation turbulence model, are coupled conservatively with rigid body motion equations of the cylinder mounted on elastic supports in order to study the amplitude and frequency response of a freely vibrating cylinder, its flow-induced motion, Vortex Street, near-wake flow structure, and unsteady loading in a moderate range of Reynolds numbers. The time accurate response of the cylinder from rest to its limit cycle is studied to explore the effects of Reynolds number on the start of large displacements, motion amplitude, and frequency. The computational results are compared with published physical experiments and numerical studies. The maximum amplitudes of displacements computed for various Reynolds numbers are smaller than the experimental values; however, the overall agreement of the results is quite satisfactory, and the upper branch of the limit-cycle displacement amplitude vs. reduced velocity response is captured, a feature that was missed by other studies. Vortex shedding modes, lock-in phenomena, frequency response, and phase angles are also in agreement with experiments.

• Journal of the Korean Society for Aviation and Aeronautics
• /
• v.17 no.2
• /
• pp.8-17
• /
• 2009

In this study, flutter analyses for supercritical airfoil have been conducted in transonic region. Advanced computational analysis system based on computational fluid dynamics (CFD) and computational structural dynamics (CSD) has been developed in order to investigate detailed static and dynamic responses of supercritical airfoil. Reynolds-averaged Navier-Stokes equations with Spalart-Allmaras (S-A) and SST ${\kappa}-{\omega}$ turbulence models are solved for unsteady flow problems. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of cascades for fluid-structure interaction (FSI) problems. Also, flow-induced vibration (FIV) analyses for various supercritical airfoil models have been conducted. Detailed flutter responses for supercritical are presented to show the physical performance and vibration characteristics in various angle of attack.

• Particle and aerosol research
• /
• v.15 no.1
• /
• pp.27-36
• /
• 2019

Numerical analysis has been conducted on three-dimensional airflow distribution in the passenger cabin of a high-speed electric train. The types of air distribution systems investigated in the present study were those of TGV and Shinkansen. The Reynolds-averaged Navier-Stokes equations governing the mass and momentum conservations of the airflow in the cabin were solved by using a finite volume method, which are coupled with the standard $k-{\varepsilon}$ turbulence model equations. Predicted velocity distributions were presented on several selected planes in the passenger cabin. The present three-dimensional simulations were found to show the overall features of the airflow in the passenger cabin fairly well. In particular, it was shown that the type of air distribution for Shinkansen was more suitable for a non-smoking cabin than that for TGV.