NUMERICAL CALCULATION OF TURBULENT OPEN-CHANNEL FLOWS BY USING A MODIFIED /g-e TURBULENCE MODEL By Iehisa NEZU and Hiroji NAKAGA WA Numerical calculat

Similar documents
<4D F736F F D208F4390B38DC58F49938A8D6595A CA90858D48985F95B F8F43959C82B382EA82BD B5F2E646F6378>

Fig. 1 Experimental apparatus.

A Higher Weissenberg Number Analysis of Die-swell Flow of Viscoelastic Fluids Using a Decoupled Finite Element Method Iwata, Shuichi * 1/Aragaki, Tsut

Study of the "Vortex of Naruto" through multilevel remote sensing. Abstract Hydrodynamic characteristics of the "Vortex of Naruto" were investigated b

teionkogaku43_527

$arrow$ $\yen$ T (Yasutala Nagano) $arrow$ $\yen$ ?,,?,., (1),, (, ).,, $\langle$2),, (3),.., (4),,,., CFD ( ),,., CFD,.,,,

IHIMU Energy-Saving Principle of the IHIMU Semicircular Duct and Its Application to the Flow Field Around Full Scale Ships IHI GHG IHIMU CFD PIV IHI M

Developement of Plastic Collocation Method Extension of Plastic Node Method by Yukio Ueda, Member Masahiko Fujikubo, Member Masahiro Miura, Member Sum

Estimation of Photovoltaic Module Temperature Rise Motonobu Yukawa, Member, Masahisa Asaoka, Non-member (Mitsubishi Electric Corp.) Keigi Takahara, Me

Optical Lenses CCD Camera Laser Sheet Wind Turbine with med Diffuser Pitot Tube PC Fig.1 Experimental facility. Transparent Diffuser Double Pulsed Nd:

Natural Convection Heat Transfer in a Horizontal Porous Enclosure with High Porosity Yasuaki SHIINA*4, Kota ISHIKAWA and Makoto HISHIDA Nuclear Applie

The Effect of the Circumferential Temperature Change on the Change in the Strain Energy of Carbon Steel during the Rotatory Bending Fatigue Test by Ch

Study on Application of the cos a Method to Neutron Stress Measurement Toshihiko SASAKI*3 and Yukio HIROSE Department of Materials Science and Enginee

Study on Throw Accuracy for Baseball Pitching Machine with Roller (Study of Seam of Ball and Roller) Shinobu SAKAI*5, Juhachi ODA, Kengo KAWATA and Yu

bosai-2002.dvi

7 OpenFOAM 6) OpenFOAM (Fujitsu PRIMERGY BX9, TFLOPS) Fluent 8) ( ) 9, 1) 11 13) OpenFOAM - realizable k-ε 1) Launder-Gibson 15) OpenFOAM 1.6 CFD ( )

Mikio Yamamoto: Dynamical Measurement of the E-effect in Iron-Cobalt Alloys. The AE-effect (change in Young's modulus of elasticity with magnetization

空力騒音シミュレータの開発

日立金属技報 Vol.34

塗装深み感の要因解析

The Evaluation of LBB Behavior and Crack Opening Displacement on Statically Indeterminate Piping System Subjected to Monotonic Load The plastic collap

Flow Around a Circular Cylinder with Tangential Blowing near a Plane Boundary (2nd Report, A Study on Unsteady Characteristics) Shimpei OKAYASU, Kotar

20 $P_{S}=v_{0}\tau_{0}/r_{0}$ (3) $v_{0}$ $r_{0}$ $l(r)$ $l(r)=p_{s}r$ $[3 $ $1+P_{s}$ $P_{s}\ll 1$ $P_{s}\gg 1$ ( ) $P_{s}$ ( ) 2 (2) (2) $t=0$ $P(t

チャネル乱流における流体線の伸長

A Study on Throw Simulation for Baseball Pitching Machine with Rollers and Its Optimization Shinobu SAKAI*5, Yuichiro KITAGAWA, Ryo KANAI and Juhachi

Tetsu-to-Hagane Vol. 87 (2001) No. 5 Table 1. Physical properties of particles. (a) side view (b) front view Fig. 1. Experimental apparatus with semic

Tetsu-to-Hagane Vol. 87 (2001) No. 1 Fig. 1. Experimental apparatus for PIV measurement. Fig. 2. Detail of A.


Fig. 2 Effect of oxygen partial pressure on interfacial tensions between molten copper and fayalite slag (Fe/Si0 2=1.23) at 1473 K. Fig. s Effect or o

(Mamoru Tanahashi) Department of Mechanical and Aerospaoe Engineering Tokyo Institute of Technology ,,., ,, $\sim$,,

$\hat{\grave{\grave{\lambda}}}$ $\grave{\neg}\backslash \backslash ^{}4$ $\approx \mathrm{t}\triangleleft\wedge$ $10^{4}$ $10^{\backslash }$ $4^{\math

Fig. 1. Schematic drawing of testing system. 71 ( 1 )

Visual Evaluation of Polka-dot Patterns Yoojin LEE and Nobuko NARUSE * Granduate School of Bunka Women's University, and * Faculty of Fashion Science,

Fig. 1 Flow diagram of experimental apparatus employed Fig. 2 Porosity change during sulfurization of reduced sample pellets

Fig. 1 Hydrostatic Thrust Bearing Fig. 2 Point loading of elastic half-space

Table 1. Reluctance equalization design. Fig. 2. Voltage vector of LSynRM. Fig. 4. Analytical model. Table 2. Specifications of analytical models. Fig

Table 1. Assumed performance of a water electrol ysis plant. Fig. 1. Structure of a proposed power generation system utilizing waste heat from factori

Rate of Oxidation of Liquid Iron by Pure Oxygen Shiro BAN-YA and Jae-Dong SHIM Synopsis: The rate of oxidation of liquid iron by oxygen gas has been s

第62巻 第1号 平成24年4月/石こうを用いた木材ペレット

1) T. L. Cottrel, A. J. Matheson, Trans. Farad. Soc., 58, 2336(1962). 2) E. N. Chesnokov, V. N. Panfilov, Teor. Eksp. Khimiya, 17, 699(1981). 3) M. Ko


Vol. 51 No (2000) Thermo-Physiological Responses of the Foot under C Thermal Conditions Fusako IWASAKI, Yuri NANAMEKI,* Tomoko KOSHIB

The Evaluation on Impact Strength of Structural Elements by Means of Drop Weight Test Elastic Response and Elastic Limit by Hiroshi Maenaka, Member Sh

T05_Nd-Fe-B磁石.indd

Fig. 4. Configuration of fatigue test specimen. Table I. Mechanical property of test materials. Table II. Full scale fatigue test conditions and test

Studies of Foot Form for Footwear Design (Part 9) : Characteristics of the Foot Form of Young and Elder Women Based on their Sizes of Ball Joint Girth

133 1.,,, [1] [2],,,,, $[3],[4]$,,,,,,,,, [5] [6],,,,,, [7], interface,,,, Navier-Stokes, $Petr\dot{o}$v-Galerkin [8], $(,)$ $()$,,

Field Measurement of Aeolian Sand Flux over a Sand Dune Slope Junaidi Keiko UDO, Shota MITSUSHIO, JUNAIDI, Shin-ichi AOKI, Shigeru KATO and Akira MANO

Corrections of the Results of Airborne Monitoring Surveys by MEXT and Ibaraki Prefecture


1. Precise Determination of BaAl2O4 Cell and Certification of the Formation of Iron Bearing Solid Solution. By Hiroshi UCHIKAWA and Koichi TSUKIYAMA (

EQUIVALENT TRANSFORMATION TECHNIQUE FOR ISLANDING DETECTION METHODS OF SYNCHRONOUS GENERATOR -REACTIVE POWER PERTURBATION METHODS USING AVR OR SVC- Ju

Fig. 1 KAMOME50-2 Table 1 Principal dimensions Fig.2 Configuration of the hydrofoils (Endurance and sprint foil) Fig. 3 Schematic view of the vortex l

Research Reports on Information Science and Electrical Engineering of Kyushu University Vol.11, No.1, March 2006 Numerical Analysis of Scattering Atom

TM

Spacecraft Propulsion Using Solar Energy Spacecraft with Magnetic Field Light from the Sun Solar Wind Thrust Mirror Solar Sail Thrust production by li

A Nutritional Study of Anemia in Pregnancy Hematologic Characteristics in Pregnancy (Part 1) Keizo Shiraki, Fumiko Hisaoka Department of Nutrition, Sc


3. ( 1 ) Linear Congruential Generator:LCG 6) (Mersenne Twister:MT ), L 1 ( 2 ) 4 4 G (i,j) < G > < G 2 > < G > 2 g (ij) i= L j= N

315 * An Experimental Study on the Characteristic of Mean Flow in Supersonic Boundary Layer Transition Shoji SAKAUE, Department of Aerospace Engineeri

Vol.7 No.2 ( ) in mm m/s 40 m/s 20 m/s m/s 20m/s 1999 US ,

16_.....E...._.I.v2006


Quantitative Relationship between SAR and Temperature Rise inside Eyeball in a Realistic Human Head Model for 1.5 GHz-Microwave Exposure Kiyofumi Taka

Experimental and Clinical Studies of Pregnant Hypertension Takashi SHIMAZU Department of Obstetrics and Gynecology, Osaka City University Medical Scho

r z m ε r ε θ z rθ

Effect of Autofrettage on Fatigue Crack Growth Rate for High Pressure Reactors For high pressure reactor vessels such as for polyethylene production,

(43) Vol.33, No.6(1977) T-239 MUTUAL DIFFUSION AND CHANGE OF THE FINE STRUCTURE OF WET SPUN ANTI-PILLING ACRYLIC FIBER DURING COAGULATION, DRAWING AND

非線形長波モデルと流体粒子法による津波シミュレータの開発 I_ m ρ v p h g a b a 2h b r ab a b Fang W r ab h 5 Wendland 1995 q= r ab /h a d W r ab h

206“ƒŁ\”ƒ-fl_“H„¤‰ZŁñ

電子部品はんだ接合部の熱疲労寿命解析

H10Masuki

2006 [3] Scratch Squeak PEN [4] PenFlowchart 2 3 PenFlowchart 4 PenFlowchart PEN xdncl PEN [5] PEN xdncl DNCL 1 1 [6] 1 PEN Fig. 1 The PEN

Relationship between Men's Commuting Clothes and their Awareness toward Work Setsuko Yajiri*, Tomoko Takaoka**, Machiko Morita*** and Shigeo Kobayashi


: u i = (2) x i Smagorinsky τ ij τ [3] ij u i u j u i u j = 2ν SGS S ij, (3) ν SGS = (C s ) 2 S (4) x i a u i ρ p P T u ν τ ij S c ν SGS S csgs

Title 混合体モデルに基づく圧縮性流体と移動する固体の熱連成計算手法 Author(s) 鳥生, 大祐 ; 牛島, 省 Citation 土木学会論文集 A2( 応用力学 ) = Journal of Japan Civil Engineers, Ser. A2 (2017), 73 Issue

( ) [1] [4] ( ) 2. [5] [6] Piano Tutor[7] [1], [2], [8], [9] Radiobaton[10] Two Finger Piano[11] Coloring-in Piano[12] ism[13] MIDI MIDI 1 Fig. 1 Syst

Fig, 1. Waveform of the short-circuit current peculiar to a metal. Fig. 2. Waveform of arc short-circuit current. 398 T. IEE Japan, Vol. 113-B, No. 4,

1) K. J. Laidler, "Reaction Kinetics", Vol. II, Pergamon Press, New York (1963) Chap. 1 ; P. G. Ashmore, "Catalysis and Inhibition of Chemical Reactio

Study on Heat Transfer and Flow Characteristics of a Bank of Tubes with Tube-to-Baffle Leakage Z 29 < Z < 55 R 0.4 < R < 0.75 RANS This study investig

食品工学.indb

Computer Simulation in Thermoplastic Injection Molding Takaaki Matsuoka Toyota Central Research and Development Laboratories, Inc. 41-1, Yokomichi, Na

レーザ誘起蛍光法( LIF法) によるピストンの油膜挙動の解析



Fig. 1 Schematic construction of a PWS vehicle Fig. 2 Main power circuit of an inverter system for two motors drive

三校永谷.indd

Journal of Geography 116 (6) Configuration of Rapid Digital Mapping System Using Tablet PC and its Application to Obtaining Ground Truth

Fig. 3 Flow diagram of image processing. Black rectangle in the photo indicates the processing area (128 x 32 pixels).

(1) 2

Instability of Aerostatic Journal Bearings with Porous Floating Bush at High Speeds Masaaki MIYATAKE *4, Shigeka YOSHIMOTO, Tomoaki CHIBA and Akira CH


3_23.dvi

2. 実 験 方 法 第 1 表 N min 1 Ns 0.15 Ns = N / 60 Q 1/2 /Had 3/4 Q ( m 3 /s ) Had ( J/kg ) 第 1 図 VLD 5 mm 2 R 2 R 2 VLD 1.38 R R 2 10 mm 1.14

DVIOUT-ajhe

Jan THE JAPANESE JOURNAL OF ANTIBIOTICS XL-1 Table 1. Outline of administering doses, routes and sampling times *: 4 ml/hr/kg Bacillus subtilis

大学における原価計算教育の現状と課題

UDC : ' : '24' : '24'26' : : A Study of Condition of Pits Formation and Their Fe

untitled

Transcription:

NUMERICAL CALCULATION OF TURBULENT OPEN-CHANNEL FLOWS BY USING A MODIFIED /g-e TURBULENCE MODEL By Iehisa NEZU and Hiroji NAKAGA WA Numerical calculation techniques of turbulent shear flows are classified into two categories:one is k-c turbulence model and the other is large eddy simulation (LES). The standard k-s model has been established at present to predict turbulent structure in closed duct flows, while LES is being developed to predict coherent eddy structure in more simple duct flows. The standard /g-e model cannot be, however, applied to open channel surface flows, because the turbulence near the free surface is more depressed than the closed duct flows. In the present study a new modified k-e model is proposed to predict reasonably the turbulent structure in open channel flows with both of low and high Reynolds numbers. The numerical calculations indicate a splendid agreement with the experimental data which were obtained by making use of hot-film and Laser Doppler anemometers. Keywords:open-channel turbulence, numerical calculation, free surface, k-e model, turbulence model

P=p9(h-y)cose+P' (4 ) (5) vt=cu.2e (6) (7)

Up1 U-*xInyp+A (16) t(-uv)ady )=1 wvt (23) (a) Experimental Values Fig. 1 Effect of free-surface damping factor on eddy viscosity,

Fig. 2 Effect of free-surface damping factor on mean velocity distribution. Fig, 3 Effect of free-surface damping factor on turbulent energy k/u. Fig. 4 Effect of free-surface damping factor on energy dissipation Eh/ U*.

Table 1 Hydraulic conditions for numerical calculations. (a) High Reynolds-number Series (b) Low Reynolds-number Series

Fig. 5 Numerical results of mean velocity distributions at high Reynolds numbers (HR). Fig. 6 Turbulent energy Ic/ U* vs y/h at high Reynolds numbers. Fig. 7 Distribution of Reynolds stress -u/u. v* Fig. 8 Comparison of calculated values of energy dissipation with experimental data of Nezu (1977).

(a) Experimental Values (Nezu 1977) High Reynolds-number Series (HR) RUN Re Rsra O 1 2.0=104 1.34=10= E 2 5.0104 2.9210' + 3 1.0104 5.34=10' X U 5.0104 2.25=10' O 5 1.0104 x 0.2210' 4 5 1.0104 3.05=10" Fig. 9 Comparison of calculated values of generationdissipation relations with experimental data.

Low Reynolds-number Series (LR) RUN Re Rsn o 1 5.0=102 3.82101 0 2 1.0x103 6.52x101 + 3 2.0x103 1.22x10 X 1 5.0x103 2.74x102 O 5 1.0x104 5.12x10 + 8 5.0x104 2.22x10 R 7 1.0=10 11.17=10' Fig. 10 Numerical results of mean velocity distributions at low Reynolds numbers (LR). Low Reynolds-number Series (LR) 0 RUN Rt Rsl. 1 5.0.10' e 2 1.010' 63.8622.1100' + 3 2.0.10' 1.220 X 4 S.0.10' 2. 74.10' o 5 1.0.10' 5. 12.10' * 6 5.0.10' 2.22.10' R 7 1.0.10' 4. I Tl0' Low Reynolds-number Series (LR) RUN Rt Rsln O 1 5.0.10' 3.82.10' e 2 1.0.10' 6.62.10' + 3 2.0.10' 1.2210' X 4 5.010' 2.7410'.O 5 1.0' 10' 5.1210' 1 6 5.0=10' 2.2210'.0 7 1.010' 10 11.1710' Fig, 11 Calculated distribution of turbulent energy in the wall region. Fig. 12 Calculated distribution of Reynolds stress in the wall region. Low Reynolds-number Series (LR) RUN Rt Rs. O 1 5.0.10' 3.82.10' e 2 1.0.10' 6.6210' + 3 2.0.1 0' 1. 2 21 0' X 4 5.0.10' 2. 74.10' O 5 1.0.10' S. 12.10' t 6 5.0.10' 2. 22x10' R 7 1.0.10' 4. 17.10' lo. Reynolds-nuwber Series (LR) O RUN R=M. 5.0x10' 1.0xt0' 2 I 3.12.10' 6.62. n + 3 2.0.10' 1.22.10' 10' x r 5.x10' 2. 7Vx10' o * S 6 1.0.10' 5. ax 10' 5.12.10' 2.22.10' x 7 I.alo4 12.10' O : E. perimentl Vlues (taufer 1954) Fig. 14 Relation between turbulent generation and dissipation in the wall region. Fig. 13 Calculated distribution of dissipation in the wall region.

2) Rodi, W.: Turbulence Models and Their Application in Hydraulics, IAHR, Delft, 1980. 3) Naot, D. and Rodi, W.: Calculation of Secondary Currents in Channel Flow, J. Hydraulics Div., ASCE, Vol. 108, HY-8, pp. 948-968, 1982. 4) Nakagawa, H., Nezu, I. and Ueda, H.: Turbulence of Open Channel Flow over Smooth and Rough Beds, Proc. JSCE, No. 241, pp. 155--168, 1975. 6) Celik, I, and Rodi, W.: Simulation of Free-Surface Effects in Turbulent Channel Flows, Physico-Chemical Hydrodynamics, Vol. 5, pp. 217-227, 1984. 7) Jones, W. P. and Launder, B. E.: The Calculation of Low-Reynolds-Number Phenomena with a Two- Equation Model of Turbulence, Int. J. Heat and Mass 10) Patel, V. C., Rodi, W. and Scheuerer, G.: Turbulence Model for Near-Wall and Low Reynolds Number Flows A Review, AIAA J., pp. 1308-1319, 1985. 11) Nezu, I. and Rodi, W.: Open-Channel Flow Measurements with a Laser Doppler Anemometer, J. Hydraulic Eng., ASCE, Vol. 112, pp. 335-355, 1986. 12) van Driest, E. R.: On Turbulent Flow Near a Wa!!, J. Aeron. Sci., Vol. 23, pp. 1007-1011, 1956. 13) Hussain, A. K. M. F, and Reynolds, W. C.: Measurements in Fully Developed Turbulent Channel Flow, J. Fluids Eng., ASME, Vol. 97, pp. 568-580, 1975. 15) Gosman, A. D., Launder, B. E, and Reece, G. J.: Computer-Aided Engineering; Heat Transfer and Fluid Flow, John Wiley & Sons, 1985. 16) Eckelmann, H.: The Structure of the Viscous Sublayer and the Adjacent Wall Region in a Turbulent Channel Flow, J. Fluid Mech., Vol. 65, pp. 439-459, 1974. 17) Laufer, J.: The Structure of Turbulence in Fully Developed Pipe Flow, NACA, TR-1174, 1954. 19) Rodi, W.: Use of Advanced Turbulence Models for Calculating the Flow and Pollutant Spreading in Rivers, 3rd Int. Symp. River Sedim., Univ. of Mississippi, pp. 1369-1382, 1986.