Source author record

Shiyi Chen

Shiyi Chen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

Enhancing large-scale motions and turbulent transport in rotating plane Poiseuille flow

Based on the properties of large-scale plume currents, an injection/suction control strategy is introduced to enhance the strength of plume currents and the turbulent transport of passive scalar in rotating plane Poiseuille flow. The control keeps the classical non-penetrative condition on the stable side, and prescribes spanwise varying wall-normal velocity on the unstable side. Direct numerical simulations show that at medium rotation numbers, very weak injection/suction (the root-mean-square velocity on the controlled side is below $1\%$ of the bulk mean velocity) with properly distributed slots could significantly enhance the large-scale plume currents which have the same spatial distribution as the injection slots. In addition, if the distance between plume currents are increased, the turbulent transport efficiency can be promoted. Therefore, taking the strength and transport efficiency into account, the distance between injection slots should be slightly larger than the intrinsic distance between plume currents to achieve maximum enhancement of turbulent transport. For comparison, the control with non-penetrative condition on the unstable side and injection/suction on the stable side is also examined with simulations. It is found that the injection/suction on the stable side has much weaker influence compared with that on the unstable side. This is because that plumes are generated on the unstable side, and that controlling the origin of plumes is more effective.

preprint2022arXiv

Rarefaction effects in head-on collision of two identical droplets

The head-on collision of two identical droplets is investigated based on the BGK-Boltzmann equation. Gauss-Hermite quadratures with different degree of precision are used to solve the kinetic equation, so that the continuum (solution truncated at the Navier-Stokes order) and non-continuum (rarefied gas dynamics) solutions can be compared. When the kinetic equation is solved with adequate accuracy, prominent variations of the vertical velocity (the collision is in the horizontal direction), the viscous stress components, and droplet morphology are observed during the formation of liquid bridge, which demonstrates the importance of the rarefaction effects and the failure of the Navier-Stokes equation. The rarefaction effects change the topology of streamlines near the droplet surface, suppress the high-magnitude vorticity concentration inside the interdroplet region, and promote the vorticity diffusion around outer droplet surface. Two physical mechanisms responsible for the local energy conversion between the free and kinetic energies are identified, namely, the total pressure-dilatation coupling effect and the interaction between the density gradient and strain rate tensor. An energy conversion analysis is performed to show that the rarefaction effects can enhance the conversion from free energy to kinetic energy and facilitate the discharge of interdroplet gas film along the vertical direction, thereby boosting droplet coalescence.

preprint2022arXiv

Study of a droplet breakup process in decaying homogeneous isotropic turbulence based on the phase-field DUGKS approach

The breakup of a spherical droplet in a decaying homogeneous isotropic turbulence is studied by solving the Cahn-Hilliard-Navier-Stokes equations, using the discrete unified gas kinetic scheme combined with the free-energy-based phase-field model. We focus on the combined effects of turbulence and surface tension on the breakup process by assuming that the two fluid phases have the same density and same viscosity. The key physical parameters of the system include the volume fraction (6.54%), the initial Weber number (21.7), and the initial Taylor microscale Reynolds number (58). Three distinct stages of droplet evolution are identified, namely, the deformation stage when the initially spherical droplet evolves into an irregular geometric shape with complex structures, the breakup stage when many daughter droplets are formed, and the restoration stage when the droplets relax towards spherical shape. These three stages are analyzed systematically from several perspectives: (1) a geometric perspective concerning the maximum equivalent diameter, the total number of droplets, total interface area, and probability distribution of droplet diameters, (2) a dynamic perspective concerning the evolution of local velocity and vorticity at the fluid-fluid interface, (3) a global perspective concerning the evolution of average kinetic energy / dissipation rate and their Fourier spectra, (4) spherical harmonics based energetics concerning simultaneous transfer of kinetic energy across different length scales and different radii relative to initial droplet center, and (5) the time evolution of global kinetic energy and free energy of the system.

preprint2020arXiv

Heterogeneous-Temporal Graph Convolutional Networks: Make the Community Detection Much Better

Community detection has long been an important yet challenging task to analyze complex networks with a focus on detecting topological structures of graph data. Essentially, real-world graph data contains various features, node and edge types which dynamically vary over time, and this invalidates most existing community detection approaches. To cope with these issues, this paper proposes the heterogeneous-temporal graph convolutional networks (HTGCN) to detect communities from hetergeneous and temporal graphs. Particularly, we first design a heterogeneous GCN component to acquire feature representations for each heterogeneous graph at each time step. Then, a residual compressed aggregation component is proposed to represent "dynamic" features for "varying" communities, which are then aggregated with "static" features extracted from current graph. Extensive experiments are evaluated on two real-world datasets, i.e., DBLP and IMDB. The promising results demonstrate that the proposed HTGCN is superior to both benchmark and the state-of-the-art approaches, e.g., GCN, GAT, GNN, LGNN, HAN and STAR, with respect to a number of evaluation criteria.

preprint2020arXiv

Obtaining the mean fields with known Reynolds stresses at steady state

With the rising of modern data science, data--driven turbulence modeling with the aid of machine learning algorithms is becoming a new promising field. Many approaches are able to achieve better Reynolds stress prediction, with much lower modeling error ($ε_M$), than traditional RANS models but they still suffer from numerical error and stability issues when the mean velocity fields are estimated using RANS equations with the predicted Reynolds stresses, illustrating that the error of solving the RANS equations ($ε_P$) is also very important. In the present work, the error $ε_P$ is studied separately by using the Reynolds stresses obtained from direct numerical simulation and we derive the sources of $ε_P$. For the implementations with known Reynolds stresses solely, we suggest to run an adjoint RANS simulation to make first guess on $ν_t^*$ and $S_{ij}^0$. With around 10 iterations, the error could be reduced by about one-order of magnitude in flow over periodic hills. The present work not only provides one robust approach to minimize $ε_P$, which may be very useful for the data-driven turbulence models, but also shows the importance of the nonlinear part of the Reynolds stresses in flow problems with flow separations.

preprint2019arXiv

Dual channels of helicity cascade in turbulent flows

Helicity, as one of only two inviscid invariants in three-dimensional turbulence, plays an important role in the generation and evolution of turbulence. From the traditional viewpoint, there exists only one channel of helicity cascade similar to that of kinetic energy cascade. Through theoretical analysis, we find that there are two channels in helicity cascade process. The first channel mainly originates from vortex twisting process, and the second channel mainly originates from vortex stretching process. By analysing the data of direct numerical simulations of typical turbulent flows, we find that these two channels behave differently. The ensemble averages of helicity flux in different channels are equal in homogeneous and isotropic turbulence, while they are different in other type of turbulent flows. The second channel is more intermittent and acts more like a scalar, especially on small scales. Besides, we find a novel mechanism of hindered even inverse energy cascade, which could be attributed to the second-channel helicity flux with large amplitude.

preprint2016arXiv

Aerodynamic Heating in Hypersonic Boundary Layers:\ Role of Dilatational Waves

The evolution of multi-mode instabilities in a hypersonic boundary layer and their effects on aerodynamic heating are investigated. Experiments are conducted in a Mach 6 wind tunnel using Rayleigh-scattering flow visualization, fast-response pressure sensors, fluorescent temperature-sensitive paint (TSP), and particle image velocimetry (PIV). Calculations are also performed based on both parabolized stability equations (PSE) and direct numerical simulations (DNS). It is found that second-mode dilatational waves, accompanied by high-frequency alternating fluid compression and expansion, produce intense aerodynamic heating in a small region that rapidly heats the fluid passing through it. As a result, the surface temperature rapidly increases and results in an overshoot over the nominal transitional value. When the dilatation waves decay downstream, the surface temperature decreases gradually until transition is completed. A theoretical analysis is provided to interpret the temperature distribution affected by the aerodynamic heating.

preprint2016arXiv

Evolution of material surfaces in the temporal transition in channel flow

We report a Lagrangian study on the evolution of material surfaces in the K-type temporal transitional channel flow. Based on the Eulerian velocity field from the DNS, a backward-particle-tracking method is applied to solve the transport equation of the Lagrangian scalar field, and then the iso-surfaces of the Lagrangian field can be extracted as material surfaces in the evolution. Three critical issues for Lagrangian investigations on the evolution of coherent structures using material surfaces are addressed. First, the initial scalar field is uniquely determined based on proposed criteria, so that the initial material surfaces can be approximated as vortex surfaces, and keep invariant in the initial laminar state. Second, the evolution of typical material surfaces initially from different wall distances is presented, and then the influential material surface with the maximum deformation is identified. Large vorticity variations with the maximum curvature growth of vortex lines are also observed on this surface. Moreover, crucial events in the transition can be characterized in a Lagrangian approach by conditional statistics on the material surfaces. Finally, the influential material surface, which is initially a vortex surface, is demonstrated as a surrogate of the vortex surface before significant topological changes of vortical structures. Therefore, this material surface can be used to elucidate the continuous evolution of vortical structures in transitional wall-bounded flows in a Lagrangian perspective. The evolution of the influential material surface is divided into three stages: the formation of a triangular bulge from an initially disturbed streamwise-spanwise sheet, rolling up of the vortex sheet near the bulge ridges with the vorticity intensification, and the generation and evolution of signature hairpin-like structures with self-induced dynamics of vortex filaments.

preprint2015arXiv

A numerical investigation on active and passive scalars in isotropic compressible turbulence

In this paper, we investigated the statistical differences between active and passive scalars in isotropic compressible turbulence. In the inertial range, the kinetic energy and scalar have the Kolmogorov spectra, and the Kolmogorov and Obukhov-Corrsin constants are 2.06 and 0.87, respectively. The local scaling exponents computed from the second-order structure functions exist plateaus for velocity and active scalar, while that for passive scalar takes first a minimum of 0.61 then a maximum of 0.73. The mixed third-order structure function of velocity and passive scalar satisfies the 4/3-law. For the scaling exponent, the one of velocity and passive scalar mixing is between those of velocity and passive scalar, while the one of velocity and active scalar mixing is below those of velocity and active scalar. At large amplitudes, the p.d.f. of active scalar fluctuations is super-Gaussian, whereas that of passive scalar fluctuations is sub-Gaussian. Moreover, the p.d.f.s of the two scalar increments are concave and convex shapes, respectively, which exhibit strong intermittency at small scales, and approach Gaussian as scale increases. The active scalar has "ramp-cliff" structures, while the passive scalar seems to be dominated by rarefaction and compression. By employing a "coarse-graining" approach, we study scalar cascades. Unlike passive scalar, the cascade of active scalar is mainly determined by the viscous dissipation at small scales and the pressure-dilatation at large scales, where the latter is substantial in the vicinity of small-scale shocklets but is negligible after space averages, because of the cancelations between rarefaction and compression regions. Finally, in the inertial range, the p.d.f.s for the subgrid-scale fluxes of scalars can collapse to the same distribution, revealing the scale-invariant feature for the statistics of active and passive scalars.

preprint2015arXiv

Effects of shock topology on temperature field in compressible turbulence

Effects of two types of shock topology, namely, small-scale shocklet and large-scale shock wave, on the statistics of temperature in compressible turbulence were investigated by simulations. The shocklet and shock wave are caused by the solenoidal and compressive modes of driven forces, respectively. Hereafter, the related two flows are called as SFT and CFT, respectively. It shows that in SFT the temperature spectrum follows the k^-5/3 power law, and the temperature field has "ramp-cliff" structures. By contrast, in CFT the temperature spectrum obeys the k^-2 power law, and the temperature field is dominated by large-scale rarefaction and compression. The power-law exponents for the p.d.f. of large negative dilatation are -2.5 in SFT and -3.5 in CFT, close to theoretical values. For the isentropic assumption of thermodynamic variables, the derivation in SFT grows with the turbulent Mach number (Mt), and for same Mt, the variables in CFT are more anisentropic. The angle statistics shows that the temperature gradient in CFT is preferentially perpendicular to the anisotropic strain rate tensor. It tends to be parallel with the first eigenvector and be orthogonal with the other two eigenvectors. For the cascade of temperature, the temperature variance is increased by the viscous dissipation and the pressure-dilatation, but is decreased by the SGS temperature flux. The statistics of pressure-dilatation proves that the negligible contribution of pressure-dilatation at small scales is due to the cancelations of high values between rarefaction and compression regions. The strongly positive skewness of the p.d.f.s of pressure-dilatation implies that the conversion from kinetic to internal energy through compression is more intense than the opposite process through rarefaction. Furthermore, in SFT the fluctuations of pressure-dilatation approximately follow the Zeman model.

preprint2015arXiv

Transition in Hypersonic Boundary Layers: Role of Dilatational Waves

Transition and turbulence production in a hypersonic boundary layer is investigated in a Mach 6 quiet wind tunnel using Rayleigh-scattering visualization, fast-response pressure measurements, and particle image velocimetry. It is found that the second instability acoustic mode is the key modulator of the transition process. The second mode experiences a rapid growth and a very fast annihilation due to the effect of bulk viscosity. The second mode interacts strongly with the first vorticity mode to directly promote a fast growth of the latter and leads to immediate transition to turbulence.