METHOD OF WATER QUALITY IMPROVEMENT BY USE OF ICE FLOE CONTROL STRUCTURE
|
|
- Thomas Quinn
- 5 years ago
- Views:
Transcription
1 METHOD OF WATER QUALITY IMPROVEMENT BY USE OF ICE FLOE CONTROL STRUCTURE Shunsuke Makita 1 and Hiroshi Saeki 1 Graduate School of Engineering, Hokkaido Universit, Sapporo, Japan Hokkaido Universit, Sapporo, Japan ABSTRACT In Hokkaido, Japan, the aquaculture has been performed in man semi-closed water areas. Ice control structures to prevent the invasion of ice floes had started to be installed at inlets of some of those areas. However, these structures have been used onl during winter. In this paper, the effective use of these structures aiming the improvement of the water qualit of these areas is proposed. B numerical analses, the became clear that it is possible to suppl the seawater to the center part of the area b use of this structure as a training dike at the time of inflow, and in the case that the water area has two or more inlets, it is possible to transfer the seawater b producing the water current in an arbitrar direction b use of this structure as a gorge at one inlet at the time of outflow. KEY WORDS: Ice control structure; Ice boom; semi-closed water area; water echange; training. INTRODUCTION In Hokkaido, more than 3 semi-closed water areas eist, and aquaculture has been performed in man of those water areas. Since semi-closed water areas connected to open sea b a narrow inlet are ver calm water area, such areas are suitable for the aquaculture of scallops, osters, and other marine products. However there is also the risk that the water qualit ma suddenl decrease with change of the surrounding environment. For permanent use for aquaculture, maintenance of the water qualit is needed. As one of the water qualit preservation methods, promotion of water echange has been pointed out. Man of these water areas are located on the coast of the Okhotsk Sea. Ice floes normall reach the coast in winter; tides cause ice floes to pass through inlets into the water area. These ice floes have caused serious damage, such as destruction and wear, to aquaculture facilities. To prevent such damage, ice control structures had started to be installed at inlets of some of those areas. These structures are composed b some piles attached to the seabed and some ice booms installed between these piles. However, ice booms are removed and these structures have not been used for periods other than winter. In this paper, the effective use of these structures aiming the improvement of the water qualit of these areas is proposed
2 METHOD OF ANALYSIS Model of numerical simulation For numerical simulation, a multi-level-model was used (Fujihara, 199). In this model, the object is divided b mesh for horizontal direction, and is divided b laer for vertical direction. The variation of the thickness of laer is considered onl for the surface laer, and the vertical flow components are considered for underling laers. Therefore this model is suitable for analsis of flow conditions in a semi-closed water area where vertical flow such as the lifting of nutrients b upwelling is important. As basic equation, the continuit equation is shown in equation (1) and the momentum equations are shown in equation () and equation (3). divu u v w = + + z (1) M k ζ = ( M kuk) ( M kvk) ( uw) + ( uw) gh k k+ 1 k t Mk Mk 1 k 1, k k, k N + N + ( τ τ ) τ f ρ ρ Nk ζ = ( Nkuk) ( Nkvk) ( vw) + ( vw) gh k k+ 1 k t Nk Nk 1 k 1, k k, k N + N + ( τ τ ) τ f ρ ρ () (3) M k, Nk ; Linear quantit, N, N ; Coefficient of horizontal vorte viscosit, ς ; Water level, hk ; Thickness of laer, ρ ; Water densit, k ; Number of laer, τ, τ ; Frictional force between laers calculated from following equations τ = ργ u u + v, τ = ργ v u + v (4) γ ; Coefficient of friction (Coefficient of friction at water surface; bottom; γ b, Coefficient of internal friction of water; i γ a, Coefficient of friction at γ ), u,v ; Water velocit The terms shown in underlined part of the momentum equations are terms of resistance forces of aquaculture facilities installed in water areas. And these resistance forces are calculated b equation (5). Coefficients of resistance are set based on two coefficients of resistance of aquaculture facilities obtained from the eperiment b authors (Makita, 4). 1 1 τ = ρcu u + v, τ = ρcv u + v (5) C ; Coefficient of resistance -138-
3 Method of evaluation of water echange The water echange was evaluated b performing mass transfer simulation for a diffusion mass initiall diffused to a uniform concentration in the semi-closed water area. Diffusion mass transfer was calculated b the adjective diffusion equation (6). The initial concentration of diffusion mass was 1. in the semi-closed water area, and. in the open sea. Therefore the area in which the concentration of a diffusion mass decreases shows the area into which the water of open sea flowed, and it is possible to evaluate the amount of outflow of water of a lagoon b the residual rate of a diffusion mass. SS SS SS SS = ( uss ) ( vss ) ( wss ) + K + K + K z t z z z (6) SS ; Concentration of diffusion mass, K, K, K ; Coefficient of horizontal diffusion Coefficients Coefficient of friction at water surface; Although the coefficient varies complel according to wind, flow, etc., the following constant is often used. 3 γ a = 1.3 Coefficient of friction at bottom; In the multi-level-model, the following constant obtained b Hasen through tidal current observation is often used. This value was used in this calculation. 3 γ b =.6 Coefficient of internal friction of water; The value that is 1/~3/4 of the coefficient of friction at bottom is commonl used. The following value was used in this calculation. 3 γ i =.5 Coefficient of horizontal diffusion; It is known that this coefficient is proportional to the scale of motion raised to the 4/3 powers, and calculated using the following equation (7). z 4 3 K = K = cl (7) c ; Constant (about.1) and L is the grid size The grid size of this simulation was m, and the following value obtained using equation (7) was used in this calculation. K = K = 4.64 (cm sec) Coefficient of vertical diffusion mass; Although the value calculated from Richardson number, is proposed, in the multi-level-model, the following constant is often used. K = z 1. (cm /sec) -139-
4 Coefficient of horizontal vorte viscosit; Although this coefficient varies according to the state and scale of motion, a value nearl equal to the coefficient of horizontal diffusion is commonl used. The value equal to the coefficient of horizontal diffusion was used in this calculation. N = N = 4.64 (cm sec) Coefficient of resistance of aquaculture facilities; Coefficients of resistance in the equation (5) are set Table 1 Coefficients of resistance based on two kinds of coefficients of resistance of 3 aquaculture facilities, the coefficient in the state where a Parallel C 1 =. flow acts parallel to the main line of the long-line formula cultivating facilit, and the coefficient in the Perpendicular C =. state where a flow acts perpendicular to it. Coefficients of resistance of aquaculture facilities are shown in Table 1. In the state that the flow is diagonal to the main wire, resistance force is calculated from these tow coefficients and the action angle of the flow based on the eperimental result b authors (Makita, 4). Object of simulation The object of this simulation is the A D Saroma Lagoon, which is a tpical semi-closed water area in Hokkaido. The B Inlet No. 1 Inlet No. C lagoon is located on the Coast of the Okhotsk Sea and is used intensivel for aquaculture of shellfish, such as scallops and osters. This lagoon is connected to the open sea b two inlets. Ice control structures have been installed at both km inlets. In this simulation, the control of the flow was made onl b the ice Figure 1 Range of the simulation control structure at the inlet No. 1. It is composed b 14 piles and 13 ice booms installed between them. The range of simulation was the Saroma Lagoon and the area 5km offshore from inlet No. 1. Figure 1 shows the range of the simulation, and the simulation was made b dividing the range into a m-mesh. As the boundar condition, Sommerfeld s open conditions were set for AB, CD and DA, as accesses leading to the open sea. Variations in eternal tide for a 4-hour period and a semi-tide range of.4m were given onl to boundar AD. In the figure, the area of hatching shows the area in which the aquaculture facilities are installed, and the direction of the line shows the direction of main line of facilities. Reliabilit of simulation In the result of the simulation, the vorte flow occurred in a position 3m from inlet No. at the time of inflow. In the local observation, it is reported that the vorte of separation occurs in the almost the same position. This can confirm that the setup of coefficients, such as the coefficient of -14-
5 vorte viscosit, is appropriate. Figure shows the vector diagram of the flow around the inlet No.. The calculated value of the flow velocit at the lagoon inlet is equivalent to the observed value. Figure 3 shows velocities at the inlet. 5 Vorte flow 3m Figure Vector diagram around inlet No. Velocit (cm/sec ) Difference of water level (m) Figure 3 Velocit at inlet Calculated value Observed value Flow control b use of ice control structure Following matters became clear from results of the analsis b authors (Makita, ). At the semi-closed water area, inflows from inlets are linear jet flows and seawater does not easil supplied to the center part of the lagoon. When the amount of outflow at one inlet decreases due to a gorge, such as an accumulation of ice floes, the amount of outflow at another inlet increases, and the water current in a direction from one inlet to another inlet is produced. In this analsis, following possibilities were eamined. The possibilit to suppl the seawater to the center part of the lagoon b use of the ice control structure as a training dike at the time of inflow, and the possibilit to transfer the seawater to the center part b producing the water current in an arbitrar direction b use of this structure as a gorge at one inlet at the time of outflow. Figure 4 shows the location of the ice control structure on the mesh. In addition, m-mesh is Inlet No. 1 nested at the gra area. 13 spans of ice boom were divided into the section A, form the 1 st to the 6 th span, and the section B, from the 7 th to the 13 th span. The calculation processing in each section is summarized in the Table. Pier Section A Section B Ice Booms Figure 4 Location of ice control structure Table Calculation processing Section A Obstruct the flow at both of the time of inflow and the time of outflow. Section B Act the resistance to the flow onl at the time of outflow
6 RESULT OF SIMULATION AND DISCUSSION Suppl of seawater The simulation was made in the state that the flow is and in the state that the flow is. The vector diagrams of the flow, at the time when the velocit of the inflow is the maimum, are shown in Figure 5 and Figure 6. It can be confirmed that the inflow from the inlet No. 1 is trained to the center part of the lagoon in the state that the flow is. The isoconcentration maps of the diffusion mass, after 1 ccle of tide, are shown in Figure 7 and Figure 8. In the state that the flow is, the area where the concentration of the diffusion mass is 1., that is, the area to which the seawater is not supplied, remains at the center part of the lagoon. In the state that the flow is, the area to which the seawater is supplied is reduced but the area where the concentration of the diffusion mass is 1. dose not remain. Figure 5 Vector diagram - Figure 6 Vector diagram Figure 7 Isoconcentration map - Figure 8 Isoconcentration map - The maimum flow velocit in each state is shown in Figure 9. In the figure, the velociti is shown b a positive value in the case of the inflow from the open sea to the lagoon, and is shown b a negative value in the case of outflow. In the state that the flow is, the velocit at the inlet No. 1 decreases and the velocit at the inlet No. increases. Velocit (m/sec) inflow outflow Inlet No. 1 Inlet No. Figure 9 Velocit at each inlet -14-
7 Production the water current The rate of quantit shown in Figure, epresses the role of each inlet in the water echange. The inlet No. 1 plas the role of inflow and the inlet No. plas the role of outflow even in the state that the flow is, both of the inlet pla each role more greatl in the state that the flow is. It is guessed that the amount the water that flows in through the inlet No. 1 and flows out through the inlet No. increases, and the water current in direction form inlet No. 1 to the inlet No. is promoted. Figure 11 shows the residual rate of the diffusion mass. The rate is.87 in the state that the flow is, and.89 in the state the flow is. It is confirmed that the amount of water echanged during 1 ccle of tide decreases. Rate of quantit Residual rate Rate of quantit= 1.5 Inlet No. 1 Inlet No. ontrolled ( quantit of inflow ) + ( quantit of outflow) ( average of the absolute value of quantit) Figure Rate of quantit Figure 11 Residual rate Water echange of ccles period of tide So far, the water echange performed during 1 ccle of tide has been eamined. Henceforth, the water echange performed during ccles is eamined. The isoconcentration maps of the diffusion mass are shown in Figure 1 and Figure 13. In the state that the flow is, the area to which the seawater is supplied is formed in the vicinit of the inlet No. 1. In the state that the flow is, it is formed in the vicinit of the center of the lagoon, and the interval of the contour line is wider in the side of the center of the lagoon than in the side of the inlet No. 1. It is guessed that the seawater is transferred b the water current in direction form inlet No. 1 to the inlet No.... Figure 1 Isoconcentration map - Figure 13 Isoconcentration map
8 The change of the amount of diffusion mass during ccles is shown in Figure 14. It is confirm that the residual rate of the diffusion mass, for ccles of tide, is higher in the state the flow is than that in the state the flow is and the difference of the rate increase at the ccle. Residual rate Ccle of tide From the above-mentioned, it is confirmed that the seawater can be Figure 14 Residual rate supplied to the center part of the lagoon b the flow control using the ice control structure, but the amount of the water echange decrease. Thus, it can be said that the utilization plan, such as limitation at the period of the instillation of the structure, is needed. CONCLUSIONS In conclusion, the possibilit of improvement of the water qualit of the semi-closed water area is summarized as follows. 1) It is possible to suppl the seawater to the center part of the area b use of the ice control structure as a training dike at the time of inflow. ) It is possible to transfer the seawater b producing the water current in an arbitrar direction b use of Ice control structure as a gorge at one inlet at the time of outflow, in the case that the water area has two or more inlets. 3) The total amount of the water echange decrease b the inference of the flow control, the utilization plan, such as limitation at the period of the instillation of the structure, is needed. REFERENCES Fujihara, M., Akada, S. and Takeuchi, T. (199), Development of multi-level densit flow model and its application to the upwelling generated b artificial structures, Tech. rept N.R.I.F.E aquacul. & fish. port 14, pp13-35 Makita, S. and Saeki, H. (4) Influence of the aquaculture facilit to the water echange at the semi-closed water area, Annual Journal of Civil Engineering in the Ocean, JSCE, vol., pp Makita, S., Saeki, H. and Furua, A. (), Effect of ice on water flow at Saroma lagoon, Proceedings of 16 th International Smposium on Ice (IAHR ), Dunedin, New Zealand, pp
EFFECT OF ICE ON WATER FLOW AT SALOMA LAGOON
Ice in the Environment: Proceedings of the 16th IAHR International Symposium on Ice Dunedin, New Zealand, 2nd 6th December 22 International Association of Hydraulic Engineering and Research EFFECT OF ICE
More informationESCI 485 Air/sea Interaction Lesson 5 Oceanic Boundary Layer
ESCI 485 Air/sea Interaction Lesson 5 Oceanic Boundar Laer References: Descriptive Phsical Oceanograph, Pickard and Emer Introductor Dnamical Oceanograph, Pond and Pickard Principles of Ocean Phsics, Apel
More informationTHE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY. V. A. Sharifulin.
THE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY 1. Introduction V. A. Sharifulin Perm State Technical Universit, Perm, Russia e-mail: sharifulin@perm.ru Water
More informationBubble curtain system for prevention of seawater intrusion in coastal aquifers
Available online at www.pelagiaresearchlibrar.com Advances in Applied Science Research, 013, 4(4):81-90 ISSN: 0976-8610 CODEN (USA): AASRFC Bubble curtain sstem for prevention of seawater intrusion in
More information7. TURBULENCE SPRING 2019
7. TRBLENCE SPRING 2019 7.1 What is turbulence? 7.2 Momentum transfer in laminar and turbulent flow 7.3 Turbulence notation 7.4 Effect of turbulence on the mean flow 7.5 Turbulence generation and transport
More informationApply mass and momentum conservation to a differential control volume. Simple classical solutions of NS equations
Module 5: Navier-Stokes Equations: Appl mass and momentum conservation to a differential control volume Derive the Navier-Stokes equations Simple classical solutions of NS equations Use dimensional analsis
More informationWater Flow in Open Channels
The Islamic Universit of Gaza Facult of Engineering Civil Engineering Department Hdraulics - ECIV 33 Chapter 6 Water Flow in Open Channels Introduction An open channel is a duct in which the liquid flows
More informationConservation of Linear Momentum for a Differential Control Volume
Conservation of Linear Momentum for a Differential Control Volume When we applied the rate-form of the conservation of mass equation to a differential control volume (open sstem in Cartesian coordinates,
More informationSTUDY ON LOCAL SCOURING AT SURAMADU BRIDGE PIERS FOR STRUCTURAL INTEGRITY MONITORING
Proceedings of the 7 th International Conference on Asian and Pacific Coasts (APAC 03) ali, Indonesia, September 4-6, 03 STUDY ON LOCAL SCOURING AT SURAMADU BRIDGE PIERS FOR STRUCTURAL INTEGRITY MONITORING
More informationPrediction of changes in tidal system and deltas at Nakdong estuary due to construction of Busan new port
Prediction of changes in tidal system and deltas at Nakdong estuary due to construction of Busan new port H. Gm1 & G.-Y. park2 l Department of Civil & Environmental Engineering, Kookmin University, Korea
More informationFluid Mechanics II. Newton s second law applied to a control volume
Fluid Mechanics II Stead flow momentum equation Newton s second law applied to a control volume Fluids, either in a static or dnamic motion state, impose forces on immersed bodies and confining boundaries.
More informationME 321: FLUID MECHANICS-I
6/07/08 ME 3: LUID MECHANI-I Dr. A.B.M. Toufique Hasan Professor Department of Mechanical Engineering Bangladesh Universit of Engineering & Technolog (BUET), Dhaka Lecture- 4/07/08 Momentum Principle teacher.buet.ac.bd/toufiquehasan/
More informationChapter 4 Transport of Pollutants
4- Introduction Phs. 645: Environmental Phsics Phsics Department Yarmouk Universit hapter 4 Transport of Pollutants - e cannot avoid the production of pollutants. hat can we do? - Transform pollutants
More informationSIGNIFICANCE OF TIDAL CHANGE ON ABRASION AREA OF STRUCTURES DUE TO SEA ICE MOVEMENT
SIGNIFICANCE OF TIDAL CHANGE ON ABRASION AREA OF STRUCTURES DUE TO SEA ICE MOVEMENT Takahiro Takeuchi, Shinji Kioka and Hiroshi Saeki 3 Hachinohe Institute of Technology, Hachinohe, Japan Civil Engineering
More informationNUMERICAL SIMULATION AND ANALYSIS OF INTERFERENCE EFFECTS ON THE PORT CONTAINER GROUP
NUMERICAL SIMULATION AND ANALYSIS OF INTERFERENCE EFFECTS ON THE PORT CONTAINER GROUP Chun-hui Zhang, Wei Sun, Xing-qian Peng Graduate student, College of Civil Engineering, Huaqiao Universit, Quanzhou,
More informationPERSPECTIVE OF COASTAL VEGETATION PATCHES WITH TOPOGRAPHY VARIATIONS FOR TSUNAMI PROTECTION IN 2D - NUMERICAL MODELING
PERSPECTIVE OF COASTAL VEGETATION PATCHES WITH TOPOGRAPHY VARIATIONS FOR TSUNAMI PROTECTION IN D - NUMERICAL MODELING N. A. K. NANDASENA ), Norio TANAKA ) and Katsutosi TANIMOTO ) Member of JSCE, B.Sc.Eng
More informationLECTURE NOTES - VIII. Prof. Dr. Atıl BULU
LECTURE NOTES - VIII «LUID MECHNICS» Istanbul Technical Universit College of Civil Engineering Civil Engineering Department Hdraulics Division CHPTER 8 DIMENSIONL NLYSIS 8. INTRODUCTION Dimensional analsis
More informationUNIT 4 HEAT TRANSFER BY CONVECTION
UNIT 4 HEAT TRANSFER BY CONVECTION 4.1 Introduction to convection 4. Convection boundar laers 4..1 Hdrodnamic boundar laer over flat plate 4.. Thermal boundar laer 4..3 Concentration boundar laer 4.3 Dimensional
More informationTitle. Author(s)Leppäranta, Matti; Shirasawa, Kunio. Issue Date Doc URL. Type. Note. File Information.
Title Influence of the ice-ocean heat flux on the ice thic Author(s)Leppäranta, Matti; Shirasawa, Kunio Issue Date 7 Doc URL http://hdl.handle.net/115/38937 Type proceedings (author version) Note 16th
More informationData Repository Hampel et al., page 1/5
GSA DATA REPOSITORY 2138 Data Repositor Hampel et al., page 1/5 SETUP OF THE FINITE-ELEMENT MODEL The finite-element models were created with the software ABAQUS and consist of a 1-km-thick lithosphere,
More informationCombined tidal and wind driven flows and residual currents
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 Combined tidal and wind driven flows and residual currents 1 Lars Erik Holmedal, Hong Wang Abstract : The effect of a residual current on the combined
More informationCOMPUTATIONAL STUDY OF THE SEPARATED FLOW STRUCTURE INDUCED BY THE SYNTHETIC JET ON A BACKWARD-FACING STEP
COMPUTATIONAL STUDY OF THE SEPARATED FLOW STRUCTURE INDUCED BY THE SYNTHETIC JET ON A BACKWARD-FACING STEP Koichi Okada, Koo Fujii and Koji Miaji ABSTRACT In order to clarif the mechanism of the snthetic
More informationTIME EVOLUTION OF SCOUR AND DEPOSITION AROUND A CYLINDRICAL PIER IN STEADY FLOW
A-7 Fourth International Conference on Scour and Erosion 28 TIME EVOLUTION OF SCOUR AND DEPOSITION AROUND A CYLINDRICAL PIER IN STEADY FLOW Shina UMEDA 1,Tatsuo YAMAZAKI 2, Hajime ISHIDA 1 Member of JSCE,
More informationCHANGES IN RADIATION PROPERTIES AND HEAT BALANCE WITH SEA ICE GROWTH IN SAROMA LAGOON AND THE GULF OF FINLAND
Ice in the Environment: Proceedings of the 16th IAHR International Symposium on Ice Dunedin, New Zealand, 2nd 6th December 22 International Association of Hydraulic Engineering and Research CHANGES IN
More informationUNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX METHOD
Proceedings of the 3rd ASME/JSME Joint Fluids Engineering Conference Jul 8-23, 999, San Francisco, California FEDSM99-8 UNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX
More informationNUMERICAL COMPUTATION OF PARTICLE TRANSPORT IN TURBULENT FLOW FIELD IN ELECTROSTATIC PRECIPITATOR
NUMERICAL COMPUTATION OF PARTICLE TRANSPORT IN TURBULENT FLOW FIELD IN ELECTROSTATIC PRECIPITATOR Jenő SUDA, István BERTA, Gergel KRISTÓF +, Tamás LAJOS # Ph.D. Student, + Assistant, # Professor, Department
More informationCONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL
Mathematical and Computational Applications,Vol. 15, No. 4, pp. 742-761, 21. c Association for Scientific Research CONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL R. Naz 1,
More informationVERTICAL ICE FORCES ON PILE STRUCTURES UNDER WATER LEVEL CHANGES
Symposium on Ice (6) VERTICAL ICE FORCES ON PILE STRUCTURES UNDER WATER LEVEL CHANGES Takashi Terahima 1, Naoki Nakazawa, Shinji Kioka, Norihiro Usami 4 and Hiroshi Saeki 5 1 Kumashiro System Frontier
More informationThe Numerical Analysis and Experiment of Shock Processing for Bouef
Int. Jnl. of Multiphsics Volume 4 Number 4 2010 329 The Numerical Analsis and Eperiment of Shock Processing for Bouef Yusuke YAMASHITA*, Asuka ODA*, Takeshi HUJII*, Hironori MAEHARA*, M Moatamedi**, Shigeru
More informationQUANTITATIVE ANALYSIS ON THE TSUNAMI INUNDATION AREA AND EROSION ALONG THE ABUKUMAGAWA RIVER AT THE GREAT EAST JAPAN EARTHQUAKE
UANTITATIVE ANALYSIS ON THE TSUNAMI INUNDATION AREA AND EROSION ALONG THE ABUKUMAGAWA RIVER AT THE GREAT EAST JAPAN EARTHUAKE Masauki Sato Graduate School of Science and Engineering Saitama Universit,
More informationSIMULATION OF A LARGE-EDDY-BREAK-UP DEVICE (LEBU) IN A MODERATE REYNOLDS NUMBER TURBULENT BOUNDARY LAYER
June - Jul 3, 15 Melbourne, Australia 9 1A-3 SIMULATION OF A LARGE-EDDY-BREAK-UP DEVICE (LEBU) IN A MODERATE REYNOLDS NUMBER TURBULENT BOUNDARY LAYER Cheng Chin Universit of Melbourne Parkville, Victoria,
More informationFast waves during transient flow in an asymmetric channel
Int. Conference on Boundar and Interior Laers BAIL 6 G. Lube, G. Rapin (Eds) c Universit of Göttingen, German, 6. Introduction Fast waves during transient flow in an asmmetric channel Dick Kachuma & Ian
More informationThe wind-driven models of Stommel and Munk employed a linearization involving a small parameter, the Rossby number, which we need to reconsider.
Equatorial twists to mid-latitude dnamics As we saw or Stommel s or Munk s wind-driven gres and or Sverdrup s balance, there was no particular problem with the equator. In act, Stommel solved his gre or
More informationA Study on Residual Flow in the Gulf of Tongking
Journal of Oceanography, Vol. 56, pp. 59 to 68. 2000 A Study on Residual Flow in the Gulf of Tongking DINH-VAN MANH 1 and TETSUO YANAGI 2 1 Department of Civil and Environmental Engineering, Ehime University,
More informationSound Propagation in Ducts
Sound Propagation in Ducts Hongbin Ju Department of Mathematics Florida State Universit, Tallahassee, FL.3306 www.aeroacoustics.info Please send comments to: hju@math.fsu.edu In this section we will discuss
More informationPerformance Optimization of Jet Actuator Arrays for Active Flow Control
Performance Optimization of Jet Actuator Arras for Active Flow Control P. Scholz, J. Ortmanns, C.J. Kähler and R. Radespiel Technische Universität Braunschweig, Institut für Strömungsmechanik, Bienroder
More informationGrowing and decaying processes and resistance of sand waves in the vicinity of the Tone River mouth
Advances in River Sediment Research Fukuoka et al. (eds) 2013 Taylor & Francis Group, London, ISBN 978-1-138-00062-9 Growing and decaying processes and resistance of sand waves in the vicinity of the Tone
More informationThe Application of a Black-Box Solver with Error Estimate to Different Systems of PDEs
The Application of a Black-Bo Solver with Error Estimate to Different Sstems of PDEs Torsten Adolph and Willi Schönauer Institut für Wissenschaftliches Rechnen, Forschungszentrum Karlsruhe torstenadolph@iwrfzkde
More informationAbstract. 1. Introduction
Factors Influencing the Coefficient of Friction between Sea Ice and Various Materials Takashi Terashima & Naoki Nakazawa Pacific Consultants, Co., Ltd., New Stage Sapporo Bldg., 2-6, Kita-7 Nishi-1, Kita-ku
More informationTurbulent Mean Flow Effects: Inclusion of Rotation
Turbulent Mean Flow Effects: Inclusion of Rotation Ocean Edd L U H Horizontal Equation of Motion Du 1 1 f u h p Dt z where D u u u ' w' kz Dt t z Added horizontal Friction (eddies) u u 1 v u u p u u u
More informationPan Pearl River Delta Physics Olympiad 2005
1 Jan. 29, 25 Morning Session (9 am 12 pm) Q1 (5 Two identical worms of length L are ling on a smooth and horizontal surface. The mass of the worms is evenl distributed along their bod length. The starting
More informationEffects of wavelength and amplitude of a wavy cylinder in cross-flow at low Reynolds numbers
J. Fluid Mech. (2009), vol. 620, pp. 195 220. c 2009 Cambridge Universit Press doi:10.1017/s0022112008004217 Printed in the United Kingdom 195 Effects of wavelength and amplitude of a wav clinder in cross-flow
More informationAnalysis of Effects of Rebounds and Aerodynamics for Trajectory of Table Tennis Ball
SICE Annual Conference 21 August 18-21, 21, The Grand Hotel, Taipei, Taiwan Analsis of Effects of Reounds and Aerodnamics for Trajector of Tale Tennis Junko Nonomura Mechanical Science and Engineering,
More informationJoule Heating Effects on MHD Natural Convection Flows in Presence of Pressure Stress Work and Viscous Dissipation from a Horizontal Circular Cylinder
Journal of Applied Fluid Mechanics, Vol. 7, No., pp. 7-3, 04. Available online at www.jafmonline.net, ISSN 735-357, EISSN 735-3645. Joule Heating Effects on MHD Natural Convection Flows in Presence of
More informationAnalysis of condenser shell side pressure drop based on the mechanical energy loss
Article Calorifics December 01 Vol.57 No.36: 4718475 doi: 10.1007/s11434-01-5558-4 Analsis of condenser shell side pressure drop based on the mechanical energ loss ZENG Hui, MENG JiAn & LI ZhiXin * Ke
More informationLarge Eddy Simulation of the flow and heat transfer in a half-corrugated channel with various wave amplitudes
Int. J Heat Mass Transfer, Vol. 76 pp. 432-446, 214 Large Edd Simulation of the flow and heat transfer in a half-corrugated channel with various wave amplitudes Maram Mirzaei, Ahmad Sohankar, Lars Davidson,
More informationVALIDATION OF CALCULATIONS OF THE VISCOUS FLOW AROUND A SHIP IN OBLIQUE MOTION
The First MARIN-NMRI Workshop in Toko (Oktober 25-26, 24) VALIDATION OF CALCULATIONS OF THE VISCOUS FLOW AROUND A SHIP IN OBLIQUE MOTION Serge TOXOPEUS Maritime Research Institute Netherlands (MARIN),
More informationThe General Circulation of the Oceans
The General Circulation of the Oceans In previous classes we discussed local balances (Inertial otion, Ekman Transport, Geostrophic Flows, etc.), but can we eplain the large-scale general circulation of
More informationResearch Article Development of a Particle Interaction Kernel Function in MPS Method for Simulating Incompressible Free Surface Flow
Journal of Applied Mathematics Volume 2, Article ID 793653, 6 pages doi:.55/2/793653 Research Article Development of a Particle Interaction Kernel Function in MPS Method for Simulating Incompressible Free
More informationDiffusion of Radioactive Materials in the Atmosphere
American Journal of Environmental Sciences (): 3-7, 009 ISSN 3-34X 009 Science Publications Diffusion of Radioactive Materials in the Atmosphere Abdul-Wali Ajlouni, Manal Abdelsalam, Hussam Al-Rabai ah
More informationAalborg Universitet. Openings in Wind Driven Natural Ventilation True, Jan Per Jensen. Publication date: 2003
Aalborg Universitet Openings in Wind Driven Natural Ventilation True, Jan Per Jensen Publication date: 003 Document Version Earl version, also known as pre-print Link to publication from Aalborg Universit
More informationAN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES *
Journal of Computational Mathematics, Vol.27, No.2-3, 29, 388 399. AN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES * Yin Yang Hunan Ke Laborator for Computation and Simulation
More informationMEASUREMENTS OF BED SHEAR STRESS, VELOCITY AND REYNOLDS STRESS AROUND ARRAYS OF VERTICAL CYLINDERS
Full Paper 11 th ISE 216, Melbourne, Australia MEASUREMENTS OF BED SHEAR STRESS, VELOCITY AND REYNOLDS STRESS AROUND ARRAYS OF VERTICAL CYLINDERS KENJIROU HAYASHI Dept. of Civil and Environmental Eng.,
More informationOcean Dynamics. Equation of motion a=σf/ρ 29/08/11. What forces might cause a parcel of water to accelerate?
Phsical oceanograph, MSCI 300 Oceanographic Processes, MSCI 5004 Dr. Ale Sen Gupta a.sengupta@unsw.e.au Ocean Dnamics Newton s Laws of Motion An object will continue to move in a straight line and at a
More informationScale-Adaptive Simulation (SAS) Turbulence Modeling. F.R. Menter, ANSYS Germany GmbH
Scale-Adaptive Simulation (SAS) Turbulence Modeling F.. Menter, ANSYS German GmbH 1 Unstead ANS Based Models UANS (Unstead enolds averaged Navier Stoes) Methods UANS gives unphsical single mode unstead
More informationLift Enhancement on Unconventional Airfoils
Lift Enhancement on nconventional Airfoils W.W.H. Yeung School of Mechanical and Aerospace Engineering Nanang Technological niversit, North Spine (N3), Level 2 50 Nanang Avenue, Singapore 639798 mwheung@ntu.edu.sg
More informationModule 1 : The equation of continuity. Lecture 4: Fourier s Law of Heat Conduction
1 Module 1 : The equation of continuit Lecture 4: Fourier s Law of Heat Conduction NPTEL, IIT Kharagpur, Prof. Saikat Chakrabort, Department of Chemical Engineering Fourier s Law of Heat Conduction According
More informationNovel Approach to Prompt Fission Neutron Investigation
Novel Approach to Prompt ission Neutron Investigation Sh. Zenalov, O. Zenalova, P.Sedshev, V. Shvetsov JIN-Joint Institute for Nuclear esearch, Dubna, ussia ISINN, Alushta Ukraine, Ma -5, 3 Motivation
More informationInternal Tides in the Bab el Mandab Strait. Ewa Jarosz and Cheryl Ann Blain Naval Research Laboratory, Stennis Space Center, MS
Internal Tides in the Bab el Mandab Strait Ewa Jarosz and Cheryl Ann Blain Naval Research Laboratory, Stennis Space Center, MS Project Objectives Numerical simulation of tides and exchange flow in the
More informationSPREADING OF OIL UNDER ICE COVERS: EFFECTS OF BOTTOM ROUGHNESS OF ICE
Ice in the Environment: Proceedings of the 16th IAHR International Symposium on Ice Dunedin, New Zealand, 2nd 6th December 2002 International Association of Hydraulic Engineering and Research SPREADING
More informationSOUND SCATTERING FROM ACOUSTICALLY TREATED BODIES USING A DOMAIN DECOMPOSITION TECHNIQUE
th AIAA/CEAS Aeroacoustics Conference (th AIAA Aeroacoustics Conference) - Ma 9, Miami, Florida th AIAA/CEAS Aeroacoustics Conference (th AIAA Aeroacoustics Conference) AIAA 9-9 AIAA-9-9 SOUND SCATTERING
More informationMechanics Departmental Exam Last updated November 2013
Mechanics Departmental Eam Last updated November 213 1. Two satellites are moving about each other in circular orbits under the influence of their mutual gravitational attractions. The satellites have
More informationAn ellipsoidal particle in tube Poiseuille flow
J. Fluid Mech. (217), vol. 822, pp. 664 688. c Cambridge Universit Press 217 doi:1.117/jfm.217.298 An ellipsoidal particle in tube Poiseuille flow 664 Haibo Huang 1, and Xi-Yun Lu 1 1 Department of Modern
More informationMATHEMATICS 200 December 2011 Final Exam Solutions
MATHEMATICS December 11 Final Eam Solutions 1. Consider the function f(, ) e +4. (a) Draw a contour map of f, showing all tpes of level curves that occur. (b) Find the equation of the tangent plane to
More informationRESPONSE OF TWO PIERS ON CONFEDERATION BRIDGE TO ICE LOADING EVENT OF APRIL 4, 2003
RESPONSE OF TWO PIERS ON CONFEDERATION BRIDGE TO ICE LOADING EVENT OF APRIL 4, 23 R. Frederking 1, I. Kubat 1 and S. Prinsenberg 2 1 National Research Council, Canadian Hydraulics Centre, Ottawa, ON, Canada
More informationImplementation of Bridge Pilings in the ADCIRC Hydrodynamic Model: Upgrade and Documentation for ADCIRC Version 34.19
Implementation of Bridge Pilings in the ADCIRC Hydrodynamic Model: Upgrade and Documentation for ADCIRC Version 34.19 Richard A. Luettich, Jr. University of North Carolina at Chapel Hill Institute of Marine
More informationOn the relationship between the vortex formation process and cylinder wake vortex patterns
J. Fluid Mech. (4), vol. 59, pp. 68. c 4 Cambridge Universit Press DOI:.7/S48 Printed in the United Kingdom 6 On the relationship between the vorte formation process and clinder wake vorte patterns B DAVID
More informationNUMERICAL INVESTIGATION AND EXPERIMENTAL VALIDATION OF HEAT TRANSFER IN A SMALL SIZE SHELL AND TUBE HEAT EXCHANGER
1 M. Kirinčić, A. rp, K. Lenić: Numerical investigation and NUMERICAL INVESIGAION AND EXPERIMENAL VALIDAION OF HEA RANSFER IN A SMALL SIZE SHELL AND UBE HEA EXCHANGER Mateo Kirinčić * Anica rp Kristian
More informationChapter 2 Basic Conservation Equations for Laminar Convection
Chapter Basic Conservation Equations for Laminar Convection Abstract In this chapter, the basic conservation equations related to laminar fluid flow conservation equations are introduced. On this basis,
More informationFluid Mechanics for International Engineers HW #4: Conservation of Linear Momentum and Conservation of Energy
2141-365 Fluid Mechanics for International Engineers 1 Problem 1 RTT and Time Rate of Change of Linear Momentum and The Corresponding Eternal Force Notation: Here a material volume (MV) is referred to
More informationSTUDY ON THERMAL FUNCTION OF IVY-COVERED WALLS
STUDY ON THERMAL FUNCTION OF IVY-COVERED WALLS Mr. Liao Zaii, Dr. Niu J.L Dept. of Building Services Engineering The Hong Kong Poltechnic Universit Hunghum, Kowloon, Hong Kong E-mail: bejlniu@polue.edu.hk
More informationEstimation Of Linearised Fluid Film Coefficients In A Rotor Bearing System Subjected To Random Excitation
Estimation Of Linearised Fluid Film Coefficients In A Rotor Bearing Sstem Subjected To Random Ecitation Arshad. Khan and Ahmad A. Khan Department of Mechanical Engineering Z.. College of Engineering &
More informationA QUALITATIVE ANALYSIS OF BREAKING LENGTH OF SHEET ICE AGAINST CONICAL STRUCTURE
A QUALITATIVE ANALYSIS OF BREAKING LENGTH OF SHEET ICE AGAINST CONICAL STRUCTURE Li Feng 1, Yue Qianjin 1, Karl N. Shkhinek 2, Tuomo Kärnä 3 1 Dalian University of Technology, China 2 St. Petersburg State
More informationUnsteady free MHD convection flow past a vertical porous plate in slip-flow regime under fluctuating thermal and mass diffusion *
CHAPTER 4 Unstead free MHD convection flow past a vertical porous plate in slip-flow regime under fluctuating thermal and mass diffusion * 4. Nomenclature A : Suction parameter C : Dimensionless species
More informationDETERMINATION OF HORIZONTAL AXIS WIND TURBINE PERFORMANCE IN YAW BY USE OF SIMPLIFIED VORTEX THEORY
5 th INTERNATIONAL MULTIDISCIPLINARY CONFERENCE DETERMINATION OF HORIZONTAL AXIS WIND TURBINE PERFORMANCE IN YAW BY USE OF SIMPLIFIED ORTEX THEORY Piotr Strzelczk, PhD Eng. Dept. Of Fluid Mechanics and
More informationWind Parameters from Scatterometry
Wind Parameters from Scatterometr Zorana Jelenak, Moira Sten-Ross and Aleksandar Jelenak Department of Phsics and Electronic Engineering Universit of Waikato Private Bag 3105 Hamilton, New Zealand Email:
More informationExperimental Uncertainty Review. Abstract. References. Measurement Uncertainties and Uncertainty Propagation
Experimental Uncertaint Review Abstract This is intended as a brief summar of the basic elements of uncertaint analsis, and a hand reference for laborator use. It provides some elementar "rules-of-thumb"
More informationSTUDY ON BRACKISH ICE IN THE GULF OF FINLAND
Ice in the Environment: Proceedings of the 16th IAHR International Symposium on Ice Dunedin, New Zealand, 2nd 6th December 2002 International Association of Hydraulic Engineering and Research STUDY ON
More informationA first contact with STAR-CCM+ Comparing analytical and finite volume solutions with STAR-CCM+ simulations
A first contact with STAR-CCM+ Comparing analtical and finite volume solutions with STAR-CCM+ simulations Michael Heer Analtical Finite volumes STAR-CCM+ What is ParisTech? ParisTech is a consortium of
More information1. Research Institute for Applied Mechanics, Kyushu University, JAPAN 2. Kagoshima University, JAPAN 3. Tokyo University of Science, JAPAN
An estimate of the tsunami-debris quantity washed ashore on the US and Canadian beaches, based on a webcam monitoring and a particle tracking model experiment Atsuhiko Isobe 1, Shinsuke Iwasaki 1, Shin
More informationCH.7. PLANE LINEAR ELASTICITY. Multimedia Course on Continuum Mechanics
CH.7. PLANE LINEAR ELASTICITY Multimedia Course on Continuum Mechanics Overview Plane Linear Elasticit Theor Plane Stress Simplifing Hpothesis Strain Field Constitutive Equation Displacement Field The
More informationAll work must be shown in this course for full credit. Unsupported answers may receive NO credit.
AP Calculus. Worksheet All work must be shown in this course for full credit. Unsupported answers ma receive NO credit.. What is the definition of a derivative?. What is the alternative definition of a
More informationCHAPTER-III CONVECTION IN A POROUS MEDIUM WITH EFFECT OF MAGNETIC FIELD, VARIABLE FLUID PROPERTIES AND VARYING WALL TEMPERATURE
CHAPER-III CONVECION IN A POROUS MEDIUM WIH EFFEC OF MAGNEIC FIELD, VARIABLE FLUID PROPERIES AND VARYING WALL EMPERAURE 3.1. INRODUCION Heat transer studies in porous media ind applications in several
More information2010 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media,
2010 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in an current or future media, including reprinting/republishing this material for advertising
More informationVortex Shedding on Combined Bodies at Incidence to a Uniform Air Stream T. Yavuz x, Y. E. Akansu xx, M. Sarıo lu xxx, and M.
Vorte Sheing on Combine Boies at Incience to a Uniform Air Stream T. Yavuz, Y. E. Akansu, M. Sarıo lu, an M. Özmert. Ba kent Universit, : Nige Universit,, : Karaeniz Technical Universit,Turke Abstract
More informationOCN/ATM/ESS 587. The wind-driven ocean circulation. Friction and stress. The Ekman layer, top and bottom. Ekman pumping, Ekman suction
OCN/ATM/ESS 587 The wind-driven ocean circulation. Friction and stress The Ekman layer, top and bottom Ekman pumping, Ekman suction Westward intensification The wind-driven ocean. The major ocean gyres
More informationPHYSICS PART II SECTION- I. Straight objective Type
PHYSICS PAT II SECTION- I Straight objective Tpe This section contains 9 multiple choice questions numbered to 1. Each question has choices,, (C) and, out of which ONLY ONE is correct.. A parallel plate
More informationGENERALISATION OF THE TWO-SCALE MOMENTUM THEORY FOR COUPLED WIND TURBINE/FARM OPTIMISATION
25 th National Symposium on Wind Engineering, Tokyo, Japan, 3-5 December 2018 第 25 回風工学シンポジウム (2018) GENERALISATION OF THE TWO-SCALE MOMENTUM THEORY FOR COUPLED WIND TURBINE/FARM OPTIMISATION Takafumi
More informationModeling the Columbia River Plume on the Oregon Shelf during Summer Upwelling. 2 Model
Modeling the Columbia River Plume on the Oregon Shelf during Summer Upwelling D. P. Fulton August 15, 2007 Abstract The effects of the Columbia River plume on circulation on the Oregon shelf are analyzed
More informationPHYSICAL REVIEW E 78, Lattice Boltzmann method for simulation of compressible flows on standard lattices
PHYSICAL REVIEW E 7, 016704 00 Lattice Boltzmann method for simulation of compressible flows on standard lattices Nikolaos I. Prasianakis 1, * and Ilia V. Karlin 1,, 1 Aerothermochemistr and Combustion
More informationDirect Numerical Simulation of Airfoil Separation Bubbles
Direct Numerical Simulation of Airfoil Separation Bubbles Ulrich Maucher, Ulrich Rist and Siegfried Wagner 1 Abstract. In the present paper results of Direct Numerical Simulations of Laminar Separation
More informationRecent Changes in Pressure Patterns and Their Regional Occurrence at Times of Heavy Snowfall and Blizzard Events in Hokkaido, Japan
Recent Changes in Pressure Patterns and Their Regional Occurrence at Times of Heavy Snowfall and Blizzard Events in Hokkaido, Japan Masaru Matsuzawa, Yusuke Harada, Satoshi Omiya, Hirotaka Takechi Civil
More informationCH. 1 FUNDAMENTAL PRINCIPLES OF MECHANICS
446.201 (Solid echanics) Professor Youn, eng Dong CH. 1 FUNDENTL PRINCIPLES OF ECHNICS Ch. 1 Fundamental Principles of echanics 1 / 14 446.201 (Solid echanics) Professor Youn, eng Dong 1.2 Generalied Procedure
More informationThe Gas-assisted Expelled Fluid Flow in the Front of a Long Bubble in a Channel
C. H. Hsu, P. C. Chen,. Y. ung, G. C. uo The Gas-assisted Epelled Fluid Flow in the Front of a Long Bubble in a Channel C.H. Hsu 1, P.C. Chen 3,.Y. ung, and G.C. uo 1 1 Department of Mechanical Engineering
More informationModeling Temporal and Spatial Dispersion of Smoke Plume Based on Digital Elevation Model (Part 1)
American Journal of Geographic Information Sstem 2013, 2(4): 82-92 DOI: 10.5923/j.ajgis.20130204.02 Modeling Temporal and Spatial Dispersion of Smoke Plume ased on Digital Elevation Model (Part 1) Carl
More informationSelected Papers from the WSEAS Conferences in Spain, September 2008 Santander, Cantabria, Spain, September 23-25, 2008
Comparisons between Velocit Profiles According to the Modified and the Navier-Stokes Equations of Motion and the Eperimental Measurements for Laminar Boundar Laer Over a Flat Plate Matthew J. Inkman, and
More informationResearch Design - - Topic 15a Introduction to Multivariate Analyses 2009 R.C. Gardner, Ph.D.
Research Design - - Topic 15a Introduction to Multivariate Analses 009 R.C. Gardner, Ph.D. Major Characteristics of Multivariate Procedures Overview of Multivariate Techniques Bivariate Regression and
More informationA Coin Detection System by Coupled Printed Spiral Inductors
A Coin Detection Sstem b Coupled Printed Spiral Inductors Junichi Fukatani, Sho Yamaguchi, Masauki Yamauchi, Kazuhisa Yoshimatsu, Hisashi Aomori and Mamoru Tanaka Department of Electrical and Electronics
More informationRECENT STUDIES OF THE SEA OF OKHOTSK
RECENT STUDIES OF THE SEA OF OKHOTSK Masaaki Wakatsuchi Institute of Low Temperature Science, Hokkaido University, Sapporo, Japan ABSTRACT We recently have had a Japan-Russia-United States International
More informationConsideration of Shock Waves in Airbag Deployment Simulations
Consideration of Shock Waves in Airbag Deploment Simulations Doris Rieger BMW Group ABSTRACT When the inflation process of a simple flat airbag was simulated with the MADYMO gas flow module, the resulting
More information