Velocity measurements of flow through a step stenosis using Magnetic Resonance Imaging

 Magnetic resonance imaging (MRI) is a versatile noninvasive tool for achieving full-field quantitative visualization of complex fluid flows. The MRI signal results from the interaction of radio-frequency (RF) pulses with nuclear spins exposed to a strong static magnetic field. The two main techniques of MRI velocimetry are time-of-flight and phase contrast techniques. Time-of- flight techniques involve tagging and tracking a material volume of fluid, whereas phase contrast techniques use magnetic field gradients to encode velocity information into the phase of the MRI signal. In this study, both techniques are used to probe the pressure-driven steady flow of water in a pipe with a step stenosis. The velocity measurements were then compared with computational results obtained using the FIDAP software package. The experiments show that the phase contrast method gives more accurate results, with 90% of the measurements within 10% of the local computational fluid dynamics (CFD) velocity predictions at Re = 100 and 94% of the measurements within 10% of the local CFD predictions at Re = 258. Although the time-of-flight experiments were not as accurate, they provide a good qualitative image of the flow field. Sources of the discrepancies between the MRI data and the CFD results are also discussed, including acceleration and spin flow-through artifacts. Received: 7 April 1999/Accepted: 20 December 1999     

Experimental investigation of integral and local flow field properties on rotating porous disc

Flow field through rotating porous disc was investigated with experimental fluid dynamics and compared with computational fluid dynamics. Open cell aluminum metal foam with 88% porosity was used. On rotating porous disc, integral measurements of static pressure difference in dependence of air volume flow rate were performed. Local measurements of velocity profiles close to disc circumference were performed with hot-wire anemometry. The airflow visualization method using smoke generator and digital camera was performed. Flow structures through porous disc were visualized at three different air volume flow rates. Numerical simulation of homogenous rotating porous disc was performed. Experimental and numerical results were compared. The results showed appropriate comparison of integral and local properties. The presented experimental approach can be used for the investigation and understanding of flow field phenomena on rotating porous materials. The proposed conclusions can be applied for variable applications on rotating porous materials.     

RANS simulations for transition and turbulent flow regimes in wire-wrapped rod bundles

Internal flows through rod assemblies are commonly found in heat exchangers, steam generators, and nuclear reactors. One of the fuel assembly designs considered for liquid metal-cooled reactors utilizes wires helically wrapped around each fuel rod as spacers. The wires keep the fuel pins separated, enhancing the turbulent mixing, and heat transfer, but also affecting the pressure drop. It is of interest the understanding of the fluid flow phenomena in the sodium-cooled fast reactor as it is one of the Generation IV advanced reactor designs and it has been a motivating topic of research for the last decade. A wire-wrapped fuel assembly replica with 61-pins has been in operation at the Thermal-Hydraulic Research Laboratory of Texas A&M University. This facility produced high-fidelity velocity and pressure drop data for validation of computational fluid dynamics codes. This study investigates the effects of geometrical features and operating conditions on the flow behavior of the 61-pin wire wrapped bundle using Reynolds-Averaged Navier-Stokes (RANS) models to predict the axial and transverse pressure drops for a range of Reynolds numbers from 1,270 to 100,000. The friction factor predictions were in satisfactory agreement with the experimental data and the Upgraded Chen and Todreas correlation. The internal subchannel velocity results were compared with experimental data and Large Eddy Simulations (LES) and found in reasonable agreement. This study demonstrates that RANS is a suitable approach in predicting velocity and pressure fields in wire-wrapped rod bundles, with a relatively low computational effort.     

Modelling Pressurized Water Reactor cores in terms of porous media

The aim of this study is to develop a tractable model of a nuclear reactor core taking the complexity of the structure (including its nonlinear behaviour) and fluid flow coupling into account. The mechanical behaviour modelling includes the dynamics of both the fuel assemblies and the fluid. Each rod bundle is modelled in the form of a deformable porous medium; then, the velocity field of the fluid and the displacement field of the structure are defined over the whole domain. The fluid and the structure are first modelled separately, before being linked together. The equations of motion for the structure are obtained using a Lagrangian approach and, to be able to link up the fluid and the structure, the equations of motion for the fluid are obtained using an arbitrary Lagragian Eulerian approach. The finite element method is applied to spatially discretize the equations. Simulations are performed to analyse the effects of the characteristics of the fluid and of the structure. Finally, the model is validated with a test involving two fuel assemblies, showing good agreement with the experimental data.     

Permeability mapping in porous media by magnetization prepared centric-scan SPRITE

The ability of porous media to transmit fluids is commonly referred to as permeability. The concept of permeability is central for hydrocarbon recovery from petroleum reservoirs and for studies of groundwater flow in aquifers. Spatially resolved measurements of permeability are of great significance for fluid dynamics studies. A convenient concept of local Darcy’s law is suggested for parallel flow systems. The product of porosity and mean velocity images in the plane across the average flow direction is directly proportional to permeability. Single Point Ramped Imaging with T ₁ Enhancement (SPRITE) permits reliable quantification of local fluid content and flow in porous media. It is particularly advantageous for reservoir rocks characterized by fast magnetic relaxation of a saturating fluid. Velocity encoding using the Cotts pulsed field gradient scheme improves the accuracy of measured flow parameters. The method is illustrated through measurements of 2D permeability maps in a capillary bundle, glass bead packs and composite sandstone samples.     

High-flux magnetorheology at elevated temperatures

Commercial applications of magnetorheological (MR) fluids often require operation at elevated temperatures as a result of surrounding environmental conditions or intense localized viscous heating. Previous experimental investigations of thermal effects on MR fluids have reported significant reductions in the magnetorheological stress with increasing temperature, exceeding the predictions made by considering the thermal variations in the individual physical properties of the fluid and solid constituents of a typical MR fluid. In the low-flux regime, designers of MR fluid actuators can alleviate this thermal reduction in stress by increasing the applied magnetic field strength. However, this is not possible in the high-flux regime because of magnetic saturation, and it becomes necessary to explore and understand the intrinsic limitations of the fluid at elevated temperature. We describe a new magnetorheological fixture, which was designed as an accessory to a commercial torsional shear rheometer, capable of applying magnetic flux densities up to 1 T and controlling the sample temperature up to 150°C. During the design of the instrument, close attention was given to the uniformity of the magnetic field applied to the sample by using numerical simulations. Incorporation of a custom-built magnetic flux sensor which matches the environmental capabilities of the fixture enables in situ measurement of the local magnetic field at each temperature. The numerical results are also validated by spatially resolved measurements of the local magnetic field. Finally, we explore the ability of a shift factor between fluid magnetization and yield strength to describe the measured variation in the MR fluid response at elevated temperatures.     

A comprehensive PIV measurement campaign on a fully equipped helicopter model

The flow field around a helicopter is characterised by its inherent complexity including effects of fluid?Cstructure interference, shock?Cboundary layer interaction, and dynamic stall. Since the advancement of computational fluid dynamics and computing capabilities has led to an increasing demand for experimental validation data, a comprehensive wind tunnel test campaign of a fully equipped and motorised generic medium transport helicopter was conducted in the framework of the GOAHEAD project. Different model configurations (with or without main/tail rotor blades) and several flight conditions were investigated. In this paper, the results of the three-component velocity field measurements around the model are surveyed. The effect of the interaction between the main rotor wake and the fuselage for cruise/tail shake flight conditions was analysed based on the flow characteristics downstream from the rotor hub and the rear fuselage hatch. The results indicated a sensible increment of the intensity of the vortex shedding from the lower part of the fuselage and a strong interaction between the blade vortex filaments and the wakes shed by the rotor hub and by the engine exhaust areas. The pitch-up phenomenon was addressed, detecting the blade tip vortices impacting on the horizontal tail plane. For high-speed forward flight, the shock wave formation on the advancing blade was detected, measuring the location on the blade chord and the intensity. Furthermore, dynamic stall on the retreating main rotor blade in high-speed forward flight was observed at r/R?=?0.5 and 0.6. The analysis of the substructures forming the dynamic stall vortex revealed an unexpected spatial concentration suggesting a rotational stabilisation of large-scale structures on the blade.     

Cascade aerodynamic gust response including steady loading effects

To predict the unsteady convected gust aerodynamic response of a cascade comprised of arbitrary thick and cambered aerofoils in an incompressible, inviscid, flow field, a complete first-order model is formulated. The flow is analysed by considering a periodic flow channel. The velocity potential is separated into steady and unsteady harmonic components, each described by a Laplace equation. The strong dependence of the unsteady aerodynamics on the steady effects of aerofoil and cascade geometry and incidence angle is manifested in the coupling of the unsteady and steady flow fields through the unsteady boundary conditions. Analytical solutions in individual grid elements of a body-fitted computational grid are then determined, with the complete solution obtained by assembly of these local solutions. The validity and capabilities of this model and solution technique are then demonstrated by analysing the steady and unsteady aerodynamics of both theoretical and experimental cascade configurations.     

Direct numerical simulation of fluid flow in a 5x5 square rod bundle

Characterization of parallel flow through rod bundles is of key importance in assessing the performance and safety of several engineering systems, including a majority of nuclear reactor concepts. Inhomogeneities in the bundle cross-section can present complex flow phenomena, including varying local conditions of turbulence. With the ever-increasing capabilities of high-performance computing, Direct Numerical Simulation (DNS) of turbulent flows is becoming more feasible. Through resolving all scales of turbulence, DNS can serve as a “numerical experiment,” and can provide substantial insight into flow physics, but at considerable computational cost. Thus to date, the DNS in open literature for rod bundle flows is relatively scarce, and largely limited to unit-cell domains. Since wall effects are important in rod bundle flows, a multiple-pin DNS study can expand understanding of rod bundle flows while providing valuable reference data for evaluating reduced-resolution techniques. In this work, DNS of a 5x5 square bare rod bundle representative of typical light water reactor fuel dimensions was performed using the spectral element code Nek5000. Turbulent microscales based on an advanced Reynolds-Averaged Navier–Stokes model were used to establish the required DNS resolution. Velocity and Reynolds stress fields are analyzed in detail, and invariant analysis is used for further investigation into flow physics. The results show stark changes in the structure of turbulence in the edge gaps, suggesting the presence of gap vortices in these regions. In addition, turbulent kinetic energy budgets are presented to more fully illustrate the various turbulent processes. These data can prove useful for rigorous evaluation of lower-fidelity turbulence modeling approaches.     

Numerical and experimental analysis of flow and heat transfer in a fuel assembly mock-up with transverse flow above the rods

Thermal-hydraulic conditions in a partially uncovered nuclear fuel assembly mock-up are studied with particular focus on the influence of the horizontal air flow above the rod bundle. The investigations are performed at the ALADIN test facility, which models a boiling water reactor fuel assembly at a 1:1 scale both axially and radially. In the scenario studied, the main heat transfer mechanisms – conduction, convection and radiation – are strongly coupled and all are of similar importance. A combination of measurements and CFD simulations serves to analyze the heat transfer processes in detail. Contrary to previous studies in this field, all heat transfer mechanisms were considered in the simulation with sophisticated models. The numerical results show a good agreement with the measurements, given the inevitable differences between the approaches. Although the successive evaporation of cooling water in a fuel assembly is a transient, multiphase process, the steady, single-phase simulation yields acceptable results. While single effects are overestimated in the simulation, the important dependencies are predicted similarly. A general result is that the maximum cladding temperature rises with decreasing water level. Further results indicate an impact of the horizontal air flow on the residual heat removal for moderate rod powers. Higher horizontal velocities above the fuel assembly lead to slightly higher temperatures inside. A characteristic flow field forms in the test facility that prevails for all studied water levels and horizontal velocities. However, it has only a minor effect on the temperature distribution in the central rod bundle. By combining experiments and numerical simulations, the study provides important information about the decisive parameters for the heat exchange in a spent fuel pool in case of an accident with loss of cooling. The exposed length of the fuel rods is of much more importance than the magnitude of the horizontal velocity above the fuel assembly.     

Mathematical modeling of brush seals

A computational fluid dynamics (CFD)-based model of brush seals has been developed and tested against other workers' experimental data. In the model, the brush is treated as an axisymmetric, anisotropic porous region with nonlinear resistance coefficients. The resistance coefficients are chosen through calibration against measurements. The CFD model gives predictions of flow rate, pressure distribution, velocity field, and bending forces on the bristles. The bristle forces are used in a separate calculation to estimate bristle bending and reaction forces on the shaft and backing plate. Bending in both the axial direction and the orthogonal plane are considered.     

结构振动对湍流近尾迹的影响

研究了圆柱绕流中流体与结构的相互作用,侧重结构振动对湍流尾迹的影响,用激光测振仪测量圆柱在升力方向的位移;用热线和LDA（二维）测量湍流的近尾迹,通过变化自由流的速度和圆柱体直径（特征尺寸）来变化雷诺数,用两个振动特性不同的（一个相对刚性,一个相对弹）圆柱来产生尾迹,研究固体结构振动对湍流近尾迹的平均速度场和湍流场的影响,结果表明,结构自由振动对湍流近尾迹场影响明显,该影响随雷诺数的变化不明显。     

Experimental PIV/V3V measurements of vortex-induced fluid–structure interaction in turbulent flow—A new benchmark FSI-PfS-2a

The investigation of the bidirectional coupling between a fluid flow and a structure motion is a growing branch of research in science and industry. Applications of the so-called fluid–structure interactions (FSI) are widespread. To improve coupled numerical FSI simulations, generic experimental benchmark studies of the fluid and the structure are necessary. In this work, the coupling of a vortex-induced periodic deformation of a flexible structure mounted behind a rigid cylinder and a fully turbulent water flow performed at a Reynolds number of Re=30 470 is experimentally investigated with a planar particle image velocimetry (PIV) and a volumetric three-component velocimetry (V3V) system. To determine the structure displacements a multiple-point laser triangulation sensor is used. The three-dimensional fluid velocity results show shedding vortices behind the structure, which reaches the second swiveling mode with a frequency of about 11.2 Hz corresponding to a Strouhal number of St=0.177. Providing phase-averaged flow and structure measurements precise experimental data for coupled computational fluid dynamics (CFD) and computational structure dynamics (CSD) validations are available for this new benchmark case denoted FSI-PfS-2a. The test case possesses four main advantages: (i) the geometry is rather simple; (ii) kinematically, the rotation of the front cylinder is avoided; (iii) the boundary conditions are well defined; (iv) nevertheless, the resulting flow features and structure displacements are challenging from the computational point of view. In addition to the flow field and displacement data a PIV-based force calculation method is used to estimate the lift and drag coefficients of the moving structure.     

Particle size distribution in CPFD modeling of gas-solid flows in a CFB riser

A computational particle fluid dynamics(CPFD) numerical method to model gas-solid flows in a circulating fluidized bed(CFB) riser was used to assess the effects of particle size distribution(PSD) on solids distribution and flow.We investigated a binary PSD and a polydisperse PSD case.Our simulations were compared with measured solids concentrations and velocity profiles from experiments,as well as with a published Eulerian-Eulerian simulation.Overall flow patterns were similar for both simulation cases,as confirmed by experimental measurements.However,our fine-mesh CPFD simulations failed to predict a dense bottom region in the riser,as seen in other numerical studies.Above this bottom region,distributions of particle volume fraction and particle vertical velocity were consistent with our experiments,and the simulated average particle diameter decreased as a power function with riser height.Interactions between particles and walls also were successfully modeled,with accurate predictions for the lateral profiles of particle vertical velocity.It was easy to implement PSD into the CPFD numerical model,and it required fewer computational resources compared with other models,especially when particles with a polydisperse PSD were present in the heterogeneous flow.     

Particle size distribution in CPFD modeling of gas–solid flows in a CFB riser

A computational particle fluid dynamics (CPFD) numerical method to model gas–solid flows in a circulating fluidized bed (CFB) riser was used to assess the effects of particle size distribution (PSD) on solids distribution and flow. We investigated a binary PSD and a polydisperse PSD case. Our simulations were compared with measured solids concentrations and velocity profiles from experiments, as well as with a published Eulerian-Eulerian simulation. Overall flow patterns were similar for both simulation cases, as confirmed by experimental measurements. However, our fine-mesh CPFD simulations failed to predict a dense bottom region in the riser, as seen in other numerical studies. Above this bottom region, distributions of particle volume fraction and particle vertical velocity were consistent with our experiments, and the simulated average particle diameter decreased as a power function with riser height. Interactions between particles and walls also were successfully modeled, with accurate predictions for the lateral profiles of particle vertical velocity. It was easy to implement PSD into the CPFD numerical model, and it required fewer computational resources compared with other models, especially when particles with a polydisperse PSD were present in the heterogeneous flow.     

Time-resolved velocity and concentration measurements in variable-viscosity turbulent jet flow

This paper addresses the ability to reliably measure the fluctuating velocity field in variable-viscosity flows (herein, a propane–air mixture), using hot-wire anemometry. Because the latter is sensitive to both velocity and concentration fluctuations, the instantaneous concentration field also needs to be inferred experimentally. To overcome this difficulty, we show that the hot-wire response becomes insensitive to the concentration of the field, when a small amount of neon is added to the air. In this way, velocity measurements can be made independently of the concentration field. Although not necessary to velocity measurements, Rayleigh light-scattering technique is also used to infer the local (fluctuating) concentration, and, therefore, the viscosity of the fluid. Velocity and concentration measurements are performed in a turbulent propane jet discharging into an air–neon co-flow, for which the density and viscosity ratios are 1.52 and 1/5.5, respectively. The Reynolds number (based on injection diameter and velocity) is 15400. These measurements are first validated: the axial decay of the mean velocity and concentration, as well as the lateral mean and RMS profiles of velocity and concentration, is in full agreement with the existing literature. The variable-viscosity flow along the axis of the round jet is then characterized and compared with a turbulent air jet discharging into still air, for which the Reynolds number (based on injection diameter and velocity) is 5400. Both flows have the same initial jet momentum. As mixing with the viscous co-flow is enhanced with increasing downstream position, the viscosity of the fluid increases rapidly for the case of the propane jet. In comparison with the air jet, the propane jet exhibits: (1) a lower local Reynolds number based on the Taylor microscale (by a factor of four); (2) a reduced range of scales present in the flow; (3) the isotropic form of the mean energy dissipation rate is first more enhanced and then drastically diminishes and (4) a progressively increasing local Schmidt number (from 1.36 to 7.5) for increasing downstream positions. Therefore, the scalar spectra exhibit an increasingly prominent Batchelor regime with a ~ k ^?1 scaling law. The experimental technique developed herein provides a reliable method for the study of variable-viscosity flows.     

面向主给水系统的流动调整器整流性能数值研究

针对某核电站主给水系统低负荷时文丘里流量计测值虚高所提出的流场改善需求,利用计算流体力学方法结合SST k–ω湍流模型,研究安装流动调整器时的流场特征和整流性能。首先分别对流场未安装和安装栅格（AMCA）调整器及不同结构参数的孔板调整器进行三维建模和计算,获得不同调整器对流场速度、压力、流场均匀度及流量计流出系数的影响,对比不同调整器对旁路支管所致流速畸变的整流性能,并确定选用孔板调整器的构型。其次改变调整器安装位置,对比安装位置对流场速度分布、阻力损失、均匀度改善及流量计流出系数的影响。三是针对满负荷大流量工况,将调整器与不同开度的主调节阀耦合建模,计算调整器所致额外压力损失,分析阀门开度与整流器压力损失的关系。获得的调整器构型和安装位置等参数的影响规律可为核电站主给水系统解决同类问题时的流动调整器选型和安装提供参考。

     

表面粗糙度对微细管内气体流动特性的影响

采用了表面粗糙度粘性系数模型对微细管内的气体流动进行数值模拟，以研究微管内壁表面粗糙度对微管内气体流动的影响。运用本文改进的表面粗糙度粘性系数模型，数值模拟与实验数据十分吻合。计算结果表明，进出口压力一定时，表面粗糙度对流场的压力、密度及温度分布的影响不大，但是对速度场的影响十分显著，表面粗糙度使气体流动速度减小，并使壁面附近的速度梯度减小，从而使通过管道的气体质量流量减小，在微管内的气体流动中，表面粗糙度的影响是不能被忽略的。     

NUMERICAL MODELLING AND PARTICLE IMAGE VELOCIMETRY MEASUREMENT OF THE LAMINAR FLOW FIELD INDUCED BY AN ENCLOSED ROTATING DISC

The fluid flow field within an enclosed cylindrical chamber with a rotating flat disc was calculated using a finite volume computational fluid dynamics (CFD) model and compared with particle image velocimetry (PIV) measurements. Two particular laminar cases near the Transitional flow regime were investigated: Reynolds number Re=2.5×1 ⁴, chamber aspect ratio G (h/R_d)=0.2 and Re=4.2×10⁴, G (h/R_d)=0.217. This enabled direct comparison with the numerical and experimental results reported by other researchers. The computational details and some major factors that affect the computed accuracy and convergence speed are also discussed in detail. PIV results containing some 4300 velocity vector points in each of seven planes for each case were obtained from the flow field parallel to the rotating disc. It was found that PIV results could be obtained in planes within the boundary layers as well as the core flow by careful use of a thin laser illumination sheet and correct choice of laser pulse separation. There was close agreement between numerical results, the present PIV measurements and other reported experimental and numerical results.