BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Wei Shyy

Wei Shyy contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.comp-phphysics.flu-dynmath.NANumerical Analysisphysics.app-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 published item(s)

preprint2022arXiv

High-order Compact Gas-kinetic Schemes for Three-dimensional Flow Simulation on Tetrahedral Mesh

A general framework for the development of high-order compact schemes has been proposed recently. The core steps of the schemes are composed of the following. 1). Based on a kinetic model equation, from a generalized initial distribution of flow variables construct a time-accurate evolution solution of gas distribution function at a cell interface ; 2). Introduce the WENO-type weighting functions into the time-derivative of the cell interface flux function in the multistage multi-derivative time stepping scheme to cope with the possible impingement of a shock wave on a cell interface within a time step; 3). Take moments of interface gas distribution function to obtain the time-accurate flow variables and the corresponding fluxes at the cell interface, and update the cell-averaged flow variables and their gradients inside each control volume; 4). Within the physical domain of dependence of the reconstructed cell, based on the cell-averaged flow variables and their gradients develop compact initial data reconstruction to get initial flow distributions at the beginning of next time step. A compact gas-kinetic scheme (GKS) up to sixth-order accuracy in space and fourth-order in time has been constructed on 2D unstructured mesh before. In this paper, the compact GKS up to fourth-order accuracy on 3D tetrahedral mesh will be further constructed with the focus on the WENO-type initial data reconstruction. Nonlinear weights are designed to achieve high-order accuracy for the smooth Navier-Stokes solution and keep super robustness in 3D computation with strong shock interactions. The fourth-order compact GKS can use a large time step with CFL number $0.6$ in the simulations from subsonic to hypersonic flow. A series of test cases are used to validate the scheme. The high-order compact GKS is ready for 3D applications with complex geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Streamwise inclination angle of wall-attached eddies in turbulent channel flows

We develop a new methodology to assess the streamwise inclination angles (SIAs) of the wall-attached eddies populating the logarithmic region with a given wall-normal height. To remove the influences originating from other scales on the SIA estimated via two-point correlation, the footprints of the targeted eddies in the vicinity of the wall and the corresponding streamwise velocity fluctuations carried by them are isolated simultaneously, by coupling the spectral stochastic estimation with the attached-eddy hypothesis. Datasets produced with direct numerical simulations spanning $Re_τ \sim O(10^2)-O(10^3)$ are dissected to study the Reynolds-number effect. The present results show, for the first time, that the SIAs of attached eddies are Reynolds-number dependent in low and medium Reynolds numbers and tend to saturate at $45^{\circ}$ as the Reynolds number increases. The mean SIA reported by vast previous experimental studies are demonstrated to be the outcomes of the additive effect contributed by multi-scale attached eddies. These findings clarify the long-term debate and perfect the picture of the attached-eddy model.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

UGKWP for three-dimensional simulation of gas-particle fluidized bed

The gas-solid particle two-phase flow in a fluidized bed shows complex physics. Following our previous work, the multi-scale framework based on gas-kinetic scheme (GKS) and unified gas-kinetic wave-particle method (UGKWP) for the gas-particle system is firstly extended to the three-dimensional simulation of the fluidized bed. For the solid particle evolution, different from the widely-used Eulerian and Lagrangian approaches, the UGKWP unifies the wave (dense particle region) and discrete particle (dilute particle region) formulation seamlessly according to a continuous variation of particle cell's Kundsen number ($Kn$). The GKS-UGKWP for the coupled gas-particle evolution system can automatically become an Eulerian-Eulerian (EE) method in the high particle collision regime and Eulerian-Lagrangian (EL) formulation in the collisionless particle regime. In the transition regime, the UGKWP can achieve a smooth transition between the Eulerian and Lagrangian limiting formulation. More importantly, the weights of mass distributions from analytical wave and discrete particle are related to the local $Kn$ by $\exp(-1/Kn)$ for wave and $(1-\exp(-1/Kn))$ for discrete particle. As a result, the UGKWP provides an optimal modeling for capturing the particle phase in terms of physical accuracy and numerical efficiency. In the numerical simulation, the UGKWP does not need any prior division of dilute/dense regions, which makes it suitable for the fluidized bed problem, where the dilute/transition/dense regions instantaneously coexist and are dynamically interconvertible. In this paper, based on the GKS-UGKWP formulation two lab-scale fluidization cases are simulated in 3D and the simulation results are compared with the experimental measurements. The typical heterogeneous flow features of the fluidized bed are well captured and the statistics are in good agreement with experiment data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Two-step multi-resolution reconstruction-based compact gas-kinetic scheme on tetrahedral mesh

In this paper, a third-order compact gas-kinetic scheme (GKS) on unstructured tetrahedral mesh is constructed for the compressible Euler and Navier-Stokes solutions. The time-dependent gas distribution function at a cell interface is used to calculate the fluxes for the updating the cell-averaged flow variables and to evaluate the time accurate cell-averaged flow variables as well for evolving the cell-averaged gradients of flow variables. With the accurate evolution model for both flow variables and their slopes, the quality of the scheme depends closely on the accuracy and reliability of the initial reconstruction of flow variables. The reconstruction scheme becomes more challenge on tetrahedral mesh, where the conventional second-order unlimited least-square reconstruction can make the scheme be linearly unstable when using cell-averaged conservative variables alone with von Neumann neighbors. Benefiting from the evolved cell-averaged slopes, on tetrahedral mesh the GKS is linearly stable from a compact third-order smooth reconstruction with a large CFL number. In order to further increase the robustness of the high-order compact GKS for capturing discontinuous solution, a new two-step multi-resolution weighted essentially non-oscillatory (WENO) reconstruction will be proposed. The novelty of the reconstruction includes the following. Firstly, it releases the stability issue from a second-order compact reconstruction through the introduction of a pre-reconstruction step. Secondly, in the third-order non-linear reconstruction, only one more large stencil is added beside those in the second-order one, which significantly simplifies the high-order reconstruction. The proposed third-order scheme shows good robustness in high speed flow computation and favorable mesh adaptability in cases with complex geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Unified gas-kinetic wave-particle method for gas-particle two phase flow from dilute to dense solid-particle limit

In this paper, a unified framework for particulate two-phase flow will be presented with a wide range of solid-particle concentration from dilute to dense limit. The two phase flow is simulated by two coupled flow solvers, i.e., the gas-kinetic scheme (GKS) for the gas phase and unified gas-kinetic wave-particle method (UGKWP) for the particle phase. The GKS is a second-order Navier-Stokes flow solver for the continuum flow. The UGKWP is a multiscale method for all flow regimes. The wave and particle decomposition in UGKWP depends on the cell's Knudsen number (Kn). At a small Kn number, the high concentrated solid particle phase will be modeled by the Eulerian hydrodynamic wave due to the intensive particle-particle collisions. At a large Kn number, the dilute solid particle will be sampled and followed by the Lagrangian particle formulation to capture the non-equilibrium transport. In the transition regime, the distribution and evolution of particle and wave in UGKWP are controlled by the local Kn number with a smooth transition between the above limits. In the current scheme, the two phase model improves the previous one in all following aspects: drag force model for different solid particle concentrations; the frictional pressure in inter-particle contacts at high solid-particle concentration; a flux limiting model to avoid solid particles' over-packing; additional non-conservative nozzle and work terms for the gas phase. Besides, the inter-particle collisions have been refined numerically for the dense particle flow through the discretization of the collision term and numerical flux function. The numerical scheme is tested in a series of typical gas-particle problems. The results validate the accuracy and reliability of the proposed method for gas-particle flow.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Unified gas-kinetic wave-particle methods VI: Disperse dilute gas-particle multiphase flow

In this paper, a unified gas-kinetic wave-particle scheme (UGKWP) for the disperse dilute gas-particle multiphase flow is proposed. The gas phase is always in the hydrodynamic regime. However, the particle phase covers different flow regimes from particle trajectory crossing to the hydrodynamic wave interaction with the variation of local particle phase Knudsen number. The UGKWP is an appropriate method for the capturing of the multiscale transport mechanism in the particle phase through its coupled wave-particle formulation. In the regime with intensive particle collision, the evolution of solid particle will be followed by the analytic wave with quasi-equilibrium distribution; while in the rarefied regime the non-equilibrium particle phase will be captured through particle tracking and collision, which plays a decisive role in recovering particle trajectory crossing behavior. The gas-kinetic scheme (GKS) is employed for the simulation of gas flow. In the highly collision regime for the particles, no particles will be sampled in UGKWP and the wave formulation for solid particle with the hydrodynamic gas phase will reduce the system to the two-fluid Eulerian model. On the other hand, in the collisionless regime for the solid particle, the free transport of solid particle will be followed in UGKWP, and coupled system will return to the Eulerian-Lagrangian formulation for the gas and particle. The scheme will be tested for in all flow regimes, which include the non-equilibrium particle trajectory crossing, the particle concentration under different Knudsen number, and the dispersion of particle flow with the variation of Stokes number. A experiment of shock-induced particle bed fluidization is simulated and the results are compared with experimental measurements. These numerical solutions validate suitability of the proposed scheme for the simulation of gas-particle multiphase flow.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A three-dimensional compact high-order gas-kinetic scheme on structured mesh

In this paper, a third-order compact gas-kinetic scheme is firstly proposed for three-dimensional computation for the compressible Euler and Navier-Stokes solutions. The scheme achieves its compactness due to the time-dependent gas distribution function in GKS, which provides not only the fluxes but also the time accurate flow variables in the next time level at a cell interface. As a result, the cell averaged first-order spatial derivatives of flow variables can be obtained naturally through the Gauss's theorem. Then, a third-order compact reconstruction involving the cell averaged values and their first-order spatial derivatives can be achieved. The trilinear interpolation is used to treat possible non-coplanar elements on general hexahedral mesh. The constrained least-square technique is applied to improve the accuracy in the smooth case. To deal with both smooth and discontinuous flows, a new HWENO reconstruction is designed in the current scheme by following the ideas in Zhu, 2018. No identification of troubled cells is needed in the current scheme. In contrast to the Riemann solver-based method, the compact scheme can achieve a third-order temporal accuracy with the two-stage two-derivative temporal discretization, instead of the three-stage Runge-Kutta method. Overall, the proposed scheme inherits the high accuracy and efficiency of the previous ones in two-dimensional case. The desired third-order accuracy can be obtained with curved boundary. The robustness of the scheme has been validated through many cases, including strong shocks in both inviscid and viscous flow computations. Quantitative comparisons for both smooth and discontinuous cases show that the current third-order scheme can give competitive results against the fifth-order non-compact GKS under the same mesh. A large CFL number around 0.5 can be used in the present scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Comparison of the performance of high-order schemes based on the gas-kinetic and HLLC fluxes

In this paper, a comparison of the performance of two high-order finite volume methods based on the gas-kinetic scheme (GKS) and HLLC fluxes is carried out in structured rectangular mesh. For both schemes, the fifth-order WENO-AO reconstruction is adopted to achieve a high-order spatial accuracy. In terms of temporal discretization, a two-stage fourth-order (S2O4) time marching strategy is adopted for WENO5-AO-GKS scheme, and the fourth-order Runge-Kutta (RK4) method is employed for WENO5-AO-HLLC scheme. For the viscous flow computation, the GKS includes both inviscid and viscous fluxes in the evolution of a single cell interface gas distribution function. While for the WENO5-AO-HLLC scheme, the inviscid flux is provided by HLLC Riemann solver, and the viscous flux is discretized by a sixth-order central difference method. Based on the tests of forward Mach step and viscous shock tube, both schemes show outstanding shock capturing property. From the Titarev-Toro and double shear layer tests, WENO5-AO-GKS scheme seems to have a better resolution than WENO5-AO-HLLC scheme. Both schemes show excellent robustness in extreme cases, such as the Le Blanc problem. From the cases of the Noh problem and the compressible isotropic turbulence, WENO5-AO-GKS scheme shows favorite robustness. In the compressible isotropic turbulence and three-dimensional Taylor-Green vortex problems, WENO-AO-GKS can use a CFL number up to 0.5, instead of 0.3 for WENO5-AO-HLLC. In terms of computational efficiency, WENO5-AO-HLLC scheme is about 27% more expensive than WENO5-AO-GKS scheme in the two-dimensional viscous flow problems, but is about 15% faster in the three-dimensional case. Due to the multi-dimensionality, WENO5-AO-GKS scheme performs better than WENO5-AO-HLLC scheme in the laminar boundary layer and the double shear layer test.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

A Compact Fourth-order Gas-kinetic Scheme for the Euler and Navier-Stokes Solutions

In this paper, a fourth-order compact gas-kinetic scheme (GKS) is developed for the compressible Euler and Navier-Stokes equations under the framework of two-stage fourth-order temporal discretization and Hermite WENO (HWENO) reconstruction. Due to the high-order gas evolution model, the GKS provides a time dependent gas distribution function at a cell interface. This time evolution solution can be used not only for the flux evaluation across a cell interface and its time derivative, but also time accurate evolution solution at a cell interface. As a result, besides updating the conservative flow variables inside each control volume, the GKS can get the cell averaged slopes inside each control volume as well through the differences of flow variables at the cell interfaces. So, with the updated flow variables and their slopes inside each cell, the HWENO reconstruction can be naturally implemented for the compact high-order reconstruction at the beginning of next step. Therefore, a compact higher-order GKS, such as the two-stages fourth-order compact scheme can be constructed. This scheme is as robust as second-order one, but more accurate solution can be obtained. In comparison with compact fourth-order DG method, the current scheme has only two stages instead of four within each time step for the fourth-order temporal accuracy, and the CFL number used here can be on the order of $0.5$ instead of $0.11$ for the DG method. Through this research, it concludes that the use of high-order time evolution model rather than the first order Riemann solution is extremely important for the design of robust, accurate, and efficient higher-order schemes for the compressible flows.

0 citations0 reviewsNo code linkedOpen access