Conceptual Design of the Flow Noise Simulator

The Flow Noise Simulator (FNS) of the 1st Research Center of TRDI/JDA (Japan Defense Agency) is a large, variable pressure, recirculating water tunnel with very low background noise level. The tunnel is 20m high and 49m long, containing 2000m3 of water. The test section has a square cross section of 2m × 2m with 10m in length. It will accept large size surface ship models of 6m, submarine models of 4m in length and full scale ship appendix models. The FNS is currently under construction and will be accomplished in 2005. It will be used for a wide variety of hydrodynamic and hydroacoustic testing of surface ships and submarines, such as propeller cavitation noise measurements and propeller-hull interaction observation, with sufficiently large scale models. Conceptual design of the FNS was started in 1996 and evaluated by following scale model studies. This paper discusses some technical issues of the FNS.

Rational Polynomial Transfer Function Approximations for Fluid Transients in Lines

This paper introduces a simple approach utilizing MATLAB® computational tools for generating rational polynomial transfer functions for fluid transients in both liquid and gas fluid transmission lines. These transfer functions are obtained by curve fitting in the frequency domain the exact solution to the distributed parameter laminar flow “Dissipative Model” for fluid transients that includes nonlinear frequency dependent viscous friction terms as well as heat transfer effects in gas lines. These transfer functions are formulated so they are applicable to arbitrary line terminations and so they can be inserted directly into SIMULINK® models for time domain simulation and analysis of a total system of which the fluid lines are only internal components. The inputs to the algorithm are the internal radius and length of the line; the kinematic viscosity, density, Prandtl number, and speed of sound of the fluid; and the maximum frequency to which an accurate curve fit of the exact solution is desired. This maximum frequency normally is equal to or greater than the bandwidth of the other components in the total system to be analyzed or the maximum frequency associated with the input. The simplicity of use and accuracy in the results of the exact solution representations are demonstrated for examples of a blocked fluid line and of a line terminating into a tank. The computational algorithms are available for download from the Author’s web site. This is the first of two papers pertaining to transfer functions for fluid transients. The second paper pertains to formulating simulation diagrams for total systems containing fluid lines represented by rational polynomial transfer functions.

Why Are Relationships Between Lagrangian and Eulerian Scales Necessary for Gas-Particle Flow Modeling?

Simulation of Gas-Particle flows can be fulfilled by Lagrangian modeling of the dispersed phase. Each type of Lagrangian method, Monte-Carlo/Eddy Interaction or Markov Chain models, needs the knowledge of Lagrangian scales associated with the turbulent flow under consideration and the type of particle dispersing in the gas carrier flow. Unfortunately, Lagrangian quantities (as well the interesting moving Eulerian time scale, that given by a sensor which would move with the mean fluid velocity) are still difficult to be obtained directly by most experimental measurement techniques (except by very recent techniques such as PIV.PTV.), contrary to Eulerian scales scales, such as those classically obtained from a fixed hot-wire or LDA control volume. It is therefore of great importance to have available accurate relationships between Eulerian and Lagrangian scales, based on fluid flow properties as well as particle characteristics.

Three-Dimensional Velocity and Droplet Size Measurements of Spray Flow (Keynote Paper)

The present investigation describes an application of a novel technique of simultaneous measurement of droplet size and three-dimensional components of velocity in a high density spray with swirl. The spray has a complicated and three-dimensional structure caused by mixing with surrounding airflow entrained by high speed fuel jet issuing from a nozzle. The breakup process of fuel film to fine-droplet-cloud, the droplet size dispersion and the velocity distribution of droplets are important factors in practical application of fuel spray for combustors. The conventional technique can be applied to local measurement of droplet speed and size. Recent methods, based on optical and image processing techniques, provide measurement of the velocity and droplet size distribution in observation area or volume. Maeda et al. proposed an excellent measurement technique of the size and the velocity distribution of droplet in spray based on interferometric laser imaging in which the fringe pattern is generated at the out of focus plane by interference between 0th order and 1st order refractions of droplet illuminated by high power laser light sheet. And also, in this technique, the separation of overlapping droplets image has been successfully done by optical method. As a practical application, the size and velocity distributions of droplets in a high density spray without swirl have been measured by this technique. In general, the droplet motion in a spray field is highly three-dimensional. Especially, a spray generated by a swirl nozzle shows complicated droplet motion in the three-dimensional field. In order to analyze the configuration of a complicated spray field, three-dimensional velocity measurement of droplets must be required. In the present paper, a combined measurement technique of the size and three velocity components of droplets in three-dimensional spray field based on doublet imaging technique of droplets and stereoscopic PIV method has been developed. And its feasibility and applicability was confirmed by practical application to measurements of spray fields induced by a swirl jet nozzle using in gas turbine.

Numerical Analysis of Pipeline Equipment Effect on Water Hammer Using Characteristic Method

In this attempt effect of pipeline equipment behavior was considered on water hammer numerically. The effect includes opening / closing of the shut off valves, loss of coefficient of the outlet bypass pipe for the air chamber, elasticity of the pipeline and loss coefficient due to friction. In order to study the behavior, mass and momentum conservation equations were solved numerically using characteristic method during transient conditions. As a water hammer phenomena accompanies with large pressure gradient, so the pipeline equipment behavior and their effect were analyzed with respect to the maximum pressure occurrence. For a pipeline of 5000 m length, 1 m diameter, 1 m3/s discharge and 100 m height between upstream and downstream, the following result were concluded: 1-If the moment of inertia of the pump impeller increases by 400 percent, the maximum pressure occurred by the water hammer will decrease by 9 percent. 2-During on and off of the shut off valve, 80 percent of pressure increase due to water hammer was created during the last 15 percent of valve closure. 3-If pressure wave velocity increases by 75 percent, then the maximum pressure generated due to the water hammer will increase by 27 percent. 4-If the loss coefficient of the by pass line of the air chamber decreases by 90 percent, then the maximum pressure due to the water hammer will decrease by 20 percent. 5-If the pipeline Moody friction coefficient increases by 92 percent, the maximum pressure due to the water hammer will increase by 66 percent.

Numerical Simulation of Water Quenching Process of an Engine Cylinder Head

Presented is a simulation of an engine cylinder head undergoing water quenching process using a recently developed approach for modeling quenching cooling of metal parts (Wang et al., 2002). The approach is based on the AVL SWIFT Eulerian two-fluid method with special emphasis on handling high mass exchange rate associated with quenching. A tetrahedral grid of 830,000 cells is generated for the computational domain, which includes the solid part of the cylinder head immersed in the fluid. Detailed vapor and temperature distributions are obtained which offer valuable information for the thermal stress analysis. It is observed that the temperature field within the cylinder head is highly non-uniform. The computed cylinder head monitoring point temperature versus time is compared with that registered by the thermal couple measurement. Reasonable agreement is observed. The simulation exercise may potentially be used to identify the cause of cracks often encountered in quenching heat treatment thereby lead to a better design of the process.

Effects of Weak Free Stream Nonuniformity on Boundary Layer Transition (Keynote Paper)

Experiments are described in which well-defined FSN (Free Stream Nonuniformity) distributions are introduced by placing fine wires upstream of the leading edge of a flat plate. Large amplitude spanwise thickness variations are present in the downstream boundary layer resulting from the interaction of the laminar wakes with the leading edge. Regions of elevated background unsteadiness appear on either side of the peak layer thickness, which share many of the characteristics of Klebanoff modes, observed at elevated Free Stream Turbulence (FST) levels. However, for the low background disturbance level of the free stream, the layer remains laminar to the end of the test section (Rx ≈ l.4×106) and there is no evidence of bursting or other phenomena associated with breakdown to turbulence. A vibrating ribbon apparatus is used to demonstrate that the deformation of the mean flow is responsible for substantial phase and amplitude distortion of Tollmien-Schlichting (TS) waves. Pseudo-flow visualization of hot-wire data shows that the breakdown of the distorted waves is more complex and occurs at a lower Reynolds number than the breakdown of the K-type secondary instability observed when the FSN is not present.

Simulation of Vortex Instability in a Pipe Trifurcation

The vortex instability in a spherical pipe trifurcation is investigated by applying a Very Large Eddy Simulation (VLES). For this approach an new adaptive turbulence model based on an extended version of the k-ε model is used. Applying a classical Reynolds-averaged Navier-Stokes-Simulation with the standard k-ε model is not able to forecast the vortex instability. However the prescribed VLES method is capable to predict this flow phenomenon. The obtained results show a reasonable agreement with measurements in a model test.

Dynamics of Bubbles Rising in Finite and Infinite Media

The dynamic behavior of single bubbles rising in quiescent liquid Suva (R134a) in a duct has been examined through the use of a high speed video system. Size, shape and velocity measurements obtained with the video system reveal a wide variety of characteristics for the bubbles as they rise in both finite and infinite media. This data, coupled with previously published data for other working fluids, has been used to assess and extend a rise velocity model given by Fan and Tsuchiya. As a result of this assessment, a new rise velocity model has been developed which maintains the physically consistent characteristics of the surface tension in the distorted bubbly regime. In addition, the model is unique in that it covers the entire range of bubble sizes contained in the spherical, distorted and planar slug regimes.

Transfer, Segregation and Deposition of Heavy Particles in Horizontal and Vertical Turbulent Channel Flows

Particle transfer in the wall region of turbulent boundary layers is dominated by the coherent structures which control the turbulence regeneration cycle. Coherent structures bring particles toward the wall and away from the wall and favour particle segregation in the viscous region giving rise to nonuniform particle distribution profiles which peak close to the wall. In this work, we focus on the transfer mechanism of different size particles and on the influence of gravity on particles deposition. By tracking O(105) particles in Direct Numerical Simulation (DNS) of a turbulent channel flow at Reτ = 150, we find that particles may reach the wall directly or may accumulate in the wall region, under the low-speed streaks. Even though low-speed streaks are ejection-like environments, particles are not re-entrained into the outer region. Particles segregated very near the wall by the trapping mechanisms we investigated in a previous work [1] are slowly driven to the wall. We find that gravity plays a role on particle distribution but, for small particles (τp+ < 3), the controlling transfer mechanism is related to near-wall turbulence structure.