Velocity Structure Functions Research Articles

A cryogenic-helium pipe flow facility with unique double-line molecular tagging velocimetry capability

Cryogenic helium-4 has extremely small kinetic viscosity, which makes it a promising material for high Reynolds (Re) number turbulence research in compact laboratory apparatus. In its superfluid phase (He II), helium has an extraordinary heat transfer capability and has been utilized in various scientific and engineering applications. In order to unlock the full potential of helium in turbulence research and to improve our understanding of the heat transfer mechanism in He II, a flow facility that allows quantitative study of helium heat-and-mass transfer processes is needed. Here, we report our work in assembling and testing a unique helium pipe-flow facility that incorporates a novel double-line molecular tagging velocimetry (DL-MTV) system. This flow facility allows us to generate turbulent pipe flows with Re above 107, and it can also be adapted to produce heat-induced counterflow in He II. The DL-MTV system, which is based on the generation and tracking of two parallel thin He2 * molecular tracer lines with an adjustable separation distance, allows us to measure not only the velocity profile but also both the transverse and longitudinal spatial velocity structure functions. We have also installed a differential pressure sensor on the flow pipe for pressure drop measurements. The testing results of the flow facility and the measuring instruments are presented. We discuss how this facility will allow us to solve some outstanding problems in the helium heat-and-mass transfer topic area.

Direct Detection of Black Hole-driven Turbulence in the Centers of Galaxy Clusters

Abstract Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this “feedback” process operates. We use high-resolution optical data to study the kinematics of multiphase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide range of scales in the centers of three nearby galaxy clusters: Perseus, A2597, and Virgo. We find that the motions of the filaments are turbulent in all three clusters studied. There is a clear correlation between features of the VSFs and the sizes of bubbles inflated by SMBH-driven jets. Our study demonstrates that SMBHs are the main driver of turbulent gas motions in the centers of relaxed galaxy clusters and suggests that this turbulence is an important channel for coupling feedback to the environment. Our measured amplitude of turbulence is in good agreement with Hitomi Doppler line broadening measurement and X-ray surface-brightness fluctuation analysis, suggesting that the motion of the cold filaments is well-coupled to that of the hot gas. The smallest scales that we probe are comparable to the mean free path in the intracluster medium. Our direct detection of turbulence on these scales provides the clearest evidence to date that isotropic viscosity is suppressed in the weakly collisional, magnetized intracluster plasma.

Projected velocity statistics of interstellar turbulence

ABSTRACT Velocity statistics is a direct probe of the dynamics of interstellar turbulence. Its observational measurements are very challenging due to the convolution between density and velocity and projection effects. We introduce the projected velocity structure function, which can be generally applied to statistical studies of both subsonic and supersonic turbulence in different interstellar phases. It recovers the turbulent velocity spectrum from the projected velocity field in different regimes, and when the thickness of a cloud is less than the driving scale of turbulence, it can also be used to determine the cloud thickness and the turbulence driving scale. By applying it to the existing core velocity dispersion measurements of the Taurus cloud, we find a transition from the Kolmogorov to the Burgers scaling of turbulent velocities with decreasing length-scales, corresponding to the large-scale solenoidal motions and small-scale compressive motions, respectively. The latter occupy a small fraction of the volume and can be selectively sampled by clusters of cores with the typical cluster size indicated by the transition scale.

Similarities between the structure functions of thermal convection and hydrodynamic turbulence

In this paper, we analyze the scaling of velocity structure functions of turbulent thermal convection. Using high-resolution numerical simulations, we show that the structure functions scale similar to those of hydrodynamic turbulence, with the scaling exponents in agreement with the predictions of She and Leveque [“Universal scaling laws in fully developed turbulence,” Phys. Rev. Lett. 72, 336–339 (1994)]. The probability distribution functions of velocity increments are non-Gaussian with wide tails in the dissipative scales and become close to Gaussian in the inertial range. The tails of the probability distribution follow a stretched exponential. We also show that in thermal convection, the energy flux in the inertial range is less than the viscous dissipation rate. This is unlike in hydrodynamic turbulence where the energy flux and the dissipation rate are equal.

Numerical simulation of two-dimensional turbulence based on Immersed Boundary Lattice Boltzmann method

The two-dimensional turbulence is numerically investigated using Immersed Boundary Lattice Boltzmann Method. The 2D turbulence should be considered as 2D channel flow where the flow is forced by the arrays of cylinders vertically placed in the inlet of 2D channel. The inverse cascade in two cases, the normal wall boundary case and rough wall boundary case, is obtained in the inertial range. The scaling behavior of energy spectrum in the inverse cascade is k−5/3, which is according with the Kraichnan theory of 2D turbulence. It is found that the time-evolving vortex number density distribution n(A)∼t−1A−1, based on the Rortex vortex definition criterion, which is in fair agreement with the theory of Burgess and Scott at intermediate scales sufficiently far from the forcing and the largest vortex scale. The energy flux of the normal wall boundary case cascades to the large scale. The energy flux of the rough wall boundary case cascades to large scale when the scale l < D. The energy flux at l > D cascades to the small scale. The energy dissipation filed εl coarse grained at the scale l has the pow-law behaviors with the scale. The intermittency measured by PDF exists in the inverse inertial range of 2D turbulence. On the basis of the vortex scaling law of Burgess and Scott and our simulation results, the universal scaling law of velocity structure function in the inverse cascade is obtained by extending the 3D turbulent S-L intermittency model to the inverse inertial range of 2D turbulence.

Circulation in High Reynolds Number Isotropic Turbulence is a Bifractal

The turbulence problem at the level of scaling exponents is hard in part because of the multifractal scaling of small scales, which demands that each moment order be treated and understood independently. This conclusion derives from studies of velocity structure functions, energy dissipation, enstrophy density (that is, square of vorticity), etc. However, it is likely that there exist other physically pertinent quantities with uncomplicated structure in the inertial range, potentially resulting in huge simplifications in the turbulence theory. We show that velocity circulation around closed loops is such a quantity. By using a large databases of isotropic turbulence, generated from numerical simulations of the Navier-Stokes equations over a wide range of Reynolds numbers, we show that circulation exhibits a bifractal behavior at the highest Reynolds number considered: space filling for moments up to order $3$ and a mono-fractal with an unchanging dimension of about $2.5$ for higher orders; this change in character roughly at the third-order moment is reminiscent of a "phase transition". We explore the possibility that circulation becomes effectively space filling at much higher Reynolds numbers even though it may technically be regarded as a bifractal. We confirm that the circulation properties depend on only the area of the loop, not its shape; and, for a figure-$8$ loop, the relevant area is the scalar sum of the two segments of the loop.

How do velocity structure functions trace gas dynamics in simulated molecular clouds?

Context. Supersonic disordered flows accompany the formation and evolution of molecular clouds (MCs). It has been argued that this is turbulence that can support against gravitational collapse and form hierarchical sub-structures. Aims. We examine the time evolution of simulated MCs to investigate: What physical process dominates the driving of turbulent flows? How can these flows be characterised? Are they consistent with uniform turbulence or gravitational collapse? Do the simulated flows agree with observations? Methods. We analysed three MCs that have formed self-consistently within kiloparsec-scale numerical simulations of the interstellar medium (ISM). The simulated ISM evolves under the influence of physical processes including self-gravity, stratification, magnetic fields, supernova-driven turbulence, and radiative heating and cooling. We characterise the flows using velocity structure functions (VSFs) with and without density weighting or a density cutoff, and computed in one or three dimensions. However, we do not include optical depth effects that can hide motions in the densest gas, limiting comparison of our results with observations. Results. In regions with sufficient resolution, the density-weighted VSFs initially appear to follow the expectations for uniform turbulence, with a first-order power-law exponent consistent with Larson’s size-velocity relationship. Supernova blast wave impacts on MCs produce short-lived coherent motions at large scales, increasing the scaling exponents for a crossing time. Gravitational contraction drives small-scale motions, producing scaling coefficients that drop or even turn negative as small scales become dominant. Removing the density weighting eliminates this effect as it emphasises the diffuse ISM. Conclusions. We conclude that two different effects coincidentally reproduce Larson’s size velocity relationship. Initially, uniform turbulence dominates, so the energy cascade produces VSFs that are consistent with Larson’s relationship. Later, contraction dominates and the density-weighted VSFs become much shallower or even inverted, but the relationship of the global average velocity dispersion of the MCs to their radius follows Larson’s relationship, reflecting virial equilibrium or free-fall collapse. The injection of energy by shocks is visible in the VSFs, but decays within a crossing time.

Towards mapping turbulence in the intra-cluster medium

Context. X-ray observations of galaxy clusters provide insights into the nature of gaseous turbulent motions, their physical scales, and the fundamental processes to which they are related. Spatially-resolved, high-resolution spectral measurements of X-ray emission lines provide diagnostics on the nature of turbulent motions in emitting atmospheres. Since they are acting on scales comparable to the size of the objects, the uncertainty on these physical parameters is limited by the number of observational measurements, through sample variance. Aims. We propose a different and complementary approach to repeating numerical simulations for the computation of sample variance (i.e. Monte-Carlo sampling) by introducing new analytical developments for lines diagnosis. Methods. We considered the model of a “turbulent gas cloud”, consisting in isotropic and uniform turbulence described by a universal Kolmogorov power-spectrum with random amplitudes and phases in an optically thin medium. Following a simple prescription for the four-term correlation of Fourier coefficients, we derived generic expressions for the sample mean and variance of line centroid shift, line broadening, and projected velocity structure function. We performed a numerical validation based on Monte-Carlo simulations for two popular models of gas emissivity based on the β-model. Results. Generic expressions for the sample variance of line centroid shifts and broadening in arbitrary apertures are derived and match the simulations within their range of applicability. Generic expressions for the mean and variance of the structure function are provided and verified against simulations. An application to the Athena/X-IFU (Advanced Telescope for High-ENergy Astrophysics/X-ray Integral Field Unit) and XRISM/Resolve (X-ray Imaging and Spectroscopy Mission) instruments forecasts the potential of sensitive, spatially-resolved spectroscopy to probe the inertial range of turbulent velocity cascades in a Coma-like galaxy cluster. Conclusions. The formulas provided are of generic relevance and can be implemented in forecasts for upcoming or current X-ray instrumentation and observing programmes.

Can We Detect Submesoscale Motions in Drifter Pair Dispersion?

AbstractA cluster of 45 drifters deployed in the Bay of Bengal is tracked for a period of four months. Pair dispersion statistics, from observed drifter trajectories and simulated trajectories based on surface geostrophic velocity, are analyzed as a function of drifter separation and time. Pair dispersion suggests nonlocal dynamics at submesoscales of 1–20 km, likely controlled by the energetic mesoscale eddies present during the observations. Second-order velocity structure functions and their Helmholtz decomposition, however, suggest local dispersion and divergent horizontal flow at scales below 20 km. This inconsistency cannot be explained by inertial oscillations alone, as has been reported in recent studies, and is likely related to other nondispersive processes that impact structure functions but do not enter pair dispersion statistics. At scales comparable to the deformation radius LD, which is approximately 60 km, we find dynamics in agreement with Richardson’s law and observe local dispersion in both pair dispersion statistics and second-order velocity structure functions.

Finite Reynolds number effect and the 4/5 law

Kolmogorov [A. N. Kolmogorov, Dokl. Akad. Nauk SSSR 30, 299 (1941)] formulated a theory of small-scale turbulence (K41), valid at extremely large Reynolds numbers, based on two similarity hypotheses and on an exact result derived from the transport equation for the second-order structure function, known as the 4/5 law. Although K41 was praised for its simplicity and elegance, Kolmogorov [A. N. Kolmogorov, J. Fluid Mech. 13, 82 (1962).] proposed a new refined similarity hypothesis (K62) mainly to account for the effect of the large scales on the small scales. It has been widely interpreted in the literature as a correction to K41 arising from the intermittency of the instantaneous energy dissipation rate ε. In this paper we argue that since K62 retains the 4/5 law, it must satisfy the same constraints as K41, viz., extremely large Reynolds number and flow stationarity. The retention of the 4/5 law is not however consistent with the presence of nonstationarity due to the effect of the large scales, as postulated by K62. A relatively extensive survey of published data shows that, indeed, the 4/5 law has not yet been observed in either experiments or simulations due to the Reynolds number not being sufficiently large. The use of the transport equation for the second-order structure function, together with an empirical model for the Kolmogorov-normalized second-order velocity structure function, confirms that the 4/5 law is established only after this structure function becomes independent of the Reynolds number.Received 25 April 2019DOI:https://doi.org/10.1103/PhysRevFluids.4.084602©2019 American Physical SocietyPhysics Subject Headings (PhySH)Research AreasTurbulenceFluid Dynamics

Characterization of Submesoscale Turbulence in the East/Japan Sea Using Geostationary Ocean Color Satellite Images

AbstractSubmesoscale processes are key in understanding physical and biological phenomena near the surface, but there remains a lack of observational evidence over large areas. We used hourly images from a geostationary satellite that can resolve variation in surface ocean color over an area of few hundred kilometers. The temporal variation in the surface chlorophyll a distribution captured by the satellite images was first used to generate a submesoscale‐permitting velocity field, from which we calculated the turbulence statistics such as kinetic energy spectra, velocity structure functions, and energy flux. Application to the April scenes in the East/Japan Sea showed that the kinetic energy spectra had a transition scale at 50 km that suggested two spectral regimes following k−3 and k−5/3, implying the coexistence of quasi‐geostrophic turbulence and surface quasi‐geostrophic turbulence. The chlorophyll a scalar spectrum suggested two spectral regimes of k−5/3 and k−1 with a transition at 3 km.

Submesoscale dispersion of surface drifters in a coastal sea near offshore wind farms

Abstract. We analysed relative dispersion of surface drifters released as pairs (six instances) or triplets (two instances) during three field experiments conducted in the German Bight in close proximity to wind farms. There is some tentative evidence that nearly exponential growth of relative dispersion (non-local dispersion) preferably occurs for drifter pairs that are most exposed to the influence of a wind farm. Kinetic energy spectra and velocity structure functions are analysed with regard to the assumption that turbulent energy could be injected by tides, possibly also via an interaction between tidal currents and wind turbine towers. Applicability of inertial range turbulence theory, however, can be doubted given distinct peaks of overtides observed in velocity power spectra. More comprehensive studies would be needed to better separate submesoscale effects of wind farms, tides and possibly baroclinic instabilities on observed drifter behaviour in a complex coastal environment.

Statistical Properties of the Velocity Field for the 3D Hydrodynamic Turbulence Onset

We study the statistical correlation functions for the three-dimensional hydrodynamic turbulence onset when the dynamics is dominated by the pancake-like high-vorticity structures. With extensive numerical simulations, we systematically examine the two-points structure functions (moments) of velocity. We observe formation of the power-law scaling for both the longitudinal and the transversal moments in the same interval of scales as for the energy spectrum. The scaling exponents for the velocity structure functions demonstrate the same key properties as for the stationary turbulence case. In particular, the exponents depend on the order of the moment non-trivially, indicating the intermittency and the anomalous scaling, and the longitudinal exponents turn out to be slightly larger than the transversal ones. When the energy spectrum has power-law scaling close to the Kolmogorov's one, the longitudinal third-order moment shows close to linear scaling with the distance, in line with the Kolmogorov's 4/5-law despite the strong anisotropy.

Artificial neural network approach to large-eddy simulation of compressible isotropic turbulence.

A subgrid-scale (SGS) model for large-eddy simulation (LES) of compressible isotropic turbulence is constructed by using a data-driven framework. An artificial neural network (ANN) based on local stencil geometry is employed to predict the unclosed SGS terms. The input features are based on the first-order and second-order derivatives of filtered velocity and temperature which appear in the second-order Taylor approximation of the SGS stress and heat flux. It is shown that the proposed ANN-7 model performs better than the gradient model in the a priori test. The correlation coefficient is larger and the relative error is smaller for ANN-7 model as compared to those of the gradient model in the a priori test. In an a posteriori analysis, the performance of ANN-7 model shows advantage over the dynamic Smagorinsky model and dynamic mixed model in the prediction of spectra and structure functions of velocity and temperature, and instantaneous flow structures. Artificial neural network is a promising tool for understanding the physical fundamentals of SGS unclosed terms with further improvement.

Cascade of turbulent energy and scalar variance in DC electrokinetic turbulence

Abstract Turbulent flow can be generated by DC electrokinetic (EK) mechanism resulted from the electric conductivity and permittivity variation in fluids, as has been demonstrated by Varshney et al (2016), where a −1.4 slope of velocity power spectrum is observed. Here, we theoretically investigate the scaling exponents of both velocity and scalar structure functions in DC EK turbulence. The scaling exponents of the second order structure functions of velocity and scalar in the so-called electric body force (EBF) dominant subrange are 2/5 (equivalent to the −7/5 slope of velocity power spectrum) and 4/5 respectively. The theory directly explains the experimental results of Varshney et al. (2016). Based on Karman–Howarth equation with forcing terms, the third order structure function of velocity which is the representative of energy transport process is investigated in the DC EK turbulence. Compared with the conventional hydrodynamic turbulence, the energy transport in DC EK turbulence has an additional scaling subrange and is governed by additional small length scales.

Resolved and subgrid dynamics of Rayleigh–Bénard convection

In this work we present and demonstrate the reliability of a theoretical framework for the study of thermally driven turbulence. It consists of scale-by-scale budget equations for the second-order velocity and temperature structure functions and their limiting cases, represented by the turbulent kinetic energy and temperature variance budgets. This framework represents an extension of the classical Kolmogorov and Yaglom equations to inhomogeneous and anisotropic flows, and allows for a novel assessment of the turbulent processes occurring at different scales and locations in the fluid domain. Two relevant characteristic scales, $\ell _{c}^{u}$ for the velocity field and $\ell _{c}^{\unicode[STIX]{x1D703}}$ for the temperature field, are identified. These variables separate the space of scales into a quasi-homogeneous range, characterized by turbulent kinetic energy and temperature variance cascades towards dissipation, and an inhomogeneity-dominated range, where the production and the transport in physical space are important. This theoretical framework is then extended to the context of large-eddy simulation to quantify the effect of a low-pass filtering operation on both resolved and subgrid dynamics of turbulent Rayleigh–Bénard convection. It consists of single-point and scale-by-scale budget equations for the filtered velocity and temperature fields. To evaluate the effect of the filter length $\ell _{F}$ on the resolved and subgrid dynamics, the velocity and temperature fields obtained from a direct numerical simulation are split into filtered and residual components using a spectral cutoff filter. It is found that when $\ell _{F}$ is smaller than the minimum values of the cross-over scales given by $\ell _{c,min}^{\unicode[STIX]{x1D703}\ast }=\ell _{c,min}^{\unicode[STIX]{x1D703}}Nu/H=0.8$, the resolved processes correspond to the exact ones, except for a depletion of viscous and thermal dissipations, and the only role of the subgrid scales is to drain turbulent kinetic energy and temperature variance to dissipate them. On the other hand, the resolved dynamics is much poorer in the near-wall region and the effects of the subgrid scales are more complex for filter lengths of the order of $\ell _{F}\approx 3\ell _{c,min}^{\unicode[STIX]{x1D703}}$ or larger. This study suggests that classic eddy-viscosity/diffusivity models employed in large-eddy simulation may suffer from some limitations for large filter lengths, and that alternative closures should be considered to account for the inhomogeneous processes at subgrid level. Moreover, the theoretical framework based on the filtered Kolmogorov and Yaglom equations may represent a valuable tool for future assessments of the subgrid-scale models.

About Universality and Thermodynamics of Turbulence.

This paper investigates the universality of the Eulerian velocity structure functions using velocity fields obtained from the stereoscopic particle image velocimetry (SPIV) technique in experiments and direct numerical simulations (DNS) of the Navier-Stokes equations. It shows that the numerical and experimental velocity structure functions up to order 9 follow a log-universality (Castaing et al. Phys. D Nonlinear Phenom. 1993); this leads to a collapse on a universal curve, when units including a logarithmic dependence on the Reynolds number are used. This paper then investigates the meaning and consequences of such log-universality, and shows that it is connected with the properties of a “multifractal free energy”, based on an analogy between multifractal and thermodynamics. It shows that in such a framework, the existence of a fluctuating dissipation scale is associated with a phase transition describing the relaminarisation of rough velocity fields with different Hölder exponents. Such a phase transition has been already observed using the Lagrangian velocity structure functions, but was so far believed to be out of reach for the Eulerian data.

An empirical expression for on the axis of a slightly heated turbulent round jet

Self-preservation analyses of the equations for the mean temperature and the second-order temperature structure function on the axis of a slightly heated turbulent round jet are exploited in an attempt to develop an analytical expression for$\unicode[STIX]{x1D716}_{\unicode[STIX]{x1D703}}$, the mean dissipation rate of$\overline{\unicode[STIX]{x1D703}^{2}}/2$, where$\overline{\unicode[STIX]{x1D703}^{2}}$is the temperature variance. The analytical approach follows that of Thiessetet al.(J. Fluid Mech., vol. 748, 2014, R2) who developed an expression for$\unicode[STIX]{x1D716}_{k}$, the mean turbulent kinetic energy dissipation rate, using the transport equation for$\overline{(\unicode[STIX]{x1D6FF}u)^{2}}$, the second-order velocity structure function. Experimental data show that complete self-preservation for all scales of motion is very well satisfied along the jet axis for streamwise distances larger than approximately 30 times the nozzle diameter. This validation of the analytical results is of particular interest as it provides justification and confidence in the analytical derivation of power laws representing the streamwise evolution of different physical quantities along the axis, such as:$\unicode[STIX]{x1D702}$,$\unicode[STIX]{x1D706}$,$\unicode[STIX]{x1D706}_{\unicode[STIX]{x1D703}}$,$R_{U}$,$R_{\unicode[STIX]{x1D6E9}}$(all representing characteristic length scales), the mean temperature excess$\unicode[STIX]{x1D6E9}_{0}$, the mixed velocity–temperature moments$\overline{u\unicode[STIX]{x1D703}^{2}}$,$\overline{v\unicode[STIX]{x1D703}^{2}}$and$\overline{\unicode[STIX]{x1D703}^{2}}$and$\unicode[STIX]{x1D716}_{\unicode[STIX]{x1D703}}$. Simple models are proposed for$\overline{u\unicode[STIX]{x1D703}^{2}}$and$\overline{v\unicode[STIX]{x1D703}^{2}}$in order to derive an analytical expression for$A_{\unicode[STIX]{x1D716}_{\unicode[STIX]{x1D703}}}$, the prefactor of the power law describing the streamwise evolution of$\unicode[STIX]{x1D716}_{\unicode[STIX]{x1D703}}$. Further, expressions are also derived for the turbulent Péclet number and the thermal-to-mechanical time scale ratio. These expressions involve global parameters that are most likely to be influenced by the initial and/or boundary conditions and are therefore expected to be flow dependent.

Surface Ocean Dispersion Observations From the Ship-Tethered Aerostat Remote Sensing System

Oil slicks and sheens reside at the air-sea interface, a region of the ocean that is notoriously difficult to measure and, therefore, little is known about the velocity field at the sea surface. The Ship-Tethered Aerostat Remote Sensing System (STARSS) was developed to measure Lagrangian velocities at the air-sea interface by tracking the transport and dispersion of hundreds of drift cards in the field of view of a high-resolution aerial imaging system. The camera had a field of view of approximately 300 m X 200 m and images were obtained every 15 seconds over periods of up to 3 hours during a series of experiments in the northern Gulf of Mexico in January-February 2016. STARSS was equipped with a GPS and inertial navigation system (INS) that was used to directly georectify the aerial images. A relative rectification technique was developed that translates and rotates the drift cards to minimize the total movement of all drift cards from one frame to the next. Rectified drift card positions were used to quantify scale-dependent dispersion by computing relative dispersion, relative diffusivity, and velocity structure functions. STARSS was part of a nested observational framework, which included deployments of large numbers of GPS-tracked surface drifters from two ships, in situ ocean measurements, and X-band radar observations of surface currents. STARSS operations were supported by weather forecasts from a high-resolution coupled atmosphere-wave-ocean model. Here we describe the STARSS system and image analysis techniques, and present results from an experiment that was conducted on a density front. To the best of our knowledge, these observations are the first of their kind and STARSS-like observations can be adopted into existing and planned oceanographic campaigns to produce a step-change in our understanding of small-scale and high-frequency variability at the air-sea interface.

The modulation of particles on coherent structure of turbulent boundary layer in dilute liquid-solid two-phase flow with PIV

Liquid-solidtwo-phase wall turbulent flow is an important flow pattern in engineering field. Studying the modulation of particles on fluid plays a significant role in revealing the mechanism of momentum, heat and mass transfer, improving production efficiency and decreasing energy consumption. The statistical characteristics and coherent structure in the flat plate turbulent boundary layer in clean water condition and near-neutrally liquid-solid two-phase flow with 200 um polystyrene particles (dp+=2.29,dp+=duτv) measured by means of PIV technique were compared and analyzed. The variation rule of turbulent statistic characteristics such as mean velocity profile, turbulence intensity and Reynolds shear stress of fluid with polystyrene particles and clean water under the same free flow velocity was investigated. Based on the locally-average velocity structure function and condition sampling technique, the coherent structure was detected, and the two-dimensional spatial topologies such as the fluctuation velocity, Reynolds shear stress and spanwise vorticity in ejection and sweep events at different wall heights under two working conditions were extracted with phase-average method. It was found that in the buffer layer the flow velocity with particles was greater than that of clean water, while there was little difference between them in the logarithmic law region. Compared with the clean water, the streamwise turbulent intensity was decreased in particle-laden flow, whereas the normal turbulent intensity, the comprehensive turbulence intensity and the Reynolds shear stress were all increased, with the expense of a slightly increase of the friction drag in the fluid with particles by 3.57%. In the bursting of coherent structure, the addition of polystyrene particles suppressed the fluid streamwise fluctuation in a turbulent boundary layer, gave a rise to the normal fluctuation and Reynolds shear stress, while inhibited the streamwise development of the spanwise vortex and promoted its normal development. The addition of particles facilitated the fluctuation of the fluid microcluster, and as well as the transverse momentum transfer in the log law region of turbulent boundary layer. The particles could increase the probability of the coherent structure ejection events and had less influence on the sweep events. The number and intensity of the ejection event was the key to improve the flow normal turbulence intensity and even to improve the energy exchange near the wall. The presence of appropriate scale particles in wall turbulence could play an important role in strengthening the transfer of momentum, heat and mass, and promoting mixing and reaction in the engineering, only by slightly increasing the fraction resistance.