BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Songze Chen

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

ResearcherUnclaimed source record
physics.flu-dynphysics.comp-ph

Catalog footprint

What is connected

7works
2topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

7 published item(s)

preprint2020arXiv

Non-ideal gas dynamics under confinement: rarefaction effect, dense effect and molecular interaction

The effects of volume exclusion and long-range intermolecular attraction are investigated by the simplified kinetic model for surface-confined inhomogeneous fluids. Gas dynamics of the ideal gas, the hard-sphere fluid and the real gas are simulated by the Boltzmann equation, the Enskog equation and the simple kinetic equation, respectively. Only the Knudsen minimum appears for the ideal gas, while both the Knudsen minimum and the Knudsen maximum occur for the hard-sphere fluid and the real gas under certain confinements, beyond which the maximum and minimum may disappear. The Boltzmann equation and the Enskog equation overestimates and underestimates the mass flow rate of the real gas dynamics under confinement, respectively, where the volume exclusion and the long-range intermolecular attractive potential among molecules are not ignorable. With the increase of the channel width, gas dynamics of the hard-sphere fluid and the real gas tends to the Boltzmann prediction gradually. The density inhomogeneity, which hinders the flow under confinement, is more obvious when the solid fraction is larger. The anomalous slip occurs for real gas under constant confinement. The flow at a smaller Knudsen number (larger solid fraction or channel width) contributes more practical amount of mass transfer, although the rarefaction effects is more prominent at larger Knudsen numbers. The temperature has no effect on density and velocity profiles of the ideal gas and the hard-sphere fluid, but the energy parameter among the real gas molecules decreases with the increasing temperature and the real gas dynamics tends to the hard-sphere ones consequently.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A well-balanced gas kinetic scheme for Navier-Stokes equations with gravitational potential

The hydrostatic equilibrium state is the consequence of the exact hydrostatic balance between hydrostatic pressure and external force. Standard finite volume or finite difference schemes cannot keep this balance exactly due to their unbalanced truncation errors. In this study, we introduce an auxiliary variable which becomes constant at isothermal hydrostatic equilibrium state and propose a well-balanced gas kinetic scheme for the Navier-Stokes equations with a global reconstruction. Through reformulating the convection term and the force term via the auxiliary variable, zero numerical flux and zero numerical source term are enforced at the hydrostatic equilibrium state instead of the balance between hydrostatic pressure and external force. Several problems are tested numerically to demonstrate the accuracy and the stability of the new scheme, and the results confirm that, the new scheme can preserve the exact hydrostatic solution. The small perturbation riding on hydrostatic equilibria can be calculated accurately. The viscous effect is also illustrated through the propagation of small perturbation and the Rayleigh-Taylor instability. More importantly, the new scheme is capable of simulating the process of converging towards hydrostatic equilibrium state from a highly non-balanced initial condition. The ultimate state of zero velocity and constant temperature is achieved up to machine accuracy. As demonstrated by the numerical experiments, the current scheme is very suitable for small amplitude perturbation and long time running under gravitational potential.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

GPU acceleration of an iterative scheme for gas-kinetic model equations with memory reduction techniques

This paper presents a Graphics Processing Units (GPUs) acceleration method of an iterative scheme for gas-kinetic model equations. Unlike the previous GPU parallelization of explicit kinetic schemes, this work features a fast converging iterative scheme. The memory reduction techniques in this method enable full three-dimensional (3D) solution of kinetic model equations in contemporary GPUs usually with a limited memory capacity that otherwise would need terabytes of memory. The GPU algorithm is validated against the DSMC simulation of the 3D lid-driven cavity flow and the supersonic rarefied gas flow past a cube with grids size up to 0.7 trillion points in the phase space. The performance of the GPU algorithm is assessed by comparing with the corresponding parallel CPU program using Message Passing Interface (MPI). The profiling on several models of GPUs shows that the algorithm has a medium to high level of utilization of the GPUs' computing and memory resources. A $190\times$ speedup can be achieved on the Tesla K40 GPUs against a single core of Intel Xeon-E5-2680v3 CPU for the 3D lid-driven cavity flow.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A simplification of the unified gas kinetic scheme

Unified gas kinetic scheme (UGKS) is an asymptotic preserving scheme for the kinetic equations. It is superior for transition flow simulations, and has been validated in the past years. However, compared to the well known discrete ordinate method (DOM) which is a classical numerical method solving the kinetic equations, the UGKS needs more computational resources. In this study, we propose a simplification of the unified gas kinetic scheme. It allows almost identical numerical cost as the DOM, but predicts numerical results as accurate as the UGKS. Based on the observation that the equilibrium part of the UGKS fluxes can be evaluated analytically, the equilibrium part in the UGKS flux is not necessary to be discretized in velocity space. In the simplified scheme, the numerical flux for the velocity distribution function and the numerical flux for the macroscopic conservative quantities are evaluated separately. The simplification is equivalent to a flux hybridization of the gas kinetic scheme for the Navier-Stokes (NS) equations and conventional discrete ordinate method. Several simplification strategies are tested, through which we can identify the key ingredient of the Navier-Stokes asymptotic preserving property. Numerical tests show that, as long as the collision effect is built into the macroscopic numerical flux, the numerical scheme is Navier-Stokes asymptotic preserving, regardless the accuracy of the microscopic numerical flux for the velocity distribution function.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cartesian Grid Method for Gas Kinetic Scheme

A Cartesian grid method combined with a simplified gas kinetic scheme is presented for subsonic and supersonic viscous flow simulation on complex geometries. Under the Cartesian mesh, the computational grid points are classified into four different categories, the fluid point, the solid point, the drop point, and the interpolation point. The boundaries are represented by a set of direction-oriented boundary points. A constrained weighted least square method is employed to evaluate the physical quantities at the interpolation points. Different boundary conditions, including isothermal boundary, adiabatic boundary, and Euler slip boundary, are presented by different interpolation strategies. We also propose a simplified gas kinetic scheme as the flux solver for both subsonic and supersonic flow computations. The methodology of constructing a simplified kinetic flux function can be extended to other flow systems. A few numerical examples are used to validate the Cartesian grid method and the simplified flux function. The reconstruction scheme for recovering the boundary conditions of compressible viscous and heat conducting flow with a Cartesian mesh can provide a smooth distribution of physical quantities at solid boundary, and present an accurate solution for the flow study with complex geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Comparative Study of an Asymptotic Preserving Scheme and Unified Gas-kinetic Scheme in Continuum Flow Limit

Asymptotic preserving (AP) schemes are targeting to simulate both continuum and rarefied flows. Many AP schemes have been developed and are capable of capturing the Euler limit in the continuum regime. However, to get accurate Navier-Stokes solutions is still challenging for many AP schemes. In order to distinguish the numerical effects of different AP schemes on the simulation results in the continuum flow limit, an implicit-explicit (IMEX) AP scheme and the unified gas kinetic scheme (UGKS) based on Bhatnagar-Gross-Krook (BGk) kinetic equation will be applied in the flow simulation in both transition and continuum flow regimes. As a benchmark test case, the lid-driven cavity flow is used for the comparison of these two AP schemes. The numerical results show that the UGKS captures the viscous solution accurately. The velocity profiles are very close to the classical benchmark solutions. However, the IMEX AP scheme seems have difficulty to get these solutions. Based on the analysis and the numerical experiments, it is realized that the dissipation of AP schemes in continuum limit is closely related to the numerical treatment of collision and transport of the kinetic equation. Numerically it becomes necessary to couple the convection and collision terms in both flux evaluation at a cell interface and the collision source term treatment inside each control volume.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Comparison and Unification of Ellipsoidal Statistical and Shakhov BGK Models

The Ellipsoidal Statistical model (ES-model) and the Shakhov model (S-model) are constructed for the correction of Prandtl number of the original BGK model through the modification of stress and heat flux. Even though in the continuum flow regime, both models can give the same Navier-Stokes equations with correct Prandtl number, their modification of the collision term may have different dynamic effect in the non-equilibrium transition flow regimes. With the introduction of one free parameter, a generalized kinetic model with the combination of the ES-model and S-model can be developed, and this new model can get the correct Navier-Stokes equations in the continuum flow regime as well, but with abundant dynamic effect through the adjustment of the new degree of freedom. In order to validate the generalized model, a numerical method based on the unified gas kinetic scheme (UGKS) has been developed for the new model. The physical performance of the new model with the variation of the free parameter has been tested, where the ES-model and S-model become the limiting cases. In transition flow regime, many physical problems, i.e., the shock structure and micro-flows, have been studied using the generalized model. With a careful choice of the free parameter, good results can be achieved for most test cases. The overall conclusion is that the S-model predicts more accurate numerical solutions in most tough test cases presented in this paper than the ES-model, while ES-model performs better in the cases when the flow is mostly driven by heat, such as a channel flow with large boundary temperature variations at high Knudsen number. The numerical study demonstrates the necessity of developing such a generalized model. With the inclusion of one more freedom, in the transition regime the new kinetic model may provide more accurate solution than the ES and Shakhov models.

0 citations0 reviewsNo code linkedOpen access