Graphite dust resuspension in an HTR-10 steam generator

Graphite dust has an important effect on the safety of high-temperature gas-cooled reactors(HTR).The flow field in the steam generator was studied by the computational fluid dynamics(CFD) method,with the results indicating that the friction velocity in the windward and the leeward of the heat transfer tubes is relatively low and is higher at the sides.Further analysis of the resuspension of graphite dust indicates that the resuspension fraction reaches nearly zero for particles with a diameter less than 1 μm,whereas it will increases as the helium velocity in the steam generator increases for particle size larger than 1 μm.Moreover,the resuspension fraction increases as the particle size increases.The results also indicate that resuspension of the particles with sizes larger than 1 μm exhibited obvious differences in different parts of the steam generator.     

Motion of a spherical particle in a fluid forced vortex by DQM and DTM

In this study,coupled equations of the motion of a particle in a fluid forced vortex were investigated using the differential transformation method(DTM) with the Pade approximation and the differential quadrature method(DQM).The significant contribution of the work is the introduction of two new,fast and efficient solutions for a spherical particle in a forced vortex that are improvements over the previous numerical results in the literature.These methods represent approximations with a high degree of accuracy and minimal computational effort for studying the particle motion in a fluid forced vortex.In addition,the velocity profiles(angular and radial) and the position trajectory of a particle in a fluid forced vortex are described in the current study.     

Numerical simulation and PIV study of compressible vortex ring evolution

Formation and evolution of a compressible vortex ring generated at the open end of a short driver section shock tube has been simulated numerically for pressure ratios (PR) of 3 and 7 in the present study. Numerical study of compressible vortex rings is essential to understand the complicated flow structure and acoustic characteristics of many high Mach number impulsive jets where simultaneously velocity, density and pressure fields are needed. The flow development, incident shock formation, shock diffraction, vortex ring formation and its evolution are simulated using the AUSM+ scheme. The main focus of the present study is to evaluate the time resolved vorticity field of the vortex ring and the shock/expansion waves in the starting jet for short driver section shock tubes—a scenario where little data are available in existing literature. An embedded shock and a vortex induced shock are observed for PR =  7. However the vortex ring remains shock free, compact and unaffected by the trailing jet for PR =  3. Numerical shadowgraph shows the evolution of embedded shock and shock/expansion waves along with their interactions. The velocity and vorticity fields obtained from simulation are validated with the particle image velocimetry results and these data match closely. The translational velocity of the vortex ring, velocity across the vortex and the centre line velocity of the jet obtained from simulation also agree well with the experimental results.     

A three‐dimensional vortex particle‐in‐cell method for vortex motions in the vicinity of a wall

A new vortex particle‐in‐cell method for the simulation of three‐dimensional unsteady incompressible viscous flow is presented. The projection of the vortex strengths onto the mesh is based on volume interpolation. The convection of vorticity is treated as a Lagrangian move operation but one where the velocity of each particle is interpolated from an Eulerian mesh solution of velocity–Poisson equations. The change in vorticity due to diffusion is also computed on the Eulerian mesh and projected back to the particles. Where diffusive fluxes cause vorticity to enter a cell not already containing any particles new particles are created. The surface vorticity and the cancellation of tangential velocity at the plate are related by the Neumann conditions. The basic framework for implementation of the procedure is also introduced where the solution update comprises a sequence of two fractional steps. The method is applied to a problem where an unsteady boundary layer develops under the impact of a vortex ring and comparison is made with the experimental and numerical literature. Copyright © 2001 John Wiley & Sons, Ltd.     

Unsteady flow calculation in a radial flow centrifugal pump with spiral casing

The prediction of the two-dimensional unsteady flow established in a radial flow centrifugal pump is considered. Assuming the fluid incompressible and inviscid, the velocity field is represented by means of source and vorticity surface distributions as well as a set of point vortices. Using this representation, a grid-free (Lagrangian) numerical method is derived based on the coupling of the boundary element and vortex particle methods. In this context the source and vorticity surface distributions are determined through the non-entry boundary condition together with the unsteady Kutta condition. In order to satisfy Kelvin's theorem, vorticity is shed at the trailing edges of the impeller blades. Then the vortex particle method is used to approximate the convection of the free vorticity distribution. Results are given for a pump configuration experimentally tested by Centre Technique des Industries Mécaniques (CETIM). Comparisons between predictions and experimental data show the capability of the proposed method to reproduce the main features of the flow considered.     

Vibrations of a cylinder in a uniform flow in the presence of a no-slip side-wall

A circular cylinder placed in a uniform flow, and that spans the entire length between two side walls, may experience either parallel or oblique vortex shedding depending on the end conditions. It was shown by Mittal and Sidharth (2014) that the spatio-temporal periodicity of the oblique vortex shedding results in constant-in-time force experienced by the cylinder. On the contrary, parallel vortex shedding leads to fluid force that fluctuates with time. The free vibrations of a circular cylinder, in the presence of a wall, are investigated. For comparison, computations with end walls, where a slip condition on velocity is specified, are also carried out. The Reynolds number, based on the diameter of the cylinder and free-stream speed of the flow, is Re=100. The initial condition for the free vibrations is the fully developed unsteady flow past a stationary cylinder with oblique shedding. It is found that as the amplitude of vibration of the cylinder builds up, the vortices shed from the cylinder align with its axis leading to parallel shedding. The response of the cylinder is associated with two branches: initial and lower. On the lower branch, the response of the cylinder is virtually identical from two- and three-dimensional computations. The flow as well as the response is different on the initial branch and outside the synchronization regime. Forced vibrations confirm the phenomena.     

Accuracy of a deterministic particle method for Navier-Stokes equations

The accuracy of a deterministic particle method in approximating the solution of the Navier-Stokes equations is investigated. The convective part is solved using a classical vortex method for inviscid fluids, and an iterative procedure is added to improve the interpolation of the vorticity function. In our examples the vorticity is radially symmetric. For a regular initial data, a discrete quadratic error on the velocity and the vorticity is considered. Otherwise, for a singular initial data, the exact and computed angular moments of the vorticity are compared.     

A combined vortex and panel method for numerical simulations of viscous flows: a comparative study of a vortex particle method and a finite volume method

This paper describes and compares two vorticity‐based integral approaches for the solution of the incompressible Navier–Stokes equations. Either a Lagrangian vortex particle method or an Eulerian finite volume scheme is implemented to solve the vorticity transport equation with a vorticity boundary condition. The Biot–Savart integral is used to compute the velocity field from a vorticity distribution over a fluid domain. The vorticity boundary condition is improved by the use of an iteration scheme connected with the well‐established panel method. In the early stages of development of flows around an impulsively started circular cylinder, and past an impulsively started foil with varying angles of attack, the computational results obtained by the Lagrangian vortex method are compared with those obtained by the Eulerian finite volume method. The comparison is performed separately for the pressure fields as well. The results obtained by the two methods are in good agreement, and give a better understanding of the vorticity‐based methods. Copyright © 2005 John Wiley & Sons, Ltd.     

Direct numerical simulation of a particle-laden low Reynolds number turbulent round jet

Direct numerical simulation method is used for the investigating of particle-laden turbulent flows in a spatially evolution of low Reynolds number axisymmetric jet, and the Eulerian–Lagrangian point-particle approach is employed in the simulation. The simulation uses an explicit coupling scheme between particles and the fluid, which considers two-way coupling between the particle and the fluid. The DNS results are compared well with experimental data with equal Reynolds number (R_e = 1700). Our objects are: (i) to investigate the correlation between the particle number density and the fluctuating of fluid streamwise velocity; (ii) to examine whether the three-dimensional vortex structures in the particle-laden jet are the same as that in the free-air jet and how the particles modulate the thee-dimensional vortex structures and turbulence properties with different Stokes number particles; (iii) to discover the particle circumferential dispersion with different Stokes number particles. Our findings: (i) all the particles, regardless of their particle size, tend to preferentially accumulate in the region with large-than-mean fluid streamwise velocity; (ii) the small Stokes number particles take an important part in the modulation of three-dimensional vortex structures, but for the intermediate and larger sized particles, this modulation effect seems not so apparent; (iii) the particle circumferential dispersion is more effective for the smaller and intermediate sized particles, especially for the intermediate sized particles.     

Experimental Investigation of Three-Dimensional Vortex Structures Downstream of Vortex Generators Over a Backward-Facing Step

An experimental investigation of vortex generators has been carried out in turbulent backward-facing step (BFS) flow. The Reynolds number, based on a freestream velocity U₀ = 10 m/s and a step height h = 30 mm, was Re_h = 2.0 × 10⁴. Low-profile wedge-type vortex generators (VGs) were implemented on the horizontal surface upstream of the step. High-resolution planar particle image velocimetry (2D-2C PIV) was used to measure the separated shear layer, recirculation region and reattachment area downstream of the BFS in a single field of view. Besides, time-resolved tomographic particle image velocimetry (TR-Tomo-PIV) was also employed to measure the flow flied of the turbulent shear layer downstream of the BFS within a three-dimensional volume of 50 × 50 × 10 mm³ at a sampling frequency of 1 kHz. The flow control result shows that time-averaged reattachment length downstream of the BFS is reduced by 29.1 % due to the application of the VGs. Meanwhile, the Reynolds shear stress downstream of the VGs is considerably increased. Proper Orthogonal Decomposition (POD) and Dynamic Mode Decomposition (DMD) have been applied to the 3D velocity vector fields to analyze the complex vortex structures in the spatial and temporal approaches, respectively. A coherent bandwidth of Strouhal number 0.3 < St_h < 0.6 is found in the VG-induced vortices, and moreover, Λ-shaped three-dimensional vortex structures at St_h = 0.37 are revealed in the energy and dynamic approaches complementarily.     

Numerical study on heavy rigid particle motion in a plane wake flow by spectral element method

Experimental particle dispersion patterns in a plane wake flow at a high Reynolds number have been predicted numerically by discrete vortex method (Phys. Fluids A 1992; 4 :2244–2251; Int. J. Multiphase Flow 2000; 26 :1583–1607). To address the particle motion at a moderate Reynolds number, spectral element method is employed to provide an instantaneous wake flow field for particle dynamics equations, which are solved to make a detail classification of the patterns in relation to the Stokes and Froude numbers. It is found that particle motion features only depend on the Stokes number at a high Froude number and depend on both numbers at a low Froude number. A ratio of the Stokes number to squared Froude number is introduced and threshold values of this parameter are evaluated that delineate the different regions of particle behavior. The parameter describes approximately the gravitational settling velocity divided by the characteristic velocity of wake flow. In order to present effects of particle density but preserve rigid sphere, hollow sphere particle dynamics in the plane wake flow is investigated. The evolution of hollow particle motion patterns for the increase of equivalent particle density corresponds to that of solid particle motion patterns for the decrease of particle size. Although the thresholds change a little, the parameter can still make a good qualitative classification of particle motion patterns as the inner diameter changes. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical study of an inviscid incompressible flow through a channel of finite length

A two‐dimensional inviscid incompressible flow in a rectilinear channel of finite length is studied numerically. Both the normal velocity and the vorticity are given at the inlet, and only the normal velocity is specified at the outlet. The flow is described in terms of the stream function and vorticity. To solve the unsteady problem numerically, we propose a version of the vortex particle method. The vorticity field is approximated using its values at a set of fluid particles. A pseudo‐symplectic integrator is employed to solve the system of ordinary differential equations governing the motion of fluid particles. The stream function is computed using the Galerkin method. Unsteady flows developing from an initial perturbation in the form of an elliptical patch of vorticity are calculated for various values of the volume flux of fluid through the channel. It is shown that if the flux of fluid is large, the initial vortex patch is washed out of the channel, and when the flux is reduced, the initial perturbation evolves to a steady flow with stagnation regions. Copyright © 2008 John Wiley & Sons, Ltd.     

Experimental investigation of the interaction between a plane wall jet and a parallel offset jet

The dual-jet flow generated by a plane wall jet and a parallel offset jet at an offset ratio of d/w = 1.0 has been investigated using Particle Image Velocimetry (PIV). The particle images are captured, processed, and subsequently used to characterize the flow in terms of the 2D velocity and vorticity distributions. Statistical characteristics of the flow are obtained through ensemble averaging of 360 instantaneous velocity fields. Also presented is a time series of instantaneous flow fields to illustrate the dynamic interaction between the two jets. Results reveal that the near field of the flow is characterized by a periodic large-scale Karman-like vortex shedding similar to what would be expected in the wake of a bluff body. The existence of the Karman-like vortices results in periodic interactions between the two jets; in addition, these vortices produce noticeable impact on the jet outer layers, i.e., the free shear layer of the offset jet and the wall boundary layer of the wall jet. A schematic of vortex/shear layer interaction is proposed to illustrate the flow pattern.     

封闭水平矩形腔内流体自然对流第一次逆叉形分岔的数值研究

本文采用二阶全展开ETG(Euler-Taylor-Galerkin)分裂步有限元方法,对长宽比为3.5(L/B=3.5,如图 1所示)的封闭矩形腔体内,三种Pr数条件下,定常层流范围内,流体自然对流叉形分岔随Rayleigh数的演化过程进行了数值模拟。研究结果表明,该矩形腔内对应三种Pr数条件下,流体的叉形分岔的演化过程中,在第二次模态Ⅱ型叉形分岔之后,均会出现两个较小尺度涡旋合并,突变为一个较大尺度涡旋的全新叉形分岔模态。即在某临界Ra数两侧,存在定常四涡结构和定常三涡结构两个定常解支,当系统控制参数Ra越过临界值,前者被后者突发性取代,这是完全不同于传统叉形分岔的逆叉形分岔。其数值预报,则采用分半法结合流动拓扑结构及典型截面处速度扩线上鞍点的变化来确定。计算结果表明,在计算的Pr数条件下,随Pr数的增加逆叉形分岔对应临界Ra数的取值也会提高。     

Predicting transport by Lagrangian coherent structures with a high-order method

Recent developments in identifying Lagrangian coherent structures from finite-time velocity data have provided a theoretical basis for understanding chaotic transport in general flows with aperiodic dependence on time. As these theoretical developments are extended and applied to more complex flows, an accurate and general numerical method for computing these structures is needed to exploit these ideas for engineering applications. We present an unstructured high-order hp/spectral-element method for solving the two-dimensional compressible form of the Navier–Stokes equations. A corresponding high-order particle tracking method is also developed for extracting the Lagrangian coherent structures from the numerically computed velocity fields. Two different techniques are used; the first computes the direct Lyapunov exponent from an unstructured initial particle distribution, providing easier resolution of structures located close to physical boundaries, whereas the second advects a small material line initialized close to a Lagrangian saddle point to delineate these structures. We demonstrate our algorithm on simulations of a bluff-body flow at a Reynolds number of Re = 150 and a Mach number of M = 0.2 with and without flow forcing. We show that, in the unforced flow, periodic vortex shedding is predicted by our numerical simulations that is in stark contrast to the aperiodic flow field in the case with forcing. An analysis of the Lagrangian structures reveals a transport barrier that inhibits cross-wake transport in the unforced flow. The transport barrier is broken with forcing, producing enhanced transport properties by chaotic advection and consequently improved mixing of advected scalars within the wake.     

The effect of air-water interface on the vortex shedding from a vertical circular cylinder

The flow past an interface piercing circular cylinder at the Reynolds number Re=2.7×10⁴ and the Froude numbers Fr=0.2 and 0.8 is investigated using large-eddy simulation. A Lagrangian dynamic subgrid-scale model and a level set based sharp interface method are used for the spatially filtered turbulence closure and the air-water interface treatment, respectively. The mean interface elevation and the rms of interface fluctuations from the simulation are in excellent agreement with the available experimental data. The organized periodic vortex shedding observed in the deep flow is attenuated and replaced by small-scale vortices at the interface. The streamwise vorticity and the outward transverse velocity generated near the edge of the separated region, which enforces the separated shear layers to deviate from each other and restrains their interaction, are primarily responsible for the devitalization of the periodic vortex shedding at the interface. The lateral gradient of the difference between the vertical and transverse Reynolds normal stresses, increasing with the Froude number, is the main source of the streamwise vorticity and the outward transverse velocity at the interface.     

Turbulent flow around a wall-mounted cube: A direct numerical simulation

An immersed-boundary method was employed to perform a direct numerical simulation (DNS) of flow around a wall-mounted cube in a fully developed turbulent channel for a Reynolds number Re = 5610, based on the bulk velocity and the channel height. Instantaneous results of the DNS of a plain channel flow were used as a fully developed inflow condition for the main channel. The results confirm the unsteadiness of the considered flow caused by the unstable interaction of a horseshoe vortex formed in front of the cube and on both its sides with an arch-type vortex behind the cube. The time-averaged data of the turbulence mean-square intensities, Reynolds shear stresses, kinetic energy and dissipation rate are presented. The negative turbulence production is predicted in the region in front of the cube where the main horseshoe vortex originates.     

Experimental investigation of the wall shear stress and the vortex dynamics in a circular impinging jet

The wall shear stress and the vortex dynamics in a circular impinging jet are investigated experimentally for Re = 1,260 and 2,450. The wall shear stress is obtained at different radial locations from the stagnation point using the polarographic method. The velocity field is given from the time resolved particle image velocimetry (TR‐PIV) technique in both the free jet region and near the wall in the impinging region. The distribution of the momentum thickness is also inspected from the jet exit toward the impinged wall. It is found that the wall shear stress is correlated with the large-scale vortex passing. Both the primary vortices and the secondary structures strongly affect the variation of the wall shear stress. The maximum mean wall shear stress is obtained just upstream from the secondary vortex generation where the primary structures impinge the wall. Spectral analysis and cross-correlations between the wall shear stress fluctuations show that the vortex passing influences the wall shear stress at different locations simultaneously. Analysis of cross-correlations between temporal fluctuations of the wall shear stress and the transverse vorticity brings out the role of different vortical structures on the wall shear stress distribution for the two Reynolds numbers.     

Simulation of a Viscous Flow Around an Oscillating Cylinder at Low Keulegan-Carpenter Numbers

A finite difference simulation method for a viscous flow around a circular cylinder sinusoidally oscillating at low Keulegan-Carpenter numbers is presented. Navier-Stokes equations in finite difference form are solved on a moving grid system, based on a time dependent coordinate transformation. Evolution with time of the flow structures induced by a circular cylinder performing sinusoidal oscillations in a fluid at rest, by means of stream lines, pressure contours and vortex shedding is studied in detail at Keulegan-Carpenter numbers, Kc = 9.4 and 14. The time dependent drag and lift are also explained.     

Controlled preparation and characterization of nano-sized hexagonal Mg（OH）2 flame retardant

Nano-sized hexagonal magnesium hydroxide （Mg（OH）2） with good dispersibility was synthesized by a double injection-hydrothermal method, utilizing polyvinylpyrrolidone （PVP） as an additive and with optimized processing parameters. SEM and BET analysis showed that the mean particle size and specific surface area of the Mg（OH）2 particles were 174 nm and 50.77 m^2/g, respectively. The FT-IR spectra and the XRD patterns showed that PVP was adsorbed on the surface of the Mg（OH）2 crystal, thus effectively limiting particle agglomeration and hindering crystal growth along the （1 01 ） plane. TGA showed a decrease in the decomposition temperature and an increase in the weight loss of the Mg（OH）2 particles due to addition of PV/.